Search NASA⌕ Search

SEARCH · Search NASA

Results for “mechanistic modeling”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Computational Fluid Dynamics (CFD) Simulations of Taylor Bubbles in Vertical and Inclined Pipes with Upward and Downward Liquid Flow

Summary Two-phase flow is a common occurrence in pipes of oil and gas developments. Current predictive tools are based on the mechanistic two-fluid model, which requires the use of closure relations to predict integral flow parameters such as liquid holdup (or void fraction) and pressure gradient. However, these closure relations carry the highest uncertainties in the model. In particular, significant discrepancies have been found between experimental data and closure relations for the Taylor bubble velocity in slug flow, which has been determined to strongly affect the mechanistic model predictions (Lizarraga-García 2016). In this work, we study the behavior of Taylor bubbles in vertical and inclined pipes with upward and downward flow using a validated 3D computational fluid dynamics (CFD) approach with level set method implemented in a commercial code. A total of 56 cases are simulated, covering a wide range of fluid properties, pipe diameters, and inclination angles: Eo ∈ [10, 700]; Mo ∈ [1×10–6, 5×103]; ReSL ∈ [–40, 10]; θ ∈ [5°, 90°]. For bubbles in vertical upward flows, the simulated distribution parameter, C0, is successfully compared with an existing model. However, the C0 values of downward and inclined slug flows where the bubble becomes asymmetric are shown to be significantly different from their respective vertical upward flow values, and no current model exists for the fluids simulated here. The main contributions of this work are (1) the relatively large 3D numerical database generated for this type of flow, (2) the study of the asymmetric nature of inclined and some vertical downward slug flows, and (3) the analysis of its impact on the distribution parameter, C0.

Engineering↗
NASA NTRS2004

Using explanatory crop models to develop simple tools for Advanced Life Support system studies

System-level analyses for Advanced Life Support require mathematical models for various processes, such as for biomass production and waste management, which would ideally be integrated into overall system models. Explanatory models (also referred to as mechanistic or process models) would provide the basis for a more robust system model, as these would be based on an understanding of specific processes. However, implementing such models at the system level may not always be practicable because of their complexity. For the area of biomass production, explanatory models were used to generate parameters and multivariable polynomial equations for basic models that are suitable for estimating the direction and magnitude of daily changes in canopy gas-exchange, harvest index, and production scheduling for both nominal and off-nominal growing conditions. c2004 COSPAR. Published by Elsevier Ltd. All rights reserved.

NASA Discipline Life Support Systems↗
NASA NTRS1992

The influence of grazing on surface climatological variables of tallgrass prairie

Mass and energy exchange between most grassland canopies and the atmosphere are mediated by grazing activities. Ambient temperatures can be increased or decreased by grazers. Data have been assembled from simulated grazing experiments on Konza Prairie Research Natural Area and observations on adjacent pastures grazed by cattle show significant changes in primary production, nutrient content, and bidirectional reflectance characteristics as a function of grazing intensity. The purpose of this research was to provide algorithms that would allow incorporation of grazing effects into models of energy budgets using remote sensing procedures. The approach involved: (1) linking empirical measurements of plant biomass and grazing intensities to remotely sensed canopy reflectance, and (2) using a higher resolution, mechanistic grazing model to derive plant ecophysiological parameters that influence reflectance and other surface climatological variables.

Seastedt, T. R.↗
NASA NTRS2006

Modeling the Stress Strain Behavior of Woven Ceramic Matrix Composites

Woven SiC fiber reinforced SiC matrix composites represent one of the most mature composite systems to date. Future components fabricated out of these woven ceramic matrix composites are expected to vary in shape, curvature, architecture, and thickness. The design of future components using woven ceramic matrix composites necessitates a modeling approach that can account for these variations which are physically controlled by local constituent contents and architecture. Research over the years supported primarily by NASA Glenn Research Center has led to the development of simple mechanistic-based models that can describe the entire stress-strain curve for composite systems fabricated with chemical vapor infiltrated matrices and melt-infiltrated matrices for a wide range of constituent content and architecture. Several examples will be presented that demonstrate the approach to modeling which incorporates a thorough understanding of the stress-dependent matrix cracking properties of the composite system.

Morscher, Gregory N.↗
DOE OSTI2022

Physics-informed Bayesian machine learning case study: Integral blade rotors

This paper provides a physics-informed Bayesian machine learning (PIBML) description and case study. The PIBML approach applies three physics-based models to establish the initial beliefs before testing to determine the probability of milling stability (or prior). These include: receptance coupling substructure analysis (RCSA) prediction for the tool tip frequency response functions; finite element software prediction of the mechanistic force model coefficients; and a spindle speed-dependent power law model for process damping. Testing was then performed to identify optimal stable machining conditions using an expected improvement in material removal rate criterion. The prior probability of stability was updated using the test results to determine the posterior probability of stability. The test results were compared to the parameter recommendations provided by the endmill manufacturer. A demonstration integral blade rotor was machined at the optimal stable machining conditions for 304 stainless steel and 6061-T6 aluminum. Finally, the disagreement between manufacturer recommendations and milling performance in both materials tested emphasizes the need for broad implementation of PIBML approaches to increase machining productivity and efficiency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Rays for Roots - Integrating Backscatter X-Ray Phenotyping, Modeling and Genetics to Increase Carbon Sequestration and Switchgrass Resource Use (Final Report)

To increase carbon (C) deposition in the soil and enhance crop resource use efficiency, characterizing root form and function is essential. Several root and soil traits have been linked to increased root-to-soil C transfer. Technology that could provide high-resolution characterization of many of these traits in field conditions would revolutionize our ability to study and understand how to increase C sequestration. In this effort, we developed an initial early prototype backscatter X-ray system for non-destructive imaging of root traits. We collected initial backscatter X-ray data in field and lab settings and carried out early analysis of these data. Along with this prototype, we also developed a suite of root phenotyping approaches including advanced minirhizotron image analysis, soil core imaging, and mesocosm imaging. Minirhizotron (MR) tubes are clear tubes inserted into the soil in the field and used to image roots and the surrounding soil. Our team has developed deep learning-based methods that can segment roots from soil that can learn from imprecise image-level labels. The ability to learn or fine-tune our deep learning algorithms from image-level labels allows easier and faster application of these approaches to new locations and new plant species. We have successfully implemented and applied our MR analysis approaches to thousands of switchgrass MR images collected across geographical regions. An advantage of MR imaging is the ability to collect root and soil images over time. Our soil core analysis included collecting hundreds of soil core samples from harvested switchgrass fields and imaging these cores with both X-ray CT and backscatter X-ray imaging. Initial segmentation approaches for the X-ray CT images of these cores have been developed and applied. An advantage of soil core analysis is that it preserves the three-dimensional structures of the roots and soil in the core collected. Our group also developed photogrammetry-based mesocosm root imaging and phenotyping approaches. In this approach, a plant was grown in a large mesocosm with a three-dimensional grid of thin supporting lines inserted throughout the mesocosm. After the plant (and, correspondingly, the root architecture is grown and established) the soil media was removed and the supporting lines approximately preserved the three-dimensional root architecture. Then, we applied photogrammetry techniques to create a three-dimensional digital representation of the root architecture for which we developed analysis algorithms including skeletonization. We carried out our phenotyping development with powerful switchgrass resources and physiological and agroecosystem modeling to deliver novel technology. This project contributes to multiple ARPA-E missions including reduction of foreign imports of energy, reduction of energy-related emissions including greenhouse gases, and ensuring that the United States maintains a technological lead in developing and deploying advanced energy technology. Furthermore, the developed tools could transform public and private plant breeding and could be broadly applicable to other crops and, potentially, other application areas. Our team of engineers, plant and soil scientists, and modelers i) developed an early prototype backscatter X-ray platform that can operate in field conditions; ii) developed a suite of root phenotyping and characterization approaches as described above; iii) developed and carried out plant biology and physiology roots studies and; iv) developed and implemented mechanistic physiological modeling.

42 ENGINEERING↗
DOE OSTI2022

Resolving the Heat of Trimethylaluminum and Water Atomic Layer Deposition Half-Reactions

We report that atomic layer deposition (ALD) is a surface synthesis technique that is characterized by self-limiting reactions between gas-phase precursors and a solid substrate. Although ALD processes have been demonstrated that span the periodic table, a greater understanding of the surface chemistry that affords ALD is necessary to enable greater precision, including area- and site-selective growth. We offer new insight into the thermodynamics and kinetics of the trimethylaluminum (TMA) and H 2 O ALD half-reactions with calibrated and time-resolved in situ pyroelectric calorimetry. The half-reactions produce 3.46 and 2.76 eV/Al heat, respectively, which is greater than the heat predicted by computational models based on crystalline Al 2 O 3 substrates and closely aligned with the heat predicted by standard heats of formation. The pyroelectric thin-film calorimeter offers submilisecond temporal resolution that uniquely and clearly resolves precursor delivery and reaction kinetics. Both half-reactions are observed to exhibit multiple kinetic rates, with average TMA half-reaction rates at least 2 orders of magnitude faster than the H 2 O half-reaction kinetics. Comparing the experimental heat with published computational literature and additional first-principles modeling highlights the need to refine our models and mechanistic understanding of even the most ubiquitous ALD reactions.

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

Mathematical model of a personalized neoantigen cancer vaccine and the human immune system

Cancer vaccines are an important component of the cancer immunotherapy toolkit enhancing immune response to malignant cells by activating CD4 + and CD8 + T cells. Multiple successful clinical applications of cancer vaccines have shown good safety and efficacy. Despite the notable progress, significant challenges remain in obtaining consistent immune responses across heterogeneous patient populations, as well as various cancers. We present a mechanistic mathematical model describing key interactions of a personalized neoantigen cancer vaccine with an individual patient’s immune system. Specifically, the model considers the vaccine concentration of tumor-specific antigen peptides and adjuvant, the patient’s major histocompatibility complexes I and II copy numbers, tumor size, T cells, and antigen presenting cells. We parametrized the model using patient-specific data from a clinical study in which individualized cancer vaccines were used to treat six melanoma patients. Model simulations predicted both immune responses, represented by T cell counts, to the vaccine as well as clinical outcome (determined as change of tumor size). This model, although complex, can be used to describe, simulate, and predict the behavior of the human immune system to a personalized cancer vaccine.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Rapidly changing high-latitude seasonality: implications for the 21st century carbon cycle in Alaska

Seasonal variations in high-latitude terrestrial carbon (C) fluxes are predominantly driven by air temperature and radiation. At present, high-latitude net C uptake is largest during the summer. Recent observations and modeling studies have demonstrated that ongoing and projected climate change will increase plant productivity, microbial respiration, and growing season lengths at high-latitudes, but impacts on high-latitude C cycle seasonality (and potential feedbacks to the climate system) remain uncertain. Here we use ecosys, a well-tested and process-rich mechanistic ecosystem model that we evaluate further in this study, to explore how climate warming under an RCP8.5 scenario will shift C cycle seasonality in Alaska throughout the 21st century. The model successfully reproduced recently reported large high-latitude C losses during the fall and winter and yet still predicts a high-latitude C sink, pointing to a resolution of the current conflict between process-model and observation-based estimates of high-latitude C balance. We find that warming will result in surprisingly large changes in net ecosystem exchange (NEE; defined as negative for uptake) seasonality, with spring net C uptake overtaking summer net C uptake by year 2100. This shift is driven by a factor of 3 relaxation of spring temperature limitation to plant productivity that results in earlier C uptake and a corresponding increase in magnitude of spring NEE from –19 to –144 gC m –2 season –1 by the end of the century. Although a similar relaxation of temperature limitation will occur in the fall, radiation limitation during those months will limit increases in C fixation. Additionally, warmer soil temperatures and increased carbon inputs from plants lead to combined fall and winter C losses (163 gC m –2 ) that are larger than summer net uptake (123 gC m –2 season –1 ) by year 2100. However, this increase in microbial activity leads to more rapid N cycling and increased plant N uptake during the fall and winter months that supports large increases in spring NPP. Due to the large increases in spring net C uptake, the high-latitude atmospheric C sink is projected to sustain throughout this century. Our analysis disentangles the effects of key environmental drivers of high-latitude seasonal C balances as climate changes over the 21st century.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Changes in precipitation and air temperature contribute comparably to permafrost degradation in a warmer climate

Abstract Surface energy budgets of high-latitude permafrost systems are poorly represented in Earth system models (ESMs), yet permafrost is rapidly degrading and these dynamics are critical to future carbon-climate feedback predictions. A potentially important factor in permafrost degradation neglected so far by ESMs is heat transfer from precipitation, although increases in soil temperature and thaw depth have been observed following increases in precipitation. Using observations and a mechanistic ecosystem model, we show here that increases in precipitation hasten active layer development beyond that caused by surface air warming across the North Slope of Alaska (NSA) under recent and 21st century climate (RCP8.5). Modeled active layer depth (ALD) in simulations that allow precipitation heat transfer agreed very well with observations from 28 Circumpolar Active Layer Monitoring sites (R 2 = 0.63; RMSE = 10 cm). Simulations that ignored precipitation heat transfer resulted in lower spatially-averaged soil temperatures and a 39 cm shallower ALD by 2100 across the NSA. The results from our sensitivity analysis show that projected increases in 21st century precipitation deepen the active layer by enhancing precipitation heat transfer and ground thermal conductivity, suggesting that precipitation is as important an environmental control on permafrost degradation as surface air temperature. We conclude that ESMs that do not account for precipitation heat transfer likely underestimate ALD rates of change, and thus likely predict biased ecosystem responses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Tritium Transport Modeling at the Pore Scale in Ceramic Breeder Materials Using TMAP8

Fusion reactors depend on the blanket material to breed and release tritium at the same rate or faster than it is consumed by the fusion reaction. Cellular ceramic breeders (CCBs) are dense materials that can maintain a high tritium breeding ratio while promoting tritium release because of highly connected pores. Assessing the tritium breeding capabilities of these materials requires a combination of extensive experimental and modeling efforts. In this work, we develop and calibrate a multiphysics model of tritium transport. This novel model accounts for ceramic and pore diffusion, trapping and detrapping, and several surface reactions at the pore surface. We perform a sensitivity analysis and calibrate the model by comparing its predictions against experimental measurements of deuterium absorption. The calibrated model is then used to model tritium absorption in samples with different pore microstructures to investigate the effect of pore interconnectivity on tritium absorption. The model is part of the development of the multiscale, multiphysics framework for tritium transport [i.e., the Tritium Migration Analysis Program (TMAP8)], which is itself built on top of the finite-element multiphysics framework multiphysics object-oriented simulation environment (MOOSE). This study demonstrates some of TMAP8’s capabilities and is the first step toward assessing the tritium breeding capabilities of ceramic breeder material designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Cellular Response to Proton Irradiation: A Simulation Study with TOPAS-nBio

The cellular response to ionizing radiation continues to be of significant research interest in cancer radiotherapy, and DNA is recognized as the critical target for most of the biologic effects of radiation. Incident particles can cause initial DNA damages through physical and chemical interactions within a short time scale. Initial DNA damages can undergo repair via different pathways available at different stages of the cell cycle. The misrepair of DNA damage results in genomic rearrangement and causes mutations and chromosome aberrations, which are drivers of cell death. This work presents an integrated study of simulating cell response after proton irradiation with energies of 0.5–500 MeV (LET of 60–0.2 keV/µm). A model of a whole nucleus with fractal DNA geometry was implemented in TOPAS-nBio for initial DNA damage simulations. The default physics and chemistry models in TOPAS-nBio were used to describe interactions of primary particles, secondary particles, and radiolysis products within the nucleus. The initial DNA double-strand break (DSB) yield was found to increase from 6.5 DSB/Gy/Gbp at low-linear energy transfer (LET) of 0.2 keV/µm to 21.2 DSB/Gy/Gbp at high LET of 60 keV/µm. A mechanistic repair model was applied to predict the characteristics of DNA damage repair and dose response of chromosome aberrations. It was found that more than 95% of the DSBs are repaired within the first 24 h and the misrepaired DSB fraction increases rapidly with LET and reaches 15.8% at 60 keV/µm with an estimated chromosome aberration detection threshold of 3 Mbp. The dicentric and acentric fragment yields and the dose response of micronuclei formation after proton irradiation were calculated and compared with experimental results.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2026

Model scripts associated with “Revisiting controls on hyporheic respiration with knowledge-guided machine learning at continental scale”

NOTE: The manuscript associated with this data package is currently in review. The data/scripts may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final scripts and additional metadata. This data package is associated with the publication “Revisiting controls on hyporheic respiration with knowledge-guided machine learning at continental scale” submitted to Environmental Science & Technology (Zheng et al. 2026). The project combines mechanistic process modeling with knowledge-guided machine learning (KGML) to evaluate how organic matter chemistry, microbial biomass, and physical substrate accessibility regulate realized respiration rates across river corridors. All data used in this paper have been previously published and can be accessed at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1729719 (Goldman et al., 2020). This data package contains 3 R-markdown (Rmd) preprocessing scripts for the previously published data and subsequent modelling workflows. The full workflow with input and output data can be found in the associated GitHub repository at https://github.com/jianqiuz/KGML-WHONDRS.

Biogeochemistry↗
DOE OSTI2022

Demonstrate Improved Ag Diffusion and Describe the Basis for Pd Penetration Modeling in SiC

In past work, an effective diffusivity coefficient was determined for Ag transport through the silicon carbide layer of a tristructural isotropic fuel particle. The effective diffusivity coefficient accounts for the microstructure of the silicon carbide and includes both bulk diffusion and grain boundary diffusion of Ag. In this report, the model has been improved by accounting for the enhanced concentration of vacancies in the bulk due to irradiation, which substantially influence bulk diffusivity at low temperatures. To improve the BISON model and make it fission rate dependent, effective diffusivity calculations have been performed that incorporate the radiation modified bulk diffusivity. The microstructure and irradiation-dependent effective diffusivity has also been implemented into BISON, and its predictions for Ag release from tristructural isotropic fuel have been successfully compared to AGR-1 post irradiation measurements. Moreover, a new feature has been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) to account for different grain boundary types. The Ag diffusivity in 5 (210)/[001] grain boundaries has been computed and was found to be greater than in random high-angle grain boundaries. The presence of the fission product Pd can also have an important effect on the properties of the silicon carbide layer in tristructural isotropic particles. The penetration of Pd into the silicon carbide layer causes a corrosion reaction that can lead to the failure of the silicon carbide layer; however, this corrosion reaction is not well understood. To enable an improved understanding of the mechanism, ab-initio molecular dynamics simulations of Pd interaction with bulk silicon carbide have been performed. The improved understanding of the reaction will form a basis for future improvements to the BISON’s Pd penetration failure model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Documentation, assessment, and improvements of fission gas modeling capabilities in BISON

Fission gas behavior critically affects nuclear fuel performance, influencing thermal conductivity, microstructure, swelling, and fuel-cladding interactions. The Simple Integrated Fission Gas Release and Swelling (Sifgrs) model, implemented within the BISON fuel performance code, provides advanced capabilities for modeling these phenomena. This manuscript reviews and assesses Sifgrs’ modular structure and submodels describing intragranular, intergranular, and release mechanisms. The modular implementation enables rapid improvements as our understanding advances and new fuel forms emerge. After an overview of capabilities, we present a comprehensive validation suite, comparing predictions to experimental data from tests under normal and transient conditions, including reproduction of the Vitanza curve. Results demonstrate Sifgrs’ accuracy, robustness, and applicability. Current limitations and ongoing efforts are discussed, including mechanistic, multiscale modeling, uncertainty quantification, and validation within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, supporting the safe and efficient use of nuclear fuels.

BISON↗
NASA NTRS1993

Laboratory and theoretical models of planetary-scale instabilities and waves

Research work is proceeding in theoretical, numerical, and experimental geophysical fluid dynamics leading up to a reflight of the GFFC (Geophysical Fluid Flow Cell Experiment) on USML-2. The work is intended not only to generate ideas for future space experiments, but to provide fundamental results concerned with nonlinear and chaotic properties of thermal convection and baroclinic waves in terrestrial and planetary atmospheres. The major efforts are focussed on thermal convection in a rapidly rotating annulus relevant to Jovian atmospheric dynamics, and on the chaotic behavior of baroclinic waves relevant to the Earth's atmosphere. The approach, in preparation for USML-2, is primarily theoretical and numerical. Mechanistic process models are solved numerically in order to identify physical mechanisms that may be observed in the GFFC, and which are important in real geophysical applications. The results from numerical simulations of geophysical fluid flow (subject to rotation and stratification) are compared with previous GFFC experiments on Spacelab-3 and with existing and proposed terrestrial laboratory experiments of various types. Pattern recognition algorithms have been employed to generate low-dimensional descriptions of the highly nonlinear and turbulent numerical simulations. Such empirically truncated descriptions provide for simplified but robust physical interpretations of the dynamics, as well as yielding highly efficient computations of these chaotic flows.

Hart, John E.↗