A lattice Boltzmann model for the compressible Euler equations with second‐order accuracy

In this paper, we propose a new lattice Boltzmann model for the compressible Euler equations. The model is based on a three‐energy‐level and three‐speed lattice Boltzmann equation by using a method of higher moments of the equilibrium distribution functions. In order to obtain second‐order accuracy, we employ the ghost field distribution functions to remove the non‐physical viscous parts. We also use the conditions of the higher moment of the ghost field equilibrium distribution functions to obtain the equilibrium distribution functions. In the numerical examples, we compare the numerical results of this scheme with those obtained by other lattice Boltzmann models for the compressible Euler equations. Copyright © 2008 John Wiley & Sons, Ltd.     

A lattice Boltzmann model for simulation of compressible flows

This paper presents a new model of lattice Boltzmann method for full compressible flows. On the basis of multi‐speed model, an extra potential energy distribution function is introduced to recover the full compressible Navier–Stokes equations with a flexible specific‐heat ratio and Prandtl number. The Chapman–Enskog expansion of the kinetic equations is performed, and the two‐dimension‐seventeen‐velocity density equilibrium distribution functions are obtained. The governing equations are discretized using the third order monotone upwind scheme for scalar conservation laws finite volume scheme. The van Albada limiter is used to avoid spurious oscillations. In order to verify the accuracy of this double‐distribution‐function model, the Riemann problems, Couette flows, and flows around a NACA0012 airfoil are simulated. It is found that the proposed lattice Boltzmann model is suitable for compressible flows, even for strong shock wave problem, which has an extremely large pressure ratio, 100,000. Copyright © 2014 John Wiley & Sons, Ltd.     

A multi-entropy-level lattice Boltzmann model for the one-dimensional compressible Euler equations

In this paper, we propose a new lattice Boltzmann model for the one-dimensional compressible Euler equations. The new model is based on a three-entropy-level and three-speed lattice Boltzmann equation by using a method of higher-order moments of the equilibrium distribution functions. In order to obtain the second-order accuracy model, we employ the ghost field distribution functions to remove the non-physical dissipation terms in the Euler equations. We also use the conditions of the higher-order moments of the ghost field equilibrium distribution functions to obtain the equilibrium distribution functions. The numerical examples show that the numerical results can be compared with those classical methods.     

Simulate natural convection in a 30 MeV cyclotron with effect of irradiation,lattice Boltzmann 3D compressible thermal model

The lattice Boltzmann equation method in three dimensions used to analyze compressible thermal flow in 30 MeV cyclotrons. Cyclotron produced radionuclides have diagnostic applications in nuclear medicine. This cyclotron has horizontal conically shaped and window air cooling in the front and water cooling of the body. Gas heated by irradiation. This scheme is a new type of simulation method for solving the time dependent Navier–Stokes equations in a compressible flow regime. The improved model is convenient to compromise the high accuracy and stability. The included dispersion term can effectively reduce the numerical oscillation at discontinuity. The lattice Boltzmann scheme with uniform mesh resolution is applied as a numerical research tool.     

Lattice Boltzmann method for the fractional sub‐diffusion equation

A lattice Boltzmann model for the fractional sub‐diffusion equation is presented. By using the Chapman–Enskog expansion and the multiscale time expansion, several higher‐order moments of equilibrium distribution functions and a series of partial differential equations in different time scales are obtained. Furthermore, the modified partial differential equation of the fractional sub‐diffusion equation with the second‐order truncation error is obtained. In the numerical simulations, comparisons between numerical results of the lattice Boltzmann models and exact solutions are given. The numerical results agree well with the classical ones. Copyright © 2015 John Wiley & Sons, Ltd.     

On coupling the Reynolds‐averaged Navier–Stokes equations with two‐equation turbulence model equations

Two methods for coupling the Reynolds‐averaged Navier–Stokes equations with the q–ω turbulence model equations on structured grid systems have been studied; namely a loosely coupled method and a strongly coupled method. The loosely coupled method first solves the Navier–Stokes equations with the turbulent viscosity fixed. In a subsequent step, the turbulence model equations are solved with all flow quantities fixed. On the other hand, the strongly coupled method solves the Reynolds‐averaged Navier–Stokes equations and the turbulence model equations simultaneously. In this paper, numerical stabilities of both methods in conjunction with the approximated factorization‐alternative direction implicit method are analysed. The effect of the turbulent kinetic energy terms in the governing equations on the convergence characteristics is also studied. The performance of the two methods is compared for several two‐ and three‐dimensional problems. Copyright © 2005 John Wiley & Sons, Ltd.     

A higher‐order accuracy lattice Boltzmann model for the wave equation

A lattice Boltzmann model with higher‐order accuracy for the wave motion is proposed. The new model is based on the technique of the higher‐order moment of equilibrium distribution functions and a series of lattice Boltzmann equations in different time scales. The forms of moments are derived from the binary wave equation by designing the higher‐order dissipation and dispersion terms. The numerical results agree well with classical ones. Copyright © 2008 John Wiley & Sons, Ltd.     

Non‐intrusive reduced‐order modelling of the Navier–Stokes equations based on RBF interpolation

We present a new non‐intrusive model reduction method for the Navier–Stokes equations. The method replaces the traditional approach of projecting the equations onto the reduced space with a radial basis function (RBF) multi‐dimensional interpolation. The main point of this method is to construct a number of multi‐dimensional interpolation functions using the RBF scatter multi‐dimensional interpolation method. The interpolation functions are used to calculate POD coefficients at each time step from POD coefficients at earlier time steps. The advantage of this method is that it does not require modifications to the source code (which would otherwise be very cumbersome), as it is independent of the governing equations of the system. Another advantage of this method is that it avoids the stability problem of POD/Galerkin. The novelty of this work lies in the application of RBF interpolation and POD to construct the reduced‐order model for the Navier–Stokes equations. Another novelty is the verification and validation of numerical examples (a lock exchange problem and a flow past a cylinder problem) using unstructured adaptive finite element ocean model. The results obtained show that CPU times are reduced by several orders of magnitude whilst the accuracy is maintained in comparison with the corresponding high‐fidelity models. Copyright © 2015 John Wiley & Sons, Ltd.     

A stabilized finite element method for the incompressible Navier–Stokes equations using a hierarchical basis

Stabilized finite element methods have been shown to yield robust, accurate numerical solutions to both the compressible and incompressible Navier–Stokes equations for laminar and turbulent flows. The present work focuses on the application of higher‐order, hierarchical basis functions to the incompressible Navier–Stokes equations using a stabilized finite element method. It is shown on a variety of problems that the most cost‐effective simulations (in terms of CPU time, memory, and disk storage) can be obtained using higher‐order basis functions when compared with the traditional linear basis. In addition, algorithms will be presented for the efficient implementation of these methods within the traditional finite element data structures. Copyright © 2001 John Wiley & Sons, Ltd.     

Fourth‐order accurate compact‐difference discretization method for Helmholtz and incompressible Navier–Stokes equations

A fourth‐order accurate solution method for the three‐dimensional Helmholtz equations is described that is based on a compact finite‐difference stencil for the Laplace operator. Similar discretization methods for the Poisson equation have been presented by various researchers for Dirichlet boundary conditions. Here, the complicated issue of imposing Neumann boundary conditions is described in detail. The method is then applied to model Helmholtz problems to verify the accuracy of the discretization method. The implementation of the solution method is also described. The Helmholtz solver is used as the basis for a fourth‐order accurate solver for the incompressible Navier–Stokes equations. Numerical results obtained with this Navier–Stokes solver for the temporal evolution of a three‐dimensional instability in a counter‐rotating vortex pair are discussed. The time‐accurate Navier–Stokes simulations show the resolving properties of the developed discretization method and the correct prediction of the initial growth rate of the three‐dimensional instability in the vortex pair. Copyright © 2004 John Wiley & Sons, Ltd.     

A direct‐forcing pressure‐based lattice Boltzmann method for solving fluid–particle interaction problems

A direct‐forcing immersed boundary‐lattice Boltzmann method (IB–LBM) is developed to simulate fluid–particle interaction problems. This method uses the pressure‐based LBM to solve the incompressible flow field and the immersed boundary method to handle the fluid–particle interactions. The pressure‐based LBM uses the pressure distribution functions instead of the density distribution functions as the independent dynamic variables. The main idea is to explicitly eliminate the compressible effect due to the density fluctuation. In the IB method, a direct‐forcing method is introduced to capture the particle motion. It directly computes an IB force density at each lattice grid from the differences between the pressure distribution functions obtained by the LBM and the equilibrium pressure distribution functions computed from the particle velocity. By applying this direct‐forcing method, the IB–LBM becomes a purely LBM version. Also, by applying the Gauss theorem, the formulas for computing the force and the torque acting on the particle from the flows are derived from the volume integrals over the particle volume instead of from the surface integrals over the particle surface. The order of accuracy of the IB–LBM is demonstrated on the errors of velocity field, wall stress, and gradients of velocity and pressure. As a demonstration of the efficiency and capabilities of the new method, sedimentation of a large number of spherical particles in an enclosure is simulated. Copyright © 2010 John Wiley & Sons, Ltd.     

A reduced‐order approach for optimal control of fluids using proper orthogonal decomposition

In this article, a reduced‐order modeling approach, suitable for active control of fluid dynamical systems, based on proper orthogonal decomposition (POD) is presented. The rationale behind the reduced‐order modeling is that numerical simulation of Navier–Stokes equations is still too costly for the purpose of optimization and control of unsteady flows. The possibility of obtaining reduced‐order models that reduce the computational complexity associated with the Navier–Stokes equations is examined while capturing the essential dynamics by using the POD. The POD allows the extraction of a reduced set of basis functions, perhaps just a few, from a computational or experimental database through an eigenvalue analysis. The solution is then obtained as a linear combination of this reduced set of basis functions by means of Galerkin projection. This makes it attractive for optimal control and estimation of systems governed by partial differential equations (PDEs). It is used here in active control of fluid flows governed by the Navier–Stokes equations. In particular, flow over a backward‐facing step is considered. Reduced‐order models/low‐dimensional dynamical models for this system are obtained using POD basis functions (global) from the finite element discretizations of the Navier–Stokes equations. Their effectiveness in flow control applications is shown on a recirculation control problem using blowing on the channel boundary. Implementational issues are discussed and numerical experiments are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

Modelling and simulation of fires in vehicle tunnels

Applying a low‐Mach asymptotic for the compressible Navier–Stokes equations, we derive a new fluid dynamics model,which should be capable to model large temperature differences in combination with the low‐Mach number limit. The model is used to simulate fires in vehicle tunnels, where the standard Boussinesq‐approximation for the incompressible Navier–Stokes seems to be inappropriate due to the high temperatures developing in the tunnel. The model is implemented using a modified finite‐difference approach for the incompressible Navier–Stokes equations and tested in some realistic fire events. Copyright © 2004 John Wiley & Sons, Ltd.     

A comparative study of direct‐forcing immersed boundary‐lattice Boltzmann methods for stationary complex boundaries

In this study, we assess several interface schemes for stationary complex boundary flows under the direct‐forcing immersed boundary‐lattice Boltzmann methods (IB‐LBM) based on a split‐forcing lattice Boltzmann equation (LBE). Our strategy is to couple various interface schemes, which were adopted in the previous direct‐forcing immersed boundary methods (IBM), with the split‐forcing LBE, which enables us to directly use the direct‐forcing concept in the lattice Boltzmann calculation algorithm with a second‐order accuracy without involving the Navier–Stokes equation. In this study, we investigate not only common diffuse interface schemes but also a sharp interface scheme. For the diffuse interface scheme, we consider explicit and implicit interface schemes. In the calculation of velocity interpolation and force distribution, we use the 2‐ and 4‐point discrete delta functions, which give the second‐order approximation. For the sharp interface scheme, we deal with the exterior sharp interface scheme, where we impose the force density on exterior (solid) nodes nearest to the boundary. All tested schemes show a second‐order overall accuracy when the simulation results of the Taylor–Green decaying vortex are compared with the analytical solutions. It is also confirmed that for stationary complex boundary flows, the sharper the interface scheme, the more accurate the results are. In the simulation of flows past a circular cylinder, the results from each interface scheme are comparable to those from other corresponding numerical schemes. Copyright © 2010 John Wiley & Sons, Ltd.     

Fourth‐order compact formulation of Navier–Stokes equations and driven cavity flow at high Reynolds numbers

A new fourth‐order compact formulation for the steady 2‐D incompressible Navier–Stokes equations is presented. The formulation is in the same form of the Navier–Stokes equations such that any numerical method that solve the Navier–Stokes equations can easily be applied to this fourth‐order compact formulation. In particular, in this work the formulation is solved with an efficient numerical method that requires the solution of tridiagonal systems using a fine grid mesh of 601 × 601. Using this formulation, the steady 2‐D incompressible flow in a driven cavity is solved up to Reynolds number with Re = 20 000 fourth‐order spatial accuracy. Detailed solutions are presented. Copyright © 2005 John Wiley & Sons, Ltd.     

A Galerkin spectral method based on helical‐wave decomposition for the incompressible Navier–Stokes equations

This paper presents a global Galerkin spectral method for solving the incompressible Navier–Stokes equations in three‐dimensional bounded domains. The method is based on helical‐wave decomposition (HWD), which uses the vector eigenfunctions of the curl operator as orthogonal basis functions. We shall first review the general theory of HWD in an arbitrary simply connected domain, along with some new developments. We then employ the HWD to construct a Galerkin spectral method. The current method innovates the existing HWD‐based spectral method by (a) adding a series of auxiliary fields to the HWD of the velocity field to fulfill the no‐slip boundary condition and to settle the convergence problem of the HWD of the curl fields, and (b) providing a pseudo‐spectral method that utilizes a fast spherical harmonic transform algorithm and Gaussian quadrature to calculate the nonlinear term in the Navier–Stokes equations. The auxiliary fields are uniquely determined by solving the Stokes and Stokes‐like equations under adequate boundary conditions. The implementation of the method under the spherical geometry is presented in detail. Several numerical examples are provided to validate the proposed method. The method can be easily extended to other domains once the helical‐wave bases, which depend only on the geometry of the domains, are available. Copyright © 2015 John Wiley & Sons, Ltd.     

A rotating reference frame‐based lattice Boltzmann flux solver for simulation of turbomachinery flows

In this paper, the newly developed lattice Boltzmann flux solver (LBFS) is developed into a version in the rotating frame of reference for simulation of turbomachinery flows. LBFS is a finite volume solver for the solution of macroscopic governing differential equations. Unlike conventional upwind or Godunov‐type flux solvers which are constructed by considering the mathematical properties of Euler equations, it evaluates numerical fluxes at the cell interface by reconstructing local solution of lattice Boltzmann equation (LBE). In other words, the numerical fluxes are physically determined rather than by some mathematical approximation. The LBE is herein expressed in a relative frame of reference in order to correctly recover the macroscopic equations, which is also the basis of LBFS. To solve the LBE, an appropriate lattice Boltzmann model needs to be established in advance. This includes both the determinations of the discrete velocity model and its associated equilibrium distribution functions. Particularly, a simple and effective D1Q4 model is adopted, and the equilibrium distribution functions could be efficiently obtained by using the direct method. The present LBFS is validated by several inviscid and viscous test cases. The numerical results demonstrate that it could be well applied to typical and complex turbomachinery flows with favorable accuracy. It is also shown that LBFS has a delicate dissipation mechanism and is thus free of some artificial fixes, which are often needed in conventional schemes. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of hypersonic flows using finite elements with Taylor–Galerkin scheme

The aim of the present work is to introduce a formulation for the numerical analysis of three‐dimensional thermochemical non‐equilibrium hypersonic flows, using the finite element method and the Taylor–Galerkin scheme and adopting Park's 2‐temperature, 5‐species (N2, O2, NO, N and O) and 17‐reaction model. Examples using Euler and Navier–Stokes equations are included and compared with experimental and numerical works presented by other authors. The results are close to those analysed by other researches and a good computational performance was obtained. Copyright © 2004 John Wiley & Sons, Ltd.     

Time‐linearized time‐harmonic 3‐D Navier–Stokes shock‐capturing schemes

In the present paper, a numerical method for the computation of time‐harmonic flows, using the time‐linearized compressible Reynolds‐averaged Navier–Stokes equations is developed and validated. The method is based on the linearization of the discretized nonlinear equations. The convective fluxes are discretized using an O(Δx) MUSCL scheme with van Leer flux‐vector‐splitting. Unsteady perturbations of the turbulent stresses are linearized using a frozen‐turbulence‐Reynolds‐number hypothesis, to approximate eddy‐viscosity perturbations. The resulting linear system is solved using a pseudo‐time‐marching implicit ADI‐AF (alternating‐directions‐implicit approximate‐factorization) procedure with local pseudo‐time‐steps, corresponding to a matrix‐successive‐underrelaxation procedure. The stability issues associated with the pseudo‐time‐marching solution of the time‐linearized Navier–Stokes equations are discussed. Comparison of computations with measurements and with time‐nonlinear computations for 3‐D shock‐wave oscillation in a square duct, for various back‐pressure fluctuation frequencies (180, 80, 20 and 10 Hz), assesses the shock‐capturing capability of the time‐linearized scheme. Copyright © 2007 John Wiley & Sons, Ltd.