Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixture 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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Scramjet Development Tests Supporting the Mach 10 Flight of the X-43

The Hyper-X Project s successful third flight of the X-43 at near Mach 10 in 2004 proved the potential for airbreathing propulsion at hypersonic speeds. The engine flowpath used in the X-43 research vehicle was developed and evaluated in a systematic series of ground tests in the NASA HyPulse Shock Tunnel at conditions duplicating Mach 10 flight using a full scale height, partial width engine model of the flight engine. Tests were conducted over a range of equivalence ratios from 0.8 to 1.6 using hydrogen and a mixture of two-percent silane in hydrogen fuels. Silane gas was used as an ignition aid during the short duration of the pulse facility tests. Variation of the engine inflow conditions, pressure, temperature, and Mach number, were parametrically varied during the test entries to broaden the database over the expected uncertainty in the flight conditions. A review of the ground test technique and comparisons of the ground test pressures along with selected flight data are presented.

Rogers, R. C.↗
NASA NTRS2001

Iron Catalyst Chemistry in High Pressure Carbon Monoxide Nanotube Reactor

The high-pressure carbon monoxide (HiPco) technique for producing single wall carbon nanotubes (SWNT) is analyzed using a chemical reaction model coupled with properties calculated along streamlines. Streamline properties for mixing jets are calculated by the FLUENT code using the k-e turbulent model for pure carbon monixide. The HiPco process introduces cold iron pentacarbonyl diluted in CO, or alternatively nitrogen, at high pressure, ca. 30 atmospheres into a conical mixing zone. Hot CO is also introduced via three jets at angles with respect to the axis of the reactor. Hot CO decomposes the Fe(CO)5 to release atomic Fe. Cluster reaction rates are from Krestinin, et aI., based on shock tube measurements. Another model is from classical cluster theory given by Girshick's team. The calculations are performed on streamlines that assume that a cold mixture of Fe(CO)5 in CO is introduced along the reactor axis. Then iron forms clusters that catalyze the formation of SWNTs from the Boudouard reaction on Fe-containing clusters by reaction with CO. To simulate the chemical process along streamlines that were calculated by the fluid dynamics code FLUENT, a time history of temperature and dilution are determined along streamlines. Alternative catalyst injection schemes are also evaluated.

Scott, Carl D.↗
NASA NTRSDate not supplied

Open-source Numerical Modeling of Solidification Cracking Susceptibility: Application to Refractory Alloy Systems

Introduction. Alloys such as aluminum, nickel-base, and austenitic stainless steels are susceptible to solidification cracking during welding and 3D printing. Compositional optimization is one method used to effectively mitigate solidification cracking of those alloy systems. With the surge in hypersonic and in-space propulsion activities, refractory metals (Nb, Mo, Ta, W, and Re) and their alloy derivatives are increasing in importance due to their extreme high melting point and retention of high-temperature strength; however, their chemistry was most typically optimized to promote ductility during mechanical operations such as drawing and forming. Welding of such alloys has been a challenge due to a number of issues including solidification cracking, atmospheric contamination (O, C, and N), as well as a shift in ductile-to-brittle transition to higher temperature following grain growth induced by welding. Compositional optimization of refractory alloys for solidification cracking resistance in particular is desirable as their usage increases with the advent of advanced manufacturing methods such as 3D printing. This work evaluates the effect of compositional variation in refractory metal systems on the solidification cracking susceptibility with the goals of optimizing existing alloys and joining process techniques, and formulating new alloys with increased solidification cracking resistance. Experimental Procedures. A python code was developed in a Jupyter notebook environment (Michael and Sowards, 2023) to facilitate the calculation of crack susceptibility index proposed by Kou (2015). Composition is entered as a single point, or as a 1-D or 2-D array. The notebook calls pycalphad (Otis and Liu, 2017 and Bocklund et al, 2020) to calculate the evolution of fraction solid as a function of temperature (under either Scheil or equilibrium assumptions) and then evaluates steepness of the fraction solid curve near the terminal stage of solidification to predict solidification cracking resistance. Open source thermodynamic databases available at online repositories are used (van de Walle). The process is setup in an automated fashion to generate plots that show variation in solidification cracking susceptibility according to composition on 1-D line plots or 2-D contour plots. The Jupyter notebook and crack susceptibility algorithm was also integrated with a widely used commercial CALPHAD code for validation and alloy exploration. Results and Discussion. The crack susceptibility model was first validated against a series of refractory alloy compositions evaluated in past work which utilized a specialized Varestraint test built inside a vacuum chamber environment (Lessman and Gold, 1971). The alloys tested in the Varestraint apparatus included T-111 (Ta-8W-2Hf), ASTAR-811C (Ta-8W-1Re-0.7Hf-0.025C), FS-85 (Nb-27Ta-10W-1Zr), T-222 (Ta-9.6W-2.4Hf-0.01C), Ta-10W, B-66 (Nb-5Mo-5V-1Zr), and SCb-291 (Nb-10W-10Ta). The initial test of the model showed a strong correlation with empirical Varestraint data, i.e., a Spearman rank correlation between model predictions and hot cracking measurements was observed to be greater than 0.8. Following the validation, a set of refractory metal binary mixtures was investigated to evaluate sensitivity of Nb, Mo, W, and Ta to C, N, and O content. A series of plots were produced that suggest ppmw ranges of C, N, and O where solidification cracking increases significantly and reaches a maximum. Also comparative ranking of each primary refractory metal to each interstitial was produced. For example C produces greater cracking response in Mo whereas O produces greater cracking response in Ta and Nb. Such compositional values have utility in setting limits on pickup of these interstitial elements during welding and printing rather than using a one-size-fits-all approach. Furthermore, the results have use in determining additive powder recycling requirements, which is especially pertinent for refractory metal powders due to their high cost compared to conventional alloys. Another application created thousands of hypothetical alloys within the nominal specified composition range of two widely used refractory alloys C103 (Nb-10Hf-1Ti) and TZM (Mo-0.5Ti-0.1Zr). The cracking index was calculated for the alloys and results were fed into machine learning regression techniques including Multiple Linear Regression, Ridge Regression, and Lasso Regression to determine relative potency each alloying element had on computed solidification cracking index. A series of linear equations were produced that relate composition of C103 and TZM to solidification cracking index. The crack susceptibility of C103 for example is described by an equation of the form: cracking index ~ O + 0.667*C + 0.635*N + 0.00037*Ta – 0.0008*Hf (in wt.%) From that equation, it is clear that O has strong propensity to induce solidification cracking. Interestingly, Hf is shown to reduce calculated cracking response. Finally, realizing the potential of this method to discover new refractory alloy formulations across the period table that have low solidification cracking sensitivity, the code was applied to new untested alloy systems including W-Zr-C, W-Ta-C, and others. Conclusions. In summary, an open source numerical method has been developed using Python code to calculate Kou’s crack susceptibility index. The method was applied to refractory metals which are inherently difficult to study from a weldability testing standpoint since inert shielding gas is not sufficient and welding is typically done in vacuum, especially in light of findings presented here where oxygen has profound influence on solidification cracking. This work revealed the effect of compositional variations on a series of refractory metals and showed the framework defined here will be useful in 1) the development of new alloys that have improved weldability and 3D printability, 2) placing compositional limits on existing alloys, and 3) ensuring adequate controls of manufacturing processes such as 3D printing where powder reuse is critical. Keywords. pycalphad; Python; refractory metals; solidification cracking. References. B. Bocklund et. al. (2020) http://doi.org/10.5281/zenodo.3630657. S. Kou. (2015) https://doi.org/10.1016/j.actamat.2015.01.034. G.G. Lessmann and R.E. Gold. Welding Journal, issue 1, pp. 1-s – 8-s (1971). F.N. Michael and J.W. Sowards. NASA/TM-20230002218 (2023). R. Otis and Z.-K. Liu. (2017) http://doi.org/10.5334/jors.140. A. Van de Wallle et. al. (2018) https://doi.org/10.1016/j.calphad.2018.04.003.

pycalphad↗
NASA NTRSDate not supplied

Effects of Mild Hypobaric Hypoxia on Response From Mild Exercise in NASA Exploration Atmosphere Tests

INTRODUCTION: The NASA Exploration Atmosphere study aims to validate a new prebreathe protocol incorporating an alternate habitat atmosphere of 56.5 kPa (8.2psia), 34% O 2 , and 66% N 2 to control and mitigate decompression sickness risk associated with spaceflight extravehicular activities. This alternate atmosphere results in a mild hypoxic environment (PIO 2 of 128 mmHg) that may influence inhabitant physiological responses to exercise. Therefore, the responses to light exercise in relation to the mild hypobaric hypoxic environment were investigated. METHODS: Eight participants (4M/4F; age=38.3±9.0 yr; weight=76.0±13.1 kg; peak aerobic capacity [VO 2 pk]= 3.1±0.7 L/min) were exposed to a mild hypobaric hypoxic environment for 11 days in NASA’s 20-foot hypobaric chamber at Johnson Space Center. Participants alternated test days between a habitat atmosphere of 56.5 kPA/34% O 2 and simulated 6-hour EVA environment, 29.6 kPa/85% O 2 . Participants completed pre-mission graded VO 2 pk tests on a LODE cycle ergometer with ParvoMedics metabolic analyzer. Submaximal aerobic tests were performed pre and during the mission (10-minutes of exercise at a workload of 40% VO 2 pk).Pre-mission, performed submaximal exercise exposed to a breathing air mixture of, 18% O 2 and balance N 2 , to simulate the mildly hypoxic environment within the 20-foot chamber. Linear mixed models (fixed effect: test day; random effects: subject, age) were performed to determine whether physiological responses (oxygen uptake [VO 2 ], carbon dioxide production [VCO 2 ], ventilation [VE], oxygen saturation [SPO 2 ], heart rate [HR], respiratory exchange ratio [RER]) to submaximal exercise performed within the 20-foot hypobaric chamber every 2 days differed from the pre-mission (18% O 2 ) testing. Results are presented as estimated marginal means with lower and upper confidence limits, with significance set to 0.05. RESULTS AND DISCUSSION: Seven of eight participants completed the chamber study. Mixed models for VO 2 , VCO 2 , SpO 2 , HR, and RER indicated negligible impact of the mild hypoxic chamber environment compared to pre-mission testing (all p>0.05).Results for VE indicated a minor impact of chamber environment compared to pre-mission testing (37.0 L/min [30.8, 43.2]), with increased VE at test day 10 (40.4 L/min [34.2, 46.6]; p= 0.029).The aerobic performance data collected suggests limited physiologic responses to mild exercise when performed at normobaric hypoxia and hypobaric hypoxia.

B Siders↗
NASA NTRS2024

Enstatite Chondrite Outgassing and Condensate Formation: Implications for Early Atmosphere Development

Early atmospheres on rocky planets where life may develop form through outgassing of their original starting blocks, likely a mixture of chondritic material (carbonaceous chondrites (CC), ordinary chondrites (OC), and enstatite chondrites (EC)). However, there is limited experimental data to inform models connecting a planet’s bulk composition to its early atmospheric properties and thus its possibility for life. Thompson et al. (2021) took a major step for-ward in exploring this knowledge gap by measuring outgassing of 3 volatile-rich CCs, providing important experimental constraints on the initial chemical com-position of early rocky planet atmospheres. These data provided novel insights into the gas chemistry released into evolving atmospheres early in a rocky planet’s history and differed from those currently assumed by many theoretical models of rocky planet atmosphere formation. In this study, we focused on the outgassing and condensation reactions of a primitive EC3, which has a lower intrinsic oxygen fugacity ( ƒ O 2 ) and lower volatile content than the CCC meteorites investigated by. We selected the EC to explore differences between EC and CC in low-pressure, high-temperature outgassing and condensation including S, Cl, Na, and C species.

B A Anzures↗
DOE OSTI2021

Potential Economic Values of Low-Vapor-Pressure Gasoline-Range Bio-Blendstocks: Property Estimation, and Blending Optimization

Adding non-traditional blendstocks derived from biomass into gasoline has gained renewed interest since it will reduce greenhouse gas emissions and potentially enhance fuel properties and boost engine efficiencies. A competitive bio-blendstock candidate should have either lower prices or superior properties or both relative to conventional gasoline blending components. Gasoline is a mixture of hundreds compounds specified by properties rather than compositions. The addition of bio-blendstocks, most likely oxygenates, will modify the property mixing rules due to the non-ideal interactions between polar and nonpolar components. This paper presents an equation-of-state-based model for predicting Reid vapor pressure, reduced order non-linear property blending models for computing key properties of final gasoline products, and a blending optimization approach to identify key economic drivers. This paper also evaluates the economic values of bio-blendstocks to petroleum refiners using those models. In additional to ethanol, six low-vapor-pressure bio-blendstock candidates were evaluated: i-propanol, n-propanol, i-butanol, diisobutylene, cyclopentanone, and a mixture of furans. Reid vapor pressure, distillation temperatures, and octane numbers were identified as the key economic drivers of adding bio-blendstock. The economic values of those bio-blendstocks range from 2.22 to 4.00 US dollar per gasoline gallon equivalent in 2013 to 2017 5-year averaged price basis. The bio-blendstocks with higher octane numbers may have higher economic values. The uncertainties in property predictions may lead to roughly 15% deviation in the economic value evaluation.

Jiang, Yuan↗
DOE OSTI2024

Hydrogen Leak Modeling for Development of Smart Distributed Monitoring Under Unintended Releases

Hydrogen is a versatile and clean energy carrier that can be produced from various renewable sources such as wind, solar, and hydropower. Hydrogen has the potential to play a crucial role in decarbonizing industrial processes that are currently reliant on fossil fuels and provide long-duration and/or seasonal energy storage to enable electricity decarbonization. Hydrogen can also be used as a fuel for fuel cell vehicles, providing a zero-emission alternative to traditional internal combustion engines. DOE launched the Hydrogen Energy Earthshot (Hydrogen Shot) in June 2021 to reduce the cost of clean hydrogen by 80% to $1 per 1 kilogram in 1 decade ("1 1 1"). While promising, Hydrogen is highly-flammable, and in the presence of oxygen, it can form explosive mixtures. . Therefore, understanding leak scenarios is essential to evaluate and mitigate the safety risks associated with potential hydrogen leaks. An increased understanding of leak behavior, and having tools to model leaks, can help assess how hydrogen would disperse in different environments, influencing emergency response plans and safety measures, and identify potential issues with materials and design systems that can withstand the challenges posed by hydrogen. Recently, researchers have attempted to study hydrogen leaks for development of risk management strategies. However, the focus has been on closed or semi-closed spaces like storage rooms, vehicles, garages, and fueling stations - all promising locations for future hydrogen infrastructure. In this presentation, the modeling environment extends the span of research further by modeling hydrogen leak in an outdoor, open space. We will present the key challenges with modeling hydrogen leaks in an uncontrollable environment, how they were handled, and how modeling results informed sensor selection and placement. A Hydrogen research facility at the National Renewable Energy Laboratory (NREL) was used as a case study to model hydrogen leaks. In the future, Hydrogen wide area detection methodologies will be developed and tested at this site to monitor for unintended and operational hydrogen releases. The data generated from modeling will be used to develop a predictive model to detect hydrogen leak location based on concentration measured by sensors in this open space. Furthermore, the facility was also chosen because controlled hydrogen releases can be performed. A computational fluid dynamics (CFD) based modeling approach was taken to model hydrogen leak. The full-scale hydrogen facility was modeled with a large ambient domain. The electrolyzer at the facility can produce a controlled release rate of 27 kg-H2/hr. Site-specific atmospheric and weather condition data such as wind direction, wind speed at various altitudes, and temperature were used as inputs to the model. To capture the variability of weather conditions, a subset of the weather conditions experienced during daytime hours without precipitation over the course of three months was generated; using established data clustering techniques, a total of 100 condition sets were chosen. The results show statistical distributions and ranges of hydrogen concentrations at locations throughout the domain. These distributions are compared to experimental data from a constant mass flow, controlled hydrogen release at the facility. The stochastic wind conditions of the release make direct validation difficult, therefore, statistical comparison approaches were used. Wind conditions are found to significantly impact the release behavior, including direction and concentration. Sensor selection and placement is proposed for the facility and is now based on release behavior predicted for the facility given its weather patterns; this is much more informed than without the modeling results. The methodology and analysis procedure can be translated to other facilities using modified geometries and site-specific weather conditions. Hydrogen holds great promise as a renewable energy fuel, but ensuring safety in its production, storage, and use is paramount. Studying potential leak scenarios in an open space will help develop sensors to detect hydrogen on a large spectrum of concentration and eventually build a smart distributed monitoring system.

CFD↗
NASA NTRSDate not supplied

A Warm Garage for a Lunar Rover

Approach: One approach to heating a rover during the lunar night is the so-called thermal wadis concept [1]. This involves heating the regolith with solar concentrators and placing the rover on the heated surface for the night. Since the regolith is heated by a relatively weak heat flux, a high thermal conductivity is required for heating a sufficiently large mass of regolith. However, lunar regolith has a low thermal conductivity. Therefore, the concept involves increasing the conductivity by sintering the regolith, which requires a significant energy input and complex procedures. Here we propose an alternative approach where the low thermal conductivity of regolith is an advantage. Specifically, we propose to use a highly exothermic combustible mixture for heat generation. The mixture pellets are placed in the surface layer of regolith and ignited. The combustion forms condensed products and releases heat, which then slowly spreads to the surrounding regolith. Heat can also be transferred, for example, by heat pipes, into radiant heating surfaces installed on the ground. A greenhouse that transmits sunlight during the day and decreases the radiative heat losses during the night can also be installed. Selection of the Heat-generating Mixture: The reactive mixture should have a high specific energy and generate only condensed products since gases could disturb the regolith layer, carry enthalpy out of the system, and lead to an explosion. There are mixtures, (sometimes called pyrolants) that possess very high specific energies. One example is magnesium-Teflon-Viton mixtures used in flares. However, they produce gases and may cause explosions. Other mixtures that include magnesium cannot be used either because of the high vapor pressure of Mg at temperatures well below the combustion temperature. Recently, mixtures that involve lithium peroxide (Li2O2) have been proposedfor using in space power systems [2]. However, they produce lithium oxide (Li2O), which boils at 2800 K at 1 atm and hence at a lower temperature in vacuum. Fortunately, there exist many mixtures that release a lot of heat and form only condensed products during the combustion. Many such mixtures have been used for self-propagating high-temperature synthesis (SHS) of various materials [3, 4]. For the application discussed here, t itanium/boron (1:2 mole ratio) mixture appears to be particularly promising. The specific energy is 4.0 MJ/kg (1.1 kWh/kg), the adiabatic flame temperature is about 3200 K, and the reaction forms solid titanium diboride (TiB2, melting point: 3500 K). The mixture can be ignited easily with a heated tungsten wire, and it has been used widely as a booster to ignite the main mixture in the SHS process.Estimates: Assuming that specific heat of regolith is 500 J/(kg∙K) [5] and all generated heat is transferred to regolith, 12.5 kg of the Ti/B mixture would increase the temperature of 1000 kg of regolith by 100 K. To evaluate the rate of heat transfer in the regolith, a spherical model was analyzed where the heat released by a 12.5 kg Ti/B core propagates by thermal conduction through a 1000 kg regolith shell with no heat loss from its outer surface. At a bulk density of 1500 kg/m3 [5], the radius of the shell was 54 cm, while the radius of the core was about 11 cm. The calculations were conducted using Thermal Desktop SINDA/FLUINT (Cullimore and Ring Technologies) software at two constant values of bulk thermal conductivityk of the regolith: 0.001 and 0.01 W/(m∙K). The results show that after 14.5 days the core lost 31% of the released heat at the lower k and 77% at the higher k. At a distance of 20 cm from the core surface, the temperature of the regolith increased by only 1 K at the lower k and by 132 K at the higher k. In reality, the regolith near the heat source will be melted, so its thermal conductivity will increase significantly. Nevertheless, the conducted estimates indicate that combustion-based heat generators, placed directly in the regolith, could provide heat during a rather long period such as the lunar night.Conclusion: Heat generators based on gasless combustion of highly energetic reactive mixtures could be installed directly in the surface layer of lunar regolith. Because of the low thermal conductivity of the regolith, such generators would keep thermal energy for days and gradually supply heat to a rover/lander.Acknowledgment: The material presented in this work is based upon the work supported by National Aeronautics and Space Administration (NASA) under Grant #80NSSC20K0293.References: [1] Balasubramaniam R. et al. (2011) J. Thermophys. Heat Trans., 25,130−139. [2] Blair R.G. and Vasu S.S. (2022) Conf. Advanced Power Systems for Deep Space Exploration. [3] Varma A. et al. (1998) Adv. Chem. Eng., 24,79−226. [4] Levashov E.A. et al. (2017) Int. Mater. Rev., 62,203−239. [5] Wood-Robinson R. et al. (2019) J. Geophys. Res. Planets, 124, 1989−2011.

lunar↗
DOE OSTI2021

Probing the Chemical Kinetics of Minimalist Functional Group Gasoline Surrogates

We report surrogate mixtures are routinely used for understanding gasoline fuel combustion in engine simulations. The general trend in surrogate formulation has been to increase the number of fuel components in a mixture to better emulate real fuel properties. Recently, a new surrogate design strategy based on functional group analysis of real gasolines was proposed using a minimal number of species [minimalist functional group (MFG)-approach]. MFG surrogates (having just one or two components) could experimentally capture the ignition delay time (IDT), threshold sooting index, and smoke point of different gasoline fuels with hundreds of components. However, other combustion characteristics were not explored, and kinetic modeling of MFG surrogates was not reported. These aspects are addressed in this paper, where the combustion behavior of MFG surrogates for various gasolines was assessed by simulating IDT, jet-stirred reactor oxidation, and premixed laminar flame speeds using chemical kinetic modeling. MFG simulations were compared with experimental data of the real gasolines as well as with the more complex multicomponent (five to nine species) surrogates. This study reveals that binary MFG surrogate mixtures are capable of accurately simulating the combustion behavior of more complex gasoline fuels with hundreds of components.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

LES and RANS Spray Combustion Analysis of OME 3-5 and n-Dodecane

Clean-burning oxygenated and synthetic fuels derived from renewable power, so-called e-fuels, are a promising pathway to decarbonize compression–ignition engines. Polyoxymethylene dimethyl ethers (PODEs or OMEs) are one candidate of such fuels with good prospects. Their lack of carbon-to-carbon bonds and high concentration of chemically bound oxygen effectively negate the emergence of polycyclic aromatic hydrocarbons (PAHs) and even their precursors like acetylene (C 2 H 2 ), enabling soot-free combustion without the soot-NO x trade-off common for diesel engines. The differences in the spray combustion process for OMEs and diesel-like reference fuels like n-dodecane and their potential implications on engine applications include discrepancies in the observed ignition delay, the stabilized flame lift-off location, and significant deviations in high-temperature flame morphology. For CFD simulations, the accurate modeling and prediction of these differences between OMEs and n-dodecane proved challenging. This study investigates the spray combustion process of an OME 3 − 5 mixture and n-dodecane with advanced optical diagnostics, Reynolds-Averaged Navier–Stokes (RANS), and Large-Eddy Simulations (LESs) within a constant-volume vessel. Cool-flame and high-temperature combustion were measured simultaneously via high-speed (50 kHz) imaging with formaldehyde (CH 2 O) planar laser-induced fluorescence (PLIF) representing the former and line-of-sight OH* chemiluminescence the latter. Both RANS and LES simulations accurately describe the cool-flame development process with the formation of CH 2 O. However, CH 2 O consumption and the onset of high-temperature reactions, signaled by the rise of OH* levels, show significant deviations between RANS, LES, and experiments as well as between n-dodecane and OME. A focus is set on the quality of the simulated results compared to the experimentally observed spatial distribution of OH*, especially in OME fuel-rich regions. The influence of the turbulence modeling is investigated for the two distinct ambient temperatures of 900 K and 1200 K within the Engine Combustion Network Spray A setup. The capabilities and limitations of the RANS simulations are demonstrated with the initial cool-flame propagation and periodic oscillations of CH 2 O formation/consumption during the quasi-steady combustion period captured by the LES.

CFD↗
DOE OSTI2025

Thermodynamic modeling of aqueous acetic acid, butyric acid, and lactic acid solutions

Based on the activity coefficient – fugacity coefficient approach, a rigorous thermodynamic modeling study is presented for accurate correlation of vapor-liquid equilibrium data of aqueous solutions of acetic acid (293 to 391 K), butyric acid (325 to 436 K), lactic acid (378 to 409 K), and acetic acid + butyric acid binary mixture (358 to 421 K). In addition, the pH data of the three aqueous, single carboxylic acid solutions were measured at 298 to 328 K and successfully correlated. Given that these aqueous carboxylic acid solutions exhibit various degrees of association behavior in both vapor and liquid phases, the thermodynamic models considered for this study include the Redlich-Kwong equation of state (RK-EoS) and the Hayden-O’Connell equation of state (HOC-EoS) for the vapor phase fugacity coefficients and the electrolyte non-random two-liquid model (eNRTL) and the association electrolyte non-random two-liquid model (AeNRTL) for the liquid phase activity coefficients. The combination of the HOC-EoS for the vapor phase and the AeNRTL model for the liquid phase is found to provide the best correlation results, consistent with the fact that the HOC-EoS and the AeNRTL model explicitly account for association behaviors in the vapor phase and the liquid phase, respectively.

09 BIOMASS FUELS↗
NASA NTRSDate not supplied

Integrated Bosch Process System Models for In-Situ Oxygen and Carbon Production

In-Situ Resource Utilization (ISRU) technology is a vital component to NASA’s mission of a sustainable presence on the Moon and Mars. Local resources can be leveraged to reduce resupply frequency and mass. Elements of the Bosch process, combined with the carbothermal reduction process, can produce oxygen on the lunar surface with minimal consumables. The Bosch process can also produce oxygen on the Martian surface by using the CO 2 -rich environment. Between both systems, adsorption pump, solar thermal energy, carbon formation reactor, and water recovery subsystems are modeled and integrated to create a functional model in MATLAB software. The model is used to simulate performance of the system and reduce mass, power, and volume requirements. This integrated system model provides a tool to scale ISRU technologies for oxygen and carbon production. The MATLAB model is created by developing a system of independent subsystem models that are solved for their quasi-steady state values which can be integrated with respect to time to determine the change in current states. A flexible time stepping method is used to ensure a high level of accuracy during periods of rapid change while still making use of a simple explicit integration method. The flexible time step is calculated for each independent subsystem and the minimum value from those is used as the overall time step. A flexible time step is calculated by dividing a resolution value, or the maximum change per time step, by the variables current rate of change. The maximum value from all points in space is used for subsystem models that contain multiple values. The process is done for every variable that is being monitored in each subsystem and the global minimum is used as that iteration’s timestep. Several assumptions used in the MATLAB model for fluid flow dynamics, such as 1-D gas flow through the sorption pump, are supported by modeling in Ansys Fluent software. The Lunar oxygen production system is outlined in Fig. 1. The carbothermal reduction subsystem uses solar energy to heat a mixture of lunar regolith and carbon powder to produce carbon monoxide. To begin, the carbon monoxide feeds to the modified Bosch subsystem along with hydrogen gas. The reactants then enter the carbon formation reactor where water and carbon powder are produced. Solar thermal energy is used to add energy to the reactor, but waste heat from the carbothermal process is another potential heat source. The water is collected and electrolyzed to produce hydrogen which reenters the Bosch subsystem, and the oxygen is stored for downstream use. The carbon powder is collected and feeds back into the carbothermal subsystem. The Martian oxygen production system uses the full Bosch process and is outlined in Fig 2. A CO 2 adsorption pump thermally cycles to scrub and pressurize CO 2 from the environment. Along with an initial supply of hydrogen, the reactants enter the Reverse Water Gas Shift Reactor (RWGSR) which produces carbon monoxide and water. Carbon monoxide and unreacted hydrogen enter the carbon formation reactor to produce water and carbon powder. The water is collected from both reactors and electrolyzed to reintroduce hydrogen and store oxygen for propellant production or life support. Carbon is removed from the carbon formation reactor and stored. The adsorption pump utilizes rapid cycle temperature swings within a stack of zeolite coated surfaces. The subsystem model solves 1-D quasi-steady conservation laws of the quasi-steady form, shown in Eq. 1, for the gas stream and heat exchange liquid to predict performance parameters such as breakthrough capacity and optimum cycle time. The source term S is used to capture interactions between the fluid flows and the sorbent. A quasi-steady-state scheme is used where no time derivatives appear in the governing equations, except for those in the source terms. This results in an autonomous system, where ∂F/∂x = ƒ(F). The fluxes F are provided at the inlet, and an explicit method is used to solve for the spatial distribution of F. The heat and mass flows to the sorbent are then extracted from the source terms. These flows are numerically integrated to produce a 1-D solution for the system’s state as a function of both time and space. The body of the adsorption pump is separated into two semi-independent models: the heat exchanger fluid flow and gas flow through the zeolite coated surfaces. Both models are solved using the above-described method to find a 1-D solution as a function of space and interact only once a timestep is taken. The interaction point is the sorbent through which all heat transfer between the two models must occur. Sorbent mass adsorption is calculated using the Lagergren model, shown in Eq. 2, where the transfer coefficient, λ D , is found by solving a system of nondimensionalized equations derived by using the heat and mass transfer analogy for transport phenomena. Using Grade 544 Type 13X zeolite as the sorbent material, the equilibrium concentration, θ eq , is calculated using the k-site Langmuir isotherm and fit parameters. Additionally, the enthalpy of adsorption used in the model is computed by interpolation of available data [1]. The subsystem model was validated using the Rapid Cycle Temperature Swing Adsorption (RC-TSA) pump. The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC). The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC).

In situ Resource Utilization↗
DOE OSTI2021

On Two-Phase Pressure and Temperature Equilibration with Mie-Grüneisen Equations of State

In the course of reactive flow modeling, the phases of high explosive (HE) reactants and products introduce additional degrees of freedom to the underlying hydrodynamic equations of compressible reactive fluid mechanics. Hydrodynamic simulations track the bulk (mixture) conservation properties of mass, momentum and energy. When multiple phases are introduced, one must decide on appropriate closure rules to determine the relative phase volumes and energy distribution between the phases. One typical closure is to assume that the phases are in pressure and temperature equilibrium. Here, the requisite equations of p-T equilibration are explored for Mie-Grüneisen Equations of State (EOS). Examples with Davis reactants and Davis products EOSs will be given. The Davis forms are used in both the WSD and AWSD reactive flow models.

42 ENGINEERING↗
DOE OSTI2022

Combustion Characteristics of a Tier II Gasoline Certification Fuel and its Surrogate with Iso-butanol: experiments and detailed numerical modeling

This presentation will describe results of an experimental and numerical study of combustion of a Tier II gasoline certification fuel and its surrogate blended with iso-butanol as a potential biofuel additive. The experiments were carried out under conditions where gas transport originated entirely by fuel evaporation with external convection (both forced and natural) were minimized. The resulting near spherically symmetric droplet flames were amenable to detailed numerical modeling. The droplet initial diameters were in the narrow range of 0.55 mm to 0.62 mm while and burning is performed under standard atmosphere. The combustion characteristics such as evolution of droplet, flame and soot shell diameters for indolene/iso-butanol and surrogate/iso-butanol mixtures are compared. The presence of iso-butanol is experimentally found to reduce the formation of soot as the iso-butanol concentration was increased. Simulations are also reported using a detailed kinetic chemistry comprised of 324 species and 17608 reactions using the CRECK POLIMI mechanism. The detailed numerical model included unsteady gas and liquid transport, soot chemistry and radiative transport. Droplet burn rates were well predicted while flame diameters showed more variability depending on how the flame was defined in the numerical simulation (peak gas temperature and OH concentration). Simulations also show the influence of iso-butanol concentration on flame temperature, greenhouse gas emissions and particulate emissions.

Guo, Songtao↗
NASA NTRS1994

Pdf modeling for premixed turbulent combustion based on the properties of iso-concentration surfaces

In premixed turbulent flames the presence of intense mixing zones located in front of and behind the flame surface leads to a requirement to study the behavior of iso-concentration surfaces defined for all values of the progress variable (equal to unity in burnt gases and to zero in fresh mixtures). To support this study, some theoretical and mathematical tools devoted to level surfaces are first developed. Then a database of direct numerical simulations of turbulent premixed flames is generated and used to investigate the internal structure of the flame brush, and a new pdf model based on the properties of iso-surfaces is proposed.

Vervisch, L.↗
DOE OSTI2025

Secondhand Exposure to Simulated Cannabis Vaping Aerosols

Emissions from cannabis vaping degrade indoor air quality and expose non-users to secondhand pollutants. We investigated how the vaping mixture composition affects indoor aerosol characteristics and exposures. Simulated cannabis vaping aerosol was produced by flash evaporation in a 20 m3 chamber of mixtures containing terpenoids, cannabinoids, cannabis extract constituents, and the adulterant vitamin E acetate (VEA). Aerosol time- and size-resolved concentrations (8 nm-2.5 μm at 1 Hz) were measured, and a dosimetry model was used to evaluate the intake of secondhand aerosols. The results showed peak particle number (PN) concentrations between 0.7 × 106 and 13 × 106 cm-3 and peak mass concentration (PM1.0) between 65 and 1191 μg m-3 at t = 5 min after emission. Concentrations decreased to 21-57% of peak PN and 33-69% of peak PM1.0 at t = 60 min. The PM1.0 yield was 0.06 for a terpenoid-only mixture, 0.22-0.36 for tetrahydrocannabinol (THC)-terpenoid mixtures, and >1 for mixtures containing high concentrations of cannabidiol (CBD) or VEA. For intake deposition, the highest aerosol fraction was deposited in the pulmonary region, followed by the tracheobronchial and head regions. Deposition increased in the presence of THC, CBD, or VEA, with aerosols <100 nm contributing the majority of particles deposited in all regions.

Tang, Xiaochen↗
NASA NTRS1962

Jet Interference Effects on a Model of a Single-Engine Four Jet V/STOL Airplane at Mach Numbers from 0.60 to 1.00

An investigation was conducted in the Langley 16-foot transonic tunnel to determine the interference from four exhaust jets on the aerodynamic characteristics of a model of a V/STOL airplane. The single- engine four-jet turbofan power plant of the airplane was simulated by inducing tunnel airflow through two large side inlets and injecting the decomposition products of hydrogen peroxide into the internal flow. The heated gas mixture was exhausted through four nozzles located on the sides of the fuselage under the wing, two near the wing leading edge and two forward of the trailing edge; the nozzles were deflected downward 1.5 deg and outward 5.0 deg to simulate cruise conditions. The wing of the model was a clipped delta with leading-edge sweep of 40 deg, aspect ratio of 3.06, taper ratio of 0.218, thickness-chord ratio of 0.09 at the root and 0.07 at the tip, and 10 deg negative dihedral. Aerodynamic and longitudinal stability coefficients were obtained for the model with the tail removed, and for horizontal-tail incidences of 0 deg and -5 deg. Data were obtained at Mach numbers from 0.60 to 1.00, angles of attack from 0 deg to 12 deg, and with jet total-pressure ratios up to 3.1. Jet operation generally caused a decrease in lift, an increase in pitching-moment coefficient, and a decrease in longitudinal stability at subsonic speeds. The jet interference effects on drag were detrimental at a Mach number of 0.60 and favorable at higher speeds for cruising-flight attitudes.

Schmeer, James W.↗