Adaptive mesh refinement method for CFD applications

The main objective of this thesis is the development of an adaptive mesh refinement (AMR) algorithm for computational fluid dynamics simulations using hexahedral and tetrahedral meshes. This numerical methodology is applied in the context of large-eddy simulations (LES) of turbulent flows and direct numerical simulations (DNS) of interfacial flows, to bring new numerical research and physical insight. For the fluid dynamics simulations, the governing equations, the spatial discretization on unstructured grids and the numerical schemes for solving Navier-Stokes equations are presented. The equations follow a discretization by conservative finite-volume on collocated meshes. For the turbulent flows formulation, the spatial discretization preserves symmetry properties of the continuous differential operators and the time integration follows a self-adaptive strategy, which has been well tested on unstructured grids. Moreover, LES model consisting of a wall adapting local-eddy-viscosity within a variational multi-scale formulation is used for the applications showed in this thesis. For the two-phase flow formulation, a conservative level-set method is applied for capturing the interface between two fluids and is implemented with a variable density projection scheme to simulate incompressible two-phase flows on unstructured meshes. The AMR algorithm developed in this thesis is based on a quad/octree data structure and keeps a relation of 1:2 between levels of refinement. In the case of tetrahedral meshes, a geometrical criterion is followed to keep the quality metric of the mesh on a reasonable basis. The parallelization strategy consists mainly in the creation of mesh elements in each sub-domain and establishes a unique global identification number, to avoid duplicate elements. Load balance is assured at each AMR iteration to keep the parallel performance of the CFD code. Moreover, a mesh multiplication algorithm (MM) is reported to create large meshes, with different kind of mesh elements, but preserving the topology from a coarser original mesh. This thesis focuses on the study of turbulent flows and two-phase flows using an AMR framework. The cases studied for LES of turbulent flows applications are the flow around one and two separated square cylinders, and the flow around a simplified car model. In this context, a physics-based refinement criterion is developed, consisting of the residual velocity calculated from a multi-scale decomposition of the instantaneous velocity. This criteria ensures grid adaptation following the main vortical structures and giving enough mesh resolution on the zones of interest, i.e., flow separation, turbulent wakes, and vortex shedding. The cases studied for the two-phase flows are the DNS of 2D and 3D gravity-driven bubble, with a particular focus on the wobbling regime. A study of rising bubbles in the wobbling regime and the effect of dimensionless numbers on the dynamic behavior of the bubbles are presented. Moreover, the use of tetrahedral AMR is applied for the numerical simulation of gravity-driven bubbles in complex domains. On this topic, the methodology is validated on bubbles rising in cylindrical channels with different topology, where the study of these cases contributed to having new numerical research and physical insight in the development of a rising bubble with wall effects.El objetivo principal de esta tesis es el desarrollo de un algoritmo adaptativo de refinamiento de malla (AMR) para simulaciones de dinámica de fluidos computacional utilizando mallas hexaédricas y tetraédricas. Esta metodología numérica se aplica en el contexto de simulaciones Large-eddie (LES) de flujos turbulentos y simulaciones numéricas directas (DNS) de flujos interfaciales, para traer nuevas investigaciones numéricas y entendimiento físicas. Para las simulaciones de dinámica de fluidos, se presentan las ecuaciones governantes, la discretización espacial en mallas no estructuradas y los esquemas numéricos para resolver las ecuaciones de Navier-Stokes. Las ecuaciones siguen una discretización conservativa por volumenes finitos en mallas colocadas. Para la formulación de flujos turbulentos, la discretización espacial preserva las propiedades de simetría de los operadores diferenciales continuos y la integración de tiempo sigue una estrategia autoadaptativa, que ha sido bien probada en mallas no estructuradas. Además, para las aplicaciones que se muestran en esta tesis, se utiliza el modelo LES que consiste en una viscosidad local que se adapta a la pared dentro de una formulación multiescala variable. Para la formulación de flujo de dos fases, se aplica un método de conjunto de niveles conservador para capturar la interfaz entre dos fluidos y se implementa con un esquema de proyección de densidad variable para simular flujos de dos fases incompresibles en mallas no estructuradas. El algoritmo AMR desarrollado en esta tesis se basa en una estructura de datos de quad / octree y mantiene una relación de 1: 2 entre los niveles de refinamiento. En el caso de las mallas tetraédricas, se sigue un criterio geométrico para mantener la calidad de la malla en una base razonable. La estrategia de paralelización consiste principalmente en la creación de elementos de malla en cada subdominio y establece un número de identificación global único, para evitar elementos duplicados. El equilibrio de carga está asegurado en cada iteración de AMR para mantener el rendimiento paralelo del código CFD. Además, se ha desarrollado un algoritmo de multiplicación de malla (MM) para crear mallas grandes, con diferentes tipos de elementos de malla, pero preservando la topología de una malla original más pequeña. Esta tesis se centra en el estudio de flujos turbulentos y flujos de dos fases utilizando un marco AMR. Los casos estudiados para aplicaciones de LES de flujos turbulentos son el flujo alrededor de uno y dos cilindros separados de sección cuadrada, y el flujo alrededor de un modelo de automóvil simplificado. En este contexto, se desarrolla un criterio de refinamiento basado en la física, que consiste en la velocidad residual calculada a partir de una descomposición de escala múltiple de la velocidad instantánea. Este criterio garantiza la adaptación de la malla siguiendo las estructuras vorticales principales y proporcionando una resolución de malla suficiente en las zonas de interés, es decir, separación de flujo, estelas turbulentas y desprendimiento de vórtices. Los casos estudiados para los flujos de dos fases son el DNS de la burbuja impulsada por la gravedad en 2D y 3D, con un enfoque particular en el régimen de oscilación. Además, el uso de AMR tetraédrico se aplica para la simulación numérica de burbujas impulsadas por la gravedad en dominios complejos. En este tema, la metodología se valida en burbujas que ascienden en canales cilíndricos con topología diferente, donde el estudio de estos casos contribuyó a tener una nueva investigación numérica y una visión física en el desarrollo de una burbuja con efectos de pared

An efficient Adaptive Mesh Refinement (AMR) algorithm for the Discontinuous Galerkin method: Applications for the computation of compressible two-phase flows

We present an Adaptive Mesh Refinement (AMR) method suitable for hybrid unstructured meshes that allows for local refinement and de-refinement of the computational grid during the evolution of the flow. The adaptive implementation of the Discontinuous Galerkin (DG) method introduced in this work (ForestDG) is based on a topological representation of the computational mesh by a hierarchical structure consisting of oct- quad- and binary trees. Adaptive mesh refinement (h-refinement) enables us to increase the spatial resolution of the computational mesh in the vicinity of the points of interest such as interfaces, geometrical features, or flow discontinuities. The local increase in the expansion order (p-refinement) at areas of high strain rates or vorticity magnitude results in an increase of the order of accuracy in the region of shear layers and vortices. A graph of unitarian-trees, representing hexahedral, prismatic and tetrahedral elements is used for the representation of the initial domain. The ancestral elements of the mesh can be split into self-similar elements allowing each tree to grow branches to an arbitrary level of refinement. The connectivity of the elements, their genealogy and their partitioning are described by linked lists of pointers. An explicit calculation of these relations, presented in this paper, facilitates the on-the-fly splitting, merging and repartitioning of the computational mesh by rearranging the links of each node of the tree with a minimal computational overhead. The modal basis used in the DG implementation facilitates the mapping of the fluxes across the non conformal faces. The AMR methodology is presented and assessed using a series of inviscid and viscous test cases. Also, the AMR methodology is used for the modelling of the interaction between droplets and the carrier phase in a two-phase flow. This approach is applied to the analysis of a spray injected into a chamber of quiescent air, using the Eulerian–Lagrangian approach. This enables us to refine the computational mesh in the vicinity of the droplet parcels and accurately resolve the coupling between the two phases

Numerical and Experimental Analysis of Transient Flow in Roots Blower

The performance of rotary positive displacement machines highly depends on the operational clearances. It is widely believed that computational fluid dynamics (CFD) can help understanding internal leakage flows. Developments of grid generating tools for analysis of leakage flows by CFD in rotary positive displacement machines have not yet been fully validated. Roots blower is a good representative of positive displacement machines and as such is convenient for optical access in order to analyse internal flows. The experimental investigation of flow in optical roots blower by phase-locked PIV (Particle Image Velocimetry) performed in the Centre for Compressor Technology at City, University of London provided the velocity field suitable for validation of the simulation model. This paper shows the results of the three-dimensional CFD transient simulation model of a Roots blower with the dynamic numerical grids generated by SCORG and flow solution solved in ANSYS CFX flow solver to obtain internal flow patterns. The velocity fields obtained by simulation agree qualitatively with the experimental results and show the correct main flow features in the working chamber. There are some differences in the velocity magnitude and vortex distribution. The flow field in roots blower is highly turbulent and three-dimensional. The axial clearances should be included, and the axial grids should be refined in the simulation method. The paper outlines some directions for future simulation and experimental work. The work described in this paper is a part of the large project set to evaluate characteristics of the internal flow in rotary positive displacement machines and to characterize leakage flow

Simulating liquids on dynamically warping grids

We introduce dynamically warping grids for adaptive liquid simulation. Our primary contributions are a strategy for dynamically deforming regular grids over the course of a simulation and a method for efficiently utilizing these deforming grids for liquid simulation. Prior work has shown that unstructured grids are very effective for adaptive fluid simulations. However, unstructured grids often lead to complicated implementations and a poor cache hit rate due to inconsistent memory access. Regular grids, on the other hand, provide a fast, fixed memory access pattern and straightforward implementation. Our method combines the advantages of both: we leverage the simplicity of regular grids while still achieving practical and controllable spatial adaptivity. We demonstrate that our method enables adaptive simulations that are fast, flexible, and robust to null-space issues. At the same time, our method is simple to implement and takes advantage of existing highly-tuned algorithms

A coupled level set and volume of fluid method for automotive exterior water management applications

Motivated by the need for practical, high fidelity, simulation of water over surface features of road vehicles a Coupled Level Set Volume of Fluid (CLSVOF) method has been implemented into a general purpose CFD code. It has been implemented such that it can be used with unstructured and non-orthogonal meshes. The interface reconstruction step needed for CLSVOF has been implemented using an iterative ‘clipping and capping’ algorithm for arbitrary cell shapes and a reinitialisation algorithm suitable for unstructured meshes is also presented. Successful verification tests of interface capturing on orthogonal and tetrahedral meshes are presented. Two macroscopic contact angle models have been implemented and the method is seen to give very good agreement with experimental data for a droplet impinging on a flat plate for both orthogonal and non-orthogonal meshes. Finally the flow of a droplet over a round edged channel is simulated in order to demonstrate the ability of the method developed to simulate surface flows over the sort of curved geometry that makes the use of a non-orthogonal grid desirable

Mesh dependence analysis for simulating unsteady cavitation around a plane convex hydrofoil

The mesh significantly influences the quality of the numerical results in computational fluid dynamics (CFD) problems due to its correlation with turbulence models. In the present paper a structured mesh was generated using Gmsh and two semi-structured meshes were generated using snappyHexMesh to determine the suitable mesh distribution for simulating unsteady cavitation around a plane convex hydrofoil. The numerical simulation was conducted by using the software OpenFOAM with the k-¿ SST SAS turbulence model and the ZwartGerber-Belamri (ZGB) cavitation model. Besides, the experimental results obtained by the Laboratory for Hydraulic Machines of École Polytechnique Fédérale de Lausanne were used to validate the numerical results. The results showed that both structured and semi-structured meshes predicted the cavitation pattern and the maximum cavity length. The semi-structured mesh with suitable refinement reproduced in detail the dynamic behavior of unsteady cavitation, while the structured mesh efficiently reproduced the phenomenon.Peer ReviewedObjectius de Desenvolupament Sostenible::7 - Energia Assequible i No ContaminantPostprint (published version