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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Experimental and modeling study of C2–C4 alcohol autoignition at intermediate temperature conditions

C2–C4 alcohols are advantageous blendstocks identified by many research groups, including the U.S. Department of Energy Co-Optima Initiative, towards enabling efficient, boosted Spark-Ignition (SI) engines. Their use in advanced engine applications requires a comprehensive understanding of their intermediate-temperature autoignition behavior. This work reports an experimental and modeling study covering their fundamental autoignition characteristics in a twin-piston rapid compression machine at pressures of 20 and 40 bar, intermediate temperatures from 750 to 980 K, and two fuel loading conditions representative of boosted SI engines. Direct comparison between these alcohols is made, where the order of reactivity is established across different thermodynamic and fuel loading conditions. Changes in preliminary exothermicity (or intermediate-temperature heat release) displayed in single-stage autoignition across different alcohols and conditions are also quantified. This provides insight into fuel-to-fuel differences, and how these could affect advanced combustion concepts such as spark-assisted compression ignition. Kinetic models are used to simulate the experiments, and reasonable agreement is obtained. Furthermore, the sensitivity analysis results demonstrate the importance of accurately capturing the autoignition kinetics, particularly H-abstraction reactions on the parent fuels by OH and HO 2 , and the branching ratio associated with these.

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

The Roles of Tidal Marshes in the Estuarine Biochemical Processes: A Numerical Modeling Study

Observations suggest that the existence of tidal marsh can alter the oxygen and nutrient dynamics in adjacent water bodies, but assessing the impacts of large tidal marshes on an estuary is challenging. In this study, we use a modeling approach to investigate the roles of tidal marshes on the estuarine biochemical processes. The marsh model, which simulates the ecological functions of marshes at seasonal and annual time-scales, is embedded inside an unstructured-grid three-dimensional hydrodynamic and eutrophication model (SCHISM-ICM). This modeling system simulates the growth and metabolism of the tidal marshes and links biological processes to nutrient dynamics in the water column and sediment. This model dynamically simulates nutrient recycling and physical transport of the materials between marshes and open water through wetting-drying processes. This coupled model system is validated and successfully applied to the York River Estuary. Model results suggest that tidal marshes influence the local diurnal dissolved oxygen (DO) cycle by exporting dissolved organic carbon and high sediment oxygen demand in the marsh system through the tidal exchange. The high deposition rates of organics and diurnal DO cycle enhance the sediment release of phosphorus. On the other hand, marshes tend to decrease dissolved inorganic nitrogen in the water column by settling particulate nutrients and enhancing the denitrification process. The study demonstrates that tidal marshes exert substantial impacts on the estuarine biochemical processes. The developed tidal marsh model enhances eutrophication modeling and advances the understanding of the feedback effects between marsh biogeochemistry and estuarine eutrophication processes on a systemic scale.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The development and use of decision support framework for informing selection of select agent toxins with modelling studies to inform permissible toxin amounts

Many countries have worked diligently to establish and implement policies and processes to regulate high consequence pathogens and toxins that could have a significant public health impact if misused. In the United States, the Antiterrorism and Effective Death Penalty Act of 1996 (Public Law 104-132, 1996), as amended by the Bioterrorism Preparedness and Response Act of 2002 (Public Law 107-188, 2002) requires that the Department of Health and Human Services (HHS) [through the Centers for Disease Control and Prevention (CDC)] establish a list of bacteria, viruses, and toxins that have the potential to pose a severe threat to public health and safety. Currently, this list is reviewed and updated on a biennial basis using input from subject matter experts (SMEs). We have developed decision support framework (DSF) approaches to facilitate selection of select toxins and, where toxicity data are known, conducted modelling studies to inform selection of toxin amounts that should be excluded from select agent regulations. Exclusion limits allow laboratories to possess toxins under an established limit to support their research or teaching activities without the requirement to register with the Federal Select Agent Program. Fact sheets capturing data from a previously vetted SME workshop convened by CDC, literature review and SME input were developed to assist in evaluating toxins using the DSF approach. The output of the DSF analysis agrees with the current select toxin designations, and no other toxins evaluated in this study were recommended for inclusion on the select agent and toxin list. To inform the selection of exclusion limits, attack scenarios were developed to estimate the amount of toxin needed to impact public health. Scenarios consisted of simulated aerosol releases of a toxin in high-population-density public facilities and the introduction of a toxin into a daily consumable product supply chain. Using published inhalation and ingestion median toxic dose (TD 50 ) and median lethal dose (LD 50 ) values, where available, a range of toxin amounts was examined to estimate the number of people exposed to these amounts in these scenarios. Based on data generated by these models, we proposed toxin exclusion values corresponding to levels below those that would trigger a significant public health response (i.e., amounts estimated to expose up to ten people by inhalation or one hundred people by ingestion to LD 50 or TD 50 levels of toxin in the modeled scenarios).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

A combined experimental and modelling study of granite hydrothermal alteration

Geochemical reactions can induce significant changes of rock reservoir porosity and permeability via mineral dissolution and precipitation processes, affecting the long-term fluid behaviour within various geological systems. Here, the understanding and quantification of these reactions rely on field and experimental studies and on the predictions of reactive transport models. The present study was aimed at assessing the extent to which current geochemical models integrating available mineral dissolution/precipitation rate equations can reproduce the experimental data obtained from 4 to 17-day long hydrothermal alteration experiments of a muscovite-biotite granite and, thus, help provide an accurate description of the evolution of geothermal systems within granitic reservoirs. The experiments were conducted at a constant temperature of 180 °C and over an aqueous fluid pH range of 2 to 8.5, using both mixed-flow and static batch reactors. Modelled major element (K, Al, Si, Ca, and Mg) concentrations were generally in satisfactory agreement with the corresponding measured elemental fluxes – the differences between modelled and experimental values were generally within a factor of 5 – and the predicted identity and mass of formed secondary phases were consistent with the microscopic observations of the reacted solids. However, larger differences between measured and modelled element concentrations were observed when significant amounts of secondary phases formed, notably at pH 2 to 3, and for longer-term batch experiments. Much of this concentration difference stems from the underestimation of the amounts of Al-phases formed at acid to near-neutral pH. Although an idealized rock composition was considered, the observed mismatch between model calculations and experimental data can be attributed to inadequate mineral precipitation reaction rates and a poor description of reactive surface areas in existing geochemical modelling codes. More accurate quantification of precipitation kinetics, including nucleation and growth, and improved descriptions of the temporal change of mineral surface area would enhance the predictive capabilities of reactive transport models and benefit, particularly, the efforts aimed at increasing the sustainability of EGS reservoirs.

58 GEOSCIENCES↗
DOE OSTI2022

Fluid – Fluid Interfacial Area and Its Impact on Relative Permeability - A Pore Network Modeling Study

Abstract Relative permeability (kr) is commonly modeled as an empirical function of phase saturation. Although current empirical models can provide a good match of one or two measured relative permeabilities using saturation alone, they are unable to predict relative permeabilities well when there is hysteresis or when physical properties such as wettability change. Further, current models often result in relative permeability discontinuities that can cause convergence and accuracy problems in simulation. To overcome these problems, recent research has modeled relative permeability as a state function of both saturation (S) and phase connectivity (X). Pore network modeling (PNM) data, however, shows small differences in relative permeability for the same S-X value when approached from a different flow direction. This paper examines the impact of one additional Minkowski parameter (Mecke and Arns, 2005), the fluid-fluid interfacial area, on relative permeability to identify if that satisfactorily explains this discrepancy. We calculate the total fluid-fluid interfacial areas (IA) during two-phase (oil/water) flow in porous media using pore network modeling. The area is calculated from PNM simulations using the areas associated with corners and throats in pore elements of different shapes. The pore network is modeled after a Bentheimer sandstone, using square, triangular, and circular pore shapes. Simulations were conducted for numerous primary drainage and imbibition cycles at a constant contact angle of 0° for the wetting phase. Simultaneous measurements of capillary pressure, relative permeability, saturation, and phase connectivity are made for each displacement. Fluid-fluid interfacial area is calculated from the PNM capillary pressure, the fluid location in the pore elements, and the pore element dimensional data. The results show that differences in the relative permeability at the same (S,X) point is explained well by differences in the fluid-fluid interfacial area (IA). That is, for a larger change in IA at these intersection points, the permeability difference is greater. That difference in relative permeability approaches zero as the difference in IA approaches zero. This confirms that relative permeability can be modeled better as a unique function of S, X, and IA. The results also show that an increase in IA restricts flow decreasing the nonwetting (oil) phase permeability. This decrease is caused by an increase in the throat area fraction compared to the corner area as the total area IA increases. The wetting phase relative permeability, however, shows the inverse trend, in that its relative permeability is greater when IA becomes larger owing to a greater fraction of the total area associated with the corners. The area IA, however, impacts the nonwetting phase relative permeability more than the wetting phase relative permeability. Corner flow improves the wetting phase relative permeability because the wetting phase is continuous there. Finally, a sensitivity analysis shows that relative permeability a is more sensitive to change in S than they are for IA for the case studied implying that if only two parameters are used to model relative permeability it is better to choose S and X.

Engineering↗
DOE OSTI2021

Detection of significant antiviral drug effects on COVID-19 with reasonable sample sizes in randomized controlled trials: A modeling study

Development of an effective antiviral drug for Coronavirus Disease 2019 (COVID-19) is a global health priority. Although several candidate drugs have been identified through in vitro and in vivo models, consistent and compelling evidence from clinical studies is limited. The lack of evidence from clinical trials may stem in part from the imperfect design of the trials. We investigated how clinical trials for antivirals need to be designed, especially focusing on the sample size in randomized controlled trials. A modeling study was conducted to help understand the reasons behind inconsistent clinical trial findings and to design better clinical trials. We first analyzed longitudinal viral load data for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) without antiviral treatment by use of a within-host virus dynamics model. The fitted viral load was categorized into 3 different groups by a clustering approach. Comparison of the estimated parameters showed that the 3 distinct groups were characterized by different virus decay rates (p-value < 0.001). The mean decay rates were 1.17 d -1 (95% CI: 1.06 to 1.27 d -1 ), 0.777 d -1 (0.716 to 0.838 d -1 ), and 0.450 d -1 (0.378 to 0.522 d -1 ) for the 3 groups, respectively. Such heterogeneity in virus dynamics could be a confounding variable if it is associated with treatment allocation in compassionate use programs (i.e., observational studies). Subsequently, we mimicked randomized controlled trials of antivirals by simulation. An antiviral effect causing a 95% to 99% reduction in viral replication was added to the model. To be realistic, we assumed that randomization and treatment are initiated with some time lag after symptom onset. Using the duration of virus shedding as an outcome, the sample size to detect a statistically significant mean difference between the treatment and placebo groups (1:1 allocation) was 13,603 and 11,670 (when the antiviral effect was 95% and 99%, respectively) per group if all patients are enrolled regardless of timing of randomization. The sample size was reduced to 584 and 458 (when the antiviral effect was 95% and 99%, respectively) if only patients who are treated within 1 day of symptom onset are enrolled. We confirmed the sample size was similarly reduced when using cumulative viral load in log scale as an outcome. We used a conventional virus dynamics model, which may not fully reflect the detailed mechanisms of viral dynamics of SARS-CoV-2. The model needs to be calibrated in terms of both parameter settings and model structure, which would yield more reliable sample size calculation. In this study, we found that estimated association in observational studies can be biased due to large heterogeneity in viral dynamics among infected individuals, and statistically significant effect in randomized controlled trials may be difficult to be detected due to small sample size. The sample size can be dramatically reduced by recruiting patients immediately after developing symptoms. We believe this is the first study investigated the study design of clinical trials for antiviral treatment using the viral dynamics model.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Coupled Experimental and High-Temperature Discrete-Element Method Modeling Studies of Aluminosilicate Particle Handling in Concentrated Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrated solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrated solar power environments. Coupled experimental and numerical studies of aluminosilicate particles in rotary kilns and dense particle curtains were performed for bulk temperatures up to 1073 K. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Temperature, spatial, and velocity profile data were extracted from experimental runs using embedded K-type thermocouple probes and particle image velocimetry techniques. Experimental and numerical studies were compared using spatial temperature profiles, velocity fields, and shape profiles of the bulk, granular flows. Numerical models were developed using commercially available discrete element method modeling software, Aspherix®. Existing Aspherix® functionality was expanded by introducing coupled radiative exchange modeling tools. The laboratory-scale rotary kiln was developed to investigate the steady-state heat and mass transfer performance of aluminosilicate particles based on particle type, bulk handling temperature, and wall roughness. The rotational speed of the rotary kiln was varied to control the relative impact of collisional and frictional effects upon the granular flow behavior. Heat and mass transfer performance was categorized based on the Froude number and the observed flow regimes of slipping, rolling, cascading, and centrifuging. Coupled discrete element method modeling studies were used to evaluate the effects of temperature-dependent, particulate mechanical properties upon bulk flow behavior and upon the relative effects of radiative, advective, and/or conductive heat transfer. A high-temperature (< 1073 K) falling particle curtain was similarly fabricated to investigate the heat and mass transfer performance of aluminosilicate particles in particle handling situations dominated by inter-particle collisions. The impact of particle type, flow preheat temperatures (< 1073K), and bulk mass flow rates were investigated upon the particle curtain shape, temperature, and velocity profiles. Coupled discrete element method modeling studies were performed to evaluate the varying impact of temperature-dependent, particulate mechanical properties on the bulk flow behavior and the temperature profile of the particle curtain.

14 SOLAR ENERGY↗
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 OSTI2021

Impact of Lateral Flow on Surface Water and Energy Budgets over the Southern Great Plains – A Modeling Study

As we see the horizontal grid spacing decrease, treatment of hydrologic processes in land surface models, such as the lateral flow of surface and subsurface flow, need to be explicitly represented. Unlike previous studies that mainly focused on the mountainous regions, in this study the offline WRF-Hydro model is employed to study the impact of lateral flow on soil moisture and energy fluxes over the relatively flat southern Great Plains (SGP). The vast amount of measurements over the SGP provide an unique opportunity to assess the model behavior. In addition, newly developed land surface properties and input forcing are ingested into the model, in an attempt to reduce uncertainties associated with the initial and boundary forcing and help to identify model deficiencies. Our results show that the more realistic inputs (parameters, soil types, forcing) lead to larger underestimation of latent heat flux and dry bias, indicating the existence of model structural uncertainty (embedded errors) in WRF-Hydro that need to be characterized to inform future model development efforts. Including lateral flow processes partly mitigates the model deficiencies in representing hydrologic processes and alleviates the dry bias. In particular, both surface and subsurface lateral flow increase soil moisture mainly over the lower elevations, except that subsurface flow also affects soil moisture over steeper terrains. Additional simulations are performed to assess the effect of routing resolution on model results. When LSM resolution is high, noticeable differences in soil moisture are produced between different routing resolutions especially over steep terrain. Whereas when LSM resolution is coarse, differences between routing resolutions become negligible, especially over flat terrain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

High-Resolution Frequency-Domain Spectroscopic and Modeling Studies of Photosystem I (PSI), PSI Mutants and PSI Supercomplexes

Photosystem I (PSI) is one of the two main pigment–protein complexes where the primary steps of oxygenic photosynthesis take place. This review describes low-temperature frequency-domain experiments (absorption, emission, circular dichroism, resonant and non-resonant hole-burned spectra) and modeling efforts reported for PSI in recent years. In particular, we focus on the spectral hole-burning studies, which are not as common in photosynthesis research as the time-domain spectroscopies. Experimental and modeling data obtained for trimeric cyanobacterial Photosystem I (PSI 3 ), PSI 3 mutants, and PSI 3 –IsiA 18 supercomplexes are analyzed to provide a more comprehensive understanding of their excitonic structure and excitation energy transfer (EET) processes. Detailed information on the excitonic structure of photosynthetic complexes is essential to determine the structure–function relationship. We will focus on the so-called “red antenna states” of cyanobacterial PSI, as these states play an important role in photochemical processes and EET pathways. The high-resolution data and modeling studies presented here provide additional information on the energetics of the lowest energy states and their chlorophyll (Chl) compositions, as well as the EET pathways and how they are altered by mutations. We present evidence that the low-energy traps observed in PSI are excitonically coupled states with significant charge-transfer (CT) character. The analysis presented for various optical spectra of PSI 3 and PSI 3 -IsiA 18 supercomplexes allowed us to make inferences about EET from the IsiA 18 ring to the PSI 3 core and demonstrate that the number of entry points varies between sample preparations studied by different groups. In our most recent samples, there most likely are three entry points for EET from the IsiA 18 ring per the PSI core monomer, with two of these entry points likely being located next to each other. Therefore, there are nine entry points from the IsiA 18 ring to the PSI 3 trimer. We anticipate that the data discussed below will stimulate further research in this area, providing even more insight into the structure-based models of these important cyanobacterial photosystems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Experimental and modeling studies of the plasma chemistry in a humid Ar radiofrequency atmospheric pressure plasma jet

Abstract While humid atmospheric pressure plasmas are extensively modeled, reaction set validation for these conditions remains limited. We present a detailed comparison of a modelling and experimental study with a goal to elucidate the plasma chemistry in a humid Ar radiofrequency-driven atmospheric pressure plasma jet. A large group of species including radicals (H, OH, O, HO 2 ) and long-lived species (H 2 , O 2 and H 2 O 2 ) in the jet effluent was experimentally quantified by molecular beam mass spectroscopy (MBMS). MBMS measurements of H 2 O 2 , OH and H were validated by direct comparison with a liquid phase colorimetric measurement, laser-induced fluorescence (LIF) and two-photon absorption LIF respectively. While an excellent agreement was found for OH and H 2 O 2 by both techniques, a significant difference was found for H and shown to be due to boundary layer effects at the MBMS sampling substrate. The measured O, OH, HO 2 and H 2 are in good agreement with the plug model while H and O 2 were underestimated and H 2 O 2 was overestimated by the model. The accuracy of both the used reaction set and the diagnostics, as well as the observed discrepancies between the modeling and experimental results, are critically assessed. The results presented in this work enable us to identify further data needs for describing H 2 O vapor chemistry in low-temperature plasmas.

Physics↗
DOE OSTI2020

Modeling study of deep direct use geothermal on the West Virginia university campus-morgantown, WV

To reduce the geothermal exploration risk, a feasibility study is performed for a deep direct-use (DDU) system proposed at the West Virginia University (WVU) Morgantown campus. This study applies numerical simulations to investigate reservoir impedance and thermal production. Because of the great depth of the geothermal reservoir, few data are available to characterize reservoir features and properties. Consequently, the study focuses on the following three aspects: 1. model choice for predicting reservoir impedance and thermal breakthrough: after investigating three potential models (one single permeability model and two dual permeability models) for flow through fractured rock, it is decided only the single permeability model is needed; 2. well placement (horizontal vs. vertical) options: horizontal well placement seems to be more robust to heterogeneity and the impedance is more acceptable; 3. Prediction uncertainty: the most influential parameters are identified using a First-Order-Second-Moment uncertainty propagation analysis, and the uncertain range of the model predictions is estimated by performing a Monte Carlo simulation. Heterogeneity has a large impact on the prediction, therefore, heterogeneity is included in the predictive model and uncertainty analysis. The numerical model results and uncertainty analysis will be used for further economic studies.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Temperature Effect on Photoelectrochemical Water Splitting: A Model Study Based on BiVO 4 Photoanodes

Photoelectrochemical (PEC) water splitting is typically studied at room temperature. In this work, the temperature effect on PEC water splitting is studied using crystalline BiVO4 thin film photoanode as a model system. Systematic temperature-dependent electrochemical study demonstrates that the PEC activity is boosted at elevated electrolyte temperatures and indicates that thermal energy plays a main role in improving charge carrier transport in the bulk of BiVO 4 . Irreversible surface reconstruction is observed after PEC reactions at elevated temperature in the presence of hole scavengers, with regularly spaced stripes emerging on BiVO4 grains. The surface-reconstructed photoanode exhibits up to 40% improvement in photocurrent densities and ~0.25 V shift of photocurrent onset to favorable direction. Detailed investigation shows the formation of amorphous layer without stoichiometric change at the reconstructed surface. Furthermore, this work provides insights of temperature effect on photoelectrode in solar water splitting and reveals the non-negligible effect of hole scavengers in photoelectrochemical measurement.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental and kinetic modeling study of tetralin: A naphtheno-aromatic fuel for gasoline, jet and diesel surrogates

Distillate fuels contain significant proportions of naphtheno-aromatic components and tetralin is a suitable surrogate component to represent this molecular moiety. The presence of aromatic and naphthyl rings makes kinetic modeling of tetralin very challenging. Primary radicals formed during the oxidation of tetralin can be aryl, benzylic or paraffinic in nature. Using available information on reaction paths and rate constants of naphthenes and alkyl-aromatics, in this study a kinetic model of tetralin has been developed with emphasis on low-temperature chemistry and high-pressure conditions. Due to the lack of high-level quantum chemical calculations on reaction pathways of tetralin, analogous rates from ab-initio studies on benzylic and paraffinic radicals have been adopted here. Some modifications to the reaction rate rules are incorporated to account for the unique characteristics of tetralin's molecular structure. Important reaction channels have been identified using reaction path and brute force sensitivity analyses. In order to investigate the model performance at low temperatures, new experiments are carried out in a rapid compression machine on blends of tetralin and 3-methylpentane. Blending of low-reactivity tetralin with a high-reactivity alkane allowed the investigation of tetralin ignition at very low temperatures (665 – 856 K). The kinetic model developed in the current study is found to predict the current experiments and literature data adequately. The new model will aid in high-fidelity surrogate predictions at engine-relevant conditions.

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

Shell-model study of calcium isotopes toward their drip line

Here, we report in this paper a study in terms of the nuclear shell model about the location of the calcium isotopes drip line. The starting point is considering the realistic two-body potential derived by Entem and Machleidt within chiral perturbation theory at next-to-next-to-next-to-leading order (N 3 LO), as well as a chiral three-body force at next-to-next-to-leading order (N 2 LO) whose structure and low-energy constants are consistent with the two-body potential. Then we construct the effective single-particle energies and residual interaction needed to diagonalize the shell-model Hamiltonian. The calculated two-neutron separation energies agree nicely with experiment until 56 Ca, which is the heaviest isotope whose mass has been measured, and do not show any sign of two-neutron emission until 70 Ca . We discuss the role of the choice of the model space in determining the neutron drip line, and also the dependence of the results on the parameters of the shell-model Hamiltonian.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

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

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

09 BIOMASS FUELS↗
DOE OSTI2022

An experimental and kinetic modeling study of the pyrolysis of isoprene, a significant biogenic hydrocarbon in naturally occurring vegetation fires

Isoprene dominates the carbon flux emitted by vegetation and constitutes 40% of non-methane biogenic emissions worldwide. Despite pyrolysis experiments at temperatures above 1000 K showing a link between isoprene combustion and aromatic species formation, comprehensive mechanistic research on isoprene is scarce in the literature. Here, in this work, we carry out an experimental and theoretical study to build, for the first time, a chemical kinetic model describing isoprene pyrolysis. The formation of polycyclic aromatic hydrocarbon (PAH) precursor species, often observed in vegetation fire plumes, is partially explained by isoprene pyrolysis experiments and theoretical modeling. Molecular dynamics (MD) simulations unveil reaction pathways from allylic isoprenyl radicals to allene and cyclopentadiene (CPD) intermediates, two relevant species detected in the experiments. Rate constants for these identified pathways are calculated using variational transition state theory to update the kinetic model, which is validated against single-pulse shock tube (SPST), and jet-stirred reactor (JSR) experimental data in the temperature range of 850–1690 K. The kinetic model presents satisfactory agreement with the SPST experimental data, and a reaction pathway analysis shows that association of propargyl radicals results in benzene formation. The JSR pathway analysis also identifies the prominent reactions for CPD, benzene, styrene, and toluene formation. Our model does not reproduce the CPD experimental profiles, indicating that additional studies are necessary. Overall, our findings advance the understanding of isoprene pyrolysis and its related atmospheric pollutants in naturally occurring vegetation fires where smoldering and oxygen-deficient combustion processes are present.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗