Organizational Unit:

Daniel Guggenheim School of Aerospace Engineering

Permanent Link

https://hdl.handle.net/1853/70742

Parent Organization

Organizational Unit

College of Engineering

Includes Organization(s)

Organizational Unit

Aerospace Design Group

Organizational Unit

Aerospace Systems Design Laboratory (ASDL)

Organizational Unit

Air Transportation Laboratory

Organizational Unit

Cognitive Engineering Center

Organizational Unit

Space Systems Design Laboratory (SSDL)

Show 1 more

ArchiveSpace Name Record

https://finding-aids.library.gatech.edu/agents/corporate_entities/106

Full item page

Publication Search Results

Now showing 1 - 10 of 60

Evaluation of Convolutional Neural Networks for Modeling Blast Propagation in Multi-room Bunkers

(Georgia Institute of Technology, 2023-12-15) Luo, Felix

The rapid evaluation of blasts in enclosed geometrically complex spaces has long eluded the design of safer blast-resistant structures. Traditional methods of determining blast responses in enclosed geometrically complex spaces oftentimes rely on the use of traditional computational fluid dynamics (CFD) solvers to compute the entire flow field of the structure. This method has an enormous computational burden, especially considering that blasts are highly transient in nature and require the transient pressure fluctuations to be determined to formulate an accurate blast response prediction. However, more efficient methods of blast evaluation are desired such that parametric sweeps or optimization processes can be performed at low cost to provide a tool for iterative design. To rectify this gap in capabilities, a convolutional neural network based (CNN) model was developed to provide rapid blast predictions for 2D structures to establish this capability to aid in the design of more blast resistant structures. This approach leverages the inherent spatial awareness of CNNs to provide predictions for peak pressures since blasts in enclosed spaces are highly dependent on the spatial relationships between blast locations and wall location. This approach provides a nearly 5,000 times speed up against CFD simulations used within this study with good convergence of errors, correlation coefficients, predicted and truth values and distributions in all situational evaluations. These computational advantages, in part, comes from using the CNN based model to directly predict peak pressures whereas traditional CFD solvers require iterations to propagate fluid flows over time. However, some limitations do exist with respect to higher errors, such as model training costs, and the capability to predict 3D structures. Nonetheless, the results provide a characterization of the capabilities CNN based models in predicting peak pressures from blasts in enclosed spaces. From these evaluations and studies, a model which can provide significant computational savings while maintaining a similar accuracy can be obtained, which enables the rapid iterative design of more blast resistant structures.
Reshock Gas Curtain Mixing Study

(Georgia Institute of Technology, 2022-06-25) Risley, Karl Robert

The current work investigates the behavior of gas curtain instabilities. A gas curtain can be visualized as an A − B − A domain, where A and B are light and heavy fluids respectively, creating a ”curtain” of heavy fluid B that is surrounded by a light fluid A. Specifically, the behavior of gas curtains following an initial shock passage and the passage of a reflected shock (reshock) through the entirety of the curtain are investigated. A gas curtain instability commonly occurs physically in a wide range of applications such as during afterburning of an explosion, inertial confinement fusion, and even supernovae explosions. Previous studies have emphasized that the physics occurring during the reshock of a gas curtain are far more complex than the behavior of a re-shock Richtmyer-Meshkov Instability, due to the interactions between the two interfaces and wave reverberations occurring. The current work attempts to understand the relationship between a gas curtain’s initial conditions and its behavior to reshock through two-dimensional numerical simulations that utilize the viscous Navier-Stokes equations. More specifically, the current work isolates the effects of the curtain’s initial thickness and shape on the post reshock mixing layer growth rate and molecular mixing of the curtain. The results for all cases indicate that the post-reshock growth rate of the curtain’s width is a function of initial thickness. The sensitivity of the curtain’s post-reshock growth rate to the initial thickness, however, depends on the curtain’s initial perturbation shape. As the initial thickness of the curtain is decreased, the interactions between the curtain’s interfaces grow in strength and impede perturbation growth, thus reducing the post reshock growth rate of the curtain’s structure width. Similarly, the results strongly suggest that a reduction in initial curtain thickness increases the late-time asymptotic molecular mixing fraction value. This result is significant, especially for reacting flows, because it indicates that faster combustion (or afterburning in an explosion) could be reached with the thinning of the gas curtain in flow systems.
UNCERTAINTY QUANTIFICATION OF DIMP PYROLYSIS KINETICS

(Georgia Institute of Technology, 2022-05-13) Patel, Pavan

To develop effective explosives and strategies for the rapid destruction of sarin stockpiles, a reliable understanding of sarin’s chemical kinetics is needed. Kinetic mechanisms of sarin simulants such as di-isopropyl methyl phosphonate (DIMP) are developed instead because they have a similar chemical structure as sarin and are less toxic. A detailed DIMP kinetics mechanism has been developed in the past; however, there is a considerable amount of uncertainty surrounding it. This uncertainty manifests through the choice of pathways, and their respective reaction rates, leading to large variations in outcomes predicted through simulations. Out of the many reaction pathways involved in the decomposition of DIMP, the initiating steps are the most crucial. Out of the two possible initiating pathways in the destruction of DIMP, the lower activation energy pathway is dominant for all temperatures. The purpose of this study is to investigate the uncertainties associated with the dominant initiating pathways of the DIMP kinetics mechanism. Propagating rate parameter uncertainties of the dominant pathway through computational models yields large uncertainties in predicting DIMP survivability at different temperatures. The prediction uncertainties are larger at lower temperatures than at high temperatures. This can significantly impact the ability to precisely predict collateral damage caused by partially destroyed DIMP in the far-field of an explosion. After reducing these rate parameter uncertainties, using Bayesian inference, the prediction uncertainties were within reasonable limits. The results here provide a reduced subspace for uncertainties associated with the first and most important step in the breakdown of DIMP, which shall enable more reliable predictions.
Modeling moderately dense to dilute multiphase reacting flows

(Georgia Institute of Technology, 2022-03-14) Panchal, Achyut

Computational modeling of multiphase flows consists of two broader sets of methods: resolved approaches where the multiphase entities (MPE) are larger and resolved on the computational grid, and dispersed approaches where the MPEs are relatively small and not resolved on the grid but they are treated as point-particles that interact with the background continuum. A consistent multiphase formulation is developed in this work that can be used to model either resolved or unresolved MPEs over a complete range of volume fractions. One phase is always modeled as a continuum Eulerian phase, whereas the other phase is modeled either as a continuum Eulerian phase, a dispersed Eulerian phase, or a dispersed Lagrangian phase. In the dispersed phase limit, a hybrid EE-EL formulation is developed from first principles, which asymptotes to well-established EE and EL methods in limiting conditions. A smooth and dynamic transition criterion and a corresponding algorithm for conversion between EE and EL are developed. To use EE and EL in their respective regions of effectiveness (dense and dilute, respectively), the transition criterion is designed as a function of the local volume fraction and the local kinetic energy of random uncorrelated motion of particles. Simulations of particle evolution in turbulence, particle dispersion in sector blast, and reactive spray jet show the method’s validity and practical relevance. In the resolved multiphase limit, the formulation limits to a compressible seven-equation diffused interface method (DIM). Novel extensions are developed for modeling surface tension, viscous effects, arbitrary EOS, multi-species, and reactions. The use of a discrete equations method (DEM) relieves the need to use conventional stiff relaxation solvers. Shock propagation through a material interface, surface tension-driven oscillating droplet, droplet acceleration in a viscous medium, and shock/detonation interaction with a deforming droplet are simulated to validate various part of the computational framework and demonstrate its applicability.
Investigating Lean Blowout of an Alternative Jet Fuel in a Gas Turbine Combustor

(Georgia Institute of Technology, 2022-01-05) Narayanan, Vijay

In the global effort to reduce the climate impact of combustion emissions, sustainable aviation fuels offer the ease and reliability of conventional petroleum-derived jet fuels without the significant pollutant effects. Ongoing research efforts include experimental testing of alternative jet fuels to identify fuel candidates that produce less pollutant combustion products and are cheaper and environmentally cleaner to source than conventional jet fuels. Fuel lean combustion already reduces the emissions of jet engines and increases fuel efficiency, but lean blowout (LBO) can occur at reduced throttle and minimum power scenarios such as descent. Lean blowout (LBO) has been identified as a critical figure of merit to ensure the stability of alternative jet fuels in the place of conventional fuels. This work aimed to further understand the LBO phenomenon, leveraging computational studies of the alternative fuel designated C-5 by the National Jet Fuel Combustion Program (NJFCP). The fuel sensitivity of LBO has been established by the NJFCP’s participants recently. In this thesis, the chemical kinetics for C-5 is first verified using zerodimensional (0-D) and one-dimensional (1-D) studies and then this is followed by three dimensional (3D) large-eddy simulations (LES). In LES to observe LBO, a direct-step and gradual equivalence ratio reduction were separately employed to assess fuel sensitivity of LBO against available experimental data. The time histories of pressure, temperature, and composition were analyzed for precursor signatures of LBO both inside and outside the flame. Localized extinction, a reduction in the vortex breakdown bubble size and magnitude, and a reduction in the exhaust velocity were all observed to occur during the LBO event.
LES of Turbulent Premixed Flame Kernel Formation and Development

(Georgia Institute of Technology, 2020-12-17) Lambert, Alexander

Spark ignition of flammable mixtures is highly sensitive to early and local conditions. Kernel formation and subsequent flame development are largely governed by turbulent conditions and interactions with igniter geometry. In order to investigate this phenomenon, the use of Large Eddy Simulation (LES) is examined for (1) modelling spherical turbulent flame development, and (2) simulating spark ignition in a channel with either laminar or turbulent inflow. A comparison between LES spherical flame simulation is made to FSD-LES results as well as experimental measurements from previous studies. For spark ignition experiments, we characterize the temporal evolution of the ignition process, and demonstrate the dependence on early velocity fluctuations and local conditions.
Entrainment, Mixing, and Ignition in Single and Multiple Jets in a Supersonic Crossflow

(Georgia Institute of Technology, 2020-08-19) Fries, Dan

Jets in crossflow are a canonical example for three-dimensional turbulent mixing. Here, non-reacting and reacting sonic jets in a supersonic crossflow are studied. The influence of injectant properties on turbulent mixing is investigated. Using pure gases, the molecular weight and specific heat ratio is varied between 4-44 g/mol and 1.24-1.66, respectively. The jets are injected into a Mach 1.71 crossflow with a stagnation temperature ~600 K. Two single jet injectors and two staged jet injectors are designed to characterize potential enhancements in turbulent mixing and combustion processes. Mixture fraction and velocity fields are determined via Mie-scattering off solid particles. Velocity vectors are obtained by processing Mie-scattering image pairs with a correlation technique (particle image velocimetry). To ignite the flow field and enable systematic variation of the ignition location a traversable laser spark system is employed. The reacting flow is probed via CH* chemiluminescence and OH planar laser induced fluorescence visualizing regions containing hot combustion products. A new trajectory scaling improves correlation between all data sets considered, suggesting that the bow shock, boundary layer and momentum flux ratio are the dominant controlling factors. Turbulent mixing rates are highest for injectants with higher molecular weight and lower specific heat ratio. The larger of two jet spacings tested yields the greater enhancement of turbulent mixing rates. Ignition locations on the symmetry plane of the flow field are evaluated for their ability to sustain chemical reactions/heat release. Most favorable ignition locations lie in the windward jet shear layer away from the regions of highest flow strain. The smallest diameter single jet with presumably more boundary layer interaction and moderate strain rates provides the best results with regard to thermal energy release after spark deposition. Trends suggest that moderate compressible strain rates and no flow expansion are advantageous to sustain thermal energy release. Implications for future research directions and opportunities are discussed.
Numerical simulations of real-gas flows with phase-equilibrium thermodynamics

(Georgia Institute of Technology, 2020-07-20) Tudisco, Principio

Motivated by the complex physics of multi-component mixtures in strongly non-ideal, real-gas (RG) conditions reported in the field of chemical engineering, this work aims to address the behavior of multi-phase thermodynamics from a broader point of view. The focus is to evaluate the differences, as well as the possible sources of errors that would arise in a computational fluid-dynamics (CFD) simulation when conventional single-phase and multi-phase equilibrium RG thermodynamics are employed: an area of research that despite the active interest in many communities (especially CFD), has not been completely understood. Knowledge of the effects that multi-phase RG thermodynamics with the assumption of vapor-liquid equilibrium (VLE) can have on a flow dynamics is important because it establishes the relevance of the fully coupled CFD-VLE solver. In fact, this relevance may go beyond the stand-alone calculation of a multi-phase state, providing important insights about the physics that may not be captured if the single-phase assumption is invoked. This work provides an extensive study of RG mixtures from a physical and numerical point of view. The difficulties associated with their modeling are discussed in detail and solutions are provided accordingly. Emphasis is given to the occurrence and suppression of numerical noise in form of pressure oscillations that can pollute the simulation to the point that it cannot be performed. Extension of existing models to eliminate such problem is achieved by incorporating the effects of VLE thermodynamics in a consistent manner, ultimately forming a new and robust tool to investigate the physics further. The resulting model is applied to non-reacting and reacting flows in canonical setups where emphasis is devoted to the discussion of the differences and sources of errors that would occur if this multi-phase behavior is not taken into account. Results show that the different thermodynamic states reached by this advanced model can have an impact on the flow physics, especially in a non-reacting (or more in general cold) regime. In particular, the strong non-linear coupling between the VLE thermodynamics and the transport properties is identified as a key element of difference with respect to the single-phase model counterpart. These differences manifest into the occurrence of localized changes in the fluid properties (such as density) that affect the flow-field in their vicinity, causing visible discrepancies even when time-averaging is performed. Concurrently, results obtained on the reacting side and carried out (for the first time) with finite-rate kinetics suggest that any VLE formation between the products and the reactants may be considered of minor importance. The latter conclusion is supported by the analysis conducted on the multi-phase field which appears to be largely composed of the vapor solution, as expected, hence limiting the analogous effect observed the non-reacting system where a broader range of phase-separation appears instead.
Ignition, topology, and growth of turbulent premixed flames in supersonic flows

(Georgia Institute of Technology, 2019-11-12) Ochs, Bradley Alan

Supersonic combustion ramjets (scramjets) are currently the most efficient combustor technology for air breathing hypersonic flight, however, lack of fundamental understanding and numerous engineering challenges hinder regular deployment of these devices. This work addresses scramjet-relevant knowledge gaps in supersonic turbulent premixed combustion, including laser ignition, numerical modeling, and flame-compressibility interaction. One of the main contributions of this work is introduction of a new turbulent premixed flame arrangement where flame-compressibility interaction can be systematically explored: flame kernels in an expanding flow field. The scramjet flow path is replaced by a simplified channel geometry with a well characterized mean flow acceleration that mimics flow field expansion typically imposed on scramjet combustors to avoid thermal choking. Spherically expanding flames are created via laser ignition and subsequent flame growth and morphology are investigated using combined physical and numerical experiments. Pressure-density misalignment due to flame-compressibility interaction produces vorticity at the flame surface through baroclinic torque, i.e. flame-compressibility interaction acts like a turbulence source. The flame ultimately evolves into a reacting vortex ring that increases the flame speed and enhances reactant consumption. To explore the relative importance of turbulence and compressibility on flame dynamics, the Mach number (M=1.5,1.75,2), equivalence ratio (φ= 1.0,0.9,0.8,0.7), and root-mean-squared turbulent velocity (u'=3.98,4.14,4.45 m/s) are varied systematically. This work also introduces flame kernels in an expanding flow field as a canonical numerical validation test case for flame-compressibility interaction. Inaccuracies in simulation results are easily identified due to high flow velocity and simplicity of the problem. The numerical setup and models are scrutinized to minimize errors. Using the appropriately verified numerical models, simulation results show very reasonable agreement with experimental data. Validated simulations are instrumental in enhancing understanding of the underlying physics of supersonic flame kernels. Laser ignition studies in supersonic flows have historically focused on ignition of non-premixed fuels within cavity flame holders. This work introduces a far simpler and more tractable problem: laser ignition of a fully premixed supersonic gas. Ignition experiments with a range of laser settings are performed to determine supersonic breakdown and ignition probabilities, length of time the ignition event influences flame growth, and Mach number influence on the ignition process. The ignition event has a long-lasting effect on kernel growth, but the influence can be minimized by properly selecting the laser energy. Mach number has a minimal impact on the ignition process, but does affect the initial kernel shape due to flow field variations with Mach number. Kernel growth matches low speed studies closely at early times, but deviates at later times due to vortex ring topology. It is not obvious how the turbulent flame speed will scale for flows with mean compressibility. Therefore, the combined physical and numerical experiments are leveraged to explore this question. The vortex ring causes significant errors in the line of sight-measured burned volume, hence correction factors to convert from line of sight to volumetric measurements are presented. Conditions for displacement and consumption speed equivalency are shown to depend heavily on the particular diagnostic used; which progress variable isocontour is measured and where it is measured within the flame brush must be considered carefully during interpretation of experimental data. Scaling with the RMS turbulent velocity cannot collapse these flame speed data, i.e. previously established flame speed scalings are inappropriate for flames interacting with compressibility. Drawing motivation from vortex ring literature, a new flame speed scaling based on the ring propagation velocity is proposed. The proposed scaling collapses the data and produces a nearly linear scaling regime, which suggests turbulence plays a secondary role to the hydrodynamic instability created by flame-compressibility interaction. In summary, flame kernels are a new and effective canonical configuration for exploring flame-compressibility interactions in supersonic flows.
Investigation of ODE-based non-equilibrium wall shear stress models for large eddy simulation

(Georgia Institute of Technology, 2019-07-30) Dzanic, Tarik

For high Reynolds number flows, wall modeling is essential for performing large eddy simulation at a reasonable computational cost. In this work, a novel low-cost ODE-based non-equilibrium wall model is introduced for wall shear stress modeling in LES. Using polynomial approximations of the pressure gradient and convective terms obtained from interpolation of the LES solution, as opposed to direct evaluation of these gradients within the wall model, the governing wall model equations reduce from coupled PDEs to uncoupled ODEs that do not require an embedded wall model grid within the LES grid. Additionally, the steady form of the wall model equations was utilized, feasible due to the spatial decoupling of the wall model equations, and the effects of the temporal evolution on the wall shear stress were modeled. The effects of polynomial degree on the accuracy of the wall shear stress predictions were explored, and an empirical lag model was built to model the unsteady effects without requiring the solution of a time-stepping problem. Wall resolved large eddy simulations of separated flow around the NASA wall mounted hump and an iced NACA 63A213 airfoil were performed and used as a reference for the comparison of the non-equilibrium wall model to a commonly used equilibrium wall model. The proposed non-equilibrium wall model was able to predict separated flow and laminar flow regions in much better agreement with the wall resolved results than the equilibrium wall model. Underpredictions in the skin friction coefficient in non-equilibrium flow regimes were reduced from 20-50% to less than 10% between the equilibrium and the non-equilibrium wall modeled approaches. Minor improvements in the pressure coefficient predictions were observed with the non-equilibrium model in the separated flow region of the iced airfoil. The results suggest that the proposed wall model can offer better predictions of separated and/or laminar flows compared to equilibrium wall models with negligible computational cost increase.