SEARCH · Engineering Papers

Results for “equilibrium 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.

SourceResearch connectionFrom yearThrough year

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

At least 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Analytical simulation of meander morphology from equilibrium to long-term evolution: Impacts of channel geometry and vegetation-induced coarsening

This study introduces an innovative approach to modeling meandering river morphology, integrating and investigating the effects of geometric characteristics and vegetation-induced channel coarsening. The developed comprehensive framework combines several advanced techniques: Genetic Programming for refining the scour factor of transverse bed slope, a Leaf Area Index (LAI)-enhanced analytical model for quantifying vegetative flow resistance, and an upstream-weighted moving average method for efficient approximation of the convolution integral in meander migration calculations. The model is validated against both an idealized Kinoshita meander and a natural bend of the Tumen River (China) in equilibrium, demonstrating its robustness across diverse scales and conditions. The model's ability to simulate the long-term evolution, including cutoff events, provides valuable insight for river management strategies. The current findings demonstrate that channel geometry, particularly width-to-depth ratio, plays a dominant role in meander evolution, with wider channels prone to more complex and rapid morphological changes. Vegetation effects are most pronounced in channels with moderate width-to-depth ratios, where they can significantly influence migration rates and bed topography. A combination of channel widening and deepening, coupled with strategic vegetation management, can effectively enhance navigability while maintaining channel stability in the studied Tumen River reach. Sensitivity analyses highlight the complex interplay between hydraulic conditions, sediment characteristics, and vegetation in shaping river morphology. This research advances understanding of the multifaceted nature of meandering river systems and offers practical tools for informed decision-making in river engineering and environmental management, particularly in the context of climate change and increasing anthropogenic pressures on fluvial ecosystems.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2004

Condensation Front Migration in a Protoplanetary Nebula

Condensation front dynamics are investigated in the mid-solar nebula region. A quasi-steady model of the evolving nebula is combined with equilibrium vapor pressure curves to determine evolutionary condensation fronts for selected species. These fronts are found to migrate inwards from the far-nebula to final positions during a period of 10(exp 7) years. The physical process governing this movement is a combination of local viscous heating and luminescent heating from the central star. Two luminescent heating models are used and their effects on the ultimate radial position of the condensation front are discussed. At first the fronts move much faster than the nebular accretion velocity, but after a time the accreting gas and dust overtakes the slowing condensation front.

Davis, Sanford S.↗

NASA NTRSDate not supplied

Predictive Modeling of Carbon Ablators

Efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. To reduce the need for extensive testing, accelerate the design cycle process, and reduce uncertainty margins applied to final designs, NASA is developing simulation and modeling tools that enable characterization of material properties and response to high-enthalpy environments. The Porous Microstructure Analysis (PuMA) code has been developed for computing macroscale (volume averaged) properties of porous materials using microscale images from micro-computed tomography (micro-CT). Microscale modeling requires a realistic representation of a material microstructure; these are obtained either synthetically during the design of the material or through X-ray micro-CT. Volume averaged properties are then used to inform macroscale material response models, such as those implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO) software, also actively developed by NASA. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale assuming local thermal equilibrium. These tools were developed to efficiently interface with other pre-existing codes such as SPARTA (direct simulation Monte Carlo), DPLR (hypersonic CFD), NEQAIR (radiative transport) and DAKOTA (uncertainty quantification and optimization). Detailed flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI]) is critical for validating these computational tools for NASA applications. Examples of modeling ablative material response using these codes will be presented including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrate the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains, through the use of massively parallel computations.

Thermal Protection Systems↗

NASA NTRSDate not supplied

Predictive Modeling of Carbon Ablators Using Micro and Macro-Scale Modeling

Thermal Protection Systems↗

DOE OSTI2020

Active Learning A Neural Network Model For Gold Clusters & Bulk From Sparse First Principles Training Data

Small metal clusters are of fundamental scientific interest and of tremendous significance in catalysis. These nanoscale clusters display diverse geometries and structural motifs depending on the cluster size; a knowledge of this size-dependent structural motifs and their dynamical evolution has been of longstanding interest. Given the high computational cost of first-principles calculations, molecular modeling and atomistic simulations such as molecular dynamics (MD) has proven to be an important complementary tool to aid this understanding. Classical MD typically employ predefined functional forms which limits their ability to capture such complex size-dependent structural and dynamical transformation. Neural Network (NN) based potentials represent flexible alternatives and in principle, well-trained NN potentials can provide high level of flexibility, transferability and accuracy on-par with the reference model used for training. A major challenge, however, is that NN models are interpolative and requires large quantities (similar to 10 4 or greater) of training data to ensure that the model adequately samples the energy landscape both near and far-from-equilibrium. A highly desirable goal is minimize the number of training data, especially if the underlying reference model is first-principles based and hence expensive. In this work, we introduce an active learning (AL) scheme that trains a NN model on-the-fly with minimal amount of first-principles based training data. Our AL workflow is initiated with a sparse training dataset (similar to 1 to 5 data points) and is updated on-the-fly via a Nested Ensemble Monte Carlo scheme that iteratively queries the energy landscape in regions of failure and updates the training pool to improve the network performance. Using a representative system of gold clusters, we demonstrate that our AL workflow can train a NN with similar to 500 total reference calculations. Using an extensive DFT test set of similar to 1100 configurations, we show that our AL-NN is able to accurately predict both the DFT energies and the forces for clusters of a myriad of different sizes. Our NN predictions are within 30 meV/atom and 40 meV/angstrom of the reference DFT calculations. Moreover, our AL-NN model also adequately captures the various size-dependent structural and dynamical properties of gold clusters in excellent agreement with DFT calculations and available experiments. We finally show that our AL-NN model also captures bulk properties reasonably well, even though they were not included in the training data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Phase-field modeling of radiation-induced segregation for multicomponent alloys

Structural alloys under irradiation are known to undergo radiation-induced solute redistribution (RIS) at grain boundaries, leading to detrimental effects such as intergranular corrosion and stress-assisted cracking. To better understand the phenomenon, improved models of RIS applicable to concentrated, multicomponent alloys, and mesoscale microstructures are needed. In this talk, we present a novel grand-potential-based phase-field model to account for the complete set of multicomponent kinetic and thermodynamic couplings between atoms and point defects in the Onsager transport equations. We demonstrate multiscale modeling capability by deriving the Onsager coefficient matrix from atomistic-based Kinetic Monte Carlo simulations. Model predictions and validations of RIS and the effect of defect production, grain boundary sink strength and density are demonstrated for a model FCC FeCrNi system. We also demonstrate the novel capability to describe RIS in the presence of equilibrium segregation described using a density-based CALPHAD thermodynamic framework. This modeling approach overcomes certain limitations in conventional RIS models and provides a step closer towards multiscale modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1981

Structure of the solar chromosphere. III - Models of the EUV brightness components of the quiet-sun

The described investigation is concerned with the solution of the non-LTE optically thick transfer equations for hydrogen, carbon, and other constituents to determine semiempirical models for six components of the quiet solar chromosphere. For a given temperature-height distribution, the solution is obtained of the equations of statistical equilibrium, radiative transfer for lines and continua, and hydrostatic equilibrium to find the ionization and excitation conditions for each atomic constituent. The emergent spectrum is calculated, and a trial and error approach is used to adjust the temperature distribution so that the emergent spectrum is in best agreement with the observed one. The relationship between semiempirical models determined in this way and theoretical models based on radiative equilibrium is discussed by Avrett (1977). Harvard Skylab EUV observations are used to determine models for a number of quiet-sun regions.

Vernazza, J. E.↗

NASA NTRS1992

Evolution of a terrestrial magma ocean: Thermodynamics, kinetics, rheology, convection, differentiation

The evolution of an initially totally molten magma ocean is constrained on the basis of analysis of various physical problems in the magma ocean. First of all an equilibrium thermodynamics of the magma ocean is developed in the melting temperature range. The equilibrium thermodynamical parameters are found as functions only of temperature and pressure and are used in the subsequent models of kinetics and convection. Kinematic processes determine the crystal size and also determine a non-equilibrium thermodynamics of the system. Rheology controls all dynamical regimes of the magma ocean. The thermal convection models for different rheological laws are developed for both the laminar convection and for turbulent convection in the case of equilibrium thermodynamics of the multiphase system. The evolution is estimated on the basis of all the above analysis.

Solomatov, V. S.↗

NASA NTRS2023

Advancing Solar and Heliospheric Science Through the Ongoing Development and Support of Atomic and Laboratory Plasma Physics

This paper outlines the necessity for the availability, accessibility, and expansion of atomic physics data and analysis tools for the meaningful interpretation of spectroscopic and polarimetric observations. As we move towards observing the Sun at higher spatio-temporal resolutions, and near-continuously at a range of wavelengths, it becomes critical to develop the appropriate atomic data and physics tools to facilitate scientific progress. We recommend the continued improvement and expansion of current databases to support the development of optically-thick/radiative transfer models, evaluate non-thermal and non-equilibrium ionization effects, and quantify uncertainties in atomic and molecular values. A critical longterm goal will require extending and strengthening collaborations across the atomic, solar/heliospheric, and laboratory plasma physics communities through the participation and training of early career scientists. We also recommend establishing funding for a centralized atomic physics resource made up of a comprehensive and user-oriented atomic database and modeling framework.

Y J Rivera↗

DOE OSTI2024

Sea-surface temperature pattern effects have slowed global warming and biased warming-based constraints on climate sensitivity

The observed rate of global warming since the 1970s has been proposed as a strong constraint on equilibrium climate sensitivity (ECS) and transient climate response (TCR)—key metrics of the global climate response to greenhouse-gas forcing. Using CMIP5/6 models, we show that the inter-model relationship between warming and these climate sensitivity metrics (the basis for the constraint) arises from a similarity in transient and equilibrium warming patterns within the models, producing an effective climate sensitivity (EffCS) governing recent warming that is comparable to the value of ECS governing long-term warming under CO 2 forcing. However, CMIP5/6 historical simulations do not reproduce observed warming patterns. When driven by observed patterns, even high ECS models produce low EffCS values consistent with the observed global warming rate. The inability of CMIP5/6 models to reproduce observed warming patterns thus results in a bias in the modeled relationship between recent global warming and climate sensitivity. Correcting for this bias means that observed warming is consistent with wide ranges of ECS and TCR extending to higher values than previously recognized. These findings are corroborated by energy balance model simulations and coupled model (CESM1-CAM5) simulations that better replicate observed patterns via tropospheric wind nudging or Antarctic meltwater fluxes. Because CMIP5/6 models fail to simulate observed warming patterns, proposed warming-based constraints on ECS, TCR, and projected global warming are biased low. The results reinforce recent findings that the unique pattern of observed warming has slowed global-mean warming over recent decades and that how the pattern will evolve in the future represents a major source of uncertainty in climate projections.

58 GEOSCIENCES↗

NASA NTRS1994

Turbulence modeling for separated flow

Two projects are described in this report. The first involves assessing turbulence models in separated flow. The second addresses the anomalous behavior of certain turbulence models in stagnation point flow. The primary motivation for developing turbulent transport models is to provide tools for computing non-equilibrium, or complex, turbulent flows. Simple flows can be analyzed using data correlations or algebraic eddy viscosities, but in more complicated flows such as a massively separated boundary layer, a more elaborate level of modeling is required. It is widely believed that at least a two-equation transport model is required in such cases. The transport equations determine the evolution of suitable velocity and time-scales of the turbulence. The present study included assessment of second-moment closures in several separated flows, including sharp edge separation; smooth wall, pressure driven separation; and unsteady vortex shedding. Flows with mean swirl are of interest for their role in enhancing mixing both by turbulent and mean motion. The swirl can have a stabilizing effect on the turbulence. An axi-symmetric extension to the INS-2D computer program was written adding the capability of computing swirling flow. High swirl can produce vortex breakdown on the centerline of the jet and it occurs in various combustors.

Durbin, Paul A.↗

NASA NTRSDate not supplied

Unified Material-Environment Interaction Model for Binary UHTC Composites

A new model for the surface chemistry of zirconium and hafnium diborides containing silicon-carbide when exposed to high temperature conditions in air is formulated. Oxidation dynamics are simplified by considering the diffusion-limited equilibrium regime within the material at elevated temperatures. This model is first assessed against experimental mass-gain data taken from UHTC samples heated to high temperature in an oxygen environment. Then, further evaluations of the new model are conducted using hypersonic CFD simulations to analyze conditions experienced by HfB2-SiC in arc jet experiments conducted at NASA Ames Research Center. Measurements of stagnation point heat transfer and pressure are then used to calibrate the simulations. Results of coupled CFD-material response simulations are then compared to the surface temperatures measured during the Ames arc jet test. Results support the use of the proposed model while highlighting the need for improved gas-phase data on the additional constituents considered.

ultra high temperature ceramic↗

DOE OSTI2023

Analytical population dynamics underlying harmful algal blooms triggered by prey avoidance

A modified version of the NPZD ecosystem model is used to analytically examine the effects of predation avoidance, a possible mechanism for triggering harmful algal blooms (HAB). To resolve HAB development caused by predation avoidance, an additional phytoplankton functional group is considered, one that has slower nutrient uptake and better predation avoidance characteristics than the non-harmful phytoplankton group used in traditional NPZD models. Because the two phytoplankton groups (one non-harmful and one HAB) compete for only one resource within the same system, steady state (equilibrium) conditions cannot occur without the presence of zooplankton; only the non-harmful phytoplankton group, which defeats the HAB group in the resource competition, can survive in the equilibrium. The presence of sufficient zooplankton effectively acts to replenish the nutrient pool by consuming the non-harmful phytoplankton. When this occurs, two equilibrium states are found: one with both phytoplankton groups coexisting, and one that only includes the HAB group. Further, the condition required for equilibrium is that the total nitrogen within the system should be larger than a threshold determined by model coefficients. The threshold and feasibility of the equilibrium are sensitive to the relative HAB predation avoidance coefficient. If the coefficient is larger than the ratio of net growth rates between the HAB and non-harmful phytoplankton group, the threshold becomes infinite, and an equilibrium is not feasible. The time scale for the system to reach an equilibrium state that includes a HAB group is determined asymptotically. The dependence of a threshold condition as a controlling factor may explain the regime shift of dominant species causing HABs. The ecosystem model is fully implemented into the Regional Ocean Modeling System and applied to an idealized coastal embayment (with depths and geometry taken from San Francisco Bay) to show numerically the dominance of prey avoidance dynamics in a natural shallow water environment that includes advection and diffusion. The analytical results improve strategies for HAB modeling and provide guidance for setting model coefficients necessary to resolve a HAB event.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Characteristics of locational uncertainty marginal price for correlated uncertainties of variable renewable generation and demands

With the rapid increase of variable renewable energy sources in power systems, how to manage and price the uncertainty of renewable resources’ power outputs is becoming an urgent issue. Current market designs considering the uncertainties are mainly based on the probabilistic scenario set of demand and renewable energy resources power outputs. This consideration makes market designs vulnerable to three significant challenges when put into practice. First, the accurate probability distribution of renewable generation is hard to obtain in real-time. Second, it is challenging to clear the market timely with many scenarios to guarantee accuracy. Third, generation cost recovery cannot be guaranteed for some scenarios. To overcome these challenges, this paper proposes a locational uncertainty marginal price model to price the uncertainty explicitly based on a scenario-free stochastic market-clearing model. Instead of using the probabilistic scenario set, the uncertainty of renewable energy sources and loads is modeled with distributionally-robust chance constraints. The correlation of uncertainties can be endogenously modeled in both the market-clearing and the locational uncertainty marginal price formation. Furthermore, this paper proves that generation cost recovery, revenue adequacy, and partial market equilibrium can be achieved using the locational uncertainty marginal price model. Numerical results from both the small and large systems simulations validate that the generation cost recovery is maintained no matter the generation participates in uncertainty mitigation or not. The transmission congestion surplus is also allocated appropriately among loads, renewable energy sources, and financial transmission right owners.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

NASA NTRSDate not supplied

Characterization of CO 4th Positive Radiation in CO2-Ar Mixtures

This paper presents measurements and simulations of CO 4th Positive spectra obtained in NASA Ames Research Center's Electric Arc Shock Tube (EAST) and CentraleSupelec's Plasma Torch. The experiments were aimed at measuring optically thin emission for the main band features of CO 4th Positive by using test gas mixtures of CO2 and Ar. By using mixtures of CO2 and Ar, the temperature of the flow (and consequently black-body curve) is increased compared to using pure CO2 (or CO2 / N2 mixtures). This temperature increase allows for emission spectroscopy of features that would nominally be black-body limited for conditions relevant to Mars entry. The CO 4th Positive system is the focus of this study as it accounts for a large percentage of the emitted radiation for Martian entry, and also due to the difficulties of obtaining experimental validation data due to the emission appearing in the Vacuum Ultraviolet (VUV) spectral range. The focus of this paper is to provide a comprehensive comparison between the EAST and plasma torch data along with various CO 4th Positive databases available in the literature. Previously reported VUV CO2 / N2 EAST data required a larger temperature than would be predicted based upon equilibrium shock relations. The need for the temperature increase implies that the EAST experiment did not fully equilibrate. Thus, it is presumed that the EAST data did not obtain an equilibrium condition. The plasma torch data, on the other hand, is very consistent with equilibrium assumptions. Thus, the equilibrium data obtained in the plasma torch is utilized to construct the CO 4th Positive emission model, which is then applied to the EAST data with a non-equilibrium interpretation.

A. M. Brandis↗

DOE OSTI2020

Phase Diagram and Structure Map of Binary Nanoparticle Superlattices from a Lennard-Jones Model

A first-principles prediction of the binary nanoparticle phase diagram assembled by solvent evaporation has eluded theoretical approaches. In this paper, we show that a binary system interacting through the Lennard-Jones (LJ) potential contains all experimental phases in which nanoparticles are effectively described as quasi hard spheres. We report a phase diagram consisting of 53 equilibrium phases, whose stability is quite insensitive to the microscopic details of the potentials, thus giving rise to some type of universality. Furthermore, we show that binary lattices may be understood as consisting of certain particle clusters, i.e., motifs, that provide a generalization of the four conventional Frank–Kasper polyhedral units. Our results show that metastable phases share the very same motifs as equilibrium phases. Finally, we discuss the connection with packing models, phase diagrams with repulsive potentials, and the prediction of likely experimental superlattices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Prediction of alkaline earth metal ion adsorption on goethite for various background electrolytes with the CD-MUSIC model

As water scarcity drives the use of more saline water sources, contaminant fate and transport models must capture the impact of high concentrations of alkaline earth metal ions (AEMs) and background electrolytes in these more complex waters. By utilizing macroscopic adsorption data from various electrolyte systems, a Charge Distribution – Multisite Complexation (CD-MUSIC) model, capable of incorporating electrolyte adsorption, was able to accurately simulate the adsorption behavior of alkaline earth metal ions onto goethite. The modeling effort was guided by previous spectroscopic and surface complexation modeling of alkaline earth metal adsorption and built on previous CD-MUSIC modeling that accounted for changes in crystal face contributions to the surface site density as a function of specific surface area. Further, the model was constrained to consider only two dominant surface complex species for each metal ion adsorption reaction. These two species were selected from 44 possible species through objective curve fitting of single-solute macroscopic adsorption data. While most of the alkaline earth metal surface complexes formed outer-sphere complexes at the goethite surface, an inner-sphere species was utilized for Mg 2+ . With the surface complex species and equilibrium constants obtained from this study, the calibrated model successfully predicted alkaline earth metal ion adsorption over a wide range of solution and surface conditions; the model predictions encompassed a wide range of pH (5–11), solute/solid ratio (1.37 × 10 -5 – 8.33 × 10 -4 mol -solute /g -solid ), ionic strengths (0.01 M – 0.7 M), and background electrolytes (Na + , Cs + , Rb + , Cl - , and NO 3 - ) using the same crystal face contribution methodology for site density, capacitance values, and surface acidity constants adopted for proton and cadmium adsorption in previous work (Han and Katz, 2019). Model simulations for a range of background water chemistries demonstrated the potential for Mg 2+ to reduce Cd 2+ adsorption to goethite in model seawater and oil- and gas-produced waters.

42 ENGINEERING↗