Tesis por compendio / [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. / [CA] 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]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/124348 / Compendio

Most realistic fluid flow problems are characterised by high Reynolds numbers and complex boundaries. Over the last ten years, immersed boundary methods (IBM) that are able to cope with realistic geometries have been applied to Lattice- Boltzmann methods (LBM). These methods, however, have normally been applied to low Reynolds number problems. In the present work, an iterative direct forcing IBM has been successfully coupled with a multi-domain cascaded LBM in order to investigate viscous flows around rigid, moving and wilfully deformed boundaries at a wide range of Reynolds numbers. The iterative force-correction immersed boundary method of (Zhang et al., 2016) has been selected due to the improved accuracy of the computation, while the cascaded LB formulation is used due to its superior stability at high Reynolds numbers. The coupling is shown to improve both the stability and numerical accuracy of the solution. The resulting solver has been applied to viscous flow (up to a Reynolds number of 100000) passed a NACA-0012 airfoil at a 10 degree angle of attack. Good agreement with results obtained using a body-fitted Navier-Stokes solver has been obtained. At moving or deformable boundary applications, emphasis should be given on the influence of the internal mass on the computation of the aerodynamic forces, focusing on deforming boundary motions where the rigid body approximation is no longer valid. Both the rigid body and the internal Lagrangian points approximations are examined. The resulting solver has been applied to viscous flows around an in-line oscillating cylinder, a pitching foil, a plunging SD7003 airfoil and a plunging and flapping NACA-0014 airfoil. Good agreement with experimental results and other numerical schemes has been obtained. It is shown that the internal Lagrangian points approximation accurately captures the internal mass effects in linear and angular motions, as well as in deforming motions, at Reynolds numbers up to 4 • 104. Finally, an expanded higher-order immersed boundary method which addresses two major drawbacks of the conventional IBM will be presented. First, an expanded velocity profile scheme has been developed, in order to compensate for the discontinuities caused by the gradient of the velocity across the boundary. Second, a numerical method derived from the Navier-Stokes equations in order to correct the pressure distribution across the boundary has been examined. The resulting hybrid solver has been applied to viscous flows around stationary and oscillating cylinders and examined the hovering flight of elliptical wings at low Reynolds numbers. It is shown that the proposed scheme smoothly expands the velocity profile across the boundary and increases the accuracy of the immersed boundary method. In addition, the pressure correction algorithm correctly expands the pressure profile across the boundary leading to very accurate pressure coefficient values along the boundary surface. The proposed numerical schemes are shown to be very efficient in terms of computational cost. The majority of the presented results are obtained within a few hours of CPU time on a 2.8 GHz Intel Core i7 MacBook Pro computer with a 16GB memory.

L'approche la plus courante en Mécanique des Fluides Numérique pour réaliser les simulations d'écoulement cardiovasculaire consiste à utiliser des méthodes numériques Body-fitted. Ces méthodes ont permis d'obtenir des simulations d'écoulement sanguin dans les artères qui sont précises et utiles. Toutefois, la génération du maillage body-fitted est une tâche qui demande beaucoup de temps et d'expertise à l'utilisateur.Les méthodes de Frontières Immergées sont des méthodes numériques alternatives qui ont l'avantage d'être plus simples d'emploi car elles ne requièrent aucune tâche de maillage de la part de l'utilisateur. Le travail présenté ici vise à évaluer le potentiel d'un méthode de Frontières Immergées à réaliser des simulations d'écoulement cardiovasculaire.Ce travail s'attache, dans un premier temps, à décrire les capacités de cette méthode numérique à rendre compte de l'imperméabilité et de la mobilité des parois sur des cas relativement simples mais représentatifs d'écoulements cardiovasculaires. Ensuite, des applications de la méthode à des cas d'écoulement cardiovasculaire plus complexes sont montrées. Il s'agira d'abord d'une simulation de l'écoulement dans un modèle rigide d'artère aorte. Puis, la simulation d'un écoulement à l'intérieur d'un ventricule cardiaque à paroi mobile sera montrée. / The most common approach in Computational Fluid Dynamics(CFD) for simulating blood flow into vessel is to make use of a body-fitted me-thod. This approach has lead to accurate and useful simulations of blood flowinto arteries. However, generation of the body-fitted grid is time consuming andrequires from the user an engineering knowledge.The Immersed Boundary Method has emerged as an alternate method whichdoes not require from the user any grid generation task. Simulations are done on astructured Cartesian grid which can be automatically generated. Here we addressthe question of the capability of an Immersed Boundary Method to cope withcardiovascular flow simulations.In particular, we assess the impermeable and moving properties of the wallwhen using the Immersed Boundary Method on simple but relevant vascular flowcases. Then, we show more complex and realistic cardiovascular flow simulations.The first application consists of blood flow simulation inside an aorta cross model.Then, the simulation of blood flow inside a cardiac ventricle with moving wall isshown.

We present a novel dimensionally split Cartesian cut cell method to compute inviscid, viscous and turbulent flows around rigid geometries. On a cut cell mesh, the existence of arbitrarily small boundary cells severely restricts the stable time step for an explicit numerical scheme. We solve this `small cell problem' when computing solutions for hyperbolic conservation laws by combining wave speed and geometric information to develop a novel stabilised cut cell flux. The convergence and stability of the developed technique are proved for the one-dimensional linear advection equation, while its multi-dimensional numerical performance is investigated through the computation of solutions to a number of test problems for the linear advection and Euler equations. This work was recently published in the Journal of Computational Physics (Gokhale et al., 2018). Subsequently, we develop the method further to be able to compute solutions for the compressible Navier-Stokes equations. The method is globally second order accurate in the L1 norm, fully conservative, and allows the use of time steps determined by the regular grid spacing. We provide a full description of the three-dimensional implementation of the method and evaluate its numerical performance by computing solutions to a wide range of test problems ranging from the nearly incompressible to the highly compressible flow regimes. This work was recently published in the Journal of Computational Physics (Gokhale et al., 2018). It is the first presentation of a dimensionally split cut cell method for the compressible Navier-Stokes equations in the literature. Finally, we also present an extension of the cut cell method to solve high Reynolds number turbulent automotive flows using a wall-modelled Large Eddy Simulation (WMLES) approach. A full description is provided of the coupling between the (implicit) LES solution and an equilibrium wall function on the cut cell mesh. The combined methodology is used to compute results for the turbulent flow over a square cylinder, and for flow over the SAE Notchback and DrivAer reference automotive geometries. We intend to publish the promising results as part of a future publication, which would be the first assessment of a WMLES Cartesian cut cell approach for computing automotive flows to be presented in the literature.

The fluid-structure interactions of wall-mounted slender structures, such as cilia, filaments, flaps, and flags, play an important role in a broad range of physical processes: from the coherent waving motion of vegetation, to the passive flow control capability of hair-like surface coatings. While these systems are ubiquitous, their coupled nonlinear response exhibits a wide variety of behaviours that is yet to be fully understood, especially when multiple structures are considered. The purpose of this work is to investigate, via numerical simulation, the fluid-structure interactions of arrays of slender structures over a range of input conditions. A direct modelling approach, whereby the individual structures and their dynamics are fully resolved, is realised via a lattice Boltzmann-immersed boundary model, which is coupled to two different structural solvers: an Euler-Bernoulli beam model, and a finite element model. Results are presented for three selected test cases - which build in scale from a single flap in a periodic array, to a small finite array of flaps, and finally to a large finite array - and the key behaviour modes are characterised and quantified. Results show a broad range of behaviours, which depend on the flow conditions and structural properties. In particular, the emergence of coherent waving motions are shown to be closely related to the natural frequency of the array. Furthermore, this behaviour is associated with a lock-in between the natural frequency of the array and the predicted frequency of the fluid instabilities. The original contributions of this work are: the development and application of a numerical tool for direct modelling of large arrays of slender structures; the characterisation of the behaviour of slender structures over a range of input conditions; and the exposition of key behaviour modes of slender structures and their relation to input conditions.

Aerothermodynamics is an integral component in the design and implementation of hypersonic transport systems. Accurate estimates of the aerodynamic forces and heat transfer rates are critical in trajectory analysis and for payload weight considerations. The present work seeks to investigate the ability of an unstructured Cartesian grid framework in modeling hypersonic viscous flows. The effectiveness of modeling viscous phenomena in hypersonic flows using the immersed boundary ghost cell methodology of this solver is analyzed. The capacity of this framework to predict the surface physics in a hypersonic non-reacting environment is investigated. High velocity argon gas flows past a 2-D cylinder are simulated for a set of freestream conditions (Reynolds numbers), and impact of the grid cell sizes on the quality of the solution is evaluated. Additionally, the formulation is verified over a series of hypersonic Mach numbers for the flow past a hemisphere, and compared to experimental results and empirical estimates. Next, a test case that involves flow separation and the interaction between a hypersonic shock wave and a boundary layer, and a separation bubble is investigated using various adaptive mesh refinement strategies. The immersed boundary ghost cell approach is tested with two temperature clipping strategies, and their impact on the overall solution accuracy and smoothness of the surface property predictions are compared. Finally, species diffusion terms in the conservation equations, and collision cross-section based transport coefficients are installed, and hypersonic flows in thermochemical nonequilibrium environments are studied, and comparisons of the off-surface flow properties and the surface physics predictions are evaluated. First, a 2-D cylinder in a hypersonic reacting air flow is tested with an adiabatic wall boundary condition. Next, the same geometry is tested to evaluate the viscous chemistry prediction capability of the solver with an isothermal wall boundary condition, and to identify the strengths and weaknesses of the immersed boundary ghost cell methodology in computing convective heating rates in such an environment.

Plunger pumps are needed for heavy duty sludge pumping at wastewater treatment facilities. America's leading pump manufacturer Wastecorp Inc. brought their plunger pump problem to us in late 2009. It was found that when the ow rate reaches a critical value, the plunger pump starts to generate a clicking noise. A one-dimensional model was built for studying the ow of a typical plunger pump operation. The velocities and pressures are calculated at certain interesting locations. Pressure jumps have been found while opening or closing the valves. The valve motion is then modeled with considerations to its geometry. The results show that as the plunger speed reaches a critical value, the valve moves more rapidly and more likely to hit the wall and generates a noise. We also provide a methodology to study the ow across the valve in higher resolution. A nite-di erence approach to the Navier-Stokes equations are presented with the immersed boundary method. / UOIT

Two immersed boundary methods (IBM) for the simulation of conjugate heat transfer problems with complex geometries are introduced: a finite element (IFEM) and a finite volume (IFVM) immersed boundary methods are discussed. In the IFEM a projection approach is presented for the coupled system of time-dependent incompressible Navier-Stokes equations (NSEs) and energy equation in conjunction with the immersed boundary method for solving fluid flow and heat transfer problems in the presence of rigid objects not represented by the underlying mesh. The IBM allows solving the flow for geometries with complex objects without the need of generating a body-fitted mesh. Dirichlet boundary constraints are satisfied applying a boundary force at the immersed body surface. Using projection and interpolation operators from the fluid volume mesh to the solid surface mesh (i.e., the “immersed” boundary) and vice versa, it is possible to impose the extra constraint to the NSEs as a Lagrange multiplier in a fashion very similar to the effect pressure has on the momentum equations to satisfy the divergence-free constraint. The IFEM approach presented shows third order accuracy in space and second order accuracy in time when the simulation results for the Taylor-Green decaying vortex are compared to the analytical solution. For the IFVM a ghost-cell approach with sharp interface scheme is used to enforce the boundary condition at the fluid/solid interface. The interpolation procedure at the immersed boundary preserves the overall second order accuracy of the base solver. The developed ghost-cell method is applied on a staggered configuration with the Semi-Implicit Method for Pressure-Linked Equations Revised algorithm. Second order accuracy in space and first order accuracy in time are obtained when the Taylor-Green decaying vortex test case is compared to the IFVM analytical solution. Computations were performed using the IFEM and IFVM approaches for the two-dimensional flow over a backward-facing step, two-dimensional flow past a stationary circular cylinder, three-dimensional flow past a sphere and two and three-dimensional natural convection in an enclosure with/without immersed body. The numerical results obtained with the discussed IFEM and IFVM were compared against other IBMs available in literature and simulations performed with the commercial computational fluid dynamics code STAR-CCM+/V7.04.006. The benchmark test cases showed that the numerical results obtained with the implemented immersed boundary methods are in good agreement with the predictions from STAR-CCM+ and the numerical data from the other IBMs. The immersed boundary method based of finite element approach is numerically more accurate than the IBM based on finite volume discretization. In contrast, the latter is computationally more efficient than the former.

In this dissertation, we explore direct-forcing immersed boundary methods (IBM) under the framework of the lattice Boltzmann method (LBM), which is called the direct-forcing immersed boundary-lattice Boltzmann method (IB-LBM). First, we derive the direct-forcing formula based on the split-forcing lattice Boltzmann equation, which recovers the Navier-Stokes equation with second-order accuracy and enables us to develop a simple and accurate formula due to its kinetic nature. Then, we assess the various interface schemes under the derived direct-forcing formula. We consider not only diffuse interface schemes but also a sharp interface scheme. All tested schemes show a second-order overall accuracy. In the simulation of stationary complex boundary flows, we can observe that the sharper the interface scheme is, the more accurate the results are. The interface schemes are also applied to moving boundary problems. The sharp interface scheme shows better accuracy than the diffuse interface schemes but generates spurious oscillation in the boundary forcing terms due to the discontinuous change of nodes for the interpolation. In contrast, the diffuse interface schemes show smooth change in the boundary forcing terms but less accurate results because of discrete delta functions. Hence, the diffuse interface scheme with a corrected radius can be adopted to obtain both accurate and smooth results. Finally, a direct-forcing immersed boundary method (IBM) for the thermal lattice Boltzmann method (TLBM) is proposed to simulate non-isothermal flows. The direct-forcing IBM formulas for thermal equations are derived based on two TLBM models: a double-population model with a simplified thermal lattice Boltzmann equation (Model 1) and a hybrid model with an advection-diffusion equation of temperature (Model 2). The proposed methods are validated through natural convection problems with stationary and moving boundaries. In terms of accuracy, the results obtained from the IBMs based on both models are comparable and show a good agreement with those from other numerical methods. In contrast, the IBM based on Model 2 is more numerically efficient than the IBM based on Model 1. Overall, this study serves to establish the feasibility of the direct-forcing IB-LBM as a viable tool for computing various complex and/or moving boundary flow problems.

Study of turbulent flow over a rough bed is highly important due to its numerous applications in the areas of sediment transport and pollutant discharge in streams, rivers and channels. Over the past few decades, many experimental studies have been conducted in this respect to understand the underlying phenomenon. However, there is a scarcity in the number of computational studies conducted on this topic. Therefore, a Large Eddy Simulation (LES) of turbulent flow over a rough channel bed was conducted to contribute further understanding of the influence of bed roughness on turbulent flow properties. For this purpose, an efficient, second order accurate 'immersed boundary method' was implemented into the LES code Hydro3d-GT, and validated for flow past bluff bodies. LES results from the present study showed excellent agreement with previous experimental studies on flow over rough beds. An in-depth analysis of time varying turbulent quantities (like the velocity fluctuations) revealed the presence of coherent structures in the flow. Also, a three dimensional visualization of the turbulent structures provided a good picture of the flow, especially in the near bed region, which is quite difficult to accomplish using experimental studies.