Compact differencing schemes for advective problems

Two-dimensional viscous flow computations on the Connecti on Machine: Unstructured meshes, upwind schemes and massively parallel computations

A three-dimensional unsteady, viscous aerodynamic analysis has been developed for the flow inside a transonic, high-through-flow, single stage compressor. The compressor stage is comprised of a low-aspect-ratio rotor and a closely coupled stator. The analysis is based on a numerical method for solving the three-dimensional Navier-Stokes equation for unsteady viscous flow through multiple turbomachinery blade rows. The method solves the fully three-dimensional Navier-Stokes equation with an implicit scheme. A two-equation turbulence model with a low-Reynolds-number modification is applied for the turbulence closure. A third-order accurate upwinding scheme is used to approximate convection terms while a second-order accurate central difference scheme is used for the discretization of the viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow field inside a transonic, high-through-flow, axial compressor stage. The numerical results are compared with available experimental data.

Laminar mixed convection over horizontal flat plate embedded in a non-Darcian porous medium with suction and injection

Numerical study of two dimensional laminar boundary layer compressible flow with pressure gradient and heat and mass transfer

The force convective heat transfer in an equilateral triangular duct of different wall heat flux configurations was analysed for the laminar hydro-dynamically developed and thermally developing flow by the use of finite volume method. Unstructured meshing was generated by multi-block technique and set of governing equations were discretized using second-order accurate up-wind scheme and numerically solved by SIMPLE Algorithm. For ensuring accuracy, grid independence study was also done. Numerical methodology was verified by comparing results with previous work and predicted results showed good agreement with them (within error of ±5 %). The different combinations of constant heat flux boundary condition were analysed and their effect on heat transfer and fluid flow for different Reynolds number was also studied. The results of different combinations were compared with the case of force convective heat transfer in the equilateral triangular duct with constant heat flux on all three walls.

On the development of a high-order compact scheme for exhibiting the switching and dissipative solution natures in the Camassa–Holm equation

This paper develops a high-order accurate gas-kinetic scheme in the framework of the finite volume method for the one- and two-dimensional flow simulations, which is an extension of the third-order accurate gas-kinetic scheme [Q.B. Li, K. Xu, and S. Fu, J. Comput. Phys., 229(2010), 6715-6731] and the second-order accurate gas-kinetic scheme [K. Xu, J. Comput. Phys., 171(2001), 289-335]. It is formed by two parts: quartic polynomial reconstruction of the macroscopic variables and fourth-order accurate flux evolution. The first part reconstructs a piecewise cell-center based quartic polynomial and a cell-vertex based quartic polynomial according to the “initial” cell average approximation of macroscopic variables to recover locally the non-equilibrium and equilibrium single particle velocity distribution functions around the cell interface. It is in view of the fact that all macroscopic variables become moments of a single particle velocity distribution function in the gas-kinetic theory. The generalized moment limiter is employed there to suppress the possible numerical oscillation. In the second part, the macroscopic flux at the cell interface is evolved in fourth-order accuracy by means of the simple particle transport mechanism in the microscopic level, i.e. free transport and the Bhatnagar-Gross-Krook (BGK) collisions. In other words, the fourth-order flux evolution is based on the solution (i.e. the particle velocity distribution function) of the BGK model for the Boltzmann equation. Several 1D and 2D test problems are numerically solved by using the proposed high-order accurate gas-kinetic scheme. By comparing with the exact solutions or the numerical solutions obtained the second-order or third-order accurate gas-kinetic scheme, the computations demonstrate that our scheme is effective and accurate for simulating invisid and viscous fluid flows, and the accuracy of the high-order GKS depends on the choice of the (numerical) collision time.

In this study, natural convection flow along a vertical wavy surface has been investigated with variable heat flux. The governing equations are transformed into dimensionless PDE by using the non-dimensional variables and then solved numerically by using an implicit finite difference scheme known as Keller Box method. The effects of the parameters amplitude of the wavy surface, ?, exponent of the variable heat flux, m, and Prandtl number on the local skin friction coefficient and local Nusselt number are shown graphically. It is found that for the negative value of exponent of the variable heat flux, m, the local skin friction coefficient increases and Nusselt number decreases but the opposite behavior is observed for the positive values of m. The comparison of limiting case with the previous study is shown through table and it is found that the solution obtained is in excellent agreement with the previous studies.

A compact fourth-order gas-kinetic scheme for the Euler and Navier–Stokes equations

In this paper, we propose and analyze a second-order accurate numerical scheme for the Poisson-Nernst-Planck equations. The proposed scheme combines a novel second-order temporal discretization with the centered finite difference method in space. By using a Taylor expansion for the logarithmic term in the chemical potential, this second-order accurate numerical scheme is able to preserve original energy dissipation. Based on the gradient flow formulation, the resulting scheme guarantees several crucial physical properties: mass conservation, positivity of ionic concentration, preservation of original energy dissipation, and steady states preservation. Remarkably, these properties are ensured without any restriction on the time step size. Furthermore, an optimal rate convergence estimate is provided for the proposed numerical scheme. Due to the non-constant mobility and the nonlinear and singular properties of the logarithm terms, a higher-order asymptotic expansion and a combination of rough and refined error estimation techniques are introduced to accomplish this analysis. A few numerical tests are provided to validate our theoretical claims.

Numerical simulation of lead-acid batteries using Keller–Box method

The two-dimensional Reynolds averaged compressible Navier-Stokes equations are solved using MacCormack's second-order accurate explicit finite difference method to simulate the separated transonic tur- bulent flowfield over an airfoil. Four different algebraic eddy viscoisity models are tested for viability to achieve turbulence closure for the class of flows considered. These models range from an unmodified boundary-layer mixing-length model to a relaxation model incorporating special considerations for the separation bubble region. Results of this study indicate the necessity for special attention to the separated flow region and suggest limits of applicability of algebraic turbulence models to these separated flowfield. each of these studies the time-dependent Reynolds averaged Navier-Stokes equations for two-dimensional compressive flow are used and tur- bulence closure is achieved by means of model equations for the Reynolds stresses. Wilcox1'2 used a first-order accurate numerical scheme and the two equation differential tur- bulence model of Saffman 12 to simulate the supersonic shock boundary-layer interaction experiment of Reda and Mur- phy 13 and the compression corner flow of Law.14 Good quan- titative agreement with the Reda and Murphy data was ob- tained, but only the qualitative features of the compression corner flow were well simulated. Using a more sophisticated second-order accurate numerical scheme, Baldwin3'4 con- sidered both the two equation differential model of Saffman and a simpler algebraic mixing-length model to simulate the hypersonic shock boundary-layer interaction experiment of Holden.15 He found the more elaborate model of Saffman to yield somewhat better results than the algebraic model, but at the cost of considerably more computing time. Good quan- titative agreement with experiment was not obtained with either model. Following Baldwin's approach all subsequent investigations have been performed using the more rigorous second-order accurate numerical scheme of Mac- Cormack.17'18 Deiwert5'6'11 considered an algebraic mixing- length model to simulate the transonic airfoil experiment of McDevitt et al. 16 while Horstman et al. 8 used a similar ap- proach to simulate their hypersonic shock boundary-layer ex- periment on an axisymmetric cylinder. In each of these studies, while qualitative features of the flows were described well, good quantitative agreement with experiment in the in- teraction regions was not obtained. Using a relaxing turbulence model Shang and Hankey7 simulated the compression corner flow of Law, and Baldwin and Rose10 simulated the flat plate flow of Reda and Murphy. In each of these studies the relaxing model was found to per- form significantly better than the simpler algebraic model and, according to Shang and Hankey, provided significantly better comparisons with measurements than were obtained by Wilcox using the two equation differential model of Saffman. In each of these studies it was essential that the full Navier- Stokes equations be considered to describe the viscous- inviscid interaction and the elliptic nature of separating-

Two-dimensional shallow-water schemes on Cartesian grids are amendable for graphics processing units and thus a convenient choice for fast flood simulations. A comparison of recent schemes and validation of important use cases is essential for developers and practitioners working with flood simulation tools. In this paper, we discuss three state-of-the-art shallow-water schemes: a first-order upwind scheme, a second-order upwind scheme, and a second-order central-upwind scheme. We analyze the advantages and disadvantages of each scheme on historical Danube river floods at three regions in Austria. We study the Lobau region as a floodplain with several small channels, the Wachau region with the meandering Danube in a steep valley, and the Marchfeld region located at the river confluence of March and Danube. The validation case studies show that the second-order schemes provide better estimates of the water levels than the first-order scheme. Still, the first order scheme is useful because it offers fast simulations and reasonable results at higher resolutions. The best trade-off between accuracy and computational effort for simulating river floods is provided by the second-order upwind scheme.

Second-order accurate upwind and centered schemes are presented in a framework that facilitates their analysis and comparison. The upwind scheme employed consists of a reconstruction step (MUSCL approach) followed by an upwind step (Roe's flux-difference splitting). The two centered schemes are of Lax-Friedrichs (L-F) type. They are the nonstaggered versions of the Nessyahu-Tadmor (N-T) scheme and the CE/SE method with epilson = 1/2. The upwind scheme is extended to the case of two spatial dimensions (2D) in a straightforward manner. The N-T and CE/SE schemes are extended in a manner similar to the 2D extensions of the CE/SE scheme by Wang and Chang for a triangular mesh and by Zhang, Yu, and Chang for a quadrilateral mesh. The slope estimates, however, are simplified. Fourier stability and accuracy analyses are carried out for these schemes for the standard 1D and the 2D quadrilateral mesh cases. In the nonstandard case of a triangular mesh, the triangles must be paired up when analyzing the upwind and N-T schemes. An observation resulting in an extended N-T scheme which is faster and uses only one-third of the storage for flow data compared with the CE/SE method is presented. Numerical results are shown. Other improvements to the schemes are discussed.

The problem of natural convection in a square cavity of a Boussinesq incompressible fluid is computed using a consistent central-difference approximation, a consistent second-order upwind scheme, and the hybrid scheme for the energy equation. It was found that, for this problem, the second-order upwind scheme yields satisfactory solutions and that an apparently low-order combination of schemes, with the second-order upwind scheme for the momentum equations and the hybrid scheme for the energy equation, can produce satisfactory predictions for this problem.