A high-order finite volume remapping scheme for nonuniform grids: The piecewise quartic method (PQM)
Introduction
Remapping is a crucial component of most arbitrary Lagrangian–Eulerian (ALE) algorithms used in computational fluid dynamics [13], [17], [14]. These algorithms involve a regridding step, whereby a new grid is generated based on some criteria, and a remapping step, whereby the variables are remapped from the old grid onto the new grid (Fig. 1). It is generally required that remapping be both conservative and monotonic in the sense that no new extrema should be created nor existing ones amplified. This is particularly important in applications where boundedness of some variables must be guaranteed.
The present study is motivated by the growing need to improve vertical coordinate systems in ocean general circulation models used for climate predictions. Over the last four decades or so, the vast majority of ocean models have used a single coordinate system in the vertical, usually aiming at a better representation of selected physical processes. These ocean models, however, have difficulties resolving physical processes for which they were not primarily designed. Hybrid coordinate ocean models have thus naturally emerged where the vertical grid is built by combining different coordinate systems in different regions [16], [2], [6], [1]. Due to the dynamical nature of the ocean, these hybrid coordinate systems are adapted in the course of the simulations, which means that vertical remapping is a key component. The accuracy of remapping is a major research issue in hybrid coordinate ocean models that prevents the hybrid framework from being truly convincing; presently, third-order reconstruction is used, at best, in these models. Hence, there is a need to explore higher-order representations.
An essential element of any remapping scheme lies in the reconstruction. Based on a set of cell averages on a given grid, the objective of a reconstruction scheme is to accurately and conservatively represent the underlying, real data with piecewise functions. Other constraints, such as monotonicity, may also apply. In this paper, we limit ourselves to piecewise polynomials but other functions can be employed, such as rational functions, e.g. [20]. Once the polynomial degree is chosen (for example, degree two leading to piecewise parabolas), a number of degrees of freedom must be determined that will produce a unique polynomial over the cell. There are basically two ways of improving a polynomial reconstruction method: First, the accuracy of the estimates for the pending degrees of freedom (e.g., edge values) may be improved. Second, when a monotonic reconstruction is required, improving the limiter may yield more accurate results. The latter domain of improvement, in particular, has drawn most of the attention in the case of the piecewise parabolic method (PPM) [4], [3], [12], [18], [11]. Besides these efforts, piecewise polynomial methods have witnessed almost no improvement regarding the edge estimates or the use of higher-degree polynomials. Only recently has the parabolic spline method been introduced but it also relies on polynomials of degree two [21]. Daru and Tenaud [5] clearly showed that successively high-order schemes lead to more accurate solutions for both limited and unlimited finite volume fluxes. The close relationship between their third-order method and PPM inspired us to seek improvement by using higher-order polynomials.
The objective of this paper is twofold: First, we introduce the piecewise quartic method (PQM) that uses piecewise polynomials of degree four and which, to our knowledge, has not been presented before. A limiter is devised that ensures monotonicity of any PQM-based remapping scheme. Second, a range of explicit and implicit schemes to estimate the edge values and slopes is investigated on the basis of accuracy and convergence analysis.
The main part (Section 2) of the paper focuses on reconstruction. We present PQM and the wide range of explicit and implicit schemes to estimate the edge values and slopes. The associated PPM and PQM schemes are evaluated in terms of accuracy. A limiter for PQM is also described in detail. The treatment of boundaries is covered in Section 3 and some comments on computational costs are made in Section 4. Conclusions are given in Section 5.
Section snippets
Piecewise polynomial reconstruction schemes
Given a nonuniform grid of cell widths and cell averages for , where N is the number of cells, the objective is to determine a piecewise polynomial reconstruction that accurately approximates the underlying data (Fig. 2). The polynomial over cell j is noted and, for convenience, use is made of a local coordinate such that the global coordinate x is given bywhere the global coordinates of the left and right cell interfaces are and
Treatment of boundaries
We now turn our attention on estimating the edge values and slopes at the boundaries. We limit our discussion to the following schemes: PPM , PPM , PQM and PQM , because they turn out to be the most effective. The concepts presented here are readily applicable to other schemes.
For PPM , edge values are estimated by using fourth-order polynomials spanning four cells. Thus, this scheme is not directly applicable to the first and last two edges for which another approach must
Computational cost
The relative computational costs of PPM and PQM schemes are now briefly investigated on the basis of remapping experiments in closed domains, as described in Fig. 11(b) (high-order boundary conditions are used). The errors and elapsed computational times are reported in Table 5 for unlimited and limited remapping experiments consisting of 20,000 cycles between a uniform 100-cell grid and a nonuniform 90-cell grid. Both the unlimited and limited versions of PQM schemes are more cost-effective
Conclusion
We have presented a hierarchy of one-dimensional high-order remapping schemes and investigated their performance with respect to accuracy and convergence rate. The schemes have also been compared based on remapping experiments in closed domains. We have introduced the new PQM scheme that is based on fifth-order accurate piecewise quartics. A limiter for this scheme has been fully described that never decreases the polynomial degree, except at the location of extrema where piecewise constants
Acknowledgments
We thank Bob Hallberg and Steve Griffies for comments regarding the topics covered in this paper. We also thank two anonymous reviewers for valuable suggestions that helped significantly improve the manuscript. Laurent White is supported by the ECCO2 project (http://ecco2.org/), “Estimating the Circulation and Climate of the Ocean, Phase II: High Resolution Global-Ocean and Sea-Ice Data Synthesis”, NASA award number NNG06GC28G. Laurent White is an honorary postdoctoral researcher with the
References (21)
- et al.
Formulation, implementation and examination of vertical coordinate choices in the Global Navy Coastal Ocean Model (NCOM)
Ocean Model
(2006) An oceanic general circulation model framed in hybrid isopycnic-Cartesian coordinates
Ocean Model
(2002)- et al.
The piecewise parabolic method (ppm) for gas-dynamical simulations
J. Comput. Phys.
(1984) - et al.
High order one-step monotonicity-preserving schemes for unsteady compressible flow calculations
J. Comput. Phys.
(2004) Evaluation of vertical coordinate and vertical mixing algorithms in the HYbrid-Coordinate Ocean Model (HYCOM)
Ocean Model
(2004)Marcelo H. Kobayashion a class of Padé finite volume methods
J. Comput. Phys.
(1999)- et al.
A finite volume formulation of compact central schemes on arbitrary structured grids
J. Comput. Phys.
(2004) - et al.
Second-order sign-preserving conservative interpolation (remapping) on general grids
J. Comput. Phys.
(2003) - et al.
A cell by cell anisotropic adaptive mesh ALE scheme for the numerical solution of the Euler equations
J. Comput. Phys.
(2007) - et al.
The repair paradigm and application to conservation laws
J. Comput. Phys.
(2004)
Cited by (38)
A vertically-Lagrangian, non-hydrostatic, multilayer model for multiscale free-surface flows
2020, Journal of Computational PhysicsCitation Excerpt :In multiple dimensions this is not a simple operation, in particular when conservation of the remapped quantities is required [42]. In the one-dimensional case we are interested in, this is much simpler and can be done using for example Piecewise Polynomial Reconstruction [43]. In our implementation we have used the quadratic piecewise remapping of the PPR library of Engwirda and Kelley (2016, https://github.com/dengwirda/PPR).
High-accurate and robust conservative remapping combining polynomial and hyperbolic tangent reconstructions
2020, Computers and FluidsCitation Excerpt :Such slope limiters may not always be trivial to design even if there exists a vast literature about slope/flux limiters, some of them are now known after their discoverers (van Leer [31–33], van Albada [34], Sweby [35], Barth-Jespersen [36], Venkatakrishnan [37], Koren [38,39], etc.) or their particularity (minmod [40], superbee [40], monotonized central [32], etc.). Relatively few attempts have been made to increase the remapping accuracy to more than second order via high order polynomial reconstructions, for several examples see [41–43]. Any attempt has to face the difficulty of designing a ‘limiter’ for non-smooth solution.
The numerics of hydrostatic structured-grid coastal ocean models: State of the art and future perspectives
2018, Ocean ModellingCitation Excerpt :Originally Hirt et al. (1974) defined the Arbitrary Lagrangian-Eulerian (ALE) method by a Lagrangian mesh movement (LVD step) followed by an optional readjustment of the mesh associated with grid-related advective transports (EVD step). The advection during the EVD step can be replaced by instantaneous “remapping” (monotone and conservative interpolation, see e.g. White and Adcroft, 2008) which relaxes the strict CFL constraint to a less severe Lipschitz type criteria of the form Δt max |∇u| < 1. Gräwe et al. (2015) demonstrated the clear advantages of adaptive coordinates compared to σ-coordinates in a realistic state-of-the-art model application with different dynamical regimes (tidal/non-tidal parts, seasonal thermocline / permanent halocline).