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.     

Balanced finite volume WENO and central WENO schemes for the shallow water and the open-channel flow equations

The goal of this work is to extend finite volume WENO and central WENO schemes to the hyperbolic balance laws with geometrical source term and spatially variable flux function. In particular, we apply proposed schemes to the shallow water and the open-channel flow equations where the source term depends on the channel geometry. For obtaining stable numerical schemes that are free of spurious oscillations, it becomes crucial to use the decomposed source term evaluation, which maintains the balancing between the flux gradient and the source term. In addition, the open-channel flow equations contain spatially variable flux function. The appropriate definitions of the terms that arise in the source term decomposition, in combination with the Roe approximate Riemann solver that includes the spatial derivative of the flux function, lead to the finite volume WENO scheme that satisfies the exact conservation property – the property of preserving the quiescent flow exactly. When the central WENO schemes are applied, additional reformulations are introduced for the transition from the staggered values to the nonstaggered ones and vice versa by using the WENO reconstruction procedure. The proposed central WENO schemes also preserve the quiescent flow, but only in prismatic channels. In various test problems the obtained balanced schemes show improvements in comparison with the standard versions of the proposed type schemes, as well as with some other first- and second-order numerical schemes.     

Generalized finite compact difference scheme for shock/complex flowfield interaction

A class of generalized high order finite compact difference schemes is proposed for shock/vortex, shock/boundary layer interaction problems. The finite compact difference scheme takes the region between two shocks as a compact stencil. The high order WENO fluxes on shock stencils are used as the internal boundary fluxes for the compact scheme. A lemma based on the property of smoothness estimators on a 5-points stencil is given to detect the shock position. There is no free parameter introduced to switch the compact scheme and the WENO scheme. Some numerical experiments are given and they demonstrate that the present scheme has low dissipation due to the compact central differencing scheme used in the smooth regions.     

Multi-domain Fourier-continuation/WENO hybrid solver for conservation laws

We introduce a multi-domain Fourier-continuation/WENO hybrid method (FC–WENO) that enables high-order and non-oscillatory solution of systems of nonlinear conservation laws, and which enjoys essentially dispersionless, spectral character away from discontinuities, as well as mild CFL constraints (comparable to those of finite difference methods). The hybrid scheme employs the expensive, shock-capturing WENO method in small regions containing discontinuities and the efficient FC method in the rest of the computational domain, yielding a highly effective overall scheme for applications with a mix of discontinuities and complex smooth structures. The smooth and discontinuous solution regions are distinguished using the multi-resolution procedure of Harten [J. Comput. Phys. 115 (1994) 319–338]. We consider WENO schemes of formal orders five and nine and a FC method of order five. The accuracy, stability and efficiency of the new hybrid method for conservation laws is investigated for problems with both smooth and non-smooth solutions. In the latter case, we solve the Euler equations for gas dynamics for the standard test case of a Mach three shock wave interacting with an entropy wave, as well as a shock wave (with Mach 1.25, three or six) interacting with a very small entropy wave and evaluate the efficiency of the hybrid FC–WENO method as compared to a purely WENO-based approach as well as alternative hybrid based techniques. We demonstrate considerable computational advantages of the new FC-based method, suggesting a potential of an order of magnitude acceleration over alternatives when extended to fully three-dimensional problems.     

A weighted-integral based scheme

A weighted-integral based scheme (WIBS) and a weighted essentially non-oscillatory (WENO)-WIBS are constructed, where the integral of the unknown function with a set of linearly independent test functions are recorded on every cell. The time evolutions of these recordings are computed with TVD Runge–Kutta method. At the boundary of every two cells, the function values are interpolated from the recordings of the neighboring cells to calculate flux and volumetric integral in the weak form.     

A hybrid numerical simulation of isotropic compressible turbulence

A novel hybrid numerical scheme with built-in hyperviscosity has been developed to address the accuracy and numerical instability in numerical simulation of isotropic compressible turbulence in a periodic domain at high turbulent Mach number. The hybrid scheme utilizes a 7th-order WENO (Weighted Essentially Non-Oscillatory) scheme for highly compressive regions (i.e., shocklet regions) and an 8th-order compact central finite difference scheme for smooth regions outside shocklets. A flux-based conservative and formally consistent formulation is developed to optimize the connection between the two schemes at the interface and to achieve a higher computational efficiency. In addition, a novel numerical hyperviscosity formulation is proposed within the context of compact finite difference scheme for the smooth regions to improve numerical stability of the hybrid method. A thorough and insightful analysis of the hyperviscosity formulation in both Fourier space and physical space is presented to show the effectiveness of the formulation in improving numerical stability, without compromising the accuracy of the hybrid method. A conservative implementation of the hyperviscosity formulation is also developed. Combining the analysis and test simulations, we have also developed a criterion to guide the specification of a numerical hyperviscosity coefficient (the only adjustable coefficient in the formulation). A series of test simulations are used to demonstrate the accuracy and numerical stability of the scheme for both decaying and forced compressible turbulence. Preliminary results for a high-resolution simulation at turbulent Mach number of 1.08 are shown. The sensitivity of the simulated flow to the detail of thermal forcing method is also briefly discussed.     

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.     

A systematic methodology for constructing high-order energy stable WENO schemes

A third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme developed by the authors of this paper [N.K. Yamaleev, M.H. Carpenter, Third-order energy stable WENO scheme, J. Comput. Phys. 228 (2009) 3025–3047] was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables “energy stable” modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes including one sixth-order central scheme; ESWENO schemes up to eighth order are presented in the Appendix. We also develop new weight functions and derive constraints on their parameters, which provide consistency, much faster convergence of the high-order ESWENO schemes to their underlying linear schemes for smooth solutions with arbitrary number of vanishing derivatives, and better resolution near strong discontinuities than the conventional counterparts.     

基于WENO重构的高阶移动网格动理学格式

提出一种基于WENO重构的高阶(至少三阶)移动网格动理学格式.利用流体力学方程的积分形式得到移动网格上离散格式,再利用自适应移动网格方法移动网格,进而得到网格速度,利用WENO重构得到高阶插值多项式,最后使用时间方向上精确的动理学数值方法构造数值通量,得到移动网格单元上新的物理量.数值实验表明这种格式同时具有高精度、高分辨率的特点.     

A stable and high-order accurate conjugate heat transfer problem

This paper analyzes well-posedness and stability of a conjugate heat transfer problem in one space dimension. We study a model problem for heat transfer between a fluid and a solid. The energy method is used to derive boundary and interface conditions that make the continuous problem well-posed and the semi-discrete problem stable. The numerical scheme is implemented using 2nd-, 3rd- and 4th-order finite difference operators on Summation-By-Parts (SBP) form. The boundary and interface conditions are implemented weakly. We investigate the spectrum of the spatial discretization to determine which type of coupling that gives attractive convergence properties. The rate of convergence is verified using the method of manufactured solutions.     

A new family of high-order compact upwind difference schemes with good spectral resolution

This paper presents a new family of high-order compact upwind difference schemes. Unknowns included in the proposed schemes are not only the values of the function but also those of its first and higher derivatives. Derivative terms in the schemes appear only on the upwind side of the stencil. One can calculate all the first derivatives exactly as one solves explicit schemes when the boundary conditions of the problem are non-periodic. When the proposed schemes are applied to periodic problems, only periodic bi-diagonal matrix inversions or periodic block-bi-diagonal matrix inversions are required. Resolution optimization is used to enhance the spectral representation of the first derivative, and this produces a scheme with the highest spectral accuracy among all known compact schemes. For non-periodic boundary conditions, boundary schemes constructed in virtue of the assistant scheme make the schemes not only possess stability for any selective length scale on every point in the computational domain but also satisfy the principle of optimal resolution. Also, an improved shock-capturing method is developed. Finally, both the effectiveness of the new hybrid method and the accuracy of the proposed schemes are verified by executing four benchmark test cases.     

A high-order cell-centered Lagrangian scheme for compressible fluid flows in two-dimensional cylindrical geometry

The goal of this paper is to present high-order cell-centered schemes for solving the equations of Lagrangian gas dynamics written in cylindrical geometry. A node-based discretization of the numerical fluxes is obtained through the computation of the time rate of change of the cell volume. It allows to derive finite volume numerical schemes that are compatible with the geometric conservation law (GCL). Two discretizations of the momentum equations are proposed depending on the form of the discrete gradient operator. The first one corresponds to the control volume scheme while the second one corresponds to the so-called area-weighted scheme. Both formulations share the same discretization for the total energy equation. In both schemes, fluxes are computed using the same nodal solver which can be viewed as a two-dimensional extension of an approximate Riemann solver. The control volume scheme is conservative for momentum, total energy and satisfies a local entropy inequality in its first-order semi-discrete form. However, it does not preserve spherical symmetry. On the other hand, the area-weighted scheme is conservative for total energy and preserves spherical symmetry for one-dimensional spherical flow on equi-angular polar grid. The two-dimensional high-order extensions of these two schemes are constructed employing the generalized Riemann problem (GRP) in the acoustic approximation. Many numerical tests are presented in order to assess these new schemes. The results obtained for various representative configurations of one and two-dimensional compressible fluid flows show the robustness and the accuracy of our new schemes.     

Perfect plane-wave source for a high-order symplectic finite-difference time-domain scheme

The method of splitting a plane-wave finite-difference time-domain (SP-FDTD) algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field (TF/SF) formulation of high-order symplectic finite-difference time-domain (SFDTD) scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain (FDTD), the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.     

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.     

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.     

A hybrid multilevel method for high-order discretization of the Euler equations on unstructured meshes

Higher order discretization has not been widely successful in industrial applications to compressible flow simulation. Among several reasons for this, one may identify the lack of tailor-suited, best-practice relaxation techniques that compare favorably to highly tuned lower order methods, such as finite-volume schemes. In this paper we investigate solution algorithms in conjunction with high-order Spectral Difference discretization for the Euler equations, using such techniques as multigrid and matrix-free implicit relaxation methods. In particular we present a novel hybrid multilevel relaxation method that combines (optionally matrix-free) implicit relaxation techniques with explicit multistage smoothing using geometric multigrid. Furthermore, we discuss efficient implementation of these concepts using such tools as automatic differentiation.     

A two-dimensional unstructured cell-centered multi-material ALE scheme using VOF interface reconstruction

We present a new cell-centered multi-material arbitrary Lagrangian–Eulerian (ALE) scheme to solve the compressible gas dynamics equations on two-dimensional unstructured grid. Our ALE method is of the explicit time-marching Lagrange plus remap type. Namely, it involves the following three phases: a Lagrangian phase wherein the flow is advanced using a cell-centered scheme; a rezone phase in which the nodes of the computational grid are moved to more optimal positions; a cell-centered remap phase which consists of interpolating conservatively the Lagrangian solution onto the rezoned grid. The multi-material modeling utilizes either concentration equations for miscible fluids or the Volume Of Fluid (VOF) capability with interface reconstruction for immiscible fluids. The main original feature of this ALE scheme lies in the introduction of a new mesh relaxation procedure which keeps the rezoned grid as close as possible to the Lagrangian one. In this formalism, the rezoned grid is defined as a convex combination between the Lagrangian grid and the grid resulting from condition number smoothing. This convex combination is constructed through the use of a scalar parameter which is a scalar function of the invariants of the Cauchy–Green tensor over the Lagrangian phase. Regarding the cell-centered remap phase, we employ two classical methods based on a partition of the rezoned cell in terms of its overlap with the Lagrangian cells. The first one is a simplified swept face-based method whereas the second one is a cell-intersection-based method. Our multi-material ALE methodology is assessed through several demanding two-dimensional tests. The corresponding numerical results provide a clear evidence of the robustness and the accuracy of this new scheme.     

A prediction-correction scheme for microchannel milling using femtosecond laser

In this paper, a prediction-correction scheme is proposed to online measure and regulate the milling depth of microchannel using an indicator of laser triggered plasma. Firstly, a prediction model, with respect to the laser fluence and feedrate, is established with several calibration tests using the least square fitting method. It is utilized to change the focal position of objective to track the depth evolution of newly generated surface. Meanwhile, a scanning path for every milling layer with an offset in Z-axis at the beginning and the end of the trajectory, is developed to drive the plasma brightness periodically changing. Then, the milling depth could be obtained when the brightness reaches to the maximum value. By doing so, an online measurement method is presented to estimate the milling depth using the trend of plasma brightness. Furthermore, a correction model is developed to iteratively adjust the feedrate with the online estimated depth. Therefore, the microchannel milling process could be monitored and controlled in a closed-loop manner, in order to accurately regulate the milling depth. Finally, an online measurement and closed-loop microchannel milling is carried out on the self-developed micro-machining center. The effectiveness and correctness of the proposed method are verified by comparing the estimated depth with the actually measured results.     

Comparison of the generalized Riemann solver and the gas-kinetic scheme for inviscid compressible flow simulations

The generalized Riemann problem (GRP) scheme for the Euler equations and gas-kinetic scheme (GKS) for the Boltzmann equation are two high resolution shock capturing schemes for fluid simulations. The difference is that one is based on the characteristics of the inviscid Euler equations and their wave interactions, and the other is based on the particle transport and collisions. The similarity between them is that both methods can use identical MUSCL-type initial reconstructions around a cell interface, and the spatial slopes on both sides of a cell interface involve in the gas evolution process and the construction of a time-dependent flux function. Although both methods have been applied successfully to the inviscid compressible flow computations, their performances have never been compared. Since both methods use the same initial reconstruction, any difference is solely coming from different underlying mechanism in their flux evaluation. Therefore, such a comparison is important to help us to understand the correspondence between physical modeling and numerical performances. Since GRP is so faithfully solving the inviscid Euler equations, the comparison can be also used to show the validity of solving the Euler equations itself. The numerical comparison shows that the GRP exhibits a slightly better computational efficiency, and has comparable accuracy with GKS for the Euler solutions in 1D case, but the GKS is more robust than GRP. For the 2D high Mach number flow simulations, the GKS is absent from the shock instability and converges to the steady state solutions faster than the GRP. The GRP has carbuncle phenomena, likes a cloud hanging over exact Riemann solvers. The GRP and GKS use different physical processes to describe the flow motion starting from a discontinuity. One is based on the assumption of equilibrium state with infinite number of particle collisions, and the other starts from the non-equilibrium free transport process to evolve into an equilibrium one through particle collisions. The different mechanism in the flux evaluation deviates their numerical performance. Through this study, we may conclude scientifically that it may NOT be valid to use the Euler equations as governing equations to construct numerical fluxes in a discretized space with limited cell resolution. To adapt the Navier–Stokes (NS) equations is NOT valid either because the NS equations describe the flow behavior on the hydrodynamic scale and have no any corresponding physics starting from a discontinuity. This fact alludes to the consistency of the Euler and Navier–Stokes equations with the continuum assumption and the necessity of a direct modeling of the physical process in the discretized space in the construction of numerical scheme when modeling very high Mach number flows. The development of numerical algorithm is similar to the modeling process in deriving the governing equations, but the control volume here cannot be shrunk to zero.