Search

For the last 140 years, the mechanisms of transport and dissipation of energy in a turbulent flow have not been completely understood due to the complexity of this phenomenon. The dissipation of energy due to turbulence is significative, and understanding turbulence physics is crucial for fighting the present climate emergency. Previous research has focused on analyzing the so-called coherent structures of the flow (Q events, streaks, and vortices), which are regions of high turbulence transport, high/low streamwise fluctuation, and rotation, respectively. However, the connection between these classically studied structures and the flow development is still uncertain. To fill this gap, here we show a data-driven methodology for objectively identifying high-importance regions in a turbulent flow. A deep-learning model is trained to predict a future state of a turbulent channel flow and the gradient-SHAP explainability algorithm is used to calculate the importance of each grid point for such a prediction. Finally, high-importance regions are computed using the SHAP data, analyzing and comparing their characteristics with those of the other coherent structures. The SHAP analysis provides an objective way to identify the regions of highest importance in the turbulent flow, which exhibit different levels of agreement with the classically studied structures.

The use of data-driven methods in fluid mechanics has surged dramatically in recent years due to their capacity to adapt to the complex and multi-scale nature of turbulent flows, as well as to detect patterns in large-scale simulations or experimental tests. In order to interpret the relationships generated in the models during the training process, numerical attributions need to be assigned to the input features. One important example are the additive-feature-attribution methods. These explainability methods link the input features with the model prediction, providing an interpretation based on a linear formulation of the models. The SHapley Additive exPlanations (SHAP values) are formulated as the only possible interpretation that offers a unique solution for understanding the model. In this manuscript, the additive-feature-attribution methods are presented, showing four common implementations in the literature: kernel SHAP, tree SHAP, gradient SHAP, and deep SHAP. Then, the main applications of the additive-feature-attribution methods are introduced, dividing them into three main groups: turbulence modeling, fluid-mechanics fundamentals, and applied problems in fluid dynamics and heat transfer. This review shows thatexplainability techniques, and in particular additive-feature-attribution methods, are crucial for implementing interpretable and physics-compliant deep-learning models in the fluid-mechanics field.

Direct numerical simulations (DNS) stand out as formidable tools in studying turbulent flows. Despite the fact that the achievable Reynolds number remains lower than those available through experimental methods, DNS offers a distinct advantage: the complete knowledge of the velocity field, facilitating the evaluation of any desired quantity. This capability should extend to compute derivatives. Among the classic functions requiring derivatives is the indicator function, $\Xi(y^+) = y^+\frac{{\rm d}\overline{U}_x^+}{{\rm d}y^+}$. This function encapsulates the wall-normal derivative of the streamwise velocity, with its value possibly influenced by mesh size and its spatial distribution. The indicator function serves as a fundamental element in unraveling the interplay between inner and outer layers in wall-bounded flows, including the overlap region, and forms a critical underpinning in turbulence modeling. Our investigation reveals a sensitivity of the indicator function on the utilized mesh distributions, prompting inquiries into the conventional mesh-sizing paradigms for DNS applications., Comment: 6 pages conference preprint. arXiv admin note: substantial text overlap with arXiv:2311.05204

Indicator functions of the streamwise normal-stress profiles (NSP), based on careful differentiation of some of the best direct numerical simulations (DNS) data from channel and pipe flows, over the range $550

The symmetry-based turbulence theory has been used to derive new scaling laws for the streamwise velocity and temperature moments of arbitrary order. For this, it has been applied to an incompressible turbulent channel flow driven by a pressure gradient with a passive scalar equation coupled in. To derive the scaling laws, symmetries of the classical Navier-Stokes and the thermal energy equations have been used together with statistical symmetries, i.e. the statistical scaling and translation symmetries of the multi-point moment equations. Specifically, the multi-point moments are built on the instantaneous velocity and temperature fields other than in the classical approach, where moments are based on the fluctuations of these fields. With this instantaneous approach, a linear system of multi-point correlation equations has been obtained, which greatly simplifies the symmetry analysis. The scaling laws have been derived in the limit of zero viscosity and heat conduction, i.e. $Re_\tau \rightarrow \infty$ and $Pr > 1$, and apply in the centre of the channel, i.e. they represent a generalization of the deficit law so herewith extending the work of done on the velocity field. The scaling laws are all power laws, with the exponent of the high moments all depending exclusively on those of the first and second moments. To validate the new scaling laws, the data from a large number of DNS for different Reynolds and Prandtl numbers have been used. The results show a very high accuracy of the scaling laws to represent the DNS data. The statistical scaling symmetry of the multi-point moment equations, which characterizes intermittency, has been the key to the new results since it generates a constant in the exponent of the final scaling law. Most important, since this constant is independent of the order of the moments, it clearly indicates anomalous scaling., Comment: 31 pages, 4 figures

Direct Numerical Simulations (DNSs) are one of the most powerful tools for studying turbulent flows. Even if achievable Reynolds numbers are lower than those obtained with experimental means, there is a clear advantage since the entire velocity field is known, and any desired quantity can be evaluated. This also includes the computation of derivatives of all relevant terms. One such derivative provides the indicator function, which is the product of wall distance by wall-normal derivative of the mean streamwise velocity. This derivative may depend on mesh spacing and distribution. However, it is extremely affected by the convergence of the simulation. The indicator function is a cornerstone to understanding inner and outer interactions in wall-bounded flows and describing the overlap region between them. We find a clear dependence of this indicator function on mesh distributions we examined, raising questions about classical mesh and convergence requirements for DNS and achievable accuracy. Within the framework of the logarithmic plus linear overlap region, coupled with a parametric study of channel flows and some pipe flows, sensitivities of extracted overlap parameters are examined, and a path is revealed to establishing their high-$Re_\tau$ or near-asymptotic values at modest Reynolds numbers accessible by high-quality DNS with reasonable ``cost''.

Despite its great scientific and technological importance, wall-bounded turbulence is an unresolved problem in classical physics that requires new perspectives to be tackled. One of the key strategies has been to study interactions among the energy-containing coherent structures in the flow. Such interactions are explored in this study for the first time using an explainable deep-learning method. The instantaneous velocity field obtained from a turbulent channel flow simulation is used to predict the velocity field in time through a U-net architecture. Based on the predicted flow, we assess the importance of each structure for this prediction using the game-theoretic algorithm of SHapley Additive exPlanations (SHAP). This work provides results in agreement with previous observations in the literature and extends them by revealing that the most important structures in the flow are not necessarily the ones with the highest contribution to the Reynolds shear stress. We also apply the method to an experimental database, where we can identify completely new structures based on their importance score. This framework has the potential to shed light on numerous fundamental phenomena of wall-bounded turbulence, including novel strategies for flow control.

The 2030 Agenda of the United Nations (UN) revolves around the Sustainable Development Goals (SDGs). A critical step towards that objective is identifying whether scientific production aligns with the SDGs' achievement. To assess this, funders and research managers need to manually estimate the impact of their funding agenda on the SDGs, focusing on accuracy, scalability, and objectiveness. With this objective in mind, in this work, we develop ASDG, an easy-to-use artificial-intelligence (AI)-based model for automatically identifying the potential impact of scientific papers on the UN SDGs. As a demonstrator of ASDG, we analyze the alignment of recent aerospace publications with the SDGs. The Aerospace data set analyzed in this paper consists of approximately 820,000 papers published in English from 2011 to 2020 and indexed in the Scopus database. The most-contributed SDGs are 7 (on clean energy), 9 (on industry), 11 (on sustainable cities) and 13 (on climate action). The establishment of the SDGs by the UN in the middle of the 2010 decade did not significantly affect the data. However, we find clear discrepancies among countries, likely indicative of different priorities. Also, different trends can be seen in the most and least cited papers, with clear differences in some SDGs. Finally, the number of abstracts the code cannot identify is decreasing with time, possibly showing the scientific community's awareness of SDG.

Turbulence is a complex phenomenon that has a chaotic nature with multiple spatio-temporal scales, making predictions of turbulent flows a challenging topic. Nowadays, an abundance of high-fidelity databases can be generated by experimental measurements and numerical simulations, but obtaining such accurate data in full-scale applications is currently not possible. This motivates utilising deep learning on subsets of the available data to reduce the required cost of reconstructing the full flow in such full-scale applications. Here, we develop a generative-adversarial-network (GAN)-based model to reconstruct the three-dimensional velocity fields from flow data represented by a cross-plane of unpaired two-dimensional velocity observations. The model could successfully reconstruct the flow fields with accurate flow structures, statistics and spectra. The results indicate that our model can be successfully utilised for reconstructing three-dimensional flows from two-dimensional experimental measurements. Consequently, a remarkable reduction of the experimental setup cost can be achieved., Comment: 21 pages, 5 figures, Nature - computational science (submitted)

High-fidelity large-eddy simulations of the flow around two rectangular obstacles are carried out at a Reynolds number of 10,000 based on the free-stream velocity and the obstacle height. The incoming flow is a developed turbulent boundary layer. Mean-velocity components, turbulence fluctuations, and the terms of the turbulent-kinetic-energy budget are analyzed for three flow regimes: skimming flow, wake interference, and isolated roughness. Three regions are identified where the flow undergoes the most significant changes: the first obstacle's wake, the region in front of the second obstacle, and that around the second obstacle. In the skimming-flow case, turbulence activity in the cavity between the obstacles is limited and mainly occurs in a small region in front of the second obstacle. In the wake-interference case, there is a strong interaction between the free-stream flow that penetrates the cavity and the wake of the first obstacle. This interaction results in more intense turbulent fluctuations between the obstacles. In the isolated-roughness case, the wake of the first obstacle is in good agreement with that of an isolated obstacle. Separation bubbles with strong turbulent fluctuations appear around the second obstacle.

Studying and interpreting the different flow patterns present in urban areas is becoming essential since they help develop new approaches to fight climate change through an improved understanding of the dynamics of the pollutants in urban environments. This study uses higher order dynamic mode decomposition (HODMD) to analyze a high-fidelity database of the turbulent flow in various simplified urban environments. The geometry simulated consists of two buildings separated by a certain distance. Three different cases have been studied, corresponding to the three different regimes identified in the bibliography. We recognize the characteristics of the well-known arch vortex forming on the leeward side of the first building and document possible generation and destruction mechanisms of this vortex based on the resulting temporal modes. These so-called vortex-generating and vortex-breaking modes are further analyzed via proper-orthogonal decomposition. We show that the arch vortex plays a prominent role in the dispersion of pollutants in urban environments, where its generation leads to an increase in their concentration; therefore, the reported mechanisms are of extreme importance in the context of urban sustainability.

Understanding flow structures in urban areas is widely recognized as a challenging concern due to its effect on urban development, air quality, and pollutant dispersion. In this study, state-of-the-art data-driven methods for modal analysis of simplified urban flows are used to study the dominant flow processes in these environments. Higher order dynamic mode decomposition (HODMD), a highly-efficient method to analyze turbulent flows, is used together with traditional techniques such as proper-orthogonal decomposition (POD) to analyze high-fidelity simulation data of a simplified urban environment. Furthermore, the spatio-temporal Koopman decomposition (STKD) will be applied to the temporal modes obtained with HODMD to perform spatial analysis. The flow interaction within the canopy influences the flow structures, particularly the arch vortex. The latter is a vortical structure generally found downstream of wall-mounted obstacles, which is generated as a consequence of flow separation. Therefore, the main objective of the present study is to characterize the mechanisms that promote these phenomena in urban areas with different geometries. Remarkably, among all the vortical structures identified by the HODMD algorithm, low- and high-frequency modes are classified according to their relation with the arch vortex. They are referred to as vortex-generator and vortex-breaker modes, respectively. This classification implies that one of the processes driving the formation and destruction of major vortical structures in between the buildings is the interaction between low- and high-frequency structures. The high energy revealed by the POD for the vortex-breaker modes points to this destruction process as the mechanism driving the flow dynamics.

The sustainability of urban environments is an increasingly relevant problem. Air pollution plays a key role in the degradation of the environment as well as the health of the citizens exposed to it. In this chapter we provide a review of the methods available to model air pollution, focusing on the application of machine-learning methods. In fact, machine-learning methods have proved to importantly increase the accuracy of traditional air-pollution approaches while limiting the development cost of the models. Machine-learning tools have opened new approaches to study air pollution, such as flow-dynamics modelling or remote-sensing methodologies.

We propose a deep probabilistic-neural-network architecture for learning a minimal and near-orthogonal set of non-linear modes from high-fidelity turbulent-flow-field data useful for flow analysis, reduced-order modeling, and flow control. Our approach is based on $\beta$-variational autoencoders ($\beta$-VAEs) and convolutional neural networks (CNNs), which allow us to extract non-linear modes from multi-scale turbulent flows while encouraging the learning of independent latent variables and penalizing the size of the latent vector. Moreover, we introduce an algorithm for ordering VAE-based modes with respect to their contribution to the reconstruction. We apply this method for non-linear mode decomposition of the turbulent flow through a simplified urban environment, where the flow-field data is obtained based on well-resolved large-eddy simulations (LESs). We demonstrate that by constraining the shape of the latent space, it is possible to motivate the orthogonality and extract a set of parsimonious modes sufficient for high-quality reconstruction. Our results show the excellent performance of the method in the reconstruction against linear-theory-based decompositions. Moreover, we compare our method with available AE-based models. We show the ability of our approach in the extraction of near-orthogonal modes that may lead to interpretability.