Central Essentially Non-oscillatory Research Articles

High‐order implicit time‐stepping with high‐order central essentially‐non‐oscillatory methods for unsteady three‐dimensional computational fluid dynamics simulations

AbstractThis article develops high‐order implicit time‐stepping methods combined with the fourth‐order central essentially‐non‐oscillatory (CENO) scheme for stiff three‐dimensional computational fluid dynamics problems having disparate characteristic time scales. Both aerodynamic and magnetohydrodynamic problems are considered on three‐dimensional multiblock body‐fitted grids with hexahedral cells. Several implicit time integration methods of third‐ and fourth‐order accuracy are considered, including the multistep backward differentiation formulas (BDF4), multistage explicitly singly diagonally implicit Runge‐Kutta (ESDIRK4), and Rosenbrock‐type methods (ROS34POW2). The resulting nonlinear algebraic system of equations is solved via a preconditioned Jacobian‐free inexact Newton–Krylov method with additive Schwarz preconditioning using block‐based incomplete LU decomposition. The performance of the high‐order implicit time‐stepping methods on smooth and stiff problems is compared with a standard fourth‐order explicit Runge‐Kutta (RK4) method. It is shown that the Rosenbrock methods, despite their advantage of only requiring the solution of linear systems, have significant drawbacks in terms of robustness issues for highly nonlinear compressible flows. The implicit BDF4 and ESDIRK4 methods are found to be much more efficient than the explicit fourth‐order RK4 method for a stiff resistive magnetohydrodynamic (MHD) problem discretized with the fourth‐order CENO method. When applied to the problem of vortex shedding governed by the Navier–Stokes equations, an A‐stable ESDIRK4 scheme proved to be the more robust and accurate implicit time‐marching scheme and was able to offer significant speedup compared with the RK4 method. Initial results are also shown for high‐order implicit time integration applied to two problems with discontinuities. The current study represents the first to achieve high‐order implicit time integration for MHD, enabling large time steps and substantial speedups for stiff MHD problems with high‐order accuracy, and it also represents the first to establish high‐order implicit time integration for high‐order CENO in space.

High-Order Finite-Volume Method with Block-Based AMR for Magnetohydrodynamics Flows

A high-order central essentially non-oscillatory (CENO) finite volume scheme combined with a block-based adaptive mesh refinement (AMR) algorithm is proposed for the solution of the ideal magnetohydrodynamics equations. The high-order CENO finite-volume scheme is implemented with fourth-order spatial accuracy within a flexible multi-block, body-fitted, hexahedral grid framework. An important feature of the high-order adaptive approach is that it allows for anisotropic refinement, which can lead to large computational savings when anisotropic flow features such as isolated propagating fronts and/or waves, shocks, shear surfaces, and current sheets are present in the flow. This approach is designed to handle complex multi-block grid configurations, including cubed-sphere grids, where some grid blocks may have degenerate edges or corners characterized by missing neighboring blocks. A procedure for building valid high-order reconstruction stencils, even at these degenerate block edges and corners, is proposed, taking into account anisotropic resolution changes in a systematic and general way. Furthermore, a non-uniform or heterogeneous block structure is used where the ghost cells of a block containing the solution content of neighboring blocks are stored directly at the resolution of the neighbors. A generalized Lagrange multiplier divergence correction technique is applied to achieve numerically divergence-free magnetic fields while preserving high-order accuracy on the anisotropic AMR grids. Parallel implementation and local grid adaptivity are achieved by using a hierarchical block-based domain partitioning strategy in which the connectivity and refinement history of grid blocks are tracked using a flexible binary tree data structure. Physics-based refinement criteria as well as the CENO smoothness indicator are both used for directing the mesh refinement. Numerical results, including solution-driven anisotropic refinement of cubed-sphere grids, are presented to demonstrate the accuracy and efficiency of the approach.

A high-order vertex-based central ENO finite-volume scheme for three-dimensional compressible flows

High-order discretization methods offer the potential to reduce the computational cost associated with modeling compressible flows. However, it is difficult to obtain accurate high-order discretizations of conservation laws that do not produce spurious oscillations near discontinuities, especially on multi-dimensional unstructured meshes. A novel, high-order, central essentially non-oscillatory (CENO) finite-volume method that does not have these difficulties is proposed for tetrahedral meshes. The proposed unstructured method is vertex-based, which differs from existing cell-based CENO formulations, and uses a hybrid reconstruction procedure that switches between two different solution representations. It applies a high-order k-exact reconstruction in smooth regions and a limited linear reconstruction when discontinuities are encountered. Both reconstructions use a single, central stencil for all variables, making the application of CENO to arbitrary unstructured meshes relatively straightforward. The new approach was applied to the conservation equations governing compressible flows and assessed in terms of accuracy and computational cost. For all problems considered, which included various function reconstructions and idealized flows, CENO demonstrated excellent reliability and robustness. Up to fifth-order accuracy was achieved in smooth regions and essentially non-oscillatory solutions were obtained near discontinuities. The high-order schemes were also more computationally efficient for high-accuracy solutions, i.e., they took less wall time than the lower-order schemes to achieve a desired level of error. In one particular case, it took a factor of 24 less wall-time to obtain a given level of error with the fourth-order CENO scheme than to obtain the same error with the second-order scheme.

A High-Order Central ENO Finite-Volume Scheme for Three-Dimensional Low-Speed Viscous Flows on Unstructured Mesh

AbstractHigh-order discretization techniques offer the potential to significantly reduce the computational costs necessary to obtain accurate predictions when compared to lower-order methods. However, efficient and universally-applicable high-order discretizations remain somewhat illusive, especially for more arbitrary unstructured meshes and for incompressible/low-speed flows. A novel, high-order, central essentially non-oscillatory (CENO), cell-centered, finite-volume scheme is proposed for the solution of the conservation equations of viscous, incompressible flows on three-dimensional unstructured meshes. Similar to finite element methods, coordinate transformations are used to maintain the scheme’s order of accuracy even when dealing with arbitrarily-shaped cells having non-planar faces. The proposed scheme is applied to the pseudo-compressibility formulation of the steady and unsteady Navier-Stokes equations and the resulting discretized equations are solved with a parallel implicit Newton-Krylov algorithm. For unsteady flows, a dual-time stepping approach is adopted and the resulting temporal derivatives are discretized using the family of high-order backward difference formulas (BDF). The proposed finite-volume scheme for fully unstructured mesh is demonstrated to provide both fast and accurate solutions for steady and unsteady viscous flows.

High-order central ENO finite-volume scheme for hyperbolic conservation laws on three-dimensional cubed-sphere grids

A fourth-order accurate finite-volume scheme for hyperbolic conservation laws on three-dimensional (3D) cubed-sphere grids is described. The approach is based on a central essentially non-oscillatory (CENO) finite-volume method that was recently introduced for two-dimensional compressible flows and is extended to 3D geometries with structured hexahedral grids. Cubed-sphere grids feature hexahedral cells with nonplanar cell surfaces, which are handled with high-order accuracy using trilinear geometry representations in the proposed approach. Varying stencil sizes and slope discontinuities in grid lines occur at the boundaries and corners of the six sectors of the cubed-sphere grid where the grid topology is unstructured, and these difficulties are handled naturally with high-order accuracy by the multidimensional least-squares based 3D CENO reconstruction with overdetermined stencils. A rotation-based mechanism is introduced to automatically select appropriate smaller stencils at degenerate block boundaries, where fewer ghost cells are available and the grid topology changes, requiring stencils to be modified. Combining these building blocks results in a finite-volume discretization for conservation laws on 3D cubed-sphere grids that is uniformly high-order accurate in all three grid directions. While solution-adaptivity is natural in the multi-block setting of our code, high-order accurate adaptive refinement on cubed-sphere grids is not pursued in this paper. The 3D CENO scheme is an accurate and robust solution method for hyperbolic conservation laws on general hexahedral grids that is attractive because it is inherently multidimensional by employing a K-exact overdetermined reconstruction scheme, and it avoids the complexity of considering multiple non-central stencil configurations that characterizes traditional ENO schemes. Extensive numerical tests demonstrate fourth-order convergence for stationary and time-dependent Euler and magnetohydrodynamic flows on cubed-sphere grids, and robustness against spurious oscillations at 3D shocks. Performance tests illustrate efficiency gains that can be potentially achieved using fourth-order schemes as compared to second-order methods for the same error level. Applications on extended cubed-sphere grids incorporating a seventh root block that discretizes the interior of the inner sphere demonstrate the versatility of the spatial discretization method.

High-order solution-adaptive central essentially non-oscillatory (CENO) method for viscous flows

A high-order, central, essentially non-oscillatory (CENO), finite-volume scheme in combination with a block-based adaptive mesh refinement (AMR) algorithm is proposed for solution of the Navier-Sto...

High-order solution-adaptive central essentially non-oscillatory (CENO) method for viscous flows

A high-order, central, essentially non-oscillatory (CENO), finite-volume scheme in combination with a block-based adaptive mesh refinement (AMR) algorithm is proposed for solution of the Navier–Stokes equations on body-fitted multi-block mesh. In contrast to other ENO schemes which require reconstruction on multiple stencils, the proposed CENO method uses a hybrid reconstruction approach based on a fixed central stencil. This feature is crucial to avoiding the complexities associated with multiple stencils of ENO schemes, providing high-order accuracy at relatively lower computational cost as well as being very well suited for extension to unstructured meshes. The spatial discretization of the inviscid (hyperbolic) fluxes combines an unlimited high-order k-exact least-squares reconstruction technique following from the optimal central stencil with a monotonicity-preserving, limited, linear, reconstruction algorithm. This hybrid reconstruction procedure retains the unlimited high-order k-exact reconstruction for cells in which the solution is fully resolved and reverts to the limited lower-order counterpart for cells with under-resolved/discontinuous solution content. Switching in the hybrid procedure is determined by a smoothness indicator. The high-order viscous (elliptic) fluxes are computed to the same order of accuracy as the hyperbolic fluxes based on a k-order accurate cell interface gradient derived from the unlimited, cell-centred, reconstruction. A somewhat novel h-refinement criterion based on the solution smoothness indicator is used to direct the steady and unsteady mesh adaptation. The proposed numerical procedure is thoroughly analyzed for advection–diffusion problems characterized by the full range of Péclet numbers, and its predictive capabilities are also demonstrated for several inviscid and laminar flows. The ability of the scheme to accurately represent solutions with smooth extrema and yet robustly handle under-resolved and/or non-smooth solution content (i.e., shocks and other discontinuities) is shown. Moreover, the ability to perform mesh refinement in regions of under-resolved and/or non-smooth solution content is also demonstrated.

High-order central ENO finite-volume scheme for ideal MHD

A high-order accurate finite-volume scheme for the compressible ideal magnetohydrodynamics (MHD) equations is proposed. The high-order MHD scheme is based on a central essentially non-oscillatory (CENO) method combined with the generalized Lagrange multiplier divergence cleaning method for MHD. The CENO method uses k-exact multidimensional reconstruction together with a monotonicity procedure that switches from a high-order reconstruction to a limited low-order reconstruction in regions of discontinuous or under-resolved solution content. Both reconstructions are performed on central stencils, and the switching procedure is based on a smoothness indicator. The proposed high-order accurate MHD scheme can be used on general polygonal grids. A highly sophisticated parallel implementation of the scheme is described that is fourth-order accurate on two-dimensional dynamically-adaptive body-fitted structured grids. The hierarchical multi-block body-fitted grid permits grid lines to conform to curved boundaries. High-order accuracy is maintained at curved domain boundaries by employing high-order spline representations and constraints at the Gauss quadrature points for flux integration. Detailed numerical results demonstrate high-order convergence for smooth flows and robustness against oscillations for problems with shocks. A new MHD extension of the well-known Shu–Osher test problem is proposed to test the ability of the high-order MHD scheme to resolve small-scale flow features in the presence of shocks. The dynamic mesh adaptation capabilities of the approach are demonstrated using adaptive time-dependent simulations of the Orszag–Tang vortex problem with high-order accuracy and unprecedented effective resolution.