Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model physics”

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 1,045 records · Page 58

Distribution System Segmented Model Simplification with Independent Dynamically Changing End-Use Loads

With more emerging technologies in power distribution systems, to support quasi-static time-series (QSTS) simulation over short time frames, accurately modeling the dynamic behavior of end-use load under various conditions is im- portant. However, time and cost can be a burden for the QSTS simulation of thousands of houses representing dynamic load behaviors on large distribution feeders. There are two main contributions of this paper in order to reduce the QSTS simulation complexity. First, the paper presents a method to stochasti- cally model the dynamic behavior of end-use load for QSTS simulations. Second, using the stochastically modeled load behavior as an input to a distribution sys- tem, and building on our previous QSTS simulation simplification/acceleration method that relies on uniform load behavior, this paper presents a segment substitution approach for the simplification of distribution system models with independent dynamically changing end-use loads. Our numerical results show that compared to GridLAB-D, a software tool that simulates end-use load be- havior through physical models, the proposed stochastic modeling approach achieved an 89.3% simulation time reduction for end-use loads, while accurately capturing the temporally changing total load behavior. Moreover, the QSTS simulation of a large test distribution system through OpenDSS showed that the proposed simplification method decreased the feeder simulation time by a factor of three at less than 0.3% voltage error. Compared to our previous method with independent dynamic loads, the voltage error is reduced by a factor of four.

Power distribution system modeling, stochastic pro↗

Raindrop oscillations

A model of the change in shape of a raindrop is presented. Raindrops measured by two orthogonal cameras were classified by shape and orientation to determine the nature of the oscillation. A physical model based on potential energy was then developed to study the amplitude variation of oscillating drops. The model results show that oscillations occur about the equilibrium axis ratio, but the time average axis ratio if significantly more spherical for large amplitudes because of asymmetry in the surface potential energy. A generalization of the model to oscillations produced by turbulence yields average axis ratios that are consistent with the camera measurements. The model results for average axis ratios were applied to rainfall studies with a dual polarized radar.

Beard, K. V.↗

Impact experiments in viscous fluid media

Available phase and group velocity data are inverted by a procedure which includes the effects of transverse anisotropy, anelastic dispersion, sphericity, and gravity. The resulting models, for average Earth, average ocean, and oceanic regions divided according to the age of the ocean floor, are quite different from previous results which ignore the above effects. The models show a low-veocity zone with age dependent anisotropy and velocities higher than derived in previous surface wave studied. The correspondence between the anisotropy variation with age and a physical model based on flow aligned olivine is suggestive. For most of the Earth SHSV in the vicinity of the low-velocity zone. Near the East Pacific Rise, however, SVSH at depth, consistent with ascending flow. Anisotropy is as important as temperature in causing radial and lateral variations in velocity. The models have a high velocity nearly isotropic layer at the top of the mantle that thickens with age. This layer defines the LID, or seismic lithosphere. In the Pacific, the LID thickens with age to a maximum thickness of about 50 km. This thickness is comparable to the thickness of the elastic lithosphere. The LID thickness is thinner than derived using isotropic or pseudo-isotropic procedures A new model for Average Earth is obtained which includes a thin LID. This model extends the fit of a P.R.E.M. type model to shorter period surface waves.

Greeley, R.↗

Image-Based Empirical Modeling of the Plasmasphere

A new suite of empirical models of plasmaspheric plasma based on remote, global images from the IMAGE EUV instrument is proposed for development. The purpose of these empirical models is to establish the statistical properties of the plasmasphere as a function of conditions. This suite of models will mark the first time the plasmaspheric plume is included in an empirical model. Development of these empirical plasmaspheric models will support synoptic studies (such as for wave propagation and growth, energetic particle loss through collisions and dust transport as influenced by charging) and serves as a benchmark against which physical models can be tested. The ability to know that a specific global density distribution occurs in response to specific magnetospheric and solar wind factors is a huge advantage over all previous in-situ based empirical models. The consequence of creating these new plasmaspheric models will be to provide much higher fidelity and much richer quantitative descriptions of the statistical properties of plasmaspheric plasma in the inner magnetosphere, whether that plasma is in the main body of the plasmasphere, nearby during recovery or in the plasmaspheric plume. Model products to be presented include statistical probabilities for being in the plasmasphere, near thermal He+ density boundaries and the complexity of its spatial structure.

Adrian, Mark L.↗

MULTIPHYSICS-MODELING OF FIRE-INDUCED URANIUM AEROSOL FORMATION – A-POSTERIORI BENCHMARKING OF EXPERIMENTS

The formulation and solution of a model that properly described the temperature profiles of uranium in fire conditions was prepared, verified and validated, to support a-posteriori benchmarking of historical experiments. Data from the multi-physics model combined with data regression from the experiments provides a useful tool for exploring sources of potential bias in historical experiments. The approaches used are augmented by visual observations and photographic evidence allowing for correlations between various phenomena, independent of whether or not that specific phenomena are part of the validated model. In this manner, the validated model is a tool for a-posteriori benchmarking, and not a model designed to replicate the experimental output exactly for any specific experiment. For experiments like the Clark (2015) experiments, this integrated model (including chemistry) demonstrated that the largest bias was associated in the metallurgy of the coupons. The model developed does serve as a basis for defending which of the historical alloys tested are the most appropriate for any other alloy being considered (e.g., one that has not yet been tested). The model identifies key reactions that impact uranium chemistry and can significantly bias experimental results. Through use of the model and energy balances, data from historical experiments by Elder and Tinkle (1980) were regressed and benchmarked. This study identifies that all uranium alloys are not equal and that there are distinct differences between the dominant alloys representing distinct metallurgical phases.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling and verification of dynamic field ionization for laser-target interactions

Integrating field ionization models into kinetic plasma simulations is required for a variety of applications, especially when field strengths vary from low to high regimes, such as in laser-target interactions. The introduction of new physics models into kinetic codes requires a rigorous verification of their accuracy through well-defined verification problems. In this work, the field ionization model that has been included in the kinetic plasma code VPIC is presented, along with the detailed approach adopted for its integration. This model includes a comprehensive range of field ionization mechanisms: multiphoton ionization, tunneling ionization, and barrier suppression ionization. New verification problems employed to evaluate the ionization model's fidelity are outlined, and the simulation parameters that affect the accuracy of simulation results are explored. Additionally, this work addresses the impact of field ionization on computational performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coupled model for liquid lithium plasma facing components

Numerical analysis provides the design choice and operating window of liquid metal Plasma Facing Components (PFC) concepts. Coupled analysis of boundary plasma together with the surrounding boundary structures is required. Here, to achieve this goal, PPPL is developing a comprehensive multi-physics model for modeling of PFCs in fusion devices. The model includes the fluid-kinetic code SOLPS-ITER and the flow and heat transfer code CFX from ANSYS. SOLPS-ITER was augmented with a liquid metal boundary condition algorithm, allowing direct two-way coupling of the plasma analysis with the two-dimensional analytical slab flow model which includes heat convection in the liquid metal PFC. The target heat flux resulting from this coupled analysis is used as a boundary condition for detailed 3D Computational Fluid Dynamics (CFD) Magneto Hydro Dynamics (MHD) and heat transfer analysis. A new formulation of MHD equations is introduced in the numerical procedure ensuring current conservation of the discretized equations. Results of the 3D analysis are used for final validation of the coupled model. A PFC design where a porous wall is used to stabilize the liquid metal surface, while MHD drive is used to push the liquid metal flow inside the PFC, will be investigated in the regimes where vapor shielding is created for enhanced volumetric plasma heat dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SAM: A Modern System Code for Advanced Non-LWR Safety Analysis

The System Analysis Module (SAM), developed at Argonne National Laboratory and by collaborators at other organizations, is for advanced non–light water reactor safety analysis. SAM aims to provide fast-running, modest-fidelity, whole-plant transient analysis capabilities that are essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. To facilitate code development, SAM utilizes the MOOSE object-oriented application framework, its underlying finite element library, and linear and nonlinear solvers to leverage modern advanced software environments and numerical methods. SAM aims to solve tightly coupled physical phenomena, including fission reaction, heat transfer, fluid dynamics, and thermal-mechanical responses in advanced reactor structures, systems, and components with high accuracy and efficiency. Finally, this paper gives an overview of the SAM code development, including goals and functional requirements, physical models, current capabilities, verification and validation, software quality assurance, and examples of simulations for advanced nuclear reactor applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fission Evaluation Tools and Analytics (FETA)

This living document presents the Python package FETA. FETA computes observables resulting from the fission process. This document provides the definition of these observables as well as the physics models that are implemented to compute them. Some of these models are used to determine the initial conditions of fission fragments, e.g., the excitation energy E* and spin distribution p(J, π) at scission for prompt decay, while others are related to the nuclear structure and decay properties of the fragments, e.g. the ground-state properties, level density and low-lying excitation spectrum, γ-strength functions and electromagnetic transitions, and neutron transmission coefficients. The end goal for FETA is to enable users to substitute every one of these models by their own files providing these quantities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Searches for Lepton Flavor Violation in Tau Decays at Belle II

Searches for lepton flavor violation in tau decays are unambiguous signatures of new physics. The branching ratios of tau leptons at the level of 10−10–10−9 can be probed using 50 ab−1 of electron-positron annihilation data being collected by the Belle II experiment at the world’s highest luminosity accelerator, the SuperKEKB, located at the High Energy Accelerator Research Organization, KEK, in Tsukuba, Japan. Searches with such expected sensitivity will either discover new physics or strongly constrain several new physics models.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling the Plasmasphere

The plasmasphere, with its accumulation of H(+), He(+), and O(+), is both a strong influence on inner magnetospheric processes and a direct indicator or large and small scale convection electric fields which develop in this region. We will present our empirical and physical modeling of these plasmas and the data that drives them.

Gallagher, Dennis L.↗

Physics-guided machine learning approaches to predict the ideal stability properties of fusion plasmas

One of the biggest challenges to achieve the goal of producing fusion energy in tokamak devices is the necessity of avoiding disruptions of the plasma current due to instabilities. The disruption event characterization and forecasting (DECAF) framework has been developed in this purpose, integrating physics models of many causal events that can lead to a disruption. Two different machine learning approaches are proposed to improve the ideal magnetohydrodynamic (MHD) no-wall limit component of the kinetic stability model included in DECAF. First, a random forest regressor (RFR), was adopted to reproduce the DCON computed change in plasma potential energy without wall effects for a large database of equilibria from the national spherical torus experiment (NSTX). This tree-based method provides an analysis of the importance of each input feature, giving an insight into the underlying physics phenomena. Secondly, a fully-connected neural network has been trained on sets of calculations with the DCON code, to get an improved closed form equation of the no-wall β limit as a function of the relevant plasma parameters indicated by the RFR. The neural network has been guided by physics theory of ideal MHD in its extension outside the domain of the NSTX experimental data. The estimated value has been incorporated into the DECAF kinetic stability model and tested against a set of experimentally stable and unstable discharges. Moreover, the neural network results were used to simulate a real-time stability assessment using only quantities available in real-time. Finally, the portability of the model was investigated, showing encouraging results by testing the NSTX-trained algorithm on the mega ampere spherical tokamak (MAST).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Highly Ordered Hierarchical Anodes for Extreme Fast Charging Batteries (Final Report)

The goal of this project is to enable extreme fast charging (XFC) of Li-ion batteries, which was accomplished through a combination of 1) rational design and manufacturing of hierarchically structured anode architectures; 2) blending graphite/hard carbon into a bulk hybrid anode; 3) engineering artificial solid-electrolyte interphase (SEI) coatings with reduced interphase impedance; 4) computational modeling of coupled transport, kinetic, and electrochemical phenomena; and 5) improved fundamental understanding of lithium plating through operando analysis. This work integrated structural, compositional, and surface modification of graphite anodes, multi-physics modeling of ion transport, electrochemical activity, and heat transfer, advanced strategies to detect Li plating, and semi-automated roll-to-roll cell assembly. The unique facilities at the University of Michigan (UM) and Sandia National Laboratories (SNL) were leveraged to manufacture, prototype, and characterize commercially relevant >2Ah and >180 Wh/kg cells with a target of <20% capacity fade over 500 XFC cycles.

25 ENERGY STORAGE↗

Using ARM Observations to Evaluate Process-Interactions in MCS Simulations Across Scales (Final Progress Report)

This project, funded by DOE Atmospheric System Research (DE-SC0020050), focused on improving the representation of mesoscale convective systems (MCSs) in numerical weather and climate models by leveraging high-resolution observations from the DOE Atmospheric Radiation Measurement (ARM) program. The research aimed to evaluate model sensitivities to grid spacing, microphysics, and planetary boundary layer (PBL) schemes, with a particular emphasis on improving convection parameterization for high-resolution modeling. Findings from this work highlight several key advancements. Model validation against ARM radar wind profiler data from the Southern Great Plains (SGP) and Manaus (MAO) sites revealed systematic biases in simulated convective mass flux profiles, leading to the development of an observationally constrained evaluation framework for diagnosing and improving model performance. Sensitivity analyses demonstrated that the representation of Amazonian MCSs was highly dependent on PBL scheme selection, while mid-latitude MCSs were more strongly influenced by microphysics parameterizations. A series of high-resolution WRF simulations, ranging from 4 km to 125 m grid spacing, provided insight into the behavior of convective drafts across scales. While updraft properties converged at sub-kilometer resolutions, biases in downdraft intensity persisted even at the finest resolution tested, emphasizing the need for further refinements in model physics. Additionally, comparisons of MCS vertical structures between mid-latitude and tropical environments revealed stronger updrafts and larger mass flux in mid-latitude MCSs, providing critical insights for improving climate model representations of storm-scale dynamics. The project’s findings have already contributed to advancing numerical modeling capabilities, particularly in WRF, MPAS, ICON, and DOE’s SCREAM model, by refining how convective processes are represented in high-resolution climate simulations. Results were disseminated through peer-reviewed publications, conference presentations, and ARM/ASR Research Highlights, engaging the broader scientific community. The project also provided valuable training opportunities for two postdoctoral researchers, who played central roles in model development, analysis, and dissemination of results. Their work contributed to several publications and conference presentations, helping prepare them for careers in atmospheric modeling. By improving the simulation of MCSs, this research directly supports the development of next-generation climate models capable of more accurately representing extreme precipitation and convective processes. The insights gained will inform future improvements in convective parameterization and guide the design of high-resolution weather and climate simulations, ultimately enhancing the reliability of climate projections and weather forecasts.

54 ENVIRONMENTAL SCIENCES↗

The effect of solute additions on the steady-state creep behavior of dispersion-strengthened aluminum.

The effect of solute additions on the steady-state creep behavior of coarse-grained dispersion-strengthened aluminum alloys was studied. Recrystallized dispersion-strengthened solid solutions were found to have stress and temperature sensitivities quite unlike those observed in single-phase solid solutions having the same composition and grain size. The addition of magnesium or copper to the matrix of a recrystallized dispersion-strengthened aluminum causes a decrease in the steady-state creep rate which is much smaller than that caused by similar amounts of solute in single-phase solid solutions. All alloys exhibited essentially a 4.0 power stress exponent in agreement with the model of Ansell and Weertman. The activation energy for steady-state creep in dispersion-strengthened Al-Mg alloys, as well as the stress dependence, was in agreement with the physical model of dislocation climb over the dispersed particles.

Reynolds, G. H.↗

Monte Carlo investigation of transient acoustic fields in partially or completely bounded medium

A simple repetitive calculation was used to investigate what happens to the field in terms of the signal paths of disturbances originating from the energy source. The computation allowed the field to be reconstructed as a function of space and time on a statistical basis. The suggested Monte Carlo method is in response to the need for a numerical method to supplement analytical methods of solution which are only valid when the boundaries have simple shapes, rather than for a medium that is bounded. For the analysis, a suitable model was created from which was developed an algorithm for the estimation of acoustic pressure variations in the region under investigation. The validity of the technique was demonstrated by analysis of simple physical models with the aid of a digital computer. The Monte Carlo method is applicable to a medium which is homogeneous and is enclosed by either rectangular or curved boundaries.

Thanedar, B. D.↗

Biological indicators for monitoring water quality of MTF canals system

Biological models, diversity indexes, were developed to predict environmental effects of NASA's Mississippi test facility (MTF) chemical operations on canal systems in the area. To predict the effects on local streams, a physical model of unpolluted streams was established. The model is fed by artesian well water free of background levels of pollutants. The species diversity and biota composition of unpolluted MTF stream was determined; resulting information will be used to form baseline data for future comparisons. Biological modeling was accomplished by adding controlled quantities or kinds of chemical pollutants and evaluating the effects of these chemicals on the biological life of the stream.

Sethi, S. L.↗

On the cosmological evolution of the X-ray emission from quasars

The average dependence of the ratio of X-ray luminosity to optical luminosity as a function of redshift and optical luminosity for quasars is derived, and it is found that the explicit dependence of this ratio is predominantly on optical luminosity. For a wide class of models for the cosmological evolution of quasars, the results imply that the evolution of the X-ray luminosity function is weaker than the evolution of the optical luminosity function. Implications for physical models of quasars and for global properties of the quasar population are indicated.

Avni, Y.↗