可压缩自由混合层中的输运平衡

为研究可压缩混合层中的输运特性 ,该文采用气动BGK格式模拟了超音速来流条件下的二维平面可压缩混合层的空间发展 ,给出了湍能、雷诺应力和动量输运方程中各项的分布情况。各平衡方程中粘性耗散项和粘性扩散项都很小。湍能主要是通过对流、湍流扩散和压力扩散进行输运 ;切应力输运方程中主要是生成项与压力扩散项和压力 -应变关联项的平衡 ;雷诺应力 (1,1)分量输运方程中以生成项和压力 -应变关联项为主 ;雷诺应力 (2 ,2 )分量输运方程中主要是压力 -应变关联项与湍流扩散项和压力扩散项之间的平衡。     

泡状流相分布及湍流结构的欧拉-拉格朗日双向耦合数值研究

使用欧拉-拉格朗日双向耦合数值模拟方法对泡状流的相分布及气泡对液相湍流产生的影响做了详细的调查.液相的速度场采用直接数值模拟方法求解,气泡的运动轨迹由牛顿运动方程计算.相间的耦合是把相间作用力作为液相动量方程的源相来实现的,所考虑的相间作用力有阻力、剪切升力、惯性力等.研究表明:阻力、浮力、剪切升力和虚拟质量力对气泡的运动具有重要影响,大量的气泡在剪切升力的作用下聚集在壁面附近;与单液体相相比,气泡的加入使液相的平均速度、湍流强度和雷诺应力均有一定幅度的降低,而且气泡的加入修改了阻力系数公式.     

旋转对平板通道内湍流流动和换热影响的直接模拟

直接模拟了展向旋转数等于0~15时平板通道内的湍流流动和换热,分析了科氏力对湍流平均特性、湍流统计特性、雷诺应力输运和湍流结构的影响.采用有限差分法离散非稳态N S方程,湍流雷诺数和普朗特数分别为150和0 71,网格数为64×128×64.计算结果表明,随着旋转数的增加,压力面附近湍流强度增大,换热增强,而吸力面附近换热减小.雷诺应力输运方程中的旋转项在各应力中的能量分配作用远大于压力分配项,且旋转项使得雷诺应力输运方程的各项分布和相对关系发生了较大的变化.由于科氏力的作用,压力面附近的条带间距减小,条带变细、变短,而吸力面流动脉动衰减,有层流化的趋势.     

基于无结构网格的有限体积法和颗粒方法相结合求解湍流的联合PDF方程

本研究着力于发展一种新的有限体积法和颗粒方法的混合算法 ,用于湍流 (冷态流以及反应流 )的速度 频率 标量联合PDF(probabilitydensityfunction)输运模型方程的数值求解 .采用一种基于无结构网格的有限体积法 ,求解质量、动量、能量平均守恒方程和平均状态方程 ,为联合PDF的输运模型方程的求解提供其所需的平均速度场和平均密度场 ;颗粒方法中 ,用MonteCarlo算法来求解PDF方程 ,为平均场的计算提供能源源项、雷诺应力、热通量、比热比等信息 .两种方法之间的耦合是一个关键 ,为了验证该算法的可行性 ,论文还对二维V型槽后的湍流冷态流场进行了数值模拟 ,并与实验、大涡模拟、k 模型的结果进行了比较 ,结果表明 ,该算法的结果与实验结果吻合得较好 .     

用颗粒相双尺度二阶矩湍流模型模拟突扩两相流动

基于将颗粒脉动分成湍流引起的大尺度脉动和颗粒间碰撞产生的小尺度脉动的概念，应用单相流动多尺度湍流模型方法，从双流体模型出发，推导和封闭了颗粒的双尺度脉动雷诺应力方程、大尺度脉动能量传递率方程、小尺度脉动耗散率方程和两相脉动关联方程，用统一的湍流模型方法构造了包含颗粒两类脉动的新的颗粒双尺度二阶矩两相湍流模型．用该模型对突扩管内的两相流动进行了数值模拟，并和单尺度的二阶矩两相湍流模型进行了比较，计算结果和实验数据吻合较好，比单尺度的二阶矩两相湍流模型的计算结果有所改进．     

矢势A的物理意义探讨

本文在前文(1)的基础上进一步讨论了矢势A的物理意义。指出矢势A是电磁场的场动量势。它和标势φ一起作为电磁场的基本场量导出电动力学方程比以E、B为基本场量的电动力学方程更为基本。这也完全符合近代物理中力的概念正在消失,能量和动量的概念才是近代物理学中最重要的概念这一情况。文中还着重说明,利用规范变换决不能消除A-B效应。同时指出,矢势和标势在理论上存在着规范变换的不确定性与场动量势和测量矢势时的唯一性是统一的。     

旋转矩形通道内湍流雷诺应力分析

采用直接数值模拟方法对旋转矩形通道内的湍流流动和换热进行了研究.非稳态N-S方程的空间离散采用二阶中心差分法,时间推进采用二阶显式Adams-Bashforth格式.通过分析近壁湍流雷诺应力输运方程各项的分布以及比较壁面附近湍流生成和耗散强度,讨论了科氏力对湍流脉动的影响.结果表明:系统旋转效应增强了压力面附近的湍流生成、黏性耗散、黏性扩散、湍流扩散及压力应变相关项和压力速度扩散项,而在吸力面附近,雷诺应力输运方程中各项因旋转而受到强烈抑制.     

沿岸海域三维斜压场的数值模拟：Ⅱ.渤海湍流能量数值计算

研究建立渤海海域三维斜压场的流能量模型。对计算网格无法分辩的湍流运动，引入湍流能量得垂向动量、质量及能量交换系统，从而形成封闭的湍流运动方程组。     

涡流发生器对3/4汽车模型风洞流场品质的影响

采用数值计算的方法对3/4开口回流式汽车风洞喷口处的涡流发生器进行了研究.对比研究了三种不同片数的涡流发生器对流场品质的影响.首先采用定常雷诺时均纳维斯托克斯方程求解了流场的定常特性,接着采用大涡模拟(LES)方法对流场的非定常特性进行了研究.并就两种计算方法给予了相应的试验验证.对比计算结果得到涡流发生器可以在一定程度上提高流场均与性,尤其是三片式涡流发生器的工况能较好地拓宽风洞测试区的高流速区域的范围,降低测试区的静压梯度,减小试验段内的湍流度和剪切层内的湍流度,延长测试段内低湍流区的长度,分散流场内的脉动能量,降低低频颤振敏感频率区的能量聚集.     

可压缩湍流封闭模型初探：湍流动能K方程

为解决可压缩流的湍流运动方程组封闭问题，对瞬时值的Ｎａｖｉｅｒ—Ｓｔｏｋｅｓ运动方程进行质量加权平均处理，根据湍流封闭理论对各脉动量关联项进行模化，得到可压缩流动的湍流动能Ｋ方程，完成了建立一个适用于三维可压缩流的湍流封闭模型的第一步工作．     

考虑井筒变质量流动的砾石充填水平井产能预测

应用质量守恒和动量守恒原理推导砾石充填水平井井筒内变质量流动压降方程,并采用拟三维思想,将地层内流体在三维空间的流动分为水平面内的向垂直裂缝流和近井区域垂直平面内的径向流,建立油藏渗流模型.提出将井筒变质量流与油藏渗流耦合的数学模型及求解方法.结果表明:利用该模型计算的水平井产能与实测结果平均误差仅为3.79%,沿井筒...     

水平井砾石充填过程中的变质量流

由于携砂液向地层的滤失以及通过筛管缝隙的流体质量交换,水平井砾石充填过程中冲筛、井筒环空两个流动系统均为复杂的变质量流。考虑变质量流特征,分别建立井筒环空、冲筛环空两个独立流动系统的砾石、携砂液的质量和动量守恒方程,以及各系统间的流动耦合方程,得到水平井砾石充填过程中的变质量流动数学模型。数值求解后可计算压力、流量、窜流量、滤失量等随时间和位置的变化规律,为水平井砾石充填的优化设计提供基础。     

非等温湍流射流的脉动关联

本文采用X热线-冷线组合探针,对具有高湍流强度的非等温湍流圆形射流的流场和温度场进行了较精确的测量,尤其是对脉动量及其关联关系的实测和分析不仅在湍流测量方法上具有一定的新意,还为湍流结构、动量与热量交换的类比关系以及相关学科的研究,提供了某些有用的依据.     

Hydrodynamic turbulence cannot transport angular momentum effectively in astrophysical disks

The most efficient energy sources known in the Universe are accretion disks. Those around black holes convert 5-40 per cent of rest-mass energy to radiation. Like water circling a drain, inflowing mass must lose angular momentum, presumably by vigorous turbulence in disks, which are essentially inviscid. The origin of the turbulence is unclear. Hot disks of electrically conducting plasma can become turbulent by way of the linear magnetorotational instability. Cool disks, such as the planet-forming disks of protostars, may be too poorly ionized for the magnetorotational instability to occur, and therefore essentially unmagnetized and linearly stable. Nonlinear hydrodynamic instability often occurs in linearly stable flows (for example, pipe flows) at sufficiently large Reynolds numbers. Although planet-forming disks have extreme Reynolds numbers, keplerian rotation enhances their linear hydrodynamic stability, so the question of whether they can be turbulent and thereby transport angular momentum effectively is controversial. Here we report a laboratory experiment, demonstrating that non-magnetic quasi-keplerian flows at Reynolds numbers up to millions are essentially steady. Scaled to accretion disks, rates of angular momentum transport lie far below astrophysical requirements. By ruling out purely hydrodynamic turbulence, our results indirectly support the magnetorotational instability as the likely cause of turbulence, even in cool disks.     

地下储气库注采井温度压力耦合分析

地下储气库进行循环注气、采气作业,引起注采井筒温度、压力分布的交替变化,准确预测井筒温度、压力分布是优化设计管柱的基础。基于气体质量、动量、能量平衡方程和地层-流体径向非稳态传热模型,考虑气体的焦耳-汤姆逊效应,得到了井筒气体温度、压力、密度和流速的耦合微分方程组,采用四阶龙格-库塔方法数值求解。通过实例分析了注采工况下井筒的温度、压力分布特征及影响因素,并与Well Cat软件计算结果相比较,结果表明注气井口温度、注气压力、注采流量和作业时间影响井筒的温度压力分布。注气中,注入压力对井底压力的影响最大,对井底温度的影响很小;采气中,井口温度随流量的增加而增大,井口压力随流量的增大而减小;作业早期温度、压力随时间变化较大。     

高Reynolds数流动的时均能量方程及熵增表达式

本文对湍流耗散过程作了详细的热力学分析,明确提出不仅在时均动量方程中,而且在时均能量方程及热力学第二定律的表达式中,湍流影响项与粘性影响项也具有相拟的张量表达式,而且与用何种湍流模型计算Reynolds 应力张量Ⅱ_r 无关.本文还引进名义滞止焓和平均流动熵的概念,推导出适合数值求解的封闭的克罗柯形式的时均流动方程组.     

岩体THM耦合模型控制方程建立

根据岩体三场耦合的机理分析,从线性动量守恒、质量守恒和能量守恒出发,理论上推导出以位移、孔隙压力和温度为未知量的THM耦合控制方程组,包括岩块变形场方程、地下水连续性方程和水、岩能量守恒方程.方程组不仅考虑了水、岩的密度随温度和压力变化,同时也考虑了水、岩温度不同具有热交换,最后还对后继研究提出了建议.     

射孔完井水平井筒变质量湍流压降规律研究

运用拟三维原理 ,提出了用于计算地层刚性渗流、弹性不稳定渗流与水平井筒内变质量湍流耦合的新方法 ,解决了在井筒变质量流体流动模拟过程中井筒边界条件难以确定、地层线性渗流与井筒内非线性湍流难以耦合的问题。同时 ,对不同射孔方式和生产制度下流体水动力学进行了分析 ,研究了射孔入流与水平井筒内流体的相互作用关系及变质量主流体的流动特性。模拟结果表明 ,与圆管湍流相比 ,井筒内变质量流具有更大的压降系数 ,且主流中心速度降低 ,速度剖面呈现内凹形式 ,水平井射孔设计应该利用井筒内流体与渗流的耦合模拟进行优化。     

拉曼耦合作用下玻色-爱因斯坦凝聚的磁化性质

研究双分量玻色气体在拉曼耦合作用下的磁化性质。对于宏观占据单一动量态的玻色?爱因斯坦凝聚, 磁化强度和磁化率分别与凝聚动量和有效质量存在普适的线性关系, 这些普适关系不依赖于拉曼耦合的强度、失谐大小及原子间的相互作用。推导出磁化率与有效质量的线性关系, 可以在经典物理图像下描述简谐势阱中凝聚体的偶极振荡, 自旋空间、动量空间及坐标空间的振幅之比可表达为有效质量和自旋磁化率的解析形式。这组振幅关系对于亚稳态附近的小幅振荡也同样适用。     

叶片泵内气液两相流的三维流动数值模型

基于双流体模型建立针对叶片泵内气液两相三维湍流流动的数学模型和数值计算方法.相间作用力考虑了阻力和附加质量力,湍流模型考虑了因相含量脉动引起的附加源项,在SIMPLEC算法基础上提出了一种气液两相流的压力修正算法.运用所建模型计算了叶片式气液混输泵在进口含气率分别为5%,15%,25%的系列工况下叶轮的内部两相流场,并作了扬程特性预测.与实验结果的对比表明,所提出的流动模型是可靠的,并具有较宽的进口含气率适用范围.