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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Enhanced thermochemical nonequilibrium computations of flow around the aeroassist flight experiment vehicle

A three-dimensional explicit, finite-rate, shock-capturing numerical algorithm is used to calculate thermochemical nonequilibrium flowfields about the Aeroassist Flight Experiment vehicle at one of its flight trajectory points. The full Navier-Stokes equations and an eleven species chemical model with the latest reaction rates are incorporated into the code. Results are compared against experimental data and other numerical solutions. The effects of changes to the physical model on the stagnation line and base region flow and surface quantities on the forebody are investigated.

Palmer, Grant↗
DOE OSTI2020

Interacting radiation after Planck and its implications for the Hubble tension

Standard cosmology predicts that prior to matter-radiation equality about 41% of the energy density was in free-streaming neutrinos. In many beyond Standard Model scenarios, however, the amount and free-streaming nature of this component is modified. For example, this occurs in models with new neutrino self-interactions or an additional dark sector with interacting light particles. Here, we consider several extensions of the standard cosmology that include a non-free-streaming radiation component as motivated by such particle physics models and use the final Planck data release to constrain them. This release contains significant improvements in the polarization likelihood which plays an important role in distinguishing free-streaming from interacting radiation species. Fixing the total amount of energy in radiation to match the expectation from standard neutrino decoupling we find that the fraction of free-streaming radiation must be ffs > 0.8 at 95% CL (combining temperature, polarization and baryon acoustic oscillation data). Allowing for arbitrary contributions of free-streaming and interacting radiation, the effective number of new non-free-streaming degrees of freedom is constrained to be Nfld < 0.6 at 95% CL. Cosmologies with additional radiation are also known to ease the discrepancy between the local measurement and CMB inference of the current expansion rate H0. We show that including a non-free-streaming radiation component allows for a larger amount of total energy density in radiation, leading to a mild improvement of the fit to cosmological data compared to previously discussed models with only a free-streaming component.

particle physics - cosmology connection↗
DOE OSTI2023

Modeling noise in global Mølmer-Sørensen interactions applied to quantum approximate optimization

Many-qubit Mølmer-Sørensen (MS) interactions applied to trapped ions offer unique capabilities for quantum information processing, with applications including quantum simulation and the quantum approximate optimization algorithm (QAOA). Here, we develop a physical model to describe many-qubit MS interactions under four sources of experimental noise: vibrational mode frequency fluctuations, laser power fluctuations, thermal initial vibrational states, and state preparation and measurement errors. The model parametrizes these errors from simple experimental measurements, without free parameters. Herein, we validate the model in comparison with experiments that implement sequences of MS interactions on two 171 Yb + ions. The model shows reasonable agreement after several MS interactions as quantified by the reduced chi-squared statistic $χ^{2}_{red}$≈2. As an application we examine MaxCut QAOA experiments on three and six ions. The experimental performance is quantified by approximation ratios that are 91% and 83% of the optimal theoretical values. Our model predicts $0.93^{+0.03}_{-0.02}$ and $0.95^{+0.04}_{-0.03}$, respectively, with disagreement in the latter value attributable to secondary noise sources beyond those considered in our analysis. With realistic experimental improvements to reduce measurement error and radial trap frequency variations, the model achieves approximation ratios that are 99% of the optimal. Incorporating these improvements into future experiments is expected to reveal new aspects of noise for future modeling and experimental improvements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Modeling Framework to Analyze Performance and Structural Reliability of Solid Oxide Electrolysis Cells

Solid oxide electrolysis cells (SOEC) have been receiving significant attention recently because of their high energy efficiency and fast hydrogen production. In this study a multi-physics model to simulate the SOEC performance and structural reliability of a state-of-the-art planar SOEC design was developed. The electrochemical reactions, fluid dynamics, species transport, electron transfer, and heat transfer were modeled in the commercial computational fluid dynamics (CFD) software STAR-CCM+. The thermomechanical analysis and the associated structural reliability evaluations were conducted using the commercial finite element analysis software ANSYS. The electrochemistry model was validated by using the experimentally obtained current-voltage (I-V) characteristics of the electrode-supported SOECs. The reliability analysis using a risk-of-rupture approach showed low failure probabilities under standard operating conditions considered in this study. For cells operated at voltages well above a thermoneutral voltage, the reliability evaluations indicated a potential risk of cell failure, but the damage was concentrated locally in specific areas of the cell which typically do not lead to total loss of cell function. The presented approach provides insights for evaluating representative cell and stack performances and structural reliability without intensive testing and for developing optimally performing and structurally reliable SOECs for efficient hydrogen generation.

25 ENERGY STORAGE↗
DOE OSTI2021

Comments on 'Do ultrafine cloud condensation nuclei invigorate deep convection?'

Here we elaborated on the deficiencies associated with the theoretical arguments and model simulations in Grawbowski and Morrison (2020; GM20), and showed that the concept of convective invigoration by aerosols can be supported by both accurate theoretical analysis and explicit physics modeling with prognostic aerosols and supersaturation. Above the freezing level, droplet freezing itself does not change mass loading and the latent heat release resulting from freezing leads to an increase in the buoyancy. The derivation in GM20 is erroneous because of omitting the liquid water sink term. They also ignored the effect of enhanced latent heating from the subsequent growth processes (deposition and riming) after freezing. Regarding aerosol effect on condensational growth, the quasi-steady assumption for supersaturation as adopted in GM20 is invalidated especially when droplet number concentration is low and updraft is strong from both theoretical analysis and our bin microphysics modeling results. The quasi-steady assumption makes condensation independent of droplet number and size and therefore reduces the overall aerosol effects on cloud dynamics. Studying aerosol effects on clouds particularly deep convective clouds requires the exact solution of supersaturation. Any assumption imposed to supersaturation such as saturation adjustment and quasi-steady approximation leads to incorrect aerosol effects on diffusional growth. The exact supersaturation equation clearly shows that condensation depends on droplet number and size, and more droplets in the polluted clouds increase condensation (thus lower supersaturation), leading to enhanced buoyancy and updraft speeds.

Convective invigoration, deep convective clouds, a↗
NASA NTRS1993

A physical optics/equivalent currents model for the RCS of trihedral corner reflectors

The scattering in the interior regions of both square and triangular trihedral corner reflectors is examined. The theoretical model presented combines geometrical and physical optics (GO and PO), used to account for reflection terms, with equivalent edge currents (EEC), used to account for first-order diffractions from the edges. First-order, second-order, and third-order reflection terms are included. Calculating the first-order reflection terms involves integrating over the entire surface of the illuminated plate. Calculating the second- and third-order reflection terms, however, is much more difficult because the illuminated area is an arbitrary polygon whose shape is dependent upon the incident angles. The method for determining the area of integration is detailed. Extensive comparisons between the high-frequency model, Finite-Difference Time-Domain (FDTD) and experimental data are used for validation of the radar cross section (RCS) of both square and triangular trihedral reflectors.

Balanis, Constantine A.↗
NASA NTRSDate not supplied

Tropospheric Ozone in the NASA GEOS Model: Effects of Satellite NO2 Data Assimilation and Improvements in Background NOx Chemistry

Tropospheric ozone is a major atmospheric oxidant, the primary source of the hydroxyl radical, and a greenhouse gas. In the troposphere, it is produced from the oxidation of carbon monoxide and volatile organic compounds in the presence of nitrogen oxides (NOx=NO+NO2) or has its origins in the stratosphere. The tropospheric ozone burdens in global atmospheric chemistry models disagree by a factor of 1.5, indicating an incomplete understanding of ozone sources and sinks. Here we examine the tropospheric ozone simulation in NASA Global Modeling and Assimilation Office’s (GMAO’s) GEOS model with the GEOS-Chem chemistry mechanism (version 14.2), using observations from satellites, aircraft, ozonesondes, and surface sites. This modeling framework is used to produce GMAO’s GEOS Composition Forecasts (GEOS-CF), and it has been updated extensively with new emission inventories, improved model physics, satellite data assimilation capability, and an up-to-date chemical mechanism that includes tropospheric halogen chemistry and NOx recycling from particulate nitrate. Earlier versions of the model using GEOS-Chem version 12.0 showed an underestimate in tropospheric ozone in the northern hemisphere. We evaluate the ozone simulation in the updated model, focusing on the effects of assimilating satellite NO2 data from the Ozone Monitoring Instrument (OMI), constraining stratospheric ozone to satellite data, and updates to the NOx chemistry. We will examine changes in the vertical profiles of tropospheric ozone over the US in support of the Tropospheric Emissions: Monitoring of Pollution (TEMPO), and discuss implications for the tropospheric ozone budget.

V. Shah↗
NASA NTRSDate not supplied

Tropospheric Ozone in the NASA GEOS Model: Effects of Satellite Data Assimilation and Improvements in Background NO x Chemistry

Tropospheric ozone is a major atmospheric oxidant, the primary source of the hydroxyl radical, and a greenhouse gas. In the troposphere, it is produced from the oxidation of carbon monoxide and volatile organic compounds in the presence of nitrogen oxides (NO x =NO+NO 2 ) or has its origins in the stratosphere. The tropospheric ozone burdens in global atmospheric chemistry models disagree by a factor of 1.5, indicating an incomplete understanding of ozone sources and sinks. Here we examine the tropospheric ozone simulation in NASA Global Modeling and Assimilation Office’s (GMAO’s) GEOS model with the GEOS-Chem chemistry mechanism (version 14.2), using observations from satellites, aircraft, ozonesondes, and surface sites. This modeling framework is used to produce GMAO’s GEOS Composition Forecasts (GEOS-CF), and it has been updated extensively with new emission inventories, improved model physics, satellite data assimilation capability, and an up-to-date chemical mechanism that includes tropospheric halogen chemistry and NOx recycling from particulate nitrate. Earlier versions of the model using GEOS-Chem version 12.0 showed an underestimate in tropospheric ozone in the northern hemisphere. We evaluate the ozone simulation in the updated model, focusing on the effects of assimilating satellite NO 2 data from the Ozone Monitoring Instrument (OMI), constraining stratospheric ozone to satellite data, and updates to the NOx chemistry. We will examine changes in the vertical profiles of tropospheric ozone over the US in support of the Tropospheric Emissions: Monitoring of Pollution (TEMPO), and discuss implications for the tropospheric ozone budget.

Viral Shah↗
DOE OSTI2024

Comparing pion production in transport simulations of heavy-ion collisions at 270⁢𝐴 MeV under controlled conditions

Within the TMEP, we present a detailed study of the performance of different transport models in Sn+Sn collisions at 270A MeV, and put particular emphasis on the production of pions and Delta resonances, which have been used as probes of the nuclear symmetry energy. We prescribe a common and rather simple physics model, and follow in detail the results of 4 BUU models and 6 QMD models. The nucleonic evolution of the collision and the nucleonic observables in these codes do not completely converge, but the differences among the codes can be understood as being due to several reasons: the basic differences between BUU and QMD models in the representation of the phase-space distributions, computational differences in the mean-field evaluation, and differences in the adopted strategies for the Pauli blocking in the collision integrals. For pionic observables, we find that a higher maximum density leads to an enhanced pion yield and a reduced pi- pi+ yield ratio, while a more effective Pauli blocking generally leads to a slightly suppressed pion yield and an enhanced pi- pi+ yield ratio. We specifically investigate the effect of the Coulomb force, and find that it increases the total pi- pi+ yield ratio but reduces the ratio at high pion energies, although differences in its implementations do not have a dominating role in the differences among the codes. Taking into account only the results of codes that strictly follow the homework specifications, we find a convergence of the codes in the final charged pion yield ratio to a 1 sigma deviation of about 5%. Furthermore, the uncertainty is expected to be reduced to about 1.6% if the same or similar strategies and ingredients, i.e., an improved Pauli blocking and calculation of the non-linear term in the mean-field potential, are similarly used in all codes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

How analytic choices can affect the extraction of electromagnetic form factors from elastic electron scattering cross section data

Scientists often try to incorporate prior knowledge into their regression algorithms, such as a particular analytic behavior or a known value at a kinematic endpoint. Unfortunately, there is often no unique way to make use of this prior knowledge, and thus, different analytic choices can lead to very different regression results from the same set of data. In this paper, to illustrate this point in the context of the proton electromagnetic form factors, we use the Mainz elastic data with its 1422 cross section points and 31 normalization parameters. Starting with a complex unbound non-linear regression, we will show how the addition of a single theory-motivated constraint removes an oscillation from the magnetic form factor and shifts the extracted proton charge radius. We then repeat both regressions using the same algorithm, but with a rebinned version of the Mainz dataset. These examples illustrate how analytic choices, such as the function that is being used or even the binning of the data, can dramatically affect the results of a complex regression. These results also demonstrate why it is critical when using regression algorithms to have either a physical model in mind or a firm mathematical basis

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

End-to-End Differentiable Modeling and Management of the Environment

Focal Area: (2) Data acquisition and assimilation enabled by machine learning, AI, and advanced methods including experimental/network design/optimization. We emphasize the importance of leveraging optimization techniques from AI/machine learning (ML) to solve challenging problems in Earth system modeling. Science Challenge: Automatic differentiation has had a transformative effect on ML by allowing the calculation of gradients of arbitrary functions in an incredibly large class of models. We can potentially realize similar improvements in parameter estimation and control for Earth system models (ESMs) by reimplementing them in computational frameworks from ML. Practitioners working with large (>10 7 parameters) models in ML and AI can obtain good predictive performance in a range of spatiotemporal tasks by making use of optimization via stochastic gradient descent and incorporating prior knowledge at multiple levels. We propose writing ESMs in open-source computational frameworks such as Torch, Tensorflow, and JAX to greatly expand the scope of environmental forecasting and management challenges, which can be addressed by leveraging automatic differentiation and gradient descent-like algorithms. We do not call for a wholesale replacement of physical models with data-driven surrogates, but rather advocate for interleaving physical and empirical equations in a manner that is most faithful to the extent of our scientific knowledge and observational data. Central to this topic is the merging of differentiable physical simulations with differentiable optimization layers, which are now both beginning to come to the forefront.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1997

Measurement and Validation of Bidirectional Reflectance of Space Shuttle and Space Station Materials for Computerized Lighting Models

Task illumination has a major impact on human performance: What a person can perceive in his environment significantly affects his ability to perform tasks, especially in space's harsh environment. Training for lighting conditions in space has long depended on physical models and simulations to emulate the effect of lighting, but such tests are expensive and time-consuming. To evaluate lighting conditions not easily simulated on Earth, personnel at NASA Johnson Space Center's (JSC) Graphics Research and Analysis Facility (GRAF) have been developing computerized simulations of various illumination conditions using the ray-tracing program, Radiance, developed by Greg Ward at Lawrence Berkeley Laboratory. Because these computer simulations are only as accurate as the data used, accurate information about the reflectance properties of materials and light distributions is needed. JSC's Lighting Environment Test Facility (LETF) personnel gathered material reflectance properties for a large number of paints, metals, and cloths used in the Space Shuttle and Space Station programs, and processed these data into reflectance parameters needed for the computer simulations. They also gathered lamp distribution data for most of the light sources used, and validated the ability to accurately simulate lighting levels by comparing predictions with measurements for several ground-based tests. The result of this study is a database of material reflectance properties for a wide variety of materials, and lighting information for most of the standard light sources used in the Shuttle/Station programs. The combination of the Radiance program and GRAF's graphics capability form a validated computerized lighting simulation capability for NASA.

Fletcher, Lauren E.↗
DOE OSTI2022

Computational models of direct and indirect X‐ray breast imaging detectors for in silico trials

Abstract Background To facilitate in silico studies that investigate digital mammography (DM) and breast tomosynthesis (DBT), models replicating the variety in imaging performance of the DM and DBT systems, observed across manufacturers are needed. Purpose The main purpose of this work is to develop generic physics models for direct and indirect detector technology used in commercially available systems, with the goal of making them available open source to manufacturers to further tweak and develop the exact in silico replicas of their systems. Methods We recently reported on an in silico version of the SIEMENS Mammomat Inspiration DM/DBT system using an open‐source GPU‐accelerated Monte Carlo x‐ray imaging simulation code (MC‐GPU). We build on the previous version of the MC‐GPU codes to mimic the imaging performances of two other Food and Drug Administration (FDA)‐approved DM/DBT systems, such as Hologic Selenia Dimensions (HSD) and the General Electric Senographe Pristina (GSP) systems. In this work, we developed a hybrid technique to model the optical spread and signal crosstalk observed in the GSP and HSD systems. MC simulations are used to track each x‐ray photon till its first interaction within the x‐ray detector. On the other hand, the signal spread in the x‐ray detectors is modeled using previously developed analytical equations. This approach allows us to preserve the modeling accuracy offered by MC methods in the patient body, while speeding up secondary carrier transport (either electron–hole pairs or optical photons) using analytical equations in the detector. The analytical optical spread model for the indirect detector includes the depth‐dependent spread and collection of optical photons and relies on a pre‐computed set of point response functions that describe the optical spread as a function of depth. To understand the capabilities of the computational x‐ray detector models, we compared image quality metrics like modulation transfer function (MTF), normalized noise power spectrum (NNPS), and detective quantum efficiency (DQE), simulated with our models against measured data. Please note that the purpose of these comparisons with measured data would be to gauge if the model developed as part of this work could replicate commercially used direct and indirect technology in general and not to achieve perfect fits with measured data. Results We found that the simulated image quality metrics such as MTF, NNPS, and DQE were in reasonable agreement with experimental data. To demonstrate the imaging performance of the three DM/DBT systems, we integrated the detector models with the VICTRE pipeline and simulated DM images of a fatty breast model containing a spiculated mass and a calcium oxalate cluster. In general, we found that the images generated using the indirect model appeared more blurred with a different noise texture and contrast as compared to the systems with direct detectors. Conclusions We have presented computational models of three commercially available FDA‐approved DM/DBT systems, which implement both direct and indirect detector technology. The updated versions of the MC‐GPU codes that can be used to replicate three systems are available in open source format through GitHub.

Sengupta, Aunnasha↗
DOE OSTI2021

Multi-objective Bayesian optimization of ferroelectric materials with interfacial control for memory and energy storage applications

Optimization of materials’ performance for specific applications often requires balancing multiple aspects of materials’ functionality. Even for the cases where a generative physical model of material behavior is known and reliable, this often requires search over multidimensional function space to identify low-dimensional manifold corresponding to the required Pareto front. In this work, we introduce the multi-objective Bayesian optimization (MOBO) workflow for the ferroelectric/antiferroelectric performance optimization for memory and energy storage applications based on the numerical solution of the Ginzburg–Landau equation with electrochemical or semiconducting boundary conditions. MOBO is a low computational cost optimization tool for expensive multi-objective functions, where we update posterior surrogate Gaussian process models from prior evaluations and then select future evaluations from maximizing an acquisition function. Using the parameters for a prototype bulk antiferroelectric (PbZrO 3 ), we first develop a physics-driven decision tree of target functions from the loop structures. We further develop a physics-driven MOBO architecture to explore multidimensional parameter space and build Pareto-frontiers by maximizing two target functions jointly—energy storage and loss. This approach allows for rapid initial materials and device parameter selection for a given application and can be further expanded toward the active experiment setting. The associated notebooks provide both the tutorial on MOBO and allow us to reproduce the reported analyses and apply them to other systems (https://github.com/arpanbiswas52/MOBO_AFI_Supplements).

36 MATERIALS SCIENCE↗
DOE OSTI2023

Optimization with Neural Network Feasibility Surrogates: Formulations and Application to Security-Constrained Optimal Power Flow

In many areas of constrained optimization, representing all possible constraints that give rise to an accurate feasible region can be difficult and computationally prohibitive for online use. Satisfying feasibility constraints becomes more challenging in high-dimensional, non-convex regimes which are common in engineering applications. A prominent example that is explored in the manuscript is the security-constrained optimal power flow (SCOPF) problem, which minimizes power generation costs, while enforcing system feasibility under contingency failures in the transmission network. In its full form, this problem has been modeled as a nonlinear two-stage stochastic programming problem. In this work, we propose a hybrid structure that incorporates and takes advantage of both a high-fidelity physical model and fast machine learning surrogates. Neural network (NN) models have been shown to classify highly non-linear functions and can be trained offline but require large training sets. In this work, we present how model-guided sampling can efficiently create datasets that are highly informative to a NN classifier for non-convex functions. We show how the resultant NN surrogates can be integrated into a non-linear program as smooth, continuous functions to simultaneously optimize the objective function and enforce feasibility using existing non-linear solvers. Overall, this allows us to optimize instances of the SCOPF problem with an order of magnitude CPU improvement over existing methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2012

Assimilation of GRACE Terrestrial Water Storage into a Land Surface Model: Evaluation 1 and Potential Value for Drought Monitoring in Western and Central Europe

A land surface model s ability to simulate states (e.g., soil moisture) and fluxes (e.g., runoff) is limited by uncertainties in meteorological forcing and parameter inputs as well as inadequacies in model physics. In this study, anomalies of terrestrial water storage (TWS) observed by the Gravity Recovery and Climate Experiment (GRACE) satellite mission were assimilated into the NASA Catchment land surface model in western and central Europe for a 7-year period, using a previously developed ensemble Kalman smoother. GRACE data assimilation led to improved runoff correlations with gauge data in 17 out of 18 hydrological basins, even in basins smaller than the effective resolution of GRACE. Improvements in root zone soil moisture were less conclusive, partly due to the shortness of the in situ data record. In addition to improving temporal correlations, GRACE data assimilation also reduced increasing trends in simulated monthly TWS and runoff associated with increasing rates of precipitation. GRACE assimilated root zone soil moisture and TWS fields exhibited significant changes in their dryness rankings relative to those without data assimilation, suggesting that GRACE data assimilation could have a substantial impact on drought monitoring. Signals of drought in GRACE TWS correlated well with MODIS Normalized Difference Vegetation Index (NDVI) data in most areas. Although they detected the same droughts during warm seasons, drought signatures in GRACE derived TWS exhibited greater persistence than those in NDVI throughout all seasons, in part due to limitations associated with the seasonality of vegetation.

Li, Bailing↗
NASA NTRS2014

Land Surface Microwave Emissivity Dynamics: Observations, Analysis and Modeling

Land surface microwave emissivity affects remote sensing of both the atmosphere and the land surface. The dynamical behavior of microwave emissivity over a very diverse sample of land surface types is studied. With seven years of satellite measurements from AMSR-E, we identified various dynamical regimes of the land surface emission. In addition, we used two radiative transfer models (RTMs), the Community Radiative Transfer Model (CRTM) and the Community Microwave Emission Modeling Platform (CMEM), to simulate land surface emissivity dynamics. With both CRTM and CMEM coupled to NASA's Land Information System, global-scale land surface microwave emissivities were simulated for five years, and evaluated against AMSR-E observations. It is found that both models have successes and failures over various types of land surfaces. Among them, the desert shows the most consistent underestimates (by approx. 70-80%), due to limitations of the physical models used, and requires a revision in both systems. Other snow-free surface types exhibit various degrees of success and it is expected that parameter tuning can improve their performances.

Modeling↗
DOE OSTI2021

A model for the fragmentation kinetics of crumpled thin sheets

As a confined thin sheet crumples, it spontaneously segments into flat facets delimited by a network of ridges. Despite the apparent disorder of this process, statistical properties of crumpled sheets exhibit striking reproducibility. Experiments have shown that the total crease length accrues logarithmically when repeatedly compacting and unfolding a sheet of paper. Here, we offer insight to this unexpected result by exploring the correspondence between crumpling and fragmentation processes. We identify a physical model for the evolution of facet area and ridge length distributions of crumpled sheets, and propose a mechanism for re-fragmentation driven by geometric frustration. This mechanism establishes a feedback loop in which the facet size distribution informs the subsequent rate of fragmentation under repeated confinement, thereby producing a new size distribution. We then demonstrate the capacity of this model to reproduce the characteristic logarithmic scaling of total crease length, thereby supplying a missing physical basis for the observed phenomenon.

42 ENGINEERING↗