No Arabic abstract
Huge data advent in high-performance computing (HPC) applications such as fluid flow simulations usually hinders the interactive processing and exploration of simulation results. Such an interactive data exploration not only allows scientiest to play with their data but also to visualise huge (distributed) data sets in both an efficient and easy way. Therefore, we propose an HPC data exploration service based on a sliding window concept, that enables researches to access remote data (available on a supercomputer or cluster) during simulation runtime without exceeding any bandwidth limitations between the HPC back-end and the user front-end.
The Python package fluidsim is introduced in this article as an extensible framework for Computational Fluid Mechanics (CFD) solvers. It is developed as a part of FluidDyn project (Augier et al., 2018), an effort to promote open-source and open-science collaboration within fluid mechanics community and intended for both educational as well as research purposes. Solvers in fluidsim are scalable, High-Performance Computing (HPC) codes which are powered under the hood by the rich, scientific Python ecosystem and the Application Programming Interfaces (API) provided by fluiddyn and fluidfft packages (Mohanan et al., 2018). The present article describes the design aspects of fluidsim, viz. use of Python as the main language; focus on the ease of use, reuse and maintenance of the code without compromising performance. The implementation details including optimization methods, modular organization of features and object-oriented approach of using classes to implement solvers are also briefly explained. Currently, fluidsim includes solvers for a variety of physical problems using different numerical methods (including finite-difference methods). However, this metapaper shall dwell only on the implementation and performance of its pseudo-spectral solvers, in particular the two- and three-dimensional Navier-Stokes solvers. We investigate the performance and scalability of fluidsim in a state of the art HPC cluster. Three similar pseudo-spectral CFD codes based on Python (Dedalus, SpectralDNS) and Fortran (NS3D) are presented and qualitatively and quantitatively compared to fluidsim. The source code is hosted at Bitbucket as a Mercurial repository bitbucket.org/fluiddyn/fluidsim and the documentation generated using Sphinx can be read online at fluidsim.readthedocs.io.
We present hidden fluid mechanics (HFM), a physics informed deep learning framework capable of encoding an important class of physical laws governing fluid motions, namely the Navier-Stokes equations. In particular, we seek to leverage the underlying conservation laws (i.e., for mass, momentum, and energy) to infer hidden quantities of interest such as velocity and pressure fields merely from spatio-temporal visualizations of a passive scaler (e.g., dye or smoke), transported in arbitrarily complex domains (e.g., in human arteries or brain aneurysms). Our approach towards solving the aforementioned data assimilation problem is unique as we design an algorithm that is agnostic to the geometry or the initial and boundary conditions. This makes HFM highly flexible in choosing the spatio-temporal domain of interest for data acquisition as well as subsequent training and predictions. Consequently, the predictions made by HFM are among those cases where a pure machine learning strategy or a mere scientific computing approach simply cannot reproduce. The proposed algorithm achieves accurate predictions of the pressure and velocity fields in both two and three dimensional flows for several benchmark problems motivated by real-world applications. Our results demonstrate that this relatively simple methodology can be used in physical and biomedical problems to extract valuable quantitative information (e.g., lift and drag forces or wall shear stresses in arteries) for which direct measurements may not be possible.
The influence of turbulent effects on a fluid flow through a (pseudo) porous media is studied by numerically solving the set of Reynolds-averaged Navier-Stokes equations with the $kappa$-$epsilon$ model for turbulence. The spatial domains are two-dimensional rectangular grids with different {it porosities} obtained by the random placing of rigid obstacles. The objective of the simulations is to access the behavior of the generalized friction factor with varying Reynolds number. A good agreement with the Forchheimers equation is observed. The flow distribution at both low and high Reynolds conditions is also analyzed.
This paper deals with simulation of flow and transport in porous media such as transport of groundwater contaminants. We first discuss how macro scale equations are derived and which terms have to be closed by models. The transport of tracers is strongly influenced by pore scale velocity structure and large scale inhomogeneities in the permeability field. The velocity structure on the pore scale is investigated by direct numerical simulations of the 3D velocity field in a random sphere pack. The velocity probability density functions are strongly skewed, including some negative velocities. The large probability for very small velocities might be the reason for non-Fickian dispersion in the initial phase of contaminant transport. We present a method to determine large scale distributions of the permeability field from point-wise velocity measurements. The adjoint-based optimisation algorithm delivers fully satisfying agreement between input and estimated permeability fields. Finally numerical methods for convection dominated tracer transports are investigated from a theoretical point of view. It is shown that high order Finite Element Methods can reduce or even eliminate non-physical oscillations in the solution without introducing additional numerical diffusivity.
This work presents an efficient numerical method based on spectral expansions for simulation of heat and moisture diffusive transfers through multilayered porous materials. Traditionally, by using the finite-difference approach, the problem is discretized in time and space domains (Method of lines) to obtain a large system of coupled Ordinary Differential Equations (ODEs), which is computationally expensive. To avoid such a cost, this paper proposes a reduced-order method that is faster and accurate, using a much smaller system of ODEs. To demonstrate the benefits of this approach, tree case studies are presented. The first one considers nonlinear heat and moisture transfer through one material layer. The second case - highly nonlinear - imposes a high moisture content gradient - simulating a rain like condition - over a two-layered domain, while the last one compares the numerical prediction against experimental data for validation purposes. Results show how the nonlinearities and the interface between materials are easily and naturally treated with the spectral reduced-order method. Concerning the reliability part, predictions show a good agreement with experimental results, which confirm robustness, calculation efficiency and high accuracy of the proposed approach for predicting the coupled heat and moisture transfer through porous materials.