Browsing by Author "Wang, Kevin Guanyuan"
Now showing 1 - 20 of 25
Results Per Page
Sort Options
- A Computational Framework for Long-Term Atomistic Analysis of Solute Diffusion in NanomaterialsSun, Xingsheng (Virginia Tech, 2018-10-04)Diffusive Molecular Dynamics (DMD) is a class of recently developed computational methods for the simulation of long-term mass transport with a full atomic fidelity. Its basic idea is to couple a discrete kinetic model for the evolution of mass transport process with a non-equilibrium thermodynamics model that governs lattice deformation and supplies the requisite driving forces for kinetics. Compared to previous atomistic models, e.g., accelerated Molecular Dynamics and on-the-fly kinetic Monte Carlo, DMD allows the use of larger time-step sizes and hence has a larger simulation time window for mass transport problems. This dissertation focuses on the development, assessment and application of a DMD computational framework for the long-term, three-dimensional, deformation-diffusion coupled analysis of solute mass transport in nanomaterials. First, a computational framework is presented, which consists mainly of: (1) a computational model for interstitial solute diffusion, which couples a nonlinear optimization problem with a first-order nonlinear ordinary differential equation; (2) two numerical methods, i.e., mean field approximation and subcycling time integration, for accelerating DMD simulations; and (3) a high-performance computational solver, which is parallelized based on Message Passing Interface (MPI) and the PETSc/TAO library for large-scale simulations. Next, the computational framework is validated and assessed in two groups of numerical experiments that simulate hydrogen mass transport in palladium. Specifically, the framework is validated against a classical lattice random walk model. Its capability to capture the atomic details in nanomaterials over a long diffusive time scale is also demonstrated. In these experiments, the effects of the proposed numerical methods on solution accuracy and computation time are assessed quantitatively. Finally, the computational framework is employed to investigate the long-term hydrogen absorption into palladium nanoparticles with different sizes and shapes. Several significant findings are shown, including the propagation of an atomistically sharp phase boundary, the dynamics of solute-induced lattice deformation and stacking faults, and the effect of lattice crystallinity on absorption rate. Specifically, the two-way interaction between phase boundary propagation and stacking fault dynamics is noteworthy. The effects of particle size and shape on both hydrogen absorption and lattice deformation are also discussed in detail.
- Computational Reconstruction and Quantification of Aerospace MaterialsLong, Matthew Thomas (Virginia Tech, 2024-05-14)Microstructure reconstruction is a necessary tool for use in multi-scale modeling, as it allows for the analysis of the microstructure of a material without the cost of measuring all of the required data for the analysis. For microstructure reconstruction to be effective, the synthetic microstructure needs to predict what a small sample of measured data would look like on a larger domain. The Markov Random Field (MRF) algorithm is a method of generating statistically similar microstructures for this process. In this work, two key factors of the MRF algorithm are analyzed. The first factor explored is how the base features of the microstructure related to orientation and grain/phase topology information influence the selection of the MRF parameters to perform the reconstruction. The second focus is on the analysis of the numerical uncertainty (epistemic uncertainty) that arises from the use of the MRF algorithm. This is done by first removing the material uncertainty (aleatoric uncertainty), which is the noise that is inherent in the original image representing the experimental data. The epistemic uncertainty that arises from the MRF algorithm is analyzed through the study of the percentage of isolated pixels and the difference in average grain sizes between the initial image and the reconstructed image. This research mainly focuses on two different microstructures, B4C-TiB2 and Ti-7Al, which are a ceramic composite and a metallic alloy, respectively. Both of them are candidate materials for many aerospace systems owing to their desirable mechanical performance under large thermo-mechanical stresses.
- Computationally-effective Modeling of Far-field Underwater Explosion for Early-stage Surface Ship DesignLu, Zhaokuan (Virginia Tech, 2020-03-23)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.
- Coupled Adjoint-based Sensitivity Analysis using a FSI Method in Time Spectral FormKim, Hyunsoon (Virginia Tech, 2019-09-26)A time spectral and coupled adjoint based sensitivity analysis of rotor blade is carried out in this study. The time spectral method is an efficient technique to solve unsteady periodic problems by transforming unsteady equation of motion to a steady state one. Due to the availability of the governing equations in the steady form, the steady form of the adjoint equations can be applied for the sensitivity analysis of the coupled fluid-structure system. An expensive computational time and memory requirement for the unsteady adjoint sensitivity analysis is thus avoided. A coupled analysis of fluid, structural, and flight dynamics is carried out through a CFD/CSD/CA coupling procedure that combines FSI analysis with enforced trim condition. Coupled sensitivity analysis results and their validations are presented and compared with aerodynamics only sensitivity analysis results. The fluid-structure coupled adjoint based sensitivity analysis will be applied to the shape optimization of a rotor blade in the future work. Minimization of required power is the objective of the optimization problem with constraints on thrust and drag of the rotor. The bump functions are considered as the design variables. Rotor blade shape changes are obtained by using the bump function on the surface of the airfoil sections along the span.
- CPU/GPU Code Acceleration on Heterogeneous Systems and Code Verification for CFD ApplicationsXue, Weicheng (Virginia Tech, 2021-01-25)Computational Fluid Dynamics (CFD) applications usually involve intensive computations, which can be accelerated through using open accelerators, especially GPUs due to their common use in the scientific computing community. In addition to code acceleration, it is important to ensure that the code and algorithm are implemented numerically correctly, which is called code verification. This dissertation focuses on accelerating research CFD codes on multi-CPUs/GPUs using MPI and OpenACC, as well as the code verification for turbulence model implementation using the method of manufactured solutions and code-to-code comparisons. First, a variety of performance optimizations both agnostic and specific to applications and platforms are developed in order to 1) improve the heterogeneous CPU/GPU compute utilization; 2) improve the memory bandwidth to the main memory; 3) reduce communication overhead between the CPU host and the GPU accelerator; and 4) reduce the tedious manual tuning work for GPU scheduling. Both finite difference and finite volume CFD codes and multiple platforms with different architectures are utilized to evaluate the performance optimizations used. A maximum speedup of over 70 is achieved on 16 V100 GPUs over 16 Xeon E5-2680v4 CPUs for multi-block test cases. In addition, systematic studies of code verification are performed for a second-order accurate finite volume research CFD code. Cross-term sinusoidal manufactured solutions are applied to verify the Spalart-Allmaras and k-omega SST model implementation, both in 2D and 3D. This dissertation shows that the spatial and temporal schemes are implemented numerically correctly.
- Cross-Sectional Stiffness Properties of Complex Drone WingsMuthirevula, Neeharika (Virginia Tech, 2017-01-05)The main purpose of this thesis is to develop a beam element in order to model the wing of a drone, made of composite materials. The proposed model consists of the framework for the structural design and analysis of long slender beam like structures, e.g., wings, wind turbine blades, and helicopter rotor blades, etc. The main feature consists of the addition of the coupling between axial and bending with torsional effects that may arise when using composite materials and the coupling stemming from the inhomogeneity in cross-sections of any arbitrary geometry. This type of modeling approach allows for an accurate yet computationally inexpensive representation of a general class of beam-like structures. The framework for beam analysis consists of main two parts, cross-sectional analysis of the beam sections and then using this section analysis to build up the finite element model. The cross-sectional analysis is performed in order to predict the structural properties for composite sections, which are used for the beam model. The thesis consists of the model to validate the convergence of the element size required for the cross-sectional analysis. This follows by the validation of the shell models of constant cross-section to assess the performance of the beam elements, including coupling terms. This framework also has the capability of calculating the strains and displacements at various points of the cross-section. Natural frequencies and mode shapes are compared for different cases of increasing complexity with those available in the papers. Then, the framework is used to analyze the wing of a drone and compare the results to a model developed in NASTRAN.
- Definition of Damage Volumes for the Rapid Prediction of Ship Vulnerability to AIREX Weapon EffectsStark, Sean Aaron (Virginia Tech, 2016-09-09)This thesis presents a damage model developed for the rapid prediction of the vulnerability of a ship concept design to AIREX weapon effects. The model uses simplified physics-based and empirical equations, threat charge size, geometry of the design, and the structure of the design as inputs. The damage volumes are customized to the design being assessed instead using of a single volume defined only by the threat charge size as in previous damage ellipsoid methods. This methodology is validated against a range of charge sizes and a library of notional threats is created. The model uses a randomized hit distribution that is generated using notional threat targeting and the geometry of the design. A Preliminary Arrangement and Vulnerability (PAandV) model is updated with this methodology and used to calculate an Overall Measure of Vulnerability (OMOV) by determining equipment failures and calculating the resulting loss of mission capabilities. A selection of baseline designs from a large design space search in a Concept and Requirements Exploration (CandRE) are assessed using this methodology.
- Discretization Error Estimation Using the Error Transport Equations for Computational Fluid Dynamics SimulationsWang, Hongyu (Virginia Tech, 2021-06-11)Computational Fluid Dynamics (CFD) has been widely used as a tool to analyze physical phenomena involving fluids. To perform a CFD simulation, the governing equations are discretized to formulate a set of nonlinear algebraic equations. Typical spatial discretization schemes include finite-difference methods, finite-volume methods, and finite-element methods. Error introduced in the discretization process is called discretization error and defined as the difference between the exact solution to the discrete equations and the exact solution to the partial differential or integral equations. For most CFD simulations, discretization error accounts for the largest portion of the numerical error in the simulation. Discretization error has a complicated nonlinear relationship with the computational grid and the discretization scheme, which makes it especially difficult to estimate. Thus, it is important to study the discretization error to characterize numerical errors in a CFD simulation. Discretization error estimation is performed using the Error Transport Equations (ETE) in this work. The original nonlinear form of the ETE can be linearized to formulate the linearized ETE. Results of the two types of the ETE are compared. This work implements the ETE for finite-volume methods and Discontinuous Galerkin (DG) finite-element methods. For finite volume methods, discretization error estimates are obtained for both steady state problems and unsteady problems. The work on steady-state problems focuses on turbulent flow modelled by the Spalart-Allmaras (SA) model and Menter's $k-omega$ SST model. Higher-order discretization error estimates are obtained for both the mean variables and the turbulence working variables. The type of pseudo temporal discretization used for the steady-state problems does not have too much influence on the final converged solution. However, the temporal discretization scheme makes a significant difference for unsteady problems. Different temporal discretizations also impact the ETE implementation. This work discusses the implementation of the ETE for the 2-step Backward Difference Formula (BDF2) and the Singly Diagonally Implicit Runge-Kutta (SDIRK) methods. Most existing work on the ETE focuses on finite-volume methods. This work also extends ETE to work with the DG methods and tests the implementation with steady state inviscid test cases. The discretization error estimates for smooth test cases achieve the expected order of accuracy. When the test case is non-smooth, the truncation error estimation scheme fails to generate an accurate truncation error estimate. This causes a reduction of the discretization error estimate to first-order accuracy. Discussions are made on how accurate truncation error estimates can be found for non-smooth test cases.
- The Dynamics of Single and Double Cavitation Bubbles and Interaction Between Bubbles and Different MaterialsZhao, Ben (Virginia Tech, 2022-09-06)We present two distinct projects in this article. In the first project, an experiment aiming to quantify the impacts of material acoustic impedance and thickness on single laser-induced cavitation bubble dynamics with measurements of exerted pressure on a specific material in order to identify the primary sources most responsible for material damages is presented in this article. Two types of major pressure sources have been identified. For bubble collapsing near a rigid wall, when standoff ratio γ < 0.6, the ring collapse is the most prominent pressure source. The jet takes the strongest effects at γ = 1.12. The pressure is minimal at γ = 0.913. After the first jet impingement, a second ring collapse will follow and input the maximum pressure to the wall. By further increasing γ, a similar pressure profile of the second collapse to the first collapse is achieved, during which the pressure for the second collapse is minimal at γ = 1.41 and for the jet is maximum at γ = 1.79. Compared with the maximum pressure dealt by the first jet, the second ring collapse and jet are increasing much faster with the bubble size and eventually overwhelm the first jet. However, the first ring collapse is still the most dominant pressure source responsible for material damages. For wall featuring smaller acoustic impedance or thickness that cannot be approximated to a rigid body, the ring collapse and jet occur at smaller standoff ratios. The cavity shrinking rate suggests the maximum pressure exerted on the wall at applicable standoff ratios should be smaller than that on a rigid wall. In the second project, a comprehensive collection of dynamics of one and two laser-induced cavitation bubbles collapsing near different boundaries is presented in this article by measuring the velocity fields using particle image velocimetry (PIV) techniques. Cases include a single bubble collapsing near the hard, medium, and soft walls characterized by acoustic impedance, free collapse of two bubbles, and two bubbles collapsing near the hard and soft walls. We implemented the most significant velocity and top velocity regions derived from each velocity field to analyze the features of these cases. Before converging to free collapse, the bubble near the hard wall experienced a significant velocity decrease before collapse, the bubble near the medium wall was severely damped at a specific standoff distance, and the bubble near the soft wall collapsed much earlier and preserved a linear velocity region at low speed. Free collapse of two same bubbles underwent a decrease of acceleration before collapse. Decreasing the size of one bubble caused a jet in the other. With the presence of a hard wall near two bubbles, the bubble closer to it may be stretched to a cavity with a high aspect ratio, leading to very mild collapse. With a bigger bubble between a smaller one and the soft wall, the merging cavity may suppress the tendency of jet formation, making the velocity stay at low levels throughout the lifetime. For configurations regarding single bubbles collapsing near a wall and free collapse of two same bubbles, we performed data scaling to study the velocity variations for different bubble sizes by controlling the standoff ratios and assessed the data quality aided by curving fitting and statistics. Results indicated measured velocity regarding a single bubble collapsing near the wall over its diameter remained the same given a standoff ratio, while measured velocity did not change given a standoff ratio for free collapse of two same bubbles within the scope of the experiment. In addition, we detailed the experimental setup and water treatment for better signal-to-noise ratios as well as validated the system from both the PIV and high speed imaging approaches using free collapse of a single bubble to ensure the reliability of this experiment.
- Evaluating the OpenACC API for Parallelization of CFD ApplicationsPickering, Brent Phillip (Virginia Tech, 2014-09-06)Directive-based programming of graphics processing units (GPUs) has recently appeared as a viable alternative to using specialized low-level languages such as CUDA C and OpenCL for general-purpose GPU programming. This technique, which uses directive or pragma statements to annotate source codes written in traditional high-level languages, is designed to permit a unified code base to serve multiple computational platforms and to simplify the transition of legacy codes to new architectures. This work analyzes the popular OpenACC programming standard, as implemented by the PGI compiler suite, in order to evaluate its utility and performance potential in computational fluid dynamics (CFD) applications. Of particular interest is the handling of stencil algorithms, which are an important component of finite-difference and finite-volume numerical methods. To this end, the process of applying the OpenACC Fortran API to a preexisting finite-difference CFD code is examined in detail, and all modifications that must be made to the original source in order to run efficiently on the GPU are noted. Optimization techniques for OpenACC are also explored, and it is demonstrated that tuning the code for a particular accelerator architecture can result in performance increases of over 30%. There are also some limitations and programming restrictions imposed by the API: it is observed that certain useful features of modern Fortran (2003/8) are effectively disabled within OpenACC regions. Finally, a combination of OpenACC and OpenMP directives is used to create a truly cross-platform Fortran code that can be compiled for either CPU or GPU hardware. The performance of the OpenACC code is measured on several contemporary NVIDIA GPU architectures, and a comparison is made between double and single precision arithmetic showing that if reduced precision can be tolerated, it can lead to significant speedups. To assess the performance gains relative to a typical CPU implementation, the execution time for a standard benchmark case (lid-driven cavity) is used as a reference. The OpenACC version is compared against the identical Fortran code recompiled to use OpenMP on multicore CPUs, as well as a highly-optimized C++ version of the code that utilizes hardware aware programming techniques to attain higher performance on the Intel Xeon platforms being tested. Low-level optimizations specific to these architectures are analyzed and it is observed that the stencil access pattern required by the structured-grid CFD code sometimes leads to performance degrading conflict misses in the hardware managed CPU caches. The GPU code, which primarily uses software managed caching, is found to be free from these issues. Overall, it is observed that the OpenACC GPU code compares favorably against even the best optimized CPU version: using a single NVIDIA K20x GPU, the Fortran+OpenACC code is seen to outperform the optimized C++ version by 20% and the Fortran+OpenMP version by more than 100% with both CPU codes running on a 16-core Xeon workstation.
- From Oscillating Flat Plate to Maneuvering Bat Flight – Role of Kinematics, Aerodynamics, and InertiaRahman, Aevelina (Virginia Tech, 2022-02-01)With the aim to understand the synergistic roles played by kinematics, aerodynamics, and inertia in flapping wing maneuvers, this thesis first investigates the plunging motion of a simple flat plate as it is a fundamental motion in the kinematics of many flying animals. A wide range of frequency (k) and amplitude (h) is investigated to account for a robust kinematic characterization in the form of plunge velocity (kh). Leading Edge Vortices (LEVs) are found to be responsible for producing thrust while Trailing Edge Vortices (TEVs) produce drag. The vortex dynamics becomes nonlinear for higher kh and three main vortex-vortex interactions (VVI) are identified in the flow-field. To estimate the sole effect of LEVs on thrust coefficient, TEVs are eliminated by introducing a splitter plate. This resulted in reduced non-linearity in VVI and facilitated a parametrization of aerodynamic thrust coefficient with key kinematic features, frequency (k) and amplitude (h) [C_T= A.k^1.4 h-B where A and B are constants]. This is followed by investigating the more direct problem of bio-inspired MAV research – the interplay of kinematics, aerodynamics, and inertia on maneuvering bat flights. At first, an ascending right turn of a H. pratti bat is investigated to elucidate on the kinematic features and aerodynamic mechanisms used to effectuate the maneuver. Deceleration in flight speed, an increase in flapping frequency, shortening of the upstroke, and thrust generation at the end of the upstroke is observed during this maneuver. The turn is initiated by the synergisytic implementation of roll and yaw rotation where the turning moments are generated by drawing the inside wing closer to the body, by introducing phase lags in force generation between the two wings and by redirecting force production to the outer part of the wing outside of the turn. Upon comparison with a similar maneuver by a H. armiger bat, some commonalities as well as differences were observed. This analysis was followed by a comparative study among different maneuvering flights (a straight flight, two ascending right turns, and a U-turn) in order to establish the complete motion dynamics of a maneuver in action. The individual effects of aerodynamics and wing inertia for maneuvering flights of a H. armiger and H. pratti are investigated. It is found that for both, translation and rotation the overall trajectory trend is mostly driven by the aerodynamic forces and moments, whereas inertial effects drive the intricate intra-cycle fluctuations as well as the vertical velocity and altitude gain during ascent. Additionally, inertial moments play a dominant role for effecting yaw rotations where the importance of the Coriolis and centrifugal moments increase with increasing acuteness of the maneuver, with the largest effect of centrifugal moments being evidenced in the U-turn.
- A Hybrid Framework of CFD Numerical Methods and its Application to the Simulation of Underwater ExplosionsSi, Nan (Virginia Tech, 2022-02-08)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.
- Hydroelasticity of High-Speed Planing Craft Subject to Slamming Events: An Experimental and Numerical Investigation of Wedge Water EntryRen, Zhongshu (Virginia Tech, 2020-08-27)High-speed planing craft operating in waves are subject to frequent water impact, or slamming, as a portion or whole of the craft exits the water and re-enters at high velocity. The global load induced by slamming can cause fatigue-related damages to structures. The local slamming can cause local damage to structures and its induced acceleration can cause damage to equipment and personnel aboard. Therefore the slamming loads in high-speed craft are critical design loads. Nowadays, due to the increasing use of composite materials in high-speed craft, the interaction between the hydrodynamic loading and structural response, or hydroelasticity, must be considered. In this work, a flexible V-shaped wedge, which vertically enters the calm water with an impact velocity, was examined experimentally and numerically to characterize the slamming of a representative cross-section of high-speed craft. Physical quantities of interest include rigid-body kinematic motions, spray root propagation, hydrodynamic loading, and structural response. In the experimental work, with varied impact velocity and flexural rigidity of the wedge bottom plate, a wide range of hydroelasticity factors were investigated. The intersection between the bottom plate and side plate is called chine. The phases before and after the spray root reached the chine are called chine-unwetted and chine-wetted phase, respectively. It was found that the maximum deflection and strain occur in the chine-unwetted phase while a structural vibration with rapidly decaying magnitude is observed in the chine-wetted phase. Furthermore, the kinematic effect of hydroelasticity changes the spray root propagation and hence the pressure, while the inertial effect elongates the natural period of the plate. Inspired by the experimental work, a computational framework was proposed to focus on the chine-unwetted phase. Several hydroelastic models can be obtained from this framework. The hydroelastic models were validated to show reasonable agreement with experiments. Various parameters were studied through the computational framework. The hydroelasticity factor was modified to account for the mass and boundary conditions. It was found that the nondimensional rigid-body kinematic motions and maximum deflection showed little dependence on the hydroelasticity factor. Hydroelastic effects increased the time it takes for the peak maximum deflection to be reached for small values of the hydroelasticity factor. Hydroelastic effects also have little influence on the magnitude of the maximum deflection. These discoveries further the understanding of hydroelastic slamming and show the potential to guide the structural optimization and design of high-speed craft.
- Induction Infrared Thermography for Non-Destructive Evaluation of Alloy SensitizationRoberts, Matthew Thomas (Virginia Tech, 2019-06-26)The sensitization of stainless steel describes the process by which a high-carbon steel alloy is heated above a certain threshold (either naturally or artificially) followed by a cooling period during which chromium (one of the elements most responsible for providing stainless steel with its corrosion-inhibiting properties) forms new compounds with the carbon present in the steel. With the chromium being taken from the parent material to form these compounds, the corrosion-resistant properties are compromised, which can lead to corrosion, cracking, and broader failure. Currently, the accepted techniques used to test for the presence of sensitization are qualitative and/or destructive in nature. Attempts have been made to non-destructively detect and characterize sensitization through various means, but all with mixed results. With the use of these high-carbon alloys in a range of industries, a comprehensive, in-place process is desirable. This thesis will focus specifically on non-destructive evaluation of sensitization seen as a result of welding steel plates using induction infrared thermography (IIRT). This process uses an induction coil to generate heat within a sample whose resulting heat signature can then be detected with an infrared (IR) camera and analyzed. Previous IIRT experimental results have shown higher levels of heating in the HAZ when sensitization is present as it modifies the original microstructure of the material. New IIRT experiments have been conducted on both welded and unwelded 440C alloy samples to establish quantitative data on the heating profiles. These results (in conjunction with the appropriate experimental parameters) were then used to create a numerical model to replicate them. Despite some limitations in populating the model with accurate parameters, the results obtained were in good agreement with the experiments and provide a foundation for future work. Future work will focus on establishing a predictive tool that can detect and quantify the level of sensitization in an arbitrary steel sample in the field.
- An Investigation of Phase Change Material (PCM)-Based Ocean Thermal Energy HarvestingWang, Guangyao (Virginia Tech, 2019-06-10)Phase change material (PCM)-based ocean thermal energy harvesting is a relatively new method, which extracts the thermal energy from the temperature gradient in the ocean thermocline. Its basic idea is to utilize the temperature variation along the ocean water depth to cyclically freeze and melt a specific kind of PCM. The volume expansion, which happens in the melting process, is used to do useful work (e.g., drive a turbine generator), thereby converting a fraction of the absorbed thermal energy into mechanical energy or electrical energy. Compared to other ocean energy technologies (e.g., wave energy converters, tidal current turbines, and ocean thermal energy conversion), the proposed PCM-based approach can be easily implemented at a small scale with a relatively simple structural system, which makes it a promising method to extend the range and service life of battery-powered devices, e.g, autonomous underwater vehicles (AUVs). This dissertation presents a combined theoretical and experimental study of the PCM-based ocean thermal energy harvesting approach, which aims at demonstrating the feasibility of the proposed approach and investigating possible methods to improve the overall performance of prototypical systems. First, a solid/liquid phase change thermodynamic model is developed, based on which a specific upperbound of the thermal efficiency is derived for the PCM-based approach. Next, a prototypical PCM-based ocean thermal energy harvesting system is designed, fabricated, and tested. To predict the performance of specific systems, a thermo-mechanical model, which couples the thermodynamic behaviors of the fluid materials and the elastic behavior of the structural system, is developed and validated based on the comparison with the experimental measurement. For the purpose of design optimization, the validated thermo-mechanical model is employed to conduct a parametric study. Based on the results of the parametric study, a new scalable and portable PCM-based ocean thermal energy harvesting system is developed and tested. In addition, the thermo-mechanical model is modified to account for the design changes. However, a combined analysis of the results from both the prototypical system and the model reveals that achieving a good performance requires maintaining a high internal pressure, which will complicate the structural design. To mitigate this issue, the idea of using a hydraulic accumulator to regulate the internal pressure is proposed, and experimentally and theoretically examined. Finally, a spatial-varying Robin transmission condition for fluid-structure coupled problems with strong added-mass effect is proposed and investigated using fluid structure interaction (FSI) model problems. This can be a potential method for the future research on the fluid-structure coupled numerical analysis of AUVs, which are integrated with and powered by the PCM-based thermal energy harvesting devices.
- Long-Pulsed Laser-Induced Cavitation: Laser-Fluid Coupling, Phase Transition, and Bubble DynamicsZhao, Xuning (Virginia Tech, 2024-02-29)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.
- Modeling Underwater Explosion (UNDEX) Shock Effects for Vulnerability Assessment in Early Stage Ship DesignMathew, Ajai Kurian (Virginia Tech, 2018-03-20)This thesis describes and assesses a simplified tool for modeling underwater explosion shock effects during early naval ship concept design. A simplified fluid model using Taylor flat-plate theory is incorporated directly into the OpenFSI module code in Nastran and used to interface with the structural solver in Nastran to simulate a far-field shockwave impacting the hull. The kick-off velocities and the shock spectra captured in this computationally efficient module is compared to results from a high-fidelity CASE (Cavitating Acoustic Spectral Element) fluid model implemented with the ABAQUS/Nastran structural solver to validate the simplified framework and assess the sufficiency of this very simple but, fast approach for early stage ship design.
- Multiphase Fluid-Material Interaction: Efficient Solution Algorithms and Shock-Dominated ApplicationsMa, Wentao (Virginia Tech, 2023-09-05)This dissertation focuses on the development and application of numerical algorithms for solving compressible multiphase fluid-material interaction problems. The first part of this dissertation is motivated by the extraordinary shock-resisting ability of elastomer coating materials (e.g., polyurea) under explosive loading conditions. Their performance, however, highly depends on their dynamic interaction with the substrate (e.g., metal) and ambient fluid (e.g., air or liquid); and the detailed interaction process is still unclear. Therefore, to certify the application of these materials, a fluid-structure coupled computational framework is needed. The first part of this dissertation developes such a framework. In particualr, the hyper-viscoelastic constitutive relation of polyurea is incorporated into a high-fidelity computational framework which couples a finite volume compressible multiphase fluid dynamics solver and a nonlinear finite element structural dynamics solver. Within this framework, the fluid-structure and liquid-gas interfaces are tracked using embedded boundary and level set methods. Then, the developed computational framework is applied to study the behavior a bilayer coating–substrate (i.e., polyurea-aluminum) system under various loading conditions. The observed two-way coupling between the structure and the bubble generated in a near-field underwater explosion motivates the next part of this dissertation. The second part of this dissertation investigates the yielding and collapse of an underwater thin-walled aluminum cylinder in near-field explosions. As the explosion intensity varies by two orders of magnitude, three different modes of collapse are discovered, including one that appears counterintuitive (i.e., one lobe extending towards the explosive charge), yet has been observed in previous laboratory experiments. Because of the transition of modes, the time it takes for the structure to reach self-contact does not decrease monotonically as the explosion intensity increases. Detailed analysis of the bubble-structure interaction suggests that, in addition to the incident shock wave, the second pressure pulse resulting from the contraction of the explosion bubble also has a significant effect on the structure's collapse. The phase difference between the structural vibration and the bubble's expansion and contraction strongly influences the structure's mode of collapse. The third part focuses on the development of efficient solution algorithms for compressible multi-material flow simulations. In these simulations, an unresolved challenge is the computation of advective fluxes across material interfaces that separate drastically different thermodynamic states and relations. A popular class of methods in this regard is to locally construct bimaterial Riemann problems, and to apply their exact solutions in flux computation, such as the one used in the preceding parts of the dissertation. For general equations of state, however, finding the exact solution of a Riemann problem is expensive as it requires nested loops. Multiplied by the large number of Riemann problems constructed during a simulation, the computational cost often becomes prohibitive. This dissertation accelerates the solution of bimaterial Riemann problems without introducing approximations or offline precomputation tasks. The basic idea is to exploit some special properties of the Riemann problem equations, and to recycle previous solutions as much as possible. Following this idea, four acceleration methods are developed. The performance of these acceleration methods is assessed using four example problems that exhibit strong shock waves, large interface deformation, contact of multiple (>2) interfaces, and interaction between gases and condensed matters. For all the problems, the solution of bimaterial Riemann problems is accelerated by 37 to 87 times. As a result, the total cost of advective flux computation, which includes the exact Riemann problem solution at material interfaces and the numerical flux calculation over the entire computational domain, is accelerated by 18 to 81 times.
- The Numerical Investigation of the Effects of Sand Ingestion on Compressor Blade ErosionCagdas, Taha Irfan (Virginia Tech, 2024-01-10)ABSTRACT The performance of aircraft engines can be significantly affected by the variety of foreign particles that are mixed into the air while operating under miscellaneous conditions. In particular, aircraft engines that operate in sandy or dusty conditions may fail within minutes of exposure to particle-laden flow due to foreign particle deposition on hot section components or erosion occurring on the compressor and turbine blades. For these reasons, the effect of sand ingestion on erosion, which may occur in the turbine and compressor blades, was studied in this master's thesis. In this master's thesis, the effect of sand ingestion on erosion on the M250 turboshaft engine's compressor blades will be investigated with the aid of numerical methods. In this study, we used the OpenFOAM software to solve the multiphase flow problem from the standpoint of finite control methods and the Eulerian-Lagrangian framework. The initial sand distribution conditions were taken from the Ph.D. thesis written by Olshefski, K. T. (2023) [1]. The compressor blade was modeled as 2D, which has a NACA 6510 profile shape, with a chord length of 63 mm. The results show that the leading edge and the suction side of the compressor, i.e. the upper half of the compressor, eroded more compared to the trailing edge, and the pressure side. Results also show that as the sand particle distribution becomes non-uniform the most eroded region shifts toward the trailing edge. In addition, for varying angles of attack, the region where the erosion occurs alters periodically. We observed that as the angle of attack increases, the eroded region shifts toward the trailing edge, but when the angle of attack is kept increasing the eroded region shifts back to the leading edge again. In conclusion, the non-uniformity of sand particle loading has a strong effect on the determination of the eroded regions. Furthermore, the variation of the angle of attack has a huge role in both the determination of eroded regions and the amount of eroded material.
- Numerical Methods for Fluid-Solid Coupled Simulations: Robin Interface Conditions and Shock-Dominated ApplicationsCao, Shunxiang (Virginia Tech, 2019-09-09)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.