基于小关联时间两相湍流色噪声的近似模型及数值模拟

从一般高斯型色噪声模型出发,通过泛函导数,应用小关联时间,近似计算多维色噪声,得到有效Fokker-Planck方程.将其应用到两相湍流中得到颗粒相的概率密度函数输运方程,从而得到颗粒相的二阶矩模型.将颗粒应力方程简化成代数方程,建立代数应力模型.将对流扩散方程的有限分析法运用到求解两相流模型中,对壁面两相射流进行数值模拟,并将求解结果与实验结果进行对比分析.     

下降管中稠密两相湍流的数值模拟

本文将统一二阶矩两相湍流模型和颗粒动力学理论结合,推导并封闭了稠密两相流考虑颗粒间碰撞的统一二阶矩两相湍流模型。该模型用颗粒动力学理论模拟颗粒之间的碰撞,用各向异性的统一二阶矩模型考虑气相和颗粒相的湍流脉动,并用输运方程描述气固两相湍流之间的相互作用。最后用该模型对下降管中的气粒两相流动进行了模拟,模拟所得的沿径向颗粒浓度分布和速度分布与实验吻合较好,预报结果也反应出了颗粒雷诺应力的各向异性。     

用视密度加权平均二阶矩模型模拟旋流两相流动

本文用视密度加权平均代替时平均,建立了视密度加权平均的统一二阶矩两相湍流模型方程组（ＭＵＳＭ）,其中用体积分数代替了数密度,用颗粒驰豫时间作为封闭两相脉动速度关联方程耗散项的时间尺度,并引入了颗粒视在的气体速度脉动的输运方程。用ＭＵＳＭ模型模拟了旋流数为０．４７的气粒两相流动。并和实验结果及时间平均的ＵＳＭ模型的模拟结果进行了对照,两种模型均能较好地预报的两相的轴向和切向速度,轴向和切向脉动速度。此外,ＭＵＳＭ模型可以减少所用方程数,节省计算量。因此视密度加权平均的统一二阶矩两相湍流模型是一种对时间平均的统一二阶矩模型的改进,今后可以进一步扩大应用。     

颗粒Reynolds正应力轨道模型及后台阶两相流动数值模拟

用随机过程理论建立气固两相耦合脉动量Lagrange方程,并建立了气固两相流随机颗粒轨道模型中颗粒Reynolds正应力的Lagrange方程.将新的模型运用于各向同性的湍流衰减的流场中,模拟颗粒的湍流扩散特性,与文[1]中的模型和实验结果作比较,并使用新模型模拟了后台阶两相流动.     

k－ε－PDF两相湍流模型和台阶后方气粒两相流动的模拟

本文提出了两相湍流的ｋ－ε－ＰＤＦ模型。ＰＤＦ模型所得的湍流两相流动的统计平均方程，与雷诺时均方程有相似的形式，但ＰＤＦ模型可精确计算出颗粒相各方程的脉动关联项。本文将该模型用于预报台阶后方湍流两相流动，与ｋ－ε－ｋｐ模型相比，它可以更合理地预报出颗粒湍流的各向异性。     

DSM-LPDF两相湍流模型及旋流两相流动的模拟

本文由流体－颗粒速度的拉氏联合概率密度函数（ＰＤＦ）输运方程出发,用Ｓｉｍｏｎｉｎ建议的Ｌａｎｇｅｖｉｎ模型封闭颗粒所遇到流体瞬时速度的条件期望项,并用Ｍｏｎｔｅ Ｃａｒｌｏ方法直接求解 ＰＤＦ输运方程,将其和求解流体雷诺应力方程模型的有限差分方法结合,建立了雷诺应力－拉氏ＰＤＦ（ＤＳＭ－ＬＰＤＦ,简称ＤＬ）两相湍流模型．用此模型模拟了旋流数为０．４７的突扩旋流气粒两相流动,并与文献中ＰＤＰＡ实验和用类似于单相流动湍流模型封闭方法的时平均统一二阶矩（ＵＳＭ）模型的预报进行了对比．     

基于Lagrange方法封闭的两相湍流场方程模型

两相湍流场方程模型采用基于Euler方法的一阶矩方程,而二阶矩方程由Lagrange方法得到,新模型的封闭不需要附加其它假设.首先基于概率密度函数给出颗粒运动的连续方程和动量方程,其次由基于平均Langevin方程的Lagrange模型推导得到颗粒二阶矩方程,最终获得封闭的二阶矩模型.将新模型用于气固两相壁面射流的数值模拟,结果表明新模型合理有效.     

湍流两相流动有燃烧颗粒相概率密度函数输运方程理论

由有燃烧的湍流气粒两相流动的瞬态方程和统计力学概率密度函数概念出发,推导了有燃烧颗粒相的质量－动量－能量联合概率密度函数（ＰＤＦ）输运方程,并对方程中条件期望项用梯度模拟概念进行了模拟封闭。封闭后的ＰＤＦ方程可作为建立颗粒拟流体模型方程和封闭二阶矩模型的基础,也可以通过Ｍｏｎｔｅ－Ｃａｒｌｏ 法求解用以直接计算颗粒雷诺应力和湍流动能,以便和二阶 矩模型的结果相对照,改善二阶矩模型。     

颗粒增强湍流的新模型和气粒流动的湍流变动

用雷诺应力方程模型和极细的网格系对单个颗粒受湍流气体绕流进行了数值模拟,研究了改变颗粒直径和气体相对速度时颗粒增强气体湍流的规律.据此构造了颗粒尾涡增强气体湍流的新模型.将此子模型加入到两相流动模型中,对竖直和水平通道内气粒两相流动进行了数值模拟,和实验结果的对照表明,考虑颗粒尾涡增强气体湍流效应得到的气体湍流脉动速度的模拟结果比不考虑此效应的模拟结果好得多.     

基于拉氏概率密度的两相湍流二阶矩模型

用颗粒运动的拉氏分析和PDF方法,改进了颗粒相的二阶矩模型。由拉氏两相运动的随机微分方程出发,采用随机过程分析和信号分析法得到湍流两相流动的PDF输运方程,双流体模型方程和两相脉动速度相关的基本模式的封闭式,和用其它方法导出的方程与封闭式的结果一致,对封闭式作了重要的改进,在分析颗粒轨道上的流体湍流作用时间时,全面地引入拉氏分析的轨道穿越效应、惯性效应、连续效应和湍流的各向异性。     

基于水平集方法二维模型的软性磨粒两相流流场特性分析方法

两相流场特性分析是软性磨粒流精密加工方法的技术关键. 针对此问题, 利用水平集方法(LSM)结构拓扑变换原理, 建立了描述磨粒流在另一流体中运动变化的二维模型. 该模型中, 用LSM捕捉两相流之间的界面, 从而对两相流的运动情况和相关参数进行模拟仿真, 深入讨论了磨粒流加工机理. 研究表明: 以LSM原理建立的二维模型, 成功地仿真模拟了两相流相变过程; 利用k-ε湍流模型与Preston方程相结合的数值分析方法, 对流道内各区域的速度等流场特性进行讨论, 为软性磨粒流加工提供了前期基础性研究. 关键词： 水平集方法 结构化表面 软性磨粒流 k-ε湍流模型')" href="#">k-ε湍流模型     

两相湍流色噪声扩维法PDF模型及数值模拟

由颗粒运动的朗之万方程出发,对流体脉动速度采用扩维方法,得到两个不同层次的PDF输运方程.通过对颗粒运动方程求解和高斯分布假设,解决PDF方程的封闭问题,获得颗粒二阶矩模型,然后将颗粒应力方程简化成代数方程,建立代数应力模型.将对流扩散方程的有限分析法运用到求解两相流模型中,对壁面两相射流进行数值模拟,对比分析数值结果与实验结果.     

Van der Waals状态方程可压缩多流体流动的PPM方法

基于流体体积分数的混合型多流体数值模型，将Piecewise Parabolic Method（PPM）方法应用于可压缩多流体流动的数值模拟，采用双波近似求解多流体van der Waals状态方程的Riemann问题．模拟高密度比且含有激波的可压缩多流体流动，典型的纯界面平移问题模拟结果表明，在接触间断的界面附近，压力和速度没有任何的振荡且界面数值耗散都被控制在2—3个网格之内；一维和二维算例表明，该数值方法可以有效地处理接触间断、激波和多维滑移线等物理问题，并能够比其它多流体数值方法更精细地模拟多流体交界面．     

Analysis of Two-Phase Cavitating Flow with Two-Fluid Model Using Integrated Boltzmann Equations

In the present work, both computational and experimental methods are employed to study the two-phase flow occurring in a model pump sump. The two-fluid model of the two-phase flow has been applied to the simulation of the three-dimensional cavitating flow. The governing equations of the two-phase cavitating flow are derived from the kinetic theory based on the Boltzmann equation. The isotropic RNG$k-\epsilon-k_{ca}$ turbulence model of two-phase flows in the form of cavity number instead of the form of cavity phase volume fraction is developed. The RNG $k-\epsilon-k_{ca}$ turbulence model, that is the RNG$k-\epsilon$ turbulence model for the liquid phase combined with the $k_{ca}$model for the cavity phase, is employed to close the governing turbulent equations of the two-phase flow. The computation of the cavitating flow through a model pump sump has been carried out with this model in three-dimensional spaces. The calculated results have been compared with the data of the PIV experiment. Good qualitative agreement has been achieved which exhibits the reliability of the numerical simulation model.     

Dynamics of the Langevin model subjected to colored noise: Functional-integral method

We have discussed the dynamics of Langevin model subjected to colored noise, by using the functional-integral method (FIM) combined with equations of motion for mean and variance of the state variable. Two sets of colored noise have been investigated: (a) one additive and one multiplicative colored noise, and (b) one additive and two multiplicative colored noise. The case (b) is examined with relevance to a recent controversy on the stationary subthreshold voltage distribution of an integrate-and-fire model including stochastic excitatory and inhibitory synapses and a noisy input. We have studied the stationary probability distribution and dynamical responses to time-dependent (pulse and sinusoidal) inputs of the linear Langevin model. Model calculations have shown that results of the FIM are in good agreement with those of direct simulations (DSs). A comparison is made among various approximate analytic solutions such as the universal colored noise approximation (UCNA). It has been pointed out that dynamical responses to pulse and sinusoidal inputs calculated by the UCNA are rather different from those of DS and the FIM, although they yield the same stationary distribution.     

α铁冲击相变的离散元方法

利用离散元方法结合基于无扩散相变的两相模型、热力学相容的自由能函数与有限速率相变动力学方程,模拟α铁的冲击相变过程.得到铁的相边界、冲击Hugoniot关系及自由面速度剖面,并对不同冲击压力下加卸载过程的波系传播和相互作用过程进行分析,获得了加卸载波速特征.     

A moment approach to non-Gaussian colored noise

The Langevin system subjected to non-Gaussian colored noise has been discussed, by using the second-order moment approach with two kinds of models for generating the noise. We have derived the effective differential equation (DE) for a variable x, from which the stationary probability distribution P(x) has been calculated with the use of the Fokker-Planck equation. The result of P(x) calculated by the moment method is compared to several expressions obtained by different methods such as the universal colored noise approximation (UCNA) [Jung and Hänggi, Phys. Rev. A 35 (1987) 4464] and the functional-integral method. It has been shown that our P(x) is in good agreement with that of direct simulations (DSs). We have also discussed dynamical properties of the model with an external input, solving DEs in the moment method.     

染料激光色立方模型的第一通过时间分布

本文首先把Ｊｕｎｇ和Ｈａｎｇｇｉ提出的统一色噪声近似推广到多维情形，从而导出了这一近似下单模染料激光色立方模型的朗之万方程和福克－普朗克方程；得到了该模型第一通过时间分布的平均值，协方差和偏斜率的近似解析式并将其和数值模拟进行比较。