Solid-fluid two-phase flows, where the solid volume fraction is large either by geometry or by population (as in slurry flows), are ubiquitous in nature and industry. The interaction between the fluid and the suspended solids, in such flows, are too strongly coupled rendering the assumption of a single-way interaction (flow influences particle motion alone but not vice-versa) invalid and inaccurate. Most commercial flow solvers do not account for twoway interactions between fluid and immersed solids. The current state-of-art is restricted to two-way coupling between spherical particles (of very small diameters, such that the particlediameter to the characteristic flow domain length scale ratio is less than 0.01) and flow. These solvers are not suitable for solving several industrial slurry flow problems such as those of hydrates which is crucial to the oil-gas industry and rheology of slurries, flows in highly constrained geometries like microchannels or sessile drops that are laden with micro-PIV beads at concentrations significant for two-way interactions to become prominent. It is therefore necessary to develop direct numerical simulation flow solvers employing rigorous two-way coupling in order to accurately characterise the flow profiles between large immersed solids and fluid. It is necessary that such a solution takes into account the full 3D governing equations of flow (Navier-Stokes and continuity equations), solid translation (Newton’s second law) and solid rotation (equation of angular momentum) while simultaneously enabling interaction at every time step between the forces in the fluid and solid domains. This thesis concerns with development and rigorous validation of a 3D solid-fluid solver based on a novel variant of immersed-boundary method (IBM). The solver takes into account full two-way fluid-solid interaction with 6 degrees-of-freedom (6DOF). The solid motion solver is seamlessly integrated into the Gerris flow solver hence called Gerris Immersed Solid Solver (GISS). The IBM developed treats both fluid and solid in the manner of “fluid fraction” such that any number of immersed solids of arbitrary geometry can be realised. Our IBM method also allows transient local mesh adaption in the fluid domain around the moving solid boundary, thereby avoiding problems caused by the mesh skewness (as seen in common mesh-adaption algorithms) and significantly improves the simulation efficiency. The solver is rigorously validated at levels of increasing complexity against theory and experiment at low to moderate flow Reynolds number. At low Reynolds numbers (Re 1) these include: the drag force and terminal settling velocities of spherical bodies (validating translational degrees of freedom), Jeffrey’s orbits tracked by elliptical solids under shear flow (validating rotational and translational degrees of freedom) and hydrodynamic interaction between a solid and wall. Studies are also carried out to understand hydrodynamic interaction between multiple solid bodies under shear flow. It is found that initial distance between bodies is crucial towards the nature of hydrodynamic interaction between them: at a distance smaller than a critical value the solid bodies cluster together (hydrodynamic attraction) and at a distance greater than this value the solid bodies travel away from each other (hydrodynamic repulsion). At moderately high flow rates (Re O(100)), the solver is validated against migratory motion of an eccentrically placed solid sphere in Poisuelle flow. Under inviscid conditions (at very high Reynolds number) the solver is validated against chaotic motion of an asymmetric solid body. These validations not only give us confidence but also demonstrate the versatility of the GISS towards tackling complex solid-fluid flows. This work demonstrates the first important step towards ultra-high resolution direct numerical simulations of solid-fluid flows. The GISS will be available as opensource code from February 2015.

This thesis presents a framework for simulating the fluid dynamical behavior of complex moving boundary problems in a high-performance computing environment. The framework is implemented in the pELAFINT3D software package. Moving boundaries are evolved in a seamless fashion through the use of distributed narrow band level set methods and the effect of moving boundaries is incorporated into the flow solution by a novel Cartesian grid method. The proposed Cartesian grid approach builds on the concept of a ghost fluid method where boundary conditions are applied through least-squares polynomial extrapolations. The method is hybridized such that computational cells adjacent to moving boundaries change discretization schemes smoothly in time to avoid the introduction of strong oscillations in the pressure field. The hybridization is shown to have minimal effect on accuracy while significantly suppressing pressure oscillations. The computational capability of the Cartesian grid approach is enhanced with a massively parallel adaptive meshing algorithm. Local mesh enrichment is effected through the use of octree refinement, and a scalable mesh pruning algorithm is used to reduce the memory footprint of the Cartesian grid for geometries which are not well bounded by a rectangular cuboid. The computational work is kept in a well-balanced state through the use of an adaptive repartitioning strategy. The numerical scheme is validated against many benchmark problems and the composite approach is demonstrated to work well on tens of thousands of computational cores. A simulation of the closure phase of a mechanical heart valve was carried out to demonstrate the ability of the pELAFINT3D package to compute high Reynolds number flows with complex moving boundaries and a wide disparity in length scales. Finally, a novel image-to-computation algorithm was implemented to demonstrate the flexibility the current method allows in designing new applications.

No description available.

[ES] La interacción de contacto entre sólidos deformables es uno de los fenómenos más complejos en el ámbito de la mecánica computacional. La resolución de este problema requiere de algoritmos robustos para el tratamiento de no linealidades geométricas. El Método de Elementos Finitos (MEF) es uno de los más utilizados para el diseño de componentes mecánicos, incluyendo la solución de problemas de contacto. En este método el coste asociado al proceso de discretización (generación de malla) está directamente vinculado a la definición del contorno a modelar, lo cual dificulta la introducción en la simulación de superficies complejas, como las superficies NURBS, cada vez más utilizadas en el diseño de componentes. Esta tesis está basada en el "Cartesian grid Finite Element Method" (cgFEM). En esta metodología, encuadrada en la categoría de métodos "Immersed Boundary", se extiende el problema a un dominio de aproximación (cuyo mallado es sencillo de generar) que contiene al dominio de análisis completamente en su interior. Al desvincular la discretización de la definición del contorno del problema se reduce drásticamente el coste de generación de malla. Es por ello que el método cgFEM es una herramienta adecuada para la resolución de problemas en los que es necesario modificar la geometría múltiples veces, como el problema de optimización de forma o la simulación de desgaste. El método cgFEM permite también crear de manera automática y eficiente modelos de Elementos Finitos a partir de imágenes médicas. La introducción de restricciones de contacto habilitaría la posibilidad de considerar los diferentes estados de integración implante-tejido en procesos de optimización personalizada de implantes. Así, en esta tesis se desarrolla una formulación para resolver problemas de contacto 3D con el método cgFEM, considerando tanto modelos de contacto sin fricción como problemas con rozamiento de Coulomb. La ausencia de nodos en el contorno en cgFEM impide la aplicación de métodos tradicionales para imponer las restricciones de contacto, por lo que se ha desarrollado una formulación estabilizada que hace uso de un campo de tensiones recuperado para asegurar la estabilidad del método. Para una mayor precisión de la solución, se ha introducido la definición analítica de las superficies en contacto en la formulación propuesta. Además, se propone la mejora de la robustez de la metodología cgFEM en dos aspectos: el control del mal condicionamiento del problema numérico mediante un método estabilizado, y la mejora del campo de tensiones recuperado, utilizado en el proceso de estimación de error. La metodología propuesta se ha validado a través de diversos ejemplos numéricos presentados en la tesis, mostrando el gran potencial de cgFEM en este tipo de problemas. / [CAT] La interacció de contacte entre sòlids deformables és un dels fenòmens més complexos en l'àmbit de la mecànica computacional. La resolució d'este problema requerix d'algoritmes robustos per al tractament de no linealitats geomètriques. El Mètode dels Elements Finits (MEF) és un dels més utilitzats per al disseny de components mecànics, incloent la solució de problemes de contacte. En este mètode el cost associat al procés de discretització (generació de malla) està directament vinculat a la definició del contorn a modelar, la qual cosa dificulta la introducció en la simulació de superfícies complexes, com les superfícies NURBS, cada vegada més utilitzades en el disseny de components. Esta tesi està basada en el "Cartesian grid Finite Element Method" (cgFEM). En esta metodologia, enquadrada en la categoria de mètodes "Immersed Boundary", s'estén el problema a un domini d'aproximació (el mallat del qual és senzill de generar) que conté al domini d'anàlisi completament en el seu interior. Al desvincular la discretització de la definició del contorn del problema es reduïx dràsticament el cost de generació de malla. És per això que el mètode cgFEM és una ferramenta adequada per a la resolució de problemes en què és necessari modificar la geometria múltiples vegades, com el problema d'optimització de forma o la simulació de desgast. El mètode cgFEM permet també crear de manera automàtica i eficient models d'Elements Finits a partir d'imatges mèdiques. La introducció de restriccions de contacte habilitaria la possibilitat de considerar els diferents estats d'integració implant-teixit en processos d'optimització personalitzada d'implants. Així, en esta tesi es desenvolupa una formulació per a resoldre problemes de contacte 3D amb el mètode cgFEM, considerant tant models de contacte sense fricció com a problemes amb fregament de Coulomb. L'absència de nodes en el contorn en cgFEM impedix l'aplicació de mètodes tradicionals per a imposar les restriccions de contacte, per la qual cosa s'ha desenvolupat una formulació estabilitzada que fa ús d'un camp de tensions recuperat per a assegurar l'estabilitat del mètode. Per a una millor precisió de la solució, s'ha introduït la definició analítica de les superfícies en contacte en la formulació proposada. A més, es proposa la millora de la robustesa de la metodologia cgFEM en dos aspectes: el control del mal condicionament del problema numèric per mitjà d'un mètode estabilitzat, i la millora del camp de tensions recuperat, utilitzat en el procés d'estimació d'error. La metodologia proposada s'ha validat a través de diversos exemples numèrics presentats en la tesi, mostrant el gran potencial de cgFEM en este tipus de problemes. / [EN] The contact interaction between elastic solids is one of the most complex phenomena in the computational mechanics research field. The solution of such problem requires robust algorithms to treat the geometrical non-linearities characteristic of the contact constrains. The Finite Element Method (FE) has become one of the most popular options for the mechanical components design, including the solution of contact problems. In this method the computational cost of the generation of the discretization (mesh generation) is directly related to the complexity of the analysis domain, namely its boundary. This complicates the introduction in the numerical simulations of complex surfaces (for example NURBS), which are being increasingly used in the CAD industry. This thesis is grounded on the Cartesian grid Finite Element Method (cgFEM). In this methodology, which belongs to the family of Immersed Boundary methods, the problem at hand is extended to an approximation domain which completely embeds the analysis domain, and its meshing is straightforward. The decoupling of the boundary definition and the discretization mesh results in a great reduction of the mesh generation's computational cost. Is for this reason that the cgFEM is a suitable tool for the solution of problems that require multiple geometry modifications, such as shape optimization problems or wear simulations. The cgFEM is also capable of automatically generating FE models from medical images without the intermediate step of generating CAD entities. The introduction of the contact interaction would open the possibility to consider different states of the union between implant and living tissue for the design of optimized implants, even in a patient-specific process. Hence, in this thesis a formulation for solving 3D contact problems with the cgFEM is presented, considering both frictionless and Coulomb's friction problems. The absence of nodes along the boundary in cgFEM prevents the enforcement of the contact constrains using the standard procedures. Thus, we develop a stabilized formulation that makes use of a recovered stress field, which ensures the stability of the method. The analytical definition of the contact surfaces (by means of NURBS) has been included in the proposed formulation in order to increase the accuracy of the solution. In addition, the robustness of the cgFEM methodology is increased in this thesis in two different aspects: the control of the numerical problem's ill-conditioning by means of a stabilized method, and the enhancement of the stress recovered field, which is used in the error estimation procedure. The proposed methodology has been validated through several numerical examples, showing the great potential of the cgFEM in these type of problems. / Navarro Jiménez, JM. (2019). Contact problem modelling using the Cartesian grid Finite Element Method [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/124348 / TESIS

No description available.

Prediction of the near wall region is crucial to the accuracy of turbulent flow computational fluid dynamics (CFD) simulation. However, sufficient near-wall resolution is often prohibitive for high Reynolds number flows with complex geometries, due to high memory and processing requirements. A common approach in these cases is to use wall functions to bridge the region from the first grid node to the wall. This thesis presents an alternative method that relaxes the near wall resolution requirement by solving one dimensional transport equations for velocity and turbulence across a locally defined subgrid contained within wall adjacent grid cells. The addition of the subgrid allows for wall adjacent primary grid sizes to vary arbitrarily from low-Re model sizing (y+ ~ 1) to wall function sizing without significant loss of accuracy or increase in computational cost.

Previous experimental and numerical studies focused on incompressible flow interactions with multi-scale fractal structures targeting the generation of turbulence at multiple scales. Depending on various flow conditions, it was found that these fractal structures are able to enhance mixing and scalar transport, and in some cases reduce flow generated sound in certain frequency ranges. The interaction of compressible flows with multi-scale fractal structures, however, did not receive attention as the focus was entirely on the incompressible regime. The objective of this study is to conduct large eddy simulations (LES) of flow interactions with a class of fractal plates in the compressible regime, and to extract and analyze different flow statistics in an attempt to determine the effect of compressibility. Immersed boundary methods (IBM) will be employed to overcome the difficulty of modeling the fractal structures via a bodyitted mesh, with adequate mesh resolution around small features of the fractal shapes.

Porosity is known to be one of the primary detrimental factors controlling fatigue life and total elongation of several cast alloy components. The two main aims of this work are to examine pore nucleation and growth effects for predicting gas microporosity and to study the physics of bifilm dynamics to gain understanding in the role of bifilms in producing defects and the mechanisms of defect creation. In the second chapter of this thesis, an innovative technique, based on the combination of a set of conservation equations that solves the transport phenomena during solidification at the macro-scale and the hydrogen diffusion into the pores at the micro-scale, was used to quantify the amount of gas microporosity in A356 alloy castings. The results were compared with published experimental data. In the reminder of this work, the Immersed Element-Free Galerkin method (IEFGM) is presented and it was used to study the physics of bifilm dynamics. The IEFGM is an extension of the Immersed Finite Element method (IFEM) developed by Zhang et al. [50] and it is an attractive technique for simulating FSI problems involving highly deformable bifilm-like solids.

A wide range of applications involve interface problems. In most of the cases, mathematical modeling of these interface problems leads to partial differential equations with non-smooth or discontinuous inputs and solutions, especially across material interfaces. Different numerical methods have been developed to solve these kinds of problems and handle the non-smooth behavior of the input data and/or the solution across the interface. The main focus of our work is the immersed finite element method to obtain optimal numerical solutions for interface problems. In this thesis, we present piecewise quadratic immersed finite element (IFE) spaces that are used with an immersed finite element (IFE) method with interior penalty (IP) for solving two-dimensional second-order elliptic interface problems without requiring the mesh to be aligned with the material interfaces. An analysis of the constructed IFE spaces and their dimensions is presented. Shape functions of Lagrange and hierarchical types are constructed for these spaces, and a proof for the existence is established. The interpolation errors in the proposed piecewise quadratic spaces yield optimal <i>O</i>(h³) and <i>O</i>(h²) convergence rates, respectively, in the L² and broken H¹ norms under mesh refinement. Furthermore, numerical results are presented to validate our theory and show the optimality of our quadratic IFE method. Our approach in this thesis is, first, to establish a theory for the simplified case of a linear interface. After that, we extend the framework to quadratic interfaces. We, then, describe a general procedure for handling arbitrary interfaces occurring in real physical practical applications and present computational examples showing the optimality of the proposed method. Furthermore, we investigate a general procedure for extending our quadratic IFE spaces to <i>p</i>-th degree and construct hierarchical shape functions for <i>p</i>=3. / Ph. D.

In this dissertation we prove the superconvergence of the minimal-dissipation local discontinuous Galerkin method for elliptic problems and construct optimal immersed finite element approximations and discontinuous immersed finite element methods for the Stokes interface problem. In the first part we present an error analysis for the minimal dissipation local discontinuous Galerkin method applied to a model elliptic problem on Cartesian meshes when polynomials of degree at most <i>k</i> and an appropriate approximation of the boundary condition are used. This special approximation allows us to achieve <i>k</i> + 1 order of convergence for both the potential and its gradient in the L² norm. Here we improve on existing estimates for the solution gradient by a factor &#8730;h. In the second part we present discontinuous immersed finite element (IFE) methods for the Stokes interface problem on Cartesian meshes that does not require the mesh to be aligned with the interface. As such, we allow unfitted meshes that are cut by the interface. Thus, elements may contain more than one fluid. On these unfitted meshes we construct an immersed Q₁/Q₀ finite element approximation that depends on the location of the interface. We discuss the basic features of the proposed Q₁/Q₀ IFE basis functions such as the unisolvent property. We present several numerical examples to demonstrate that the proposed IFE approximations applied to solve interface Stokes problems maintain the optimal approximation capability of their standard counterpart applied to solve the homogeneous Stokes problem. Similarly, we also show that discontinuous Galerkin IFE solutions of the Stokes interface problem maintain the optimal convergence rates in both L² and broken H¹ norms. Furthermore, we extend our method to solve the axisymmetric Stokes interface problem with a moving interface and test the proposed method by solving several benchmark problems from the literature. / Ph. D.