Search

This paper describes the numerical implementation in a high-performance computing environment of an open-source library for model order reduction in fluid dynamics. This library, called pyLOM, contains the algorithms of proper orthogonal decomposition (POD), dynamic mode decomposition (DMD) and spectral proper orthogonal decomposition (SPOD), as well as, efficient SVD and matrix-matrix multiplication, all of them tailored for supercomputers. The library is profiled in detail under the MareNostrum IV supercomputer. The bottleneck is found to be in the QR factorization, which has been solved by an efficient binary tree communications pattern. Strong and weak scalability benchmarks reveal that the serial part (i.e., the part of the code that cannot be parallelized) of these algorithms is under 10% for the strong scaling and under 0.7% for the weak scaling. Using pyLOM, a POD of a dataset containing 1.14 x 108 gridpoints and 1808 snapshots that takes 6.3Tb of memory can be computed in 81.08 seconds using 10368 CPUs. Additioally, the algorithms are validated using the datasets of a flow around a circular cylinder at ReD = 100 and ReD = 1 x 104, as well as the flow in the Stanford diffuser at Reh = 1 x 104.

The control efficacy of classical periodic forcing and deep reinforcement learning (DRL) is assessed for a turbulent separation bubble (TSB) at $Re_\tau=180$ on the upstream region before separation occurs. The TSB can resemble a separation phenomenon naturally arising in wings, and a successful reduction of the TSB can have practical implications in the reduction of the aviation carbon footprint. We find that the classical zero-net-mas-flux (ZNMF) periodic control is able to reduce the TSB by 15.7%. On the other hand, the DRL-based control achieves 25.3% reduction and provides a smoother control strategy while also being ZNMF. To the best of our knowledge, the current test case is the highest Reynolds-number flow that has been successfully controlled using DRL to this date. In future work, these results will be scaled to well-resolved large-eddy simulation grids. Furthermore, we provide details of our open-source CFD-DRL framework suited for the next generation of exascale computing machines., Comment: 19 pages, 14 figures, 3 tables

This paper presents for the first time successful results of active flow control with multiple independently controlled zero-net-mass-flux synthetic jets. The jets are placed on a three-dimensional cylinder along its span with the aim of reducing the drag coefficient. The method is based on a deep-reinforcement-learning framework that couples a computational-fluid-dynamics solver with an agent using the proximal-policy-optimization algorithm. We implement a multi-agent reinforcement-learning framework which offers numerous advantages: it exploits local invariants, makes the control adaptable to different geometries, facilitates transfer learning and cross-application of agents and results in significant training speedup. In this contribution we report significant drag reduction after applying the DRL-based control in three different configurations of the problem., Comment: ETMM14 2023 conference proceeding paper

Deep artificial neural networks (ANNs) used together with deep reinforcement learning (DRL) are receiving growing attention due to their capabilities to control complex problems. This technique has been recently used to solve problems related to flow control. In this work, an ANN trained through a DRL agent is used to perform active flow control. Two-dimensional simulations of the flow around a cylinder are conducted and an active control based on two jets located on the walls of the cylinder is considered. By gathering information from the flow surrounding the cylinder, the ANN agent is able to learn effective control strategies for the jets, leading to a significant drag reduction. In the present work, a Reynolds-number range beyond those previously considered is studied and compared with results obtained using classical flow-control methods. Significantly different nature in the control strategies is identified by the DRL as the Reynolds number Re increases. For Re <= 1000 the classical control strategy based on an opposition control relative to the wake oscillation is obtained. For Re = 2000 the new strategy consists of an energisation of the boundary layers and the separation area, which modulate the flow separation and reduce drag in a fashion similar to that of the drag crisis, through a high frequency actuation. A cross-application of agents is performed for a flow at Re = 2000, obtaining similar results in terms of drag reduction with the agents trained at Re = 1000 and 2000. The fact that two different strategies yield the same performance make us question whether this Reynolds number regime (Re = 2000) belongs to a transition towards a nature-different flow which would only admit a high-frequency actuation strategy to obtain drag reduction. This finding allows the application of ANNs trained at lower Reynolds numbers but comparable in nature, saving computational resources., Comment: 26 pages, 5 figures, submitted to MDPI

In this paper, we describe an upgrade of the Alya code with up-to-date parallel linear solvers capable of achieving reliability, efficiency and scalability in the computation of the pressure field at each time step of the numerical procedure for solving a Large Eddy Simulation formulation of the incompressible Navier-Stokes equations. We developed a software module in Alya's kernel to interface the libraries included in the current version of PSCToolkit, a framework for the iterative solution of sparse linear systems on parallel distributed-memory computers by Krylov methods coupled to Algebraic MultiGrid preconditioners. The Toolkit has undergone various extensions within the EoCoE-II project with the primary goal of facing the exascale challenge. Results on a realistic benchmark for airflow simulations in wind farm applications show that the PSCToolkit solvers significantly outperform the original versions of the Conjugate Gradient method available in the Alya's kernel in terms of scalability and parallel efficiency and represent a very promising software layer to move the Alya code towards exascale.

Communities around the world are increasingly interested in nature-based solutions to mitigation of coastal risks like coastal forests, but it remains unclear how much protective benefits vegetation provides, particularly in the limit of highly energetic flows after tsunami impact. The current study, using a three-dimensional incompressible computational fluid dynamics model with a fluid-structure interaction approach, aims to quantify how energy reflection and dissipation vary with different degrees of rigidity and vegetation density of a coastal forest. We represent tree trunks as cylinders and use the elastic modulus of hardwood trees such as pine or oak to characterize the rigidity of these cylinders. The numerical results show that energy reflection increases with rigidity only for a single cylinder. In the presence of multiple cylinders, the difference in energy reflection created by varying rigidity diminishes as the number of cylinders increases. Instead of rigidity, we find that the blockage area created by the presence of multiple tree trunks dominates energy reflection. As tree trunks are deformed by the hydrodynamic forces, they alter the flow field around them, causing turbulent kinetic energy generation in the wake region. As a consequence, trees dissipate flow energy, highlighting coastal forests reducing the onshore energy flux of tsunamis by means of both reflection and dissipation.

This work assesses Lagrangian droplet evaporation models frequently used in spray combustion simulations, with the purpose of identifying the influence of modeling decisions on the single droplet behavior. Besides more simplistic models, the evaluated strategies include a simple method to incorporate Stefan flow effects in the heat transfer (Bird's correction), a method to consider the interaction of Stefan flow with the heat and mass transfer films (Abramzon-Sirignano model), and a method to incorporate non-equilibrium thermodynamics (Langmuir-Knudsen model). The importance of each phenomena is quantified analytically and numerically under various conditions. Evaporation models ignoring Stefan flow are found to be invalid under the studied conditions. The Langmuir-Knudsen model is also deemed inadequate for high temperature evaporation, while Bird's correction and the Abramzon-Sirignano model are identified as the most relevant for numerical studies of spray combustion systems. Latter is the most elaborate model studied here, as it considers Reynolds number effects beyond the empirical correlation of Ranz and Marshall derived for low-transfer rates. Thus, the Abramzon-Sirignano model is identified as the state of the art alternative in the scope of this study., Comment: 46 pages, 16 figures