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 379 records · Page 21

Development of a kinetic-thermodynamic model for lime-stabilization of Na-bentonite

This study presents the first kinetic model to predict the solid and pore solution composition of Na-bentonite clay reacting with slaked lime over a period of 720 days. The model successfully accounts for most experimental data using a single kinetic rate constant. The following sequence of reactions was predicted by the model: initial rapid dissolution of portlandite within the first 7 days, leading to a decrease in pH and dissolved calcium, and concurrent formation of calcium silicate hydrates (C-S-H: jennite), calcium aluminate hydrate (C-A-H: C₄AH₁₃), calcium aluminosilicate hydrates (stratlingite) and hydrotalcite. After 7 days, jennite and stratlingite are predicted to transform into tobermorite-II, contributing to strength development up to 28 days. From 28 to 90 days, continued montmorillonite dissolution is predicted, along with minor formation of ettringite, partial tobermorite-II dissolution, and precipitation of secondary phases such as albite and talc. Experimentally, portlandite dissolution was confirmed by TGA and XRD and found to be complete within 7 days, in agreement with model predictions. However, other predicted solid-phase transformations (e.g., tobermorite-II formation and dissolution, ettringite, albite, and talc formation) could not be conclusively verified through experimental techniques. Aqueous phase measurements confirmed that the pH and Ca trends in solution, and that equilibrium was reached by 90 days.

Chemical kinetics↗

Use of Shallow Ice Core Measurements to Evaluate and Constrain 1980–1990 Global Reanalyses of Ice Sheet Precipitation Rates

Abstract Sea‐level rise (SLR) projections by Earth System Models (ESMs) depend on ice sheet surface mass balances. Accurate, global atmosphere reanalyses would be ideal for providing equilibrated ice sheet model initial conditions in fully coupled ESM simulations. Here we present the first evaluation of 1980–1990 global reanalysis precipitation over Greenland and Antarctica that uses independent observations of net accumulation rates derived from shallow ice cores. Precipitation distributions from both the European Centre for Medium‐Range Weather Forecast's Reanalysis (ERA5) and the Modern‐Era Retrospective Analysis for Research and Applications (MERRA‐2) are highly correlated with contemporaneous co‐located net accumulation rates from Greenland ( r 2 > 0.95) and West Antarctica ( r 2 > 0.7). Three other commonly used reanalyses (WFDE5, CRUNCEP, and GSWP3) exhibit significantly weaker correlations on one or both ice sheets. Our findings imply that ESMs should use ERA5 or MERRA‐2 in data‐forced simulations to validate ice sheet model dynamics and precondition firn for SLR projections.

Schneider, Adam↗

Improving an Acoustic Vehicle Detector Using an Iterative Self-Supervision Procedure

In many non-canonical data science scenarios, obtaining, detecting, attributing, and annotating enough high-quality training data is the primary barrier to developing highly effective models. Moreover, in many problems that are not sufficiently defined or constrained, manually developing a training dataset can often overlook interesting phenomena that should be included. To this end, we have developed and demonstrated an iterative self-supervised learning procedure, whereby models are successfully trained and applied to new data to extract new training examples that are added to the corpus of training data. Successive generations of classifiers are then trained on this augmented corpus. Using low-frequency acoustic data collected by a network of infrasound sensors deployed around the High Flux Isotope Reactor and Radiochemical Engineering Development Center at Oak Ridge National Laboratory, we test the viability of our proposed approach to develop a powerful classifier with the goal of identifying vehicles from continuously streamed data and differentiating these from other sources of noise such as tools, people, airplanes, and wind. Using a small collection of exhaustively manually labeled data, we test several implementation details of the procedure and demonstrate its success regardless of the fidelity of the initial model used to seed the iterative procedure. Finally, we demonstrate the method’s ability to update a model to accommodate changes in the data-generating distribution encountered during long-term persistent data collection.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Experimental Simulations of Lunar Magma Ocean Crystallization: The Plot (But Not the Crust) Thickens

Numerical models of differentiation of a global-scale lunar magma ocean (LMO) have raised as many questions as they have answered. Recent orbital missions and sample studies have provided new context for a large range of lithologies, from the comparatively magnesian "purest anorthosite" reported by to Si-rich domes and spinel-rich clasts with widespread areal distributions. In addition, the GRAIL mission provided strong constraints on lunar crustal density and average thickness. Can this increasingly complex geology be accounted for via the formation and evolution of the LMO? We have in recent years been conducting extensive sets of petrologic experiments designed to fully simulate LMO crystallization, which had not been attempted previously. Here we review the key results from these experiments, which show that LMO differentiation is more complex than initial models suggested. Several important features expected from LMO crystallization models have yet to be reproduced experimentally; combined modelling and experimental work by our group is ongoing.

Draper, D. S.↗

Data assimilation and model evaluation experiment datasets

The Institute for Naval Oceanography, in cooperation with Naval Research Laboratories and universities, executed the Data Assimilation and Model Evaluation Experiment (DAMEE) for the Gulf Stream region during fiscal years 1991-1993. Enormous effort has gone into the preparation of several high-quality and consistent datasets for model initialization and verification. This paper describes the preparation process, the temporal and spatial scopes, the contents, the structure, etc., of these datasets. The goal of DAMEE and the need of data for the four phases of experiment are briefly stated. The preparation of DAMEE datasets consisted of a series of processes: (1) collection of observational data; (2) analysis and interpretation; (3) interpolation using the Optimum Thermal Interpolation System package; (4) quality control and re-analysis; and (5) data archiving and software documentation. The data products from these processes included a time series of 3D fields of temperature and salinity, 2D fields of surface dynamic height and mixed-layer depth, analysis of the Gulf Stream and rings system, and bathythermograph profiles. To date, these are the most detailed and high-quality data for mesoscale ocean modeling, data assimilation, and forecasting research. Feedback from ocean modeling groups who tested this data was incorporated into its refinement. Suggestions for DAMEE data usages include (1) ocean modeling and data assimilation studies, (2) diagnosis and theoretical studies, and (3) comparisons with locally detailed observations.

Lai, Chung-Cheng A.↗

n-Propylbenzene Droplet Combustion: Experiments and Numerical Modeling

n-Propylbenzene (PB, C9H12, 432K boiling point) is important as a constituent in surrogates for some transportation fuels but little is understood about its droplet burning characteristics. We present an experimental and numerical study of a fuel droplet burning under conditions where the sole means of gas transport is fuel evaporation at the droplet surface thereby creating a one-dimensional gas transport that is amenable to detailed numerical modeling. The one-dimensional gas transport that results is amenable to detailed numerical modeling. Initial droplet diameters (Do) ranged from 0.3 mm to about 5 mm to better understand the mechanisms that control burning over such a wide range of Do. For Do < 1 mm droplets were burned while anchored to 14 m SiC fibers within a sealed container in free-fall during which time the droplets were ignited by spark discharge and the burning process recorded by video cameras. Larger and free-floating droplets were studied in the spaced-based platform of the orbiting International Space Station for Do > 1.0 mm which provided unlimited times for experimental observation. Measurements of droplet, flame and soot shell diameters were made and are discussed in the presentation. The experimental results showed formation of optically thick soot clouds developing during burning that required a robust image analysis algorithm to enable measurement of droplet diameter. The measured droplet diameters were well predicted by a numerical model that included soot formation, a PB kinetic mechanism comprising 576 species and 27369 reactions, unsteady gas and liquid transport, variable properties and radiative transport. Flame diameters showed more variability compared to the simulations, depending on how the flame was defined (peak gas temperature and OH concentration). Droplets burned slower as Do was increased and the PB sooting propensity increased with Do. Flame extinction was observed for Do > 3.3 mm indicating a strong effect of radiation as Do was increased compared to the smaller droplet diameters examined where extinction was not experimentally observed. For the largest droplets examined the possibility that the observed extinction is indicative of a transition to cool flame burning is still under investigation.

Guo, Songtao↗

Updated Trends of the Stratospheric Ozone Vertical Distribution in the 60°S-60°N Latitude Range Based on the LOTUS Regression Model

This study presents an updated evaluation of stratospheric ozone profile trends in the 60°S - 60°N latitude range over the 2000 - 2020 period using an updated version of the Long-term Ozone Trends and Uncertainties in the Stratosphere (LOTUS) regression model that was used to evaluate such trends up to 2016 for the last WMO Ozone Assessment (2018). In addition to the derivation of detailed trends as a function of latitude and vertical coordinates, the regressions are performed with the data sets averaged over broad latitude bands, i.e., 60°S–35°S, 20°S–20°N and 35°N–60°N. The same methodology as in the last Assessment is applied to combine trends in these broad latitude bands in order to compare the results with the previous studies. Longitudinally resolved merged satellite records are also considered in order to provide a better comparison with trends retrieved from ground-based records, e.g., lidar, ozone sondes, Umkehr, microwave and Fourier Transform Infrared (FTIR) spectrometers at selected stations where long-term time series are available. The study includes a comparison with trends derived from the REF-C2 simulations of the Chemistry Climate Model Initiative (CCMI-1). This work confirms past results showing an ozone increase in the upper stratosphere, which is now significant in the three broad latitude bands. The increase is largest in the northern and southern hemisphere midlatitudes, with ~2.2%/decade at ~2.1 hPa, and ~2.1%/decade at ~3.2 hPa respectively, compared to ~1.6%/decade at ~2.6 hPa in the tropics. New trend signals have emerged from the records, such as a significant decrease of ozone in the tropics around 35 hPa and a non-significant increase of ozone in the southern midlatitudes at about 20 hPa. Non-significant negative ozone trends are derived in the lowermost stratosphere, with the most pronounced trends in the tropics. While a very good agreement is obtained between trends from merged satellite records and the CCMI-1 REF-C2 simulation in the upper stratosphere, observed negative trends in the lower stratosphere are not reproduced by models at southern and, in particular, at northern midlatitudes, where models report an ozone increase. However, the lower stratospheric trend uncertainties are quite large, for both measured and modelled trends. Finally, 2000-2020 stratospheric ozone trends derived from the ground-based and longitudinally resolved satellite records are in reasonable agreement over the European Alpine and tropical regions, while at the Lauder station in the southern hemisphere mid-latitudes they show some differences.

Stratospheric ozone trends↗

Computational Inference of Vibratory System with Incomplete Modal Information Using Parallel, Interactive and Adaptive Markov Chains

Inverse analysis of vibratory system is an important subject in fault identification, model updating, and robust design and control. It is challenging subject because 1) the problem is oftentimes underdetermined while the measurements are limited and/or incomplete; 2) many combinations of parameters may yield results that are similar with respect to actual response measurements; and 3) uncertainties inevitably exist. The aim of this research is to leverage upon computational intelligence through statistical inference to facilitate an enhanced, probabilistic framework using incomplete modal response measurement. This new framework is built upon efficient inverse identification through optimization, whereas Bayesian inference is employed to account for the effect of uncertainties. To overcome the computational cost barrier, we adopt Markov chain Monte Carlo (MCMC) to characterize the target function/distribution. Instead of using single Markov chain in conventional Bayesian approach, we develop a new sampling theory with multiple parallel, interactive and adaptive Markov chains and incorporate it into Bayesian inference. This can harness the collective power of these Markov chains to realize the concurrent search of multiple local optima. The number of required Markov chains and their respective initial model parameters are automatically determined via Monte Carlo simulation-based sample pre-screening followed by K-means clustering analysis. These enhancements can effectively address the aforementioned challenges in finite element inverse analysis. The validity of this framework is systematically demonstrated through case studies.

K Zhou↗

An ensemble data assimilation modeling system for operational outdoor microalgae growth forecasting

Microalgae have received increasing attention as a potential feedstock for biofuel or biobased products. Forecasting the microalgae growth is beneficial for managers in planning pond operations and harvesting decisions. This study proposed a biomass forecasting system comprised of the Huesemann Algae Biomass Growth Model (BGM), the Modular Aquatic Simulation System in Two Dimensions (MASS2), ensemble data assimilation (DA), and numerical weather prediction Global Ensemble Forecast System (GEFS) ensemble meteorological forecasts. The novelty of this study is to seek the use of ensemble DA to improve both BGM and MASS2 model initial conditions with the assimilation of biomass and water temperature measurements and consequently improve short-term biomass forecasting skills. This study introduces the theory behind the proposed integrated biomass forecasting system, with an application undertaken in pseudo-real-time in three outdoor ponds cultured with Chlorella sorokiniana in Delhi, California, United States. Results from all three case studies demonstrate that the biomass forecasting system improved the short-term (i.e., 7-day) biomass forecasting skills by about 60% on average, comparing to forecasts without using the ensemble DA method. Given the satisfactory performances achieved in this study, it is probable that the integrated BGM-MASS2-DA forecasting system can be used operationally to inform managers in making pond operation and harvesting planning decisions.

59 BASIC BIOLOGICAL SCIENCES↗

Initialization and assimilation of cloud and rainwater in a regional model

The initialization and assimilation of cloud and rainwater quantities in a mesoscale regional model was examined. Forecasts of explicit cloud and rainwater are made using conservation equations. The physical processes include condensation, evaporation, autoconversion, accretion, and the removal of rainwater by fallout. These physical processes, some of which are parameterized, represent source and sink in terms in the conservation equations. The question of how to initialize the explicit liquid water calculations in numerical models and how to retain information about precipitation processes during the 4-D assimilation cycle are important issues that are addressed.

Raymond, William H.↗

Implementation of Triply Periodic Minimal Surfaces (TPMS) as surface objects in OpenMC

Triply Periodic Minimal Surfaces (TPMS) represent a promising geometry for future fuel designs due to their significant surface-to-volume ratio, which facilitates efficient cooling of nuclear fuel, a crucial factor for safety and efficiency. Demonstrating the remarkable capabilities of TPMS fuel requires initial modeling and simulation. This paper presents an implementation of TPMS in the Monte Carlo code OpenMC, enabling reactor physics modeling of TPMS. Here, the primary advantages over traditional methods using CAD files include reduced memory requirements for computations and high-fidelity implementation. This implementation has been tested against CAD files loaded in Serpent2, yielding promising results with low biases in the $k_{\textrm{eff}}$, comparable to biases in the material balance sheet. The implementation presented in this work will be used in future reactor physics computations related to new reactor designs involving TPMS-based fuels.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Monitoring spatiotemporal evolution of fractures during hydraulic stimulations at the first EGS collab testbed using anisotropic elastic-waveform inversion

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. Here we use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

15 GEOTHERMAL ENERGY↗

The local structure of 0.5Ba(Zr 0.2 Ti 0.8 )O 3 -0.5(Ba 0.7 Ca 0.3 )TiO 3 from neutron total scattering measurements and multi-edge X-ray absorption analysis

Neutron total scattering measurements and multi-edge X-ray absorption analysis were performed on a 0.5Ba(Zr 0.2 Ti 0.8 )O 3 -0.5(Ba 0.7 Ca 0.3 )TiO 3 powder sample. The TiIII-II edges showed that titanium is distorted in the rhombohedral [111] direction (in a BaTiO 3 -like chemical and bonding environment). Extended X-ray fine structure analysis (EXAFS) was performed on the zirconium K-edge and it was found that zirconium resides in the center of the ZrO 6 octahedra in cubic symmetry. Furthermore, the local structure was initially modelled with the orthorhombic Amm2 space group with unsatisfactory results. To better model the local structure, the results from the EXAFS analysis were incorporated into the small-box pair distribution function (PDF) refinements and the best fit of the observed data was a combination of barium calcium titanate (BCT) and barium calcium zirconate (BCZ) phases.

36 MATERIALS SCIENCE↗

The Rise of Neural Networks for Materials and Chemical Dynamics

Machine learning (ML) is quickly becoming a premier tool for modeling chemical processes and materials. ML-based force fields, trained on large data sets of high-quality electron structure calculations, are particularly attractive due their unique combination of computational efficiency and physical accuracy. This Perspective summarizes some recent advances in the development of neural network-based interatomic potentials. Designing high-quality training data sets is crucial to overall model accuracy. One strategy is active learning, in which new data are automatically collected for atomic configurations that produce large ML uncertainties. Another strategy is to use the highest levels of quantum theory possible. Transfer learning allows training to a data set of mixed fidelity. A model initially trained to a large data set of density functional theory calculations can be significantly improved by retraining to a relatively small data set of expensive coupled cluster theory calculations. These advances are exemplified by applications to molecules and materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Strong Upper-Plate Heterogeneity at the Hikurangi Subduction Margin (North Island, New Zealand) Imaged by Adjoint Tomography

We use earthquake-based adjoint tomography to invert for three-dimensional structure of the North Island, New Zealand, and the adjacent Hikurangi subduction zone. The study area, having a shallow depth to the plate interface below the North Island, offers a rare opportunity for imaging material properties at an active subduction zone using land-based measurements. Starting from an initial model derived using ray tomography, we perform iterative model updates using spectral element and adjoint simulations to fit waveforms with periods ranging from 4–30 s. We perform 28 model updates using an L-BFGS optimization algorithm, improving data fit and introducing P- and S-wave velocity changes of up to ±30%. Resolution analysis using point spread functions show that our measurements are most sensitive to heterogeneities in the upper 30 km. The most striking velocity changes coincide with areas related to the active Hikurangi subduction zone. Lateral velocity structures in the upper 5 km correlate well with New Zealand geology. Here, the inversion reveals increased along-strike heterogeneity on the margin. In Cook Strait we observe a low-velocity zone interpreted as deep sedimentary basins. In the central North Island, low-velocity anomalies are linked to surface geology, and we relate velocity structures at depth to crustal magmatic activity below the Taupō Volcanic Zone.

58 GEOSCIENCES↗

Changes in Tropical Cyclones Undergoing Extratropical Transition in a Warming Climate: Quasi-Idealized Numerical Experiments of North Atlantic Landfalling Events

The current study extends earlier work that demonstrated future extratropical transition (ET) events will feature greater intensity and heavier precipitation to specifically consider potential changes in the impacts of landfalling ET events in a warming climate. A quasi-idealized modeling framework allows comparison of highly similar present-day and future event simulations; the model initial conditions are based on observational composites, increasing representativeness of the results. The future composite ET event features substantially more impactful weather conditions in coastal areas, with heavier precipitation and greater storm intensity. Specifically, a Category 2 present-day storm attained Category 4 Saffir-Simpson intensity in the future simulation and maintained greater intensity throughout the entire life cycle, although the storm undergoes less reintensification during the post-ET process, a result of reduced baroclinic conversion. These findings suggest increased potential for coastal hazards due to stronger tropical cyclone winds and heavier rainfall, leading to more severe coastal flooding and storm surge.

54 ENVIRONMENTAL SCIENCES↗

Corrugation-driven symmetry breaking in magic-angle twisted bilayer graphene

Abstract The discovery of unconventional superconductivity in magic-angle twisted bilayer graphene (tBLG) supported the twist-angle-induced flat band structure predictions made a decade earlier. Numerous physical properties have since been linked to the interlayer twist angle using the flat band prediction as a guideline. However, some key observations like the nematic phase and striped charge order behind the superconductivity are missing in this initial model. Here we show that a thermodynamically stable large out-of-plane displacement, or corrugation of the bilayer, induced by the interlayer twist, demonstrates partially filled states of the flat band structure, accompanied by a broken symmetry, in the magic-angle regime and the presence of symmetry breaking associated with the superconductivity in tBLG. The distinction between low and high corrugation can also explain the observed evolution of the vibrational spectra of tBLG as a function of twist angle. Our observation that large out-of-plane deformation modes enable partial filling of states near the Fermi energy may lead to a strategy for offsetting the effects of disorder in the local twist angle, which suppresses unconventional superconductivity and correlated insulator behavior in magic-angle tBLG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Tailoring diesel bioblendstock from integrated catalytic upgrading of carboxylic acids: a “fuel property first” approach

Abstract Lignocellulosic biomass offers the potential to produce renewable fuels at a scale commensurate with petroleum consumption. Hybrid approaches that combine biological and chemocatalytic processes have garnered increasing attention due to their flexibility for feedstock utilization and diversity of potential products. Of note, lignocellulosic sugars can be converted biologically to short-chain carboxylic acids, while subsequent chemocatalytic upgrading can elongate the carbon backbone and remove oxygen from the structure to produce drop-in hydrocarbon fuels. However, hybrid conversion processes are typically not designed with the fuel properties in mind a priori . In this work, we apply a “fuel property first” design approach to produce a tailored hydrocarbon bioblendstock with lower intrinsic sooting and drop-in diesel fuel potential. Initially, model predictions for six fuel properties critical to diesel applications (physicochemical requirements, energy content, safety considerations, autoignition ability, and sooting tendency) were used to screen an array of hydrocarbons accessible from upgrading individual and mixed C 2 /C 4 acids. This screening step allowed for down-selection to a non-cyclic branched C 14 hydrocarbon (5-ethyl-4-propylnonane) that can be synthesized from butyric acid through sequential catalytic reactions of acid ketonization, ketone condensation, and hydrodeoxygenation. Following evaluation of each conversion step with model compounds, butyric acid was then converted through an integrated catalytic process scheme to achieve >80% overall carbon yield to a hydrocarbon mixture product containing >60% of the target C 14 hydrocarbon. The potential of this conversion strategy to produce a hydrocarbon diesel bioblendstock from lignocellulosic biomass was then demonstrated using corn stover-derived butyric acid produced from Clostridium butyricum fermentation. Experimental fuel property testing of the purified C 14 blendstock validated the majority of the fuel property model predictions, including <50% of the intrinsic sooting tendency when compared to conventional diesel. Meanwhile, the crude conversion product met fuel property target metrics, validating conversion process development. When the C 14 bioblendstock was blended into a petroleum diesel at 20 vol.%, the blend maintained low cloud point, high energy density, and cetane number. Notably, the blend reduced sooting tendency by more than 10%, highlighting the potential of the tailored bioblendstock to reduce particulate emissions.

09 BIOMASS FUELS↗