用Level Set方法追踪运动界面

首先介绍了近年来发展起来的界面追踪技术Level Set方法,然后采用五阶WENO格式和积分平均型TVD格式计算Level Set方程,用修正的Godunov方法求解重新初始化的Level Set方程,数值求解了同心圆在常数流场,圆和矩形界面在剪切流场,缺口圆在旋转流场中的界面变形和还原效果,比较了时间导数离散精度和Level Set函数有无重新初始化对界面追踪效果的影响．最后,通过和其它界面处理方法的比较可以看出,Level Set方法不仅能够比较准确地追踪运动界面,而且无须进行复杂繁琐的界面重构技术,容易编程,具有较大的通用性．     

三维非牛顿熔体前沿界面的Level Set/Ghost/SIMPLEC模拟

建立了三维粘性不可压非牛顿流体流动的控制方程,采用Level Set/Ghost/SIMPLEC方法模拟了注塑成型充模阶段的三维流动过程;追踪到了不同时刻的熔体前沿界面,预测并分析了流动过程中不同时刻的压力、速度等重要的流动特征参数,并与牛顿流体相应的流动特征参数做了对比.研究结果表明:Lovol Set/Ghost/SIMPLEC方法可以准确追踪非牛顿熔体前沿界面:幂律熔体在流动过程中的压差明显大于牛顿熔体的压差,沿横截面的速度分布也有明显的差别.     

Level Set追踪等温非牛顿熔全前沿界面

应用Level Set方法追踪薄壁型腔内Hele-Shaw熔体流动前沿界面,采用5阶加权本质无振荡格式耦合中心差分格式实现了充填阶段的动态模拟.准确追踪到了不同时刻熔体前沿界面,并得到了对应的压力等值线分布,数值结果表明Level Set方法是准确追踪注塑成型熔体前沿界面的一种行之有效的方法.     

注塑成型充模过程的温度场计算

对注塑成型过程中充填阶段的非等温效应的数值模拟进行了研究 ,模拟了非牛顿流体在任意形状型腔内的流动。这里 ,数学模型是基于非牛顿流体在非等温状态下的广义 Hele-Shaw流动 ,利用有限元 /有限差分混合数值方法求解流场中的压力和温度分布     

气液两相流动与固壁相互作用耦合求解的研究

气液两相流动与固壁相互作用的研究是液滴撞击壁面运动研究的重要基础.以结合了VOF和Level Set两种方法优点的用于气液相界面追踪的复合Level Set-VOF方法和利用唯象分析方法建立的能够反映接触角滞后性及壁面性质对润湿过程影响的壁面润湿模型为基础,提出了气液两相流动与固壁相互作用耦合求解流程,给出了气液两相流动与固壁相互作用耦合求解过程中接触线速度的计算方法及边界条件的确定方法.通过与已有实验结果的对比,对提出的气液两相流动与固壁相互作用耦合求解方法的有效性进行了验证.     

Level Set方法在地面效应数值模拟中的应用

利用Level Set方法隐式捕捉界面并能处理复杂的物质界面及其拓扑结构变化（如合并、交叉、破碎等）情况的优势,基于可压缩N-S方程的预处理技术,采用基于界面两侧介质黎曼关系的修正虚拟流体方法对介质边界进行处理;将只能应用于可压缩流动界面模拟的修正虚拟流体方法推广应用于不可压缩流动问题的模拟,数值模拟了气-水运动界面与近水面飞行器的相互影响规律。通过与将水面视为对称面的简化模型对比发现：利用自由运动界面模型计算的NACA6409翼型升力较小,能更准确地反映飞行器流场对水面形状的影响,其结果更为符合实际;文中还对比了翼型距水面不同高度时的结果,得到了翼型受水面影响的强弱与翼型距水面的高度成反比的结论。     

用迭代法模拟注塑成型的三维流动

Hele-Shaw模型是模拟注塑成型过程的常用模型,它的主要缺点是不能模拟一些重要的物理现象及引入实际制件中并不存在的“中面”概念,为了消除这些不足,本文开发了真三维的流动分析程序。建立了粘性、不可压缩的非牛顿流体流动的控制方程,为了避免同时求解耦合的压力场、速度场,本文引入拟稳定技术独立地求解这些方程,用迭代法求耦合方程的解。这种方法可以减少内存并提高数值方法的稳定性。算例表明数值结果与实验结果吻合较好,这种方法成功地模拟了注塑成型流动过程中的重要特征。     

RKDG有限元方法计算流体与刚体耦合

用Level set方法配合Runge-Kutta discontinuous Galerkin (RKDG)有限元方法求解流体与刚体耦合问题。用RKDG有限元方法求解欧拉方程,通过求解Level set方程对界面进行追踪,并用推广的Ghost fluid方法对流刚界面进行处理。数值实验表明,该方法具有较高的分辨率。由于该方法不需要对移动网格进行处理,因此可以处理任意形状的拓扑问题,并且很容易推广到三维。     

基于Level Set的多维双曲守恒律标量方程的激波追踪方法

结合四阶CWENO(Cemral Weighted Essentially Non-Oscillatory)格式、四阶NCE(Natural Continuous Extensions)Runge-Kutta法和Level Set方法，很好地处理了一维双曲守恒律标量方程的激波追踪问题。针对二维双曲守恒律标量方程，成功地用五阶WENO格式、非TVD格式的四阶Runge-Kutta方法和Level Set方法进行激波追踪。将所得的数值解与标准的高阶激波捕捉方法所得的数值解进行比较，说明基于Level Set的激波追踪方法的有效性与逐点收敛性。     

钝感炸药爆轰波传播和驱动的数值模拟研究

为实现对非理想爆轰波传播和爆轰驱动问题的精确计算,研究了与贴体坐标下Level Set方法相结合的爆轰冲击波动力学(DSD)计算方法;根据Hamilton-Jacobi方程的Godunov差分格式,提出了非正交的贴体坐标系中Level Set函数方程的差分格式及其相应的数值方法;并将该DSD方法与二维流体动力学程序TD...     

用Level Set方法求解具有自由面的流动问题

为采用Ｌｅｖｅｌ Ｓｅｔ方法来计算有自由的流动问题提出了一种方案,把自由水面视为水和空气的交界面,两种介质用统一的Ｎ－Ｓ方法求解,在自由面两侧分别采用各自的密度和粘性,并在自由面上给以适当的光滑;采用边界元法求解双调和方程来确定距离函数;Ｎ－Ｓ方程用投影法求解,文中给出了二维水池水面振荡和瞬时溃坝问题的算例,可以看出用ＬｅｖｅｌＳｅｔ方法求解有自由面流动问题是有效的。     

Analysis of a benchmark solution for non-Newtonian radial displacement in porous media

We present an analytical formulation useful to interpret the key phenomena involved in non-Newtonian displacement in porous media and an analysis of the results obtained by considering the uncertainty associated with relevant problem parameters. To derive a benchmark solution, we consider the radial dynamics of a moving stable interface in a porous domain saturated by two fluids, displacing and displaced, both non-Newtonian of shear-thinning power-law behavior, assuming the pressure and velocity to be continuous at the interface, and constant initial pressure. The flow law for both fluids is a modified Darcy’s law. Coupling the nonlinear flow law with the continuity equation, and taking into account compressibility effects, yields a set of nonlinear second-order partial differential equations. Considering two fluids with the same flow behavior index n allows transformation of the latter equations via a self-similar variable; further transformation of the equations incorporating the conditions at the interface shows for n<1 the existence of a compression front ahead of the moving interface. Solving the resulting set of nonlinear equations yields the positions of the moving interface and compression front, and the pressure distributions; the latter are derived in closed form for any value of n. A sensitivity analysis of the model responses is conducted both in a deterministic and a stochastic framework. In the latter case, Global Sensitivity Analysis (GSA) of the benchmark analytical model is adopted to study how the effects of uncertainty affecting selected parameters: (a) the fluids flow behavior index, (b) the relative total compressibility and mobility in the displaced and displacing fluid domains, and (c) the domain permeability and porosity, propagate to state variables. The relative influence of input parameters on model outputs is evaluated by means of associated Sobol indices, calculated via the Polynomial Chaos Expansion (PCE) technique. The goodness of the results obtained by the PCE is assessed by comparison against a traditional Monte Carlo (MC) approach.     

A free surface finite element model for low Froude number mould filling problems on fixed meshes

The simulation of low Froude number mould filling problems on fixed meshes presents significant difficulties. As the Froude number decreases, the coupling between the position of the interface and the resulting flow field increases. The usual two‐phase flow model provides poor results for such flow. In order to overcome the difficulties, a free surface model that applies boundary conditions at the interface accurately is used. Moreover, the use of wall laws on curved boundaries also fails in the case of low Froude number flows. To solve this second problem, we combine wall laws with ‘do nothing’ boundary conditions. A special feature of our approach is that ‘do nothing’ boundary conditions are only applied in the normal direction. These two key ingredients together with the Level Set method allow us to simulate three‐dimensional mould filling problems borrowed directly from the foundry. Copyright © 2010 John Wiley & Sons, Ltd.     

土体颗粒物流动物质点法模拟的弹塑性和非牛顿流体本构模型比较研究

土体颗粒物流动是一种典型的大变形破坏,具有非牛顿流体的流动特征。准确模拟土体颗粒物的流动及冲击过程,对滑坡和泥石流等地质灾害的防治具有重要意义。物质点法是一种无网格粒子类方法,已在各类大变形问题中得到了广泛应用。以往土体颗粒物流动的模拟,通常采用弹塑性本构模型,但缺乏对非牛顿本构模型的模拟分析。本文引入非牛顿本构模型的模拟分析,旨在为土体颗粒物流动模拟提供一种新的方法与思路。非牛顿本构模型的模拟分析是将非牛顿广义Cross模型引入三维物质点法,通过人工阻尼力模拟颗粒间的摩擦力,对土体颗粒物的坍塌、沿斜面滑动以及冲击障碍物等问题进行了动态模拟,研究了其运动全过程,并与弹塑性本构模型的模拟结果进行了对比验证。结果表明,基于非牛顿流体本构模型的物质点法可以较好地模拟土体颗粒物加速、减速到再次稳定的流动全过程及其对障碍物的冲击效应。     

Lattice Boltzmann method for the simulation of viscoelastic fluid flows

The simulation of viscoelastic fluids is a challenging task from the theoretical and numerical points of view. This class of fluids has been extensively studied with the help of classical numerical methods. In this paper we propose a new approach based on the lattice Boltzmann method in order to simulate linear and non-linear viscoelastic fluids and in particular those described by the Oldroyd-B and FENE-P constitutive equations. We study the accuracy and stability of our model on three different benchmarks: the 3D Taylor–Green vortex decay, the simplified 2D four-rolls mill, and the 2D Poiseuille flow. To our knowledge, the methodology described in this work is a first attempt for the simulation of non-trivial flows of viscoelastic fluids using the lattice Boltzmann method to discretize the constitutive and conservation equations.     

Interaction of a solitary wave and a front step simulated by level set method

ＩｎｔｒｏｄｕｃｔｉｏｎＬｅｖｅｌＳｅｔｍｅｔｈｏｄｈａｓｂｅｅｎｄｅｖｅｌｏｐｅｄｆｏｒｃｏｍｐｕｔｉｎｇｔｈｅｍｏｔｉｏｎｏｆｔｗｏ＿ｐｈａｓｅｆｌｏｗｂｙＯｓｈｅｒｅｔａｌ.[1,2 ].Ｉｎｔｈｉｓｗｏｒｋ ,ｔｈｅｉｎｃｏｍｐｒｅｓｓｉｂｌｅＮ＿Ｓｅｑｕａｔｉｏｎｓａｒｅｓｏｌｖｅｄｉｎｔｈｅｗｈｏｌｅｄｏｍａｉｎ ,ａｎｄｔｈｅｉｎｔｅｒｆａｃｅｉｓｔｒａｃｋｅｄｔｈｒｏｕｇｈａＬｅｖｅｌＳｅｔｆｕｎｃｔｉｏｎ .ＴｈｅＬｅｖｅｌＳｅｔｆｕｎｃｔｉｏｎｉｓｄｅｆｉｎｅｄａｓａｄｉｓｔａｎｃｅｆｕｎ…     

Simulation of the two-fluid model on incompressible flow with Fractional Step method for both resolved and unresolved scale interfaces

In the present paper, the Fractional Step method usually used in single fluid flow is here extended and applied for the two-fluid model resolution using the finite volume discretization. The use of a projection method resolution instead of the usual pressure-correction method for multi-fluid flow, successfully avoids iteration processes. On the other hand, the main weakness of the two fluid model used for simulations of free surface flows, which is the numerical diffusion of the interface, is also solved by means of the conservative Level Set method (interface sharpening) (Strubelj et al., 2009). Moreover, the use of the algorithm proposed has allowed presenting different free-surface cases with or without Level Set implementation even under coarse meshes under a wide range of density ratios. Thus, the numerical results presented, numerically verified, experimentally validated and converged under high density ratios, shows the capability and reliability of this resolution method for both mixed and unmixed flows.     

Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau rheological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (ρ₁/ρ_g up to 10³) and high viscosity (η₁/η_g up to 10⁴). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids.     

基于四阶半离散中心迎风格式的虚拟流方法的应用

给出了求解多维无粘可压Euler方程组的四阶半离散中心迎风格式,该格式根据非线性波在网格单元边界上传播的局部速度来更准确地估计局部Riemann的宽度,避免了计算网格的交错,降低了格式的数值粘性。同时,考虑到Level Set函数能隐式地追踪到界面的位置,而虚拟流的构造能隐式地捕捉到界面的边界条件,因此再将新的四阶半离散中心迎风格式与Level Set方法以及虚拟流方法相结合,成功地处理了非反应激波和多介质流中爆轰间断的追踪问题。     

Transient free convection of non-Newtonian fluids along a wavy vertical plate including the magnetic field effect

A Prandtl transformation method is applied to study the transient free convection of non-Newtonian fluids along a wavy vertical plate in the presence of a magnetic field. A simple transformation is proposed to transform the governing equations into the boundary-layer equations and solved numerically by the cubic spline approximation. A simple coordinate transformation is employed to transform the complex wavy surface to a vertical flat plate for a constant wall temperature by the numerical method. The effects of the magnetic field parameter, the wavy geometry and the non-Newtonian nature of the fluids on the flow characteristics and heat transfer are discussed in detail.