Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic modeling”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Role of copper in nitrous oxide accumulation in wetland soils from Illinois and South Carolina and stream sediments from Tennessee

Denitrification is microbially-mediated through enzymes containing metal cofactors. Laboratory studies of pure cultures have highlighted that the availability of copper (Cu), required for the multicopper enzyme nitrous oxide reductase, can limit nitrous oxide (N2O) reduction. However, in natural aquatic systems, such as wetlands and hyporheic zones in stream beds, the role of Cu in controlling denitrification remains incompletely understood. In this study, we collected soils and sediments from three natural environments -- riparian wetlands, marsh wetlands, and a stream -- to investigate their nitrogen species transformation activity at background Cu levels and different supplemented Cu loadings. All of the systems contained solid-phase associated Cu below or around geological levels (40–280 nmol g-1) and exhibited low dissolved Cu (3–50 nM), which made them appropriate sites for evaluating the effect of limited Cu availability on denitrification.The dataset contains the variation in the nutrient (nitrate (NO3-), nitrite (NO2-), and ammonium (NH4+)) and N2O concentrations during the incubation experiments to evaluate the effect of Cu. We have also reported the variation in metal concentrations (Cu, iron(Fe) and manganese (Mn)) and dissolved organic carbon (DOC) concentrations during the incubations. The parameters obtained using a kinetic model to quantitatively report the effect of Cu on nitrogen species conversion in natural aquatic systems are also added in the dataset. The dataset also contains the labile concentration of Cu estimated using a speciation model in MINTEQ. All data are provided in text-based CSV format with header sections indicating the data contained in each file and the corresponding units.Our study suggests that high concentrations of N2O accumulated in all microcosms lacking Cu amendment except for one stream sediment sample. With Cu added to provide dissolved concentrations at trace levels (10–300 nM), the reduction rate of N2O to N2 in the wetland soils and stream sediments was enhanced. A kinetic model could account for the trends in nitrogen species by combining the reactions for microbial reduction of NO3- to NO2-/N2O/N2 and abiotic reduction of NO2 to nitrogen (N2_. The model revealed that the rate of N2O to N2 conversion increased significantly in the presence of Cu. For riparian wetland soils and stream sediments, the kinetic model also suggested that overall denitrification is driven by abiotic reduction of NO2- in the presence of inorganic electron donors. This study demonstrated that natural aquatic systems containing Cu at concentrations less than or equal to crustal abundances may display incomplete reduction of N2O to N2 that would cause N2O accumulation and release to the atmosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

An experimental and modeling study on autoignition of 2-phenylethanol and its blends with n-heptane

2-Phenylethanol (2-PE) is an aromatic alcohol with high research octane number, high octane sensitivity, and a potential to be produced using biomass. Considering that 2-PE can be used as a fuel additive for boosting the anti-knocking quality of gasoline in spark-ignition engines and as the low reactivity fuel or fuel component in dual-fuel reactivity controlled compression ignition (RCCI) engines, it is of fundamental and practical interest to understand the autoignition chemistry of 2-PE, especially at low-to-intermediate temperatures (<1000 K). Based upon the experimental ignition delay time (IDT) results of neat 2-PE obtained from our previous rapid compression machine (RCM) investigation and the literature shock tube study, a detailed chemical kinetic model of 2-PE is developed herein, covering low-to-high temperature regimes. Besides, RCM experiments using binary fuel blends of 2-PE and n-heptane (nC7) are conducted in this work to investigate the nC7/2-PE blending effects, as they represent a dual-fuel system for RCCI operations. Furthermore, the newly developed 2-PE model is merged with a well-validated nC7 kinetic model to generate the current nC7/2-PE binary blend model. Overall, the consolidated model reasonably predicts the experimental IDT data of neat 2-PE and nC7/2-PE blends, as well as captures the experimental effects of pressure, equivalence ratio, and blending ratio on autoignition. Finally, model-based chemical kinetic analyses are carried out to understand and identify the controlling chemistry accounting for the observed blending effects in RCM experiments. In conclusion, the analyses reveal that nC7 enhances 2-PE autoignition via providing extra $\dot{O}$H radicals to the shared radical pool, while the diminished nC7 promoting effect on 2-PE autoignition with increasing temperature is due to the negative temperature coefficient characteristics of nC7.

2-phenlyethanol↗
DOE OSTI2020

Development of a high temperature pyrolysis mechanism for cyclopentanone, a potential biofuel derived from biomass

Cyclopentanone (CPO) is a promising biofuel for spark-ignition engines due to its ring strain and high auto-ignition resistance. The Department of Energy has identified CPO as a potential biofuel candidate because of its unique chemical properties. Understanding CPO decomposition is crucial to building a high-temperature combustion model. In this study, we present a comprehensive kinetic model for high-temperature pyrolysis of CPO with verified results from high-pressure shock tube (HPST) measurements. The time histories of carbon monoxide (CO), ethylene (C2H4), and CPO absorbances over the temperature range of 1156-1416 K and 8.53-10.06 atm were measured. A corresponding kinetic model was generated using the Reaction Mechanism Generator (RMG) with the parameters of dominant pathways from either previous studies or quantum calculation within current work. The final model containing 821 species and 79,859 reactions exhibited a good agreement with the experimental results. In this study, the C2H4/CO ratio was used as an essential factor to validate the model and to prove that radical-involved bimolecular pathways were as significant as unimolecular decomposition of CPO. The rate of production (ROP) analysis showed that H radicals play major roles in both primary reactions and product pyrolysis. The insights from this work can be used to generate a better CPO combustion model and help evaluate CPO as a biofuel. Such understanding can also be helpful when proposing pyrolysis or combustion models for other biofuels.

cyclopentanone, reaction mechanism, pyrolysis, sho↗
DOE OSTI2025

Experimental Measurements of Soot and Soot Precursors of Sustainable Aviation Fuels

Sustainable Aviation Fuels (SAFs) are promising for decarbonizing the challenging-to-electrify aerospace sector and reducing yearly greenhouse gas emissions by thousands of tons. Developing predictive detailed chemical kinetic models for SAF (and traditional) fuel mixtures necessitates a robust experimental database to elucidate the intricate pyrolysis and oxidation chemistry of SAF fuel blend components. Two critical components of SAF surrogate mixtures are iso-dodecane isomers and trans-decalin. This report characterizes the structure of three slightly sooting, laminar, non-premixed, nitrogen-diluted Planar Mixing Layer Flames (PMLFs) fueled by either pure ethylene, ethylene doped with 1500ppm of either 2,2,4,6,6-pentamethyl-heptane or trans-decalin. The flames have the same stoichiometric mixture fraction and total mole fraction of fuels in the fuel stream (X F,F = X C2H4,F + X Dopant,F = 0.26), resulting in approximately the same maximum temperature (T max ≈ 1800 K) to highlight the effect of doping on soot and its gaseous precursors. The PMLF is unconfined and established between adjacent planar jets of the fuel and oxidizer streams surrounded by annular shielding nitrogen. The flow is stabilized by exhausting the hot buoyant flame products through a slot in a downstream plate onto which the cold flow impinges. Importantly, any horizontal cross-section of the PMLF has a self-similar thermochemical structure of a reactive boundary layer that grows at increasing Height Above the Burner (HAB). Therefore, the horizontal cross-sections of a PMLF can be modeled as a One-Dimensional Counterflow Flame (1D-CF) with very low strain rates which are unachievable in traditional CFs and yield a structure that is several millimeters thick. In brief, the PMLF is ideal for performing spatially resolved measurements with minimal effects of the sampling-induced perturbations, to develop and validate detailed chemical kinetic models on time scales of a few tens of milliseconds that are relevant in aviation applications. The horizontal profiles of C0-C18 gaseous species at HAB = 50 mm and soot volume fraction at 25 and 50 mm are measured using capillary sampling followed by GC-MS analyses and Laser Induced Emission Spectroscopy (LIES), respectively. Additionally, the Elastic Laser Light Scattering (E-LLS) coefficient is measured at HAB = 50 mm and 25 mm to determine the strain rates of the 1D-CFs equivalent to the characterized PMLF horizontal cross sections and the profile of the E-LLS equivalent diameter of soot. The comparison of the results in the base PMLF fueled only by ethylene with those in either doped PMLF in which 1500 ppm of ethylene are replaced with 2,2,4,6,6-pentamethyl-heptane or trans-decalin, reveals that the doping increases the mole fraction of several Polycyclic Aromatic Hydrocarbons (PAHs) at both HABs by a factor of approximately 2.5, and the soot volume fraction at both HABs by factors of ~1.5 for 2,2,4,6,6-pentamethyl-heptane and ~1.25 for trans-decalin doping, respectively. Conversely, doping causes a near doubling of the soot equivalent diameter only at HAB = 50 mm, without significant effects at 25 mm. The experimental results partially validate the tested state-of-the-art detailed chemical kinetic model and also point toward the further improvement of its predictive capabilities.

10 SYNTHETIC FUELS↗
DOE OSTI2022

A self-consistent hybrid model of kinetic striations in low-current argon discharges

A self-consistent hybrid model of standing and moving striations was developed for low-current DC discharges in noble gases. We introduced the concept of surface diffusion in phase space ( r , u ) (where u denotes the electron kinetic energy) described by a tensor diffusion in the nonlocal Fokker–Planck kinetic equation for electrons in the collisional plasma. Electrons diffuse along surfaces of constant total energy ε = u - eφ ( r ) between energy jumps in inelastic collisions with atoms. Numerical solutions of the 1d1u kinetic equation for electrons were obtained by two methods and coupled to ion transport and Poisson solver. We studied the dynamics of striation formation in Townsend and glow discharges in argon gas at low discharge currents using a two-level excitation-ionization model and a ‘full-chemistry’ model, which includes stepwise and Penning ionization. Standing striations appeared in Townsend and glow discharges at low currents, and moving striations were obtained for the discharge currents exceeding a critical value. These waves originate at the anode and propagate towards the cathode. We have seen two types of moving striations with the two-level and full-chemistry models, which resemble the s and p striations previously observed in the experiments. Simulations indicate that processes in the anode region could control moving striations in the positive column plasma. The developed model helps clarify the nature of standing and moving striations in DC discharges of noble gases at low discharge currents and low gas pressures.

Physics↗
DOE OSTI2026

Data and code from: Selective oxidation of 3-hydroxypropionic acid to malonic acid over Pd/C: Mechanistic and kinetic study

Malonic acid (MA) is a high-value dicarboxylic acid with strong industrial demand, yet its current production heavily relies on petrochemical feedstocks. Here, we report the first systematic study for the sustainable production of MA via oxidation of 3-hydroxypropionic acid (3-HP) with a Pd/Carbon catalyst. The effects of oxidant type (O2 and H2O2), pH, and temperature on the reaction chemistry were comprehensively evaluated. Guided by experimental observation and DFT-calculated thermochemical energetics, reaction networks for 3-HP oxidation with both oxidants are proposed and validated through kinetic modeling. MA was identified as the primary oxidation product, while further oxidation yields acetic acid and oxalic acid. The kinetic model validated the network, displaying excellent agreement (R2 > 0.95). Kinetic observations also enabled a direct comparison between O2 and H2O2 and revealed their distinct behaviors. The model was further developed into a temperature-time map, providing insight into conditions that maximize MA production. Malonic acid selectivity of 56.9 % and yield of 50.5 % were achieved at 3 bar oxygen, equimolar NaOH:3-HP ratio, and 50 °C.

3-Hydroxypropionic Acid↗
NASA NTRS2010

Multi-Scale Modeling of Global of Magnetospheric Dynamics

To understand the role of magnetic reconnection in global evolution of magnetosphere and to place spacecraft observations into global context it is essential to perform global simulations with physically motivated model of dissipation that is capable to reproduce reconnection rates predicted by kinetic models. In our efforts to bridge the gap between small scale kinetic modeling and global simulations we introduced an approach that allows to quantify the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions near reconnection sites. We utilized the high resolution global MHD code BATSRUS and incorporate primary mechanism controlling the dissipation in the vicinity of reconnection sites in terms of kinetic corrections to induction and energy equations. One of the key elements of the multiscale modeling of magnetic reconnection is identification of reconnection sites and boundaries of surrounding diffusion regions where non-MHD corrections are required. Reconnection site search in the equatorial plane implemented in our previous studies is extended to cusp and magnetopause reconnection, as well as for magnetotail reconnection in realistic asymmetric configurations. The role of feedback between the non-ideal effects in diffusion regions and global magnetosphere structure and dynamics will be discussed.

Kuznetsova, M. M.↗
DOE OSTI2026

Selective oxidation of 3-hydroxypropionic acid to malonic acid over Pd/C: Mechanistic and kinetic study

Malonic acid (MA) is a high-value dicarboxylic acid with strong industrial demand, yet its current production heavily relies on petrochemical feedstocks. Here, we report the first systematic study for sustainable production of MA via oxidation of 3-hydroxypropionic acid (3-HP) with a Pd/Carbon catalyst. The effects of oxidant type (O 2 and H 2 O 2 ), pH, and temperature on the reaction chemistry were comprehensively evaluated. Guided by experimental observation and DFT-calculated thermochemical energetics, reaction networks for 3-HP oxidation with both oxidants are proposed and validated through kinetic modeling. MA was identified as the primary oxidation product, while further oxidation yields acetic acid and oxalic acid. The kinetic model validated the network, displaying excellent agreement (R 2 > 0.95). Kinetic observations also enabled a direct comparison between O 2 and H 2 O 2 and revealed their distinct behaviors. The model was further developed into a temperature-time map, providing insight into conditions that maximize MA production. Malonic acid selectivity of 56.9 % and yield of 50.5 % were achieved at 3 bar oxygen, equimolar NaOH:3-HP ratio, and 50 °C.

3-hydroxypropionic acid↗
DOE OSTI2024

Toward UV models of kinetic mixing and portal matter. VI. A more complex dark matter sector?

Portal matter (PM), having both Standard Model (SM) and dark sector charges, can induce kinetic mixing between the 𝑈⁢(1) 𝐷 dark photon and the SM gauge fields at the 1-loop level offering an attractive mechanism by which light (≲1 GeV) thermal dark matter (DM) can interact with visible matter and obtain its observed relic density. In doing so, if the DM is fermionic, the CMB and other astrophysical observations inform us that it must be Majorana/pseudo-Dirac in nature to avoid velocity/temperature-independent 𝑠-wave annihilation to SM final states. How does this idea fit into a more UV-complete picture also including the SM interactions? There are some reasons to believe that at least a first step along this path may not lie too far away in energy due to the renormalization group equations running of the dark gauge coupling, which for a significant range of parameters, becomes nonperturbative at/before the ∼10’s of TeV energy range. This implies that 𝑈⁢(1) 𝐷 must become embedded in an asymptotically free, non-Abelian group, 𝐺 𝐷 , before this can occur. The breaking of this larger group then produces the masses for the PM and the additional gauge fields associated with 𝐺 𝐷 then can lead to new interactions between the SM and the dark sector. Following several bottom-up approaches, we have examined a set of distinctive and testable phenomenological features associated with this general setup, based upon a number of simplifying assumptions. Clearly, it behooves us to explore the impact of these specific assumptions on these predictions for the array of possible experimental tests of this class of models. In most past analyses it has been assumed that DM is a vectorlike, complex singlet under the group 𝐺 𝐷 . If this assumption is relaxed, the dark sector must be augmented by additional fermion(s) and the associated scalar fields needed to break the gauge symmetries while generating the needed Majorana-like mass terms for the DM. In this paper, we analyze the simplest extension of this kind wherein the DM lies in a vectorlike doublet of 𝐺 𝐷 , which we take to have the structure 𝑆⁢𝑈⁢(2) 𝐼 ×𝑈⁢(1) 𝑌 𝐼 as in earlier work, leading to new phenomenological implications. We find, for example, that given the current LHC search constraints on the masses of heavy gauge bosons, the production of these new dark states with large rates is unlikely to occur at colliders unless they are produced singly in 𝑔⁡𝑔 fusion or their pair production cross sections are resonantly enhanced. Here, we also find that an additional mechanism arises to generate hierarchal neutrino masses in such a setup.

Extensions of Higgs sector↗
NASA NTRS1993

Kinetic theory model predictions compared with low-thrust axisymmetric nozzle plume data

A system of nonlinear integral equations equivalent to the steady-state Krook kinetic equation was used to model the flow from a low-thrust axisymmetric nozzle. The mathematical model was used to numerically calculate the number density, temperature, and velocity of a simple gas as it expands into a near vacuum. With these quantities the gas pressure and flow directions of the gas near the exit plane were calculated and compared with experimental values for a low-thrust nozzle of the same geometry and mass flow rate.

Riley, B. R.↗
DOE OSTI2020

Final Report for Tech-X Corporations’s contribution to Enhancing Understanding of High Energy Density Plasmas from Wire Array and Solid Liner Implosions Using Fluid Modeling with Kinetic Closures

Recent results from experiments and simulation of magnetically driven pulsed power liners have explored the role of the early-time electrothermal instability in the evolution of the magneto-Rayleigh-Taylor instability. Our focus will be on understanding the development of such instabilities and the potential stabilization mechanisms via electron thermal conduction, viscosity, and large magnetic fields which we expect could play a significant role in supporting the success of the MagLIF program. Experiments have shown that studies of high-energy density plasmas from wire-array implosions require physics modeling that goes well beyond simple models such as ideal magnetohydrodynamics. The goal of this work is to provide increased understanding of these experiments by employing simulations with a multi-fluid extended-MHD model which uses kinetic closures for thermal conduction, resistivity and viscosity. We will use codes easily available to the wider research community, including university students, with a secondary goal of providing the community with well-benchmarked tools capable of advanced modeling of high-energy-density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Cellular repair/misrepair track model

A repair/misrepair cell kinetics model is superimposed onto the track structure model of Katz to provide for a repair mechanism. The model is tested on the repair-dependent data of Yang et al. and provides an adequate description of that data. The misrepair rate determines the maximum relative biological effectiveness (RBE), but similar results could arise from indirect X-ray lethality not include in the present model.

Wilson, John W.↗
DOE OSTI2025

Corrosion of 316 L stainless steel under the natural circulation of molten NaCl-MgCl 2 salt

The corrosion behavior of 316 L stainless steel (SS) was studied via the natural circulation of molten eutectic NaCl-MgCl 2 salt through a microloop. The post-corrosion tested 316 L SS microloop sections were characterized with microscopy techniques to determine the microstructural and microchemical changes that occurred at the alloy/salt interface. It was found that 316 L SS showed heterogeneous dissolution at the hot-leg, whereas deposition of corrosion products occurred at the cold-leg. For the first time, experimentally obtained molten salt flow-induced corrosion of 316 L SS results were combined with computational thermodynamic-kinetic models to validate the dissolution and deposition in terms of elemental compositional changes at the alloy/salt interface. The thermodynamic-kinetic modeling predicted that the heterogeneous dissolution of Cr and Fe from the hot-leg section of 316 L SS persisted throughout the salt circulation. The model also estimated that, despite Cr deposition starting earlier than Fe, the total redeposition of Fe is expected to be significantly greater than that of Cr over the circulation of salt. Furthermore, the modeling accurately predicted the subsurface enrichment of Mo which is attributed to the reduced Cr activity and the relatively higher diffusion rate of Mo within the alloy matrix. Here, the agreement between modeling and experimental results confirms that Fe chlorides dissolve at the hot-leg and subsequently deposit at the cold-leg due to activity changes driven by the thermal gradient. By contrast, Cr was not detected in the cold-leg deposits, which is attributed to its weaker temperature dependence on activity, limiting its redeposition under these conditions.

36 - MATERIALS SCIENCE↗
NASA NTRS1987

Core plasma in the magnetosphere

The paper presents a review of new findings related to core plasma (which includes ions and electrons with energies less than 50 eV) identified in studies of magnetosphere during 1983-1986. These findings include the cleft ion fountain; the plasmaspheric heavy ion torus; core molecular ions; high core-plasma densities in the plasma-sheet boundary layer; intense transverse heating of core plasmas in the equatorial regions; the supersonic polar wind (H/+/ and He/+/); toroidal or ring heavy-ion distribution functions in the auroral region; conic-to-field-aligned ion signatures of inverted-V auroral events; and the large-scale electron density structure of the magnetosphere. The progress made on modeling core plasma includes simulation of thermal helium heating at the equator; kinetic models of the plasma transport from localized ionospheric regions; time-dependent hydrodynamic models of the polar wind and plasmasphere-ionosphere coupling; and kinetic models of plasmasphere refilling.

Horwitz, J. L.↗
NASA NTRS2008

Gaseous Surrogate Hydrocarbons for a Hifire Scramjet that Mimic Opposed Jet Extinction Limits for Cracked JP Fuels

This paper describes, first, the top-down methodology used to define simple gaseous surrogate hydrocarbon (HC) fuel mixtures for a hypersonic scramjet combustion subtask of the HiFIRE program. It then presents new and updated Opposed Jet Burner (OJB) extinction-limit Flame Strength (FS) data obtained from laminar non-premixed HC vs. air counterflow diffusion flames at 1-atm, which follow from earlier investigations. FS represents a strain-induced extinction limit based on cross-section-average air jet velocity, U(sub air), that sustains combustion of a counter jet of gaseous fuel just before extinction. FS uniquely characterizes a kinetically limited fuel combustion rate. More generally, Applied Stress Rates (ASRs) at extinction (U(sub air) normalized by nozzle or tube diameter, D(sub n or t) can directly be compared with extinction limits determined numerically using either a 1-D or (preferably) a 2-D Navier Stokes simulation with detailed transport and finite rate chemistry. The FS results help to characterize and define three candidate surrogate HC fuel mixtures that exhibit a common FS 70% greater than for vaporized JP-7 fuel. These include a binary fuel mixture of 64% ethylene + 36% methane, which is our primary recommendation. It is intended to mimic the critical flameholding limit of a thermally- or catalytically-cracked JP-7 like fuel in HiFIRE scramjet combustion tests. Our supporting experimental results include: (1) An idealized kinetically-limited ASR reactivity scale, which represents maximum strength non-premixed flames for several gaseous and vaporized liquid HCs; (2) FS characterizations of Colket and Spadaccini s suggested ternary surrogate, of 60% ethylene + 30% methane + 10% n-heptane, which matches the ignition delay of a typical cracked JP fuel; (3) Data showing how our recommended binary surrogate, of 64% ethylene + 36% methane, has an identical FS; (4) Data that characterize an alternate surrogate of 44% ethylene + 56% ethane with identical FS and nearly equal molecular weights; this could be useful when systematically varying the fuel composition. However, the mixture liquefies at much lower pressure, which limits on-board storage of gaseous fuel; (5) Dynamic Flame Weakening results that show how oscillations in OJB input flow (and composition) can weaken (extinguish) surrogate flames up to 200 Hz, but the weakening is 2.5x smaller compared to pure methane; and finally, (6) FS limits at 1-atm that compare with three published 1-D numerical OJB extinction results using four chemical kinetic models. The methane kinetics generally agree closely at 1-atm, whereas, the various ethylene models predict extinction limits that average ~ 45% high, which represents a significant problem for numerical simulation of surrogate-based flameholding in a scramjet cavity. Finally, we continue advocating the FS approach as more direct and fundamental for assessing idealized scramjet flameholding potentials than measurements of "unstrained" premixed laminar burning velocity or blowout in a Perfectly Stirred Reactor.

Pellett, Gerald L.↗
DOE OSTI2024

Insight into premixed diethoxymethane flames: Laminar burning velocities, temperatures, and emissions behaviour

Diethoxymethane ((CH 3 CH 2 O) 2 CH 2 , DEM) is a promising carbon-neutral fuel. DEM is a diether or acetal with a molecular structure similar to oxymethylene ethers (CH 3 O–(CH 2 O)n–CH 3 , OME n ). Thus, DEM can be expected to have a similar combustion behavior to OMEs, reducing harmful emissions such as NO x and particulate matter (PM) in internal combustion engines. From both experimental and kinetic modeling, fundamental studies on DEM are scarce in the literature. More studies are required to gain a detailed insight into the oxidation kinetics of DEM. Laminar burning velocity (LBV) is a critical property that allows a detailed assessment of the potential application of DEM in combustion devices. Unfortunately, the literature on the LBV of DEM is limited. Therefore, in this study we have investigated the LBV of DEM using two reactors for the first time, namely a heat flux burner and a combustion chamber. The experimental data is reported for equivalence ratio between 0.7 and 1.7, initial temperatures of 368–423 K, and initial pressure of 1–5 bar. In addition, we developed a detailed kinetic model extending our recent work of Shrestha et al. (Combust. Flame. 246 (2022) 112,426) to characterize the combustion behavior of DEM utilizing the new experimental data from this work and the literature data. Our model performs remarkably well in capturing the newly measured LBV experimental data over various experimental conditions. We found that DEM and dimethoxy methane (DMM) have similar values of LBVs (within ±1.5 cm/s) for a given condition, which indicates that intermediate chemistry governs the flame chemistry. Despite DEM being a larger molecule that is expected to have slightly lower LBVs than DMM, its effect on the measured values of LBVs is negligible. Finally, we experimentally measured NO x formation in DEM flame for the first time. The stochiometric flame has the highest NO x formation. The proposed model predicted the equivalence ratio dependence of NO x nicely. However, it overestimates the NO x formation for stoichiometric DEM/air mixtures by ~30 %. The model suggests that the thermal NO formation route is favored at lean and stochiometric conditions. In contrast, the prompt NO formation route is enhanced for rich mixtures.

10 SYNTHETIC FUELS↗
DOE OSTI2023

Fuel Flow Modeling under Thermal Fluidics Considerations within Operational Domains of Molten Salt Reactors

Molten Salt Reactors (MSRs) offer significant versatility supported by their design flexibility and safety advantages due to high temperatures, design customization capabilities, strong negative reactivity feedback effects, and salt chemistry characteristics. Different control schedules and different mass flow rates in the MSR primary systems impact neutronics and thermal hydraulics coupling. The result is two-fold. First, it leads to adaptability in various applications when the customization of a MSR unit is desired. Second, it calls for the quantification of dynamics and safety characteristics to assess the operational domain and make sure reactor stability and inherent safety characteristics are maintained. This paper provides an in-depth transient analysis of the operational implications of MSRs having heat-generating fuel and heat-transporting salt mixed forming liquid fuel salts contained within primary systems. A coupling model was developed to simulate fuel flow in molten salt reactors. A zero-dimensional reactor kinetics model was used with a one-dimensional heat transfer model in the reactor core. The temperature reactivity feedbacks of the fuel salt resulting from the thermal-hydraulics model were used in the reactor kinetics model to complete the coupled code. The model was applied to simulate the Molten Salt Breeder Reactor (MSBR) response at steady state and transients. The developed model was shown to be a suitable tool for the dynamic analysis of molten salt reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1995

Convective and radiative heating for vehicle return from the Moon and Mars

The aerothermal environment is examined for two vehicle forebodies near the peak heating points of lunar and martian return-to-earth trajectories at several nominal entry velocities. The first vehicle forebody is that of a 70 deg aerobrake for entry into earth orbit; the second, a capsule of Apollo configuration for direct entry into the earth's atmosphere. The configurations and trajectories are considered likely candidates for such missions. Two-temperature, thermochemical nonequilibrium models are used in the flow field analyses. In addition to Park's empirical model for dissociation under conditions of thermal nonequilibrium, the Gordiets kinetic model for the homonuclear dissociation of N2 and O2 is also considered. Temperature and emission profiles indicate nonequilibrium effects in a 2 to 5 cm post shock region. Substantial portions of the shock layer flow appear to be in equilibrium. The shock layer over an aerobrake for return from the moon exhibits the largest extent of nonequilibrium effects of all considered missions. Differences between the Gordiets and Parks kinetic model were generally very small for the lunar return aerobrake case, the greatest difference of 6.1 percent occurring in the radiative heating levels.

Greendyke, Robert B.↗