Multi‐material closure model for high‐order finite element Lagrangian hydrodynamics

We present a new closure model for single fluid, multi‐material Lagrangian hydrodynamics and its application to high‐order finite element discretizations of these equations 1 . The model is general with respect to the number of materials, dimension and space and time discretizations. Knowledge about exact material interfaces is not required. Material indicator functions are evolved by a closure computation at each quadrature point of mixed cells, which can be viewed as a high‐order variational generalization of the method of Tipton 2 . This computation is defined by the notion of partial non‐instantaneous pressure equilibration, while the full pressure equilibration is achieved by both the closure model and the hydrodynamic motion. Exchange of internal energy between materials is derived through entropy considerations, that is, every material produces positive entropy, and the total entropy production is maximized in compression and minimized in expansion. Results are presented for standard one‐dimensional two‐material problems, followed by two‐dimensional and three‐dimensional multi‐material high‐velocity impact arbitrary Lagrangian–Eulerian calculations. Published 2016. This article is a U.S. Government work and is in the public domain in the USA.     

A compressible single‐temperature conservative two‐phase model with phase transitions

A model for multidimensional compressible two‐phase flow with pressure and velocity relaxations based on the theory of thermodynamically compatible system is extended to study liquid–gas flows with cavitation. The model assumes for each phase its own pressure and velocity, while a common temperature is considered. The governing equations form a hyperbolic system in conservative form and are derived through the theory of a thermodynamically compatible system. The phase pressure‐equalizing process and the interfacial friction are taken into account in the balance laws for the volume fractions of one phase and for the relative velocity by adding two relaxation source terms, while the phase transition is introduced into the model considering in the balance equation for the mass of one phase the relaxation of the Gibbs free energies of the two phases. A modification of the central finite‐volume Kurganov–Noelle–Petrova method is adopted in this work to solve the homogeneous hyperbolic part, while the relaxation source terms are treated implicitly. In order to investigate the effect of the mass transfer in the solution, a 1D cavitation tube problem is presented. In addition, two 2D numerical simulations regarding cavitation problem are also studied: a cavitating Richtmyer–Meshkov instability and a laser‐induced cavitation problem. Copyright © 2014 John Wiley & Sons, Ltd.     

Two‐level finite element method with a stabilizing subgrid for the incompressible Navier–Stokes equations

We consider the Galerkin finite element method for the incompressible Navier–Stokes equations in two dimensions. The domain is discretized into a set of regular triangular elements and the finite‐dimensional spaces employed consist of piecewise continuous linear interpolants enriched with the residual‐free bubble functions. To find the bubble part of the solution, a two‐level finite element method with a stabilizing subgrid of a single node is described, and its application to the Navier–Stokes equation is displayed. Numerical approximations employing the proposed algorithm are presented for three benchmark problems. The results show that the proper choice of the subgrid node is crucial in obtaining stable and accurate numerical approximations consistent with the physical configuration of the problem at a cheap computational cost. Copyright © 2008 John Wiley & Sons, Ltd.     

A higher resolution edge‐based finite volume method for the simulation of the oil–water displacement in heterogeneous and anisotropic porous media using a modified IMPES method

In this article, we present a higher‐order finite volume method with a ‘Modified Implicit Pressure Explicit Saturation’ (MIMPES) formulation to model the 2D incompressible and immiscible two‐phase flow of oil and water in heterogeneous and anisotropic porous media. We used a median‐dual vertex‐centered finite volume method with an edge‐based data structure to discretize both, the elliptic pressure and the hyperbolic saturation equations. In the classical IMPES approach, first, the pressure equation is solved implicitly from an initial saturation distribution; then, the velocity field is computed explicitly from the pressure field, and finally, the saturation equation is solved explicitly. This saturation field is then used to re‐compute the pressure field, and the process follows until the end of the simulation is reached. Because of the explicit solution of the saturation equation, severe time restrictions are imposed on the simulation. In order to circumvent this problem, an edge‐based implementation of the MIMPES method of Hurtado and co‐workers was developed. In the MIMPES approach, the pressure equation is solved, and the velocity field is computed less frequently than the saturation field, using the fact that, usually, the velocity field varies slowly throughout the simulation. The solution of the pressure equation is performed using a modification of Crumpton's two‐step approach, which was designed to handle material discontinuity properly. The saturation equation is solved explicitly using an edge‐based implementation of a modified second‐order monotonic upstream scheme for conservation laws type method. Some examples are presented in order to validate the proposed formulation. Our results match quite well with others found in literature. Copyright © 2016 John Wiley & Sons, Ltd.     

Two‐level finite element method with a stabilizing subgrid for the incompressible MHD equations

We consider the Galerkin finite element method (FEM) for the incompressible magnetohydrodynamic (MHD) equations in two dimension. The domain is discretized into a set of regular triangular elements and the finite‐dimensional spaces employed consist of piecewise continuous linear interpolants enriched with the residual‐free bubble functions. To find the bubble part of the solution, a two‐level FEM with a stabilizing subgrid of a single node is described and its application to the MHD equations is displayed. Numerical approximations employing the proposed algorithm are presented for three benchmark problems including the MHD cavity flow and the MHD flow over a step. The results show that the proper choice of the subgrid node is crucial to get stable and accurate numerical approximations consistent with the physical configuration of the problem at a cheap computational cost. Furthermore, the approximate solutions obtained show the well‐known characteristics of the MHD flow. Copyright © 2009 John Wiley & Sons, Ltd.     

Large eddy simulation (2D) using diffusion–velocity method and vortex‐in‐cell

A large eddy Simulation based on the diffusion‐velocity method and the discrete vortex method is presented. The vorticity‐based and eddy viscosity type subgrid scale model simulating the enstrophy transfer between the large and small scale appears as a convective term in the diffusion‐velocity formulation. The methodology has been tested on a spatially growing mixing layer using the two‐dimensional vortex‐in‐cell method and the Smagorinsky subgrid scale model. The effects on the vorticity contours, momemtum thickness, mean streamwise velocity profiles, root‐mean‐square velocity and vorticity fluctuations and negative cross‐stream correlation are discussed. Comparison is made with experiment and numerical work where diffusion is simulated using random walk. Copyright © 2004 John Wiley & Sons, Ltd.     

Large eddy simulation (2D) of spatially developing mixing layer using vortex‐in‐cell for flow field and filtered probability density function for scalar field

A large eddy simulation based on filtered vorticity transport equation has been coupled with filtered probability density function transport equation for scalar field, to predict the velocity and passive scalar fields. The filtered vorticity transport has been formulated using diffusion‐velocity method and then solved using the vortex method. The methodology has been tested on a spatially growing mixing layer using the two‐dimensional vortex‐in‐cell method in conjunction with both Smagorinsky and dynamic eddy viscosity subgrid scale models for an anisotropic flow. The transport equation for filtered probability density function is solved using the Lagrangian Monte‐Carlo method. The unresolved subgrid scale convective term in filtered density function transport is modelled using the gradient diffusion model. The unresolved subgrid scale mixing term is modelled using the modified Curl model. The effects of subgrid scale models on the vorticity contours, mean streamwise velocity profiles, root‐mean‐square velocity and vorticity fluctuations profiles and negative cross‐stream correlations are discussed. Also the characteristics of the passive scalar, i.e. mean concentration profiles, root‐mean‐square concentration fluctuations profiles and filtered probability density function are presented and compared with previous experimental and numerical works. The sensitivity of the results to the Schmidt number, constant in mixing frequency and inflow boundary conditions are discussed. Copyright © 2005 John Wiley & Sons, Ltd.     

Application of a thermodynamically compatible two‐phase flow model to the high‐resolution simulations of compressible gas–magma flow

This paper reports on the application and development of a fully hyperbolic and fully conservative two‐phase flow model for the simulation of gas and magma flow within volcanic processes. The model solves a set of mixture conservation equations for the gas and magma two‐phase flow with velocity non‐equilibrium. In this model, the effect of the relative velocity is introduced by a kinetic constitutive equation with other equations for volume and mass fractions of the gas phase. The model is examined numerically by the widely used finite volume Godunov methods of centered‐type. Using the Riemann problem, we numerically simulate wave propagation and the development of shocks and rarefactions in volcanic eruptions. These simulations are of magma fragmentation type where the relative velocity continues to dominate. A series of test cases whose solution contains features relevant to gas–magma mixtures are conducted. In particular, numerical results indicate that the model implementation predicts key features of the relative velocity within volcanic processes without any mathematical or physical simplifications. Simulation results are sharply and accurately provided without any spurious oscillations in all of the flow variables. The numerical methods and results are also compared with other numerical methods available in the literature. It is found that the provided resolutions are more accurate for the considered test cases. Copyright © 2014 John Wiley & Sons, Ltd.     

Development of temporal and spatial high‐order schemes for two‐fluid seven‐equation two‐pressure model and its applications in two‐phase flow benchmark problems

Current existing main nuclear thermal‐hydraulics (T‐H) system analysis codes, such as RALAP5, TRACE, and CATHARE, play a crucial role in the nuclear engineering field for the design and safety analysis of nuclear reactor systems. However, two‐fluid model used in these T‐H system analysis codes is ill posed, easily leading to numerical oscillations, and the classical first‐order methods for temporal and special discretization are widely employed for numerical simulations, yielding excessive numerical diffusion. Two‐fluid seven‐equation two‐pressure model is of particular interest due to the inherent well‐posed advantage. Moreover, high‐order accuracy schemes have also attracted great attention to overcome the challenge of serious numerical diffusion induced by low‐order time and space schemes for accurately simulating nuclear T‐H problems. In this paper, the semi‐implicit solution algorithm with high‐order accuracy in space and time is developed for this well‐posed two‐fluid model and the robustness and accuracy are verified and assessed against several important two‐phase flow benchmark tests in the nuclear engineering T‐H field, which include two linear advection problems, the oscillation problem of the liquid column, the Ransom water faucet problem, the reversed water faucet problem, and the two‐phase shock tube problem. The following conclusions are achieved. (1) The proposed semi‐implicit solution algorithm is robust in solving two‐phase flows, even for fast transients and discontinuous solutions. (2) High‐order schemes in both time and space could prevent excessive numerical diffusion effectively and the numerical simulation results are more accurate than those of first‐order time and space schemes, which demonstrates the advantage of using high‐order schemes.     

A Riemann solver and upwind methods for a two‐phase flow model in non‐conservative form

We present a theoretical solution for the Riemann problem for the five‐equation two‐phase non‐conservative model of Saurel and Abgrall. This solution is then utilized in the construction of upwind non‐conservative methods to solve the general initial‐boundary value problem for the two‐phase flow model in non‐conservative form. The basic upwind scheme constructed is the non‐conservative analogue of the Godunov first‐order upwind method. Second‐order methods in space and time are then constructed via the MUSCL and ADER approaches. The methods are systematically assessed via a series of test problems with theoretical solutions. Copyright © 2005 John Wiley & Sons, Ltd.     

Large eddy simulations of wall‐bounded flows using a simplified immersed boundary method and high‐order compact schemes

This paper presents a solution algorithm based on an immersed boundary (IB) method that can be easily implemented in high‐order codes for incompressible flows. The time integration is performed using a predictor‐corrector approach, and the projection method is used for pressure‐velocity coupling. Spatial discretization is based on compact difference schemes and is performed on half‐staggered meshes. A basic algorithm for body‐fitted meshes using the aforementioned solution method was developed by A. Tyliszczak (see article “A high‐order compact difference algorithm for half‐staggered grids for laminar and turbulent incompressible flows” in Journal of Computational Physics) and proved to be very accurate. In this paper, the formulated algorithm is adapted for use with the IB method in the framework of large eddy simulations. The IB method is implemented using its simplified variant without the interpolation (stepwise approach). The computations are performed for a laminar flow around a 2D cylinder, a turbulent flow in a channel with a wavy wall, and around a sphere. Comparisons with literature data confirm that the proposed method can be successfully applied for complex flow problems. The results are verified using the classical approach with body‐fitted meshes and show very good agreement both in laminar and turbulent regimes. The mean (velocity and turbulent kinetic energy profiles and drag coefficients) and time‐dependent (Strouhal number based on the drag coefficient) quantities are analyzed, and they agree well with reference solutions. Two subfilter models are compared, ie, the model of Vreman (see article “An eddy‐viscosity subgrid‐scale model for turbulent shear flow: algebraic theory and applications” in Physics and Fluids) and σ model (Nicoud et al, see article “Using singular values to build a subgrid‐scale model for large eddy simulations” in Physics and Fluids). The tests did not reveal evident advantages of any of these models, and from the point of view of solution accuracy, the quality of the computational meshes turned out to be much more important than the subfilter modeling.     

The numerical solution of the transient two‐phase flow in rigid pipelines

Consideration is given in this paper to the numerical solution of the transient two‐phase flow in rigid pipelines. The governing equations for such flows are two coupled, non‐linear, hyperbolic, partial differential equations with pressure dependent coefficients. The fluid pressure and velocity are considered as two principle dependent variables. The fluid is a homogeneous gas–liquid mixture for which the density is defined by an expression averaging the two‐component densities where a polytropic process of the gaseous phase is admitted. Instead of the void fraction, which varies with the pressure, the gas–fluid mass ratio (or the quality) is assumed to be constant, and is used in the mathematical formulation. The problem has been solved by the method of non‐linear characteristics and the finite difference conservative scheme. To verify their validity, the computed results of the two numerical techniques are compared for different values of the quality, in the case where the liquid compressibility and the pipe wall elasticity are neglected. Copyright © 1999 John Wiley & Sons, Ltd.     

Accurate and consistent particle tracking on unstructured grids

A new numerical method for particle tracking (Lagrangian particle advection) on 2‐D unstructured grids with triangular cells is presented and tested. This method combines key attributes of published methods, including streamline closure for steady flows and local mass conservation (uniformity preservation). The subgrid‐scale velocity reconstruction is linear, and this linear velocity field is integrated analytically to obtain particle trajectories. A complete analytic solution to the 2‐D system of ordinary differential equations (ODEs) governing particle trajectories within a grid cell is provided. The analytic solution to the linear system of locally mass‐conserving constraints that must be enforced to obtain the coefficients in the ODEs is also provided. Numerical experiments are performed to demonstrate that the new method has substantial advantages in accuracy over previously published methods and that it does not suffer from unphysical particle clustering. The method can be used not only in particle‐tracking applications but also as part of a semi‐Lagrangian advection scheme.Copyright © 2015 John Wiley & Sons, Ltd.     

Large eddy simulation using a dynamic mixing length subgrid‐scale model

A novel dynamic mixing length (DML) subgrid‐scale model for large eddy simulations is proposed in this work to improve the cutoff length of the Smagorinsky model. The characteristic mixing length (or the characteristic wave number) is dynamically estimated for the subgrid‐scale fluctuation of turbulence by the cutoff wave‐number, k_c, and the dissipation wave‐number, k_d. The dissipation wave number is derived from the kinetic energy spectrum equation and the dissipation spectrum equation. To prove the promise of the DML model, this model is used to simulate the lid‐driven cubical cavity with max‐velocity‐based Reynolds numbers 8850 and 12,000, the channel flows with friction‐velocity‐based Reynolds numbers 180, 395, 590, and 950, and the turbulent flow past a square cylinder at the higher Reynolds number 21,400, respectively, compared with the Smagorinsky model and Germano et al.'s dynamic Smagorinsky model. Different numerical experiments with different Reynolds numbers show that the DML model can be used in simulations of flows with a wide range of Reynolds numbers without the occurrence of singular values. The DML model can alleviate the dissipation of the Smagorinsky model without the loss of its robustness. The DML model shows some advantages over Germano et al.'s dynamic Smagorinsky model in its high stability and simplicity of calculation because the coefficient of the DML model always stays positive. The characteristic mixing length in the DML model reflects the subgrid‐scale fluctuation of turbulence in nature and thus the characteristic mixing length has a spatial and temporal distribution in turbulent flow. Copyright © 2011 John Wiley & Sons, Ltd.     

Pressure‐stabilized maximum‐entropy methods for incompressible Stokes

We present a parameter‐free stable maximum‐entropy method for incompressible Stokes flow. Derived from a least‐biased optimization inspired by information theory, the meshfree maximum‐entropy method appears as an interesting alternative to classical approximation schemes like the finite element method. Especially compared with other meshfree methods, e.g. the moving least‐squares method, it allows for a straightforward imposition of boundary conditions. However, no Eulerian approach has yet been presented for real incompressible flow, encountering the convective and pressure instabilities. In this paper, we exclusively address the pressure instabilities caused by the mixed velocity‐pressure formulation of incompressible Stokes flow. In a preparatory discussion, existing stable and stabilized methods are investigated and evaluated. This is used to develop different approaches towards a stable maximum‐entropy formulation. We show results for two analytical tests, including a presentation of the convergence behavior. As a typical benchmark problem, results are also shown for the leaky lid‐driven cavity. The already presented information‐flux method for convection‐dominated problems in mind, we see this as the last step towards a maximum‐entropy method capable of simulating full incompressible flow problems. Copyright © 2015 John Wiley & Sons, Ltd.     

An unfitted interior penalty discontinuous Galerkin method for incompressible Navier–Stokes two‐phase flow

A discontinuous Galerkin method for the solution of the immiscible and incompressible two‐phase flow problem based on the nonsymmetric interior penalty method is presented. Therefore, the incompressible Navier–Stokes equation is solved for a domain decomposed into two subdomains with different values of viscosity and density as well as a singular surface tension force. On the basis of a piecewise linear approximation of the interface, meshes for both phases are cut out of a structured mesh. The discontinuous finite elements are defined on the resulting Cartesian cut‐cell mesh and may therefore approximate the discontinuities of the pressure and the velocity derivatives across the interface with high accuracy. As the mesh resolves the interface, regularization of the density and viscosity jumps across the interface is not required. This preserves the local conservation property of the velocity field even in the vicinity of the interface and constitutes a significant advantage compared with standard methods that require regularization of these discontinuities and cannot represent the jumps and kinks in pressure and velocity. A powerful subtessellation algorithm is incorporated to allow the usage of standard time integrators (such as Crank–Nicholson) on the time‐dependent mesh. The presented discretization is applicable to both the two‐dimensional and three‐dimensional cases. The performance of our approach is demonstrated by application to a two‐dimensional benchmark problem, allowing for a thorough comparison with other numerical methods. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical investigation of the aerodynamics of the near‐slot film cooling

Fluid injection from slot or holes into cross‐flow produces highly complicated flow fields. Physical situations encountering the above problem range from turbine blade cooling to waste discharge into rivers. In this paper, the flow field created by a two‐dimensional slot cooling geometry is examined using the finite volume approach with a second‐order upwind differencing scheme. The time‐averaged Navier–Stokes equations were solved on a collocated Cartesian grid with a two‐equation model of turbulence. Attempting to solve the flow field by assuming a uniform velocity profile at the slot exit leads to inaccurate results, while extending the solution domain improves significantly the results, but proves to be costly, both in memory and in computing time (particularly in the case of multiple holes). A pressure‐type boundary condition, based on uniform total pressure, is developed for the slot exit (easily applied to a three‐dimensional geometry), which yields more accurate results than the widely used uniform velocity assumption. It is also found that the implementation of low Reynolds number turbulence models on this geometry provides no significant differences from the standard k–ε model. Copyright © 2000 John Wiley & Sons, Ltd.     

Semi‐implicit subgrid modelling of three‐dimensional free‐surface flows

In this paper, a semi‐implicit numerical model for two‐ and three‐dimensional free‐surface flows will be formulated in such a fashion as to intrinsically account for subgrid bathymetric details. It will be shown that with the proposed subgrid approach the model accuracy can be substantially improved without increasing the corresponding computational effort. Copyright © 2010 John Wiley & Sons, Ltd.     

Implementation of CLEAR algorithm on non‐orthogonal curvilinear co‐ordinates for solution of incompressible flow and heat transfer

In this paper, the CLEAR (coupled and linked equations algorithm revised) algorithm is extended to non‐orthogonal curvilinear collocated grids. The CLEAR algorithm does not introduce pressure correction in order to obtain an incompressible flow field which satisfies the mass conservation law. Rather, it improves the intermediate velocity by solving an improved pressure equation to make the algorithm fully implicit since there is no term omitted in the derivation process. In the extension of CLEAR algorithm from a staggered grid system in Cartesian coordinates to collocated grids in non‐orthogonal curvilinear coordinates, three important issues are appropriately treated so that the extended CLEAR can lead to a unique solution without oscillation of pressure field and with high robustness. These three issues are (1) solution independency on the under‐relaxation factor; (2) strong coupling between velocity and pressure; and (3) treatment of the cross pressure gradient terms. The flow and heat transfer problems in a rectangular enclosure with an internal eccentric circle and the flow in a lid‐driven inclined cavity are computed by using the extended CLEAR. The results show that the extended CLEAR can guarantee the solution independency on the under‐relaxation factor, the smoothness of pressure profile even at very small under‐relaxation factor and good robustness which leads to a converged solution for the small inclined angle of 5^° only with 5‐point computational molecule while the extended SIMPLE‐series algorithm usually can get a converged solution for the inclined angle larger than 30^° under the same condition. Copyright © 2006 John Wiley & Sons, Ltd.     

A solution adaptive simulation of compressible multi‐fluid flows with general equation of state

The unstructured quadrilateral mesh‐based solution adaptive method is proposed in this article for simulation of compressible multi‐fluid flows with a general form of equation of state (EOS). The five equation model (J. Comput. Phys. 2002; 118 :577–616) is employed to describe the compressible multi‐fluid flows. To preserve the oscillation‐free property of velocity and pressure across the interface, the non‐conservative transport equation is discretized in a compatible way of the HLLC scheme for the conservative Euler equations on the unstructured quadrilateral cell‐based adaptive mesh. Five numerical examples, including an interface translation problem, a shock tube problem with two fluids, a solid impact problem, a two‐dimensional Riemann problem and a bubble explosion under free surface, are used to examine its performance in solving the various compressible multi‐fluid flow problems with either the same types of EOS or different types of EOS. The results are compared with those calculated by the following methods: the method with ROE scheme (J. Comput. Phys. 2002; 118 :577–616), the seven equation model (J. Comput. Phys. 1999; 150 :425–467), Shyue's fluid‐mixture model (J. Comput. Phys. 2001; 171 :678–707) or the method in Liu et al. (Comp. Fluids 2001; 30 :315–337). The comparisons for the test problems show that the proposed method seems to be more accurate than the method in Allaire et al. (J. Comput. Phys. 2002; 118 :577–616) or the seven‐equation model (J. Comput. Phys. 1999; 150 :425–467). They also show that it can adaptively and accurately solve these compressible multi‐fluid problems and preserve the oscillation‐free property of pressure and velocity across the material interface. Copyright © 2010 John Wiley & Sons, Ltd.