The thesis develops singly-periodic mathematical models for complex laminar wakes which are formed behind vortex-shedding bluff bodies. These wake structures exhibit a variety of patterns as the bodies oscillate or are in close proximity of one another. The most well-known formation comprises two counter-rotating vortices in each shedding cycle and is popularly known as the vk vortex street. Of the more complex configurations, as a specific example, this thesis investigates one of the most commonly occurring wake arrangements, which consists of two pairs of vortices in each shedding period. The paired vortices are, in general, counter-rotating and belong to a more general definition of the 2P mode, which involves periodic release of four vortices into the flow. The 2P arrangement can, primarily, be sub-classed into two types: one with a symmetric orientation of the two vortex pairs about the streamwise direction in a periodic domain and the other in which the two vortex pairs per period are placed in a staggered geometry about the wake centerline. The thesis explores the governing dynamics of such wakes and characterizes the corresponding relative vortex motion. In general, for both the symmetric as well as the staggered four vortex periodic arrangements, the thesis develops two-dimensional potential flow models (consisting of an integrable Hamiltonian system of point vortices) that consider spatially periodic arrays of four vortices with their strengths being +/-1 and +/-2. Vortex formations observed in the experiments inspire the assumed spatial symmetry. The models demonstrate a number of dynamic modes that are classified using a bifurcation analysis of the phase space topology, consisting of level curves of the Hamiltonian. Despite the vortex strengths in each pair being unequal in magnitude, some initial conditions lead to relative equilibrium when the vortex configuration moves with invariant size and shape. The scaled comparisons of the model results with experiments conducted in a flowing soap film with an airfoil, which was imparted with forced oscillations, are satisfactory and validate the reduced order modeling framework. The experiments have been performed by a collaborator group at the Department of Physics and Fluid Dynamics at the Technical University of Denmark (DTU), led by Dr. Anders Andersen. Similar experiments have also been run at Virginia Tech as part of this dissertation and the preliminary results are included in this treatise. The thesis also employs the same dynamical systems techniques, which have been applied to study the 2P regime dynamics, to develop a mathematical model for the P+S mode vortex wakes, with three vortices present in each shedding cycle. The model results have also been compared favorably with an experiment and the predictions regarding the vortex circulation data match well with the previous results from literature. Finally, the thesis introduces a novel concept of clean and renewable energy extraction from vortex-induced vibrations of bluff bodies. The slow-moving currents in the off-shore marine environments and riverine flows are beyond the operational capabilities of the more established hydrokinetic energy converters and the discussed technology promises to be a significant tool to generate useful power from these copiously available but previously untapped sources. / Ph. D.

Numerical simulations of fluid-structure interaction (FSI) problems were performed in an effort to verify and validate a commercially available FSI tool. This tool uses an iterative partitioned coupling scheme between CD-adapco's STAR-CCM+ finite volume fluid solver and Simulia's Abaqus finite element structural solver to simulate the FSI response of a system. Preliminary verification and validation work (VandV) was carried out to understand the numerical behavior of the codes individually and together as a FSI tool. Verification and Validation work that was completed included code order verification of the respective fluid and structural solvers with Couette-Pouiselle flow and Euler-Bernoulli beam theory. These results confirmed the 2nd order accuracy of the spatial discretizations used. Following that, a mixture of solution verifications and model calibrations was performed with the inclusion of the physics models implemented in the solution of the FSI problems. Solution verifications were completed for fluid and structural stand-alone models as well as for the coupled FSI solutions. These results re-confirmed the spatial order of accuracy but for more complex flows and physics models as well as the order of accuracy of the temporal discretizations. In lieu of a good material definition, model calibration is performed to reproduce the experimental results. This work used model calibration for both instances of hyperelastic materials which were presented in the literature as validation cases because these materials were defined as linear elastic. Calibrated, three dimensional models of the bow seal on the University of Michigan bow seal test platform showed the ability to reproduce the experimental results qualitatively through averaging of the forces and seal displacements. These simulations represent the only current 3D results for this case. One significant result of this study is the ability to visualize the flow around the seal and to directly measure the seal resistances at varying cushion pressures, seal immersions, forward speeds, and different seal materials. SES design analysis could greatly benefit from the inclusion of flexible seals in simulations, and this work is a positive step in that direction. In future work, the inclusion of more complex seal geometries and contact will further enhance the capability of this tool. / Ph. D.

This thesis investigates the characteristics of a close proximity underwater explosion and its effect on a ship-like structure. Finite element model tests are conducted to verify and validate the propagation of a pressure wave generated by an underwater explosion through a fluid medium, and the transmission of the pressure wave in the fluid to a structure using the Multi-Material Arbitrary Lagrangian/Eulerian method. A one dimensional case modeling the detonation of a spherical TNT charge underwater is investigated. Three dimensional cases modeling the detonation of an underwater spherical TNT charge, and US Navy Blast Test cases modeling a shape charge and a circular steel plate, and a shape charge and a Sandwich Plate System (SPS) are also investigated. This thesis provides evidence that existing tools and methodologies have some capability for predicting early-time/close proximity underwater explosion effects, but are insufficient for analyses beyond the arrival of the initial shock wave. This thesis shows that a true infinite boundary condition, a modified Gruneisen equation of state near the charge, and the ability to capture shock without a very small element size is needed in order to provide a sufficient means for predicting early-time/close proximity underwater explosion effects beyond the arrival of the initial shock wave. / Master of Science

<p dir="ltr">The focus is on the complex nature of primary blast injury (PBI) and employs advanced simulation techniques to model the physiological impacts using a TBI-on-a-chip system. This study involves a two-way Fluid-Structure Interaction (FSI) model in ANSYS, coupling Transient Structural and Fluent modules to simulate the effects of a blast wave on brain tissue. The research explores the creation and validation of boundary conditions, such as fixed support and varying strain rates, to ensure the reliability of the experimental setup. Key findings include the non-uniform distribution of strain, which has significant implications for understanding injury mechanisms and inflammatory marker analysis. The project also provides a detailed workflow for FSI simulations, highlighting the advantages of uniform mechanical loading and its impact on experimental accuracy.</p>

Cardiovascular diseases, including heart attacks and strokes, are leading causes of mortality worldwide. A significant contributor to these conditions is the formation of blood clots obstructing blood flow to critical organs. This study looks specifically at the blood flow in the left ventricle (LV) using simulations with Finite Element Method (FEM). Computational simulations of blood flow offer valuable insights into potential risks, aiding in pre-surgical assessments. However, interpreting simulation results, often presented as ambiguous velocity fields, can be difficult. This study examines two different post-processing measurements, triple decomposition of the velocity gradient tensor and E-wave propagation index (EPI), that can turn these complexities into simpler scalar fields and values. This study confirms that both measurements are computable and most likely useful given a velocity vector field. Further studies, particularly by medical professionals, need to be done to analyse the effects of the measurements and verify its usefulness by comparison to the same measurements in real, not simulated, patients. Some specific studies are proposed to address this. With the triple decomposition implemented some basic analysis was done. High levels of shear were observed in certain areas, particularly around the papillary muscles, where strain was also present.

This thesis gives an overview of the process that led to the development of a novel semi-implicit fluid-structure interaction model. The thesis is dedicated to the creation of a new numerical model that allows to study ship generated waves and ship manoeuvers in waterways for various vessel characteristics and speeds in different external current situations. A model like this requires a coupling between the fluid and the solid to generate the waves and the hydrodynamic forces on the hull. Since the horizontal dimensions are significantly larger than the vertical dimension, we started by employing the shallow water equations, which are based on the assumption of hydrostatic pressure. The discretization was carried out taking only the nonlinear advective terms explicitly while the pressure terms are discretized implicitly, which makes the CFL condition milder. The price to pay for this semi-implicit discretization is an increase in the algorithm complexity compared to a fully-explicit method, but it is still much simpler than a fully-implicit discretization of the governing equations. Indeed, the mass and momentum equations couple, and finding the unknowns involves solving a system of equations with dimensions equal to the number of cells. The grid supporting the discretization is staggered, overlapping and Cartesian. Since the aimed application domain is inland waterways, it is paramount to allow wetting and drying of the cells. This was achieved by acting on the depth function, the relationship between the free-surface elevation and the water depth in the cell. The main novelty of this research project is the two-way coupling of the PDE system for the water flow with the ODE system for the rigid body motion of the ship. The hull defines the ship region, and its shape can range from a simple box to an STL file of a real 3D ship geometry. Where the hull is in contact with the water, the cells are pressurized. This pressurized group of cells generates waves as it moves, and its motion is influenced by incoming external waves. This result is obtained by imposing an upper bound to the depth function, so that the water depth does not increase when it reaches the hull elevation, while the pressure is allowed to increase. This upper bound increases the nonlinearity of the system, which may have dry cells, wet free-surface cells and pressurized cells. The solution of this system is found by a single nested-Newton iterative solver of Casulli and Zanolli [36], in which with two separate linearizations the system is written in a sparse, symmetric, positive semi-definite form. This particular form allows us to employ a matrix-free conjugate gradient method, and efficiently get the unknown pressure. The integral of the pressure over the hull is applied for the hydrodynamic force and torque acting on the ship. After adding the skin friction and other external forces from the propeller or the rudder, the total force is inserted in the equation of motion of the rigid body. The ODE system is discretized with a second-order Taylor method, and it is solved for the six degrees of freedom (3 coordinates for the position vector of the barycenter and 3 rotation angles), providing the next position and orientation of the ship. The vertical translation of the rigid body is governed by the gravitational force and the restoring force from Archimedes' principle. As the ship oscillates up and down, the gravitational potential energy is partially transferred to the radiated free-surface water waves, damping and eventually stopping the motion. Also, the ship pushes and pulls the water around it, inducing the added mass force. All these elements constitute the ODE that was used for the verification of the vertical degree of freedom. The numerical simulation gave the expected results for the vertical motion. The horizontal translation, important for the manoeuvers, presented a numerical instability unseen in our previous test cases, which is connected to the relative motion between the ship and the grid. In each time step in which the ship enters a new cell, the pressure sharply increases and decreases at the ship bow. An oscillation can build up in time and create an unphysical void below the vessel. We implemented a few ideas to attenuate the oscillations. At the heart of all the following techniques is the reduction of the time derivative of the water depth, especially for those cells transitioning to a pressurized state. All these modifications were effective at controlling the oscillations, each with a different intensity, and simulations with a horizontal motion are much more stable than without these techniques. With the collaboration of the BAW research institute, we worked on the model validation. We used data from two separate experiments to compare the measurements with the numerical results. Specifically, we focused on the ship-generated wave height and the hydrodynamic forces on the hull. The comparison is satisfactory for the wave height. The force and torque prediction is plausible but underestimated compared to the measurements. The model seems to displace the water volume correctly during the ship passage, while the force and torque response might need additional work to be trusted in applications. Even though the hydrostatic assumption is mostly correct in our range of applications, the presence and the motion of a ship could generate strong vertical accelerations of the flow, which may not be negligible. For this reason, we implemented an algorithm that corrects the velocity field, introducing also dispersive effects due to a non-hydrostatic pressure. The correction consists of a higher-order Bousinnesq-type term in the momentum equation and the solution of the resulting system. The non-hydrostatic update has a small influence on the wave generation, while it alters significantly the reaction forces. The subgrid method implementation allowed to benefit from high-resolution bottom descriptions while keeping the grid size coarse. The same subgrid can also be used for a refined definition of the hull, which makes the volume computations more accurate. Furthermore, the subgrid introduces new possible states for the cells, as they can be partially dry or partially pressurized. These intermediate states translate into smoother transitions from one state to the other when the free-surface is close to the bathymetry or to the hull. Concerning the software implementation of the developed scheme, in order to improve the execution performance of the prototype script formulated initially in Matlab, the numerical method was rewritten as a Fortran program. Also, thanks to the domain decomposition technique and the MPI standard, each simulation can run in parallel on multiple CPUs, leveraging the computational power of supercomputers. The coupling of the PDE and ODE system, together with an appropriate redefinition of the depth function, proved to be a valuable method for studying fluid-structure interaction problems. The combination of efficient numerical techniques led to the development of a tool with a potential to be applied in the practice for the simulation of floating objects in wide domains.

Air cushion vehicle flexible seals operate in a complex and chaotic environment dominated by fluid-structure interaction. An efficient means to explore interdependencies between various governing parameters that affect performance is through high fidelity numerical simulation. As previous numerical efforts have employed separate iterative partitioned solvers, or have implemented simplified physics, the approaches have been complex, computationally expensive, or of limited utility. This research effort performs numerical simulations to verify and validate the commercial multi-physics tool STAR-CCM+ as a stand-alone partitioned approach for fluid-structure interaction problems with or without a free surface. A dimensional analysis is first conducted to identify potential non-dimensional forms of parameters related to seal resistance. Then, an implicit, Reynolds-averaged Navier-Stokes finite volume fluid solver is coupled to an implicit, nonlinear finite element structural solver to successfully replicate benchmark results for an elastic beam in unsteady laminar flow. To validate the implementation as a seal parameter exploratory tool, a planer bow seal model is developed and results are obtained for various cushion pressures and inflow speeds. Previous numerical and experimental results for deflection and resistance are compared, showing good agreement. An uncertainty analysis for inflow velocity reveals an inversely proportional resistance dependency. Using Abaqus/Explicit, methodologies are also developed for a two-way, loosely coupled explicit approach to large deformation fluid-structure interaction problems, with and without a free surface. Following numerous verification and validation problems, Abaqus is ultimately abandoned due to the inability to converge the fluid pressure field and achieve steady state. This work is a stepping stone for future researchers having interests in ACV seal design and other large deformation, fluid-structure interaction problems. By modeling all necessary physics within a verified and validated stand-alone approach, a designer's ability to comprehensively investigate seal geometries and interactions has never been more promising. / Ph. D. / Air cushion vehicles are specialized marine craft that utilize flexible seals to enable improved performance and fully amphibious operation. An efficient means to explore interdependencies between various seal design parameters that affect performance is through computer modeling of the fluid-structure interaction between the seal and the sea. This research effort performs numerical simulations to verify and validate the commercial multi-physics tool STAR-CCM+ as a single computer program for fluid-structure interaction problems occurring on the water surface. A dimensional analysis is first conducted to identify parameters related to seal resistance. Then, a fluid model is coupled to a structural model to successfully replicate benchmark results for a flexible beam in an oscillating fluid flow. To validate the implementation as a seal parameter exploratory tool, a model of an ACV bow seal is developed and results are obtained for various operational conditions and inflow speeds. Previous numerical and experimental results for seal deflection and seal resistance are compared, showing good agreement. This work is a stepping stone for future researchers having interests in ACV seal design and other large deformation, fluid-structure interaction problems. By modeling all necessary physics within a verified and validated stand-alone computer program, a designer’s ability to comprehensively investigate seal geometries and interactions has never been more promising.

In the context of interior noise reduction, the present work aims at proposing Finite Element (FE) solution strategies for interior structural-acoustic applications including 3D modelling of homogeneous and isotropic poroelastic materials, under timeharmonic excitations, and in the low frequency range. A model based on the Biot-Allard theory is used for the poroelastic materials, which is known to be very costly in terms of computational resources. Reduced models offer the possibility to enhance the resolution of such complex problems. However, their applicability to porous materials remained to be demonstrated.First, this thesis presents FE resolutions of poro-elasto-acoustic coupled problems using modal-based approaches both for the acoustic and porous domains. The original modal approach proposed for porous media, together with a dedicated mode selection and truncation procedure, are validated on 1D to 3D applications.In a second part, modal-reduced models are combined with a Padé approximants reconstruction scheme in order to further improve the efficiency.A concluding chapter presents a comparison and a combination of the proposed methods on a 3D academic application, showing promising performances. Conclusions are then drawn to provide indications for future research and tests to be conducted in order to further enhance the methodologies proposed in this thesis. / Dans le contexte de lutte contre les nuisances sonores, cette thèse porte sur le développement de méthodes de résolution efficaces par éléments finis, pour des problèmes de vibroacoustique interne avec interfaces dissipatives, dans le domaine des basses fréquences. L’étude se limite à l’utilisation de solutions passives telles que l’intégration de matériaux poreux homogènes et isotropes, modélisés par une approche fondée sur la théorie de Biot-Allard. Ces modèles étant coûteux en terme de résolution, un des objectifs de cette thèse est de proposer une approche modale pour la réduction du problème poroélastique, bien que l’adéquation d’une telle approche avec le comportement dynamique des matériaux poreux soit à démontrer.Dans un premier temps, la résolution de problèmes couplés élasto-poro-acoustiques par sous-structuration dynamique des domaines acoustiques et poreux est établie. L’approche modale originale proposée pour les milieux poroélastiques, ainsi qu’une procédure de sélection des modes significatifs, sont validées sur des exemples 1D à 3D.Une deuxième partie présente une méthode combinant l’utilisation des modèles réduits précédemment établis avec une procédure d’approximation de solution par approximants de Padé. Il est montré qu’une telle combinaison offre la possibilité d’accroître les performances de la résolution (allocation mémoire et ressources en temps de calcul).Un chapitre dédié aux applications permet d’évaluer et comparer les approches sur un problème académique 3D, mettant en valeur leurs performances encourageantes. Afin d’améliorer les méthodes établies dans cette thèse, des perspectives à ces travaux de recherche sont apportées en conclusion. / <p>QC 20120224</p> / FP6 Marie-Curie Smart Structures / FP7 Marie-Curie Mid-Frequency

Dans le contexte de lutte contre les nuisances sonores, cette thèse porte sur le développement de méthodes de résolution efficaces par éléments finis, pour des problèmes de vibroacoustique interne avec interfaces dissipatives, dans le domaine des basses fréquences. L’étude se limite à l’utilisation de solutions passives telles que l’intégration de matériaux poreux homogènes et isotropes, modélisés par une approche fondée sur la théorie de Biot-Allard. Ces modèles étant coûteux en terme de résolution, un des objectifs de cette thèse est de proposer une approche modale pour la réduction du problème poroélastique, bien que l’adéquation d’une telle approche avec le comportement dynamique des matériaux poreux soit à démontrer. Dans un premier temps, la résolution de problèmes couplés élasto-poro-acoustiques par sous-structuration dynamique des domaines acoustiques et poreux est établie. L’approche modale originale proposée pour les milieux poroélastiques, ainsi qu’une procédure de sélection des modes significatifs, sont validées sur des exemples 1D à 3D. Une deuxième partie présente une méthode combinant l’utilisation des modèles réduits précédemment établis avec une procédure d’approximation de solution par approximants de Padé. Il est montré qu’une telle combinaison offre la possibilité d’accroître les performances de la résolution (allocation mémoire et ressources en temps de calcul). Un chapitre dédié aux applications permet d’évaluer et comparer les approches sur un problème académique 3D, mettant en valeur leurs performances encourageantes. Afin d’améliorer les méthodes établies dans cette thèse, des perspectives à ces travaux de recherche sont apportées en conclusion. / In the context of interior noise reduction, the present work aims at proposing Finite Element (FE) solution strategies for interior structural-acoustic applications including 3D modelling of homogeneous and isotropic poroelastic materials, under timeharmonic excitations, and in the low frequency range. A model based on the Biot-Allard theory is used for the poroelastic materials, which is known to be very costly in terms of computational resources. Reduced models offer the possibility to enhance the resolution of such complex problems. However, their applicability to porous materials remained to be demonstrated.First, this thesis presents FE resolutions of poro-elasto-acoustic coupled problems using modal-based approaches both for the acoustic and porous domains. The original modal approach proposed for porous media, together with a dedicated mode selection and truncation procedure, are validated on 1D to 3D applications.In a second part, modal-reduced models are combined with a Padé approximants reconstruction scheme in order to further improve the efficiency.A concluding chapter presents a comparison and a combination of the proposed methods on a 3D academic application, showing promising performances. Conclusions are then drawn to provide indications for future research and tests to be conducted in order to further enhance the methodologies proposed in this thesis.

In this doctoral work, we adress various problems arising when dealing with multi-physical simulations using a segregated (non-monolithic) approach. We concentrate on a few specific problems and focus on the solution of aeroelastic <p>flutter for linear elastic structures in compressible fl<p>ows, conjugate heat transfer for re-entry vehicles including thermo-chemical reactions and finally, industrial electro-chemical plating processes which often include<p>stiff source terms. These problems are often solved using specifically developed<p>solvers, but these cannot easily be reused for different purposes. We have therefore considered the development of a <p>flexible and reusable software platform for the simulation of multi-physics problems. We have based this<p>development on the COOLFluiD framework developed at the von Karman Institute in collaboration with a group of partner institutions.<p>For the solution of fl<p>uid fl<p>ow problems involving compressible <p>flows, we have used the Finite Volume method and we have focused on the application of the method to moving and deforming computational domains using the Arbitrary Lagrangian Eulerian formulation. Validation on a series of testcases (including turbulent flows) is shown. In parallel, novel time integration<p>methods have been derived from two popular time discretization methods.<p>They allow to reduce the computational effort needed for unsteady fl<p>ow computations.<p>Good numerical properties have been obtained for both methods.<p>For the computations on deforming domains, a series of mesh deformation techniques are described and compared. In particular, the effect of the stiffness definition is analyzed for the Solid material analogy technique. Using<p>the techniques developed, large movements can be obtained while preserving a good mesh quality. In order to account for very large movements for which mesh deformation techniques lead to badly behaved meshes, remeshing is also considered.<p>We also focus on the numerical discretization of a class of physical models that are often associated with <p>fluid fl<p>ows in coupled problems. For the elliptic problems considered here (elasticity, heat conduction and electrochemical<p>potential problems), the implementation of a Finite Element solver is presented. Standard techniques are described and applied for a variety of problems, both steady and unsteady.<p>Finally, we discuss the coupling of the <p>fluid flow solver with the finite element solver for a series of applications. We concentrate only on loosely and strongly coupled algorithms and the issues associated with their use and implementation. The treatment of non-conformal meshes at the interface between two coupled computational domains is discussed and the problem<p>of the conservation of global quantities is analyzed. The software development of a <p>flexible multi-physics framework is also detailed. Then, several coupling algorithms are described and assessed for testcases in aeroelasticity and conjugate heat transfer showing the integration of the <p>fluid and solid solvers within a multi-physics framework. A novel strongly coupled algorithm, based on a Jacobian-Free Newton-Krylov method is also presented and applied to stiff coupled electrochemical potential problems. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished