Within OpenFOAM, we develop a pressure-based solver for the Euler equations written in conservative form using density, momentum, and total energy as variables. Under simplifying assumptions, these equations are used to describe non-hydrostatic atmospheric flow. For the stabilization of the Euler equations and to capture sub-grid processes, we consider two Large Eddy Simulation models: the classical Smagorinsky model and the one equation eddy-viscosity model. To achieve high computational efficiency, our solver uses a splitting scheme that decouples the computation of each variable. The numerical results obtained with our solver are validated against numerical data available in the literature for two classical benchmarks: the rising thermal bubble and the density current. Through qualitative and quantitative comparisons, we show that our approach is accurate. This paper is meant to lay the foundation for a new open-source package specifically created for the quick assessment of new computational approaches for the simulation of atmospheric flows at the mesoscale level.

We present a stabilized POD–Galerkin reduced order method (ROM) for a Leray model. For the implementation of the model, we combine a two-step algorithm called Evolve-Filter (EF) with a computationally efficient finite volume method. In both steps of the EF algorithm, velocity and pressure fields are approximated using different POD basis and coefficients. To achieve pressure stabilization, we consider and compare two strategies: the pressure Poisson equation and the supremizer enrichment of the velocity space. We show that the evolve and filtered velocity spaces have to be enriched with the supremizer solutions related to both evolve and filter pressure fields in order to obtain stable and accurate solutions with the supremizer enrichment method. We test our ROM approach on a 2D unsteady flow past a cylinder at Reynolds number $0\le Re\le 100$. We find that both stabilization strategies produce comparable errors in the reconstruction of the lift and drag coefficients, with the pressure Poisson equation method being more computationally efficient.

The goal of this paper is to test solids4Foam, the fluid-structure interaction (FSI) toolbox developed for foam-extend (a branch of OpenFOAM), and assess its flexibility in handling more complex flows. For this purpose, we consider the interaction of an incompressible fluid described by a Leray model with a hyperelastic structure modeled as a Saint Venant-Kirchho material. We focus on a strongly coupled, partitioned fluid-structure interaction (FSI) solver in a finite volume environment, combined with an arbitrary Lagrangian-Eulerian approach to deal with the motion of the fluid domain. For the implementation of the Leray model, which features a nonlinear differential low-pass filter, we adopt a three-step algorithm called Evolve-Filter-Relax. We validate our approach against numerical data available in the literature for the 3D cross flow past a cantilever beam at Reynolds number 100 and 400.

A data-driven reduced order model (ROM) based on a proper orthogonal decomposition-radial basis function (POD-RBF) approach is adopted in this paper for the analysis of blood flow dynamics in a patient-specific case of atrial fibrillation (AF). The full order model (FOM) is represented by incompressible Navier–Stokes equations, discretized with a finite volume (FV) approach. Both the Newtonian and the Casson’s constitutive laws are employed. The aim is to build a computational tool able to efficiently and accurately reconstruct the patterns of relevant hemodynamics indices related to the stasis of the blood in a physical parametrization framework including the cardiac output in the Newtonian case and also the plasma viscosity and the hematocrit in the non-Newtonian one. Many FOM-ROM comparisons are shown to analyze the performance of our approach as regards errors and computational speed-up.

Reduced order modeling has gained considerable attention in recent decades owing to the advantages offered in reduced computational times and multiple solutions for parametric problems. The focus of this manuscript is the application of model order reduction techniques in various engineering and scientific applications including but not limited to mechanical, naval and aeronautical engineering. The focus here is kept limited to computational fluid mechanics and related applications. The advances in the reduced order modeling with proper orthogonal decomposition and reduced basis method are presented as well as a brief discussion of dynamic mode decomposition and also some present advances in the parameter space reduction. Here, an overview of the challenges faced and possible solutions are presented with examples from various problems.