High-order Finite Volume Research Articles

Comparison between a priori and a posteriori slope limiters for high-order finite volume schemes

Comparison between a priori and a posteriori slope limiters for high-order finite volume schemes

A fourth-order compact finite volume method on unstructured grids for simulation of two-dimensional incompressible flow

This paper introduces a novel and efficient compact reconstruction procedure for high-order finite volume methods applied to unstructured grids. In this procedure, we establish a set of constitutive relations that ensure the continuity of the reconstruction polynomial and its normal derivatives between adjacent elements at control points. The paper delves into the details of the fourth-order compact reconstruction method specifically designed for two-dimensional triangular grids. This method can be considered an extension of the one-dimensional compact scheme and cubic spline interpolation methods to two-dimensional triangular unstructured grids. In the cubic polynomial reconstruction, a two-dimensional cubic polynomial is reconstructed using the cell averages of elements in the standard stencil and the function values at control points. The determination of function values at control points is achieved by adhering to the constraint of normal derivative continuity. This reconstruction is solved using a relaxation iteration method and has been verified to be solvable for triangular grids. Compared to other fourth-order implicit compact polynomial reconstructions, our method requires solving a smaller number of unknowns, which means less computational cost in reconstruction. By applying this method to solve the Poisson equation, the linear convection equation, and incompressible flow benchmarks, it demonstrates that the proposed method exhibits the expected high-order accuracy and performs well in incompressible flow problems.

High-order finite volume methods for solving compressible multicomponent flows

High-order finite volume methods for solving compressible multicomponent flows

High-order finite volume method for solving compressible multicomponent flows with Mie–Grüneisen equation of state

In this paper, we propose a new high-order finite volume method for solving the multicomponent fluids problem with Mie–Grüneisen EOS. Firstly, based on the cell averages of conservative variables, we develop a procedure to reconstruct the cell averages of the primitive variables in a high-order manner. Secondly, the high-order reconstructions employed in computing numerical fluxes are implemented in a characteristic-wise manner to reduce numerical oscillations as much as possible and obtain high-resolution results. Thirdly, advection equation within the governing system is rewritten in a conservative form with a source term to enhance the scheme’s performance. We utilize integration by parts and high-order numerical integration techniques to handle the source terms. Finally, all variables are evolved by using Runge–Kutta time discretization. All steps are carefully designed to maintain the equilibrium of pressure and velocity for the interface-only problem, which is crucial in designing a high-resolution scheme and adapting to more complex multicomponent problems. We have performed extensive numerical tests for both one- and two-dimensional problems to verify our scheme’s high resolution and accuracy.

A shock capturing artificial viscosity scheme in consistent with the compact high-order finite volume methods

This paper presents a shock capturing artificial viscosity scheme for the compact high-order finite volume methods in terms of the variational reconstructions on unstructured grids. The key for the design of the present artificial viscosity is the smoothness indicator, which is based on the concept of interfacial jump integration, measuring the discontinuities of the reconstruction polynomial and its spatial derivatives across a cell interface. Since the variational reconstruction is carried out by minimizing the functional in terms of the interfacial jump integration, the present smoothness indicator gives a discretization-consistent measurement of the smoothness of the flow fields that is sufficiently large in the region near discontinuities, and is in the same order of magnitude as the spatial truncation error of the finite volume scheme in smooth regions. These properties ensure that the newly developed artificial viscosity scheme has the problem-independent capability to suppress non-physical oscillations near discontinuities and preserve the theoretical order of accuracy for smooth flow. The shock capturing capability of the proposed artificial viscosity scheme has been demonstrated by a number of numerical examples confirming its essentially non-oscillatory and high-resolution properties. Additionally, the proposed artificial viscosity scheme exhibits higher computational efficiency than the approach based on a traditional limiter.

A robust and well-balanced finite volume solver for investigating the effects of tides on water renewal timescale in the Nador lagoon, Morocco

Coastal lagoons, particularly those not well-connected to the sea, are highly susceptible to continuous water retention, adversely affecting water quality and ecosystems. Understanding and quantifying the processes governing water renewal in such semi-enclosed regions is crucial. Using the principles of constituent-oriented age and residence time theory (www.climate.be/cart), we calculated timescales to explain the water renewal process in semi-enclosed areas. To improve our understanding and define the dynamics of the Nador lagoon, we use two distinct time scales: residence time and exposure time. Residence time indicates the length of time a parcel of water leaves the region of interest for the first time, while exposure time represents the overall length of time a parcel of water remains in that region, including any subsequent inflows. The modeling system adopts an Eulerian approach, including two interlinked model components: the shallow-water equations incorporating bottom friction, eddy viscosity, wind shear stresses, momentum diffusion, and Coriolis forces for modeling hydrodynamics, and a transport-diffusion equation utilized to simulate the advection and diffusion of the passive tracer concentration. The entire system is integrated into a high-order finite volume solver that employs unstructured meshes, upwind numerical fluxes, and slope limiters to achieve precise resolution of steep bathymetric gradients that may arise in the approximate solution. The scheme is non-oscillatory and preserves the conservation properties. In this research, we explore various situations related to the connection between the Mediterranean Sea and the Nador lagoon. Specifically, we suggest novel entry points that have the potential to significantly decrease the water residence within the lagoon.

Exploring the potential of TENO and WENO schemes for simulating under-resolved turbulent flows in the atmosphere using Euler equations

This paper focuses on the design and analysis of very high-order finite volume methods for the computation of simplified meso- and micro-scale atmospheric flows. In a dry atmosphere, these flows can be represented by the Euler equations with a gravitational source term. Two different approaches are considered here. While one of the approaches is fully conservative for the total energy, the other is formulated in a non-conservative form. The main focus of the paper is to analyze the performance of such models in combination with the traditional WENO reconstruction and the novel TENO reconstruction by examining the spectral properties of these reconstruction methods. The overarching goal is to determine whether the combination of these models and numerical schemes can be used to build an implicit Large Eddy Simulation framework, shedding light on their potential advantages or limitations in representing under-resolved atmospheric processes in the meso- and micro-scales.

Quinpi: Integrating Stiff Hyperbolic Systems with Implicit High Order Finite Volume Schemes

Quinpi: Integrating Stiff Hyperbolic Systems with Implicit High Order Finite Volume Schemes

A third-order compact finite volume scheme on unstructured grid for fluid flows

The large reconstruction stencil is one of the biggest obstacles in design and implementation of high-order finite volume schemes on unstructured grid. To tackle this problem, a third-order compact reconstruction scheme on unstructured grid is developed on the basis of compact least-squares (CLS) finite volume method. In building the reconstruction polynomial, the general boundary condition is adopted to ensure the unified third-order accuracy in both the internal and boundary cells. The reconstruction polynomial is solved in a semi-implicit sense to ensure high computational efficiency without solving the global linear algebraic equations. The theoretical accuracy order and stability limit are firstly verified by the Fourier analysis. The numerical accuracy and efficiency are then verified and validated with several incompressible and compressible flow cases. This work provides a reference to design and analysis of high-order compact finite volume schemes on unstructured grid with unified accuracy.

KFVM-WENO: A High-order Accurate Kernel-based Finite Volume Method for Compressible Hydrodynamics

This paper presents a fully multidimensional kernel-based reconstruction scheme for finite volume methods applied to systems of hyperbolic conservation laws, with a particular emphasis on the compressible Euler equations. Nonoscillatory reconstruction is achieved through an adaptive-order weighted essentially nonoscillatory (WENO) method cast into a form suited to multidimensional reconstruction. A kernel-based approach inspired by radial basis functions and Gaussian process modeling, which we call kernel-based finite volume method with WENO, is presented here. This approach allows the creation of a scheme of arbitrary order of accuracy with simply defined multidimensional stencils and substencils. Furthermore, the fully multidimensional nature of the reconstruction allows for a more straightforward extension to higher spatial dimensions and removes the need for complicated boundary conditions on intermediate quantities in modified dimension-by-dimension methods. In addition, a new simple yet effective set of reconstruction variables is introduced, which could be useful in existing schemes with little modification. The proposed scheme is applied to a suite of stringent and informative benchmark problems to demonstrate its efficacy and utility. A highly parallel multi-GPU implementation using Kokkos and the message-passing interface is also provided.

Combination of Advanced Actuator Line/Disk Model and High-Order Unstructured Finite Volume Solver for Helicopter Rotors

In the research field of rotorcraft aerodynamics, there are two fundamental challenges: resolving the complex vortex structures in rotor wakes and representing the moving rotor blades in the ambient airflow. In this paper, we address the first challenge by utilizing a third-order unstructured finite volume solver, which exhibits lower numerical dissipation than its second-order counterpart. This allows for sufficient resolution of small vortex structures on relatively coarse meshes. With this flow solver, the second challenge is addressed by modeling each rotor as an actuator disk (i.e., the actuator disk model (ADM)) or modeling each blade as an actuator line (i.e., the actuator line model (ALM)). Both of the two models are equipped with an improved tip loss correction, which is introduced in detail in the methodology section. In the section of numerical experiments, the numerical convergence properties of the two types of solvers have been compared in the case of two-dimensional infinite wing. In addition, the relationship between the ALM and the lifting line theory is discussed in the cases of fixed-wing calculations. Another goal of these cases is to validate the tip loss correction presented. The validation of the ALM/ADM and comparisons of computational efficiency are also demonstrated in simulations involving both hover and forward flight rotors. It was found that the combination of the third-order finite volume solver and the ALM/ADM with the improved tip loss correction presents an efficient way of performing the aerodynamic analysis of rotor-induced downwash flow.

A general positivity-preserving algorithm for implicit high-order finite volume schemes solving the Euler and Navier-Stokes equations

This paper presents a general positivity-preserving algorithm for implicit high-order finite volume schemes that solve compressible Euler and Navier-Stokes equations to ensure the positivity of density and internal energy (or pressure). Previous positivity-preserving algorithms are mainly based on the slope limiting or flux limiting technique, which rely on the existence of low-order positivity-preserving schemes. This dependency poses serious restrictions on extending these algorithms to temporally implicit schemes since it is difficult to know if a low-order implicit scheme is positivity-preserving. In the present paper, a new positivity-preserving algorithm is proposed in terms of the flux correction technique. And the factors of the flux correction are determined by a residual correction procedure. For a finite volume scheme that is capable of achieving a converged solution, we show that the correction factors are in the order of unity with additional high-order terms corresponding to the spatial and temporal rates of convergence. Therefore, the proposed positivity-preserving algorithm is accuracy-reserving and asymptotically consistent. The notable advantage of this method is that it does not rely on the existence of low-order positivity-preserving baseline schemes. Therefore, it can be applied to the implicit schemes solving Euler and especially Navier-Stokes equations. In the present paper, the proposed technique is applied to an implicit dual time-stepping finite volume scheme with temporal second-order and spatial high-order accuracy. The present positivity-preserving algorithm is implemented in an iterative manner to ensure that the dual time-stepping iteration will converge to the positivity-preserving solution. Another similar correction technique is also proposed to ensure that the solution remains positivity-preserving at each sub-iteration. Numerical results demonstrate that the proposed algorithm preserves positive density and internal energy in all test cases and significantly improves the robustness of the numerical schemes.

High-order genuinely multidimensional finite volume methods via kernel-based WENO

In this paper a family of fully multidimensional kernel-based reconstruction schemes for use in finite volume methods (FVMs) will be presented. These methods are intended for use in shock dominated problems, and stability is achieved through a suitable adaptation of the Adaptive Order Weighted Essentially Non-Oscillatory (WENO-AO) method to the proposed kernel-based reconstruction schemes. There are a number of key difficulties in the design of high-order finite volume schemes which will be discussed and addressed. High (4th and 6th) order convergence will be demonstrated on smooth exact solutions of the ideal MHD equations. The very same scheme will then be applied to extremely stringent astrophysical benchmark problems.

High-order finite volume multi-resolution WENO schemes with adaptive linear weights on triangular meshes

This paper presents high-order finite volume multi-resolution weighted essentially non-oscillatory schemes with adaptive linear weights to solve hyperbolic conservation laws on triangular meshes. They are abbreviated as the ALW-MR-WENO schemes. The novel third-order, fourth-order, and fifth-order ALW-MR-WENO schemes are designed by applying two unequal-sized hierarchical central stencils in comparison to the classical WENO schemes which utilize many equal-sized upwind biased/central stencils. With the application of one simple condition, only two linear weights are automatically adjusted to be positive values on condition that their summation is one. The novel finite volume ALW-MR-WENO schemes could maintain the designed order of accuracy in smooth areas and reduce to the first-order accuracy so as to keep essentially non-oscillatory properties around strong discontinuities. So it is the first time that any high-order WENO schemes with the application of only two unequal-sized stencils are obtained on triangular meshes. And the major benefits are their efficiency, compactness, and simplicity in large scale engineering applications on unstructured meshes. Finally, several tests are used to indicate the effectiveness of these new finite volume WENO schemes.

A New Discretely Divergence-Free Positivity-Preserving High-Order Finite Volume Method for Ideal MHD Equations

A New Discretely Divergence-Free Positivity-Preserving High-Order Finite Volume Method for Ideal MHD Equations

Improved higher-order finite volume scheme and its convergence analysis for collisional breakage equation

A new volume and number consistent finite volume scheme for the numerical solution of a collisional nonlinear breakage problem is introduced. The number consistency is achieved by introducing a single weight function in the flux formulation of finite volume scheme, whereas existing schemes for a linear fragmentation equation [Kumar et al. SIAM J. Numer. Anal. 53 (4), 1672-1689] and standard collisional nonlinear breakage equation [Das et al. SIAM J. Sci. Comp. 42 (6), B1570-B1598] require two weights for preserving both volume and number of particles. The higher efficiency and robustness of the proposed scheme allow it to be easily coupled with computational fluid dynamics (CFD) softwares such as COMSOL, Ansys and gPROMS, which is currently one of the predominant topics of discussion in particle technology. Consistency and stability via Lipschitz criterion are studied in detail to demonstrate second order convergence rate for the proposed scheme irrespective of both breakage kernel and nature of grids. Several benchmark problems are solved and validated against its analytical solution to analyze the accuracy of the new scheme.

Modelling of the water on deck problem for an offshore supply vessel running in regular astern waves

The motion of trapped water on the weather deck of an offshore supply vessel running in regular astern seas is modelled using shallow water equations. The inflow and outflow through the weather deck openings and over the bulwarks are modelled using boundary conditions in the flow domain. The hyperbolic flow equations are numerically solved using a higher-order finite volume method. The prescribed motion data for the weather deck is obtained directly from systematic model experiments. Simulated results of the deck water motion are compared against monitored data from the experiments, recorded using a system of pressure sensors and cameras placed on the weather deck. Additionally, the moment caused by the motion of the trapped water and its variation with the ship speed are investigated.

A new type of increasingly higher order finite difference and finite volume MR-WENO schemes with adaptive linear weights for hyperbolic conservation laws

In this paper, the new fifth-order, seventh-order, and ninth-order finite difference and finite volume multi-resolution weighted essentially non-oscillatory (MR-WENO) schemes with adaptive linear weights are presented for hyperbolic conservation laws on structured meshes. They are termed as high-order finite difference and finite volume ALW-WENO schemes. These ALW-WENO schemes only apply one small stencil and one large stencil in reconstruction processes, which could achieve the desired accuracy in the region of smoothness and non-oscillatory properties in the region of containing strong shocks. The linear weights that sum to one can be automatically adjusted to any positive numbers. This is the first time that arbitrary high-order finite difference and finite volume WENO schemes are designed by using only two unequal-sized central spatial stencils. The structure of these novel WENO schemes is simple, so it is easier for obtaining high-order accuracy and solving multi-dimensional problems in large scale engineering applications. Compared to traditional MR-WENO schemes with same order, the computational efficiency can be further improved. Some benchmark tests indicate that these new ALW-WENO schemes have good robustness and performance.

On high order positivity-preserving well-balanced finite volume methods for the Euler equations with gravitation

This paper studies three high-order structure-preserving finite volume weighted essentially non-oscillatory (WENO) methods, which are not only well balanced (WB) for a general known hydrostatic equilibrium state but also preserve the positivity of density and pressure, for the compressible Euler equations under gravitational fields. These methods are built on a simple local scaling positivity-preserving (PP) limiter and a modified WENO-ZQ reconstruction exactly preserving the cell average value and scaling invariance. The WB properties of these three methods are achieved based on suitable numerical fluxes and approximation to the gravitational source terms. Based on some convex decomposition techniques as well as several critical properties of the admissible states and numerical flux, we carry out rigorous positivity-preserving analyses for these three WB schemes. We rigorously prove that the three WB methods, coupled with the PP limiter and a strong-stability-preserving time discretization, are always PP under suitable Courant-Friedrichs-Lewy conditions. Extensive numerical examples are provided to confirm WB and PP properties of three methods.

Imposing slip conditions on curved boundaries for 3D incompressible flows with a very high-order accurate finite volume scheme on polygonal meshes

The conventional no-slip boundary condition does not always hold in several fluid flow applications and must be replaced with appropriate slip conditions according to the wall and fluid properties. However, not only slip boundary conditions are still a subject of discussion among fluid dynamicists, but also their numerical treatment is far from being well-developed, particularly in the context of very high-order accurate methods. The complexity of these conditions significantly increases when the boundary is not aligned with the chosen coordinate system and, even more challenging, when the fluid slips along a curved boundary. The present work proposes a simple and efficient numerical treatment of general slip boundary conditions on arbitrary curved boundaries for three-dimensional fluid flow problems governed by the incompressible Navier–Stokes equations. In that regard, two critical challenges arise: (i) achieving very high-order of convergence with arbitrary curved boundaries for the classical no-slip boundary conditions and (ii) extending the developed numerical techniques to impose general slip boundary conditions. The conventional treatment of curved boundaries relies on generating curved meshes to eliminate the geometrical mismatch between the physical and computational boundaries and achieve high-order of convergence. However, such an approach requires sophisticated meshing algorithms, cumbersome quadrature rules on curved elements, and complex non-linear transformations. In contrast, the reconstruction for off-site data (ROD) method handles arbitrary curved boundaries approximated with linear piecewise elements while employing polynomial reconstructions with specific linear constraints to fulfil the prescribed boundary conditions. For that purpose, the general slip boundary conditions are reformulated on a local orthonormal basis to allow a straightforward application of the ROD method with scalar boundary conditions. The Navier–Stokes equations are discretised with a staggered finite volume method, and the numerical fluxes are computed solely on the polygonal mesh elements. Several benchmark test cases of fluid flow problems in non-trivial three-dimensional curved domains are addressed and confirm that the proposed method effectively achieves up to the eighth-order of convergence.