Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model initialization”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗
DOE OSTI2023

Numerical simulations of the Los Alamos gapstick experiment

The gapstick is a high explosive (HE) sensitivity test recently developed at Los Alamos. The experiment design was motivated by the traditional gap test and consists of a series of HE and inert pellets in a rate stick configuration. The inert pellets are made increasingly longer so that eventually the attenuated shock is unable to initiate the next HE pellet. While the gapstick is a mechanically simple experiment it poses several challenges for numerical simulation. In particular, accurately modeling the HE initiation and detonation phenomena requires a reactive burn model with sufficient mesh resolution to capture the reaction scales. In this work, the Scaled Uniform Reactive Front (SURF) and Arrhenius Wescott-Stewart-Davis (AWSD) burn models are used for simulations of a PBX 9501 (95 wt% HMX, 5% binder) gapstick. Furthermore, the material model for the inert pellet material, Polyvinylidene Fluoride (PVDF), has a direct influence on the shock propagation and a new equation of state (EOS) calibration for PVDF is developed using available Hugoniot data. Despite some challenges, the simulations are able to reasonably predict shock transit velocities and detonation failure in the gapstick.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Technical Report on Waveform Fit Metrics for Global Models

The new WAVEFORMS Initiative in the Ground-based Nuclear Detonation Detection (GNDD) program includes an increased emphasis on the development of Earth models and methods to predict entire seismic and acoustic waveforms more accurately. In general, this increased emphasis is predicated on the need to better characterize seismic events and provide improved model-based discrimination between event types including earthquakes and explosions. More specifically, while current moment tensor inversion methods tend to work well for larger events (M>~4) using tuned 1-D Earth models, the development of state-of-the-art 3-D models and methods is required for the prediction of shorter period waves over large areas for discrimination of smaller events. There is no standard metric for model-based waveform prediction accuracy used in the waveform modeling/inversion community. However, there are several popular waveform misfit definitions; and minimizing the corresponding objective functions is the goal of waveform inversion. Some example misfit definitions employed for adjoint waveform tomography include measures of simple travel time differences (e.g. Tape et al., 2010), cross-correlation travel time differences (e.g. Luo and Schuster, 1991), multi-taper frequency dependent methods (e.g. Lei et al., 2020), time-frequency phase misfit functions (e.g. Fichtner 2010; Rodgers et al., 2022), normalized cross-correlation methods (e.g. Tao et al., 2018), and others. In some cases, these misfit definitions also involve complicated weighting schemes and summations over multiple frequency bands making it difficult to duplicate the misfit measurement with alternative models and datasets. Although each of the misfit definitions mentioned above are useful for developing waveform models, the actual misfit values are not usually meaningful outside of a given project, model, and/or dataset. Therefore, it is difficult to understand and communicate model performance for predicting waveforms and comparing to other models and/or new model iterations with a different dataset. Therefore, there is a need for a generalized method for evaluating overall model performance that is independent from the specific misfit chosen to develop the waveform models that is also intuitive and meaningful. In this report, we describe a new metric we refer to as ‘Percent of Correlated Signal’. The following sections describe and demonstrate the metric with a case study event and a more rigorous test using a random selection of globally distributed events. While the focus here is on global tomography models, the metric is meant to applicable to regional ‘wiggle-for-wiggle’ waveform models/studies as well.

58 GEOSCIENCES↗
DOE OSTI2026

Data from: "Towards CONUS-Wide ML-Augmented Conceptually-Interpretable Modeling of Catchment-Scale Precipitation-Storage-Runoff Dynamics"

This data package was generated to support the manuscript “Towards CONUS-Wide Machine Learning-Augmented Conceptually Interpretable Modeling of Catchment-Scale Precipitation-Storage-Runoff Dynamics.” It provides input files, model outputs, plotting data, scripts, notebooks, and documentation used to develop, evaluate, and reproduce Mass-Conserving Perceptron (MCP)-based hydrologic modeling experiments across 513 selected Catchment Attributes and Meteorology for Large-sample Studies in the United States (CAMELS-US) basins. The files are organized by modeling component and analysis purpose, including rainfall–runoff experiments, snow module experiments, coupled hydrologic-snow experiments, Long Short-Term Memory (LSTM) benchmark results, model skill metrics, initialization and epoch records, cell-state normalization files, Akaike Information Criterion (AIC)-based model comparison files, and data used to generate manuscript figures. Tabular files can be opened using standard spreadsheet software or Python/R data-analysis tools. Python scripts, Jupyter notebooks, and selected MATLAB scripts are included for model execution, postprocessing, plotting, and statistical analysis. Quality assurance and quality control were conducted through the source-data selection and modeling workflow. Meteorological forcing, streamflow, and static catchment attributes were derived from the CAMELS-US dataset, and snow water equivalent data were derived from the University of Arizona (UA) Snow Water Equivalent dataset. Selected basins and time periods were screened during the associated research workflow to avoid missing observations or poor-quality cases. Static geospatial features were processed primarily using Quantum Geographic Information System (QGIS) and Geospatial Data Abstraction Library (GDAL) workflows. Additional details are provided in the associated manuscript and documentation.

ESS-DIVE CSV File Formatting Guidelines Reporting ↗
DOE OSTI2024

gcamfaostat: An R package to prepare, process, and synthesize FAOSTAT data for global agroeconomic and multisector dynamic modeling

The gcamfaostat R package is designed for the preparation, processing, and synthesis of the Food and Agriculture Organization (FAO) Statistics (FAOSTAT) agroeconomic data. The primary purpose is to facilitate FAOSTAT data use in global economic and multisector dynamic models while ensuring transparency, traceability, and reproducibility. Here, we provide an overview of the development of gcamfaostat (v1.0) and demonstrate its capabilities in generating and maintaining agroeconomic data required for the Global Change Analysis Model (GCAM). Our initiative seeks to enhance the quality and accessibility of data for the global agroeconomic modeling community, with the aim of fostering more robust and harmonized outcomes in a collaborative, efficient, and open-source framework. One of the important features of the package is the possibility to construct the FAO Food Balance Sheets at the disaggregated commodity level (with over 500 commodities), which provides a comprehensive and detailed data input for a variety of analytical and modeling applications. The processed data and visualizations offered by gcamfaostat can also be valuable to a broader audience interested in gaining insights into the intricacies of global agriculture.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Field Demonstration of The Krauklis Seismic Wave In A Novel Mva Method For Geologic CO 2 Storage

The Energy & Environmental Research Center (EERC) conducted a feasibility study of a new, low-impact geophysics monitoring method designed to incrementally track the saturation front of CO 2 injected into a reservoir. The method leverages a new way of transmitting energy from the surface to the reservoir and employs a new subsurface signal called the Krauklis wave (K-wave) and other guided waves for seismic monitoring. The EERC partnered with Seismos, Inc., the company that holds the exclusive technology license for this method, to advance development of K-wave technology, and Denbury Onshore, LLC, to field-test the system. The system was field-tested at Denbury’s Bell Creek oil field, currently undergoing CO 2 -based enhanced oil recovery (EOR) in southeastern Montana. The K-wave system uses a wellhead-mounted source to propagate energy down the wellbore in the form of a tube wave and through the well perforations into the reservoir. Part of the energy transmitted through the perforations into the reservoir converts to slow moving K-waves and/or other guided waves that are expected to propagate through the formation to neighboring wells. The propagated signal would induce tube waves in the neighboring wells that return to the surface to be recorded with a wellhead-mounted receiver. The propagation time of the signal is directly affected by the presence of CO 2 between the wells, so the difference in transit time of the returning signal between pre- and postinjection data sets could be used to help define the lateral extent of a CO 2 saturation front. Objectives of the study were to demonstrate, validate, and evaluate the K-wave method’s ability to monitor the morphology of injected CO 2 and estimate saturation distribution. Results from the demonstration were to be compared to conventional 4-D surface seismic methods, dynamic reservoir simulations, and CO 2 breakthrough data to validate the technology and provide useful information that may result in improvements. Prior to mobilizing for the full-scale K-wave baseline survey, a field test was conducted to better understand wellhead hardware configurations and test the response from the Bell Creek reservoir using three wells. Useful operational information resulted but signal from the first-generation gas release/gas injection source was not recorded in neighboring wells. After the initial field test, numerical modeling of the K-wave method was initiated using parameters specific to the Bell Creek Field. Modeling results indicated that the fundamental concepts behind the K-wave method exhibited problems and revealed a significantly greater attenuation of the signal at the wellbore perforations than expected, totaling on the order of 10-8 from source well to receiver well over a 400-m well spacing. Modeling also showed that greater source frequency bandwidth was needed to produce guided waves within a 30-foot-thick reservoir layer such as that at Bell Creek and that reservoir layer seismic velocity would need to be slower than the velocities of the bounding layers in order to act as a wave guide opposite to what exists at Bell Creek. Modeling results prompted a second field test with a stronger, second-generation displacement source developed in the interim. Concurrent with the second field test, a 3-D seismic survey was acquired to serve as a baseline for project validation and to serve as a monitor data set for the portion of the project area where CO 2 injection had already commenced. The new source was tested using a four-well geometry and failed to produce a detectable signal, which appeared to confirm modeling results and indicated that the Bell Creek Field is not an ideal site for the K-wave test. It was concluded that a source with sufficient strength and bandwidth could not be built to overcome the factors discovered by the modeling within a reasonable amount of time. Several other sites undergoing CO 2 -based EOR were also evaluated; however, none of these sites was amenable because of geology, limitations of source and receiver hardware, or the advanced state of the CO 2 flood. In addition, two other configurations of the source and receiver layout were considered; however, these design modifications either changed the scope of the project or did not meet the goals of the Carbon Storage Program. In consultation with the U.S. Department of Energy (DOE) Project Manager and Technology Manager, a no-go decision was invoked to forego additional expenditures as the technology would not meet current project objectives. Project activities to date have moved the technology forward by developing modeling equations that can be applied elsewhere, advanced the state of wellhead sensors and sources, and spurred improvements to field operations. The 3-D seismic survey in the project area and the dynamic reservoir simulations done as part of this study have provided value for ongoing and future synergistic activities at Bell Creek related to DOE carbon storage goals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Hydrodynamic Test Requirements Process Improvements

Hydrodynamic testing at Los Alamos National Laboratory would benefit from a process improvement for the requirements process. Cameo was used as a digital solution for requirements management to allow Lead Engineers to track requirements more effectively. This was identified as a process improvement throughout this Capstone project. This report includes a project proposal, business case, literature review, methodology, project plan, data analysis, decision-making report, financial analysis report, discussion, and conclusion. Initially, this project focused on figuring out a solution for the hydrotest requirement process improvements. The scope narrowed to focus on the use of Cameo for requirement capture and management. During this Capstone, four tests had digital models produced for requirements management in Cameo. The initial model was the baseline, with core requirements used across the tests. Commonalities in tests were used and the core requirements allowed for process efficiencies. In the data analysis, it was seen that overall, the implementation of using Cameo for requirements resulted in a decreasing trend for both schedule and normalized cost. Tests have different complexity levels which is also a factor in how long the requirements process will take. Additional data is needed to continue analyzing process improvements. Through decisionmaking and financial analysis, the recommendation was to use Cameo for requirements process improvement. Multiple experts provided feedback for requirements that were then captured within models. Numerous tangible and intangible benefits were identified with this process improvement. For return on investment, the metric of success was schedule reduction, which was overall seen. Four tests were analyzed, so future analysis will be needed. There is also not a great financial risk because the main cost would be purchasing more licenses. Individuals must generate requirements whether using this software or not. Overall, training is needed to help improve the skillsets of Lead Engineers but is already being supported on a regular frequency. A desktop guide was started associated with explaining the process, however, it is a work in progress. The team plans on adding additional information in the digital models to help status when requirements are met using verification methods and artifacts. From working on this project an improved understanding of Cameo and requirements was the result. There are future opportunities to extend the usage for requirements management and progress will continue after this project.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2019

Evolving R Coronae Borealis stars with mesa

The R Coronae Borealis (RCB) stars are rare hydrogen-deficient, carbon-rich supergiants. They undergo extreme, irregular declines in brightness of many magnitudes due to the formation of thick clouds of carbon dust. It is thought that RCB stars result from the mergers of CO/He white dwarf (WD) binaries. We constructed post-merger spherically symmetric models computed with the mesa code, and then followed the evolution into the region of the Hertzsprung-Russell (H-R) diagram where the RCB stars are located. We also investigated nucleosynthesis in the dynamically accreting material of CO/He WD mergers which may provide a suitable environment for significant production of 18 O and the very low 16 O/ 18 O values observed. We have also discovered that the N abundance depends sensitively on the peak temperature in the He-burning shell. Our mesa modelling consists of engineering the star by adding He-WD material to an initial CO-WD model, and then following the post-merger evolution using a nuclear-reaction network to match the observed RCB abundances as it expands and cools to become an RCB star. These new models are more physical because they include rotation, mixing, mass-loss, and nucleosynthesis within mesa . We follow the later evolution beyond the RCB phase to determine the stars’ likely lifetimes. The relative numbers of known RCB and extreme helium stars correspond well to the lifetimes predicted from the mesa models. In addition, most of computed abundances agree very well with the observed range of abundances for the RCB class.

Astronomy & Astrophysics↗
DOE OSTI2022

Intermountain West Energy Sustainability & Transitions Initiative: CO 2 Transport and Geologic Storage Modeling Results

This report and its associated Excel spreadsheet file provide the results of technoeconomic modeling of carbon dioxide (CO 2 ) pipeline transport modeling, CO 2 storage in deep saline subsurface formations, and CO 2 enhanced oil recovery (EOR). This report provides the methodology, describes the models and summarizes results for the three aspects of the analysis done for the Intermountain West Energy Sustainability and Transition (I-WEST) Initiative.

42 ENGINEERING↗
DOE OSTI2026

Triso Analysis Tool For Predictive Source Terms

Source term modeling for TRi-structural ISOtropic (TRISO) fuel has been performed for previous reactor designs, but few are available in the open literature. Thus, there is a need to develop a simple, versatile, and mechanistic model of fission product release and transport in gas reactor cores that can be applied to a variety of reactors through user inputs and reactor-specific radionuclide inventories. To meet this need, the TRISO Analysis Tool for Predictive Source terms (TRISO-ATOPS) was developed. This model calculates the release of the key safety-important fission products by diffusion through the kernel, silicon carbide and graphite based on fuel and graphite temperatures in the reactor under normal operation. These releases from the fuel enter the coolant where they can plate-out on cooler surfaces. A clean-up model is included for designs with a coolant purification system for removing fission gases. This initial distribution of fission products in the reactor serves as an initial condition for potential releases under postulated accident conditions. From this initial condition, the model calculates the fission product release for any transient temperature profile, and the fission product releases can then be used to assess radiological dose to the workers and the public using conventional radiological dose tools. Data on the diffusion of fission products is based on historic German TRISO experiments and the more current Department of Energy Advanced Gas Reactor TRISO fuel development program. The example cases in this work demonstrate the flexibility of the model

Stoyer, Benjamin [Idaho National Laboratory (INL),↗
DOE OSTI2024

Lower length scale model for palladium attack of silicon carbide in TRISO fuel

TRistructural ISOtropic (TRISO) particle fuels rely on silicon carbide (SiC) as the primary barrier for metallic fission product (FP) release. Palladium (Pd) generated by fission degrades the SiC layer, resulting in the formation of lamellar layers of palladium silicides (PdxSi) and carbon (C) perpendicular to the direction of attack. The Pd attack has been hypothesized to be responsible for failure of the SiC layer and enhance FP release. To better understand and quantify Pd attack of SiC in TRISO particles, a multiscale, mechanistic model of Pd transport is being developed by the NEAMS program. Previous work provided an initial hypothesis for modeling lamellar microstructure formation in SiC due to Pd attack using a phase-field model. The work described in this report builds on the previous model by using molecular dynamics (MD) simulations to parameterize the phase-field model kinetics, and build a reduced order model in BISON using the improved mesoscale Pd penetration model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Stochastic Framework for Optimal Control of Planetary Reentry Trajectories Under Multilevel Uncertainties

We present a novel stochastic optimal control framework that accounts for various types of uncertainties, with application to reentry trajectory planning. The formulation of the optimal trajectory control problem is presented in the context of an indirect method where a functional objective associated with the terminal vehicle speed is to be minimized. Uncertain input parameters in the optimal trajectory control model, including aerodynamic parameters and initial and terminal conditions, are modeled as aleatory random variables, while the statistical parameters of these aleatory distributions are themselves random variables. The parametric and model uncertainties are simultaneously propagated through an extended polynomial chaos expansion (EPCE) formalism. Several metrics are described to evaluate response statistics and presented as insightful tools for robust decision making. Specifically, the response probability density function (PDF) reflecting influence of both epistemic and aleatory uncertainties is obtained. By sampling over the random variables representing model error, an ensemble of response PDFs is generated and the associated failure probability is estimated as a random variable with its own polynomial chaos expansion. Besides, the sensitivity index functions of response PDF with respect to the statistical parameters are evaluated. Coupling parametric and model uncertainties within the EPCE framework leads to a robust and efficient paradigm for multilevel uncertainty propagation and PDF characterization in general optimal control problems.

Engineering↗
DOE OSTI2020

Lidar measurements of yawed-wind-turbine wakes: characterization and validation of analytical models

Abstract. Wake measurements of a scanning Doppler lidar mounted on the nacelle of a full-scale wind turbine during a wake-steering experiment were used for the characterization of the wake flow, the evaluation of the wake-steering set-up, and the validation of analytical wake models. Inflow-scanning Doppler lidars, a meteorological mast, and the supervisory control and data acquisition (SCADA) system of the wind turbine complemented the set-up. Results from the wake-scanning Doppler lidar showed an increase in the wake deflection with the yaw angle and that the wake deflection was not in all cases beneficial for the power output of a downstream turbine due to a bias of the inflow wind direction perceived by the yawed wind turbine and the wake-steering design implemented. Both observations could be reproduced with an analytical model that was initialized with the inflow measurements. Error propagation from the inflow measurements that were used as model input and the power coefficient of a waked wind turbine contributed significantly to the model uncertainty. Lastly, the span-wise cross section of the wake was strongly affected by wind veer, masking the effects of the yawed wind turbine on the wake cross sections.

17 WIND ENERGY↗
DOE OSTI2023

Development and modeling for a small-scale, rapidly heated high explosives initiation time (HEIT) experiment

Most small-scale assessments of explosive sensitivity, including the popular drop-weight impact test, convolute thermal and mechanical phenomena to the extent that it has been extremely challenging to decipher how an explosive ignites and propagates reactions. For instance, an impact generates heat through a range of dissipation mechanisms, which can in turn, depending on the reaction rates of the explosive, lead to chemical decomposition. To deconvolute the various contributions to the sub-shock initiation and propagation of explosive reactions, we describe the development and modeling of the High Explosives Initiation Time (HEIT) test - a new, small-scale, high-throughput experiment designed to rapidly heat milligram quantities of energetic materials confined within small diameter steel needles. Specifically, we have modeled and designed a 250 joule pulsed power system capable of rapidly delivering electrical current to the needles, resulting in rapid heat delivery to the sample. In conclusion, the energy deposition rate into the sample is controlled by different transmission line topologies. Modeling in COMSOL is performed to understand the energy required to heat up the explosive sample.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Combustion of 10–100 μm aluminum droplets in detonation products gases

Here we describe a two-phase model of combustion effects in aluminized high explosive (HE) charges. It is based on: (i) a Gas Dynamic Model of the expansion of the detonation product gases and their turbulent combustion with air; and (ii) a Heterogeneous Continuum Model of aluminum (Al) droplets and their combustion with the detonation product gases. Initial conditions are based on an analytical similarity solution for a cylindrical Chapman-Jouguet (CJ) detonation propagating at the CJ detonation velocity. The CJ jump conditions are computed at the thermodynamic equilibrium state by the Cheetah code, assuming the Al droplets are inert. We assume that the Al is 10 % of the charge mass and occurs as droplets at the CJ state. Different initial droplet diameters, ranging from 10 to 100 microns, are studied. A hydrodynamic combustion model based on large Damköhler numbers is employed in this study, -3σ(1 + 0.276 $\sqrt{Re}$)/($Kd^2_w$). It has a square-root dependence on the Reynolds number (Re) and inversed-squared dependence on the droplet diameter (d w ). The burnout time (t B ) of the Al droplets has a three-halves dependence on the droplet diameter, t B ~ ($d^0_w$) 3/2 . After burnout, the detonation products act as detonation products of the HE charge with active Al. They turbulently mix with air and form a combustion layer on the outer edge of the fireball. Details of the two-phase model, initial conditions and evolution of the flow field will be described.

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

SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD

We present a new capability of the ice sheet model SICOPOLIS that enables flexible adjoint code generation via source transformation using the open-source algorithmic differentiation (AD) tool OpenAD.The adjoint code enables efficient calculation of the sensitivities of a scalar-valued objective function or quantity of interest (QoI) to a range of important, often spatially varying and uncertain model input variables, including initial and boundary conditions, as well as model parameters.Compared to earlier work on the adjoint code generation of SICOPOLIS, our work makes several important advances:(i) it is embedded within the up-to-date trunk of the SICOPOLIS repository – accounting for 1.5 decades of code development and improvements – and is readily available to the wider community;(ii) the AD tool used, OpenAD, is an open-source tool;(iii) the adjoint code developed is applicable to both Greenland and Antarctica, including grounded ice as well as floating ice shelves, with an extended choice of thermodynamical representations.A number of code refactorization steps were required. They are discussed in detail in an Appendix as they hold lessons for the application of AD to legacy codes at large.As an example application, we examine the sensitivity of the total Antarctic Ice Sheet volume to changes in initial ice thickness, austral summer precipitation, and basal and surface temperatures across the ice sheet.Simulations of Antarctica with floating ice shelves show that over 100 years of simulation the sensitivity of total ice sheet volume to the initial ice thickness and precipitation is almost uniformly positive, while the sensitivities to surface and basal temperature are almost uniformly negative. Sensitivity to austral summer precipitation is largest on floating ice shelves from Queen Maud to Queen Mary Land. The largest sensitivity to initial ice thickness is at outlet glaciers around Antarctica. Comparison between total ice sheet volume sensitivities to surface and basal temperature shows that surface temperature sensitivities are higher broadly across the floating ice shelves, while basal temperature sensitivities are highest at the grounding lines of floating ice shelves and outlet glaciers. A uniformly perturbed region of East Antarctica reveals that, among the four control variables tested here, total ice sheet volume is the most sensitive to variations in austral summer precipitation as formulated in SICOPOLIS.Comparison between adjoint- and finite-difference-derived sensitivities shows good agreement, lending confidence that the AD tool is producing correct adjoint code.The new modeling infrastructure is freely available at http://www.sicopolis.net (last access: 2 April 2020) under the development trunk.

58 GEOSCIENCES↗
DOE OSTI2025

Model Intercomparison of the Impacts of Varying Cloud Droplet–Nucleating Aerosols on the Life Cycle and Microphysics of Isolated Deep Convection

The microphysical impacts of aerosol particles on scattered isolated deep convective cells near Houston, Texas, on 19 June 2013, are examined using multiple cloud-resolving model (CRM) simulations initialized with vertical profiles of low and high concentrations of cloud droplet–nucleating aerosols. These simulations formed part of the Model Intercomparison Project (MIP) conducted by the Deep Convective Working Group of the Aerosol, Cloud, Precipitation and Climate (ACPC) initiative. Each CRM generated a field of convective cells representing those observed during the case study with varying degrees of accuracy. The Tracking and Object-Based Analysis of Clouds (tobac) cell-tracking algorithm was applied to each MIP CRM simulation to track relatively long-lived convective cells (20–60 min). Most of the CRMs produced similar aerosol loading impacts on the warm phase of tracked cell properties with reduced autoconversion and accretion growth of rain, increased cloud water, reduced rainfall, and reduced near-surface evaporation of rain. The sign of aerosol impacts on the warm-phase properties of the convective cells was also quite consistent over cell lifetimes with the greatest magnitude of influence in the first half of the life cycle in most CRMs. In contrast, the ice-phase response to aerosol loading was highly variable among CRMs and included increases or decreases in ice amounts at inconsistent stages of the cell life cycle and midlevel versus upper-level changes in ice. This intermodel variability in ice is indicative both of the complex indirect interactions between aerosols and ice-phase processes in deep convection and their associated parameterizations.

Aerosol-cloud interaction↗
DOE OSTI2019

The paleoclimatic footprint in the soil carbon stock of the Tibetan permafrost region

Tibetan permafrost largely formed during the late Pleistocene glacial period and shrank in the Holocene Thermal Maximum period. Quantifying the impacts of paleoclimatic extremes on soil carbon stock can shed light on the vulnerability of permafrost carbon in the future. Here, we synthesize data from 1114 sites across the Tibetan permafrost region to report that paleoclimate is more important than modern climate in shaping current permafrost carbon distribution, and its importance increases with soil depth, mainly through forming the soil's physiochemical properties. We derive a new estimate of modern soil carbon stock to 3 m depth by including the paleoclimate effects, and find that the stock ( 36 .6 &#x2212;<!-- - --> 2 .4 + 2 .3 PgC) is triple that predicted by ecosystem models (11.5 ± 4.2 s.e.m PgC), which use pre-industrial climate to initialize the soil carbon pool. The discrepancy highlights the urgent need to incorporate paleoclimate information into model initialization for simulating permafrost soil carbon stocks.

54 ENVIRONMENTAL SCIENCES↗