SEARCH · Engineering Papers

Results for “model parameterization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A data-driven multiscale model for reactive wetting simulations

Here, we describe a data-driven, multiscale technique to model reactive wetting of a silver–aluminum alloy on a Kovar™ (Fe-Ni-Co alloy) surface. We employ molecular dynamics simulations to elucidate the dependence of surface tension and wetting angle on the drop’s composition and temperature. A design of computational experiments is used to efficiently generate training data of surface tension and wetting angle from a limited number of molecular dynamics simulations. The simulation results are used to parameterize models of the material’s wetting properties and compute the uncertainty in the models due to limited data. The data-driven models are incorporated into an engineering-scale (continuum) model of a silver–aluminum sessile drop on a Kovar™ substrate. Model predictions of the wetting angle are compared with experiments of pure silver spreading on Kovar™ to quantify the model-form errors introduced by the limited training data versus the simplifications inherent in the molecular dynamics simulations. The paper presents innovations in the determination of “convergence” of noisy MD simulations before they are used to extract the wetting angle and surface tension, and the construction of their models which approximate physio-chemical processes that are left unresolved by the engineering-scale model. Together, these constitute a multiscale approach that integrates molecular-scale information into continuum scale models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2013

Towards Estimation of HIV-1 Date of Infection: A Time-Continuous IgG-Model Shows That Seroconversion Does Not Occur at the Midpoint between Negative and Positive Tests

Estimating date of infection for HIV-1-infected patients is vital for disease tracking and informed public health decisions, but is difficult to obtain because most patients have an established infection of unknown duration at diagnosis. Previous studies have used HIV-1-specific immunoglobulin G (IgG) levels as measured by the IgG capture BED enzyme immunoassay (BED assay) to indicate if a patient was infected recently, but a time-continuous model has not been available. Therefore, we developed a logistic model of IgG production over time. We used previously published metadata from 792 patients for whom the HIV-1-specific IgG levels had been longitudinally measured using the BED assay. To account for patient variability, we used mixed effects modeling to estimate general population parameters. The typical patient IgG production rate was estimated at r = 6.72[approximate 95% CI 6.17,7.33]x10 -3 OD-n units day -1 , and the carrying capacity at K = 1.84[1.75,1.95] OD-n units, predicting how recently patients seroconverted in the interval Λ t = (31,711) days. Final model selection and validation was performed on new BED data from a population of 819 Swedish HIV-1 patients diagnosed in 2002–2010. On an appropriate subset of 350 patients, the best model parameterization had an accuracy of 94% finding a realistic seroconversion date. We found that seroconversion on average is at the midpoint between last negative and first positive HIV-1 test for patients diagnosed in prospective/cohort studies such as those included in the training dataset. In contrast, seroconversion is strongly skewed towards the first positive sample for patients identified by regular public health diagnostic testing as illustrated in the validation dataset. Our model opens the door to more accurate estimates of date of infection for HIV-1 patients, which may facilitate a better understanding of HIV-1 epidemiology on a population level and individualized prevention, such as guidance during contact tracing.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1989

Superclustering in the explosion scenario

A simple toy model is used to study the spatial distribution of rich clusters in a generic type of explosion scenario. The model, parameterized by the distribution of shell radii and the filling factor, places spherical shells at random and identifies each 'knot' as a cluster. The resulting cluster correlation function is close to a power law extending to the diameter of the largest spheres. Richer clusters form at the intersections of bigger shells and so have stronger correlations. Typical shell radii and filling factors are required to produce the observed number density of clusters. Models with a power-law radius distribution also reproduce the richness distribution of clusters in the Abell catalog. Supercluster multiplicity functions, void probabilities, number counts, topology statistics, and velocity correlations confirm the presence of strong superclustering and quantify the non-Gaussian nature of the model.

Weinberg, David H.↗

DOE OSTI2019

Uncertainty quantification and propagation in CALPHAD modeling

Design is about making decisions bounded by a quantifiable degree of certainty. In the context of alloy design, Integrated Computational Materials Engineering (ICME) provides the framework whereby performance requirements are ultimately transformed into alloy/processing specifications through the combination of (complex) computational models connecting process-structure-property-performance relationships and experiments. Most ICME approaches consider the models used as deterministic and thus do not provide the means to make alloy design decisions with proper confidence measures. At the root of ICME lie CALPHAD models that describe the thermodynamics and phase stability of phases under specific thermodynamic boundary conditions. To date, the vast majority of efforts within the CALPHAD community have been deterministic in that thermodynamic models and the resulting thermodynamic properties and phase diagram features do not explicitly account for the uncertainties inherent in the model formulation or in the experimental/computational data used. In this contribution, we provide an overview of the state of the field. We review major efforts thus far and we then provide a (brief) tutorial on basic concepts of uncertainty quantification and propagation (UQ/UP) in CALPHAD. We discuss the major features of frequentist and Bayesian interpretations of uncertainty and proceed with a discussion of recent case studies in which UQ has been used to parameterize models for the thermodynamic properties of phases. We follow our discussion by presenting frameworks and demonstrating the propagation of uncertainty in thermodynamic properties and phase diagram predictions and briefly discuss how we can use Bayesian frameworks for rigorous model selection as well as for model fusion. Furthermore, we close our contribution by providing context for what has been done and what remains to be accomplished in order to fully embrace the management of uncertainty in CALPHAD modeling, a foundational element of ICME.

42 ENGINEERING↗

NASA NTRS2022

Potential Link Between Ice Nucleation and Climate Model Spread in Arctic Amplification

Arctic amplification (AA) is simulated by all global climate models, however the spread in the degree of projected warming is large and the underlying mechanisms driving it are poorly understood. The impact of the temperature dependence of immersion freezing on cloud feedbacks and AA is studied using NASA's GEOS-5 model. Parameterizations that exhibit low ice-nucleating particle (INP) concentrations in the high Arctic during summer are found to weaken the cloud-phase feedback. This allows sunlight to readily melt sea-ice in the summer, which decreases the stability of the lower troposphere, causing a decrease in wintertime cloud fraction over open ocean. Arctic amplification was found to span from ∼1.4 to >2.6, which spans 30% of the spread in AA in the coupled model intercomparison project models, depending on the temperature dependence of immersion freezing. These results suggest that summertime INP concentrations may provide an observational constraint on AA.

Ice Nucleation↗

DOE OSTI2023

Quantification of early nonpharmaceutical interventions aimed at slowing transmission of Coronavirus Disease 2019 in the Navajo Nation and surrounding states (Arizona, Colorado, New Mexico, and Utah)

During an early period of the Coronavirus Disease 2019 (COVID-19) pandemic, the Navajo Nation, much like New York City, experienced a relatively high rate of disease transmission. Yet, between January and October 2020, it experienced only a single period of growth in new COVID-19 cases, which ended when cases peaked in May 2020. The daily number of new cases slowly decayed in the summer of 2020 until late September 2020. In contrast, the surrounding states of Arizona, Colorado, New Mexico, and Utah all experienced at least two periods of growth in the same time frame, with second surges beginning in late May to early June. Here, we investigated these differences in disease transmission dynamics with the objective of quantifying the contributions of non-pharmaceutical interventions (NPIs) (e.g., behaviors that limit disease transmission). We considered a compartmental model accounting for distinct periods of NPIs to analyze the epidemic in each of the five regions. We used Bayesian inference to estimate region-specific model parameters from regional surveillance data (daily reports of new COVID-19 cases) and to quantify uncertainty in parameter estimates and model predictions. Our results suggest that NPIs in the Navajo Nation were sustained over the period of interest, whereas in the surrounding states, NPIs were relaxed, which allowed for subsequent surges in cases. Our region-specific model parameterizations allow us to quantify the impacts of NPIs on disease incidence in the regions of interest.

60 APPLIED LIFE SCIENCES↗

NASA NTRSDate not supplied

Implications of using MERRA-2 cloud water in aerosol wet scavenging for Pb-210 surface concentrations and vertical profiles in GEOS-Chem

Precipitation scavenging is the dominant loss process for a whole suite of aerosols but model parameterizations of this process are highly uncertain, substantially contributing to large uncertainties in the simulated loadings and radiative forcing of aerosols. Lead-210 (210Pb, radioactive half-life of 22.3 years) is produced by radioactive decay of soil-emitted gaseous 222Rn. It attaches to ambient submicron aerosols and is subject to precipitation scavenging processes. Liu et al. [2001] estimated the global mean lifetime of tropospheric 210Pb aerosols to be ~9 days using the GEOS-Chem model. Luo et al. [2019, 2020] provided an option to use MERRA-2 cloud water in the model aerosol scavenging scheme. In this study, we examine the impact of this revision on 210Pb simulations and compare results with climatological observations of 210Pb surface and UT/LS concentrations, as well as aircraft profile measurements from eleven NASA field missions.

Bo Zhang↗

NASA NTRS2000

NASA/GSFC Research Activities for the Global Ocean Carbon Cycle: A Prospectus for the 21st Century

There are increasing concerns that anthropogenic inputs of carbon dioxide into the Earth system have the potential for climate change. In response to these concerns, the GSFC Laboratory for Hydrospheric Processes has formed the Ocean Carbon Science Team (OCST) to contribute to greater understanding of the global ocean carbon cycle. The overall goals of the OCST are to: 1) detect changes in biological components of the ocean carbon cycle through remote sensing of biooptical properties, 2) refine understanding of ocean carbon uptake and sequestration through application of basic research results, new satellite algorithms, and improved model parameterizations, 3) develop and implement new sensors providing critical missing environmental information related to the oceanic carbon cycle and the flux of CO2 across the air-sea interface. The specific objectives of the OCST are to: 1) establish a 20-year time series of ocean color, 2) develop new remote sensing technologies, 3) validate ocean remote sensing observations, 4) conduct ocean carbon cycle scientific investigations directly related to remote sensing data, emphasizing physiological, empirical and coupled physical/biological models, satellite algorithm development and improvement, and analysis of satellite data sets. These research and mission objectives are intended to improve our understanding of global ocean carbon cycling and contribute to national goals by maximizing the use of remote sensing data.

Gregg, W. W.↗

DOE OSTI2022

Parameterization of generic positive sequence models to represent behavior of inverter based resources in low short circuit scenarios

The connection of large inverter based resources (IBR) in transmission systems is often located geographically and electrically far away from load centers. This, coupled with the displacement of synchronous machine plants, results in a reduction of the network short circuit strength at the point of connection. Under these conditions, state-of-the-art positive sequence simulation platforms and models can have difficulties maintaining numerical stability and/or providing an accurate representation of IBR plant dynamic behavior. As a result, computationally heavy time domain electromagnetic transient (EMT) simulations may be required to evaluate these systems. In this paper, a recently developed improved generic positive sequence model is parameterized to provide a representation of IBR behavior in low short circuit networks. Further, comparisons against generic and detailed EMT models demonstrate the suitability of the improved positive sequence model to study practical stability issues experienced with presently in-service plants. Such a model can provide some of the accuracy of an EMT representation with a much lower computational burden. The performance of the positive sequence model is validated against the behavior shown by both open white box and closed black box EMT domain models from around the world.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Calibrating constitutive models with full‐field data via physics informed neural networks

Abstract The calibration of solid constitutive models with full‐field experimental data is a long‐standing challenge, especially in materials that undergo large deformations. In this paper, we propose a physics‐informed deep‐learning framework for the discovery of hyperelastic constitutive model parameterizations given full‐field surface displacement data and global force‐displacement data. Contrary to the majority of recent literature in this field, we work with the weak form of the governing equations rather than the strong form to impose physical constraints upon the neural network predictions. The approach presented in this paper is computationally efficient, suitable for irregular geometric domains, and readily ingests displacement data without the need for interpolation onto a computational grid. A selection of canonical hyperelastic material models suitable for different material classes is considered including the Neo–Hookean, Gent, and Blatz–Ko constitutive models as exemplars for general non‐linear elastic behaviour, elastomer behaviour with finite strain lock‐up, and compressible foam behaviour, respectively. We demonstrate that physics informed machine learning is an enabling technology and may shift the paradigm of how full‐field experimental data are utilized to calibrate constitutive models under finite deformations.

Hamel, Craig M.↗

NASA NTRS1998

FIRE Arctic Clouds Experiment

An overview is given of the First ISCCP Regional Experiment (FIRE) Arctic Clouds Experiment that was conducted in the Arctic during April through July, 1998. The principal goal of the field experiment was to gather the data needed to examine the impact of arctic clouds on the radiation exchange between the surface, atmosphere, and space, and to study how the surface influences the evolution of boundary layer clouds. The observations will be used to evaluate and improve climate model parameterizations of cloud and radiation processes, satellite remote sensing of cloud and surface characteristics, and understanding of cloud-radiation feedbacks in the Arctic. The experiment utilized four research aircraft that flew over surface-based observational sites in the Arctic Ocean and Barrow, Alaska. In this paper we describe the programmatic and science objectives of the project, the experimental design (including research platforms and instrumentation), conditions that were encountered during the field experiment, and some highlights of preliminary observations, modelling, and satellite remote sensing studies.

Curry, J. A.↗

DOE OSTI2017

Modeling Virus Coinfection to Inform Management of Maize Lethal Necrosis in Kenya

Maize lethal necrosis (MLN) has emerged as a serious threat to food security in sub-Saharan Africa. MLN is caused by coinfection with two viruses, Maize chlorotic mottle virus and a potyvirus, often Sugarcane mosaic virus . To better understand the dynamics of MLN and to provide insight into disease management, we modeled the spread of the viruses causing MLN within and between growing seasons. The model allows for transmission via vectors, soil, and seed, as well as exogenous sources of infection. Following model parameterization, we predict how management affects disease prevalence and crop performance over multiple seasons. Resource-rich farmers with large holdings can achieve good control by combining clean seed and insect control. However, crop rotation is often required to effect full control. Resource-poor farmers with smaller holdings must rely on rotation and roguing, and achieve more limited control. For both types of farmer, unless management is synchronized over large areas, exogenous sources of infection can thwart control. As well as providing practical guidance, our modeling framework is potentially informative for other cropping systems in which coinfection has devastating effects. Finally, our work also emphasizes how mathematical modeling can inform management of an emerging disease even when epidemiological information remains scanty.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2015

Alteration of forest succession and carbon cycling under elevated CO 2

Abstract Regenerating forests influence the global carbon (C) cycle, and understanding how climate change will affect patterns of regeneration and C storage is necessary to predict the rate of atmospheric carbon dioxide (CO 2 ) increase in future decades. While experimental elevation of CO 2 has revealed that young forests respond with increased productivity, there remains considerable uncertainty as to how the long‐term dynamics of forest regrowth are shaped by elevated CO 2 (eCO 2 ). Here, we use the mechanistic size‐ and age‐ structured Ecosystem Demography model to investigate the effects of CO 2 enrichment on forest regeneration, using data from the Duke Forest Free‐Air Carbon dioxide Enrichment (FACE) experiment, a forest chronosequence, and an eddy‐covariance tower for model parameterization and evaluation. We find that the dynamics of forest regeneration are accelerated, and stands consistently hit a variety of developmental benchmarks earlier under eCO 2 . Because responses to eCO 2 varied by plant functional type, successional pathways, and mature forest composition differed under eCO 2 , with mid‐ and late‐successional hardwood functional types experiencing greater increases in biomass compared to early‐successional functional types and the pine canopy. Over the simulation period, eCO 2 led to an increase in total ecosystem C storage of 9.7 Mg C ha ‐1 . Model predictions of mature forest biomass and ecosystem–atmosphere exchange of CO 2 and H 2 O were sensitive to assumptions about nitrogen limitation; both the magnitude and persistence of the ecosystem response to eCO 2 were reduced under N limitation. In summary, our simulations demonstrate that eCO 2 can result in a general acceleration of forest regeneration while altering the course of successional change and having a lasting impact on forest ecosystems.

Miller, Adam D.↗

DOE OSTI2015

Comparison of the kinetics of different Markov models for ligand binding under varying conditions

We recently derived a Markov model for macromolecular ligand binding dynamics from few physical assumptions and showed that its stationary distribution is the grand canonical ensemble [J. W. R. Martini, M. Habeck, and M. Schlather, J. Math. Chem. 52, 665 (2014)]. The transition probabilities of the proposed Markov process define a particular Glauber dynamics and have some similarity to the Metropolis-Hastings algorithm. Here, we illustrate that this model is the stochastic analog of (pseudo) rate equations and the corresponding system of differential equations. Moreover, it can be viewed as a limiting case of general stochastic simulations of chemical kinetics. Thus, the model links stochastic and deterministic approaches as well as kinetics and equilibrium described by the grand canonical ensemble. We demonstrate that the family of transition matrices of our model, parameterized by temperature and ligand activity, generates ligand binding kinetics that respond to changes in these parameters in a qualitatively similar way as experimentally observed kinetics. In contrast, neither the Metropolis-Hastings algorithm nor the Glauber heat bath reflects changes in the external conditions correctly. Both converge rapidly to the stationary distribution, which is advantageous when the major interest is in the equilibrium state, but fail to describe the kinetics of ligand binding realistically. To simulate cellular processes that involve the reversible stochastic binding of multiple factors, our pseudo rate equation model should therefore be preferred to the Metropolis-Hastings algorithm and the Glauber heat bath, if the stationary distribution is not of only interest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

The physiological basis for estimating photosynthesis from Chl>em>a fluorescence

We report solar-induced Chl fluorescence (SIF) offers the potential to curb large uncertainties in the estimation of photosynthesis across biomes and climates, and at different spatiotemporal scales. However, it remains unclear how SIF should be used to mechanistically estimate photosynthesis. In this study, we built a quantitative framework for the estimation of photosynthesis, based on a mechanistic light reaction model with the Chla fluorescence of Photosystem II (SIFPSII) as an input (MLR-SIF). Utilizing 29 C 3 and C 4 plant species that are representative of major plant biomes across the globe, we confirmed the validity of this framework at the leaf level. The MLR-SIF model is capable of accurately reproducing photosynthesis for all C 3 and C 4 species under diverse light, temperature, and CO 2 conditions. We further tested the robustness of the MLR-SIF model using Monte Carlo simulations, and found that photosynthesis estimates were much less sensitive to parameter uncertainties relative to the conventional Farquhar, von Caemmerer, Berry (FvCB) model because of the additional independent information contained in SIF PSII . Once inferred from direct observables of SIF, SIF PSII provides ‘parameter savings’ to the MLR-SIF model, compared to the mechanistically equivalent FvCB model, and thus avoids the uncertainties arising as a result of imperfect model parameterization. Our findings set the stage for future efforts to employ SIF mechanistically to improve photosynthesis estimates across a variety of scales, functional groups, and environmental conditions.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2003

Regional-Scale Modeling at NASA Goddard Space Flight Center

Over the past decade, the Goddard Mesoscale Modeling and Dynamics Group has used a popular regional scale model, MM5, to study precipitation processes. Our group is making contributions to the MM5 by incorporating the following physical and numerical packages: improved Goddard cloud processes, a land processes model (Parameterization for Land-Atmosphere-Cloud Exchange - PLACE), efficient but sophisticated radiative processes, conservation of hydrometeor mass (water budget), four-dimensional data assimilation for rainfall, and better computational methods for trace gas transport. At NASA Goddard, the MM5 has been used to study: (1) the impact of initial conditions, assimilation of satellite-derived rainfall, and cumulus parameterizations on rapidly intensifying oceanic cyclones, hurricanes and typhoons, (2) the dynamic and thermodynamic processes associated with the development of narrow cold frontal rainbands, (3) regional climate and water cycles, (4) the impact of vertical transport by clouds and lightning on trace gas distributiodproduction associated with South and North American mesoscale convective systems, (5) the development of a westerly wind burst (WWB) that occurred during the TOGA COARE and the diurnal variation of precipitation in the tropics, (6) a Florida sea breeze convective event and a Mid-US flood event using a sophisticated land surface model, (7) the influence of soil heterogeneity on land surface energy balance in the southwest GCIP region, (8) explicit simulations (with 1.33 to 4 km horizontal resolution) of hurricanes Bob (1991) and Bonnie (1998), (9) a heavy precipitation event over Taiwan, and (10) to make real time forecasts for a major NASA field program. In this paper, the modifications and simulated cases will be described and discussed.

Tao, W.-K.↗

NASA NTRS2001

Risk analysis: divergent models and convergent interpretations

Material presented at a NASA-sponsored workshop on risk models for exposure conditions relevant to prolonged space flight are described in this paper. Analyses used mortality data from experiments conducted at Argonne National Laboratory on the long-term effects of external whole-body irradiation on B6CF1 mice by 60Co gamma rays and fission neutrons delivered as a single exposure or protracted over either 24 or 60 once-weekly exposures. The maximum dose considered was restricted to 1 Gy for neutrons and 10 Gy for gamma rays. Proportional hazard models were used to investigate the shape of the dose response at these lower doses for deaths caused by solid-tissue tumors and tumors of either connective or epithelial tissue origin. For protracted exposures, a significant mortality effect was detected at a neutron dose of 14 cGy and a gamma-ray dose of 3 Gy. For single exposures, radiation-induced mortality for neutrons also occurred within the range of 10-20 cGy, but dropped to 86 cGy for gamma rays. Plots of risk relative to control estimated for each observed dose gave a visual impression of nonlinearity for both neutrons and gamma rays. At least for solid-tissue tumors, male and female mortality was nearly identical for gamma-ray exposures, but mortality risks for females were higher than for males for neutron exposures. As expected, protracting the gamma-ray dose reduced mortality risks. Although curvature consistent with that observed visually could be detected by a model parameterized to detect curvature, a relative risk term containing only a simple term for total dose was usually sufficient to describe the dose response. Although detectable mortality for the three pathology end points considered typically occurred at the same level of dose, the highest risks were almost always associated with deaths caused by tumors of epithelial tissue origin.

Non-NASA Center↗

DOE OSTI2020

Signatures of Velocity-dependent Dark Matter Self-interactions in Milky Way-mass Halos

In this work, we explore the impact of elastic, anisotropic, velocity-dependent dark matter (DM) self-interactions on the host halo and subhalos of Milky Way (MW)--mass systems. We consider a generic self-interacting dark matter (SIDM) model parameterized by the masses of a light mediator and the DM particle. The ratio of these masses, w, sets the velocity scale above which momentum transfer due to DM self-interactions becomes inefficient. We perform high-resolution zoom-in simulations of an MW-mass halo for values of w that span scenarios in which self-interactions either between the host and its subhalos or only within subhalos efficiently transfer momentum, and we study the effects of self-interactions on the host halo and on the abundance, radial distribution, orbital dynamics, and density profiles of subhalos in each case. The abundance and properties of surviving subhalos are consistent with being determined primarily by subhalo--host halo interactions. In particular, subhalos on radial orbits in models with larger values of the cross section at the host halo velocity scale are more susceptible to tidal disruption owing to mass loss from ram pressure stripping caused by self-interactions with the host. This mechanism suppresses the abundance of surviving subhalos relative to collisionless DM simulations, with stronger suppression for larger values of w. Thus, probes of subhalo abundance around MW-mass hosts can be used to place upper limits on the self-interaction cross section at velocity scales of ~200 km s-1, and combining these measurements with the orbital properties and internal dynamics of subhalos may break degeneracies among velocity-dependent SIDM models.

79 ASTRONOMY AND ASTROPHYSICS↗