一维含化学反应流体力学方程组的ALE间断有限元方法

研究一维含化学反应流体力学方程组的数值模拟方法.结合理想气体状态方程并利用HLLC解法器在各个单元边界处的数值通量，给出ALE间断有限元方法.高阶计算时，使用TVD斜率限制器对数值解可能产生的非物理振荡进行抑制.结果表明：该算法能够保持物理量的守恒性和高精度，并能够清晰地捕捉爆轰波的结构特征.     

A FINITE DIFFERENCE SCHEME FOR SOLVING THE NONLINEAR POISSON-BOLTZMANN EQUATION MODELING CHARGED SPHERES

In this work, we propose an efficient numerical method for computing the electrostaticinteraction between two like-charged spherical particles which is governed by the nonlinearPoisson-Boltzmann equation. The nonlinear problem is solved by a monotone iterativemethod which leads to a sequence of linearized equations. A modified central finite differ-ence scheme is developed to solve the linearized equations on an exterior irregular domainusing a uniform Cartesian grid. With uniform grids, the method is simple, and as aconsequence, multigrid solvers can be employed to speed up the convergence. Numericalexperiments on cases with two isolated spheres and two spheres confined in a chargedcylindrical pore are carried out using the proposed method. Our numerical schemes arefound efficient and the numerical results are found in good agreement with the previouspublished results.     

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 Schur–Fourier decomposition for the efficient solution of high‐order Poisson equations on loosely coupled parallel computers

In this paper a parallel direct Schur–Fourier decomposition (DSFD) algorithm for the direct solution of arbitrary order discrete Poisson equations on parallel computers is proposed. It is based on a combination of a Direct Schur method and a Fourier decomposition and allows to solve each Poisson equation almost to machine accuracy using only one communication episode. Thus, it is well suited for loosely coupled parallel computers, that have a high network latency compared with the CPU performance. Several three‐dimensional direct numerical simulations (DNS) of wall‐bounded turbulent incompressible flows have been carried out using the DSFD algorithm. Numerical examples illustrating the robustness and scalability of the method on a PC cluster with a conventional 100 Mbits/s network are also presented. Copyright © 2005 John Wiley & Sons, Ltd.     

Navier–Stokes computation of transonic vortices over a round leading edge delta wing

A 3D Navier–Stokes solver has been developed to simulate laminar compressible flow over quadrilateral wings. The finite volume technique is employed for spatial discretization with a novel variant for the viscous fluxes. An explicit three-stage Runge–Kutta scheme is used for time integration, taking local time steps according to the linear stability condition derived for application to the Navier–Stokes equations. The code is applied to compute primary and secondary separation vortices at transonic speeds over a 65° swept delta wing with round leading edges and cropped tips. The results are compared with experimental data and Euler solutions, and Reynolds number effects are investigated.     

A finite volume method for viscous incompressible flows using a consistent flux reconstruction scheme

An incompressible Navier–Stokes solver using curvilinear body‐fitted collocated grid has been developed to solve unconfined flow past arbitrary two‐dimensional body geometries. In this solver, the full Navier–Stokes equations have been solved numerically in the physical plane itself without using any transformation to the computational plane. For the proper coupling of pressure and velocity field on collocated grid, a new scheme, designated ‘consistent flux reconstruction’ (CFR) scheme, has been developed. In this scheme, the cell face centre velocities are obtained explicitly by solving the momentum equations at the centre of the cell faces. The velocities at the cell centres are also updated explicitly by solving the momentum equations at the cell centres. By resorting to such a fully explicit treatment considerable simplification has been achieved compared to earlier approaches. In the present investigation the solver has been applied to unconfined flow past a square cylinder at zero and non‐zero incidence at low and moderate Reynolds numbers and reasonably good agreement has been obtained with results available from literature. Copyright © 2006 John Wiley & Sons, Ltd.     

A novel fully implicit finite volume method applied to the lid‐driven cavity problem—Part I: High Reynolds number flow calculations

A novel implicit cell‐vertex finite volume method is described for the solution of the Navier–Stokes equations at high Reynolds numbers. The key idea is the elimination of the pressure term from the momentum equation by multiplying the momentum equation with the unit normal vector to a control volume boundary and integrating thereafter around this boundary. The resulting equations are expressed solely in terms of the velocity components. Thus any difficulties with pressure or vorticity boundary conditions are circumvented and the number of primary variables that need to be determined equals the number of space dimensions. The method is applied to both the steady and unsteady two‐dimensional lid‐driven cavity problem at Reynolds numbers up to 10000. Results are compared with those in the literature and show excellent agreement. Copyright © 2003 John Wiley & Sons, Ltd.     

Multiple grid methods for classical molecular dynamics

Presented in the context of classical molecular mechanics and dynamics are multilevel summation methods for the fast calculation of energies/forces for pairwise interactions, which are based on the hierarchical interpolation of interaction potentials on multiple grids. The concepts and details underlying multigrid interpolation are described. For integration of molecular dynamics the use of different time steps for different interactions allows longer time steps for many of the interactions, and this can be combined with multiple grids in space. Comparison is made to the fast multipole method, and evidence is presented suggesting that for molecular simulations multigrid methods may be superior to the fast multipole method and other tree methods.     

Implementation of all‐Mach Roe‐type schemes in fully implicit CFD solvers – demonstration for wind turbine flows

This paper presents the implementation of all‐Mach Roe‐type schemes in a fully implicit CFD solver. Simple 2D cases, such as the flow around inviscid and viscous aerofoils, were used for an initial validation of these methods, along with more challenging computations consisting of the 3D flow around the Model Experiments in Controlled Conditions wind turbine, in parked and rotating conditions. This work is motivated by the increased interest of the wind turbine industry in larger diameter wind turbines where compressibility effects near the blade tips may be important. Instead of using an incompressible flow solver, this paper explores the option of modifying an existing, efficient, compressible flow solver for use at lower Mach numbers. The good performance of the Roe solver and its popularity influenced the selection of schemes for this work. The results suggest that effective all‐Mach solutions are possible with implicit solvers, and the paper defines the implementation of the new fluxes and Jacobian, including an investigation of some numerical parameters, using as platform the Helicopter Multi‐Block solver of Liverpool University. Copyright © 2013 John Wiley & Sons, Ltd.     

A robust kernel-based solver for variable-order time fractional PDEs under 2D/3D irregular domains

This study presents a robust kernel-based collocation method (KBCM) for solving multi-term variable-order time fractional partial differential equations (VOTFPDEs). In the proposed method, Radial basis functions (RBFs) and Muntz polynomials basis (MPB) are implemented to discretize the spatial and temporal derivative terms in the VOTFPDEs, respectively. Due to the properties of the RBFs, the spatial discretization in the proposed method is mathematically simple and truly meshless, which avoids troublesome mesh generation for high-dimensional problems involving irregular geometries. Due to the properties of the MPB, only few temporal discretization is required to achieve the satisfactory accuracy. Numerical efficiency of the proposed method is investigated under several typical examples.