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.

This dissertation develops a computational method for simulating laser-induced cavitation and investigates the mechanism behind the formation of non-spherical bubbles induced by long-pulsed lasers. The proposed computational method accounts for the laser emission and absorption, phase transition, and the dynamics and thermodynamics of a two-phase fluid flow. In this new method, the model combines the Navier-Stokes (NS) equations for a compressible inviscid two-phase fluid flow, a new laser radiation equation, and a novel local thermodynamic model of phase transition. The Navier-Stokes equations are solved using the FInite Volume method with Exact two-phase Riemann solvers (FIVER). Following this method, numerical fluxes across phase boundaries are computed by constructing and solving one-dimensional bi-material Riemann problems. The new laser radiation equation is derived by customizing the radiative transfer equation (RTE) using the special properties of laser, including monochromaticity, directionality, high intensity, and a measurable focusing or diverging angle. An embedded boundary finite volume method is developed to solve the laser radiation equation on the same mesh created for the NS equations. The fluid mesh usually does not resolve the boundary and propagation directions of the laser beam, leading to the challenges of imposing the boundary conditions on the laser domain. To overcome this challenge, ghost nodes outside the laser domain are populated by mirroring and interpolation techniques. The existence and uniqueness of the solution are proved for the two-dimensional case, leveraging the special geometry of the laser domain. The method is up to second-order accuracy, which is also proved, and verified using numerical tests. A method of latent heat reservoir is developed to predict the onset of vaporization, which accounts for the accumulation and release of latent heat. In this work, the localized level set method is employed to track the bubble surface. Furthermore, the continuation of phase transition is possible in laser-induced cavitation problems, especially for long-pulsed lasers. A method of local correction and reinitialization is developed to account for continuous phase transitions. Several numerical tests are presented to verify the convergence of these methods. This multiphase laser-fluid coupled computational model is employed to simulate the formation and expansion of bubbles with different shapes induced by different long-pulsed lasers. The simulation results show that the computational method can capture the key phenomena in the laser-induced cavitation problems, including non-spherical bubble expansion, shock waves, and the ``Moses effect''. Additionally, the observed complex non-spherical shapes of vapor bubbles generated by long-pulsed laser reflect some characteristics (e.g., direction, width) of the laser beam. The dissertation also investigates the relation between bubble shapes and laser parameters and explores the transition between two commonly observed shapes -- namely, a rounded pear-like shape and an elongated conical shape -- using the proposed computational model. Two laboratory experiments are simulated, in which Holmium:YAG and Thulium fiber lasers are used respectively to generate bubbles of different shapes. In both cases, the predicted bubble nucleation and morphology agree reasonably well with the experimental observation. The full-field results of laser radiance, temperature, velocity, and pressure are analyzed to explain bubble dynamics and energy transmission. It is found that due to the lasting energy input, the vapor bubble's dynamics is driven not only by advection, but also by the continued vaporization at its surface. Vaporization lasts less than 1 microsecond in the case of the pear-shaped bubble, compared to over 50 microseconds for the elongated bubble. It is thus hypothesized that the bubble's morphology is determined by a competition between the speed of bubble growth due to advection and continuous vaporization. When the speed of advection is higher than that of vaporization, the bubble tends to grow spherically. Otherwise, it elongates along the laser beam direction. To test this hypothesis, the two speeds are defined analytically using a model problem and then estimated for the experiments using simulation results. The results support the hypothesis and also suggest that when the laser's power is fixed, a higher laser absorption coefficient and a narrower beam facilitate bubble elongation. / Doctor of Philosophy / Laser-induced cavitation is a process where laser beams create bubbles in a liquid. This phenomenon is widely applied in research and microfluidic applications for precise control of bubble dynamics. It also naturally occurs in various laser-based processes involving liquid environments. Understanding laser-induced cavitation is important for enhancing the effectiveness and safety of related technologies. However, experimental studies encounter limitations, highlighting the development of numerical methods to advance the understanding of laser-induced cavitation. The laser-induced cavitation can be roughly described as localized boiling through thermal radiation. The detailed physics involves the absorption of laser light by a liquid, the formation of vapor bubbles due to localized heating, and the dynamics of both the bubbles and the surrounding liquid. The first part of the dissertation introduces a new computational method for modeling these phenomena. The dynamics of the two-phase flow are modeled by the Navier-Stokes equations, which are solved using the FInite Volume method with Exact two-phase Riemann solvers (FIVER). The absorption of the laser light is modeled by a new laser radiation equation, which is derived from laser energy conservation and special properties of the laser. An embedded boundary finite volume method is developed to solve this equation on the same mesh created for the NS equations. Additionally, a method of latent heat reservoir is developed to predict the onset of vaporization. In this work, the level set method is employed to track the bubble surface, and a method of local correction and reinitialization is developed to account for possible continuous phase transitions. After developing this new method, several test cases are simulated. The simulation results show that the method can capture the key phenomena in the laser-induced cavitation problems, including the absorption of laser light, non-spherical bubble expansion, and shock waves. When the laser pulse is comparable to or longer than the acoustic time scale (long-pulsed laser), vapor bubbles generated often have complex non-spherical shapes. The bubble shapes reflect some characteristics (e.g., direction, width) of the laser beam. The second part of the dissertation investigates the relation between bubble shapes and laser parameters. Two laboratory experiments are simulated, in which two different lasers are used to generate bubbles of different shapes, namely, a rounded pear-like shape and an elongated conical shape. In both cases, the simulated bubbles exhibit shapes and sizes that reasonably match the experimental results. The simulation results of temperature, pressure, and velocity fields are analyzed to explain bubble dynamics and energy transmission. The analysis shows that the expansion of bubbles induced by long-pulsed lasers is determined not only by advection but also by the continued vaporization at its surface. Vaporization lasts less than $1$ microsecond in the case of the pear-shaped bubble, compared to over $50$ microseconds for the elongated bubble. It is thus hypothesized that the bubble expansion is determined by a competition between the speed of bubble growth due to advection and continuous vaporization. When the speed of advection is higher than that of vaporization, the bubble tends to grow spherically. Otherwise, it elongates along the laser beam direction. To test this hypothesis, the two speeds are defined analytically using a model problem and then estimated for the experiments using simulation results. The results support the hypothesis and also suggest that when the laser's power is fixed, a higher laser absorption coefficient and a narrower beam facilitate bubble elongation.

La méthode de projection et l'approche variationnelle de Sasaki sont deux techniques permettant d'obtenir un champ vectoriel à divergence nulle à partir d'un champ initial quelconque. Pour une vitesse d'un vent en haute altitude, un champ de vitesse sur une grille décalée est généré au-dessus d'une topographie donnée par une fonction analytique. L'approche cartésienne nommée Embedded Boundary Method est utilisée pour résoudre une équation de Poisson découlant de la projection sur un domaine irrégulier avec des conditions aux limites mixtes. La solution obtenue permet de corriger le champ initial afin d'obtenir un champ respectant la loi de conservation de la masse et prenant également en compte les effets dûs à la géométrie du terrain. Le champ de vitesse ainsi généré permettra de propager un feu de forêt sur la topographie à l'aide de la méthode iso-niveaux. L'algorithme est décrit pour le cas en deux et trois dimensions et des tests de convergence sont effectués. / The Projection method and Sasaki's variational technique are two methods allowing one to extract a divergence-free vector field from any vector field. From a high altitude wind speed, a velocity field is generated on a staggered grid over a topography given by an analytical function. The Cartesian grid Embedded Boundary method is used for solving a Poisson equation, obtained from the projection, on an irregular domain with mixed boundary conditions. The solution of this equation gives the correction for the initial velocity field to make it such that it satisfies the conservation of mass and takes into account the effects of the terrain. The incompressible velocity field will be used to spread a wildfire over the topography with the Level Set Method. The algorithm is described for the two and three dimensional cases and convergence tests are done.

La méthode de projection et l'approche variationnelle de Sasaki sont deux techniques permettant d'obtenir un champ vectoriel à divergence nulle à partir d'un champ initial quelconque. Pour une vitesse d'un vent en haute altitude, un champ de vitesse sur une grille décalée est généré au-dessus d'une topographie donnée par une fonction analytique. L'approche cartésienne nommée Embedded Boundary Method est utilisée pour résoudre une équation de Poisson découlant de la projection sur un domaine irrégulier avec des conditions aux limites mixtes. La solution obtenue permet de corriger le champ initial afin d'obtenir un champ respectant la loi de conservation de la masse et prenant également en compte les effets dûs à la géométrie du terrain. Le champ de vitesse ainsi généré permettra de propager un feu de forêt sur la topographie à l'aide de la méthode iso-niveaux. L'algorithme est décrit pour le cas en deux et trois dimensions et des tests de convergence sont effectués. / The Projection method and Sasaki's variational technique are two methods allowing one to extract a divergence-free vector field from any vector field. From a high altitude wind speed, a velocity field is generated on a staggered grid over a topography given by an analytical function. The Cartesian grid Embedded Boundary method is used for solving a Poisson equation, obtained from the projection, on an irregular domain with mixed boundary conditions. The solution of this equation gives the correction for the initial velocity field to make it such that it satisfies the conservation of mass and takes into account the effects of the terrain. The incompressible velocity field will be used to spread a wildfire over the topography with the Level Set Method. The algorithm is described for the two and three dimensional cases and convergence tests are done.

Underwater explosions (UNDEX) and a ship's vulnerability to them are problems of interest in early-stage ship design. A series of events occur sequentially in an UNDEX scenario in both the fluid and structural domains and these events happen over a wide range of time and spatial scales. Because of the complexity of the physics involved, it is a common practice to separate the description of UNDEX into early-time and late-time, and far-field and near-field. The research described in this dissertation is focused on the simulation of near-field and early-time UNDEX. It assembles a hybrid framework of algorithms to provide results while maintaining computational efficiency. These algorithms include Runge-Kutta, Discontinuous Galerkin, Level Set, Direct Ghost Fluid and Embedded Boundary methods. Computational fluid dynamics (CFD) solvers are developed using this framework of algorithms to demonstrate the computational methods and their ability to effectively and efficiently solve UNDEX problems. Contributions, made in the process of satisfying the objective of this research include: the derivation of eigenvectors of flux Jacobians and their application to the implementation of the slope limiter in the fluid discretization; the three-dimensional extension of Direct Ghost Fluid Method and its application to the multi-fluid treatment in UNDEX flows; the enforcement of an improved non-reflecting boundary condition and its application to UNDEX simulations; and an improvement to the projection-based embedded boundary method and its application to fluid-structure interaction simulations of UNDEX problems. / Doctor of Philosophy / Underwater explosions (UNDEX) and a ship's vulnerability to them are problems of interest in early-stage ship design. A series of events occur sequentially in an UNDEX scenario in both the fluid and structural domains and these events happen over a wide range of time and spatial scales. Because of the complexity of the physics involved, it is a common practice to separate the description of UNDEX into early-time and late-time, and far-field and near-field. The research described in this dissertation is focused on the simulation of near-field and early-time UNDEX. It assembles a hybrid framework of algorithms to provide results while maintaining computational efficiency. These algorithms include Runge-Kutta, Discontinuous Galerkin, Level Set, Direct Ghost Fluid and Embedded Boundary methods. Computational fluid dynamics (CFD) solvers are developed using this framework of algorithms to demonstrate these computational methods and their ability to effectively and efficiently solve UNDEX problems.