具有H2引燃的CH4，煤油超声速混合的三维数值研究

用双流体模型和分区算法,气相多组分N－S方程和液相Euler方程分别用迎风TVD格式和NND格式进行数值求解,相间相互作用项方程用二阶Runge-Kutta法求解,并首次对H2引燃的CH4和煤油横向喷射,混合问题进行了数值研究,结果表明：与喷射加碘空气的PLIF结果相比,本文和实验结果符合得很好,对直通道,后面喷嘴的穿透深度大,煤油较气态 穿透深度大,但难以被卷吸到喷嘴前的回流区内,流场出现无煤油区,H2和CH4均可扩散到回流区内,CH4和煤油引燃机理不同,对后台阶直通道,存在两个回流区,H2可扩散到这两个回流区,但CH4只能被卷吸到其喷嘴前的回流区,后台阶直通道增强混合的效果和引燃点火的可靠性更好。     

Lindelf的工作与Chaplygin的工作的互补性

运用Vakonomic模型导出Lindel f方程 ,表明Lindel f的工作与Vakonomic模型相吻合 ;运用Chetaev模型导出Chaplygin方程 ,表明Chaplygin的工作与Chetaev模型相吻合·　在此基础上 ,通过改进Chaplygin方程和Lindel f方程的表示形式 ,实现了从Lindel f方程向Chaplygin方程的合理过渡和从Chaplygin方程向Lindel f方程的合理的过渡·　最后 ,给出一个典型实例·　结果表明 ,正如Vako nomic模型与Chetaev模型是互补的一样 ,Lindel f的工作与Chaplygin的工作也是互补的·     

分数阶扩散方程的隐差分格式

根据移位的Grnwald方法,得到求解分数阶扩散方程的三类隐差分格式.利用分数阶von Neumann方法,证明了求解亚扩散方程的两类差分格式是无条件稳定的,而求解超扩散方程的差分格式是条件稳定的,同时也给出了相应差分格式的局部截断误差估计.最后,通过两个数值例子证实了所提出的差分格式的正确性和有效性.     

激波在异种气体中传播及诱导的剪切混合研究

利用二阶迎风TVD格式求解多组分,层流全N－S方程,针对直通道和突扩直通道,研究了马赫数为2和4的激波在H2和空气界面上的传播及诱导的燃料剪切混合,计算结果表明：（1）直通道中,剪切层中的激波阵面要发生畸变,存在对混合起主要作用的卷吸涡,激波马赫数不同,卷吸涡结构和横向混合的尺寸也不同,激波马赫数低,剪切混合效果好,（2）在突扩直通道中,马赫数为2和4的激波在H2中产生不同强度激波,在剪切层中都能产生顺时针,尺度较大的卷吸涡,后台阶增强了剪切层的混合。     

关于一阶剪切板的Kármán型精化理论

依据Mindlin_Reissner理论 ,着重研究一阶剪切板的K rm n型精化理论 ,并推导出仅以挠度和应力函数为未知变量的广义K rm n型大挠度方程 ,适用于复合材料上复合构造剪切板的非线性分析·　在当前板的精化理论中 ,消去的两个转角以隐含形式作用于板的整体变形·　这一工程理论适用于各种计及横向剪切复合材料板、正交各向异性中厚板和夹层板等的线性和非线性分析·　容易发现 ,针对具体的工程应用 ,由该方程可直接获得相应的退化形式 ,并且与文献所载的一致·     

一个求解Euler方程的特殊矩阵分裂格式

§1.引言 自[1]提出矢通量分裂格式以来,在求解气动方程方面得到广泛应用。矢通量分裂格式是一种求解守恒型双曲方程组的方法,它将方程中代表质量、动量和能量的矢通量按照矢通量Jacobian矩阵正负特征值分裂为两个亚矢通量项,目的在于改进显式格式和隐式格式的计算效率和提高求解时的稳定性。在求解方法上,对于二维问题,需要求解以4×4块矩阵为矩阵元的上三角矩阵和下三角矩阵,比中心差分格式需要求解两个块三     

求解对流扩散反应方程的四阶混合紧致差分方法

针对一维对流扩散反应方程,基于对流扩散方程的四阶指数型紧致差分格式,以及一阶导数的四阶Padé公式,发展了一种高效求解对流扩散反应方程的混合型四阶紧致差分格式.数值实验结果验证了格式对于边界层问题或大雷诺数或大Pelect数的大梯度问题的求解的高精度和鲁棒性的优点.     

三维含弥散可压核废料污染问题的交替方向有限元方法

考虑到数值求解三维可压缩核废料污染问题计算量大,利用块有限元逼近技术提出了交替方向有限元格式,将三维问题化为一系列一维问题逐次求解,大大降低了计算量.由于考虑的模型问题为可压缩且同时包含分子扩散和弥散的一般情形,这为误差分析带来困难,本文采用对误差方程进行差商(dt)处理的技巧,证明了格式的最优H1-模误差估计.     

求解一维非定常对流扩散反应方程的高精度紧致差分格式

提出了数值求解一维非定常对流扩散反应方程的一种高精度紧致隐式差分格式,其截断误差为O(τ~4+τ~2h~2+h~4),即格式整体具有四阶精度.差分方程在每一时间层上只用到了三个网格节点,所形成的代数方程组为三对角型,可采用追赶法进行求解,最后通过数值算例验证了格式的精确性和可靠性.     

非定常不可压缩Navier-Stokes方程的旋度形式压力校正格式分析

将求解不可压缩流动的旋度形式压力校正格式从Stokes方程延拓到非定常不可压缩Navier-Stokes方程.在第一步需要求解一个线性的对流-扩散方程,在第二步求解一个Stokes问题.首先给出格式的稳定性分析,然后采用P_2—P_1元分别使用标准形式的压力校正格式和旋度形式的压力校正格式进行了数值模拟,模拟结果表明,对于速度的L~2,H~1误差,两种格式几乎一样,对于压力的L~2误差,旋度形式的压力校正格式略有改进.     

Two-phase laminar axisymmetric jet flow: Explicit,implicit, and split-operator approximations

An hybrid Eulerian-Lagrangian numerical scheme is developed for a two-phase problem and four finite-difference schemes are compared. For this purpose, the problem of hydrodynamic and thermal interactions between a fuel spray and a mixing region of two laminar, unconfined axisymmetric jets is formulated in terms of a set of parabolic differential equations for the gas phase and a set of Lagrangian ordinary differential equations for the condensed phase. Consistent, second-order accurate hybrid numerical schemes, with the exception of the explicit scheme with an accuracy between linear and quadratic, are used to solve these equations. The subset of gas-phase equations has been solved by four different numerical methods: a predictor-corrector explicit method, a sequential implicit method, a block implicit method, and a symmetric operator-splitting method. The subsystem of liquid-phase equations is solved along the droplet trajectories by a second-order Runge-Kutta scheme. The computations have been made to predict the hydro-dynamic and thermal mixing regions of the gas phase as well as the trajectories of each individual group of droplets. In addition, the size, velocity and temperature associated with each group are predicted along these trajectories. The relative merits of the above four difference-schemes are discussed by constructing effectiveness curves. At low error tolerances, the sequential implicit method gives the best results, where for large error tolerances, the explicit and operator splitting give better results. The block implicit scheme is the least effective at all accuracy requirements.     

一类非线性对流扩散方程两重网格特征有限元方法及误差估计

主要研究了一类非线性对流扩散方程的全离散特征有限元方法的两重网格算法及其误差估计.首先在网格步长为H的粗网格上计算一个较小的非线性问题,然后利用一阶牛顿迭代和粗网格解将网格步长为h的细网格上的非线性问题转化为线性问题求解.由于非线性问题的求解仅在粗网格上进行,该两重网格算法可以节省大量的计算工作量,同时具有较高的精度,证明了该两重网格算法L~2模先验误差估计结果为O(△t+h~2+H~(4-d/2)),其中d为空间维数.     

3D numerical simulation of ambient discharge of buoyant water

Waste heat and wastewater are frequently discharged into ambient water and become intermittent sources of buoyancy. In order to control and reduce the environmental impact of these discharges, the mixing characteristics of such discharge in ambient flow should be determined. In this work the transport, mixing and turbulence characteristics of intermittent discharge of buoyant fluid in ambient flow are simulated by a 3D numerical model incorporating a buoyancy extended k–ε model for turbulence. In the numerical model the governing equations are split into three parts in the finite difference solution: advection, dispersion and propagation. The advection part is solved by a characteristics-based scheme. The dispersion part is solved by the central difference method and the propagation part is solved implicitly by using the Gauss–Seidel iteration method. The model has been applied to cases of instantaneous and continuous discharges of buoyancy in ambient water with or without current. Dimensional analysis is used to estimate the initial values. The estimated range of values are found not sensitive to the solution. Satisfactorily comparison between computed results and the experimental results is achieved for the trajectories and lateral widths of the buoyant discharge. The engineering applicability of the model is thus ascertained.     

带有横向射流的三维超声速湍流流场分析

应用显式的五阶WENO格式,结合k-ε湍流模型,求解三维Favre平均N-S方程,计算了从方孔横向喷出的声速气流与马赫数为30的超声速气流的干扰流场。结果表明,在射流上游,射流的阻碍使超声速气流产生分离,形成两个主要的回流区域,主回流导致在方孔射流两侧形成马蹄涡区域,射流下游存在低压区域,形成较小的回流以及一对螺旋形旋涡。     

圆形垂直浮力射流的稳定性与混合特性研究

建立浅水静止环境中圆形轴对称垂直浮力射流的k-ε模型,采用混合有限分析方法进行了数值计算．针对两种不同的流动形态:近区的混合流体以浮力表面层的形式沿径向扩散的稳定排放;近区产生旋涡,浮力热水对混合热水形成二次挟带的非稳定排放,并对稳定性判据进行了验证,最后对两种不同流动结构下的远区的混合特性进行了数值模拟,结果同Lee和Jirka的试验和理论资料均十分吻合．     

Non-isothermal,compressible gas flow for the simulation of an enhanced gas recovery application

In this work, we present a framework for numerical modeling of CO₂ injection into porous media for enhanced gas recovery (EGR) from depleted reservoirs. Physically, we have to deal with non-isothermal, compressible gas flows resulting in a system of coupled non-linear PDEs. We describe the mathematical framework for the underlying balance equations as well as the equations of state for mixing gases. We use an object-oriented finite element method implemented in C++. The numerical model has been tested against an analytical solution for a simplified problem and then applied to CO₂ injection into a real reservoir. Numerical modeling allows to investigate physical phenomena and to predict reservoir pressures as well as temperatures depending on injection scenarios and is therefore a useful tool for applied numerical analysis.     

交错网格紧致差分格式和满足等价性的压力Poisson方程

1．引言随着电子计算机的发展，越来越多的实际问题数值模拟成为现实，但还有很多非线性问题数值计算时间太长，内存要求过大．数值方法的改进可使计算量和存储量大大减少，例如，对二维非定常问题，要使误差达到N－4量级，二阶格式计算点数为（N2）3，（包括时间方向），而四阶格式计算点数仅为N3，差N3倍！而计算量差的倍数更多．当N＝16时N3＝4096，当N＝256时，N31678万．紧致差分格式具有精度高，差分式基点少，＜线性）稳定性好，对高频波分辨率高，边界差分点少等优点【’，‘，’，’。’，’，‘’］，本文中的格式基点数为3，…     

An implicit finite volume method for the solution of 3D low Mach number viscous flows using a local preconditioning technique

This paper presents a cell-centered high order finite volume scheme for the solution of the three-dimensional (3D) Navier–Stokes equations with low Mach number. The system of non-linear equations is solved by means of a fully implicit pseudo-transient scheme. Each pseudo-time step is solved by a Newton-GMRes procedure. A local preconditioning technique is used to scale the speed of sound and to improve the system condition number for low Mach number and low cell Reynolds number. This preconditioning is applied to the AUSM+up flux vector splitting function. The method is tested on 2D and 3D low Mach number laminar flows.