Search

We present a filter stabilization technique for the mildly compressible Euler equations that relies on a linear or nonlinear indicator function to identify the regions of the domain where artificial viscosity is needed and determine its amount. For the realization of this technique, we adopt a three step algorithm called Evolve-Filter-Relax (EFR), which at every time step evolves the solution (i.e., solves the Euler equations on a coarse mesh), then filters the computed solution, and finally performs a relaxation step to combine the filtered and non-filtered solutions. We show that the EFR algorithm is equivalent to an eddy-viscosity model in Large Eddy Simulation. Three indicator functions are considered: a constant function (leading to a linear filter), a function proportional to the norm of the velocity gradient (recovering a Smagorinsky-like model), and a function based on approximate deconvolution operators. Through well-known benchmarks for atmospheric flow, we show that the deconvolution-based filter yields stable solutions that are much less dissipative than the linear filter and the Samgorinsky-like model and we highlight the efficiency of the EFR algorithm.

The investigation of fluid-solid systems is very important in a lot of industrial processes. From a computational point of view, the simulation of such systems is very expensive, especially when a huge number of parametric configurations needs to be studied. In this context, we develop a non-intrusive data-driven reduced order model (ROM) built using the proper orthogonal decomposition with interpolation (PODI) method for Computational Fluid Dynamics (CFD) -- Discrete Element Method (DEM) simulations. The main novelties of the proposed approach rely in (i) the combination of ROM and FV methods, (ii) a numerical sensitivity analysis of the ROM accuracy with respect to the number of POD modes and to the cardinality of the training set and (iii) a parametric study with respect to the Stokes number. We test our ROM on the fluidized bed benchmark problem. The accuracy of the ROM is assessed against results obtained with the FOM both for Eulerian (the fluid volume fraction) and Lagrangian (position and velocity of the particles) quantities. We also discuss the efficiency of our ROM approach.

Tumble dryers offer a fast and convenient way of drying textiles independent of weather conditions and therefore are frequently used in ordinary households. However, artificial drying of textiles consumes considerable amounts of energy, approximately 8.2 percent of the residential electricity consumption is for drying of textiles in northern European countries (Cranston et al., 2019). Several authors have investigated the aspects of the clothes drying cycle with experimental and numerical methods to understand and improve the process. The first turning point study on understanding the physics of evaporation for tumble dryers was presented by Lambert et al. (1991) in the early 90s. With the aid of Chilton_Colburn analogy, they introduced the concept of area-mass transfer coefficient to address evaporation rate. Afterwards, several experimental or numerical studies were published based on this concept, and furthermore, the model was then developed into 0-dimensional (Deans, 2001) and 1-dimensional (Wei et al., 2017) to gain more accuracy. The evaporation rate is considered to be the main system parameter for dryers with which other performance parameters including drying time, effectiveness, moisture content and efficiency can be estimated. More recent literature focused on utilizing dimensional analysis or image processing techniques to correlate drying indices with system parameters. However, the validity of these regressed models is machine-specific, and hence, cannot be generalized yet. All the previous models for estimating the evaporation rate in tumble dryers are discussed. The review of the related literature showed that all of the previous models for the prediction of the evaporation rate in the clothes dryers have some limitations in terms of accuracy and applicability., Comment: 25 Pages

The simulation of atmospheric flows by means of traditional discretization methods remains computationally intensive, hindering the achievement of high forecasting accuracy in short time frames. In this paper, we apply three reduced order models that have successfully reduced the computational time for different applications in computational fluid dynamics while preserving accuracy: Dynamic Mode Decomposition (DMD), Hankel Dynamic Mode Decomposition (HDMD), and Proper Orthogonal Decomposition with Interpolation (PODI). The three methods are compared in terms of computational time and accuracy in the simulation of two well-known benchmarks for mesoscale flow. The accuracy of the DMD and HDMD solutions deteriorates rather quickly as the forecast time window expands, although these methods are designed to predict the dynamics of a system. The reason is likely the strong nonlinearity in the benchmark flows. The PODI solution is accurate for the entire duration of the time interval of interest thanks to the use of interpolation with radial basis functions. This holds true also when the model features a physical parameter expected to vary in a given range, as is typically the case in weather prediction.

We introduce GEA (Geophysical and Environmental Applications), a new open-source atmosphere and ocean modeling framework within the finite volume C++ library OpenFOAM. Here, we present the development of a non-hydrostatic atmospheric model consisting of a pressure-based solver for the Euler equations written in conservative form using density, momentum, and total energy as variables. We validate the solver for two idealized test cases involving buoyancy driven flows: smooth and non-smooth rising thermal bubble. Through qualitative and quantitative comparisons against numerical data available in the literature, we show that our approach is accurate., Comment: 8 pages, 2 figures, 2 tables. arXiv admin note: substantial text overlap with arXiv:2302.04836

In this work a machine learning-based Reduced Order Model (ROM) is developed to investigate in a rapid and reliable way the hemodynamic patterns in a patient-specific configuration of Coronary Artery Bypass Graft (CABG). The computational domain is composed by the left branches of coronary arteries when a stenosis of the Left Main Coronary Artery (LMCA) occurs. A reduced basis space is extracted from a collection of Finite Volume (FV) solutions of the incompressible Navier-Stokes equations by using the Proper Orthogonal Decomposition (POD) algorithm. Artificial Neural Networks (ANNs) are employed to compute the modal coefficients. Stenosis is introduced by morphing the volume meshes with a Free Form Deformation (FFD) by means of a Non-Uniform Rational Basis Spline (NURBS) volumetric parameterization.

Partial differential equations can be used to model many problems in several fields of application including, e.g., fluid mechanics, heat and mass transfer, and electromagnetism. Accurate discretization methods (e.g., finite element or finite volume methods, the so-called full order models) are widely used to numerically solve these problems. However, when many physical and/or geometrical parameters are involved, the computational cost required by full order models becomes prohibitively expensive and this is not acceptable for real-time computations that are becoming more and more popular for rapid prototyping. Therefore, there is the need to introduce reduced order methods (also referred to as reduced basis methods) able to provide, as the input parameters change, fast and reliable solutions at a reduced computational cost.

In this work, the Immersed Boundary Method (IBM) with feedback forcing introduced by Goldstein et al. (1993) and often referred in the literature as the Virtual Boundary Method (VBM), is addressed. The VBM has been extensively applied both within a Spectral and a Finite Difference (FD) framework. Here, we propose to combine the VBM with a computationally efficient Finite Volume (FV) method. We will show that for similar computational configurations, FV and FD methods provide significantly different results. Furthermore, we propose to modify the standard feedback forcing scheme, based on a Proportional-Integral (PI) controller, with the introduction of a derivative action, in order to obtain a Proportial-Integral-Derivative (PID) controller. The stability analysis for the Backward Differentiation Formula of order 1 (BDF1) time scheme is modified accordingly, and extended to the Backward Differentiation Formula of order 2 (BDF2) time scheme. We will show that, for the BDF2 time scheme, the derivative action allows to improve the stability characteristics of the system. Our approach is validated against numerical data available in the literature for a stationary/rigidly moving 2D circular cylinder in several configurations. Finally, a Fluid-Structure Interaction (FSI) benchmark, related to the frequency response of a cantilever beam coupled with a fluid, is presented: we numerically demonstrate that the introduction of the derivative action plays an important role in order to properly detect the fluid-structure interaction coupling., Comment: 23 pages, 18 figures, 9 tables