Search NASA⌕ Search

SEARCH · Search NASA

Results for “modelling studies”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

An experimental and chemical kinetic modeling study of 4-butoxyheptane combustion

Here, the combustion kinetics of a novel oxygenated bioblendstock for diesel, 4-butoxyheptane (4-BH), was investigated experimentally using a flow reactor and a heated, high-pressure shock tube. The flow reactor experiments employed oxygen as the oxidizer and helium as the diluent with oxidation conducted at atmospheric pressure and 10 bar for temperatures from 400 to 1000 K at 20-K intervals. The fuel, oxidizer, and diluent flow rates were varied at different temperatures to maintain a constant initial fuel mole fraction of 1000 ppm, with stoichiometric equivalence ratio, and a residence time of 2.0 s. The reacted gas was fed to two separate GC systems that could qualitatively and quantitatively detect product species. Additionally, real fuel-air ignition delay time (IDT) data were collected using a heated, high-pressure shock-tube facility. Fuel lean (φ = 0.5) and stoichiometric (φ = 1.0) mixtures were investigated at 10 atm as well as at 30 atm for the fuel lean case for temperatures between 847 and 1259 K. A detailed chemical kinetics mechanism was developed to model the product distribution from the flow reactor and IDTs from the shock tube. The proposed model was able to predict the double NTC behavior in flow reactor experiments reasonably well. Model predictions at low temperatures were observed to be highly sensitive to the rate constants of ketohydroperoxide (KHP) decomposition in the case of the OOH group in α position which were modeled based on existing literature studies on ethers. It was noted that in the absence of theoretical or experimental studies, the rate constants for KHP decomposition used in the literature were empirically set. Additional studies are required to address the gap in model prediction obtained in this study and to reduce the uncertainty in kinetics models for ether oxidation. Predicted product concentrations and IDTs showed some quantitative agreement with experimental data, but the overall reactivity of the IDTs is underpredicted. Additionally, significant deviation is observed for the IDT results at 10 atm for the stoichiometric case with minor deviations for the other cases. The reaction pathways to the missing products were then further analyzed theoretically through quantum-mechanical calculations.

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

Effects of groundwater pumping on ground surface temperature: A regional modeling study in the North China Plain

Over-exploitation of groundwater (GW) in the North China Plain (NCP) since the 1960s has many environmental consequences. However, understanding of the dominant mechanisms remains limited, particularly at regional scale. In this study, the coupled ParFlow.CLM model representing subsurface and land-surface processes and their interactions was applied in the NCP at high spatio-temporal resolutions. The model was validated using the water and energy fluxes reported in previous studies and from the JRA-55 reanalysis. Numerical experiments were designed to examine the relative impacts of GW pumping and irrigation on the ground surface temperature (GST). Results show significant effects of GW pumping on GST in the NCP. Generally, the subsurface acts as a buffer to temporal variations in heat fluxes at the land-surface, but long-term pumping can gradually weaken this buffer, resulting in increases in the spatio-temporal variability of GST, as exemplified by hotter summers and colder winters. Considering that changes of water table depth (WTD) can significantly affect land surface heat fluxes when WTD ranges roughly between 0.01–10 m, the 0.5 m/year increase of WTD simulated by the model due to pumping can continue to increase the average WTD and hence, disturb GST, for about 20 years from the pre-pumping WTD in the NCP, before the WTD exceeds 10 m. The variations of GST are expected to be faster initially and gradually slow down due to the nonlinear behaviors of GST with WTD. The findings from this study in the NCP may also have implications for other regions with GW depletion.

Yang, Chen↗
DOE OSTI2023

Isomeric effects on the reactivity of branched alkenes: An experimental and kinetic modeling study of methylbutenes

Here, a detailed experimental study of the low-to-intermediate temperature combustion of methylbutene isomers, i.e., branched C 5 alkenes, has been undertaken with multiple experimental facilities. Ignition delay times were measured at equivalence ratios 0.5–2.0, 685–1020 K and up to 45 bar condition from two rapid compression machines and showed slight deviation from an Arrhenius behavior for all three isomers, while their reactivity order differs as temperature changes. Sampled intermediates formed during the oxidation process of mixtures at 900–1150 K and 0.82 bar from a flow reactor and at 730 K and 20 bar from a rapid compression machine were analyzed using gas chromatography techniques. Trends in the formation and consumption of sampled intermediates were modeled using a kinetic model developed in this work for all three isomers. Rate of production and sensitivity analyses emphasize the role of double bond-specific reactions governing the global reactivity of these fuels. Additional studies of the addition reactions of HO 2 radicals to the double bond and to allylic radicals may improve the model performance.

2-Methyl-1-butene↗
DOE OSTI2022

Multistage Surface-Heated Vacuum Membrane Distillation Process Enables High Water Recovery and Excellent Heat Utilization: A Modeling Study

Surface-heated membrane distillation (MD) enhances the energy efficiency of desalination by mitigating temperature polarization (TP). However, systematic investigations of larger scale, multistage, surface-heated MD system with high water recovery and heat recycling are limited. Here, we explore the design and performance of a multistage surface-heated vacuum MD (SHVMD) with heat recovery through a comprehensive finite difference model. In this process, the latent heat of condensation is recovered through an internal heat exchanger (HX) using the retentate from one stage as the condensing fluid for the next stage and an external HX using the feed as the condensing fluid. Model results show that surface heating enhances the performance compared to conventional vacuum MD (VMD). Specifically, in a six-stage SHVMD process, 54.44% water recovery and a gained output ratio (GOR) of 3.28 are achieved with a surface heat density of 2000 W m –2 , whereas a similar six-stage VMD process only reaches 18.19% water recovery and a GOR of 2.15. Mass and energy balances suggest that by mitigating TP, surface heating increases the latent heat trapped in vapor. The internal and external HXs capture and reuse the additional heat, which enhances the GOR values. We show for SHVMD that the hybrid internal/external heat recovery design can have GOR value 1.44 times higher than that of systems with only internal or external heat recovery. Furthermore, by only increasing six stages to eight stages, a GOR value as high as 4.35 is achieved. The results further show that surface heating can reduce the energy consumption of MD for brine concentration. As a result, the multistage SHVMD technology exhibits a promising potential for the management of brine from industrial plants.

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

Styrene Thermal Decomposition and Its Reaction with Acetylene under Shock Tube Pyrolysis Conditions: an Experimental and Kinetic Modeling Study

Styrene is an important compound for polymer production and a key intermediate in gas-phase kinetics of polycyclic aromatic hydrocarbons (PAHs). For the first time, the pyrolysis of styrene with and without the presence of acetylene is investigated in a single-pulse shock tube coupled to gas chromatography and mass spectrometry. For each reaction system, quantitative speciation profiles are probed within the temperature range of 1100-1730 K, nominal pressure of 20 bar, and reaction duration of similar to 4 ms. A kinetic model is built to simulate the results. The model explains how styrene is consumed under high-pressure pyrolytic conditions, how the secondary chemistry of intermediate products affect subsequent PAH formation, and how acetylene addition alters the reaction pathways. Throughout the temperature range, styrene breakdown is dominated by the bimolecular interaction between styrene and hydrogen atom, which produces benzene+vinyl or phenyl and ethylene through the stabilization of 2-phenylethyl and its subsequent dissociation. As a result, large amounts of phenyl accumulate, which react with styrene to form C14H12 species while simultaneously releasing H atoms through addition/elimination reactions. The reactivity of fuel consumption is preserved by the regeneration of H atoms as chain carriers. Several C14H10 compounds are formed as a result of the following breakdown of the C14H12 isomers, particularly stilbene, 1,1-diphenyl ethylene, and 9-methyl-9H-fluorene. The presence of acetylene as a co-reactant with styrene allows the Hydrogen-Abstraction-Acetylene-Addition (HACA) pathway to proceed from phenyl radical to enhance the production of phenylacetylene at very low temperatures and acenaphthylene. This hinders the formation of C14H12 isomers, exclusive products from pure styrene dissociation by competing with the styrene+phenyl routes.

Kinetic Modeling↗
DOE OSTI2024

Generating synthetic signaling networks for in silico modeling studies

Predictive models of signaling pathways have proven to be difficult to develop. Reasons include the uncertainty in the number of species, the complexity in species’ interactions, and the sparseness and uncertainty in experimental data. Traditional approaches to developing mechanistic models rely on collecting experimental data and fitting a single model to that data. This approach works for simple systems but has proven unreliable for complex systems such as biological signaling networks. For example, uncertainty and sparseness of the data often result in overfitted models that have little predictive value beyond recapitulating the experimental data itself. Thus, there is a need to develop new approaches to create predictive mechanistic models of complex systems. However, to determine the effectiveness of any new algorithm, a baseline model is needed to test its performance. To meet this need, we developed a method for generating artificial synthetic networks that are reasonably realistic and thus can be treated as ground truth models. These synthetic models can then be used to generate synthetic data for developing and testing algorithms designed to recover the underlying network topology and associated parameters. Here, we describe a simple approach for generating synthetic signaling networks that can be used for this purpose.

42 ENGINEERING↗
DOE OSTI2022

Flowability of Crumbler Rotary Shear Size-Reduced Granular Biomass: An Experiment-Informed Modeling Study on the Angle of Repose

Biomass has potential as a carbon-neutral alternative to petroleum for chemical and energy products. However, complete replacement of fossil fuel is contingent upon efficient processes to eliminate undesirable characteristics of biomass, e.g., low bulk density, variability, and storage-induced quality problems. Mechanical size reduction via comminution is a processing operation to engineer favorable biomass flowability in handling. Crumbler rotary shear mill has been empirically demonstrated to produce more uniformly shaped particles with higher flowability than hammermilled biomass. This study combines modeling and experimentation to unveil fundamental understandings of the relation between granular particle characteristics and biomass flow behavior, which elucidate underlying mechanisms and guide selection of critical processing parameters. For this purpose, the impact of critical material attributes, including particle size (2–6 mm), particle shape (briquette, chip, clumped-sphere, cube, etc.), and surface roughness, on the angle of repose (AOR) of milled pine chips were investigated using discrete element method (DEM) simulations. Forest Concepts Crumbler rotary shear system is used to produce milled pine particles within the same size range considered in DEM simulations. AOR of different sets of these particles were measured experimentally to benchmark DEM results against experimental data. Specific energy consumption for the comminution of biomass with different particle size and moisture content are measured for technoeconomic analysis. Our results show that the smaller size (2 mm) of pine particle achieves better followability (i.e., smaller AOR) while the energy cost of comminution is significantly higher and bulk density is almost the same as the 6-mm pine particles. For the 2-mm particle size, Crumbles from veneer have better flow properties than Crumbles from chips. Contrarily, no significant difference was observed between the AOR of the two materials for the 6-mm particle size. Furthermore, from DEM simulations, mechanical interlocking between particles was found as a dominant factor in determining AOR of complex-shaped particles such as milled pine, which cannot be accurately captured by using simple particle shapes (e.g., mono-sphere) with a rolling resistance model. Conversely, clumped-sphere model alleviates this limitation without increasing computational cost significantly and can be used for accurate representation of biomass granular particles when simulating free-flow behavior.

09 BIOMASS FUELS↗
DOE OSTI2022

Investigating the kinetic effect of prenol on iso -octane auto-ignition by means of an experimental and modeling study

A detailed experimental and kinetic modeling study was dedicated to understand the reported octane hyperboosting effect of prenol, by means of the measurement of the ignition delay times of its blends with iso -octane, and measurement of the mole fraction profiles of the fuels and intermediates inside the ULille rapid compression machine. Here, these results show that prenol addition leads to a reduction of the first-stage ignition phenomena and negative temperature coefficient behavior, which is only qualitatively captured by the model and is consistent with knock resistance improvement. It is suggested that this behavior is caused by two different factors. The first originates from gas-phase reactivity of prenol, and spans from the formation of unreactive unsaturated species through resonance-stabilized radicals, thereby constituting a competitive pathway for the radical pool generated by iso -octane. The second is of catalytic nature and cannot be captured by means of gas-phase kinetic modeling, but could also play an important role in the behavior of prenol in internal combustion engines.

42 ENGINEERING↗
DOE OSTI2020

Electro-reduction of organics on metal cathodes: A multiscale-modeling study of benzaldehyde on Au (111)

Here, we present a multiscale modeling study on the electro-reduction of benzaldehyde (BZY) on Au(111) under realistic electrochemical conditions. To model the electrochemical cell, we adopt a capacitor model in which complex solvents are confined between a Au cathode and a carbon anode. Classical molecular dynamics simulations reveal that electrode charge density and the presence of alcohol show strong effects on the density, adsorption geometry and dynamics of benzaldehyde on the Au electrode. Under charging conditions, the surface concentration of benzaldehyde on the Au electrode decreases, while the content of other species increases. Finally, we proposed a scheme that correlates the electric current running through the Au/solvent interface with the applied bias. This study provides a molecular level understanding of how solvent composition, in this case water/alcohol content, controls the activity of electrocatalytic hydrogenation

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

A Model for Hybrid Systems for Production Cost Modeling Studies Considering Ancillary Services

This paper introduces a model for simulating hybrid plants participating in energy and ancillary services for bulk power system studies. The model considers a hybrid plant comprised of a renewable energy source, a thermal power unit, a storage unit, fixed power loads, or any combination of these technologies. The model focuses on Production Cost Modeling (PCM) studies under the assumption of centralized dispatch. We present an example case study to illustrate the use of the model in a single-stage production cost model similar to those conducted by planning agencies. We explore the allocation of behind-the-meter ancillary services products and total energy participation to hybrid plant sub-assets and the resulting impacts on the system's ancillary service allocation. The model is implemented and simulated in a unit commitment problem in the RTS test system, which was modified to include a hybrid plant asset.

ancillary services↗
DOE OSTI2024

A Model for Hybrid Systems for Production Cost Modeling Studies Considering Ancillary Services: Preprint

This paper introduces a model for simulating hybrid plants participating in energy and ancillary services for bulk power system studies. The model considers a hybrid plant comprised of a renewable energy source, a thermal power unit, a storage unit, fixed power loads, or any combination of these technologies. The model focuses on Production Cost Modeling (PCM) studies under the assumption of centralized dispatch. We present an example case study to illustrate the use of the model in a single-stage production cost model similar to those conducted by planning agencies. We explore the allocation of behind-themeter ancillary services products and total energy participation to hybrid plant sub-assets and the resulting impacts on the system's ancillary service allocation. The model is implemented and simulated in a unit commitment problem in the RTS test system, which was modified to include a hybrid plant asset.

ancillary services↗
DOE OSTI2026

Theoretical and kinetic modeling study of hydrazine oxidation

The present work constitutes the first theoretical and kinetic modeling study of hydrazine oxidation, which may be important for burnout in ammonia-fueled combustion. The kinetics of the oxidation of N 2 H 4 , N 2 H 3 and tHNNH by molecular oxygen were investigated via a quantum chemistry/canonical transition state theory approach. Geometries and anharmonic frequencies were obtained with density functional theory, and energies from coupled cluster calculations (CCSD(T)) extrapolated to the infinite basis set limit, with corrections for core-valence electron correlation, scalar relativistic effects, and higher level correlation up to lambda coupled cluster, CCSDT(Q) Λ . The key reactions occurred on the N 2 H 4 O 2 potential energy surface, where the results indicated a fast reaction of N 2 H 3 with HO 2 via singlet adducts to yield tHNNH + H 2 O 2 and HNN(H)O + H 2 O, while reaction on the triplet surface proceeds via a bound complex followed by a tight, submerged barrier to yield N 2 H 4 + O 2 . The results were incorporated in a detailed reaction mechanism, which was used to interpret the shock tube results from Michel and Wagner (1965) on the effect of O 2 on hydrazine conversion at 1100–1400 K. The kinetic model captured qualitatively the observed behavior, but underestimated the reaction rate under oxidizing conditions. The hydrazine pyrolysis chemistry dominated conversion at reducing conditions and/or high temperature. At oxidizing conditions and intermediate temperatures (≲ 1400 K), reactions of N 2 -amines with HO 2 and O 2 were important for the oxidation rate.

Ab initio calculations↗
DOE OSTI2021

An experimental and kinetic modeling study of NO x sensitization on methane autoignition and oxidation

An experimental and kinetic modeling study of the influence of NO x (i.e. NO 2 , NO and N 2 O) addition on the ignition behavior of methane/‘air’ mixtures is performed. Ignition delay time measurements are taken in a rapid compression machine (RCM) and in a shock tube (ST) at temperatures and pressures ranging from 900–1500 K and 1.5–3.0 MPa, respectively for equivalence ratios of 0.5–2.0 in ‘air’. The conditions chosen are relevant to spark ignition and homogeneous charge compression ignition engine operating conditions where exhaust gas recirculation can potentially add NO x to the premixed charge. The RCM measurements show that the addition of 200 ppm NO 2 to the stoichiometric CH 4 /oxidizer mixture results in a factor of three increase in reactivity compared to the baseline case without NOx for temperatures in the range 600–1000 K. However, adding up to 1000 ppm N 2 O does not show any appreciable effect on the measurements. The promoting effect of NO 2 was found to increase with temperature in the range 950–1150 K, while the sensitization effect decreases at higher pressures. The experimental results measured are simulated using NUIGMech1.2 comprising an updated NO x sub-chemistry in this work. A kinetic analysis indicates that the competition between the reactions ĊH 3 + NO 2 ↔ CH 3 Ȯ + NO and ĊH 3 + NO 2 (+M) ↔ CH 3 NO 2 (+M), the former being a propagation reaction and the latter being a termination reaction governs NO x sensitization on CH 4 ignition. Recent calculations by Matsugi and Shiina (A. Matsugi, H. Shiina, J. Phys. Chem. A. 121 (2017) 4218–4224) for the nitromethane formation reaction CH 3 + NO 2 (+M) ↔ CH 3 NO 2 (+M), together with the recently calculated rate constants for HONO/HNO 2 reactions significantly improve ignition delay time predictions in the temperature range 600–1000 K. Furthermore, the experiments with NO addition reveal a non-monotonous sensitization impact on CH 4 ignition at lower temperatures with NO initially acting as an inhibitor at low NO concentrations and then as a promoter as NO concentrations increase in the mixture. This non-monotonous trend is attributed to the role of the chain-termination reaction ĊH 3 + NO 2 (+M) ↔ CH 3 NO 2 (+M) and the impact of NO on the transition to the chain-branching steps CH 2 O + HȮ 2 ↔ HĊO + H 2 O 2 , H 2 O 2 (+M) ↔ ȮH + ȮH (+M), HĊO ↔ CO + Ḣ followed by CO + O 2 ↔ CO 2 + Ö and Ḣ + O 2 ↔ Ö + ȮH. NUIGMech1.2 is systematically validated against the new ignition delay measurements taken here together with species measurements and high temperature ignition delay time data available in the literature for CH 4 /oxidizer mixtures diluted with NO 2 /N 2 O/NO and is observed to accurately capture the sensitization trends.

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

Pyrocumulonimbus Events Over British Columbia in 2017: An Ensemble Model Study of Parameter Sensitivities and Climate Impacts of Wildfire Smoke in the Stratosphere

Abstract The Pyrocumulonimbus (pyroCb) events over British Columbia in 2017 were observed in the lower stratosphere for about 8–10 months after the smoke injections. Several previous studies used global climate models to investigate the physical parameters for the 2017 pyroCb events, but the conclusions show strong model dependency. In this study, we use the Energy Exascale Earth System Model (E3SM) and complete an ensemble of runs exploring three injection parameters: smoke aerosol mass, the percentage of black carbon within the smoke aerosols, and plume injection height. Additionally, we consider the heterogeneous reaction of ozone and primary organic matter. According to the satellite daily observed aerosol optical depth (AOD), we find that the best ensemble member is the simulation with 0.4 Tg of smoke, 3% of which is black carbon, a 13.5 km smoke injection height, and a 10 −5 probability factor of the heterogeneous reaction. Besides AOD, we examine the ensemble score based on the metrics used in previous studies: the metrics of the extinction coefficient at 18 km altitude and the maximum plume height. The conclusion of the best estimate of the injection parameters for 2017 pyroCb events shows strong not only model but also evaluation metric dependency. We use the Random Forest machine learning technique to quantify the relative importance of each parameter in accurately simulating the 2017 pyroCb events and find that the injection height is the most critical feature, no matter which metric is used to score the ensemble members.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

An experimental and detailed kinetic modeling study of the pyrolysis and oxidation of allene and propyne over a wide range of conditions

Allene and propyne are important intermediates in the pyrolysis and oxidation of higher hydrocarbon fuels, and they are also a major source of propargyl radical formation, which can recombine into different C 6 H 6 isomers and finally produce soot. In a prior work (Panigrahy et al., “A comprehensive experimental and improved kinetic modeling study on the pyrolysis and oxidation of propyne”, Proc. Combust. Inst 38 (2021)), the pyrolysis, ignition , and laminar flame speed of propyne were investigated. To understand the kinetic features of initial fuel breakdown and oxidation of the two C 3 H 4 isomers, new measurements for allene pyrolysis and oxidation are conducted in the present paper at the same operating conditions as those studied previously for propyne. Ignition delay times of allene are measured using a high-pressure shock tube and a heated twin-opposed piston rapid compression machine in the temperature range 690– 1450 K at equivalence ratios of 0.5, 1.0 and 2.0 in ‘air’, and at pressures of 10 and 30 bar. Pyrolysis species measurements of allene and propyne are also performed using a gas chromatography integrated singlepulse shock tube in the temperature range 1000–1700 K at pressure of 2 and 5 bar. Furthermore, laminar flame speeds of allene are measured at elevated gas temperatures of 373 K at pressures of 1 and 2 bar for a wide range of equivalence ratios from 0.6 to 1.5. A newly updated kinetic mechanism developed for this study is the first model that can well reproduce all of the experimental results for both allene and propyne. It is observed that in the pyrolysis process, allene dissociates faster than propyne. Both isomers exhibit similar ignition delay times at high temperatures (>1000 K), while, at intermediate temperatures (770–1000 K) propyne is the faster to ignite, and at lower temperatures (< 770 K) allene becomes more reactive. Furthermore, laminar flame speeds for propyne are found to be slightly faster than those for allene under the conditions studied in this work

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

EXPERIMENTAL AND MODELING STUDY OF MOLECULAR TRANSPORT PHENOMENA IN COMPLEX LIQUID MIXTURES USING COVALENT ORGANIC FRAMEWORK MEMBRANES

Currently, 10-15% of the world’s energy is consumed by chemical separations, and more than 80% of that is used to purify organic liquids and recover rare earth elements and critical minerals. Membrane nanofiltration presents an energy-efficient, cost-effective, and eco-friendly alternative to current separation technologies. These complex liquid environments require high- performance, chemically-resistant membrane materials. Recently discovered Covalent Organic Frameworks (COFs) possess desirable properties as membrane materials for complex liquid separations. COFs are highly crystalline, chemically and thermally stable, with tunable size and charge. COFs, especially two-dimensional highly crystalline COFs, possess narrow pore-size distribution and controlled porous structures. Two-dimensional COFs are also resistant to swelling, another feature beneficial in complex liquid separations. Molecular interactions in complex liquid systems often dictate final membrane performance, resulting in a discrepancy between designed and apparent COF properties. Understanding and resolving this discrepancy is critical in utilizing COF membranes as a viable alternative to energy-intensive separations. The primary objective of this thesis was to use a commercially available COF TpPa-1 as a platform to investigate how mixed solvent environment affects COF membrane performance via experimental and molecular modeling studies. Specifically, the effects of solution pH, solvent- solvent-solute interactions in mixed solvents, and COF chemistry were carefully examined on apparent TpPa-1 pore size and permeability and targeted solute rejection of COF membranes for neat and mixed solvents. Experimental COF filtration performance was compared and corroborated with modeling predication. Findings from this study will provide guidance for future COF design, synthesis, and desired functional COF membrane performance.

Barnes, Anastasia↗