运动激波与气泡串相互作用的多介质数值模拟

通过对激波和流体界面相互作用而诱导的大变形界面演化的数值模拟,验证了Level set 方法精确模拟多个流体界面的有效性.采用间断有限元Galerkin方法求解欧拉方程得到流场解,采用5阶WENO格式求解Level set方程追踪多流体界面,界面附近的边界条件由虚拟流体方法处理.对运动激波和两个气泡相互作用过程进行了数值模拟,得到了不同时刻的压力和密度等值线分布,并分析了计算域中两个气泡同是氦气泡,以及一个是氦气泡,一个是R22气泡情况下的计算结果.计算结果表明:利用多界面Level set方程可高质量地捕捉多个流体界面,处理3种多介质流场数值模拟问题.     

Stiffened状态方程下γ-model方法跟踪二维多介质流动

在Stiffened状态方程下，运用Level—set方程跟踪界面运动变化，把界面捕捉的等效方程、Level-set函数和欧拉方程组耦合，求解包含比热比和材料参数的耦合形式的流体力学方程组。计算方法采用二阶精度Wave-propagation算法，通过对气体-液体两种流体的一维Riemann问题、激波和气泡相互作用以及气液两种流体Richtmyer—Meshkov界面不稳定性问题进行数值模拟，抑制了不同介质界面两侧的非物理振荡。     

用虚拟区域法结合虚拟边界条件求解对称及管壁流动问题

分析了用基于Lagrange乘子的虚拟区域法数值求解时,翼型后缘存在的奇异问题.为了保持直角坐标网格结构简单及可利用快速算法的特点,提出引入虚拟边界条件控制虚拟流动的方法处理奇异问题.用虚拟区域法数值模拟了翼型不可压对称绕流;并把该方法推广应用到管壁流动的数值模拟,给出了喷管内流场的算例.数值结果表明,虚拟边界条件的引入有助于控制Lagrange乘子的变化,从而改善数值模拟的精度.     

基于预处理的低速小展弦比机翼数值模拟

以雷诺平均的Navier-Stokes 方程为基础, 应用非结构混合网格技术,基于有限体积的时间推进预处理方法,对小展弦比机翼的低速流动进行了数值模拟.其中对流项的离散应用了基于Roe的Riemann近似解迎风格式,时间离散应用了基于LU-SGS的隐式格式.为描述湍流流动,采用了Spalart-Allmaras一方程湍流模型.为了验证本文方法对低速小展弦比机翼粘性流动的数值分析效果, 对几个不同展弦比的机翼进行了数值模拟.本文的数值结果与文献的实验数据基本一致, 表明本文方法能够准确模拟小展弦比机翼在低速流动时机翼的展弦比对其气动特性的影响.     

边界条件对三维空腔流动振荡的影响

基于全隐式无分裂方法,数值求解三维非定常雷诺平均的N av ier-S tokes方程,紊流模型采用修正的B-L代数模型。运用此方法,本文数值模拟了三维开式空腔的超声速、跨声速、亚声速流动。研究了空腔口处剪切层的非定常运动现象,着重探讨了入口边界条件对三维空腔流自持性振荡的影响。数值计算结果表明,与层流速度型相比,入口条件为紊流速度型和特征边界型得到的结果能够吻合。     

竖直上升通道内微气泡湍流流动的数值模拟

为了提高气泡浓度分布预测精度,考虑流体和气泡之间的双向耦合作用,对竖直上升通道内含有微小气泡的湍流流动进行了数值模拟研究。采用直接数值模拟方法模拟流体流动,并用拉格朗日法跟踪气泡受流体和重力作用所产生的运动轨迹。气泡受力包括相间阻力、虚质量力、升力、壁面升力、压力梯度和重力等。数值模拟结果表明,通过加入壁面升力和采用双向耦合方法,可以较准确地预测气泡浓度分布和揭示气泡流动对流体湍流结构的影响。     

关于粒子图像测速仪精度的数值测试

利用数值模拟的方法对粒子图像测速仪的精度问题进行了数值实验。速度场是通过对一个片光源上两次曝光的图像进行互相关处理来得到的。为了研究流体运动的不均匀性对于一个诊断区内速度矢量测量结果的影响，针对五种典型的不均匀性进行讨论。数值实验的结果表明有几种流体运动不均匀性会显著地影响这种粒子图像测速方法的精度。     

DES方法模拟空腔流动及噪声分析

空腔流动中包含多种复杂的物理现象,一直是流体力学中的一个重要问题。本文研究了三维非定常雷诺平均N-S方程和分离涡模拟(Detached eddy simulation,DES)方法在空腔流动及空腔噪音问题的应用。利用两种方法对三维空腔流动及噪音进行了数值计算并与实验数据进行了对比,通过分析比较发现DES方法得到的结果更加准确,特别是通过声压级(Sound pressure level,SPL)分析发现DES方法能够较好地捕捉空腔流动中的压强脉动及噪声水平,说明DES方法能够有效处理空腔流动问题。     

一维非均匀变截面结构瞬态响应的传递函数渐近解法

分析了具有单结构参数变化的一维非均匀变截面结构。引入状态变量将运动微分方程及边界条件写成状态空间形式。定义小参数，通过Laplace变换并利用摄动方法，得到常系数微分方程，从而得到问题的摄动解。通过反Laplace变换得到时域的响应。给出了数值算例，验证了方法的可行性。     

电磁力作用下微型管道中流体流动与混合的实验研究

设计了一种使微管道中导电液体得到混合的方案,利用电场和磁场耦合作用下产生的电磁力作用,使流体产生往复运动及流体界面的弯曲延伸,不同流体的接触面积大大增加,从而提高混合效率.在此基础上搭建了实验平台,利用Micro-PIV系统进行了详细的实验研究,得到了微管道中流体流动的速度场,并对不同电极排列方式下的流场进行分析比较.在实验基础上进行数值验证,并对不同工况下的混合效率进行了分析.结果表明:通过电磁力的扰动作用,确实能有效地提高微管道中流体的混合效率.     

NUMERICAL STUDY OF TWO INTERFACE METHODS FOR COMPRESSIBLE MULTI-FLUID FLOWS

INTRODUCTIONThe direct numerical simulation of com-pressible multi- fluid flows is an important re-search topic with various key applications.Therelatively dominant difficulty of simulation lies intreatment of the moving material interface. Ingeneral,there are two basic approaches to treatcontact discontinuities.The first,front trackingmethod[1,2 ] ,explicitly characterize the interfacemotion by the motion of individual points on theinterface. The other approach is front capturingmethod[3…     

Lattice Boltzmann Flux Solver： An Efficient Approach for Numerical Simulation of Fluid Flows

A lattice Boltzmann flux solver （LBFS） is presented for simulation of fluid flows. Like the conventional computational fluid dynamics （CFD） solvers, the new solver also applies the finite volume method to discretize the governing differential equations, but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers. Instead, it is evaluated from local solution of lattice Boltzmann equation （LBE） at cell interface. Two versions of LBFS are presented in this paper. One is to locally apply one-dimensional compressible lattice Boltzmann （LB） model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves. The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows. The present solver removes the drawbacks of conventional lattice Boltzmann method （LBM） such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.     

GHOST FLUID METHOD FOR RICHTMYER-MESHKOV INSTABILITY SIMULATION

INTRODUCTIONThe energy of the Richtmyer- Meshkov in-stability is excited by the interaction between ashock wave and a contact interface across thatfluids are different.Itis very important thatcap-ture the interface in multi- material flows exact-ly.Here we will apply the ghost fluid method tosimulate two dimensional Richtmyer- Meshkovinstability.The mathematical equation used to simulatethe Ritchmyer- Meshkov instability is the com-pressible Euler equation for two ideal gases.ρρuρvE…     

NUMERICAL SIMULATION ON MULTI-FLUID INTERFACES

INTRODUCTIONThe study on phenomenon of multi- fluid in-terface driven by shock waves include:( 1 ) totrack the motion of multimaterial interfaces;( 2 )to simulate the chemical reaction on the inter-faces. As the pioneering work,chemical reactionis neglected in this paper. The objective is totrack the interfaces. Of course,the fluids aretreated as compressible medium under such ahigh impact.Although Eulerian schemes workwell for most gas flows,they have been shownto be too dissipative nea…     

边界层稳定性的非平行性问题研究

采用抛物化稳定性方程的新方法研究稳定性的非平行问题，其控制方程从涡量形式的N－S方程导出，为提高计算精度和收敛速度，在法向采用高精度的高阶谱方法，结合边界层特征安排边界条件的方程，以及对无限区域数值问题作高效的代数变换处理，通过对预估校正迭代过程中推进步长有效的控制，以满足正规化条件，研究了不同类型压力梯度对稳定性的作用，由系列算例得到的增长率变化曲线和扰动中性曲线等，精确地给出了非平行性对流动稳定性的影响，计算的结果与相关数据符合甚好。     

应用Plemelj积分生成二维通道内的正交网格

首先选用的解析函数由速度模以及流向角组成，在已知通道边界上的流向角时，求得通道边界上的速度。由势函数和流函数的定义，求出壁面上的势函数及流函数。其次将解析函数用势函数表示，运用Ｐｌｅｍｅｌｊ边界积分关系的内点公式，求得通道内任意点的势函数和流函数。由此，从等势线和等流函数值线，可形成流体力学数值计算用的正交网格。文中给出两个算例，该方法无需迭代。     

撞击流混合器内速度脉动的间歇性分析

应用激光多普勒测速仪(LDA),对撞击流混合器的流场速度进行了实验测量.利用小波变换模极大值方法与分形理论结合,对撞击流混合器撞击区的速度脉动信号进行多重分形分析,计算其轴向和径向速度信号的间歇性指数,并与相关流动特性进行关联分析,探讨撞击流混合器撞击区非线性、非均匀性和混沌特性机理.通过分析同一转速不同位置和同一位置不同转速下的信号得出:撞击流混合器内部流体流动具有多重分形的特性,通过间歇性参数随位置变化出现的最低拐点可以判断混合器内部流体运动特征,得出流型的转变区域;转速对撞击区的影响大于对回流区的影响.为揭示撞击流混合器的流动机制提供了一种新的方法.     

改进的CFD/CSD耦合系统设计方法及其计算

论述了界面上的数据交换方法和时间推进格式在流体动力学（CFD）和结构动力学（CSD）耦合求解非线性气动弹性问题中的关键作用;改进了CFD／CSD耦合系统．包括基于边界元方法设计了CFD和CSD耦合界面的数据转换方法,该方法可在同一映射矩阵处理结构响应和非定常气动载荷转换,保证了耦合边界上的能量守恒;改进了一种松耦合方法流程,在保持模块化基础上提高了计算精度和效率。最后将本文方法应用于二维位移外插算例和A—GARD445．6机翼的动响应分析中,结果表明本文方法具有较高精度。     

一类逆风ENO线方法的构造及其对Euler方程的应用

本文利用ENO(essentially-non-oscillatory)插值思想,对非线性双曲型守恒律=0提出了一类逆风ENO线方法,使用这种格式其计算量比C.Shu在文[6]中提出的通量分裂方法明显减少,进而同Runge-kutta时间离散相结合,提出了时空高精度差分格式。本文还把这类方法推广到方程组并用于计算Euler方程组的黎曼问题和拟一维喷管问题,得到了比较满意的数值结果。