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

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

两介质流界面-激波相互作用RKDG 方法应用分析

为精确模拟多介质流界面运动现象,采用RKDG方法结合虚拟流体方法对气-气、气-液和液-气等多种界面-激波相互作用问题展开研究。数值结果表明,RKDG方法的时空高精度特征使其能够精确、稳健地求解各种复杂界面运动问题。最后,对水下激波自由面折射问题用多种DG格式限制器进行了计算,对比了它们的间断捕捉能力。     

An eigenvector‐based linear reconstruction approach for time stepping in discontinuous Galerkin scheme used to solve shallow water equations

Discontinuous Galerkin (DG) methods have shown promising results for solving the two‐dimensional shallow water equations. In this paper, the classical Runge–Kutta (RK) time discretisation is replaced by the eigenvector‐based reconstruction (EVR) that allows the second‐order time accuracy to be achieved within a single time‐stepping procedure. Moreover, the EVRDG approach yields stable solutions near drying and wetting fronts, whereas the classical RKDG approach yields instabilities. The proposed EVRDG technique is compared with the original RKDG approach on various test cases with analytical solutions. The EVRDG solutions are shown to be as accurate as those obtained with the RKDG scheme. Besides, the EVRDG scheme is 1.6 times faster than the RKDG method. Simulating dambreaks involving dry beds confirms that EVRDG scheme gives correct solutions, whereas the RKDG method yields instabilities. Copyright © 2012 John Wiley & Sons, Ltd.     

Adaptive Runge–Kutta discontinuous Galerkin method for complex geometry problems on Cartesian grid

A Cartesian grid method using immersed boundary technique to simulate the impact of body in fluid has become an important research topic in computational fluid dynamics because of its simplification, automation of grid generation, and accuracy of results. In the frame of Cartesian grid, one often uses finite volume method with second order accuracy or finite difference method. In this paper, an h‐adaptive Runge–Kutta discontinuous Galerkin (RKDG) method on Cartesian grid with ghost cell immersed boundary method for arbitrarily complex geometries is developed. A ghost cell immersed boundary treatment with the modification of normal velocity is presented. The method is validated versus well documented test problems involving both steady and unsteady compressible flows through complex bodies over a wide range of Mach numbers. The numerical results show that the present boundary treatment to some extent reduces the error of entropy and demonstrate the efficiency, robustness, and versatility of the proposed approach. Copyright © 2013 John Wiley & Sons, Ltd.     

用改进的ghost fluid方法模拟强激波与界面的相互作用

为了解决原来的ghost fluid方法在计算强激波和界面相互作用时界面附近出现的速度和压力振荡问题,对原来的ghost fluid方法进行了改进,通过在界面处构造Riemann问题并求出界面的压力和速度,ghost fluid流体的压力和速度分别用界面的压力和速度代替,ghost流体的密度通过熵常数外推得到。改进的ghost fluid保持了原来的ghost fluid的简单性,对一维强激波与气-气、气-液界面的相互作用问题以及射流问题进行了数值计算,得到了分辨率较高的计算结果。     

A practical implementation of high‐order RKDG models for the 1D open‐channel flow equations

This paper comprises an implementation of a fourth‐order Runge–Kutta discontinuous Galerkin (RKDG4) scheme for computing the open‐channel flow equations. The main features of the proposed methodology are simplicity and easiness in the implementation, which may be of possible interest to water resources numerical modellers. A version of the RKDG4 is blended with the Roe Riemann solver, an adaptive high‐order slope limiting procedure, and high‐order source terms approximations. A comparison of the performance of the proposed method with lower‐order RKDG models is performed showing a benefit of considering the RKDG4 model. The scheme is applied to computerize the Saint Venant system by considering benchmark tests that have exact solutions. Finally, numerical results are illustrated discussing the performance of the proposed high‐order model. Copyright © 2008 John Wiley & Sons, Ltd.     

A robust algorithm for moving interface of multi-material fluids based on Riemann solutions

In the paper, the numerical simulation of interface problems for multiple material fluids is studied. The level set function is designed to capture the location of the material interface. For multi-dimensional and multi-material fluids, the modified ghost fluid method needs a Riemann solution to renew the variable states near the interface. Here we present a new convenient and effective algorithm for solving the Riemann problem in the normal direction. The extrapolated variables are populated by Taylor series expansions in the direction. The anti-diffusive high order WENO difference scheme with the limiter is adopted for the numerical simulation. Finally we implement a series of numerical experiments of multi-material flows. The obtained results are satisfying, compared to those by other methods.The English text was polished by Boyi Wang.     

Well‐balancing issues related to the RKDG2 scheme for the shallow water equations

Discontinuous Galerkin (DG) finite element methods have salient features that are mainly highlighted by their locality, their easiness in balancing the flux and source term gradients and their component‐wise structure. In the light of this, this paper aims to provide insights into the well‐balancing property of a second‐order Runge–Kutta Discontinuous Galerkin (RKDG2) method. For this purpose, a Godunov‐type RKDG2 method is presented for solving the shallow water equations. The scheme is based on local DG linear approximations and does not entail any special treatment of the source terms in order to achieve well‐balanced numerical results. The performance of the present RKDG2 scheme in reproducing conserved solutions for both free surface and discharge over strongly irregular topography is demonstrated by applying to several hydraulic benchmarks. Meanwhile, the effects of different slope limiting procedures on the well‐balancing property are investigated and discussed. This work may provide useful guidelines for developing a well‐balanced RKDG2 numerical scheme for shallow water flow simulation. Copyright © 2009 John Wiley & Sons, Ltd.     

Simulations of chemical transport and reaction in a suspension of cells I: an augmented forcing point method for the stationary case

A novel augmented forcing point method is presented to solve the problem of chemical transport in the fluid outside of a collection of suspended cells coupled with chemical reactions on the surfaces of the cells. In this method, the chemical concentrations and the forcing function values are determined simultaneously from an augmented system of equations. The method is more stable and accurate than predictor‐corrector‐type forcing point methods, yields the same solution as a corresponding ghost cell method with much less computational cost, and provides solutions that are pointwise second‐order accurate in space and time, even for closely‐spaced cells. Copyright © 2011 John Wiley & Sons, Ltd.     

Towards numerical simulations of supersonic liquid jets using ghost fluid method

A computational tool based on the ghost fluid method (GFM) is developed to study supersonic liquid jets involving strong shocks and contact discontinuities with high density ratios. The solver utilizes constrained reinitialization method and is capable of switching between the exact and approximate Riemann solvers to increase the robustness. The numerical methodology is validated through several benchmark test problems; these include one-dimensional multiphase shock tube problem, shock–bubble interaction, air cavity collapse in water, and underwater-explosion. A comparison between our results and numerical and experimental observations indicate that the developed solver performs well investigating these problems. The code is then used to simulate the emergence of a supersonic liquid jet into a quiescent gaseous medium, which is the very first time to be studied by a ghost fluid method. The results of simulations are in good agreement with the experimental investigations. Also some of the famous flow characteristics, like the propagation of pressure-waves from the liquid jet interface and dependence of the Mach cone structure on the inlet Mach number, are reproduced numerically. The numerical simulations conducted here suggest that the ghost fluid method is an affordable and reliable scheme to study complicated interfacial evolutions in complex multiphase systems such as supersonic liquid jets.     

求解无粘可压Euler方程组的虚拟流方法

首先将三阶Godunov型半离散中心迎风格式推广到四阶,之后再将该新的四阶半离散中心迎风格式与Level Set方法以及虚拟流方法结合起来,成功地处理了非反应激波问题和多介质流中的爆轰间断问题。由于Level Set函数能隐式地追踪到界面的位置,而虚拟流的构造能隐式地捕捉到界面的边界条件,故而本文的方法可以很自然地推广到多维情况。     

MGFM在强激波与物质界面作用中的应用

原始的虚拟流方法(GFM)在计算强激波和物质界面作用时无法得到正确的计算结果,而改进虚拟流方法(MGFM)处理这类问题的能力大大提高.本文用Level set捕获物质界面,用MGFM方法定义虚拟流节点参数,Euler方程采用HLLC格式离散求解,完成了强激波和物质界面作用的一维和二维数值实验.结果表明改进虚拟流方法在强激波与物质界面作用中的应用是成功的.     

An improvement of classical slope limiters for high‐order discontinuous Galerkin method

In this paper, we describe some existing slope limiters (Cockburn and Shu's slope limiter and Hoteit's slope limiter) for the two‐dimensional Runge–Kutta discontinuous Galerkin (RKDG) method on arbitrary unstructured triangular grids. We describe the strategies for detecting discontinuities and for limiting spurious oscillations near such discontinuities, when solving hyperbolic systems of conservation laws by high‐order discontinuous Galerkin methods. The disadvantage of these slope limiters is that they depend on a positive constant, which is, for specific hydraulic problems, difficult to estimate in order to eliminate oscillations near discontinuities without decreasing the high‐order accuracy of the scheme in the smooth regions. We introduce the idea of a simple modification of Cockburn and Shu's slope limiter to avoid the use of this constant number. This modification consists in: slopes are limited so that the solution at the integration points is in the range spanned by the neighboring solution averages. Numerical results are presented for a nonlinear system: the shallow water equations. Four hydraulic problems of discontinuous solutions of two‐dimensional shallow water are presented. The idealized dam break problem, the oblique hydraulic jump problem, flow in a channel with concave bed and the dam break problem in a converging–diverging channel are solved by using the different slope limiters. Numerical comparisons on unstructured meshes show a superior accuracy with the modified slope limiter. Moreover, it does not require the choice of any constant number for the limiter condition. Copyright © 2008 John Wiley & Sons, Ltd.     

A comparative study of lattice Boltzmann methods using bounce‐back schemes and immersed boundary ones for flow acoustic problems

In order to find applicable treatments of moving boundary conditions based on the lattice Boltzmann method in flow acoustic problems, three bounce‐back (BB) methods and four kinds of immersed boundary (IB) methods are compared. We focused on fluid–solid boundary conditions for flow acoustic problems especially the simulations of sound waves from moving boundaries. BB methods include link bounce‐back, interpolation bounce‐back and unified interpolation bounce‐back methods. Five IB methods are explicit and implicit direct‐forcing (Explicit‐IB and Implicit‐IB), two kinds of partially saturated computational methods and ghost fluid method. In order to reduce the spurious pressure generated by the fresh grid node changing from solid domain to fluid domain for BB methods and sharp IB methods, we proposed two new kinds of treatments and compared them with two existing ones. Simulations of the benchmark problems prove that the local evolutionary iteration (LI) is the best one in treatments of the fresh nodes. In addition, for standing boundary problems, although BB methods have a little higher accuracy, all the methods have similar accuracy. However, for moving boundary problems, IB methods are more appropriate than BB methods, because IB methods' smooth interpolation of pressure eld produces less disturbing spurious pressure waves. With improved treatments of fresh nodes, BB methods are also acceptable for moving boundary acoustic problems. In comparative tests in respective type, unified interpolation bounce‐back with LI, Implicit‐IB, and ghost fluid with LI are the best choices. Copyright © 2013 John Wiley & Sons, Ltd.     

High-resolution algorithms of sharp interface treatment for compressible two-phase flows

In this paper, a kind of arbitrary high order derivatives (ADER) scheme based on the generalised Riemann problem is proposed to simulate multi-material flows by a coupling ghost fluid method. The states at cell interfaces are reconstructed by interpolating polynomials which are piece-wise smooth functions. The states are treated as the equivalent of the left and right states of the Riemann problem. The contact solvers are extrapolated in the vicinity of contact points to facilitate ghost fluids. The numerical method is applied to compressible flows with sharp discontinuities, such as the collision of two fluids of different physical states and gas–liquid two-phase flows. The numerical results demonstrate that unexpected physical oscillations through the contact discontinuities can be prevented effectively and the sharp interface can be captured efficiently.     

三维复杂流动的间断有限元方法模拟

本文基于三维可压缩Euler方程,采用基于Runge-Kutta时间离散的间断有限元方法(RKDG方法),对三维前台阶、三维Riemann问题和球Riemann等问题进行了模拟。结果表明,本文的RKDG方法能够在很少的网格内清晰地捕捉到三维复杂流场中的激波和接触间断;同时,将球Riemann问题中z=0.4平面压强沿到对称轴距离的分布与文献中的近似精确解相比,吻合较好,这也验证了本文的RKDG方法不仅能够进行三维复杂流场的定性描述,也能够应用于三维复杂流场的定量计算。     

SPH方法的运动边界条件研究Research on the moving boundary condition in SPH method

提出一种设置运动边界条件的方法,研究边界附近流体粒子积分截断和非物理穿透边界的问题。边界外的虚粒子在每个时间步由边界附近流体粒子对称生成,赋予相应的物理量,并在虚粒子中引入排斥力,利用拉格朗日形式的N‐S方程自编SPH程序,参考一维激波管的精确解验证边界方法的适用性,研究运动边界条件在计算模型中应用。激波管的模拟结果与精确解基本一致,且在运动边界模型中也计算获得合理的结果。文中提出的运动边界条件,避免了边界附近流体粒子积分截断问题,阻止流体粒子在边界处发生非物理穿透的现象。     

Coupled ghost fluid/two‐phase level set method for curvilinear body‐fitted grids

A coupled ghost fluid/two‐phase level set method to simulate air/water turbulent flow for complex geometries using curvilinear body‐fitted grids is presented. The proposed method is intended to treat ship hydrodynamics problems. The original level set method for moving interface flows was based on Heaviside functions to smooth all fluid properties across the interface. We call this the Heaviside function method (HFM). The HFM requires fine grids across the interface. The ghost fluid method (GFM) has been designed to explicitly enforce the interfacial jump conditions, but the implementation of the jump conditions in curvilinear grids is intricate. To overcome these difficulties a coupled GFM/HFM method was developed in which approximate jump conditions are derived for piezometric pressure and velocity and pressure gradients based on exact continuous velocity and stress and jump in momentum conditions with the jump in density maintained but continuity of the molecular and turbulent viscosities imposed. The implementation of the ghost points is such that no duplication of memory storage is necessary. The level set method is adopted to locate the air/water interface, and a fast marching method was implemented in curvilinear grids to reinitialize the level set function. Validations are performed for three tests: super‐ and sub‐critical flow without wave breaking and an impulsive plunging wave breaking over 2D submerged bumps, and the flow around surface combatant model DTMB 5512. Comparisons are made against experimental data, HFM and single‐phase level set computations. The proposed method performed very well and shows great potential to treat complicated turbulent flows related to ship flows. Copyright © 2007 John Wiley & Sons, Ltd.     

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

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

修正的虚拟介质方法用于水下爆炸冲击薄板的模拟

近期发展的修正的虚拟介质方法(MGFM)已经成功应用于求解多介质流动问题.本文研究了近体水下爆炸中强激波对薄板冲击的演变及薄板的动力响应过程.为了真实地考虑不同介质在界面处的非线性相互作用,本文将MGFM推广应用于处理流体与弹性薄板的相互作用.研究发现尽管流体和薄板结构的模拟基于各自不同的求解方法,但是用修正的虚拟介质...