Characteristic boundary conditions for simulations of compressible reacting flows with multi-dimensional,viscous and reaction effects

A generalized formulation of the characteristic boundary conditions for compressible reacting flows is proposed. The new and improved approach resolves a number of lingering issues of spurious solution behaviour encountered in turbulent reacting flow simulations in the past. This is accomplished (a) by accounting for all the relevant terms in the determination of the characteristic wave amplitudes and (b) by accommodating a relaxation treatment for the transverse gradient terms with the relaxation coefficient properly determined by the low Mach number asymptotic expansion. The new boundary conditions are applied to a comprehensive set of test problems including: vortex-convection; turbulent inflow; ignition front propagation; non-reacting and reacting Poiseuille flows; and counterflow cases. It is demonstrated that the improved boundary conditions perform consistently superior to existing approaches, and result in robust and accurate solutions with minimal acoustic wave interactions at the boundary in hostile turbulent combustion simulation conditions.     

An adaptive version of ghost-cell immersed boundary method for incompressible flows with complex stationary and moving boundaries

An adaptive version of immersed boundary method for simulating flows with complex stationary and moving boundaries is presented.The method employs a ghost-cell methodology which allows for a sharp representation of the immersed boundary.To simplify the implementation of the methodology,a volume-of-fluid method is introduced to identify the immersed boundary.In addition,the domain is spatially discretized using a tree-based discretization which is relatively simple to implement a fully flexible adaptive refi...     

A discontinuous Galerkin method for inviscid low Mach number flows

In this work we extend the high-order discontinuous Galerkin (DG) finite element method to inviscid low Mach number flows. The method here presented is designed to improve the accuracy and efficiency of the solution at low Mach numbers using both explicit and implicit schemes for the temporal discretization of the compressible Euler equations. The algorithm is based on a classical preconditioning technique that in general entails modifying both the instationary term of the governing equations and the dissipative term of the numerical flux function (full preconditioning approach). In the paper we show that full preconditioning is beneficial for explicit time integration while the implicit scheme turns out to be efficient and accurate using just the modified numerical flux function. Thus the implicit scheme could also be used for time accurate computations. The performance of the method is demonstrated by solving an inviscid flow past a NACA0012 airfoil at different low Mach numbers using various degrees of polynomial approximations. Computations with and without preconditioning are performed on different grid topologies to analyze the influence of the spatial discretization on the accuracy of the DG solutions at low Mach numbers.     

Analysis of an immersed boundary method for three-dimensional flows in vorticity formulation

This article presents numerical analysis and practical considerations for three-dimensional flow computation using an implicit immersed boundary method. The Euler equations, or half a step of the Navier–Stokes equations when using fractional step algorithms, are investigated in their vorticity formulation. The context of flow computation around an arbitrarily shaped body is especially investigated.In conventional immersed boundary methods using vorticity, singular vortex are dispatched over the body surface. In the present study, one prefers using sources of potential velocity field, dispatched on the body, whose nature is not vorticity. Such a formulation is compatible to the Euler equations. In practice, these sources of potential flow produce a velocity through this surface, aiming in practice at cancelling a flow-through velocity.This article focuses on the use of the source-to-flow-through linear application, its properties being the key points for fast convergence. Its self-adjointness, or lack thereof, conditioning and preconditioning aspects are investigated. It follows that computing a velocity field with no-flow-through conditions in complex geometry, when using the source-to-flow-through linear application, can be achieved for 4/3 of the computational cost of standard Poisson equation in a Cartesian box.The robustness of immersed boundaries is especially interesting when used together with vortex-in-cell methods, well known for their robustness in time and their ability to compute accurately convective effects. A few examples, based on real-world geometries, illustrate the method capabilities.     

A front-tracking/ghost-fluid method for fluid interfaces in compressible flows

A front-tracking/ghost-fluid method is introduced for simulations of fluid interfaces in compressible flows. The new method captures fluid interfaces using explicit front-tracking and defines interface conditions with the ghost-fluid method. Several examples of multiphase flow simulations, including a shock–bubble interaction, the Richtmyer–Meshkov instability, the Rayleigh–Taylor instability, the collapse of an air bubble in water and the breakup of a water drop in air, using the Euler or the Navier–Stokes equations, are performed in order to demonstrate the accuracy and capability of the new method. The computational results are compared with experiments and earlier computational studies. The results show that the new method can simulate interface dynamics accurately, including the effect of surface tension. Results for compressible gas–water systems show that the new method can be used for simulations of fluid interface with large density differences.     

A Brinkman penalization method for compressible flows in complex geometries

To simulate flows around solid obstacles of complex geometries, various immersed boundary methods had been developed. Their main advantage is the efficient implementation for stationary or moving solid boundaries of arbitrary complexity on fixed non-body conformal Cartesian grids. The Brinkman penalization method was proposed for incompressible viscous flows by penalizing the momentum equations. Its main idea is to model solid obstacles as porous media with porosity, , and viscous permeability approaching zero. It has the pronounced advantages of mathematical proof of error bound, strong convergence, and ease of numerical implementation with the volume penalization technique. In this paper, it is extended to compressible flows. The straightforward extension of penalizing momentum and energy equations using Brinkman penalization with respective normalized viscous, η, and thermal, η_T, permeabilities produces unsatisfactory results, mostly due to nonphysical wave transmissions into obstacles, resulting in considerable energy and mass losses in reflected waves. The objective of this paper is to extend the Brinkman penalization technique to compressible flows based on a physically sound mathematical model for compressible flows through porous media. In addition to penalizing momentum and energy equations, the continuity equation for porous media is considered inside obstacles. In this model, the penalized porous region acts as a high impedance medium, resulting in negligible wave transmissions. The asymptotic analysis reveals that the proposed Brinkman penalization technique results in the amplitude and phase errors of order O((η)^1/2) and O((η/η_T)^1/43/4), when the boundary layer within the porous media is respectively resolved or unresolved. The proposed method is tested using 1- and 2-D benchmark problems. The results of direct numerical simulation are in excellent agreement with the analytical solutions. The numerical simulations verify the accuracy and convergence rates.     

A Monte Carlo model for seeded atomic flows in the transition regime

A simple model for the numerical determination of separation effects in seeded atomic gas flows is presented. The model is based on the known possibility to provide a statistically convergent estimate of the exact solution for a linear transport equation using the test particle Monte Carlo method. Accordingly, the flow field of the main gas is preliminary calculated and as a second step the linear transport equations obtained by fixing the target distribution in the collision term of the Boltzmann equation for both main and minority components are solved. Both solutions are based on appropriately devised test particle Monte Carlo methods. The second step, the critical one in evaluating the separation effects, is exact and thereby completely free of numerical diffusion. The model is described in details and illustrated by 2D test cases of atomic separation in shock fronts.     

Characteristic boundary conditions for direct simulations of turbulent counterflow flames

Improved Navier–Stokes characteristic boundary conditions (NSCBC) are formulated for the direct numerical simulations (DNS) of laminar and turbulent counterflow flame configurations with a compressible flow formulation. The new boundary scheme properly accounts for multi-dimensional flow effects and provides nonreflecting inflow and outflow conditions that maintain the mean imposed velocity and pressure, while substantially eliminating spurious acoustic wave reflections. Applications to various counterflow configurations demonstrate that the proposed boundary conditions yield accurate and robust solutions over a wide range of flow and scalar variables, allowing high fidelity in detailed numerical studies of turbulent counterflow flames.     

Large eddy simulation using a new set of sixth order schemes for compressible viscous terms

A set of conservative sixth order central differencing schemes is suggested for compressible flows with variable viscosity coefficient. This new set of central differencing schemes has the stencil width matching that of the seventh order weighted essentially non-oscillatory scheme (WENO). This feature is important to maintain the compactness of the seventh order WENO scheme and facilitate boundary condition treatment. As an application example, a large eddy simulation (LES) is conducted for a cylinder flow using the seventh order WENO scheme for the convective terms and the new set of sixth order central differencing scheme for the viscous terms. The results are compared with those from other research groups and those obtained using the fifth order WENO scheme and fourth order central differencing.     

Implementation of the GRP scheme for computing radially symmetric compressible fluid flows

The study of radially symmetric compressible fluid flows is interesting both from the theoretical and numerical points of view. Spherical explosion and implosion in air, water and other media are well-known problems in application. Typical difficulties lie in the treatment of singularity in the geometrical source and the imposition of boundary conditions at the symmetric center, in addition to the resolution of classical discontinuities (shocks and contact discontinuities). In the present paper we present the implementation of direct generalized Riemann problem (GRP) scheme to resolve this issue. The scheme is obtained directly by the time integration of the fluid flows. Our new contribution is to show rigorously that the singularity is removable and derive the updating formulae for mass and energy at the center. Together with the vanishing of the momentum, we obtain new numerical boundary conditions at the center, which are then incorporated into the GRP scheme. The main ingredient is the passage from the Cartesian coordinates to the radially symmetric coordinates.     

A massively parallel multi-block hybrid compact–WENO scheme for compressible flows

A new multi-block hybrid compact–WENO finite-difference method for the massively parallel computation of compressible flows is presented. In contrast to earlier methods, our approach breaks the global dependence of compact methods by using explicit finite-difference methods at block interfaces and is fully conservative. The resulting method is fifth- and sixth-order accurate for the convective and diffusive fluxes, respectively. The impact of the explicit interface treatment on the stability and accuracy of the multi-block method is quantified for the advection and diffusion equations. Numerical errors increase slightly as the number of blocks is increased. It is also found that the maximum allowable time steps increase with the number of blocks. The method demonstrates excellent scalability on up to 1264 processors.     

A high order ENO conservative Lagrangian type scheme for the compressible Euler equations

We develop a class of Lagrangian type schemes for solving the Euler equations of compressible gas dynamics both in the Cartesian and in the cylindrical coordinates. The schemes are based on high order essentially non-oscillatory (ENO) reconstruction. They are conservative for the density, momentum and total energy, can maintain formal high order accuracy both in space and time and can achieve at least uniformly second-order accuracy with moving and distorted Lagrangian meshes, are essentially non-oscillatory, and have no parameters to be tuned for individual test cases. One and two-dimensional numerical examples in the Cartesian and cylindrical coordinates are presented to demonstrate the performance of the schemes in terms of accuracy, resolution for discontinuities, and non-oscillatory properties.     

A moving Kriging interpolation-based boundary node method for two-dimensional potential problems

In this paper,a meshfree boundary integral equation(BIE) method,called the moving Kriging interpolationbased boundary node method(MKIBNM),is developed for solving two-dimensional potential problems.This study combines the BIE method with the moving Kriging interpolation to present a boundary-type meshfree method,and the corresponding formulae of the MKIBNM are derived.In the present method,the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker’s delta property,then the boundary conditions can be imposed directly and easily.To verify the accuracy and stability of the present formulation,three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.     

Analytic solutions of simple flows and analysis of nonslip boundary conditions for the lattice Boltzmann BGK model

In this paper we analytically solve the velocity of the lattice Boltzmann BGK equation (LBGK) for several simple flows. The analysis provides a framework to theoretically analyze various boundary conditions. In particular, the analysis is used to derive the slip velocities generated by various schemes for the nonslip boundary condition. We find that the slip velocity is zero as long as fe=0 at boundaries, no matter what combination of distributions is chosen. The schemes proposed by Nobleet al. and by Inamuroet al. yield the correct zeroslip velocity, while some other schemes, such as the bounce-back scheme and the equilibrium distribution scheme, would inevitably generate a nonzero slip velocity. The bounce-back scheme with the wall located halfway between a flow node and a bounce-back node is also studied for the simple flows considered and is shown to produce results of second-order accuracy. The momentum exchange at boundaries seems to be highly related to the slip velocity at boundaries. To be specific, the slip velocity is zero only when the momentum dissipated by boundaries is equal to the stress provided by fluids.     

A numerical method for simulating the dynamics of 3D axisymmetric vesicles suspended in viscous flows

We extend [Shravan K. Veerapaneni, Denis Gueyffier, Denis Zorin, George Biros, A boundary integral method for simulating the dynamics of inextensible vesicles suspended in a viscous fluid in 2D, Journal of Computational Physics 228(7) (2009) 2334–2353] to the case of three-dimensional axisymmetric vesicles of spherical or toroidal topology immersed in viscous flows. Although the main components of the algorithm are similar in spirit to the 2D case—spectral approximation in space, semi-implicit time-stepping scheme—the main differences are that the bending and viscous force require new analysis, the linearization for the semi-implicit schemes must be rederived, a fully implicit scheme must be used for the toroidal topology to eliminate a CFL-type restriction and a novel numerical scheme for the evaluation of the 3D Stokes single layer potential on an axisymmetric surface is necessary to speed up the calculations. By introducing these novel components, we obtain a time-scheme that experimentally is unconditionally stable, has low cost per time step, and is third-order accurate in time. We present numerical results to analyze the cost and convergence rates of the scheme. To verify the solver, we compare it to a constrained variational approach to compute equilibrium shapes that does not involve interactions with a viscous fluid. To illustrate the applicability of method, we consider a few vesicle-flow interaction problems: the sedimentation of a vesicle, interactions of one and three vesicles with a background Poiseuille flow.     

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.     

Inverse Lax-Wendroff procedure for numerical boundary conditions of conservation laws

We develop a high order finite difference numerical boundary condition for solving hyperbolic conservation laws on a Cartesian mesh. The challenge results from the wide stencil of the interior high order scheme and the fact that the boundary intersects the grids in an arbitrary fashion. Our method is based on an inverse Lax-Wendroff procedure for the inflow boundary conditions. We repeatedly use the partial differential equation to write the normal derivatives to the inflow boundary in terms of the time derivatives and the tangential derivatives. With these normal derivatives, we can then impose accurate values of ghost points near the boundary by a Taylor expansion. At outflow boundaries, we use Lagrange extrapolation or least squares extrapolation if the solution is smooth, or a weighted essentially non-oscillatory (WENO) type extrapolation if a shock is close to the boundary. Extensive numerical examples are provided to illustrate that our method is high order accurate and has good performance when applied to one and two-dimensional scalar or system cases with the physical boundary not aligned with the grids and with various boundary conditions including the solid wall boundary condition. Additional numerical cost due to our boundary treatment is discussed in some of the examples.