Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physics Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,495 records · Page 83

An Integrated Reliability and Physics-Based Risk Modeling Approach for Assessing Human Spaceflight Systems

This paper presents an integrated reliability and physics-based risk modeling approach for assessing human spaceflight systems. The approach is demonstrated using an example, end-to-end risk assessment of a generic-crewed space transportation system during a reference mission to the International Space Station. The behavior of the system is modeled using analysis techniques from multiple disciplines in order to properly capture the dynamic time- and state- dependent consequences of failures encountered in different mission phases. We discuss how to combine traditional reliability analyses with Monte Carlo simulation methods and physics-based engineering models to produce loss-of- mission and loss-of-crew risk estimates supporting risk-based decision-making and requirement verification. This approach facilitates risk-informed design by providing more realistic representation of system failures and interactions; identifying key risk-driving sensitivities, dependencies, and assumptions; and tracking multiple figures of merit within a single, responsive assessment framework that can readily incorporate evolving design information throughout system development.

Risk assessment↗

EVIDENCE FOR LOW BLACK HOLE SPIN AND PHYSICALLY MOTIVATED ACCRETION MODELS FROM MILLIMETER-VLBI OBSERVATIONS OF SAGITTARIUS A*

Millimeter very long baseline interferometry (mm-VLBI) provides the novel capacity to probe the emission region of a handful of supermassive black holes on sub-horizon scales. For Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, this provides access to the region in the immediate vicinity of the horizon. Broderick et al. have already shown that by leveraging spectral and polarization information as well as accretion theory, it is possible to extract accretion-model parameters (including black hole spin) from mm-VLBI experiments containing only a handful of telescopes. Here we repeat this analysis with the most recent mm-VLBI data, considering a class of aligned, radiatively inefficient accretion flow (RIAF) models. We find that the combined data set rules out symmetric models for Sgr A*'s flux distribution at the 3.9{sigma} level, strongly favoring length-to-width ratios of roughly 2.4:1. More importantly, we find that physically motivated accretion flow models provide a significantly better fit to the mm-VLBI observations than phenomenological models, at the 2.9{sigma} level. This implies that not only is mm-VLBI presently capable of distinguishing between potential physical models for Sgr A*'s emission, but further that it is sensitive to the strong gravitational lensing associated with the propagation of photons near the black hole. Based upon this analysis we find that the most probable magnitude, viewing angle, and position angle for the black hole spin are a = 0.0{sup +0.64+0.86}, {theta}=68{sup o+5o+9o}{sub -20}{sup o}{sub -28}{sup o}, and {xi}=-52{sup o+17o+33o}{sub -15}{sup o}{sub -24}{sup o} east of north, where the errors quoted are the 1{sigma} and 2{sigma} uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Detection of thermal SZ-CMB lensing cross-correlation in Planck nominal mission data

The nominal mission maps from the Planck satellite contain a wealth of information about secondary anisotropies in the cosmic microwave background (CMB), including those induced by the thermal Sunyaev-Zel'dovich (tSZ) effect and gravitational lensing. As both the tSZ and CMB lensing signals trace the large-scale matter density field, the anisotropies sourced by these processes are expected to be correlated. We report the first detection of this cross-correlation signal, which we measure at 6.2σ significance using the Planck data. We take advantage of Planck's multifrequency coverage to construct a tSZ map using internal linear combination techniques, which we subsequently cross-correlate with the publicly-released Planck CMB lensing potential map. The cross-correlation is subject to contamination from the cosmic infrared background (CIB), which is known to correlate strongly with CMB lensing. We correct for this contamination via cross-correlating our tSZ map with the Planck 857 GHz map and confirm the robustness of our measurement using several null tests. We interpret the signal using halo model calculations, which indicate that the tSZ-CMB lensing cross-correlation is a unique probe of the physics of intracluster gas in high-redshift, low-mass groups and clusters. Our results are consistent with extrapolations of existing gas physics models to this previously unexplored regime and show clear evidence for contributions from both the one- and two-halo terms, but no statistically significant evidence for contributions from diffuse, unbound gas outside of collapsed halos. We also show that the amplitude of the signal depends rather sensitively on the amplitude of fluctuations (σ{sub 8}) and the matter density (Ω{sub m}), scaling as σ{sub 8}{sup 6.1}Ω{sub m}{sup 1.5} at ℓ = 1000. We constrain the degenerate combination σ{sub 8}(Ω{sub m}/0.282){sup 0.26} = 0.824±0.029, a result that is in less tension with primordial CMB constraints than some recent tSZ analyses. We also combine our measurement with the Planck measurement of the tSZ auto-power spectrum to demonstrate a technique that can in principle constrain both cosmology and the physics of intracluster gas simultaneously. Our detection is a direct confirmation that hot, ionized gas traces the dark matter distribution over a wide range of scales in the universe ( ∼ 0.1–50 Mpc/h)

79 ASTRONOMY AND ASTROPHYSICS↗

Tunable thermoelectric transport in nanomeshes via elastic strain engineering

Recent experimental explorations of silicon nanomeshes have shown that the unique metastructures exhibit reduced thermal conductivity while preserving bulk electrical conductivity via feature sizes between relevant phonon and electron mean free paths, aiding in the continued promise that nanometer-scale engineering may further enhance thermoelectric behavior. Here, we introduce a strategy for tuning thermoelectric transport phenomena in semiconductor nanomeshes via heterogeneous elastic strain engineering, using silicon as a model material for demonstration of the concept. By combining analytical models for electron mobility in uniformly stressed silicon with finite element analysis of strained silicon nanomeshes in a lumped physical model, we show that the nonuniform and multiaxial strain fields defined by the nanomesh geometry give rise to spatially varying band shifts and warping, which in aggregate accelerate electron transport along directions of applied stress. This allows for global electrical conductivity and Seebeck enhancements beyond those of homogenous samples under equivalent far-field stresses, ultimately increasing thermoelectric power factor nearly 50% over unstrained samples. The proposed concept and structures—generic to a wide class of materials with large dynamic ranges of elastic strain in nanoscale volumes—may enable a new pathway for active and tunable control of transport properties relevant to waste heat scavenging and thermal management.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Development of Physics-Based Transition Models for Unstructured-Mesh CFD Codes Using Deep Learning Models

Predicting transition locations over a vehicle surface is of fundamental importance for many engineering applications. With the transition information, the Reynolds-averaged Navier-Stokes (RANS) computations can turn on the turbulence model at the right locations so that drag, lift and other aerodynamic quantities can be accurately predicted. In contrast to the popularity of RANS-based transition modeling in which transition onset is governed by the turbulence equations, physics-based transition models that account for instability waves within the boundary layer, thus more compliant to flow physics, only gained more attention in recent years. This paper describes the development of a new physics-based transition model based on either the linear stability theory (LST) or parabolized stability equations (PSE). The model is designed to communicate with a structured or unstructured-mesh RANS solver back and forth in order to more accurately compute transition fronts over a three-dimensional body. In the developed model, the Python suite of interface codes in conjunction with the LASTRAC software can be executed autonomously to produce transition onset locations for a given laminar or RANS-computed transitional state. In addition, as a proof of concept, the tool set consists of a deep learning neural network model that has been designed and trained to predict instability wave evolutions inside the boundary layer for various instability wave mechanisms across a selected speed range. A machine-learned intelligent profile interpolation model has also been devised to enable reliable instability-wave spectra predictions with just a few points in the mean flow profiles.

Transition Modeling↗

Simulation of Impedance Changes with Aging in Lithium Titanate-based Cells Using Physics-Based Dimensionless Modeling

Quantifying aging effects in lithium-ion cells with chemistries that have a flat open circuit potential is challenging. We implement a physics-based electrochemical model to track changes in the electrochemical impedance response of lithium titanate-based cells. Frequency domain equations of a pseudo two-dimensional model are made dimensionless, and the corresponding non-dimensional parameters are estimated using a Levenberg-Marquardt routine. The model weighs the relative contributions of changes in diffusion, ionic conduction within the electrolyte phase against solid phase electronic conduction towards cell aging. Solid-phase diffusion, charge transfer resistance and double layer capacitance at the solid-liquid interface are accounted for in the particle impedance. The estimation routine tracks dimensionless parameters using accelerated cycling data from full cells over 1000 cycles. The model can be deployed within a short time for state estimation using physics-based models without requiring prior knowledge of the battery chemistry, format, or capacity.

25 ENERGY STORAGE↗

Modeling Warm Dense Matter Experiments using the 3D ALE-AMR Code and the Move Toward Exascale Computing

The Neutralized Drift Compression Experiment II (NDCX II) is an induction accelerator planned for initial commissioning in 2012. The final design calls for a 3 MeV, Li+ ion beam, delivered in a bunch with characteristic pulse duration of 1 ns, and transverse dimension of order 1 mm. The NDCX II will be used in studies of material in the warm dense matter (WDM) regime, and ion beam/hydrodynamic coupling experiments relevant to heavy ion based inertial fusion energy. We discuss recent efforts to adapt the 3D ALE-AMR code to model WDM experiments on NDCX II. The code, which combines Arbitrary Lagrangian Eulerian (ALE) hydrodynamics with Adaptive Mesh Refinement (AMR), has physics models that include ion deposition, radiation hydrodynamics, thermal diffusion, anisotropic material strength with material time history, and advanced models for fragmentation. Experiments at NDCX-II will explore the process of bubble and droplet formation (two-phase expansion) of superheated metal solids using ion beams. Experiments at higher temperatures will explore equation of state and heavy ion fusion beam-to-target energy coupling efficiency. Ion beams allow precise control of local beam energy deposition providing uniform volumetric heating on a timescale shorter than that of hydrodynamic expansion. The ALE-AMR code does not have any export control restrictions and is currently running at the National Energy Research Scientific Computing Center (NERSC) at LBNL and has been shown to scale well to thousands of CPUs. New surface tension models that are being implemented and applied to WDM experiments. Some of the approaches use a diffuse interface surface tension model that is based on the advective Cahn-Hilliard equations, which allows for droplet breakup in divergent velocity fields without the need for imposed perturbations. Other methods require seeding or other methods for droplet breakup. We also briefly discuss the effects of the move to exascale computing and related computational changes on general modeling codes in fusion energy.

42 ENGINEERING↗

Recent Improvements to the Calibration Models for RXTE/PCA

We are updating the calibration of the PCA to correct for slow variations, primarily in energy to channel relationship. We have also improved the physical model in the vicinity of the Xe K-edge, which should increase the reliability of continuum fits above 20 keV. The improvements to the matrix are especially important to simultaneous observations, where the PCA is often used to constrain the continuum while other higher resolution spectrometers are used to study the shape of lines and edges associated with Iron.

Jahoda, K.↗

A composite position independent monitor of reactor fuel irradiation using Pu, Cs, and Ba isotope ratios

When post-irradiation materials from the nuclear fuel cycle are released to the environment, certain isotopes of actinides and fission products carry signatures of irradiation history that can potentially aid a nuclear forensic investigation into the material's provenance. In this study, combinations of Pu, Cs, and Ba isotope ratios that produce position (in the reactor core) independent monitors of irradiation history in spent light water reactor fuel are identified and explored. These position independent monitors (PIMs) are modeled for various irradiation scenarios using automated depletion codes as well as ordinary differential equation solutions to approximate nuclear physics models. Experimental validation was performed using irradiated low enriched uranium oxide fuel from a light water reactor, which was sampled at 8 axial positions from a single rod. Plutonium, barium and cesium were chemically separated and isotope ratio measurements of the separated solutions were made by quadrupole and multi-collector inductively coupled mass spectrometry (Cs and Pu, respectively) and thermal ionization mass spectrometry (Ba). The effect of axial variations in neutron fluence and energy spectrum are evident in the measured isotope ratios. Furthermore, two versions of a combined Pu and Cs based PIM are developed. A linear PIM model, which can be used to solve for irradiation time is found to work well for natural U fuel with <10% 240 Pu and known or short cooling times. A non-linear PIM model, which cannot be solved explicitly for irradiation time without additional information, can nonetheless still group samples by irradiation history, including high burnup LEU fuel with unknown cooling time. 137 Ba/ 138 Ba is also observed to act as a position independent monitor; it is nearly single valued across the sampled fuel rod, indicating that samples sharing an irradiation history (same irradiation time and cooling time) in a reactor despite experiencing different neutron fluxes will have a common 137 Ba/ 138 Ba ratio. Modeling of this Ba PIM shows it increases monotonically with irradiation and cooling time, and a confirmatory first order analytical solution is also presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Ensemble Spread Behavior in Coupled Climate Models: Insights From the Energy Exascale Earth System Model Version 1 Large Ensemble

AbstractAssessing uncertainty in future climate projections requires understanding both internal climate variability and external forcing. For this reason, single‐model initial condition large ensembles (SMILEs) run with Earth System Models (ESMs) have recently become popular. Here we present a new 20‐member SMILE with the Energy Exascale Earth System Model version 1 (E3SMv1‐LE), which uses a “macro” initialization strategy choosing coupled atmosphere/ocean states based on inter‐basin contrasts in ocean heat content (OHC). The E3SMv1‐LE simulates tropical climate variability well, albeit with a muted warming trend over the twentieth century due to overly strong aerosol forcing. The E3SMv1‐LE's initial climate spread is comparable to other (larger) SMILEs, suggesting that maximizing inter‐basin ocean heat contrasts may be an efficient method of generating ensemble spread. We also compare different ensemble spread across multiple SMILEs, using surface air temperature and OHC. The Community Earth system Model version 1, the only ensemble which utilizes a “micro” initialization approach perturbing only atmospheric initial conditions, yields lower spread in the first ∼30 years. The E3SMv1‐LE exhibits a relatively large spread, with some evidence for anthropogenic forcing influencing spread in the late twentieth century. However, systematic effects of differing “macro” initialization strategies are difficult to detect, possibly resulting from differing model physics or responses to external forcing. Notably, the method of standardizing results affects ensemble spread: control simulations for most models have either large background trends or multi‐centennial variability in OHC. This spurious disequlibrium behavior is a substantial roadblock to understanding both internal climate variability and its response to forcing.

Stevenson, Samantha↗

New Methods of Uncertainty Quantification for Mixed Discrete-Continuous Variable Models

The scale and complexity of problems such as designing power grids or planning for climate change is growing rapidly, driving the development of complicated computer models. More complex models have longer run times and incorporate larger numbers of inputs, both continuous and discrete. For example, a detailed physics model may have continuous variables such as temperature, height or pressure along with discrete variables that indicate the choice of a material for a particular piece or the model to be used to calculate air flow. A power grid design model may have continuous variables such as generation capacity, power flow or demand along with discrete variables such as number of generators, number of transmission lines or binary variables to indicate whether or not a node is chosen for generation expansion. A growing awareness of uncertainty and the desire to make risk-informed decisions is causing uncertainty quantification (UQ) to be more routine and often required. UQ provides the underpinnings necessary to establish confidence in models and their use; therefore, much time and effort is being invested in creating efficient approaches for UQ. However, these efforts have been focused on models that take continuous variables as inputs. When discrete inputs are thrown into the mix, the basic approach is to repeat the UQ analysis for each combination of discrete inputs or some subset thereof; this rapidly becomes intractable. Because of the computational complexity inherent in mixed discrete-continuous models, researchers will focus on the uncertainty in their particular problem finding ways to take advantage of symmetries, simplifications or structures. For example, uncertainty propagation in certain dynamical systems can be efficiently carried out after various decomposition steps or uncertainty propagation in stochastic programming is confined to scenario generation. Unfortunately models are not always available for such machinations: models may be embedded in legacy codes, may utilize commercial off the shelf codes or may be created by stringing a series of codes together. It is also time consuming to start each problem from scratch; worse there may not be any simplifications or symmetries to take advantage of. For these situations a UQ method developed for any black box function is necessary. This report documents a new conceptual model for performing UQ for mixed discrete-continuous models which not only applies to any simulator function, but allows the use of the efficient UQ methods that have been developed for continuous inputs only. The conceptual model is presented and an estimation procedure is fleshed out for one class of problems. This is applied to variations of a mixed discrete-continuous optimization test problem. This procedure provides comparable results to a benchmark solution with fewer function evaluations.

Bauman, Lara E.↗

Development and Validation of a Two-Phase Thermal-Hydraulic CFD Code NEK-2P

A project is underway to develop, verify and validate an advanced two-phase flow modeling capability for the highly-scalable, high-performance Computational Fluid Dynamics (CFD) code NEK5000. The goal of this work is to verify and validate the two-phase version of the NEK5000 code, named NEK-2P, to simulate the two-phase flow and heat transfer phenomena that occur in a Boiling Water Reactor (BWR) fuel bundle under various operating conditions. The NEK-2P two-phase flow models follow the approach used for the Extended Boiling Framework (EBF) previously developed at Argonne but include more fundamental physical models of boiling phenomena and advanced numerical algorithms for improved computational accuracy, robustness, and computational speed. The development of the NEK-2P two-phase solver and the implementation of the Extended Boiling Framework two-phase models were initially supported by Argonne National Laboratory (Argonne) through a Laboratory Directed Research and Development (LDRD) project during FY2014-2016. The development and validation of the two-phase models through analyses of selected two-phase boiling flow experiments was supported by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY2017-2020. This report focuses on verification and validation of the water-steam boiling model NEK-2P Two-Phase, CFD code. The NEK-2P was validated with Nuclear Power Engineering Corporation (NUPEC) Pressurized Water Reactor (PWR) Sub-channel and Bundle Test (PSBT) void distribution benchmark. Three different simulations were performed and analyzed for various operating conditions such as wall-heat flux and sub-cooled inlet temperatures. Reasonably good agreement with measured data was obtained in predicting the measured void distributions. Simulations were performed for Virginia Tech. (VT) 3x3 rod bundle geometry with and without spacers. The preliminary results were presented for Simplified Spacer Grid (SSG). In addition, the implementation of interface reconstruction model was tested with one of the Becker benchmark Critical Heat Flux (CHF) experiments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A study of liquid propellant autoignition

Data and theory pertinent to the autoignition of liquid oxygen/liquid hydrogen and liquid oxygen propellants were reviewed. Physical models of the processes supporting or contributing to autoignition were developed. Emphasis was placed on the description of the physical environment and its relationship to the autoignition phenomenon.

Lester, D. H.↗

Quantitative NDE applied to composites and metals

Research at the NASA/Langley Research Center concerning quantitative NDE of composites and metals is reviewed. The relationship between ultrasonics and polymer cure is outlined. NDE models are presented, which can be used to develop measurement technologies for characterizing the curing of a polymer system for composite materials. The models can be used to determine the glass transition temperature, the degree of cure, and the cure rate. The application of the model to control autoclave processing of composite materials is noted. Consideration is given to the use of thermal diffusion models combined with controlled thermal input measurements to determine the thermal diffusivity of materials. Also, a two-dimensional physical model is described that permits delaminations in samples of Space Shuttle Solid Rocket Motors to be detected in thermograms in the presence of cooling effects and uneven heating.

Heyman, Joseph S.↗

Modeling Carbon and Hydrocarbon Molecular Structures in EZTB

A software module that models the electronic and mechanical aspects of hydrocarbon molecules and carbon molecular structures on the basis of first principles has been written for incorporation into, and execution within, the Easy (Modular) Tight-Binding (EZTB) software infrastructure, which is summarized briefly in the immediately preceding article. Of particular interest, this module can model carbon crystals and nanotubes characterized by various coordinates and containing defects, without need to adjust parameters of the physical model. The module has been used to study the changes in electronic properties of carbon nanotubes, caused by bending of the nanotubes, for potential utility as the basis of a nonvolatile, electriccharge- free memory devices. For example, in one application of the module, it was found that an initially 50-nmlong carbon, (10,10)-chirality nanotube, which is a metallic conductor when straight, becomes a semiconductor with an energy gap of .3 meV when bent to a lateral displacement of 4 nm at the middle.

Lee, Seungwon↗

Flavor-violating Higgs and Z boson decays at a future circular lepton collider

Recent advances in b , c , and s quark tagging coupled with novel statistical analysis techniques will allow future high energy and high statistics electron-positron colliders, such as the FCC-ee, to place phenomenologically relevant bounds on flavor violating Higgs and Z decays to quarks. We assess the FCC-ee reach for Z / h → b s , c u decays as a function of jet tagging performance. We also update the standard model (SM) predictions for the corresponding branching ratios, as well as the indirect constraints on the flavor violating Higgs and Z couplings to quarks. Using the type III two Higgs doublet model as an example of beyond the standard model physics, we show that the searches for h → b s , c u decays at FCC-ee can probe new parameter space not excluded by indirect searches. We also reinterpret the FCC-ee reach for Z → b s , c u in terms of the constraints on models with vectorlike quarks. Published by the American Physical Society 2024

Kamenik, Jernej F. (ORCID:0000000151838637)↗

Solving key challenges in collider physics with foundation models

Foundation models are neural networks that are capable of simultaneously solving many problems. Large language foundation models like ChatGPT have revolutionized many aspects of daily life, but their impact for science is not yet clear. In this paper, we use a new foundation model for hadronic jets to solve three key challenges in collider physics. In particular, we show how experiments can (1) save significant computing power when developing reconstruction algorithms, (2) perform a complete uncertainty quantification for high-dimensional measurements, and (3) search for new physics with model agnostic methods using low-level inputs. In each case, there are significant computational or methodological challenges with current methods that limit the science potential of deep learning algorithms. By solving each problem, we take jet foundation models beyond proof-of-principle studies and into the toolkit of practitioners.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

White dwarfs constraints on dark sector models with light particles

The white dwarf luminosity function is well understood in terms of standard model physics and leaves little room for exotic cooling mechanisms related to the possible existence of new weakly interacting light particles. This puts significant constraints on the parameter space of models that contain a massive dark photon and light dark sector particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗