Ghost Fluid方法与双介质可压缩流动计算

应用带有Isobaric修正的GhostFluid方法配合LevelSet方法计算可压缩双介质无粘流动.该方法可以消除计算流体界面时所产生的数值跳动和耗散,且编程上比界面跟踪法简单.应用WENO格式数值求解欧拉方程和LevelSet方程,对由刚性气体状态方程所支配的一二维双介质流动进行数值计算,得到了分辨率较高的计算结果.     

基于近似Riemann解的有限体积ALE方法

研究二维平面坐标系和二维轴对称坐标系中四边形网格上可压缩流体力学的有限体积ALE(Arbitrary Lagrangian Eulerian)方法.数值方法采用节点中心有限体积法,数值通量采用适用于任意状态方程的HLLC(Harten-Lax-Van Leer-Collela)通量.空间二阶精度通过用WENO(weighted essentially non-oscillatory)方法对原始变量进行重构获得,时间离散采用两步显式Runge-Kutta格式.数值例子显示,方法具有良好的激波分辨能力和高精度的数值逼近能力.     

An Eulerian method for multi-component problems in non-linear elasticity with sliding interfaces

This paper is devoted to developing a multi-material numerical scheme for non-linear elastic solids, with emphasis on the inclusion of interfacial boundary conditions. In particular for colliding solid objects it is desirable to allow large deformations and relative slide, whilst employing fixed grids and maintaining sharp interfaces. Existing schemes utilising interface tracking methods such as volume-of-fluid typically introduce erroneous transport of tangential momentum across material boundaries. Aside from combatting these difficulties one can also make improvements in a numerical scheme for multiple compressible solids by utilising governing models that facilitate application of high-order shock capturing methods developed for hydrodynamics. A numerical scheme that simultaneously allows for sliding boundaries and utilises such high-order shock capturing methods has not yet been demonstrated. A scheme is proposed here that directly addresses these challenges by extending a ghost cell method for gas-dynamics to solid mechanics, by using a first-order model for elastic materials in conservative form. Interface interactions are captured using the solution of a multi-material Riemann problem which is derived in detail. Several different boundary conditions are considered including solid/solid and solid/vacuum contact problems. Interfaces are tracked using level-set functions. The underlying single material numerical method includes a characteristic based Riemann solver and high-order WENO reconstruction. Numerical solutions of example multi-material problems are provided in comparison to exact solutions for the one-dimensional augmented system, and for a two-dimensional friction experiment.     

Van der Waals状态方程可压缩多流体流动的PPM方法

基于流体体积分数的混合型多流体数值模型，将Piecewise Parabolic Method（PPM）方法应用于可压缩多流体流动的数值模拟，采用双波近似求解多流体van der Waals状态方程的Riemann问题．模拟高密度比且含有激波的可压缩多流体流动，典型的纯界面平移问题模拟结果表明，在接触间断的界面附近，压力和速度没有任何的振荡且界面数值耗散都被控制在2—3个网格之内；一维和二维算例表明，该数值方法可以有效地处理接触间断、激波和多维滑移线等物理问题，并能够比其它多流体数值方法更精细地模拟多流体交界面．     

A fully discrete,kinetic energy consistent finite-volume scheme for compressible flows

A robust, implicit, low-dissipation method suitable for LES/DNS of compressible turbulent flows is discussed. The scheme is designed such that the discrete flux of kinetic energy and its rate of change are consistent with those predicted by the momentum and continuity equations. The resulting spatial fluxes are similar to those derived using the so-called skew-symmetric formulation of the convective terms. Enforcing consistency for the time derivative results in a novel density weighted Crank–Nicolson type scheme. The method is stable without the addition of any explicit dissipation terms at very high Reynolds numbers for flows without shocks. Shock capturing is achieved by switching on a dissipative flux term which tends to zero in smooth regions of the flow. Numerical examples include a one-dimensional shock tube problem, the Taylor–Green problem, simulations of isotropic turbulence, hypersonic flow over a double-cone geometry, and compressible turbulent channel flow.     

Simple and efficient relaxation methods for interfaces separating compressible fluids,cavitating flows and shocks in multiphase mixtures

Numerical approximation of the five-equation two-phase flow of Kapila et al. [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10) (2001) 3002–3024] is examined. This model has shown excellent capabilities for the numerical resolution of interfaces separating compressible fluids as well as wave propagation in compressible mixtures [A. Murrone, H. Guillard, A five equation reduced model for compressible two phase flow problems, Journal of Computational Physics 202(2) (2005) 664–698; R. Abgrall, V. Perrier, Asymptotic expansion of a multiscale numerical scheme for compressible multiphase flows, SIAM Journal of Multiscale and Modeling and Simulation (5) (2006) 84–115; F. Petitpas, E. Franquet, R. Saurel, O. Le Metayer, A relaxation-projection method for compressible flows. Part II. The artificial heat exchange for multiphase shocks, Journal of Computational Physics 225(2) (2007) 2214–2248]. However, its numerical approximation poses some serious difficulties. Among them, the non-monotonic behavior of the sound speed causes inaccuracies in wave’s transmission across interfaces. Moreover, volume fraction variation across acoustic waves results in difficulties for the Riemann problem resolution, and in particular for the derivation of approximate solvers. Volume fraction positivity in the presence of shocks or strong expansion waves is another issue resulting in lack of robustness. To circumvent these difficulties, the pressure equilibrium assumption is relaxed and a pressure non-equilibrium model is developed. It results in a single velocity, non-conservative hyperbolic model with two energy equations involving relaxation terms. It fulfills the equation of state and energy conservation on both sides of interfaces and guarantees correct transmission of shocks across them. This formulation considerably simplifies numerical resolution. Following a strategy developed previously for another flow model [R. Saurel, R. Abgrall, A multiphase Godunov method for multifluid and multiphase flows, Journal of Computational Physics 150 (1999) 425–467], the hyperbolic part is first solved without relaxation terms with a simple, fast and robust algorithm, valid for unstructured meshes. Second, stiff relaxation terms are solved with a Newton method that also guarantees positivity and robustness. The algorithm and model are compared to exact solutions of the Euler equations as well as solutions of the five-equation model under extreme flow conditions, for interface computation and cavitating flows involving dynamics appearance of interfaces. In order to deal with correct dynamic of shock waves propagating through multiphase mixtures, the artificial heat exchange method of Petitpas et al. [F. Petitpas, E. Franquet, R. Saurel, O. Le Metayer, A relaxation-projection method for compressible flows. Part II. The artificial heat exchange for multiphase shocks, Journal of Computational Physics 225(2) (2007) 2214–2248] is adapted to the present formulation.     

三角形网格上多介质流动问题界面处理方法

多介质流动问题的求解一般是在结构网格上实现,而三角形网格对于复杂计算区域具有更好的适应性,本文结合rGFM方法,给出三角形网格上多介质流动问题界面处理方法.利用level-set方法跟踪界面,在界面处构造Riemann问题,得到界面处流体准确的流动状态.通过定义界面边界条件,将多介质流动问题转化为单介质流动问题,利用高精度RKDG方法求解.采用多个算例验证该方法的稳健性和有效性,结果表明该方法能准确捕捉界面和激波的位置,保持界面清晰.     

A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows

We propose a new model and a solution method for two-phase compressible flows. The model involves six equations obtained from conservation principles applied to each phase, completed by a seventh equation for the evolution of the volume fraction. This equation is necessary to close the overall system. The model is valid for fluid mixtures, as well as for pure fluids. The system of partial differential equations is hyperbolic. Hyperbolicity is obtained because each phase is considered to be compressible. Two difficulties arise for the solution: one of the equations is written in non-conservative form; non-conservative terms exist in the momentum and energy equations. We propose robust and accurate discretisation of these terms. The method solves the same system at each mesh point with the same algorithm. It allows the simulation of interface problems between pure fluids as well as multiphase mixtures. Several test cases where fluids have compressible behavior are shown as well as some other test problems where one of the phases is incompressible. The method provides reliable results, is able to compute strong shock waves, and deals with complex equations of state.     

Preventing numerical errors generated by interface-capturing schemes in compressible multi-material flows

In the present work, errors generated in computations of compressible multi-material flows using shock-capturing schemes are examined, specifically pressure oscillations (when the specific heats ratio is variable), but also temperature spikes and species conservation errors. These numerical errors are generated at material discontinuities due to an inconsistent treatment of the convective terms. Though temperature errors are irrelevant to solutions to the Euler equations, it is shown that they have the potential to lead to problems when physical diffusion is included, i.e., for the Navier–Stokes equations. These errors are studied analytically and numerically by considering the one-dimensional advection of isolated material discontinuities. A methodology preventing such errors for weighted essentially non-oscillatory (WENO) schemes is presented, in which modified WENO weights are used to solve the transport equation for mass fraction in conservative form to prevent temperature and species conservation errors. Pressure errors are prevented by solving an additional transport equation for a given function of the ratio of specific heats. Several multi-dimensional problems with various discontinuities (shocks, material interfaces and contact discontinuities), including the single-mode Richtmyer–Meshkov instability, and turbulence are considered to test the method.     

Assessment of localized artificial diffusivity scheme for large-eddy simulation of compressible turbulent flows

The localized artificial diffusivity method is investigated in the context of large-eddy simulation of compressible turbulent flows. The performance of different artificial bulk viscosity models are evaluated through detailed results from the evolution of decaying compressible isotropic turbulence with eddy shocklets and supersonic turbulent boundary layer. Effects of subgrid-scale (SGS) models and implicit time-integration scheme/time-step size are also investigated within the framework of the numerical scheme used. The use of a shock sensor along with artificial bulk viscosity significantly improves the scheme for simulating turbulent flows involving shocks while retaining the shock-capturing capability. The proposed combination of Ducros-type sensor with a negative dilatation sensor removes unnecessary bulk viscosity within expansion and weakly compressible turbulence regions without shocks and allows it to localize near the shocks. It also eliminates the need for a wall-damping function for the bulk viscosity while simulating wall-bounded turbulent flows. For the numerical schemes used, better results are obtained without adding an explicit SGS model than with SGS model at moderate Reynolds number. Inclusion of a SGS model in addition to the low-pass filtering and artificial bulk viscosity results in additional damping of the resolved turbulence. However, investigations at higher Reynolds numbers suggest the need for an explicit SGS model. The flow statistics obtained using the second-order implicit time-integration scheme with three sub-iterations closely agrees with the explicit scheme if the maximum Courant–Friedrichs–Lewy is kept near unity.     

A shock-capturing methodology based on adaptative spatial filtering for high-order non-linear computations

A shock-capturing methodology is developed for non-linear computations using low-dissipation schemes and centered finite differences. It consists in applying an adaptative second-order filtering to handle discontinuities in combination with a background selective filtering to remove grid-to-grid oscillations. The shock-capturing filtering is written in its conservative form, and its magnitude is determined dynamically from the flow solutions. A shock-detection procedure based on a Jameson-like shock sensor is derived so as to apply the shock-capturing filtering only around shocks. A second-order filter with reduced errors in the Fourier space with respect to the standard second-order filter is also designed. Linear and non-linear 1D and 2D problems are solved to show that the methodology is capable of capturing shocks without providing dissipation outside shocks. The shock detection allows in particular to distinguish shocks from linear waves, and from vortices when it is performed from dilatation rather than from pressure. Finally the methodology is simple to implement and reasonable in terms of computational cost.     

FIVER: A finite volume method based on exact two-phase Riemann problems and sparse grids for multi-material flows with large density jumps

A robust finite volume method for the solution of high-speed compressible flows in multi-material domains involving arbitrary equations of state and large density jumps is presented. The global domain of interest can include a moving or deformable subdomain that furthermore may undergo topological changes due to, for example, crack propagation. The key components of the proposed method include: (a) the definition of a discrete surrogate material interface, (b) the computation of a reliable approximation of the fluid state vector on each side of a discrete material interface via the construction and solution of a local, exact, two-phase Riemann problem, (c) the algebraic solution of this auxiliary problem when the equation of state allows it, and (d) the solution of this two-phase Riemann problem using sparse grid tabulations otherwise. The proposed computational method is illustrated with the three-dimensional simulation of the dynamics of an underwater explosion bubble.     

Assessment of high-resolution methods for numerical simulations of compressible turbulence with shock waves

Flows in which shock waves and turbulence are present and interact dynamically occur in a wide range of applications, including inertial confinement fusion, supernovae explosion, and scramjet propulsion. Accurate simulations of such problems are challenging because of the contradictory requirements of numerical methods used to simulate turbulence, which must minimize any numerical dissipation that would otherwise overwhelm the small scales, and shock-capturing schemes, which introduce numerical dissipation to stabilize the solution. The objective of the present work is to evaluate the performance of several numerical methods capable of simultaneously handling turbulence and shock waves. A comprehensive range of high-resolution methods (WENO, hybrid WENO/central difference, artificial diffusivity, adaptive characteristic-based filter, and shock fitting) and suite of test cases (Taylor–Green vortex, Shu–Osher problem, shock-vorticity/entropy wave interaction, Noh problem, compressible isotropic turbulence) relevant to problems with shocks and turbulence are considered. The results indicate that the WENO methods provide sharp shock profiles, but overwhelm the physical dissipation. The hybrid method is minimally dissipative and leads to sharp shocks and well-resolved broadband turbulence, but relies on an appropriate shock sensor. Artificial diffusivity methods in which the artificial bulk viscosity is based on the magnitude of the strain-rate tensor resolve vortical structures well but damp dilatational modes in compressible turbulence; dilatation-based artificial bulk viscosity methods significantly improve this behavior. For well-defined shocks, the shock fitting approach yields good results.     

A low numerical dissipation immersed interface method for the compressible Navier–Stokes equations

A numerical method to solve the compressible Navier–Stokes equations around objects of arbitrary shape using Cartesian grids is described. The approach considered here uses an embedded geometry representation of the objects and approximate the governing equations with a low numerical dissipation centered finite-difference discretization. The method is suitable for compressible flows without shocks and can be classified as an immersed interface method. The objects are sharply captured by the Cartesian mesh by appropriately adapting the discretization stencils around the irregular grid nodes, located around the boundary. In contrast with available methods, no jump conditions are used or explicitly derived from the boundary conditions, although a number of elements are adopted from previous immersed interface approaches. A new element in the present approach is the use of the summation-by-parts formalism to develop stable non-stiff first-order derivative approximations at the irregular grid points. Second-order derivative approximations, as those appearing in the transport terms, can be stiff when irregular grid points are located too close to the boundary. This is addressed using a semi-implicit time integration method. Moreover, it is shown that the resulting implicit equations can be solved explicitly in the case of constant transport properties. Convergence studies are performed for a rotating cylinder and vortex shedding behind objects of varying shapes at different Mach and Reynolds numbers.     

Relativistic Hydrodynamics and Essentially Non-oscillatory Shock Capturing Schemes

Numerical solutions of relativistic hydrodynamic equations are obtained with essentially non-oscillatory (ENO) finite differencing schemes. The method is explicit, conservative, consistent with the entropy condition, and high order accurate in space and time. The present implementation is applicable to the most general, three-dimensional problems with an arbitrary equation of state. Numerical experiments, including computations of multi-dimensional flows, demonstrate that the method delivers sharp, non-oscillatory shock transitions without sacrificing high resolution of the smooth regions. This extends results already established for the Euler gas dynamics to the relativistic regime, suggesting the usefulness of ENO schemes for modelling relativistic nuclear collisions.     

Two-phase flow model for energetic proton beam induced pressure waves in mercury target systems in the planned European Spallation Source

Two-phase flow calculations are presented to investigate the thermo-hydraulical effects of the interaction between 2 ms long 1.3 GeV proton pulses with a closed mercury loop which can be considered as a model system of the target of the planned European Spallation Source (ESS) facility. The two-fluid model consists of six first-order partial differential equations that present one dimensional mass, momentum and energy balances for mercury vapor and liquid phases are capable to describe quick transients like cavitation effects or shock waves. The absorption of the proton beam is represented as instantaneous heat source in the energy balance equations. Densities and internal energies of the mercury liquid-vapor system is calculated from the van der Waals equation, but general method how to obtain such properties using arbitrary equation of state is also presented. A second order accurate high-resolution shock-capturing numerical scheme is applied with different kind of limiters in the numerical calculations. Our analysis show that even 75 degree temperature heat shocks cannot cause considerable cavitation effects in mercury.     

A continuum model of gas flows with localized density variations

We discuss the kinetic representation of gases and the derivation of macroscopic equations governing the thermomechanical behavior of a dilute gas viewed at the macroscopic level as a continuous medium. We introduce an approach to kinetic theory where spatial distributions of the molecules are incorporated through a mean-free-volume argument. The new kinetic equation derived contains an extra term involving the evolution of this volume, which we attribute to changes in the thermodynamic properties of the medium. Our kinetic equation leads to a macroscopic set of continuum equations in which the gradients of thermodynamic properties, in particular density gradients, impact on diffusive fluxes. New transport terms bearing both convective and diffusive natures arise and are interpreted as purely macroscopic expansion or compression. Our new model is useful for describing gas flows that display non-local-thermodynamic-equilibrium (rarefied gas flows), flows with relatively large variations of macroscopic properties, and/or highly compressible fluid flows.     

Comparison of the generalized Riemann solver and the gas-kinetic scheme for inviscid compressible flow simulations

The generalized Riemann problem (GRP) scheme for the Euler equations and gas-kinetic scheme (GKS) for the Boltzmann equation are two high resolution shock capturing schemes for fluid simulations. The difference is that one is based on the characteristics of the inviscid Euler equations and their wave interactions, and the other is based on the particle transport and collisions. The similarity between them is that both methods can use identical MUSCL-type initial reconstructions around a cell interface, and the spatial slopes on both sides of a cell interface involve in the gas evolution process and the construction of a time-dependent flux function. Although both methods have been applied successfully to the inviscid compressible flow computations, their performances have never been compared. Since both methods use the same initial reconstruction, any difference is solely coming from different underlying mechanism in their flux evaluation. Therefore, such a comparison is important to help us to understand the correspondence between physical modeling and numerical performances. Since GRP is so faithfully solving the inviscid Euler equations, the comparison can be also used to show the validity of solving the Euler equations itself. The numerical comparison shows that the GRP exhibits a slightly better computational efficiency, and has comparable accuracy with GKS for the Euler solutions in 1D case, but the GKS is more robust than GRP. For the 2D high Mach number flow simulations, the GKS is absent from the shock instability and converges to the steady state solutions faster than the GRP. The GRP has carbuncle phenomena, likes a cloud hanging over exact Riemann solvers. The GRP and GKS use different physical processes to describe the flow motion starting from a discontinuity. One is based on the assumption of equilibrium state with infinite number of particle collisions, and the other starts from the non-equilibrium free transport process to evolve into an equilibrium one through particle collisions. The different mechanism in the flux evaluation deviates their numerical performance. Through this study, we may conclude scientifically that it may NOT be valid to use the Euler equations as governing equations to construct numerical fluxes in a discretized space with limited cell resolution. To adapt the Navier–Stokes (NS) equations is NOT valid either because the NS equations describe the flow behavior on the hydrodynamic scale and have no any corresponding physics starting from a discontinuity. This fact alludes to the consistency of the Euler and Navier–Stokes equations with the continuum assumption and the necessity of a direct modeling of the physical process in the discretized space in the construction of numerical scheme when modeling very high Mach number flows. The development of numerical algorithm is similar to the modeling process in deriving the governing equations, but the control volume here cannot be shrunk to zero.     

High-order shock-fitting methods for direct numerical simulation of hypersonic flow with chemical and thermal nonequilibrium

In recent years, much progress has been made in the direct numerical simulation of laminar-turbulent transition of hypersonic boundary layer flow. However, most of the efforts at the direct numerical simulation of transition previously have been focused on the idealized perfect gas flow or “cold” hypersonic flows. For practical problems in hypersonic flows, high-temperature effects of thermal and chemical nonequilibrium are important and cannot be modeled by a perfect gas model. Therefore, it is necessary to include the real gas models in the numerical simulation of hypersonic boundary layer transition in order to accurately predict flow field parameters. Currently most numerical methods for hypersonic flow with thermo-chemical nonequilibrium are based on shock-capturing approach at relatively low order of accuracy. Shock capturing schemes reduce to first-order accuracy near the shock and have been shown to produce spurious oscillations behind curved strong shocks. There is a need to develop new methods capable of simulating nonequilibrium hypersonic flow fields with uniformly high-order accuracy and avoid spurious oscillations near the shock. This paper presents a fifth-order shock-fitting method for numerical simulation of thermal and chemical nonequilibrium in hypersonic flows. The method is developed based on the state-of-the-art real gas models for thermo-chemical nonequilibrium and transport phenomena. Shock-fitting approach is used because it has the advantage of capturing the entire flow field with high-order accuracy and without any oscillations near the shock. The new method has been tested and validated for a number of test cases over a wide span of free stream conditions. The developed method is applied for the study of receptivity of free stream acoustic waves over a blunt cone for hypervelocity flow. Some preliminary results of the computations of the high order shock fitting method for the above mentioned study have also been presented.     

非结构网格下非定常流场双时间步长的加权ENO-强紧致杂交高分辨率格式

针对三维非定常、可压缩流场的Navier-Stokes方程组,本文提出一种新的双时间步长高精度快速迭代格式。该格式在时间上具有二阶精度,在空间离散上不低于三阶。在对流项与粘性项的处理上,本格式分别采用了加权ENO-强紧致格式与紧致四阶精度格式的思想。几个典型算例的实践表明:计算结果与相关实验数据比较吻合,初步表明了该算法可以在非结构网格下具有高效率与高分辨率的特征。