用改进的耦合型Level Set方法计算一维双介质可压缩流动

用带有虚拟流体(Ghost Fluid)修正的Level Set方法计算了一维可压缩双介质流动,把描述流动的Euler方程和描述流体界面运动的Level Set方程耦合起来,得到一个整体的守恒律系统,应用高分辨率差分格式求解;为了解决流体界面附近的数值跳动问题,在界面附近引入了虚拟流体方法的Isobaric修正,并给出了算例.     

移动网格与Level Set耦合方法在气-液两相流数值模拟中的应用

本文采用自适应移动网格与Level Set函数相耦合的方法来实现气-液两相流的数值模拟与计算.作为自适应网格方法的一种,移动网格方法主要是为了解决发展方程的计算问题而设计的方法.文中给出了移动网格的生成方程,并针对方程的非线性,给出了一种半隐式的离散方法用于进行求解.本文将移动网格方法与Level Set方法相耦合,将控制流体运动的Navier-Stokes方程以及追踪相界面的Level Set方程转换到曲线坐标下,应用一套曲线坐标方程组来同时描述气、液两相流的运动规律,成功实现了对气-液两相流问题的数值模拟.通过对顶盖驱动流的计算以及对液滴沉降现象的模拟计算,验证了本文方法的可靠性.本文对常重力与微重力下两气泡融合的发展规律进行了数值模拟,通过分析对比,得到了重力对两气泡融合变形的影响规律.     

聚合物充填过程中裹气现象的数值模拟

针对聚合物充填过程中的裹气现象,采用一种有限元（FEM）-间断有限元（DG）耦合算法对其进行数值模拟。对于自由运动界面,采用水平集（Level Set）方法进行捕捉;用XPP（eXtended Pom-Pom）本构模型来描述黏弹性流体的流变行为。采用有限元-间断有限元耦合算法求解统一的流场方程,并采用隐式间断有限元求解XPP本构方程、Level Set及其重新初始化方程。数值结果与文献中的实验结果及模拟结果吻合较好,验证了数值方法的稳定性及准确性。分析带有非规则嵌件型腔内,注射速度及浇口尺寸对裹气现象的影响,裹气容易出现在较高注射速度及较小浇口的情形。     

二维多介质可压缩流的RKDG有限元方法

应用RKDG(Runge-Kutta Discontinuous Galerkin)有限元方法、Level Set方法和Ghost Fluid方法数值模拟二维多介质可压缩流,其中Euler方程组、Level Set方程和重新初始化方程的空间离散采用DG(Discontinuous Galerkin)有限元方法,时间离散采用Runge-Kutta方法.对二维的气-气和气-液两相流进行了数值计算,得到了分辨率较高的计算结果.     

一维多介质可压缩Euler方程的高精度RKDG有限元方法

采用RKDG有限元目的、Level Set目的和改进的带"Isentropic"修正的Ghost Fluid目的模拟了一维多介质可压缩Euler方程,其中Euler方程、Level Set方程和重新初始化方程都采用了三阶精度的RKDG有限元目的进行离散,并对一维两种介质可压缩流体进行了数值实验,得到了较高分辨率的计算结果.     

重新初始化的LS方法跟踪二维可压缩多介质流界面运动

 用非耦合求解方法计算Level Set函数方程与流体力学方程组，应用重新初始化的Level Set函数确保距离函数性质，流体力学方程组采用二阶精度多介质流波传播差分格式计算，重新初始化方程采用五阶WENO格式计算。并给出了二维可压缩多介质流界面运动的计算结果。     

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

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

贴体坐标系中非理想爆轰波阵面传播计算的一个问题

 研究在贴体坐标系中采用Level Set方法计算多维非理想爆轰波阵面传播问题,根据Hamiltom-Jacobi方程的Godunov差分格式,提出了非正交的贴体坐标系中Level Set函数方程的差分格式及其相应的数值方法。最后给出了几个典型数值算例,并且与解析结果进行了比较。     

基于Level Set方法的双层流体热毛细对流的数值研究

基于Level Set方法建立了双层流体热毛细对流的数学模型,通过变密度二阶投影法求解控制方程,C-N隐式技术用于扩散项更新,三阶龙格库塔技术用于对流项的更新,采用连续表面张力模型(CSF)模拟Marangoni效应。三维数值模拟了微重力环境下双层流体系统中交界面变形的热毛细对流,结果显示,在Marangoni效应的作用下,交界面在热端凸起,在冷端凹陷;随着Marangoni数增大,双层流体交界面的变形率随之增大,对流强度也随之增大;交界面与壁面的接触条件会影响热毛细对流的流场和温度场。     

在直角坐标系用NGFM模拟复杂计算域水下高压气泡膨胀问题

采用NGFM(New version of Ghost Fluid Method)处理复杂计算域的固壁边界,用RGFM(Real Ghost Fluid Method)求解气-水界面附近网格节点的状态参数,从而在直角坐标系下对复杂计算域的水下高压气泡膨胀问题进行数值模拟。流场控制方程选用Euler方程,用五阶WENO格式离散空间导数项,二阶Runge-Kutta法离散时间导数项;气-水界面追踪使用Level Set方法,对Level Set方程,用五阶HJ-WENO(Hamilton-Jacobi WENO)和三阶Runge-Kutta法求解。将计算结果与任意坐标系下的结果进行对比,验证了NGFM在笛卡尔网格中处理复杂形状固壁边界的可行性。得到了水下流场压力等值线图、高压气泡的演变过程以及特定点处的压力-时间曲线。计算结果表明,高压气泡在固壁反射激波的作用下,膨胀过程受到抑制;强激波在固壁的反射会导致固壁附近出现大范围的空化流动。     

湿蒸汽两相流动问题的耦合求解

控制湿蒸汽两相流动(汽相和液相)的方程个数多,其数值模拟通常采用将两相分开、两相间用参数传递的方法来处理。本文首次尝试将湿蒸汽两相流动汽相和液相的8个控制方程统一求解,推导了这一方程组的特征值及对应的特征向量,结合控制体积方法,采用TVD格式和时间推进法数值求解这一耦合方程组,研发了相应计算软件。对二维喷管中的自发凝结和非均质凝结算例进行的数值模拟和分析,表明了本文方法和计算软件的有效性和可靠性。     

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.     

非定常粘性空化流动模型及其数值计算

基于气液两相当地均相介质模型,本文给出了一种模拟非定常粘性空化流动的计算模型。认为空化绕流流场中流动介质是一种当地均匀的气液混相物。控制方程采用了应用两相流模型的N-S方程。基于液相和气相的状态方程推导了混合介质密度的表达式。为了保证数值计算的稳定性,控制方程的数值求解采用了TVD-MacCormack格式。为了评价计算模型的可靠性,分别计算了绕台阶和管道水锤的空化流动,所得结果是合理的,说明该方法可以用于空化流动的数值计算。     

A high order kinetic flux-vector splitting method for the reduced five-equation model of compressible two-fluid flows

We present a high order kinetic flux-vector splitting (KFVS) scheme for the numerical solution of a conservative interface-capturing five-equation model of compressible two-fluid flows. This model was initially introduced by Wackers and Koren (2004) [21]. The flow equations are the bulk equations, combined with mass and energy equations for one of the two fluids. The latter equation contains a source term in order to account for the energy exchange. We numerically investigate both one- and two-dimensional flow models. The proposed numerical scheme is based on the direct splitting of macroscopic flux functions of the system of equations. In two space dimensions the scheme is derived in a usual dimensionally split manner. The second order accuracy of the scheme is achieved by using MUSCL-type initial reconstruction and Runge–Kutta time stepping method. For validation, the results of our scheme are compared with those from the high resolution central scheme of Nessyahu and Tadmor [14]. The accuracy, efficiency and simplicity of the KFVS scheme demonstrate its potential for modeling two-phase flows.     

Lattice Boltzmann Simulation of Free-Surface Temperature Dispersion in Shallow Water Flows

We develop a lattice Boltzmann method for modeling free-surface temperature dispersion in the shallow water flows. The governing equations are derived from the incompressible Navier-Stokes equations with assumptions of shallow water flows including bed frictions, eddy viscosity, wind shear stresses and Coriolis forces. The thermal effects are incorporated in the momentum equation by using a Boussinesq approximation. The dispersion of free-surface temperature is modelled by an advection-diffusion equation. Two distribution functions are used in the lattice Boltzmann method to recover the flow and temperature variables using the same lattice structure. Neither upwind discretization procedures nor Riemann problem solvers are needed in discretizing the shallow water equations. In addition, the source terms are straightforwardly included in the model without relying on well-balanced techniques to treat flux gradients and source terms. We validate the model for a class of problems with known analytical solutions and we also present numerical results for sea-surface temperature distribution in the Strait of Gibraltar.     

Hybrid Flux-Splitting Schemes for a Two-Phase Flow Model

In this paper we deal with the construction of hybrid flux-vector-splitting (FVS) schemes and flux-difference-splitting (FDS) schemes for a two-phase model for one-dimensional flow. The model consists of two mass conservation equations (one for each phase) and a common momentum equation. The complexity of this model, as far as numerical computation is concerned, is related to the fact that the flux cannot be expressed in terms of its conservative variables. This is the motivation for studying numerical schemes which are not based on (approximate) Riemann solvers and/or calculations of Jacobian matrix. This work concerns the extension of an FVS type scheme, a Van Leer type scheme, and an advection upstream splitting method (AUSM) type scheme to the current two-phase model. Our schemes are obtained through natural extensions of corresponding schemes studied by Y. Wada and M.-S. Liou (1997, SIAM J. Sci. Comput.18, 633–657) for Euler equations. We explore the various schemes for flow cases which involve both fast and slow transients. In particular, we demonstrate that the FVS scheme is able to capture fast-propagating acoustic waves in a monotone way, while it introduces an excessive numerical dissipation at volume fraction contact (steady and moving) discontinuities. On the other hand, the AUSM scheme gives accurate resolution of contact discontinuities but produces oscillatory approximations of acoustic waves. This motivates us to propose other hybrid FVS/FDS schemes obtained by removing numerical dissipation at contact discontinuities in the FVS and Van Leer schemes.     

旋流排气管的一维非定常流动计算

动力机械装置中广泛存在着非定常旋流流动现象．本文根据质量、动量、能量和旋流动量矩守恒方程，建立了管内非定常旋流流动的一维计算模型，并应用特征线方法推导出了其数值计算格式，是管内非定常一维流动计算的扩展．应用于一台四缸涡轮增压柴油机旋流排气管的计算，通过与实测压力波的比较，表明计算模型有较好的计算精度．     

Explicit and implicit finite difference schemes for fractional Cattaneo equation

In this paper, the numerical solution of fractional (non-integer)-order Cattaneo equation for describing anomalous diffusion has been investigated. Two finite difference schemes namely an explicit predictor–corrector and totally implicit schemes have been developed. In developing each scheme, a separate formulation approach for the governing equations has been considered. The explicit predictor–corrector scheme is the fractional generalization of well-known MacCormack scheme and has been called Generalized MacCormack scheme. This scheme solves two coupled low-order equations and simultaneously computes the flux term with the main variable. Fully implicit scheme however solves a single high-order undecomposed equation. For Generalized MacCormack scheme, stability analysis has been studied through Fourier method. Through a numerical test, the experimental order of convergency of both schemes has been found. Then, the domain of applicability and some numerical properties of each scheme have been discussed.     

高密度和高粘度比率下气液两相流动的数值模拟

本文在VOF方法的基础上,采用粗细两套网格对高密度和高粘度比率下的气液两相流动模拟进行了研究分析.在细网格中求解流体体积函数方程,在粗网格中采用交错网格求解动量方程和压力修正方程,通过粗细网格间的数据传递获得求解动量方程时需要的准确的界面密度和粘度及控制体密度,克服了高密度和高粘度比率下通过插值方法计算界面密度和粘度及控制体密度带来较大误差的困难,保证了质量和动量同时守恒.高密度和高粘度比率下气液两相流动中气液交界面处密度、速度和压力急剧变化,为了保证格式的有界性和稳定性,采用稳定的有界高阶组合格式STOIC.最后模拟了不同工况下气泡在液体中的运动,并通过实验和模拟结果验证了方法的可行性及准确性.     

跨声速流函数方程强隐式解及确定密度场的新方案

基于人工可压缩性对密度的修正,本文用强隐式格式快速求解了非正交曲线坐标系下跨声速流函数方程;在流函数场解出后,通过求解一个由动量方程、能量方程和连续方程组合而成的关于密度的一阶偏微分方程来获得密度场,因此流函数解法中常遇到的密度双值问题在这里已不存在;通常所讲的完成流函数场{Ψ}与密度场{ρ}间的迭代在本文便体现在流函数主方程与这个新推出的一阶微分方程间的迭代计算上;几个典型算例表明了本方法的有效性。