A dam-break flood simulation model in curvilinear coordinates

A dam-break flood model based on a contravariant integral form of the shallow water equations is presented. The numerical integration of the equations of motion is carried out by means of a finite volumefinite difference numerical scheme that involves an exact Riemann solver and which is based on a high-order WENO reconstruction procedure. An original scheme for the simulation of the wet front progress on the dry bed is adopted. The proposed model capacity to correctly simulate the wet front progress velocity is tested by numerically reproducing the dry bed dam-break problem. The model is adopted for the real case study of the Rio Fucino lake-dam collapse and subsequent flood wave propagation, downstream of the Campotosto reservoir (Italy)

FullSWOF_Paral: Comparison of two parallelization strategies (MPI and SKELGIS) on a software designed for hydrology applications

In this paper, we perform a comparison of two approaches for the parallelization of an existing, free software, FullSWOF 2D (http://www. univ-orleans.fr/mapmo/soft/FullSWOF/ that solves shallow water equations for applications in hydrology) based on a domain decomposition strategy. The first approach is based on the classical MPI library while the second approach uses Parallel Algorithmic Skeletons and more precisely a library named SkelGIS (Skeletons for Geographical Information Systems). The first results presented in this article show that the two approaches are similar in terms of performance and scalability. The two implementation strategies are however very different and we discuss the advantages of each one.Comment: 27 page

DNS of multiphase flows: study of atomization and free-surface phenomena

The present thesis focuses on the numerical analysis of some diverse physical set-up that involve the interaction of two -or three immiscible and incompressible phases. The simulations are carried out by means of finite-volume algorithms developed on the in-house Computational Fluid Dynamics platform TermoFluids, developed by the Heat and Mass Transfer Technological Center (CTTC). They are intended to give detailed insights on the physics of the analyzed phenomena by carrying out Direct Numerical Simulations (DNS). In the context of multiphase flows, namely, Computational Multi-Fluid Dynamics (CMFD) field, DNS means that all the interfacial and turbulent scales of the phenomenon must be fully resolved. In the Introduction, a general overview of the engineering applications and the computational methods related to multiphase flows is proposed. The various types of physics analyzed in this work and the numerical approaches applied here to carry out efficient simulations are introduced. In Chapter 2, a low-dissipation convection scheme for the stable discretization of multiphase flow by means of interface-capturing schemes is analyzed. The hybrid form of the convective operator proposed incorporates localized low-dispersion characteristics to limit the growth of spurious flow solutions. Moreover, in comparison to pure-dissipative schemes, the discretization aims at minimizing the differences in kinetic energy preservation with respect to the continuous governing equations. This property plays a fundamental role in the case of flows presenting significant levels of turbulence. The simulation of a turbulent 2D coaxial jet with the low-dissipation convection scheme demonstrates its capability of solving correctly the two-phase turbulent problems. In Chapter 3, all the work carried out on the simulation of two-phase flow with the aid of Adaptive Mesh Refinement (AMR) strategies is described. The model is globally addressed at improving the representation of interfacial and turbulent scales in general multiphase flows. It is first applied to the simulation of simple multiphase phenomena, as 2D and 3D rising bubbles, to demonstrate the convergence of the method and the important computational savings in comparison to static mesh computations. However, its adoption becames essential in the simulation of instability and break-up phenomena, where the necessity of representing accurately the complex structures that appear at the interface, as ligaments and droplets, make the simulation particularly expensive in terms of computational cost. In Chapter 4, we analyze in detail the simulations of 3-D atomizing phenomena, including the coaxial jet case, characterized by the parallel injection of high speed liquid and gas fluxes, and the liquid spray case, characterized by the injection of a high speed liquid inside a still air chamber. In Chapter 5, an original single-phase scheme for the DNS of free-surface problems on 3-D unstructured meshes is presented. The scheme is based on a novel treatment of the interface for the deactivation of the light-phase, allowing an optimization of the classic two-phase model for the cases in which the influence of the lighter phase is negligible. Consequently, the model is particularly addressed at analyzing problems involving the movement of free-surfaces, as the evolution of waves on the sea, and their interaction with fixed and moving obstacles. Some practical cases of application are proposed, as the evaluation of stresses on an object due to the action of a dam-break event, and the interaction of linear waves with an oscillating water column device. In the same Chapter we describe the procedure to couple the single-phase model to the Immersed Boundary Method. The method is aimed at representing the interaction of a solid moving with prescribed velocity and the two-phase flow. The most significant example consists in the simulation of a sliding wedge into a liquid basin.Esta tesis se focaliza en la simulación numérica de algunos set-up físicos que involucran la interacción entre dos o tres fluidos incompresibles y immiscibles. Las simulaciones se realizan por medio de algoritmos de volúmenes-finitos desarrollados en la plataforma propia de Fluido-Dinámica Computacional (CFD) denominada TermoFluids, desarrollada en el Centro Tecnologico de Trasferencia de Calor (CTTC). Las simulaciones quieren estudiar en detalle la física de los fenómenos analizados, realizando su Simulación Numérica Directa (DNS). En el contexto de los flujos multifase, DNS significa que todas las escalas interfaciales y turbulentas del fenómeno han de ser totalmente resueltas. En la Introducción, se propone una panorámica general de las aplicaciones de ingeniería y de los métodos computacionales relacionados con flujos multifases. Se introducen los varios tipos de física analizados en este trabajo y las estrategias numéricas aplicadas aquí para efectuar su simulación de manera eficiente. En el Capitulo 2 se analiza un esquema convectivo de baja-disipación para la discretization de flujo multifase por medio de métodos de interface-capturing. La forma híbrida del operador convectivo propuesto incorpora la característica de una baja dispersión localizada, focalizada en limitar el crecimiento de soluciones numéricas espurias. Además, en comparación con métodos disipativos puros, la discretización apunta a minimizar las diferencias en la conservación de energía cinética en respeto a las ecuaciones continuas que gobiernan el flujo. Esta propiedad juega un papel fundamental en el caso de flujo caracterizado por un alto nivel de turbulencia. La simulación de un jet 2D coaxial turbulento con el método convectivo de baja disipación demuestra su capacidad de resolver correctamente un flujo de dos fases turbulentos. En el Capitulo 3 se reporta todo el trabajo realizado sobre la simulación de flujo multifase con el auxilio de técnicas de refinamiento adaptativo de malla (AMR). El modelo es globalmente dirigido a la mejora de la representación de las escalas turbulentas y interfaciales en flujos multifases en general. Se aplica inicialmente a la simulación de flujos sencillos, como unos casos de burbujas flotantes 2D y 3D, demostrando la convergencia del método y los importantes ahorros computacional en comparación con los cálculos de mallas estáticas. La adopción de la técnica se hace esencial en la simulación de fenómenos de inestabilidad y de ruptura, donde la necesidad de representar sacramentalmente las estructuras complejas que aparecen en la interfaz, como ligamentos o pequeñas gotas, hacen que la simulación sea particularmente pesada en términos de coste computacional. En el Capitulo 4 se reportan en detalle las simulaciones de fenómenos de atomización 3D. Esas incluyen el caso del jet coaxial, caracterizado por la inyección paralela de flujos de aire y liquido de altas velocidades, y el caso del spray liquido, que consiste en la inyección de un liquido dentro de una cámara de aire. En el Capitulo 5 se presenta un esquema de single-phase original, para el DNS de problemas de superficie libre en mallas 3D no-estructuradas. El esquema se basa en un nuevo tratamiento de la interfase para la desactivación de la fase ligera, permitiendo la optimización del solver clásico de dos fases para los casos en que la influencia de la fase mas ligera sea despreciable. En consecuencia, el modelo es particularmente indicado para la análisis de problemas que involucran el movimiento de superficies libres, como la evolución de olas en la superficie marina y su interacción con obstáculos fijos o muebles. Se proponen algunos casos prácticos de aplicación, como la evaluación de las fuerzas sobre un objeto debidos a un episodio de dam-break, o el estudio de las olas generadas por el impacto de un solido deslizante (representado integrando la tecnica de Immersed Boundary con el presente metodo de single-phase) con un embalse de agua.Postprint (published version

DNS of multiphase flows: study of atomization and free-surface phenomena

The present thesis focuses on the numerical analysis of some diverse physical set-up that involve the interaction of two -or three immiscible and incompressible phases. The simulations are carried out by means of finite-volume algorithms developed on the in-house Computational Fluid Dynamics platform TermoFluids, developed by the Heat and Mass Transfer Technological Center (CTTC). They are intended to give detailed insights on the physics of the analyzed phenomena by carrying out Direct Numerical Simulations (DNS). In the context of multiphase flows, namely, Computational Multi-Fluid Dynamics (CMFD) field, DNS means that all the interfacial and turbulent scales of the phenomenon must be fully resolved. In the Introduction, a general overview of the engineering applications and the computational methods related to multiphase flows is proposed. The various types of physics analyzed in this work and the numerical approaches applied here to carry out efficient simulations are introduced. In Chapter 2, a low-dissipation convection scheme for the stable discretization of multiphase flow by means of interface-capturing schemes is analyzed. The hybrid form of the convective operator proposed incorporates localized low-dispersion characteristics to limit the growth of spurious flow solutions. Moreover, in comparison to pure-dissipative schemes, the discretization aims at minimizing the differences in kinetic energy preservation with respect to the continuous governing equations. This property plays a fundamental role in the case of flows presenting significant levels of turbulence. The simulation of a turbulent 2D coaxial jet with the low-dissipation convection scheme demonstrates its capability of solving correctly the two-phase turbulent problems. In Chapter 3, all the work carried out on the simulation of two-phase flow with the aid of Adaptive Mesh Refinement (AMR) strategies is described. The model is globally addressed at improving the representation of interfacial and turbulent scales in general multiphase flows. It is first applied to the simulation of simple multiphase phenomena, as 2D and 3D rising bubbles, to demonstrate the convergence of the method and the important computational savings in comparison to static mesh computations. However, its adoption becames essential in the simulation of instability and break-up phenomena, where the necessity of representing accurately the complex structures that appear at the interface, as ligaments and droplets, make the simulation particularly expensive in terms of computational cost. In Chapter 4, we analyze in detail the simulations of 3-D atomizing phenomena, including the coaxial jet case, characterized by the parallel injection of high speed liquid and gas fluxes, and the liquid spray case, characterized by the injection of a high speed liquid inside a still air chamber. In Chapter 5, an original single-phase scheme for the DNS of free-surface problems on 3-D unstructured meshes is presented. The scheme is based on a novel treatment of the interface for the deactivation of the light-phase, allowing an optimization of the classic two-phase model for the cases in which the influence of the lighter phase is negligible. Consequently, the model is particularly addressed at analyzing problems involving the movement of free-surfaces, as the evolution of waves on the sea, and their interaction with fixed and moving obstacles. Some practical cases of application are proposed, as the evaluation of stresses on an object due to the action of a dam-break event, and the interaction of linear waves with an oscillating water column device. In the same Chapter we describe the procedure to couple the single-phase model to the Immersed Boundary Method. The method is aimed at representing the interaction of a solid moving with prescribed velocity and the two-phase flow. The most significant example consists in the simulation of a sliding wedge into a liquid basin.Esta tesis se focaliza en la simulación numérica de algunos set-up físicos que involucran la interacción entre dos o tres fluidos incompresibles y immiscibles. Las simulaciones se realizan por medio de algoritmos de volúmenes-finitos desarrollados en la plataforma propia de Fluido-Dinámica Computacional (CFD) denominada TermoFluids, desarrollada en el Centro Tecnologico de Trasferencia de Calor (CTTC). Las simulaciones quieren estudiar en detalle la física de los fenómenos analizados, realizando su Simulación Numérica Directa (DNS). En el contexto de los flujos multifase, DNS significa que todas las escalas interfaciales y turbulentas del fenómeno han de ser totalmente resueltas. En la Introducción, se propone una panorámica general de las aplicaciones de ingeniería y de los métodos computacionales relacionados con flujos multifases. Se introducen los varios tipos de física analizados en este trabajo y las estrategias numéricas aplicadas aquí para efectuar su simulación de manera eficiente. En el Capitulo 2 se analiza un esquema convectivo de baja-disipación para la discretization de flujo multifase por medio de métodos de interface-capturing. La forma híbrida del operador convectivo propuesto incorpora la característica de una baja dispersión localizada, focalizada en limitar el crecimiento de soluciones numéricas espurias. Además, en comparación con métodos disipativos puros, la discretización apunta a minimizar las diferencias en la conservación de energía cinética en respeto a las ecuaciones continuas que gobiernan el flujo. Esta propiedad juega un papel fundamental en el caso de flujo caracterizado por un alto nivel de turbulencia. La simulación de un jet 2D coaxial turbulento con el método convectivo de baja disipación demuestra su capacidad de resolver correctamente un flujo de dos fases turbulentos. En el Capitulo 3 se reporta todo el trabajo realizado sobre la simulación de flujo multifase con el auxilio de técnicas de refinamiento adaptativo de malla (AMR). El modelo es globalmente dirigido a la mejora de la representación de las escalas turbulentas y interfaciales en flujos multifases en general. Se aplica inicialmente a la simulación de flujos sencillos, como unos casos de burbujas flotantes 2D y 3D, demostrando la convergencia del método y los importantes ahorros computacional en comparación con los cálculos de mallas estáticas. La adopción de la técnica se hace esencial en la simulación de fenómenos de inestabilidad y de ruptura, donde la necesidad de representar sacramentalmente las estructuras complejas que aparecen en la interfaz, como ligamentos o pequeñas gotas, hacen que la simulación sea particularmente pesada en términos de coste computacional. En el Capitulo 4 se reportan en detalle las simulaciones de fenómenos de atomización 3D. Esas incluyen el caso del jet coaxial, caracterizado por la inyección paralela de flujos de aire y liquido de altas velocidades, y el caso del spray liquido, que consiste en la inyección de un liquido dentro de una cámara de aire. En el Capitulo 5 se presenta un esquema de single-phase original, para el DNS de problemas de superficie libre en mallas 3D no-estructuradas. El esquema se basa en un nuevo tratamiento de la interfase para la desactivación de la fase ligera, permitiendo la optimización del solver clásico de dos fases para los casos en que la influencia de la fase mas ligera sea despreciable. En consecuencia, el modelo es particularmente indicado para la análisis de problemas que involucran el movimiento de superficies libres, como la evolución de olas en la superficie marina y su interacción con obstáculos fijos o muebles. Se proponen algunos casos prácticos de aplicación, como la evaluación de las fuerzas sobre un objeto debidos a un episodio de dam-break, o el estudio de las olas generadas por el impacto de un solido deslizante (representado integrando la tecnica de Immersed Boundary con el presente metodo de single-phase) con un embalse de agua

Deep Fluids: A Generative Network for Parameterized Fluid Simulations

This paper presents a novel generative model to synthesize fluid simulations from a set of reduced parameters. A convolutional neural network is trained on a collection of discrete, parameterizable fluid simulation velocity fields. Due to the capability of deep learning architectures to learn representative features of the data, our generative model is able to accurately approximate the training data set, while providing plausible interpolated in-betweens. The proposed generative model is optimized for fluids by a novel loss function that guarantees divergence-free velocity fields at all times. In addition, we demonstrate that we can handle complex parameterizations in reduced spaces, and advance simulations in time by integrating in the latent space with a second network. Our method models a wide variety of fluid behaviors, thus enabling applications such as fast construction of simulations, interpolation of fluids with different parameters, time re-sampling, latent space simulations, and compression of fluid simulation data. Reconstructed velocity fields are generated up to 700x faster than re-simulating the data with the underlying CPU solver, while achieving compression rates of up to 1300x.Comment: Computer Graphics Forum (Proceedings of EUROGRAPHICS 2019), additional materials: http://www.byungsoo.me/project/deep-fluids

A "well-balanced" finite volume scheme for blood flow simulation

We are interested in simulating blood flow in arteries with a one dimensional model. Thanks to recent developments in the analysis of hyperbolic system of conservation laws (in the Saint-Venant/ shallow water equations context) we will perform a simple finite volume scheme. We focus on conservation properties of this scheme which were not previously considered. To emphasize the necessity of this scheme, we present how a too simple numerical scheme may induce spurious flows when the basic static shape of the radius changes. On contrary, the proposed scheme is "well-balanced": it preserves equilibria of Q = 0. Then examples of analytical or linearized solutions with and without viscous damping are presented to validate the calculations. The influence of abrupt change of basic radius is emphasized in the case of an aneurism.Comment: 36 page

Modelling Shallow Water Flows by a High Resolution Riemann Solver

Shock-capturing methods originally developed for compressible gas dynamics has been applied to many other non-linear hyperbolic systems of conservation equations, like reactive flows, two-phase flow, porous media flow and finally shallow water flows, with applications to rivers, estuaries, dam break and in particular tsunami propagation resulting from earthquake and landslide. The basic ingredient of the proposed model for solving the shallow water equations is a finite volume Flux Vector Splitting Riemann Solver with a second order resolution scheme and implicit treatment of the source terms. The model includes entropy fix and step bed treatments, and finally shore line tracking, giving the ability to capture local discontinuities - like shock waves - and reducing numerical diffusion and unphysical viscosity effects which dominates in all finite difference methods. The numerical model has been validated in respect to different numerical test cases and comparisons with the exact solution of the Riemann problems are presented.JRC.G.7-Traceability and vulnerability assessmen