Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stefan model”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Verification of Data-Driven Models of Physical Phenomena using Interpretable Approximation

Machine-learned models, specifically neural networks, are increasingly used as “closures” or “constitutive models” in engineering simulators to represent fine-scale physical phenomena that are too computationally expensive to resolve explicitly. However, these neural net models of unresolved physical phenomena tend to fail unpredictably and are therefore not used in mission-critical simulations. In this report, we describe new methods to authenticate them, i.e., to determine the (physical) information content of their training datasets, qualify the scenarios where they may be used and to verify that the neural net, as trained, adhere to physics theory. We demonstrate these methods with neural net closure of turbulent phenomena used in Reynolds Averaged Navier-Stokes equations. We show the types of turbulent physics extant in our training datasets, and, using a test flow of an impinging jet, identify the exact locations where the neural network would be extrapolating i.e., where it would be used outside the feature-space where it was trained. Using Generalized Linear Mixed Models, we also generate explanations of the neural net (à la Local Interpretable Model agnostic Explanations) at prototypes placed in the training data and compare them with approximate analytical models from turbulence theory. Finally, we verify our findings by reproducing them using two different methods.

42 ENGINEERING↗
DOE OSTI2025

The impact of non-local parallel electron transport on plasma-impurity reaction rates in tokamak scrape-off layer plasmas

Abstract Plasma-impurity reaction rates are a crucial part of modelling tokamak scrape-off layer (SOL) plasmas. To avoid calculating the full set of rates for the large number of important processes involved, a set of effective rates are typically derived which assume Maxwellian electrons. However, non-local parallel electron transport may result in non-Maxwellian electrons, particularly close to divertor targets. Here, the validity of using Maxwellian-averaged rates in this context is investigated by computing the full set of rate equations for a fixed plasma background from kinetic and fluid SOL simulations. We consider the effect of the electron distribution as well as the impact of the electron transport model on plasma profiles. Results are presented for lithium, beryllium, carbon, nitrogen, neon and argon. It is found that electron distributions with enhanced high-energy tails can result in significant modifications to the ionisation balance and radiative power loss rates from excitation, on the order of 50%–75% for the latter. Fluid electron models with Spitzer-Härm or flux-limited Spitzer-Härm thermal conductivity, combined with Maxwellian electrons for rate calculations, can increase or decrease this error, depending on the impurity species and plasma conditions. Based on these results, we also discuss some approaches to experimentally observing non-local electron transport in SOL plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2017

Classifying Multi-Model Wheat Yield Impact Response Surfaces Showing Sensitivity to Temperature and Precipitation Change

Crop growth simulation models can differ greatly in their treatment of key processes and hence in their response to environmental conditions. Here, we used an ensemble of 26 process-based wheat models applied at sites across a European transect to compare their sensitivity to changes in temperature (minus 2 to plus 9 degrees Centigrade) and precipitation (minus 50 to plus 50 percent). Model results were analysed by plotting them as impact response surfaces (IRSs), classifying the IRS patterns of individual model simulations, describing these classes and analysing factors that may explain the major differences in model responses. The model ensemble was used to simulate yields of winter and spring wheat at four sites in Finland, Germany and Spain. Results were plotted as IRSs that show changes in yields relative to the baseline with respect to temperature and precipitation. IRSs of 30-year means and selected extreme years were classified using two approaches describing their pattern. The expert diagnostic approach (EDA) combines two aspects of IRS patterns: location of the maximum yield (nine classes) and strength of the yield response with respect to climate (four classes), resulting in a total of 36 combined classes defined using criteria pre-specified by experts. The statistical diagnostic approach (SDA) groups IRSs by comparing their pattern and magnitude, without attempting to interpret these features. It applies a hierarchical clustering method, grouping response patterns using a distance metric that combines the spatial correlation and Euclidian distance between IRS pairs. The two approaches were used to investigate whether different patterns of yield response could be related to different properties of the crop models, specifically their genealogy, calibration and process description. Although no single model property across a large model ensemble was found to explain the integrated yield response to temperature and precipitation perturbations, the application of the EDA and SDA approaches revealed their capability to distinguish: (i) stronger yield responses to precipitation for winter wheat than spring wheat; (ii) differing strengths of response to climate changes for years with anomalous weather conditions compared to period-average conditions; (iii) the influence of site conditions on yield patterns; (iv) similarities in IRS patterns among models with related genealogy; (v) similarities in IRS patterns for models with simpler process descriptions of root growth and water uptake compared to those with more complex descriptions; and (vi) a closer correspondence of IRS patterns in models using partitioning schemes to represent yield formation than in those using a harvest index. Such results can inform future crop modelling studies that seek to exploit the diversity of multi-model ensembles, by distinguishing ensemble members that span a wide range of responses as well as those that display implausible behaviour or strong mutual similarities.

Classification; Climate change; Crop model; Ensemb↗
DOE OSTI2020

Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model

Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.

58 GEOSCIENCES↗
DOE OSTI2020

AATTENUATION—The Atmospheric Attenuation Model for CSP Tower Plants: A Look-Up Table for Operational Implementation

Attenuation of solar radiation between the receiver and the heliostat field in concentrated solar power (CSP) tower plants can reduce the overall system performance significantly. The attenuation varies strongly with time and the average attenuation at different sites might also vary strongly from each other. If no site specific attenuation data is available, the optimal plant design cannot be determined and rough estimations of the attenuation effect are required leading to high uncertainties of yield analysis calculations. The attenuation is caused mainly by water vapor content and aerosol particles in the lower atmospheric layer above ground. Although several on-site measurement systems have been developed during recent years, attenuation data sets are usually not available to be included during the plant project development. An Atmospheric Attenuation (AATTENUATION) model to derive the atmospheric transmittance between a heliostat and receiver on the basis of common direct normal irradiance (DNI), temperature, relative humidity, and barometric pressure measurements was developed and validated by the authors earlier. The model allows the accurate estimation of attenuation for sites with low attenuation and gives an estimation of the attenuation for less clear sites. However, the site-dependent coefficients of the AATTENUATION model had to be developed individually for each site of interest, which required time-consuming radiative transfer simulations, considering the exact location and altitude, as well as the pre-dominant aerosol type at the location. This strongly limited the application of the model despite its typically available input data. In this manuscript, a look-up table (LUT) is presented which enables the application of the AATTENUATION model at the site of interest without the necessity to perform the according complex radiative transfer calculations for each site individually. This enables the application of the AATTENUATION model for virtually all resource assessments for tower plants and in an operational mode in real time within plant monitoring systems around the world. The LUT also facilitates the generation of solar attenuation maps on the basis of long-term meteorological data sets which can be considered during resource assessment for CSP tower plant projects. The LUTs are provided together with this manuscript as supplementary files. The LUT for the AATTENUATION model was developed for a solar zenith angle (SZA) grid of 1°, an altitude grid of 100 m, 7 different standard aerosol types and the standard AFGL atmospheres for mid-latitudes and the tropics. The LUT was tested against the original version of the AATTENUATION model at 4 sites in Morocco and Spain, and it was found that the additional uncertainty introduced by the application of the LUT is negligible. With the information of latitude, longitude, altitude above mean sea level, DNI, relative humidity (RH), ambient temperature (Tair), and barometric pressure (bp), the attenuation can be now derived easily for each site of interest.

14 SOLAR ENERGY↗
DOE OSTI2023

Warming of the Willamette River, 1850–present: the effects of climate change and river system alterations

Abstract. Using archival research methods, we recovered and combined data from multiple sources to produce a unique, 140-year record of daily water temperature (Tw) in the lower Willamette River, Oregon (1881–1890, 1941–present). Additional daily weather and river flow records from the 1850s onwards are used to develop and validate a statistical regression model of Tw for 1850–2020. The model simulates the time-lagged response of Tw to air temperature and river flow and is calibrated for three distinct time periods: the late 19th, mid-20th, and early 21st centuries. Results show that Tw has trended upwards at 1.1 ∘C per century since the mid-19th century, with the largest shift in January and February (1.3 ∘C per century) and the smallest in May and June (∼ 0.8 ∘C per century). The duration that the river exceeds the ecologically important threshold of 20 ∘C has increased by about 20 d since the 1800s, to about 60 d yr−1. Moreover, cold-water days below 2 ∘C have virtually disappeared, and the river no longer freezes. Since 1900, changes are primarily correlated with increases in air temperature (Tw increase of 0.81 ± 0.25 ∘C) but also occur due to alterations in the river system such as depth increases from reservoirs (0.34 ± 0.12 ∘C). Managed release of water affects Tw seasonally, with an average reduction of up to 0.56 ∘C estimated for September. River system changes have decreased variability (σ) in daily minimum Tw by 0.44 ∘C, increased thermal memory, reduced interannual variability, and reduced the response to short-term meteorological forcing (e.g., heat waves). These changes fundamentally alter the response of Tw to climate change, posing additional stressors on fauna.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2021

Methodology to Assess the Changing Risk of Yield Failure Due to Heat and Drought Stress Under Climate Change

While the understanding of average impacts of climate change on crop yields is improving, few assessments have quantified expected impacts on yield distributions and the risk of yield failures. Here we present the relative distribution as a method to assess how the risk of yield failure due to heat and drought stress (measured in terms of return period between yields falling 15% below previous five year Olympic average yield) responds to changes of the underlying yield distributions under climate change. Relative distributions are used to capture differences in the entire yield distribution between baseline and climate change scenarios, and to further decompose them into changes in the location and shape of the distribution. The methodology is applied here for the case of rainfed wheat and grain maize across Europe using an ensemble of crop models under three climate change scenarios with simulations conducted at 25 km resolution. Under climate change, maize generally displayed shorter return periods of yield failures (with changes under RCP 4.5 between −0.3 and 0 years compared to the baseline scenario) associated with a shift of the yield distribution towards lower values and changes in shape of the distribution that further reduced the frequency of high yields. This response was prominent in the areas characterized in the baseline scenario by high yields and relatively long return periods of failure. Conversely, for wheat, yield failures were projected to become less frequent under future scenarios (with changes in the return period of −0.1 to +0.4 years under RCP 4.5) and were associated with a shift of the distribution towards higher values and a change in shape increasing the frequency of extreme yields at both ends. Our study offers an approach to quantify the changes in yield distributions that drive crop yield failures. Actual risk assessments additionally require models that capture the variety of drivers determining crop yield variability and scenario climate input data that samples the range of probable climate variation.

Climate risk assessment↗
NASA NTRS1990

Retrieval of cirrus cloud properties from comparative analyses of aircraft and satellite measurements made during the 1986 FIRE IFO

Results are presented of a comparison of cirrus cloud properties obtained from aircraft measurements made during the FIRE Intensive Field Observations and Landsat 5 and NOAA-9 observations carried out during several near coincident overpasses by these satellites. Results of the analyses of these measurements and a comparison between measurements and model calculations were used to determine particle sizes within cirrus clouds and the IR optical depths.

Hammer, Philip D.↗
NASA NTRS2010

A Model-Based Probabilistic Inversion Framework for Wire Fault Detection Using TDR

Time-domain reflectometry (TDR) is one of the standard methods for diagnosing faults in electrical wiring and interconnect systems, with a long-standing history focused mainly on hardware development of both high-fidelity systems for laboratory use and portable hand-held devices for field deployment. While these devices can easily assess distance to hard faults such as sustained opens or shorts, their ability to assess subtle but important degradation such as chafing remains an open question. This paper presents a unified framework for TDR-based chafing fault detection in lossy coaxial cables by combining an S-parameter based forward modeling approach with a probabilistic (Bayesian) inference algorithm. Results are presented for the estimation of nominal and faulty cable parameters from laboratory data.

Schuet, Stefan R.↗
NASA NTRSDate not supplied

Hybrid Modeling of Unmanned Aerial Vehicle Electric Powertrain for Fault Detection and Diagnostics

This paper shows the application of hybrid physics-informed machine learning to a representative electric powertrain for unmanned aerial vehicles. The model is composed of physics-derived principles and empirical equations, as well as fully connected networks that are strategically placed within the model to substitute equations that are subject to large uncertainty. Polynomial fit driven by heuristics or empirical observations can be substituted by more flexible networks that can minimize the error between model predictions and observations without being restricted to a predefined functional form. This modeling strategy allows training of networks deep inside the model and unknown parameters in a single learning stage. It has already been applied to Li-ion batteries in the past, and in this work, we extend the applications to other components of an electric powertrain, namely electronic speed controller with pulse-width modulation, and brushless DC motor with connected propeller. Training and testing of the model is carried out using experimental data from Li-ion battery discharge and powertrain testing in a laboratory environment.

Physics-Informed Machine Learning↗
DOE OSTI2025

Designing a Framework for Solving Multiobjective Simulation Optimization Problems

Multiobjective simulation optimization (MOSO) problems are optimization problems with multiple conflicting objectives, where evaluation of at least one of the objectives depends on a black-box numerical code or real-world experiment, which we refer to as a simulation. Whereas an extensive body of research is dedicated to developing new algorithms and methods for solving these and related problems, it is challenging and time-consuming to integrate these techniques into real-world production-ready solvers. This is partly because of the diversity and complexity of modern state-of-the-art MOSO algorithms and methods and partly because of the complexity and specificity of many real-world problems and their corresponding computing environments. The complexity of this problem is only compounded when introducing potentially complex and/or domain-specific surrogate-modeling techniques, problem formulations, design spaces, and data acquisition functions. Here, this paper carefully surveys the current state of the art in MOSO algorithms, techniques, and solvers, as well as problem types and computational environments where MOSO is commonly applied. We then present several key challenges in the design of a parallel multiobjective simulation optimization framework (ParMOO) and how they have been addressed. Finally, we provide two case studies demonstrating how customized ParMOO solvers can be quickly built and deployed to solve real-world MOSO problems.

engineering design optimization↗
DOE OSTI2026

Effect of soot concentration on radiative transport in a 2 m JP-8 pool fire

Hubbard et al. (2024) investigated conjugate heat transfer in a 7.9 m JP-8 pool fire, utilizing both heavily sooting and non-sooting fuel models. Their findings indicated that all models produced pool radiative heat fluxes within the measurement uncertainty for large-scale outdoor experiments. This study aims to further investigate the influence of soot concentration on radiative transport. To this end, Large Eddy Simulations (LES) of a 2 m diameter JP-8 pool fire were conducted, employing four turbulent combustion chemistry models: the Eddy Dissipation Concept (EDC) and three variants of flamelet model. Additionally, three variations of the one-equation subgrid scale turbulence kinetic energy model, and multiple mesh resolutions were utilized. The primary validation metrics were pool radiative heat flux and plume radiative fraction, with simulation data compared to experimental results published by Blanchat et al. (2010). The EDC model produced high soot concentrations near the pool surface, resulting in a significant shielding effect. In contrast, the flamelet combustion model has different kinetic limiters for soot production, leading to soot formation higher in the flame, where it can be more readily oxidized. Flamelet models incorporating both acetylene and benzene contributions to soot growth demonstrated better agreement with experimental data. To assess the sensitivity of pool radiative heat flux to soot mass concentration, we used a one-dimensional approximation over a wide range of soot concentrations and compared it to three-dimensional solutions obtained using Discrete Ordinates. Furthermore, this analysis illustrates the impacts of soot shielding and Turbulence Radiation Interactions (TRI), with TRI emission enhancement factors predicted to be in the range of 3–5 near the edges of the plume.

Computational fluid dynamics↗
DOE OSTI2022

Thermodynamics-guided alloy and process design for additive manufacturing

In conventional processing, metals go through multiple manufacturing steps including casting, plastic deformation, and heat treatment to achieve the desired property. In additive manufacturing (AM) the same target must be reached in one fabrication process, involving solidification and cyclic remelting. The thermodynamic and kinetic differences between the solid and liquid phases lead to constitutional undercooling, local variations in the solidification interval, and unexpected precipitation of secondary phases. These features may cause many undesired defects, one of which is the so-called hot cracking. The response of the thermodynamic and kinetic nature of these phenomena to high cooling rates provides access to the knowledge-based and tailored design of alloys for AM. Here, we illustrate such an approach by solving the hot cracking problem, using the commercially important IN738LC superalloy as a model material. The same approach could also be applied to adapt other hot-cracking susceptible alloy systems for AM.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Structure and flexibility of the yeast NuA4 histone acetyltransferase complex

The NuA4 protein complex acetylates histones H4 and H2A to activate both transcription and DNA repair. We report the 3.1-Å resolution cryo-electron microscopy structure of the central hub of NuA4, which flexibly tethers the histone acetyltransferase (HAT) and Trimer Independent of NuA4 involved in Transcription Interactions with Nucleosomes (TINTIN) modules. The hub contains the large Tra1 subunit and a core that includes Swc4, Arp4, Act1, Eaf1, and the C-terminal region of Epl1. Eaf1 stands out as the primary scaffolding factor that interacts with the Tra1, Swc4, and Epl1 subunits and contributes the conserved HSA helix to the Arp module. Using nucleosome-binding assays, we find that the HAT module, which is anchored to the core through Epl1, recognizes H3K4me3 nucleosomes with hyperacetylated H3 tails, while the TINTIN module, anchored to the core via Eaf1, recognizes nucleosomes that have hyperacetylated H2A and H4 tails. Together with the known interaction of Tra1 with site-specific transcription factors, our data suggest a model in which Tra1 recruits NuA4 to specific genomic sites then allowing the flexible HAT and TINTIN modules to select nearby nucleosomes for acetylation.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1992

Cirrus microphysics and radiative transfer - Cloud field study on 28 October 1986

An analysis of remote-sensing measurements which had detected an inhomogeneous cloud structure, at 6-11 km altitude, in the data acquired during the FIRE 86 experiment in a 75 x 50 km cirrus cloud field, has derived fluxes whose comparison with modeled fluxes imply a modeling-based underestimation of both solar reflectivity/attenuation and downward IR fluxes. Reconciliation of model results with measurements can be achieved either by adding large concentrations of ice crystals or by altering the backscattering properties of ice crystals.

Kinne, Stefan↗