Second-order Finite Volume Research Articles

Wind tunnel testing and modelling of a rod-airfoil setup: validation and 3D acoustic imaging

This work revisits the rod-airfoil benchmark experiment from Jacob et al. (2004) with the aim of performing a detailed numerical validation in the nearfield and far field as well as to perform 3D acoustic imaging. The setup consists of a rod near the leading edge of a NACA-0012 airfoil and is placed in a semi-open wind tunnel in anechoic conditions. In addition to a far-field microphone array, the setup is equipped with two planar microphone arrays as well as wall-mounted microphones on the closed wind tunnel section and on the airfoil. The setup is modelled using a second-order finite volume method, and particular attention is devoted to including installation effects and sensor response properties in order to guarantee simulation accuracy. Validation examples are provided, including single microphone spectra as well as the acoustic field reconstructed on the test object by means of 3D acoustic imaging using the planar microphone arrays. Finally, the experimentally-validated numerical model is used as a means to gain further insight into the tested object. For instance, the model allows to separate acoustical and turbulent effects present in the sensors located within the flow section.

A well-balanced finite volume solver for the 2D shallow water magnetohydrodynamic equations with topography

In this paper, a second-order finite volume Non-Homogeneous Riemann Solver is used to obtain an approximate solution for the two-dimensional shallow water magnetohydrodynamic (SWMHD) equations considering non-flat bottom topography. We investigate the stability of a perturbed steady state, as well as the stability of energy in these equations after a perturbation of a steady state using a dispersive analysis. To address the elliptic constraint ∇⋅hB=0, the GLM (Generalized Lagrange Multiplier) method designed specifically for finite volume schemes, is used. The proposed solver is implemented on unstructured meshes and verifies the exact conservation property. Several numerical results are presented to validate the high accuracy of our schemes, the well-balanced, and the ability to resolve smooth and discontinuous solutions. The developed finite volume Non-Homogeneous Riemann Solver and the GLM method offer a reliable approach for solving the SWMHD equations, preserving numerical and physical equilibrium, and ensuring stability in the presence of perturbations.

An arbitrary Lagrangian–Eulerian positivity-preserving finite volume scheme for radiation hydrodynamics equations in the equilibrium-diffusion limit

In this paper, an arbitrary Lagrangian–Eulerian (ALE) positivity-preserving finite volume scheme is constructed for radiation hydrodynamics equations (RHE) in the equilibrium-diffusion limit. Firstly, the integral form equations of RHE in the ALE framework are presented. Then, the operator-splitting method is applied to divide the equations into hyperbolic part and parabolic part. In addition, a second-order positivity-preserving finite volume scheme is constructed for the hyperbolic part based on MUSCL reconstruction. The vertex velocity is obtained by the predictor–corrector strategy. The Winslow method is applied to improve the quality of the Lagrangian mesh. Furthermore, a nonlinear positivity-preserving finite volume scheme suitable for distorted mesh is proposed for the parabolic part. Finally, some numerical examples are given to show the accuracy and reliability of the numerical scheme.

A well-balanced second-order finite volume approximation for a coupled system of granular flow

A well-balanced second-order finite volume scheme is proposed and analyzed for a 2×2 system of non-linear partial differential equations which describes the dynamics of growing sandpiles created by a vertical source on a flat, bounded rectangular table in multiple dimensions. To derive a second-order scheme, we combine a MUSCL type spatial reconstruction with strong stability preserving Runge-Kutta time stepping method. The resulting scheme is ensured to be well-balanced through a modified limiting approach that allows the scheme to revert to a well-balanced first-order scheme if it fails to capture or preserve the discrete steady state solution during its time evolution. Further, the scheme maintains the second-order accuracy away from the steady state. The well-balance property of the scheme is proven analytically in one dimension and demonstrated numerically in two dimensions. Additionally, numerical experiments reveal that the second-order scheme reduces finite time oscillations and provides better resolutions of the physical properties of the state variables, as compared to the existing first-order schemes of the literature.

Energy stable and structure-preserving schemes for the stochastic Galerkin shallow water equations

The shallow water flow model is widely used to describe water flows in rivers, lakes, and coastal areas. Accounting for uncertainty in the corresponding transport-dominated nonlinear PDE models presents theoretical and numerical challenges that motivate the central advances of this paper. Starting with a spatially one-dimensional hyperbolicity-preserving, positivity-preserving stochastic Galerkin formulation of the parametric/uncertain shallow water equations, we derive an entropy-entropy flux pair for the system. We exploit this entropy-entropy flux pair to construct structure-preserving second-order energy conservative, and first- and second-order energy stable finite volume schemes for the stochastic Galerkin shallow water system. The performance of the methods is illustrated on several numerical experiments.

A solver for subsonic flow around airfoils based on physics-informed neural networks and mesh transformation

Physics-informed neural networks (PINNs) have recently become a new popular method for solving forward and inverse problems governed by partial differential equations. However, in the flow around airfoils, the fluid is greatly accelerated near the leading edge, resulting in a local sharper transition, which is difficult to capture by PINNs. Therefore, PINNs are still rarely used to solve the flow around airfoils. In this study, we combine physical-informed neural networks with mesh transformation, using a neural network to learn the flow in the uniform computational space instead of physical space. Mesh transformation avoids the network from capturing the local sharper transition and learning flow with internal boundary (wall boundary). We successfully solve inviscid flow and provide an open-source subsonic flow solver for arbitrary airfoils. Our results show that the solver exhibits higher-order attributes, achieving nearly an order of magnitude error reduction over second-order finite volume method (FVM) on very sparse meshes. Limited by the learning ability and optimization difficulties of the neural network, the accuracy of this solver will not improve significantly with mesh refinement. Nevertheless, it achieves comparable accuracy and efficiency to second-order FVM on fine meshes. Finally, we highlight the significant advantage of the solver in solving parametric problems, as it can efficiently obtain solutions in the continuous parameter space about the angle of attack.

Erratum to : “A family of second-order dissipative finite volume schemes for hyperbolic systems of conservation laws”

Erratum to : “A family of second-order dissipative finite volume schemes for hyperbolic systems of conservation laws”

A robust unstructured finite volume scheme for urban flash flood simulation: Duhok city as a case study

ABSTRACT A robust well-balanced finite volume scheme is proposed for simulating the flash floods due to rainfall-runoff process. The solver’s implementation on unstructured triangular meshes is based on a Godunov-type second-order upwind finite volume formulation, and Roe’s average scheme is used to determine the fluxes. The bed slope, friction and rainfall source terms are discretized using an upwind approach. The suggested algorithm uses the slope limiter and Hancock’s predictor-corrector method in order to secure numerical stability. The proposed model overwhelms the drawbacks of many commercial codes that use the finite difference, finite elements and finite volume in their scheme, especially in dealing with rapid change in topography. The prepared scheme solves the two-dimensional shallow water equations over a real topography of Hishkaro basin located in Duhok city in Kurdistan of Iraq where the variation of bed elevation reaches 717 m. The scheme is validated against real observed data in Duhok city provided by the Sewage directory. The data of five observation points over the study area were used to validate the code. The results showed that the max water level reaches 3.67, 2.12, 1.20, 0.67 and 0.3 m at observed points 3, 2, 4, 1 and 5, respectively.

A conservative sharp interface method for two-dimensional incompressible two-phase flows with phase change

A conservative sharp interface method is proposed in this work to simulate two-dimensional/axisymmetric incompressible two-phase flows with phase change. In this method, we use the cut cell method to generate unstructured meshes near the interface, of which the cell edges overlap with the interface at each time step. On such mesh, the mass and heat transfer during phase change and all the jump conditions can be incorporated into the calculation of fluxes at the cell edges, to ensure that they are strictly satisfied at the interface in a sharp manner. The governing equations, including the incompressible Navier–Stokes equations, heat equation, and vapor mass fraction equation, are discretized by a second-order finite volume method in the arbitrary Lagrangian–Eulerian framework. To well couple the mass, heat, momentum, and interface evolution, the solution procedure is carefully designed and performed with several techniques. In such a way, the sharp discontinuity of the velocity, stress, temperature gradient, and vapor fraction, caused by the mass/heat transfer during phase change, can be simulated accurately and robustly. The performance of this method is systematically examined by cases of phase change at or below the saturated temperature, including vapor bubble in superheated liquid, film boiling, droplet evaporation at different relative humidity conditions, droplet evaporation under gravity, and droplet evaporation under forced convection. The applicability of the present method for incompressible two-phase flows with phase change is well demonstrated by comparing the numerical results with the benchmark, theoretical or experimental ones.

A Finite Volume Method for a 2D Dam-Break Simulation on a Wet Bed Using a Modified HLLC Scheme

This study proposes a numerical model for depth-averaged Reynolds equations (shallow-water equations) to investigate a dam-break problem, based upon a two-dimensional (2D) second-order upwind cell-centre finite volume method. The transportation terms were modelled using a modified approximate HLLC Riemann solver with the first-order accuracy. The proposed 2D model was assessed and validated through experimental data and analytical solutions for several dam-break cases on a wet and dry bed. The results showed that the error values of the model are lower than those of existing numerical methods at different points. Our findings also revealed that the dimensionless error parameters decrease as the wave propagates downstream. In general, the new model can model the dam-break problem and captures the shock wave superbly.

Theoretical Study of Non-Isothermal Gradient Elution Liquid Chromatography.

A two-component model of gradient elution chromatography is investigated to theoretically study the effects of simultaneous variations in temperature and solvent strength on the retention behaviors of elution profiles in thermally insulated liquid chromatographic columns. The gradient elution technique is based on the gradual increase or decrease in eluent strength during the chromatographic operation by varying the composition of the mobile phase. The enthalpy of adsorption is primarily responsible for internal temperature variations inside the column, as heat adsorbs during the adsorption process and releases in the desorption phase. Both types of variations change the propagation speeds of moving pulses inside the column which can lead to better separation of the components and a reduction in the recycling time for the next injection. The equilibrium dispersive model (EDM) coupled with the energy balance equation for temperature and transport equation for the volume fraction of the solvent is utilized to simulate this complex process. The resulting nonlinear model equations are approximated by applying a semi-discrete second-order finite volume scheme. The numerical solutions are used to study the impact of a gradient starting and ending times, volume-fraction of the solvent, solvent strength parameter, the slope of gradient, enthalpy of adsorption, injection temperature, and the ratio of specific heats on the concentration profiles.

Simulation of laminar to transitional wakes past cylinders with a discontinuous Galerkin inviscid shallow water model

Laminar to transitional wakes occur in slow, quasi-steady flows past cylinders at low cylinder Reynolds numbers (R ed ≤ 250). Inviscid numerical solvers of the depth-averaged shallow water equations (SWE) introduce numerical dissipation that, depending on Red , may imitate the mechanisms of viscous turbulent models. However, the numerical dissipation rate in a second-order finite volume (FV2) SWE solver is so large at a practical resolution that this can instead hide these mechanisms. The extra numerical complexity of the second-order discontinuous Galerkin (DG2) SWE solver results in a lower dissipation rate, making it a potential alternative to the FV2 solver to reproduce cylinder wakes. This paper compares the DG2 and FV2 solvers, initially for wake formation behind one cylinder. The findings confirm that DG2 can reproduce the expected wake formations, which FV2 fails to capture, even at a 10-fold finer resolution. It is further demonstrated that DG2 is capable of reproducing key features of the flow fields observed in a laboratory random cylinder array.

Analysis of Liquid Chromatography Considering a Linear Single-Component Heterogeneous-Type Reactive General Rate Model.

This paper elaborates on the significance of liquid chromatography for a single-component reactive linear general rate model. The model equations consist of a set of two coupled partial differential equations, which include diffusion, interfacial mass transfer, axial dispersion, external and intraparticle pore diffusivity, and heterogeneous chemical reaction of the first order with two sets of boundary conditions. The model equations are solved by the Laplace transformation. The actual time domain solution is obtained by numerical Laplace inversion, as analytical inversion cannot be obtained. The graphical sketch of different physical parameters is presented to analyze the dynamics of the elution profiles. The result indicates that the chromatographic reactor works more efficiently on increasing the value of the heterogeneous-type first-order reaction. To check the analytical results, a second-order high-resolution finite volume scheme is used. Both results are in good agreement and indicate the correctness of the numerical scheme. The current work is also compared with the previously available numerical schemes, which shows that the proposed numerical scheme is better for elaborating the chromatographic reactor performance. A comparison table is also presented to compute error analysis and computational run time for analyzing the efficiency of the reactor. A graphical sketch of the numerical temporal moment analysis is also presented, which gives significant information about the performance and the shape of the concentration profiles.

A Second-Order Finite Volume Method for Field-Scale Reservoir Simulation

Subsurface reservoirs are large complex systems. Reservoir flow models are defined on complex grids that follow geology with relatively large block sizes to make consistent simulations feasible. Reservoir engineers rely on established reservoir simulation software to model fluid flow. Nevertheless, fluid front position inaccuracies and front smearing on large grids may cause significant errors and make it hard to predict hydrocarbon production efficiency. We investigate higher-order methods that reduce these undesired effects without refining the grid, thus making reservoir simulation more accurate and robust. For this paper, we implemented a second-order finite volume method with linear programming (LP) reconstruction in the open-source industry-grade reservoir simulator OPM Flow (part of the open porous media initiative, OPM). We benchmark it against the first-order method on full-scale cases with standard coarse and refined grids. We prepared open refined-grid models of a synthetic reservoir with an unstructured grid and refined Norne field example. Our results confirm that the LP method predicts front positions as accurately as the first-order method on the refined grid for problems dominated by transport. These include the water alternating gas scenario on the synthetic reservoir and piston-type injection on the Norne field. Moreover, we study the gains from the LP method for CO_{2} injection problems on the Norne field with full multi-phase complexity beyond transport. We observe the relevant difference between the first- and the second-order methods in these cases. However, in some configurations, the reservoir complexity overshadows the gains from the second-order methods.

Transient solidification and melting numerical simulation of lauric acid PCM filled stepped solar still basin used in water desalination process

It is crucial to produce potable water from salt water in arid places using a stepped solar still desalination process. This work presents the temperature distribution and phase transformation of lauric acid PCM in inclined stepped solar still desalination system. A second-order finite volume method is employed for discretizing the two-dimensional geometry, while energy, continuity, and momentum equations are solved and then coupled using the second-order PRESTO algorithm. The solidification and melting CFD models have been used to capture the melting stages of lauric acid PCM. Four steps with an inclination of 30° to ground and two different cases with saline water levels of 10 mm and 40 mm have been reported. The results reveal that the temperature and mass fraction of lauric acid PCM strongly depends on the saline water level at each step, the number of steps, the angle of inclination, and the location of the solar still step. The maximum percentage increase in PCM melting fraction is 50%, and a PCM temperature rise of 4 °C is observed in the case of SWL-max compared to SWL-min. The total enthalpy of lauric acid PCM is increased by 30% for the SWL-max water level compared to the SWL-min water level.

Finite Volume Method for Transient Pipe Flow with an Air Cushion Surge Chamber Considering Unsteady Friction and Experimental Validation

In various water transmission systems such as long-distance water transfer projects and hydropower stations, accurate simulation of water hammer is extremely important for safe and stable operation and the realization of intelligent operations. Previous water hammer calculations usually consider only steady-state friction, underestimating the decay of transient pressure. A second-order Finite Volume Method (FVM) considering the effect of unsteady friction factor is developed to simulate the water hammer and the dynamic behavior of air cushion surge chamber in a water pipeline system, while an experimental pipe system is conducted to validate the proposed numerical model. Two unsteady friction models, Brunone and TVB models, were incorporated into the water hammer equations, which are solved by the MUSCL–Hancock method. One virtual boundary method was proposed to realize the FVM simulation of Air Cushion Surge Chamber. Comparisons with water hammer experimental results show that, while the steady friction model only accurately predicts the first pressure peak, it seriously underestimates pressure attenuation in later stages. Incorporating an unsteady friction factor can better predict the entire pressure attenuation process; in particular, the TVB unsteady friction model more accurately reproduces the pressure peaks and the whole pressure oscillation periods. For water pipeline systems with an air cushion surge chamber, energy attenuation of the elastic pipe water hammer is primarily due to pipe friction and the air cushion. The experimental results with the air cushion surge chamber demonstrate that the proposed FVM model with the TVB unsteady friction model and the air chamber polytropic exponent near 1.0 can well reproduce the experimental pressure oscillations.

A three-dimensional high-order flux reconstruction lattice boltzmann flux solver for incompressible laminar and turbulent flows

A three-dimensional (3D) high-order flux reconstruction lattice Boltzmann flux solver (FR-LBFS) is developed in this paper for accurate and efficient simulations of incompressible laminar and turbulent flows. The original lattice Boltzmann flux solver (LBFS) is a second-order finite volume method (FVM) which combines the advantages of conventional Navier-Stokes (N-S) solvers and lattice Boltzmann equation (LBE) solvers. But it is not convenient to construct high-order FVM, and the compactness is always a bottleneck. The present work separates LBFS from FVM and combines LBFS with high-order FR scheme. Thus, the FR-LBFS inherits the superiority of LBFS and FR and shares the following attractive features: (i) a fully-explicit arbitrary order compact method, (ii) weakly compressible (WC) model for unsteady incompressible flows and (iii) evaluating inviscid and viscous fluxes simultaneously and no particular technique requirement, such as local discontinuous Galerkin method (LDG), for viscous term discretization. To validate the current method, various 3D incompressible laminar isothermal and thermal numerical experiments are presented first. Good agreements are achieved between the current results and the previous experimental and numerical data whilst the present high-order solver adopt coarse meshes. Furthermore, the ability of the proposed method to perform accurate and stable computations of incompressible decaying homogeneous isotropic turbulent flows and turbulent channel flows with different mesh resolutions and accuracy orders are investigated. The results show that the present 3D FR-LBFS is a promising tool for under-resolved or implicit large eddy simulation (ILES) of incompressible turbulent flows.

Second-order accurate implicit finite volume method for two-dimensional modeling of PFAS transport in unsaturated porous media with variable surface tension

Per- and polyfluoroalkyl substances (PFASs) have become emerging contaminants of critical concern. Comprehensive understanding of the transport and fate of PFAS in the vadose-zone, a type of water-unsaturated porous media, is key to determination of the risks of the PFAS contamination in the subsurface and to the development of the effective remediation strategies. Accurate modeling of the PFAS transport in the unsaturated porous media is still a challenge due to the variable surface tension induced by the adsorption of PFAS to the air-water interfaces. In an effort to address this challenge, we propose a multidimensional modeling framework for the transient PFAS transport in the unsaturated porous media based on the second order accuracy finite volume method. In the modeling, the adsorption of PFAS to the solid surfaces and to the air-water interfaces is described by the two-domain sorption kinetics model, i.e., both the instantaneous and the rate limited PFAS adsorptions are taken into account. The diffusive and convective terms in the governing equations for the PFAS transport and the water flow are discretized by the central difference and the quadratic upstream interpolation for convective kinetics schemes, respectively. We investigate the effects of the convergent criteria, coupling method, and variation of the surface tension on the average and the local PFAS concentration and water content in the computational domain. We find that the convergent criteria should be chosen carefully so as to get the accurate results. The differences between the different coupling methods are affected by the boundary conditions. The variation in the surface tension due to the variation of the PFAS concentration cannot be neglected. The dimensionless parameters relevant to the properties of porous media as well as the relation between the surface tension and the PFAS concentration play an important role in transport of PFAS in the unsaturated porous media. The effects of the Péclet number, Damköhler values, and fraction of instantaneous sorption are not significant in the ranges studied in this work. These findings provide a better understanding of transport of PFAS in the vadose zone.

Numerical Dissipation Rate Analysis of Finite-Volume and Continuous-Galerkin Methods for LES of Combustor Flow-Field

A growing body of literature indicates that element-based high-order methods can exhibit considerable accuracy/cost benefit over conventional second-order finite-volume (FV) methods for large-eddy simulations (LES). This may even hold true for complex configurations relevant to industry that involve under-resolving unstructured grids. However, it is not often clear whether the accuracy/cost benefit stems from the low-dissipative nature of the high-order numerical schemes or from using a different LES approach (implicit/explicit), or a combination of the two. The present paper employs a numerical dissipation rate analysis technique due to Schranner et al. (Comput Fluids 114:84–97, 2015) to better understand the reasons for the high-order benefit seen previously on a complex LES test case related to gas-turbine combustors. It is established that a high(fifth)-order LES run provides better accuracy than its second-order FV counterpart at the same computational cost primarily because of lower numerical dissipation and the LES model dissipation has a secondary role to play. The numerical dissipation is found to contribute 60–90% of the total (numerical and LES model) dissipation.

Global entropy stability for a class of unlimited second-order schemes for 2D hyperbolic systems of conservation laws on unstructured meshes

We introduce a class of second-order accurate finite volume schemes defined on two dimensional unstructured meshes. They satisfy a global discrete entropy stability for an arbitrary strictly convex entropy provided the second-order discretization verifies a dissipative inequality. These schemes do not require usual limitation techniques. Numerical experiments are performed to confirm the accuracy and the shock capturing efficiency of the schemes.