Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reaction Control Propulsion System”

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.

12 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A wide range experimental and kinetic modeling study of the oxidation of 2,3-dimethyl-2-butene: Part 1

2,3-Dimethyl-2-butene (TME) is a potential fuel additive with high research octane number (RON) and octane sensitivity (S), which can improve internal combustion engine performance and efficiency. However, the combustion characteristics of TME have not been comprehensively investigated. Thus, it is essential to study the combustion characteristics of TME and construct a detailed chemical kinetic model to describe its combustion. In this paper, two high-pressure shock tubes and a constant-volume reactor are used to measure ignition delay times and laminar flame speeds of TME oxidation. Further, the ignition delay times were measured at equivalence ratios of 0.5, 1.0, and 2.0 in “air”, at pressures of 5 and 10 bar, in the temperature range of 950 – 1500 K. Flame speeds of the TME/ “air” mixtures were measured at atmospheric pressure, at a temperature of 325 K, for equivalence ratios ranging from 0.78 to 1.31. Two detailed kinetic mechanisms were constructed independently using different methodologies; the KAUST TME mechanism was constructed based on NUIGMech1.1, and the MIT TME mechanism was built using the Reaction Mechanism Generator (RMG). Both mechanisms were used to simulate the experimental results using Chemkin Pro. In the present work, reaction flux and sensitivity analyses were performed using the KAUST mechanism to determine the critical reactions controlling TME oxidation at the conditions studied.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Assessment of the impact of reactor residence time distribution on non-equilibrium product selectivity of polypropylene pyrolysis using reactive molecular dynamics simulations

Mass residence time distribution (RTD) is considered to be an important factor controlling the product selectivity in the pyrolysis of biomass and plastic wastes along with the pyrolysis chemistry. However, due to the complex pyrolysis chemistry of biomass and plastic waste, the coupling between the reaction chemistry, RTD, and product selectivity is challenging to understand. Here, we introduce a reaction molecular dynamics-based method to examine pyrolysis chemistry and species timescales to assess the impact of RTD on product selectivity and yield. To validate this method, reactive molecular dynamics simulations were conducted for polypropylene pyrolysis and its non-equilibrium product selectivity using 6 different RTDs. We find that the RTD and the reaction chemistry control the peak non-equilibrium product concentrations. The peak monomer (C 3 H 6 ) concentration during pyrolysis can be increased by up to 25 % by using a narrow RTD in the case of polypropylene pyrolysis. We also find that product selectivity is strongly affected by the average residence time and RTD. This coupling between the reaction chemistry, RTD, and product selectivity highlights the need to understand detailed reaction chemistry to control RTD and optimize non-equilibrium product selectivity during polymer and biomass pyrolysis. The present method provides a new way to design RTD for reactors to reach maximized product selectivity of plastic waste and biomass.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

A comprehensive experimental and kinetic modeling study of di-isobutylene isomers: Part 1

We report di-isobutylene has received significant attention as a promising fuel blendstock, as it can be synthesized via biological routes and is a short-listed molecule from the Co-Optima initiative. Di-isobutylene is also popularly used as an alkene representative in multi-component surrogate models for engine studies of gasoline fuels. However, there is limited experimental data available in the literature for neat di-isobutylene under engine-like conditions. Hence, most existing di-isobutylene models have not been extensively validated, particularly at lower temperatures (< 1000 K). Most gasoline surrogate models include the di-isobutylene sub-mechanism published by Metcalfe et al. with little or no modification. The current study is undertaken to develop a detailed kinetic model for di-isobutylene and validate the model using a wide range of relevant experimental data. Part 1 of this study exclusively focuses on the low- to intermediate temperature kinetics of di-isobutylene. An upcoming part 2 discusses the high-temperature model development and validation of the relevant experimental targets. Ignition delay time measurements for the di-isobutylene isomers were performed at pressures ranging from 15 - 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 diluted in air and in the temperature range 650 - 900 K using two independent rapid compression machine facilities. In addition, measurements of species identified during the oxidation of these isomers were performed in a jet-stirred reactor and in a rapid compression machine. A detailed kinetic model for the di-isobutylene isomers is developed to capture the wide range of new experimental targets. For the first time, a comprehensive low-temperature chemistry submodel is included. The differences in the important reaction pathways for the accurate prediction of the oxidation of the two DIB isomers are compared using reaction path analysis. In conclusion, the most sensitive reactions controlling the ignition delay times of the DIB isomers under the pressure and temperature conditions necessary for autoignition in engines are identified.

09 BIOMASS FUELS↗
DOE OSTI2022

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

NO x Control and Measurement Technology for Heavy-Duty Diesel Engines (CRADA Final Report)

The Catalyst CRADA focused on lean-NO x trap (LNT) and selective catalytic reduction (SCR) automotive emission control catalyst technology. It particularly had the goal of developing and applying advanced diagnostics capable of elucidating the network and sequence of catalyst reactions, as well as how they are affected by real-world operation and aging. Another overarching goal and action of this CRADA section was to develop methods to quantify and monitor catalyst states and support the development of full and reduced models for design and in-use control applications. William Partridge was the ORNL principal investigator (PI) throughout the Catalyst CRADA; Neal Currier was the Cummins PI through circa 2017, and Saurabh Joshi was the subsequent Cummins PI.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Optimization of WAAM Process to Produce AUSC Components with Increased Service Life

Additive manufacturing has the potential to revolutionize industrial hardware and unlock efficiency gains through the fabrication of geometries and architectures not possible by conventional processing. Wire Arc Additive manufacturing (WAAM) process is a class of directed energy deposition process enabling higher build rate and allowing custom wire feedstock allowing spatial variation of microstructure. The larger spot size and lower speed creates a larger melt pool, reducing residual stress and often time creating directional/columnar microstructure for Ni-superalloy. The use of flexible platform provides freedom in deposition strategy, which can accommodate complex substrates, including feature addition onto existing structures, non-flat layers, and repair methods. However, the certification of the final component needs to match the strength requirement. Hence, the quality of the component produced is stringently monitored to avoid buildup of residual stress, cracks, porosity and to reduce detrimental segregated phases commonly observed during alloy solidification. Simulation plays a huge role in predicting the melt pool dimension and can be used to optimize the process parameter. Similar development is also required to perform physics-based modeling of microstructural development in WAAM that can predict the microstructural features during solidification and can be used to optimize the process more effectively. To move toward this goal, Raytheon Technologies Research Center together with Siemens worked to create a set of computational tools to control the process parameters, enable on-line measurements and acquisition with feedback to the optimized process parameters, and eventually track material evolution through each step of the additive process. Computational fluid dynamics is used for accurate prediction and calibration of the thermal field during WAAM process. and phase field models for microstructure evolution as a function of processing parameters to establish a connection between additive parameters and the final microstructure. Here we report cellular automata (CA) model development to predict the dendritic microstructure evolution with surface and bulk nuclei for single track and multiple layers. The CA model was developed to account for secondary element addition and predict segregation, local melting, and latent heat release as well as prediction of Euler angles from orientation information and validated against experiments. This framework was utilized to tailor spatially-varying composition in a part by appropriately controlling the microstructure evolution during the additive process. Functionally graded Haynes 282 alloy with high Cr content at the surface was tested for oxidation and mechanical properties. An advanced physics-based reaction-diffusion model predicting the simultaneous creation of chromium oxide and alumina is developed and validated to extend life expectancy of the WAAM manufactured high temperature part. A machine-learning data-driven framework establishing the process-structure relationship from a dataset of real microstructure images and corresponding process history data has been developed and implemented in the NX Siemens design system. The digital twin configuration along with the tool path generation enabled prediction of WAAM component buildup time and the techno-economic analysis provided a favorable option for all 4 cases with 15-40% cost reduction.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Detailed chemistry modeling of rotating detonations with dilute n -heptane sprays and preheated air

Utilization of liquid fuels is crucial to enabling commercialization of rotating detonation engines (RDE) in the near future. In this study, Eulerian-Lagrangian simulations are conducted for rotating detonative combustion with dilute n-heptane sprays and preheated air. Two-dimensional flattened configuration is used and a skeletal chemical mechanism with 44 species and 112 elementary reactions for n-heptane combustion is adopted. The flow structure, droplet distribution, and thermochemical parameters in the refill zone are first analyzed. It is shown that the mixture in the refill zone is heterogeneous, including evaporating droplets, vapor, and air. When the total temperature is below 950 K, the average equivalence ratio increases with the total temperature. When it is higher than 950 K, the average equivalence ratio is almost constant. Subsequently, the chemical explosive mode analysis is applied to identify the controlling reactions and dominant combustion modes in the fuel refill zone and reaction fronts. Results demonstrate that the initiation reaction (R104: n-C 7 H 16 + O 2 2-C 7 H 15 + HO 2 ) and low-temperature reaction (R107: RO 2 R'O 2 H) are dominant in the upstream and downstream of the refill zone, respectively. The intermediate species from low-temperature chemistry, R'O 2 H, is found to be important for the chemical explosive mode in the undetonated mixture. The influence of species diffusion and dispersed droplets is further analyzed. Results show that vapor autoignition facilitated by droplet evaporation occurs in the refill zone. Finally, the effects of the air total temperature on the detonation propagation speed and RDE propulsion performance are investigated. Further, it is found that the detonation propagation speed and specific impulse increase with air total temperature. The total pressure ratio first increases and then decreases as the air total temperature increases. Moreover, when the total temperature of the preheated air is above 1,300 K, the effects of low-temperature chemistry are negligible.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Full-Loop Simulation of the Combustion of Biomass in a Circulating Fluidized Bed Combustor

Interest in circulating fluidized bed (CFB) boilers as a power generation technology has skyrocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing, which results in higher combustion efficiency, and the ability to use lower rank fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions, while SO 2 emissions can be conveniently controlled through the addition of Ca-based sulfur sorbents within the combustor. This report summarizes the development of a comprehensive computational fluid dynamics (CFD) model for biomass combustion in a 50 kW th CFB combustor employing the multiphase particle-in-cell (PIC) approach in the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite at the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL). The experimental combustor was designed, built, and operated at CanmetENERGY, Natural Resources Canada in Ottawa, Canada. This work is the culmination of a multi-year collaboration between NETL and CanmetENERGY. The hydrodynamics in the full loop are validated against cold-flow experiments conducted at CanmetENERGY. To mimic the dense packing in the standpipe in the absence of a collisional model in the MFiX-PIC approach, the ram valve in the return leg is partially closed to obstruct the flow. A combustion reaction scheme was developed that considers the volatiles released by the fuel as a lumped species that fully converts to CO 2 and H 2 O. The so-called enhanced simplified combustion scheme produces an excellent match with experimental data in terms of the species concentrations of CO 2 and O 2 at the outlet and the freeboard, demonstrating the utility of the comprehensive model for a range of flow conditions and operating temperatures. Furthermore, the model can be readily adapted to simulate combustion of a variety of biomass fuels as well as co-fired systems under air and oxy-fuel conditions.

09 BIOMASS FUELS↗
DOE OSTI2024

When and Where Lithium Plating Occurs, Its Correlation with Microstructure Heterogeneity, and the Mechanisms That Initiate and Self-Regulate Electrochemical Heterogeneity (A02-0444)

A microstructure scale electrochemical LIB model was used to investigate lithium plating onset, material non-uniform utilization, and in-plane heterogeneities for an NMC-graphite full cell. Model predicts active material particle surface roughness and size distribution (respectively, non-uniform curvature within and between particles) initiate in-plane heterogeneity, and that particle size heterogeneity at the separator interface controls the lithium plating preferential deposition ("Where"). These in-plane heterogeneities are then exacerbated by through-plane heterogeneities induced at fast charge as electrolyte depletion occurs and concentrates intercalation reaction near the anode-separator interface. Also, magnitude and occurrence of lithium plating is controlled by effective, or macroscale, microstructure parameters ("When"). As local states of charge start to diverge between nearby active material regions, overpotential differences induced by OCP difference kick in and contribute to reduce these SOC local heterogeneities. However, for staged materials such as graphite, with OCP profile alternating between plateaus and varying regions, this balancing mechanism is, respectively, inactive and active. This leads to a dynamic, non-monotonic, in-plane heterogeneity time evolution for state of charge and Faraday current density, for which their respective in-plane heterogeneity magnitude alternates. Such behavior has been modeled both for the whole electrode at the microstructure scale and at the particle scale. In-plane heterogeneities are usually considered to be detrimental, as they result in material non-uniform utilization (i.e., under and over stressed regions) and earlier degradations. However, this work provides a more granular approach as it discriminates between a harmful in-plane heterogeneity (non-uniform curvature) that triggers SOC in-plane heterogeneity, and a beneficial in-plane heterogeneity (Faraday current density) that contributes to reduce SOC in-plane heterogeneity. This work comprehensively explains the mechanisms that initiate, exacerbate, and regulate heterogeneity at the microstructure scale, while providing some design suggestions to reduce both in-plane and through-plane heterogeneities, as summarized in the graphical abstract.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Asynchronicity in opposed-piston RCMs: Does it matter?

Rapid Compression Machines (RCMs) are widely utilized to study combustion phenomena at engine-relevant conditions, and significant efforts are typically made to create a quiescent environment, particularly for investigations of autoignition chemistry. Opposed-piston configurations can be advantageous due to shorter compression times and reduced surface area to volume ratios. Each side must be actuated simultaneously, but this can be challenging in practice. These devices, like most RCMs, utilize hydraulics for actuation, speed control and arrestation of the piston at the end of the stroke; there is no mechanical control or linkage of the two piston trajectories. To quantify the magnitudes and effects of piston asynchronous behavior, this work employs both detailed experimental measurements and, for the first time, high-fidelity, Direct Numerical Simulation (DNS). The boundary conditions are carefully considered applying insight from high-resolution linear variable differential transformer (LVDT) measurements of the piston trajectory and a zero-dimensional kinematics model of the piston-shaft assembly. Sufficient resolution in the piston crevice region is used. The complicated fluid dynamical behavior that can evolve during piston compression and the ensuing delay processes due to offset timings from t offset = 0-10 ms is elucidated. It is found that near t offset = 6 ms and beyond, the boundary layer on the face of the first-seating piston can be sufficiently perturbed, due initially to reemergence of gas from the crevice of the firstseating piston, so that the adiabatic core can become degraded at long ignition delay times. Substantial mixing of colder gas into the interior of the reaction chamber can alter the measurements, similar to effects previously observed for improper piston crevice configuration. In conclusion, experimental techniques to mitigate asynchronous behavior are discussed and demonstrated.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Three-Dimensional CFD Investigation of Pre-Spark Heat Release in a Boosted SI Engine

Low-temperature heat release (LTHR) in spark-ignited internal combustion engines is a critical step toward the occurrence of auto-ignition, which can lead to an undesirable phenomenon known as engine knock. As such, correct predictions of LTHR are of utmost importance to improve the understanding of knock and enable techniques aimed at controlling it. While LTHR is typically obscured by the deflagration following the spark ignition, extremely late ignition timings can lead to LTHR occurrence prior to the spark, i.e., pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection SI engine was numerically investigated using three-dimensional computational fluid dynamics (CFD). A hybrid approach was used, based on the G-equation model for representing the turbulent flame front and the multi-zone well-stirred reactor model for tracking the chemical reactions within the unburnt region. A recently developed best practice was also employed which keeps the well-stirred reactor model active throughout the entire simulation. This allowed for correct predictions of the previous cycle trapped residuals which have a considerable effect on the onset of LTHR. Multi-cycle simulations were conducted using Co-Optima alkylate and E30 fuels. The predicted in-cylinder pressure and heat release rate agreed well with the experimental data and served as validation for the CFD model. Following the initial validation, the dynamics of PSHR was discussed and a series of parameters of interest were assessed. First, the effect of exhaust valve temperature was investigated, qualitatively highlighting the importance of boundary conditions uncertainty. Further analyses were carried out on the effects of fuel properties, including laminar flame speed (LFS) and heat of vaporization (HOV). The results indicate that PSHR phasing is slightly advanced with lower LFS as more trapped unburnt fuel is made available for the next cycle. A similar trend was observed with lower HOV as less intense spray cooling led to higher mixture temperatures, i.e., higher mixture reactivity. Finally, a comparison of Co-Optima alkylate and E30 fuels was made using the pressure-temperature trajectory framework. It was shown that the differences between Co-Optima alkylate and E30 in PSHR tendency are correlated with both HOV and chemical effects.

33 ADVANCED PROPULSION SYSTEMS↗