Deforming composite grids for solving fluid structure problems

We describe a mixed Eulerian–Lagrangian approach for solving fluid–structure interaction (FSI) problems. The technique, which uses deforming composite grids (DCG), is applied to FSI problems that couple high speed compressible flow with elastic solids. The fluid and solid domains are discretized with composite overlapping grids. Curvilinear grids are aligned with each interface and these grids deform as the interface evolves. The majority of grid points in the fluid domain generally belong to background Cartesian grids which do not move during a simulation. The FSI-DCG approach allows large displacements of the interfaces while retaining high quality grids. Efficiency is obtained through the use of structured grids and Cartesian grids. The governing equations in the fluid and solid domains are evolved in a partitioned approach. We solve the compressible Euler equations in the fluid domains using a high-order Godunov finite-volume scheme. We solve the linear elastodynamic equations in the solid domains using a second-order upwind scheme. We develop interface approximations based on the solution of a fluid–solid Riemann problem that results in a stable scheme even for the difficult case of light solids coupled to heavy fluids. The FSI-DCG approach is verified for three problems with known solutions, an elastic-piston problem, the superseismic shock problem and a deforming diffuser. In addition, a self convergence study is performed for an elastic shock hitting a fluid filled cavity. The overall FSI-DCG scheme is shown to be second-order accurate in the max-norm for smooth solutions, and robust and stable for problems with discontinuous solutions for a wide range of constitutive parameters.     

A third-order finite-volume residual-based scheme for the 2D Euler equations on unstructured grids

A residual-based (RB) scheme relies on the vanishing of residual at the steady-state to design a transient first-order dissipation, which becomes high-order at steady-state. Initially designed within a finite-difference framework for computations of compressible flows on structured grids, the RB schemes displayed good convergence, accuracy and shock-capturing properties which motivated their extension to unstructured grids using a finite volume (FV) method. A second-order formulation of the FV–RB scheme for compressible flows on general unstructured grids was presented in a previous paper. The present paper describes the derivation of a third-order FV–RB scheme and its application to hyperbolic model problems as well as subsonic, transonic and supersonic internal and external inviscid flows.     

Comparison of Several Spatial Discretizations for the Navier–Stokes Equations

Grid convergence studies for subsonic and transonic flows over airfoils are presented in order to compare the accuracy of several spatial discretizations for the compressible Navier–Stokes equations. The discretizations include the following schemes for the inviscid fluxes: (1) second-order-accurate centered differences with third-order matrix numerical dissipation, (2) the second-order convective upstream split pressure scheme (CUSP), (3) third-order upwind-biased differencing with Roe's flux-difference splitting, and (4) fourth-order centered differences with third-order matrix numerical dissipation. The first three are combined with second-order differencing for the grid metrics and viscous terms. The fourth discretization uses fourth-order differencing for the grid metrics and viscous terms, as well as higher-order approximations near boundaries and for the numerical integration used to calculate forces and moments. The results indicate that the discretization using higher-order approximations for all terms is substantially more accurate than the others, producing less than two percent numerical error in lift and drag components on grids with less than 13,000 nodes for subsonic cases and less than 18,000 nodes for transonic cases. Since the cost per grid node of all of the discretizations studied is comparable, the higher-order discretization produces solutions of a given accuracy much more efficiently than the others.     

分区拼接网格算法数值模拟超声速复杂流场

以分区拼接网格数值模拟方法为研究对象,对二阶Godunov方法在拼接网格中流场的数值模拟进行了研究,发展了适用于Godunov格式的通量守恒型算法,结合二阶Godunov有限体积法离散非定常Euler方程,数值模拟了捆绑火箭及航天飞机的超声速流场,计算结果正确描述了流场中的激波相交、反射等干扰特性.     

Block structured adaptive mesh and time refinement for hybrid,hyperbolic + N-body systems

We present a new numerical algorithm for the solution of coupled collisional and collisionless systems, based on the block structured adaptive mesh and time refinement strategy (AMR). We describe the issues associated with the discretization of the system equations and the synchronization of the numerical solution on the hierarchy of grid levels. We implement a code based on a higher order, conservative and directionally unsplit Godunov’s method for hydrodynamics; a symmetric, time centered modified symplectic scheme for collisionless component; and a multilevel, multigrid relaxation algorithm for the elliptic equation coupling the two components. Numerical results that illustrate the accuracy of the code and the relative merit of various implemented schemes are also presented.     

四阶紧致差分格式对静力模式中垂直网格计算频散性的影响

为了说明四阶紧致差分格式在大气和海洋数值模式中的潜在价值,提出一种通用方法,推导静力线性斜压适应方程组在微分和差分情况下的频散关系,水平尺度分100 km,10 km和1 km三种情况,从频率、水平群速和垂直群速方面,对采用二阶中央差和四阶紧致差分格式情况下,非跳点网格(N网格)、Lorenz网格(L网格)、Charney-Phillips网格(CP网格)、Lorenz时间跳点网格(LTS网格)和Charney-Phillips时间跳点网格(CPTS网格)的计算特性进行比较,发现采用高精度的四阶紧致差分格式总体上可以明显减少上述三种水平尺度波动在N网格、CP网格、L网格和CPTS网格上的频率、水平群速和垂直群速误差,但对LTS网格,采用四阶紧致差分格式,会使得计算水平群速和垂直群速误差变大.     

A simple multigrid scheme for solving the Poisson equation with arbitrary domain boundaries

We present a new multigrid scheme for solving the Poisson equation with Dirichlet boundary conditions on a Cartesian grid with irregular domain boundaries. This scheme was developed in the context of the Adaptive Mesh Refinement (AMR) schemes based on a graded-octree data structure. The Poisson equation is solved on a level-by-level basis, using a “one-way interface” scheme in which boundary conditions are interpolated from the previous coarser level solution. Such a scheme is particularly well suited for self-gravitating astrophysical flows requiring an adaptive time stepping strategy. By constructing a multigrid hierarchy covering the active cells of each AMR level, we have designed a memory-efficient algorithm that can benefit fully from the multigrid acceleration. We present a simple method for capturing the boundary conditions across the multigrid hierarchy, based on a second-order accurate reconstruction of the boundaries of the multigrid levels. In case of very complex boundaries, small scale features become smaller than the discretization cell size of coarse multigrid levels and convergence problems arise. We propose a simple solution to address these issues. Using our scheme, the convergence rate usually depends on the grid size for complex grids, but good linear convergence is maintained. The proposed method was successfully implemented on distributed memory architectures in the RAMSES code, for which we present and discuss convergence and accuracy properties as well as timing performances.     

Upwind schemes for the wave equation in second-order form

We develop new high-order accurate upwind schemes for the wave equation in second-order form. These schemes are developed directly for the equations in second-order form, as opposed to transforming the equations to a first-order hyperbolic system. The schemes are based on the solution to a local Riemann-type problem that uses d’Alembert’s exact solution. We construct conservative finite difference approximations, although finite volume approximations are also possible. High-order accuracy is obtained using a space-time procedure which requires only two discrete time levels. The advantages of our approach include efficiency in both memory and speed together with accuracy and robustness. The stability and accuracy of the approximations in one and two space dimensions are studied through normal-mode analysis. The form of the dissipation and dispersion introduced by the schemes is elucidated from the modified equations. Upwind schemes are implemented and verified in one dimension for approximations up to sixth-order accuracy, and in two dimensions for approximations up to fourth-order accuracy. Numerical computations demonstrate the attractive properties of the approach for solutions with varying degrees of smoothness.     

An evaluation of the FCT method for high-speed flows on structured overlapping grids

This study considers the development and assessment of a flux-corrected transport (FCT) algorithm for simulating high-speed flows on structured overlapping grids. This class of algorithm shows promise for solving some difficult highly-nonlinear problems where robustness and control of certain features, such as maintaining positive densities, is important. Complex, possibly moving, geometry is treated through the use of structured overlapping grids. Adaptive mesh refinement (AMR) is employed to ensure sharp resolution of discontinuities in an efficient manner. Improvements to the FCT algorithm are proposed for the treatment of strong rarefaction waves as well as rarefaction waves containing a sonic point. Simulation results are obtained for a set of test problems and the convergence characteristics are demonstrated and compared to a high-resolution Godunov method. The problems considered are an isolated shock, an isolated contact, a modified Sod shock tube problem, a two-shock Riemann problem, the Shu–Osher test problem, shock impingement on single cylinder, and irregular Mach reflection of a strong shock striking an inclined plane.     

守恒律方程组的大时间步长叠波格式

大时间步长叠波格式最初思想为LeVeque提出的大时间步长Godunov格式,通过叠加间断分解发出的强波来构造数值格式.原方法只给出了间断强波的穿越叠加方法,文章对其进行了完善,并推广到多维.针对膨胀波提出了一种网格单元分解法可以自动满足熵条件,避免出现非物理解.给出了格式的具体计算公式,并用单个守恒律方程、一维/多维Euler方程组进行了数值计算.计算结果表明,新格式除了可以采用大时间步长的优点外,在一定范围内随CFL数增加其耗散反而更低,因而对激波接触间断膨胀波的分辨率更高.     

Finite volume treatment of dispersion-relation-preserving and optimized prefactored compact schemes for wave propagation

In developing suitable numerical techniques for computational aero-acoustics, the dispersion-relation-preserving (DRP) scheme by Tam and co-workers and the optimized prefactored compact (OPC) scheme by Ashcroft and Zhang have shown desirable properties of reducing both dissipative and dispersive errors. These schemes, originally based on the finite difference, attempt to optimize the coefficients for better resolution of short waves with respect to the computational grid while maintaining pre-determined formal orders of accuracy. In the present study, finite volume formulations of both schemes are presented to better handle the nonlinearity and complex geometry encountered in many engineering applications. Linear and nonlinear wave equations, with and without viscous dissipation, have been adopted as the test problems. Highlighting the principal characteristics of the schemes and utilizing linear and nonlinear wave equations with different wavelengths as the test cases, the performance of these approaches is documented. For the linear wave equation, there is no major difference between the DRP and OPC schemes. For the nonlinear wave equations, the finite volume version of both DRP and OPC schemes offers substantially better solutions in regions of high gradient or discontinuity.     

Adaptive mesh refinement based on high order finite difference WENO scheme for multi-scale simulations

In this paper, we propose a finite difference AMR-WENO method for hyperbolic conservation laws. The proposed method combines the adaptive mesh refinement (AMR) framework  and  with the high order finite difference weighted essentially non-oscillatory (WENO) method in space and the total variation diminishing (TVD) Runge–Kutta (RK) method in time (WENO-RK)  and  by a high order coupling. Our goal is to realize mesh adaptivity in the AMR framework, while maintaining very high (higher than second) order accuracy of the WENO-RK method in the finite difference setting. The high order coupling of AMR and WENO-RK is accomplished by high order prolongation in both space (WENO interpolation) and time (Hermite interpolation) from coarse to fine grid solutions, and at ghost points. The resulting AMR-WENO method is accurate, robust and efficient, due to the mesh adaptivity and very high order spatial and temporal accuracy. We have experimented with both the third and the fifth order AMR-WENO schemes. We demonstrate the accuracy of the proposed scheme using smooth test problems, and their quality and efficiency using several 1D and 2D nonlinear hyperbolic problems with very challenging initial conditions. The AMR solutions are observed to perform as well as, and in some cases even better than, the corresponding uniform fine grid solutions. We conclude that there is significant improvement of the fifth order AMR-WENO over the third order one, not only in accuracy for smooth problems, but also in its ability in resolving complicated solution structures, due to the very low numerical diffusion of high order schemes. In our work, we found that it is difficult to design a robust AMR-WENO scheme that is both conservative and high order (higher than second order), due to the mass inconsistency of coarse and fine grid solutions at the initial stage in a finite difference scheme. Resolving these issues as well as conducting comprehensive evaluation of computational efficiency constitute our future work.     

A high-resolution gas-kinetic scheme with minimized dispersion and controllable dissipation reconstruction

In order to simulate multiscale problems such as turbulent flows effectively, the high-order accurate reconstruction based on mini- mized dispersion and controllable dissipation (MDCD) is implemented in the second-order accurate gas-kinetic scheme (GKS) to improve the accuracy and resolution. MDCD is firstly extended to non-uniform grids through the modification of dissipation and dispersion coefficients for uniform grids based on the local stretch ratio. Remarkable improvements in accuracy and resolution are achieved on general grids. Then a new scheme, MDCD-GKS is constructed, with the help of MDCD reconstruction, not only for conservative variables, but also for their gradients. MDCD-GKS shows good accuracy and efficiency in typical numerical tests. MDCD-GKS is also coupled with the improved delayed detached-eddy simulation (IDDES) hybrid model and applied in the fine simulation of turbulent flow around a cylinder, and the prediction is in good agreement with experiments when using the relatively coarse grid. The high accuracy and resolution of the developed GKS guarantee its high efficiency in practical applications.     

A class of finite difference schemes with low dispersion and controllable dissipation for DNS of compressible turbulence

In this paper, a class of finite difference schemes which achieves low dispersion and controllable dissipation in smooth region and robust shock-capturing capabilities in the vicinity of discontinuities is presented. Firstly, a sufficient condition for semi-discrete finite difference schemes to have independent dispersion and dissipation is derived. This condition enables a novel approach to separately optimize the dissipation and dispersion properties of finite difference schemes and a class of schemes with minimized dispersion and controllable dissipation is thus obtained. Secondly, for the purpose of shock-capturing, one of these schemes is used as the linear part of the WENO scheme with symmetrical stencils to constructed an improved WENO scheme. At last, the improved WENO scheme is blended with its linear counterpart to form a new hybrid scheme for practical applications. The proposed scheme is accurate, flexible and robust. The accuracy and resolution of the proposed scheme are tested by the solutions of several benchmark test cases. The performance of this scheme is further demonstrated by its application in the direct numerical simulation of compressible turbulent channel flow between isothermal walls.     

求解Euler方程的高精度的空间--时间守恒格式

空间-时间守恒(STC)格式是近年来发展出的一种计算格式,在现有的STC格式构造过程中,流动变量在解元中的分布都用其一阶Taylor展开式来表示.STC格式的精度与所采用的Taylor展开式的阶数有关.该文采用流动变量的二阶Taylor展开式来表示其在解元上的分布、构造出了求解一维Euler方程的STC格式.用该格式对几个问题进行了计算,将计算结果与精确解进行了比较,比较表明该格式有较高的精度.     

Numerical wave propagation on non-uniform one-dimensional staggered grids

The wave propagation behaviour of centered difference schemes on one-dimensional non-uniform staggered grids is investigated. Previous results for the linear advection equation are extended to the case of the shallow water equations on staggered grids. For waves of a given frequency, the wave field is decomposed into right- and left-propagating components, and a wave energy conservation law is derived in terms of these components. For slowly varying grids, separate evolution equations for the right- and left-propagating components are derived, leading to the result that there is asymptotically no reflection in the limit of a slowly varying grid, provided that waves of that frequency are resolvable. However, there will be reflection from any location at which the wave group velocity goes to zero. The possibility for wave energy to tunnel through a narrow region of the grid too coarse for propagation is noted. Grids with an abrupt jump in resolution are also investigated. It is possible to tailor the scheme at the jump to minimize spurious wave reflection over a range of frequencies provided the waves are resolvable on both sides of the jump. However, it does not appear possible to avoid complete reflection, except by introducing extra dissipation terms, if the waves are not resolvable on one side of the jump. An example is presented of a second-order accurate scheme that spontaneously radiates waves from the resolution jump.     

A kernel-free boundary integral method for elliptic boundary value problems

This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.     

A semi-discrete central scheme for magnetohydrodynamics on orthogonal–curvilinear grids

This work describes a novel scheme for the equations of magnetohydrodynamics on orthogonal–curvilinear grids within a finite-volume framework. The scheme is based on a combination of central-upwind techniques for hyperbolic conservation laws and projection–evolution methods originally developed for Hamilton–Jacobi equations. The scheme is derived in semi-discrete form, and a full-fledged version is obtained by applying any stable and accurate solver for integration in time. The divergence-free condition of the magnetic field is a built-in property of the scheme by virtue of a constrained-transport ansatz for the induction equation. From the general formulation second-order accurate schemes for cylindrical grids and spherical grids are introduced in some more detail pointing out their potential importance in many applications. Special emphasis in this context is put to a treatment of the geometric axis implying severe complications because of the presence of coordinate singularities and associated grid degeneracy. An attempt to tackle these problems is presented. Numerical experiments illustrate the overall robustness and performance of the scheme for a small suite of tests.     

On a new defect of shock-capturing methods

This paper proposes an explanation and a cure (or avoidance) to the new defect found of Eulerian shock-capturing methods in “A note on the conservative schemes for the Euler equations” by Tang and Liu [H. Tang, Tiegang Liu, A note on the conservative schemes for the Euler equations, J. Comput. Phys. 218 (2006) 451–459]. The latter gives a numerical investigation using several popular high resolution conservative schemes applied to Riemann problems of inviscid, compressible, perfect gas flows in Eulerian and Lagrangian coordinates with an initial high density ratio as well as a high pressure ratio. The results show that these methods work very inefficiently when applied to such problems and may give inaccurate numerical results, especially in shock location (or speed), even with a very fine grid.We have found that in problems of this type a strong rarefaction wave (SRW) is present adjacent to a contact line. Godunov averaging over the wave then produces large errors which, when the wave is strong, also persist for a long time. The cumulative error is thus very large which violates the strength of the contact line adjacent to it which, in turn, affects the speed and hence the location of the shock on the other side of the contact. We confirm this numerically using a method based on the unified coordinates with the shock-adaptive Godunov scheme plus contact strength preserving. The method, when applied to the Examples 2.1 and 2.2 of Tang and Liu [H. Tang, Tiegang Liu, A note on the conservative schemes for the Euler equations, J. Comput. Phys. 218 (2006) 451–459], produces high quality results even for comparatively coarse grids.     

An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws

In this article we develop an improved version of the classical fifth-order weighted essentially non-oscillatory finite difference scheme of [G.S. Jiang, C.W. Shu, Efficient implementation of weighted ENO schemes, J. Comput. Phys. 126 (1996) 202–228] (WENO-JS) for hyperbolic conservation laws. Through the novel use of a linear combination of the low order smoothness indicators already present in the framework of WENO-JS, a new smoothness indicator of higher order is devised and new non-oscillatory weights are built, providing a new WENO scheme (WENO-Z) with less dissipation and higher resolution than the classical WENO. This new scheme generates solutions that are sharp as the ones of the mapped WENO scheme (WENO-M) of Henrick et al. [A.K. Henrick, T.D. Aslam, J.M. Powers, Mapped weighted essentially non-oscillatory schemes: achieving optimal order near critical points, J. Comput. Phys. 207 (2005) 542–567], however with a 25% reduction in CPU costs, since no mapping is necessary. We also provide a detailed analysis of the convergence of the WENO-Z scheme at critical points of smooth solutions and show that the solution enhancements of WENO-Z and WENO-M at problems with shocks comes from their ability to assign substantially larger weights to discontinuous stencils than the WENO-JS scheme, not from their superior order of convergence at critical points. Numerical solutions of the linear advection of discontinuous functions and nonlinear hyperbolic conservation laws as the one dimensional Euler equations with Riemann initial value problems, the Mach 3 shock–density wave interaction and the blastwave problems are compared with the ones generated by the WENO-JS and WENO-M schemes. The good performance of the WENO-Z scheme is also demonstrated in the simulation of two dimensional problems as the shock–vortex interaction and a Mach 4.46 Richtmyer–Meshkov Instability (RMI) modeled via the two dimensional Euler equations.