Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model based development”

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 127 records · Page 7

Best practices in software development for robust and reproducible geoscientific models based on insights from the Global Carbon Budget's dynamic vegetation models

Computational models play an increasingly vital role in scientific research by enabling the numerical simulation of complex processes. Such models are also fundamental in geosciences. For instance, they offer critical insights into the impacts of global change on the Earth system today and in the future. Beyond their value as research tools, models are also software products and should therefore adhere to certain established software engineering standards. However, scientists are rarely trained as software developers, which can lead to potential deficiencies in software quality like unreadable, inefficient, or erroneous code. The complexity of models, coupled with their integration into broader workflows, also often makes it challenging to reproduce results, evaluate processes, and build upon them. In this paper, we review the state and current practices of the development processes of the state-of-the-art land surface models used by the Global Carbon Budget. We combine the experience of modelers from the respective research groups with the expertise of software engineers from tech companies to outline key principles and tools for improving software quality in research. We explore four main areas: (1) model testing and validation, (2) scientific, technical, and user documentation, (3) version control, continuous integration, and code review, and (4) the portability and reproducibility of workflows. Our review reveals that while modeling communities are incorporating many best practices, significant room for improvement remains in areas such as automated testing, automated documentation, and reproducibility. Therefore, we here identify and promote essential software engineering practices, including numerous examples of practices from within the community that can serve as guidelines for other models and could help streamline processes across the entire community. We conclude with an open-source example implementation of these principles, demonstrating portable and reproducible data flows, a continuous integration setup, and web-based visualizations. This example may serve as a practical resource for model developers, users, and all scientists engaged in scientific programming.

Gregor, Konstantin [Technical Univ. of Munich (Ger↗

Localization landscape theory of disorder in semiconductors. II. Urbach tails of disordered quantum well layers

Urbach tails in semiconductors are often associated to effects of compositional disorder. The Urbach tail observed in InGaN alloy quantum wells of solar cells and LEDs by biased photocurrent spectroscopy is shown to be characteristic of the ternary alloy disorder. The broadening of the absorption edge observed for quantum wells emitting from violet to green (indium content ranging from 0% to 28%) corresponds to a typical Urbach energy of 20 meV. A three-dimensional absorption model is developed based on a recent theory of disorder-induced localization which provides the effective potential seen by the localized carriers without having to resort to the solution of the Schrödinger equation in a disordered potential. This model incorporating compositional disorder accounts well for the experimental broadening of the Urbach tail of the absorption edge. For energies below the Urbach tail of the InGaN quantum wells, type-II well-to-barrier transitions are observed and modeled. This contribution to the below-band-gap absorption is particularly efficient in near-ultraviolet emitting quantum wells.When reverse biasing the device, the well-to-barrier below-band-gap absorption exhibits a red-shift, while the Urbach tail corresponding to the absorption within the quantum wells is blue-shifted, due to the partial compensation of the internal piezoelectric fields by the external bias. The good agreement between the measured Urbach tail and its modeling by the localization theory thus demonstrates the applicability of the latter to compositional disorder effects in nitride semiconductors.

14 SOLAR ENERGY↗

X-ray-induced piezoresponse during X-ray photon correlation spectroscopy of PbMg 1/3 Nb 2/3 O 3

X-ray photon correlation spectroscopy (XPCS) holds strong promise for observing atomic-scale dynamics in materials, both at equilibrium and during non-equilibrium transitions. Here an in situ XPCS study of the relaxor ferroelectric PbMg 1/3 Nb 2/3 O 3 (PMN) is reported. A weak applied AC electric field generates strong response in the speckle of the diffuse scattering from the polar nanodomains, which is captured using the two-time correlation function. Correlated motions of the Bragg peak are also observed, which indicate dynamic tilting of the illuminated volume. This tilting quantitatively accounts for the observed two-time speckle correlations. The magnitude of the tilting would not be expected solely from the modest applied field, since PMN is an electrostrictive material with no linear strain response to the field. A model is developed based on non-uniform static charging of the illuminated surface spot by the incident micrometre-scale X-ray beam and the electrostrictive material response to the combination of static and dynamic fields. The model qualitatively explains the direction and magnitude of the observed tilting, and predicts that X-ray-induced piezoresponse could be an important factor in correctly interpreting results from XPCS and nanodiffraction studies of other insulating materials under applied AC field or varying X-ray illumination.

36 MATERIALS SCIENCE↗

An analysis of parameters affecting ampacity in aircraft bipolar MVDC power cables via coupled electrical, thermal, and computational fluid dynamic modelling

Abstract The next generation of aircraft, including more electric aircraft and all‐electric aircraft (AEA), requires electric power systems with high power density and low system mass specifications. Increasing the voltage of the system to the range of a few kV, medium voltage (MV), is a reasonable approach to achieving high‐power‐density and low‐system‐mass EPSs for aircraft applications. Higher voltages, however, pose many challenges for aviation MV power cables such as arcs and arc tracking, partial discharges (PDs), and thermal management. In this regard, thermal management is more challenging since heat transfer by convection is greatly reduced at wide‐body aircraft's cruising altitudes due to the reduced air pressure. In this paper, a finite element method (FEM) model is developed in COMSOL Multiphysics for an aircraft bipolar MVDC (±5 kV) power cable. Using the model, the maximum permissible cable current at a low pressure of 18.8 kPa (at an altitude of 12.2 km from sea level, the usual cruising altitude for wide‐body aircraft) is calculated. Also, an analytical model is developed based on analytical and proven empirical correlations governing conductive, radiative, and convective heat transfers at the steady state to estimate the ampacity of the bipolar cable system at reduced pressure. It was shown that the proposed analytical model can be used for atmospheric pressure and systems with a larger number of poles, expanding its range of applications. The results of the FEM and analytical models correlate at wide ranges of parameters such as ambient temperature, duct size, distance between the positive and negative pole cables, and the overall diameter of the cables. The influence of horizontal and vertical arrangement of poles is included in the analytical model. The results of this study can be used to design bipolar MVDC power cable systems for the envisaged wide‐body AEA.

Azizi, Arian↗

ROP model for PDC bits in Geothermal drilling

Geothermal energy is a renewable source of energy, where heat extraction is preferentially balanced with the reservoir's natural heat recharge rate. The objective of this paper is to present and validate a novel rate of penetration (ROP) model for drilling hard and abrasive formations including granite formations for polycrystalline diamond compact (PDC) bits. The ROP model was developed based on a derived relationship of a threshold weight on cutter (WOC) and its corresponding depth of cut (DOC) for a single cutter. Laboratory data was used to scale the derived single cutter relationship to a full-hole ROP model for PDC bits. The ROP model includes a non-linear correlation for Phase I (inefficient drilling due to low WOB values) and a linear Phase II (efficient drilling) ROP response to WOB. The ROP model was verified using measured drilling parameter data from Utah FORGE well# 58-32 and data from Chocolate Mountains well # 17-8 in Southern California. When compared to oil and gas well drilling, geothermal drilling in granitic formations can be more difficult and complicated due to rock hardness and high temperatures. PDC bits can increase ROP and optimize drilling for these types of hard formations. This paper provides novel insight into the ROP response of PDC bits to drilling operational parameters.

15 GEOTHERMAL ENERGY↗

A Predictive Transient Model of the TREAT-SIRIUS Experiments

To support experiments modeling and transient analysis of the NASA-sponsored SIRIUS experiments in the Transient Reactor Test Facility (TREAT), an innovative predictive transient model has been developed to simulate various transient experiments. The developed model is based on a new approach that utilizes steady-state Monte Carlo calculations along with a surrogate model based on polynomial regression to determine the core reactivity. The reactivity is supplied to an exact point kinetics model that determines the reactor power, and it is coupled to an adiabatic feedback model to determine the fuel average temperature. In this work, the newly developed model is introduced along with validation test results considering SIRIUS-1 experiments at different power levels. The initial test results show a very good agreement with experimental data at different power levels, which indicates that the model can be used to predict the power of other transient tests.

42 ENGINEERING↗

cabbi-bio/BECT-ABM

The emergence of the cellulosic bioeconomy requires not only adequate technological, economic, and policy advances, but also effective communication and coordination among the multiple stakeholders in the bioeconomy community. Aiming to facilitate stakeholder communication and collaboration for cellulosic bioeconomy growth, an agent-based model is developed as a community communication tool in this study. The tool simulates the behaviors and interactions of key stakeholders, as well as the responses and interactions between the stakeholders and biophysical, socioeconomic, and environmental conditions. The agent behavior rules are developed and validated based on behavior and economic theories and information attained from a stakeholder focus group study and a land use survey.

Yang, Pan↗

AI-Batt (Autonomous Identification of Battery Life Models) [SWR 21-36]

Autonomous Identification of Battery Life Models (AI-Batt) AI-Batt is a MATLAB code base for developing lifetime models for batteries from accelerated aging data. The code base provides many functions for processing, visualizing, and modeling battery aging data, making the data processing, exploration, and modeling workflow substantially faster. These tools are tailored for working with battery aging data sets, which usually consist of many separate time-series for each cell, with many test conditions and possible replicates at each condition, which makes it difficult to simply process or visualize the data set. Complex modeling tasks, such as cross-validation, sensitivity analysis, and uncertainty quantification have been implemented to enable thorough statistical investigation of model predictions. Additionally, several machine-learning algorithms are implemented to autonomously identify suitable models via symbolic regression. Data processing functions automatically cast data from the struct data type, which is commonly used to store experimental data, but is not an acceptable input for most algorithms, to the table data type, which can be easily used as input to any optimization algorithm. Also, the data can be separated into time-invariant and time-variant data tables, which is helpful for exploring the data set as well as developing separate models for time-variant and time-invariant aging mechanisms. For example, in aging tests with constant temperature, temperature is a time-invariant experimental condition. Visualization tools enable plotting of data, model fits, and model simulations possible with single-line function calls, empowering data exploration of complex data sets with both time-varying and time-invariant trends. Plots can be automatically generated for the whole data set, or separated by data group (groups of test replicates) or individual data series. Data points or data series can be automatically colored by the value of a variable with a variety of color maps, and model predictions can also be colored by the value of a fit statistic. Comparisons between data sets and the predictions/simulations of different models on the same data set can be easily plotted as well. Distributions of parameter values from bootstrap resampling can be plotted to visualize the reliability of parameter estimation, or determine any correlations between parameters. Modeling tools handle the complex task of creating and parsing symbolic equations for modeling battery lifetime. Equations are parsed to grab relevant data variables, parameter values, or specified sub-models for input into optimization, evaluation, or simulation functions. Models can be optimized locally (one set of parameters for each data series), bi-level (some parameters shared across the data set), or globally (single set of parameters for all data). Functions implementing symbolic regression algorithms help users to discover effective model equations, even in poorly sampled, high-dimensional data.

Smith, Kandler↗

Empirical Dynamic Modeling Reveals Complexity of Methane Fluxes in a Temperate Salt Marsh

Abstract Methane dynamics within salt marshes are complex because vegetation types, temperature, oscillating water levels, and changes in salinity and redox conditions influence CH 4 production, consumption, oxidation, and emissions. These non‐linear and complex interactions among variables affect the traditionally expected functional relationships and present challenges for interpreting and developing process‐based models. We employed empirical dynamic modeling (EDM) and convergent cross mapping (CCM) as a novel approach for characterizing seasonal/multiday and diurnal CH 4 dynamics by inferring causal variables, lags, and interconnections among multiple biophysical variables within a temperate salt marsh using 5 years of eddy covariance data. EDM/CCM is a nonparametric approach capable of quantifying the coupling between variables while determining time scales where variable interactions are the most relevant. We found that gross primary productivity, tidal creek dissolved oxygen, and temperature were important for seasonal/multiday dynamics (rho = 0.73–0.80), while water level was most important for diurnal dynamics during both the growing and dormancy phenoperiods (rho = 0.72 and 0.56, respectively). Lags for the top‐ranked variables (i.e., gross primary productivity, dissolved oxygen, temperature, water level) occurred between 1 and 5 weeks at the seasonal scale and 1–24 hr at the diurnal scale. The EDM had high prediction capabilities for intra‐/inter‐seasonal patterns and annual CH 4 sums but had limitations in representing large, infrequent fluxes. Results highlight the importance of non‐linearity, drivers, lag times, and interconnections among multiple biophysical variables that regulate CH 4 fluxes in tidal wetlands. This research introduces a novel approach to examining CH 4 fluxes, which will aid in evaluating current paradigms in wetlands and other ecosystems.

Environmental Sciences & Ecology↗

Model-based optimizations of packaged rooftop air conditioners using low global warming potential refrigerants

Based on laboratory investigations for R-22 and R-410A alternative low GWP refrigerants in two baseline rooftop air conditioners (RTU), the DOE/ORNL Heat Pump Design Model was used to model the two RTUs and the models were calibrated against the experimental data. We compared the compressor efficiencies and heat exchanger performances. An efficiency-based compressor mapping method was developed. Extensive model-based optimizations were conducted to provide a fair comparison between all the low GWP candidates by selecting optimal configurations. The results illustrate that all the R-22 low GWP refrigerants will lead to slightly lower COPs. ARM-20B appears to be the best R-22 replacement at normal conditions. At higher ambient temperatures, ARM-20A exhibits better performance. All R-410A low GWP candidates will result in similar or better efficiencies than R-410A. R-32 has the best COP while requiring the smallest compressor. Finally, R-452B uses the closest compressor displacement volume and achieves the same efficiency as R-410A.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Dynamic Differential Image Circle Diameter Measurement Precision Assessment: Application to Burning Droplets

Dynamic measurement precision assessment has been achieved for a differential circle measurement application. Differential circle diameter measurement, in image analysis, typically requires fitting a circle model that optimizes for image distortions, defects or occlusions. The differential task occurs when precise measurements of diameter change are required given object size variation with time. An automated system was designed to provide diameter measurements and associated measurement precision of images of a fuel droplet undergoing combustion in zero gravity for the FLEX-2 dataset. An image gradient-based, least-squares boundary point fitting method to a circle or ellipse model is used for diameter measurement. The presence of soot aggregates poses significant challenges for diameter measurements when it occludes part of the droplet boundary. The precision of the diameter measurements depends upon the image quality. Using synthetic image simulations that model the soot behavior, we developed a model based on image quality measures that assesses the measurement precision for each individual diameter measurement. Thus, diameter measurements with precision assessments were made available for follow-up scientific analysis. As a result, the algorithm's success rate for measurable runs was 98%. In cases of limited occlusion, a measurement precision of ±0.2 pixels for the FLEX-2 dataset was achieved.

42 ENGINEERING↗

Rate of MgO Pickup in Alumina Inclusions in Aluminum-Killed Steel

This work aims to clarify the rate and mechanism of MgO pickup by alumina inclusions and the effect of oxide impurities in MgO crucibles on this transformation. Two MgO crucibles from different batches from the same supplier were used in laboratory experiments with Al-killed steel. A kinetic model was developed, based on mass-transfer control in the liquid steel. Rate constants were fitted using inclusion analysis. The rate of magnesium transfer from the two types of crucibles was found to differ by a factor of 20; faster magnesium transfer was associated with formation of a slag layer on the inner surface of the crucible wall (rather than a solid spinel product layer). The kinetic model was also used to simulate industrial scale ladle refining (1) to illustrate the effects of total oxygen concentration and (2) to evaluate the contribution of steel-refractory reaction (in addition to steel-slag reaction) on the rate of MgO pickup in alumina inclusion. The rate of MgO pickup was higher with a lower inclusion concentration. For ladle desulfurization, the extent of MgO pickup in inclusions is directly linked to the extent of desulfurization; both reactions are controlled by the oxygen potential at the steel-slag interface.

36 MATERIALS SCIENCE↗

Prediction of Hydroxymethylfurfural Yield in Glucose Conversion through Investigation of Lewis Acid and Organic Solvent Effects

Hydroxymethylfurfural (HMF) is one of the important renewable platform compounds that can be obtained from biomass feedstocks through glucose conversion catalyzed by Brønsted and Lewis acids. However, it is challenging to enhance the HMF yield due to side reactions. In this study, a systematic approach combining theory and experiment was performed to investigate the influence of Lewis acids and organic solvents on the HMF yield. For the Lewis acid effect, a relationship between chemical hardness and experimental HMF yields was found in the rate-limiting step of glucose-to-fructose isomerization for six metal chlorides; HMF production was promoted when the metal chloride and a substrate had a similar chemical hardness. To study the organic solvent effect, a multivariate model was developed based on the insights gained from the mechanistic study of fructose dehydration, to predict HMF yields in a given water-organic cosolvent system. It showed a reliable accuracy in evaluating HMF yields with a mean absolute error (MAE) of 3.0% with respect to experimental HMF yields for 13 solvents, and also predicted HMF yields with a MAE of 10.7% for four new solvents. Chemical interpretation of the model revealed that it is desirable to use a solvent capable of stabilizing the carbocation intermediates with low proton transfer activity and high hydrogen bond basicity, to maximize the HMF yield. This multivariate model informs experimentalists about rational selection of solvents with very low computational costs needed to calculate only six variables for each solvent. It can be expanded to other catalytic systems such as heterogeneous Brønsted–Lewis bifunctional catalysts and enables optimization of reaction conditions to obtain other useful platform molecules through biomass conversion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling of deposit formation in mesoporous substrates via atomic layer deposition: Insights from pore‐scale simulation

Atomic layer deposition (ALD) has been a promising technique in fabricating membranes and tuning their properties with a precision at the atomic level. Fabrication of zeolitic imidazolate framework (ZIF) membranes using the ligand-induced permselectivation (LIPS) method starts with the formation of an oxide in a mesoporous substrate by ALD and is followed by the transformation of this oxide to ZIF using imidazolate vapor treatment. The objective of the ALD step is to block the mesopores with a thin deposit, that is, one with small penetration depth and small thickness on the top surface of the substrate. Unlike typical ALD on nonporous substrates, where all available sites react per ALD cycle, thin deposit formation in a mesoporous substrate requires that only a small fraction of the available deposition sites (i.e., close to the substrate surface) is subjected to ALD. Consequently, reactant dosing and duration of pulses are important process variables which, together with diffusion and reaction kinetics determine the deposit structure. Quantitative understanding of the interplay of these variables and phenomena can enable the rational design of ALD within mesoporous substrates. Here, we extend our earlier modeling effort considering the coexistence of ALD both inside the pores and on the external surface of the substrate. Finite-volume based models were developed and validated to simulate the two distinct modes of deposition cycle by cycle. The total mass uptake of the substrate with ALD cycles can be predicted using the combined surface deposition and pore reaction–diffusion models as affirmed by in situ quartz crystal microbalance experimental data. The ALD reactor model combined with the deposition model can accurately capture the number of ALD cycles needed to block the pores of the substrate. Based on the model, we designed a modified ALD process and examined the performance of the corresponding LIPS membranes. Furthermore, the present modeling work provides a new understanding of the deposit formation via ALD within mesoporous substrates for a variety of membrane applications.

atomic layer deposition↗

Comparing Cycling Characteristics of Symmetric Lithium-Polymer-Lithium Cells with Theoretical Predictions

We develop a model based on concentrated solution theory for predicting the cycling characteristics of a lithium-polymer-lithium symmetric cell containing an electrolyte with known transport properties. The electrolytes used in this study are mixtures of polyethylene oxide (PEO) and lithium bis(trifluoromethanesulfonyl) imide (LiTFSI) salt, prepared over a wide range of salt concentrations. The transport properties of PEO/LiTFSI previously reported in the literature are used as inputs for our model. We calculate salt concentration and potential profiles, which develop in these electrolytes under a constant dc polarization, as a function of current density, electrolyte thickness, and salt concentration. These profiles are nonlinear at steady-state due to the strong concentration dependence of the transport properties of this electrolyte. The effect of this nonlinearity on limiting current is demonstrated. Cycling characteristics of a series of lithium symmetric cells were measured to test the validity of our model, without resorting to any adjustable parameters. The time-dependence and steady-state value of the potential measured during cycling experiments were in excellent agreement with model predictions.

25 ENERGY STORAGE↗

Effect of disturbances on MED-TVC plant characteristics: Dynamic modeling and simulation

This study describes a detailed steady-state analysis and transient operation of a multi-effect desalination plant with a thermal vapor compressor unit (MED-TVC). The model developed is based on solving the basic conservation equations of mass, energy, and salt. It can predict steady-state operation and transient behavior when the plant is subject to abrupt changes in the main input parameters including motive steam flow rate (m ms ), cooling seawater flow rate (m cw ) and cooling seawater temperature (T cw ). These sudden changes simulate the conditions a MED-TVC plant may be exposed to. Both the steady-state and dynamic aspects of the model are validated against published experimental data. Monitoring the variation of four state variables (vapor temperature, brine level, brine and vapor flow rates) in each effect shows that disturbances in the motive steam flow rate and the cooling seawater temperature have a significant effect on the plant performance. On the other hand, the disturbance in the cooling seawater mass flow rate has only a moderate effect. For all applied disturbances, the change in the brine level is the slowest compared to the changes in vapor temperature, and brine and vapor flow rates. Furthermore, the effect of changing seawater salinity shows a slight effect on the total distillate production and the specific heat transfer area.

42 ENGINEERING↗

Optimal CO 2 storage management considering safety constraints in multi-stakeholder multi-site GCS projects: A Markov game perspective

Geological carbon storage (GCS) projects could involve a diverse array of stakeholders or players from public, private, and regulatory sectors, each with different objectives and responsibilities. Given the complexity, scale, and long-term nature of GCS operations, determining whether individual stakeholders can independently optimize their interests — or whether collaborative coalition agreements are needed — remains a central question for effective GCS project planning and management. To access large, high-quality storage resources, future GCS deployment may increasingly occur in geologically connected sites, where shared geological features such as pressure space and reservoir pore capacity can lead to competitive behavior among stakeholders. In this work, we propose a paradigm based on Markov games to quantitatively investigate how different coalition structures affect the goals of stakeholders. We frame this multi-stakeholder multi-site problem as a multi-agent reinforcement learning problem with safety constraints. Our approach enables agents to learn optimal strategies while complying with safety regulations. We present an example where multiple operators are injecting CO 2 into their respective project areas in a geologically connected basin. To address the high computational cost of repeated simulations of high fidelity models, a previously developed surrogate model based on the Embed-to-Control (E2C) framework is employed. Our results demonstrate the effectiveness of the proposed framework in addressing optimal management of CO 2 storage when multiple stakeholders with different objectives and goals are involved.

58 GEOSCIENCES↗