Corrugated pipes are widely used in offshore gas and oil fields for their flexibility while offering local rigidity. However, self-sustained pressure pulsations associated with the flow in corrugated pipes results in a noisy environment, high running costs and eventually structure fatigue failure upon long exposure. Recent literature has addressed either the flow over a single cavity or the global oscillations. The current research aims at understanding the flow over multiple cavities as a first step to correlate the rich single cavity literature and the actual corrugated pipe problem with the ultimate goal of predicting oscillations amplitude in corrugate pipes. The standing wave method (SWM), which is an efficient experimental tool, has been successfully adapted in the first phase of the project to measure the source of multiple cavity configurations. One, two and three-cavity configurations have been investigated by means of the SWM. The source non-linearly becomes more pronounced as the number of cavities increases. The cavity length (L) is still found to be the appropriate length scale to define the oscillation dimensionless frequency (the Strouhal number). The measured source data have been successfully employed in a semi-empirical model to predict the amplitude of the self-excited oscillations. Accurate model performance is achieved for the single, double and triple cavity configurations. Including the absorption losses at the cavity corners has been found to be crucial for the model prediction accuracy. The separation distance (Lp) effect on the generated source is investigated for two and three-cavity configurations using the SWM over a practical range of spacing ratios. At extremum spacing ratios of (Lp/L) 0.5 and 1.375, constructive hydrodynamic interference associated with strong sources has been observed. At high excitation levels the source consistently becomes weaker upon increasing the spacing ratio. The reported trends are consistent for both the double and triple cavity configurations. However, the destructive interference spacing ratio is found to depend on the number of cavities indicating a relatively more complicated interaction mechanism. The different interaction patterns have been analytically interpreted based on the synchronization of the hydrodynamic cycle of the cavity shear layer and the disturbance convection along the pipe spacing between the cavities. Moreover, the three-cavity constructive interference cases have been visualized using Particle Image Velocimetry (PIV). The source evaluated based on the PIV data and applying Howe’s analogy revealed each cavity share of the global source, which fairly agrees with the SWM measured source. The source contribution due to gradually increasing the number of cavities is investigated using the SWM up to a six-cavity configuration. The source contribution reaches asymptotically a consistent value starting from the fourth cavity. This persistent contribution defines a building unit cavity source which is representative of a general cavity in a long corrugated pipe. The building unit source fairly agrees with the ninth-cavity source in a twelve-cavity configuration extracted by means of the PIV technique. Finally, a predication model, based on the building unit source, successfully predicts the oscillations amplitude of a twelve-cavity configuration, which serves as a model for a corrugated pipe section. / Thesis / Doctor of Philosophy (PhD)

Recently, numerous flexible structures have been employed in various fields of industry. Loading conditions sustained by these flexible structures are often not described well enough for engineering analyses even though these conditions are important. Here, a flexible tent with an interior Collective Protection System, which is subjected to an explosion, is analyzed. The tent protects personnel from biological and chemical agents with a pressurized liner inside the tent as an environmental barrier. Field tests showed unexpected damage to the liner, and most of the damage occurred on tent's leeward side. To solve this problem, various tests and analyses have been performed, involving material characteristics of the liner, canvas, and zip seals, modeling of the blast loading over the tent and inside the tent, and structural response of the tent to the blast loading as collaborative research works with others. It was found that the blast loading and the structural response can not be analyzed separately due to the interaction between the flexible structure and the dynamic pressure loading. In this dissertation, the dynamic loadings imposed on both the interior and the exterior sides of the tent structure due to the airblasts and the resulting dynamic responses were studied. First, the blast loadings were obtained by a newly proposed theoretical method of analytical/empirical models which was developed into a FORTRAN program. Then, a numerical method of an iterative Fluid-Structure Interaction using Computational Fluid Dynamics and Computational Structural Dynamics was employed to simulate the blast wave propagation inside and outside the flexible structure and to calculate the dynamic loads on it. All the results were compared with the field test data conducted by the Air Force Research Laboratory. The experimental pressure data were gathered from pressure gauges attached to the tent surfaces at different locations. The comparison showed that the proposed methods can be a good design tool to analyze the loading conditions for rigid or flexible structures under explosive loads. In particular, the causes of the failure of the liner on the leeward were explained. Also, the results showed that the effect of fluid-structure interaction should be considered in the pressure load calculation on the structure where the structural deflection rate can influence the solution of the flow field surrounding the structure. / Ph. D.

Arch-shaped microelectromechanical systems (MEMS) have been used as mechanical memories, micro-sensors, micro-actuators, and micro-valves. A bi-stable structure, such as an arch, is characterized by a multivalued load deflection curve. Here we study the symmetry breaking, the snap-through instability, and the pull-in instability of bi-stable arch shaped MEMS under steady and transient electric loads. We analyze transient finite electroelastodynamic deformations of perfect electrically conducting clamped-clamped beams and arches suspended over a flat rigid semi-infinite perfect conductor. The coupled nonlinear partial differential equations (PDEs) for mechanical deformations are solved numerically by the finite element method (FEM) and those for the electrical problem by the boundary element method. The coupled nonlinear PDE governing transient deformations of the arch based on the Euler-Bernoulli beam theory is solved numerically using the Galerkin method, mode shapes for a beam as basis functions, and integrated numerically with respect to time. For the static problem, the displacement control and the pseudo-arc length continuation (PALC) methods are used to obtain the bifurcation curve of arch's deflection versus the electric potential. The displacement control method fails to compute arch's asymmetric deformations that are found by the PALC method. For the dynamic problem, two distinct mechanisms of the snap-through instability are found. It is shown that critical loads and geometric parameters for instabilities of an arch with and without the consideration of mechanical inertia effects are quite different. A phase diagram between a critical load parameter and the arch height is constructed to delineate different regions of instabilities. The local water slamming refers to the impact of a part of a ship hull on stationary water for a short duration during which high local pressures occur. We simulate slamming impact of rigid and deformable hull bottom panels by using the coupled Lagrangian and Eulerian formulation in the commercial FE software LS-DYNA. The Lagrangian formulation is used to describe planestrain deformations of the wedge and the Eulerian description of motion for deformations of the water. A penalty contact algorithm couples the wedge with the water surface. Damage and delamination induced, respectively, in a fiber reinforced composite panel and a sandwich composite panel and due to hydroelastic pressure are studied. / Ph. D.

The vulnerability of a ship to the impact of underwater explosions (UNDEX) and how to incorporate this factor into early-stage ship design is an important aspect in the ship survivability study. In this dissertation, attention is focused on the cost-efficient simulation of the ship response to a far-field UNDEX which involves fluid shock waves, cavitation, and fluid-structural interaction. Traditional fluid numerical simulation approaches using the Finite Element Method to track wave propagation and cavitation requires a high-level of mesh refinement to prevent numerical dispersion from discontinuities. Computation also becomes quite expensive for full ship-related problems due to the large fluid domain necessary to envelop the ship. The burden is aggravated by the need to generate a fluid mesh around the irregular ship hull geometry, which typically requires significant manual intervention. To accelerate the design process and enable the consideration of far-field UNDEX vulnerability, several contributions are made in this dissertation to make the simulation more efficient. First, a Cavitating Acoustic Spectral Element approach which has shown computational advantages in UNDEX problems, but not systematically assessed in total ship application, is used to model the fluid. The use of spectral elements shows greater structural response accuracy and lower computational cost than the traditional FEM. Second, a novel fully automatic all-hexahedral mesh generation scheme is applied to generate the fluid mesh. Along with the spectral element, the all-hex mesh shows greater accuracy than the all-tetrahedral finite element mesh which is typically used. This new meshing approach significantly saves time for mesh generation and allows the spectral element, which is confined to the hexahedral element, to be applied in practical ship problems. A further contribution of this dissertation is the development of a surrogate non-numerical approach to predict structural peak responses based on the shock factor concept. The regression analysis reveals a reasonably strong linear relationship between the structural peak response and the shock factor. The shock factor can be conveniently employed in the design aspects where the peak response is sufficient, using much less computational resources than numerical solvers. / Doctor of Philosophy / The vulnerability of a ship to the impact of underwater explosions (UNDEX) and how to incorporate this factor into early-stage ship design is an important aspect in the ship survivability study. In this dissertation, attention is focused on the cost-efficient simulation of the ship response to a far-field UNDEX which involves fluid shock waves, cavitation, and fluid-structural interaction. Traditional fluid numerical simulation approaches using the Finite Element Method to track wave propagation and cavitation requires a highly refined mesh to deal with large numerical errors. Computation also becomes quite expensive for full ship-related problems due to the large fluid domain necessary to envelop the ship. The burden is aggravated by the need to generate a fluid mesh around the irregular ship hull geometry, which typically requires significant manual intervention. To accelerate the design process and enable the consideration of far-field UNDEX vulnerability, several contributions are made in this dissertation to make the simulation more efficient. First, a Cavitating Acoustic Spectral Element approach, which has shown computational advantages in UNDEX problems but not systematically assessed in total ship application, is used to model the fluid. The use of spectral elements shows greater structural response accuracy and lower computational cost than the traditional FEM. Second, a novel fully automatic all-hexahedral mesh generation scheme is applied to generate the fluid mesh. Along with the spectral element, the all-hex mesh shows greater accuracy than the all-tetrahedral finite element mesh which is typically used. A further contribution of this dissertation is the development of a non-numerical approach which can approximate peak structural responses comparable to the numerical solution with far less computational effort.

This dissertation investigates the development of numerical algorithms for coupling computational fluid dynamics (CFD) and computational solid dynamics (CSD) solvers, and the use of these solvers for simulating fluid-solid interaction (FSI) problems involving large deformation, shock waves, and multiphase flow. The dissertation consists of two parts. The first part investigates the use of Robin interface conditions to resolve the well-known numerical added-mass instability, which affects partitioned coupling procedures for solving problems with incompressible flow and strong added-mass effect. First, a one-parameter Robin interface condition is developed by linearly combining the conventional Dirichlet and Neumann interface conditions. Next, a numerical algorithm is developed to implement the Robin interface condition in an embedded boundary method for coupling a parallel, projection-based incompressible viscous flow solver with a nonlinear finite element solid solver. Both an analytical study and a numerical study reveal that the new algorithm can clearly outperform conventional Dirichlet-Neumann procedures in terms of both stability and accuracy, when the parameter value is carefully selected. Moreover, the studies also indicate that the optimal parameter value depends on the materials and geometry of the problem. Therefore, to efficiently solve FSI problems involving non-uniform structures, a generalized Robin interface condition is presented, in which the constant parameter is replaced by a spatially varying function that depends on the local material and geometric properties of the structure. Numerical experiments using two benchmark problems show that the spatially varying Robin interface condition can clearly improve numerical accuracy compared to the constant- parameter version with the same computational cost. The second part of this dissertation focuses on simulating complex FSI problems featuring shock waves, multiphase flow (e.g., bubbles), and shock-induced material damage and fracture. A recently developed three-dimensional computational framework is employed, which couples a multiphase, compressible CFD solver and a nonlinear finite element CSD solver using an embedded boundary method and a partitioned procedure. In particular, the CFD solver applies a level-set method to capture the evolution of the bubble surface, and the CSD solver utilizes a continuum damage mechanics model and an element erosion method to simulate the dynamic fracture of the material. Two computational studies are presented. The first one investigates the dynamic response and failure of a brittle material exposed to a prescribed shock wave. The predictive capability of the computational framework is first demonstrated by simulating a series of laboratory experiments in the context of shock wave lithotripsy. Then, a parametric study is conducted to elucidate the significant effects of the shock wave's profile on material damage. In the second study, the computational framework is applied to simulate shock-induced bubble collapse near various solid and soft materials. The reciprocal effect of the material's properties (e.g., acoustic impedance, Young's modulus) on bubble dynamics is discussed in detail. / Doctor of Philosophy / Numerical simulations that couple computational fluid dynamics (CFD) solvers and computational solid dynamics (CSD) solvers have been widely used in the solution of nonlinear fluid-solid interaction (FSI) problems underlying many engineering applications. This is primarily because they are based on partitioned solutions of fluid and solid subsystems, which facilitates the use of existing numerical methods and computational codes developed for each subsystem. The first part of this dissertation focuses on developing advanced numerical algorithms for coupling the two subsystems. The aim is to resolve a major numerical instability issue that occurs when solving problems involving incompressible, heavy fluids and thin, lightweight structures. Specifically, this work first presents a new coupling algorithm based on a one-parameter Robin interface condition. An embedded boundary method is developed to enforce the Robin interface condition, which can be advantageous in solving problems involving complex geometry and large deformation. The new coupling algorithm has been shown to significantly improve numerical stability when the constant parameter is carefully selected. Next, the constant parameter is generalized into a spatially varying function whose local value is determined by the local material and geometric properties of the structure. Numerical studies show that when solving FSI problems involving non-uniform structures, using this spatially varying Robin interface condition can outperform the constant-parameter version in both stability and accuracy under the same computational cost. In the second part of this dissertation, a recently developed three-dimensional multiphase CFD - CSD coupled solver is extended to simulate complex FSI problems featuring shock wave, bubbles, and material damage and fracture. The aim is to understand the material’s response to loading induced by a shock wave and the collapse of nearby bubbles, which is important for advancing the beneficial use of shock wave and bubble collapse for material modification. Two computational studies are presented. The first one investigates the dynamic response and failure of a brittle material exposed to a prescribed shock wave. The causal relationship between shock loading and material failure, and the effects of the shock wave’s profile on material damage are discussed. The second study investigates the shock-induced bubble collapse near various solid and soft materials. The two-way interaction between bubble dynamics and materials response, and the reciprocal effects of the material’s properties are discussed in detail.

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.

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.

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

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.