Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “robust parameter design”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37

$\mathrm{SageNet}$: Fast Neural Network Emulation of the Stiff-amplified Gravitational Waves from Inflation

Accurate modeling of the inflationary gravitational waves (GWs) requires time-consuming, iterative numerical integrations of differential equations to take into account their backreaction on the expansion history. To improve computational efficiency while preserving accuracy, we present the Stiff-amplified Gravitational-wave Emulator Network (SageNet), a deep learning framework designed to replace conventional numerical solvers (code available at https://github.com/YifangLuo/SageNet). SageNet employs a long short-term memory architecture to emulate the present-day energy density spectrum of the inflationary GWs with possible stiff amplification, Ω GW (f). Trained on a data set of 25,689 numerically generated solutions, SageNet allows accurate reconstructions of Ω GW (f) and generalizes well to a wide range of cosmological parameters; 90.9% of the test emulations with randomly distributed parameters exhibit errors of under 4%. In addition, SageNet demonstrates its ability to learn and reproduce the artificial, adaptive sampling patterns in numerical calculations, which implement denser sampling of frequencies around changes in spectral indices in Ω GW (f). The dual capability of learning both physical and artificial features of the numerical GW spectra establishes SageNet as a robust alternative to exact numerical methods. Finally, our benchmark tests show that SageNet reduces the computation time from tens of seconds to milliseconds, achieving a speedup of ∼10 4 times over standard CPU-based numerical solvers with the potential for further acceleration on GPU hardware. These capabilities make SageNet a powerful tool for accelerating Bayesian inference procedures for extended cosmological models. In a broad sense, the SageNet framework offers a fast, accurate, and generalizable solution to modeling cosmological observables whose theoretical predictions demand costly differential equation solvers.

Astronomy data modeling↗

Multilevel Error Estimation and Adaptive h-Refinement for Cartesian Meshes with Embedded Boundaries

This paper presents the development of a mesh adaptation module for a multilevel Cartesian solver. While the module allows mesh refinement to be driven by a variety of different refinement parameters, a central feature in its design is the incorporation of a multilevel error estimator based upon direct estimates of the local truncation error using tau-extrapolation. This error indicator exploits the fact that in regions of uniform Cartesian mesh, the spatial operator is exactly the same on the fine and coarse grids, and local truncation error estimates can be constructed by evaluating the residual on the coarse grid of the restricted solution from the fine grid. A new strategy for adaptive h-refinement is also developed to prevent errors in smooth regions of the flow from being masked by shocks and other discontinuous features. For certain classes of error histograms, this strategy is optimal for achieving equidistribution of the refinement parameters on hierarchical meshes, and therefore ensures grid converged solutions will be achieved for appropriately chosen refinement parameters. The robustness and accuracy of the adaptation module is demonstrated using both simple model problems and complex three dimensional examples using meshes with from 10(exp 6), to 10(exp 7) cells.

Aftosmis, M. J.↗

Environmental exposure to industrial air pollution is associated with decreased male fertility

Objective: To understand how chronic exposure to industrial air pollution is associated with male fertility through semen parameters. Design: Retrospective cohort study. Subjects: Men in the Subfertility, Health and Assisted Reproduction cohort who underwent a semen analysis 2005-2017 with ≥1 measured semen parameter (N=21,563). Intervention(s): Residential histories for each man were constructed using locations from administrative records linked through the Utah Population Database. Industrial facilities with air emissions of nine endocrine disrupting compound chemical classes were identified from the Environmental Protection Agency Risk-Screening Environmental Indicators microdata. Chemical levels were linked with residential histories for the 5 years prior to each semen analysis. Main Outcome Measures: Semen analyses were classified as azoospermic or oligozoospermic (< 15 M/mL) using World Health Organization cutoffs for concentration. Bulk semen parameters such as concentration, total count, ejaculate volume, total motility, total motile count, and total progressive motile count were also measured. Multivariable regression models with robust standard errors were used to associate exposure quartiles for each of the nine chemical classes with each semen parameter, adjusting for age, race, and ethnicity, as well as neighborhood socioeconomic disadvantage. Results: After adjustment for demographic covariates, several chemical classes were associated with azoospermia and decreased total motility and volume. For exposure in the 4th relative to 1st quartile, significant associations were observed for acrylonitrile (β total motility = -0.87 pp), aromatic hydrocarbons (odds ratio [OR]azoospermia = 1.53; β volume = -0.14 mL), dioxins (OR azoospermia = 1.31; β volume = -0.09 mL; β total motility = -2.65 pp), heavy metals (β total motility = -2.78pp), organic solvents (OR azoospermia = 1.75; β volume = -0.10 mL), organochlorines (OR azoospermia = 2.09; β volume = -0.12 mL), phthalates (OR azoospermia = 1.44; β volume = -0.09 mL; β total motility = -1.21 pp), and silver particles (OR azoospermia = 1.64; β volume = -0.11 mL). All semen parameters significantly decreased with increasing socioeconomic disadvantage. Men who lived in the most disadvantaged areas had concentration, volume, and total motility of 6.70 M/mL, 0.13 mL, and 1.79 pp lower, respectively. Count, motile count, and total progressive motile count all decreased by 30–34 M.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Jipole: A Differentiable ipole-based Code for Radiative Transfer in Curved Spacetimes

Recent imaging of supermassive black holes by the Event Horizon Telescope has relied on exhaustive parameter-space searches, matching observations to large, precomputed libraries of theoretical models. As observational data become increasingly precise, the limitations of this computationally expensive approach grow more acute, creating a pressing need for more efficient methods. In this work, we present Jipole, an automatically differentiable (AD), ipole-based code for radiative transfer in curved spacetimes, designed to compute image gradients with respect to underlying model parameters. These gradients quantify how parameter changes—such as the black hole’s spin or the observer’s inclination—affect the image, enabling more efficient parameter estimation and reducing the number of required images. We validate Jipole against ipole in two analytical tests and then compare pixelwise intensity derivatives from AD with those from finite-difference methods. We then demonstrate the utility of these gradients by performing parameter recovery for an analytical model in three increasingly complex cases for the injected image: ideal, blurred, and blurred with added noise. In most cases, high-accuracy fits are obtained in only a few optimization steps, failing only in cases with extremely low signal-to-noise ratios. These results highlight the potential of AD-based methods to accelerate robust, high-fidelity model-data comparisons in current and future black hole imaging efforts.

79 ASTRONOMY AND ASTROPHYSICS↗

The Third Flight of the Marshall Grazing Incidence X-Ray Spectrometer (MaGIXS-3)

The Marshall Grazing Incidence X-ray Spectrometer (MaGIXS) is the first X-ray slitless imaging spectrograph sounding rocket instrument designed to observe spectrally dispersed soft X-ray images of the solar corona over a wide field-of-view. During the first flight of MaGIXS (MaGIXS-1), occurred on 30 July 2021, several emission lines from coronal structures including X-ray bright points were observed. Further, MaGIXS-1 analysis also demonstrated the successful inversion of overlappograms using robust unfolding algorithms. Given the demonstrated success of MaGIXS-1, the second flight of the instrument with a simplified optical design, MaGIXS-2, is scheduled for 2024 to observe high temperature diagnostic emission lines. Results from MaGIXS-1 discovered dominant missing emission lines near 15A, arising from relatively cool plasma that peaks around 2 MK. This wavelength region hosts several closely spaced satellite lines of Fe XVII, Fe XVI and Fe XV ions, which are expected to be enhanced at lower temperatures and are currently unmodeled in the CHIANTI atomic database. This wavelength region offers one of the most unique diagnostics to measure electron temperature, as well as signatures of equilibrium state of the plasma in active regions, which has not been studied so far. The goal of MaGIXS-3 mission is the to determine coronal heating parameters, such as the spatial and temporal properties of coronal heating events, by measuring discriminating observations, such as electron temperature, effective temperature, density, abundance, and departures from thermal equilibrium of the plasma, by observing the Sun in the SXR wavelength range. To meet this goal, the MaGIXS instrument will be upgraded with a new X-ray telescope mirror that will provide higher spatial and spectral resolution and throughput, allowing for spectral lines to be observed at the relevant spatial and temporal scales. In addition, MaGIXS-3 will also carry The Resolving Inversion Context X-ray Spectrometer (TRICXS), a high dispersion Bragg crystal spectrometer to spectrally resolve the lines near 14.9 to 15.9A, critical to unlock the full diagnostic potential of this wavelength range. Here we will present the preliminary concept design of MaGIXS-3 and discuss the potential observations.

X-ray Imaging↗

Topological valley Hall polariton condensation

A photonic topological insulator features robust directional propagation and immunity to defect perturbations of the edge/surface state. Exciton-polaritons, that is, the hybrid quasiparticles of excitons and photons in semiconductor microcavities, have been proposed as a tunable nonlinear platform for emulating topological phenomena. However, mainly due to excitonic material limitations, experimental observations so far have not been able to enter the nonlinear condensation regime or only show localized condensation in one dimension. Here, in this study, we show a topological propagating edge state with polariton condensation at room temperature and without any external magnetic field. We overcome material limitations by using excitonic CsPbCl 3 halide perovskites with a valley Hall lattice design. The polariton lattice features a large bandgap of 18.8 meV and exhibits strong nonlinear polariton condensation with clear long-range spatial coherence across the critical pumping density. The geometric parameters and material composition of our nonlinear many-body photonic system platform can in principle be tailored to study topological phenomena of other interquasiparticle interactions. Polariton condensation of topological propagating edge states is demonstrated with halide perovskite microcavities.

42 ENGINEERING↗

A Metallic Fuel Performance Benchmark Problem Based on the IFR-1 Experiment

Metallic nuclear fuels are an active area of research and development for use in advanced reactors. Robust and accurate metallic fuel performance models are necessary for the design, safety analysis, and licensing of such reactors. However, metallic fuel performance models are not as mature as ceramic fuel performance models and require additional development. In order to better gauge the accuracy of existing models, identify high-priority models for development, and quantify the improvements made by future model development, ORNL has developed a benchmark problem based on the IFR-1 experiment. This work gathers all relevant information on the IFR-1 experiment and uses it to develop the benchmark problem. It simulates the problem using the fuel performance code BISON and compare the results to post-irradiation examination data from the experiment. Then a sensitivity analysis is performed on the BISON model to determine the benchmark problem’s sensitivity to uncertainty in the input parameters. The results suggest that BISON’s mechanical models require additional development. Plastic deformation is under-predicted in the cladding and axial swelling is over-predicted in the fuel. These problems may be related. Additionally, there is a bias in the temperature which may be a result of uncertainty in the input parameters rather than an issue with the fuel performance models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Real Time Data Analysis With the ATLAS Trigger at the LHC in Run-2

The trigger selection capabilities of the ATLAS detector have been significantly enhanced for the LHC Run- 2 in order to cope with the higher event rates and with the large number of simultaneous interactions (pile-up) per protonproton bunch crossing. A new hardware system, designed to analyse real time event-topologies at Level-1 came to full use in 2017. A hardware-based track reconstruction system, expected to be used real-time in 2018, is designed to provide track information to the high-level software trigger at its full input rate. The high-level trigger selections are largely relying on offline-like reconstruction techniques, and in some cases multivariate analysis methods. Despite the sudden change in LHC operations during the second half of 2017, which caused an increase in pile-up and therefore also in CPU usage of the trigger algorithms, the set of triggers (so called trigger menu) running online has undergone only minor modifications thanks to the robustness and redundancy of the trigger system, and the use of a levelling luminosity scheme in agreement with LHC and other experiments. Here, this presentation gives a brief yet comprehensive review of the real-time performance of the ATLAS trigger system in 2017. Considerations will be presented on the most relevant parameters of the trigger (efficiency to collect signal and output data rate) as well as details on some aspects of the algorithms which are run real-time on the High Level Trigger CPU farm.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Design and Analysis of Convolutional Neural Network for RF Signal Modulation Classification for In-Orbit Deployment

To effectively transmit data to and from satellites requires a complex and robust RF communication system. Commonly, several different types of signal modulations may be required to maximize satellite efficiency depending on a variety of unexpected channel impairments. We propose a neural network algorithm capable of learning these RF signal modulations using a supervised learning technique designed for low power, high-efficiency in-orbit deployment. The work presented demonstrates a convolutional neural network (CNN) capable of learning and recognizing a set of modulation schemes commonly used to transmit RF information. We are capable of recognizing the modulation scheme from the I and Q data channels directly, with no preprocessing or data conversion required other than breaking the incoming signal into a set of uniform normalized samples. We perform a network design and size analysis, showing that reasonably high accuracy can be obtained using networks with a relatively low number of trainable parameters. Given that a user of a system such as this may wish to receive a signal using a modulation scheme that the network has not previously learned, we demonstrate that transfer learning can learn new modulation schemes by retraining only the fully connected layers in the CNN. Thus, this type of network would excel in outer space deployment using high-efficiency transfer learning hardware. Modulation recognition can be performed through rapid feedforward computation, and the CNN training process is significantly simplified when learning new modulations is required.

CNN↗

From Atom to Engine: Understanding Fundamental Effects of Structure on Combustion Using Tandem Experiment and Computation

The ability to tie structural features to a fuel candidate’s ignition properties provides a path for rational design of advantaged fuels that facilitate higher-efficiency, lower-emitting combustion in engines. The time evolution of a fuel’s radical population directly correlates with global combustion parameters, and because the generated radicals depend on the initial structure and reaction pathways of the fuel molecule, it follows that the molecular structure of a fuel has a direct impact on this cascade, and therefore salient ignition-tied fuel properties. This work aims to understand the impact of structural features in the context of global combustion properties for direct application to emerging engine technologies. Specifically, we highlight the ability to identify potential fuel blendstocks for specific engine strategies by examining and elucidating the radical cascade tied to ignition. A representative subset of alcohols possessing a variety of structural features (e.g. branching, chain length) was selected for study using a robust combination of experimental and computational methods. The approach begins with identifying short-lived radicals and likely abstraction sites upon pyrolysis in a microreactor coupled with photoionization mass spectrometry for direct detection of product species. In tandem, electronic structure calculations were performed to calculate the relevant potential energy surfaces for fuel pyrolysis as well as the subsequent O2 addition-isomerization reactions, which govern ignition behavior. Furthermore, the theoretical energetic and ro-vibrational data are used to calculate rate constants for direct ignition modeling to explore how structural effects impact NTC behavior. The microreactor data combined with electronic structure theory calculations provide a well-informed picture of the reactions relevant to ignition for each fuel in this study, as well as the most likely radicals generated under engine conditions. This thorough fundamental picture of fuel decomposition and ignition chemistry is used to correlate the observed ignition behavior to molecular structural effects.

advanced biofuels↗

Integrating Safety, Security, and Nuclear Operations for Advanced Reactors

The traditional separation between safety, security, and operations teams has created significant barriers to achieving optimal outcomes. When security considerations are introduced late in the design process, they often conflict with already-established architectural, operational, or engineering parameters. Retrofitting security measures can lead to increased costs, schedule delays, and compromises in security effectiveness. For instance, the need to retrofit physical barriers or surveillance systems often results in trade-offs that could have been avoided with earlier input from security professionals. Delayed integration can also affect regulatory processes and result in licensing delays. Security reviews conducted at later stages frequently identify gaps that necessitate significant redesign efforts, impacting not only scope, schedule, and budget, but also adding risk and lowering stakeholder confidence in the project. This paper aims to address these challenges by identifying practical opportunities for integrating security considerations seamlessly with design and operations teams throughout the entire lifecycle of nuclear facilities—from conceptual design to commissioning and beyond. The research emphasizes the value of early and continuous collaboration among stakeholders to ensure that security measures are robust, operationally effective, and cost-efficient. By examining case studies, analyzing past incidents, and leveraging best practices from other high-security industries, this study highlights actionable strategies for bridging the gap between safety, security, and operations teams.

Zineddin, Dr. Z. [ORNL] (ORCID:0009000848740725)↗

Adaptive control in the presence of unmodeled dynamics

Stability and robustness properties of a wide class of adaptive control algorithms in the presence of unmodeled dynamics and output disturbances were investigated. The class of adaptive algorithms considered are those commonly referred to as model reference adaptive control algorithms, self-tuning controllers, and dead beat adaptive controllers, developed for both continuous-time systems and discrete-time systems. A unified analytical approach was developed to examine the class of existing adaptive algorithms. It was discovered that all existing algorithms contain an infinite gain operator in the dynamic system that defines command reference errors and parameter errors; it is argued that such an infinite gain operator appears to be generic to all adaptive algorithms, whether they exhibit explicit or implicit parameter identification. It is concluded that none of the adaptive algorithms considered can be used with confidence in a practical control system design, because instability will set in with a high probability.

Rohrs, C. E.↗

Computational Modeling of Molten Salt Infiltration and Oxidation in Nuclear Graphite

Graphite is utilized as a moderator and reflector in advanced nuclear reactor designs due to its high thermal conductivity, neutron moderation properties, and resistance to radiation damage. However, its longterm performance and reliability are challenged by degradation mechanisms such as molten salt infiltration in molten salt reactors (MSRs) and oxidation in gas-cooled reactors (GCRs). These mechanisms can compromise the structural integrity and operational lifetime of graphite components, necessitating a more detailed assessment of their physical behavior. This report focuses on the development of computational models for molten salt infiltration and oxidation of graphite to aid the design and performance analysis of graphite components. For molten salt infiltration, a computational framework is developed that couples incompressible Navier-Stokes and phase-field model to simulate the penetration of molten salt into graphite?s interconnected pore structure. Initial model verification is performed using two-phase flows in two dimensions, demonstrating the models ability to capture fundamental physical behavior and agree with analytical solution. This framework is then applied to a realistic IG110 nuclear graphite , where a computed tomography extracted pore geometry is used to analyse the infiltration behavior of FLiNaK molten salt. This model provides insights into how the microstructure and other relevant parameters influence the transport pathways of molten salt into graphite, potentially offering a means to rapidly evaluate a graphite grade?s resistance to infiltration. For oxidation, the report details pore-scale mass and heat transport models, describing the diffusion of gases, reaction kinetics, and thermal effects. Additionally, this report highlights inconsistencies in the existing volume-averaged macroscopic model, particularly in upscaling of reaction kinetics and flux terms, and surface to volume transformations. These inconsistencies suggest that current formulations may not accurately capture the experimentally observed graphite oxidation process, highlighting the need for improved model development. This work advances the development of physics-based computational models for graphite degradation, contributing to improved predictive models for next-generation nuclear reactor designs. Future efforts will focus on refining the infiltration model to address non-physical behaviors and enhance its robustness. Additionally, for oxidation, further studies will employ the principles of volume averaging to rigorously derive the upscaled equations, potentially in collaboration with subject matter experts.

Computational Modeling of Molten Salt Infiltration↗

Evaluation of Irradiation Creep Effects in HT9 Cladding for FAST Experiments

The push for advanced reactor fuels for improved reactor safety and efficiency had led to a renewed interest in metallic fuel for nuclear reactor applications. Experimental investigation is necessary to ensure a robust understanding of the thermomechanical properties of new metallic fuel designs. Unfortunately, with the current experimental facilities, thoroughly investigating the responses of metallic fuel burnup would take a prohibitively long time. To alleviate this, the Fission Accelerated Steady State Test (FAST) was developed to accelerate the irradiation testing while simultaneously decreasing the sensitivity to fabrication tolerances by reducing the fuel diameter and scaling the experiment. This method successfully scales the radiation effects on the fuel, but the HT9 cladding is not exposed to prototypic radiation conditions. This raises questions on whether the FAST experiment results are truly indicative of the HT9 cladding performance due to radiation induced creep effects not being appropriately accounted for. Using BISON fuel performance code, the simulated FAST cladding strain is compared to simulated EBR-II cladding strain. This is done through a sensitivity study of input parameters and scaling of neutron fluence on the cladding. This allows a parametric comparison of physical phenomena on the effective difference between cladding strains between FAST and equivalent burnup EBR-II fuel pins. The results show that the irradiation induced deformation (creep or swelling) is insignificant compared to the thermal-mechanical deformation. Therefore, the difference between the FAST experiment cladding and the EBR-II experiment cladding is negligible and comparison of fuel system performance between the two experiments is appropriate.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Approach for Input Uncertainty Propagation and Robust Design in CFD Using Sensitivity Derivatives

An implementation of the approximate statistical moment method for uncertainty propagation and robust optimization for quasi 3-D Euler CFD code is presented. Given uncertainties in statistically independent, random, normally distributed input variables, first- and second-order statistical moment procedures are performed to approximate the uncertainty in the CFD output. Efficient calculation of both first- and second-order sensitivity derivatives is required. In order to assess the validity of the approximations, these moments are compared with statistical moments generated through Monte Carlo simulations. The uncertainties in the CFD input variables are also incorporated into a robust optimization procedure. For this optimization, statistical moments involving first-order sensitivity derivatives appear in the objective function and system constraints. Second-order sensitivity derivatives are used in a gradient-based search to successfully execute a robust optimization. The approximate methods used throughout the analyses are found to be valid when considering robustness about input parameter mean values.

Putko, Michele M.↗

Sensitivity of Optimal Midcourse Correction Scheduling for Robust Cislunar Trajectory Design

A new approach to optimal trajectory design is the determination of optimal trajectories that are robust to initial trajectory dispersions, navigation errors, maneuver execution errors, and environment modeling errors. This paper investigates the sensitivity of cislunar robust optimal trajectory design to launch date, duration of navigation measurement passes, and navigation measurement frequency. For a given cislunar trajectory from translunar injection (TLI) to lunar orbit insertion (LOI), the optimal locations of midcourse corrections, also known as trajectory correction maneuvers (TCM) are determined by minimizing the final 3-σ ∆v subject to a final 3-σ position dispersion constraint for a given launch date, specified measurement pass duration prior to each maneuver, and measurement frequency. Optical navigation (OpNav) is assumed, and OpNav field-of-view (FOV) and lighting constraints are employed. These constraints turn out to be important elements of the problem. The sensitivity of the optimal TCM locations are then investigated by varying the launch date, duration of OpNav measurement passes, and the OpNav measurement frequency, and then re-optimizing the locations of the TCMs. Given the problem parameters provided herein, results show that while the optimal TCM locations with respect to TLI vary greatly from one launch date to another, their locations with respect to LOI are nearly invariant over a 2-month launch window. Results also show that in all cases the optimal location of the last TCM is found to be at the point where the OpNav lunar FOV constraint is first violated. For all other TCMs, OpNav measurement pass duration and measurement frequency can have a moderate to large affect on the optimal TCM locations.

Linear Covariance Analysis↗

Deep Learning Based Approach to Integrate MyShake's Trigger Data with ShakeAlert for Faster and Robust EEW Alerts

Earthquake Early Warning (EEW) systems depend on the dense seismic network to make fast and accurate earthquake detections to issue useful early alerts. The recent development and implementation of the ShakeAlert system is relying on the west coast seismic network to detect and estimate the earthquake parameters in real-time. While working on making improvements on the current system, one potential approach is to include more seismic recordings from various sources, such as the low-cost smartphone seismic network, such as the MyShake network, or the Google Android smartphones, since these smartphone seismic networks have much more portable stations (phones) to potentially provide useful data to the system. This report summarizes the initial exploration of using a deep learning approach to combine the data from both traditional seismic stations and the smartphone data in California. Due to the heterogeneity nature of the data, we aggregate data from the phones (using simulation data), as well as that from traditional seismic stations to grid cells. By generating real-time triggering ratio grid cell maps, the designed deep learning algorithm can process the data from multiple sources and detect the earthquake faster than only using that from a traditional seismic network.

58 GEOSCIENCES↗

CFD Analysis of the 24-inch JIRAD Hybrid Rocket Motor

A series of multispecies, multiphase computational fluid dynamics (CFD) analyses of the 24-inch diameter joint government industry industrial research and development (JIRAD) hybrid rocket motor is described. The 24-inch JIRAD hybrid motor operates by injection of liquid oxygen (LOX) into a vaporization plenum chamber upstream of ports in the hydroxyl-terminated polybutadiene (HTPB) solid fuel. The injector spray pattern had a strong influence on combustion stability of the JIRAD motor so a CFD study was initiated to define the injector end flow field under different oxidizer spray patterns and operating conditions. By using CFD to gain a clear picture of the flow field and temperature distribution within the JIRAD motor, it is hoped that the fundamental mechanisms of hybrid combustion instability may be identified and then suppressed by simple alterations to the oxidizer injection parameters such as injection angle and velocity. The simulations in this study were carried out using the General Algorithm for Analysis of Combustion SYstems (GALACSY) multiphase combustion codes. GALACSY consists of a comprehensive set of droplet dynamic submodels (atomization, evaporation, etc.) and a computationally efficient hydrocarbon chemistry package built around a robust Navier-Stokes solver optimized for low Mach number flows. Lagrangian tracking of dispersed particles describes a closely coupled spray phase. The CFD cases described in this paper represent various levels of simplification of the problem. They include: (A) gaseous oxygen with combusting fuel vapor blowing off the walls at various oxidizer injection angles and velocities, (B) gaseous oxygen with combusting fuel vapor blowing off the walls, and (C) liquid oxygen with combusting fuel vapor blowing off the walls. The study used an axisymmetric model and the results indicate that the injector design significantly effects the flow field in the injector end of the motor. Markedly different recirculation patterns are observed in the vaporization chamber as the oxygen velocity and/or spray pattern is varied. The ability of these recirculation patterns to stabilize the diffusion flame above the surface of the solid fuel gives a plausible explanation for the experimentally determined combustion stability characteristics of the JIRAD motor, and suggests how combustion stability can be assured by modifications to the injector design.

Liang, Pak-Yan↗