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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A seamless multiscale operator neural network for inferring bubble dynamics

Modelling multiscale systems from nanoscale to macroscale requires the use of atomistic and continuum methods and, correspondingly, different computer codes. Here, we develop a seamless method based on DeepONet, which is a composite deep neural network (a branch and a trunk network) for regressing operators. In particular, we consider bubble growth dynamics, and we model tiny bubbles of initial size from 100 nm to 10 $\mathrm {\mu }\textrm {m}$ , modelled by the Rayleigh–Plesset equation in the continuum regime above 1 $\mathrm {\mu }\textrm {m}$ and the dissipative particle dynamics method for bubbles below 1 $\mathrm {\mu }\textrm {m}$ in the atomistic regime. After an offline training based on data from both regimes, DeepONet can make accurate predictions of bubble growth on-the-fly (within a fraction of a second) across four orders of magnitude difference in spatial scales and two orders of magnitude in temporal scales. The framework of DeepONet is general and can be used for unifying physical models of different scales in diverse multiscale applications.

Mechanics↗
DOE OSTI2023

Accurate estimation of angular power spectra for maps with correlated masks

A common procedure when analyzing maps of the cosmic microwave background (CMB) or other cosmological signals is the need to remove ("mask") regions of the maps that are heavily contaminated, e.g., by non-cosmological foreground emission. After applying such a mask, one must account for its effect when inferring statistical properties of interest, such as the angular power spectrum of the field in the original map. A widely used approach to correct for such mask-induced effects was presented by Hivon et al. (2002), now widely known as the "MASTER" formalism. However, it is often the case that the map and mask are correlated in some way, such as point source masks used in CMB analyses, which have nonzero correlation with CMB secondary anisotropy fields and other mm-wave sky signals. In such situations, the MASTER approach gives biased results, as it assumes that the unmasked map and mask have zero correlation. While such effects have been discussed before with regard to specific physical models, here we derive a completely general formalism for any case where the map and mask are correlated. We show that our result ("reMASTERed") reconstructs ensemble-averaged angular power spectra to effectively exact precision, with significant improvements over traditional estimators for cases where the map and mask are correlated. An important consequence of our result is that for maps with correlated masks, it is no longer possible to invert a simple equation to obtain the true power spectrum from the observed (masked) power spectrum. Instead, our result necessitates the use of forward modeling from theory space into the observable domain of the masked power spectrum. We publicly release our software implementation of these results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Integration of seismic-pressure-petrophysics inversion of continuous active-seismic monitoring data for monitoring and quantifying CO 2 plume (Final Report)

The overall objective of this project is to develop and validate an integrated package of joint seismic-pressure-petrophysics inversion (jSPPI) of continuous active-source seismic monitoring dataset capable of providing real-time monitoring of CO 2 plume during geologic carbon sequestration (GCS). The three specific developments include: (a) the methodologies for fast seismic full waveform inversion of continuous active source seismic monitoring, (CASSM) datasets for simultaneously estimating velocity and attenuation, and with data assimilation; (b) joint Bayesian petrophysical inversion of seismic models and pressure data for providing and updating CO 2 saturation models; (c) the methods using multiple datasets including (Crainfield and Frio-II borehole) synthetic, laboratory, and field CASSM datasets. The outcomes of jSPPI include (a) a workflow for processing CASSM data, (b) Bayesian inversion algorithms using CASSM data and pressure response data, and (c) integration with data assimilation algorithms for continuously updating site-specific models used for prediction and reservoir management. The validation of joint FWI will be conducted using synthetic models based on the Cranfield and Frio experiments as well as field CASSM datasets collected as part of the Frio-II pilot injection. To quantify and map the mass and distribution of CO 2 (saturation), we will jointly invert velocity and attenuation measurements from the FWI with a Bayesian approach using a rock physics model for attenuation (e.g., White’s attenuation model with two selected patch sizes (White, 1976; Dutta and Seriff, 1979)). The Bayesian inversion will be applied to each time step in the CASSM survey in an updating scheme, which integrates with an ensemble of reservoir simulations at each step. A more complete experimental validation dataset will be collected as part of a mesoscale (2-3 m) gas-CO 2 injection experiment utilizing a higher frequency version of the CASSM system developed for laboratory studies; the integrated inversion will be demonstrated using this dataset which will provide both a dense geometry as well as more precise secondary confirmation measurements (e.g. saturation) typically not available in the field. The resulting real-time map of CO 2 saturation is able to provide a deeper scientific understanding of the complex, time-varying dynamics of subsurface fluid flow migration path as well as the rapid detection of CO 2 leakage hazards.

25 ENERGY STORAGE↗
DOE OSTI2022

Pulse-level noisy quantum circuits with QuTiP

The study of the impact of noise on quantum circuits is especially relevant to guide the progress of Noisy Intermediate-Scale Quantum (NISQ) computing. In this paper, we address the pulse-level simulation of noisy quantum circuits with the Quantum Toolbox in Python (QuTiP). We introduce new tools in qutip-qip, QuTiP's quantum information processing package. These tools simulate quantum circuits at the pulse level, leveraging QuTiP's quantum dynamics solvers and control optimization features. We show how quantum circuits can be compiled on simulated processors, with control pulses acting on a target Hamiltonian that describes the unitary evolution of the physical qubits. Various types of noise can be introduced based on the physical model, e.g., by simulating the Lindblad density-matrix dynamics or Monte Carlo quantum trajectories. In particular, the user can define environment-induced decoherence at the processor level and include noise simulation at the level of control pulses. We illustrate how the Deutsch-Jozsa algorithm is compiled and executed on a superconducting-qubit-based processor, on a spin-chain-based processor and using control optimization algorithms. We also show how to easily reproduce experimental results on cross-talk noise in an ion-based processor, and how a Ramsey experiment can be modeled with Lindblad dynamics. Finally, we illustrate how to integrate these features with other software frameworks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A High-Fidelity Model of the Peach Bottom 2 Turbine-Trip Benchmark Using VERA

This work presents a high-fidelity simulation of the Peach Bottom turbine trip (PBTT) benchmark using the Virtual Environment for Reactor Applications (VERA), a multiphysics reactor modeling tool developed by the U.S. Department of Energy’s Consortium for Advanced Simulation of Light Water Reactors energy innovation hub. The PBTT benchmark, based on a 1977 transient event at the end of cycle 2 in a General Electric Type-4 boiling water reactor (BWR), is a critical test case for validating core physics models with thermal feedback during rapid reactivity events. VERA was employed to perform end-to-end, pin-resolved simulations from conditions at the beginning of cycle 1 through the turbine-trip transient, incorporating detailed neutron transport, fuel depletion, and subchannel thermal hydraulics. The simulation reproduced key benchmark observables with high accuracy: the peak power excursion occurred at 0.75 s, matching the scram time and closely aligning with the benchmark average of 0.742 s; the simulated maximum power spike was approximately 7600 MW, which is within 3% of the benchmark average of 7400 MW; and void-collapse dynamics were consistent with benchmark expectations. Reactivity predictions during cycles 1 and 2 remained within 1500 pcm and 400 pcm of criticality, respectively. These results confirm VERA’s ability to model complex coupled neutronic and thermal hydraulic behavior in a BWR turbine-trip transient, which will support its use in future studies of modeling dryout, fuel performance, and uncertainty quantification for transients of this type.

BWR↗
NASA NTRS2001

New Third-Order Moments for the CBL

Turbulent convection is an inherently non-local phenomenon and a primary condition for a successful treatment of the CBL (convective boundary layer) is a reliable model of non-locality. In the dynamic equations governing the convective flux, the turbulent kinetic energy, etc., non-locality is represented by the third-order moments, TOMs. Since the simplest form, the so-called down gradient approximation (DGA), severely underestimates the TOMs (up to an order of magnitude), a more physical model is needed. In 1994, an analytical model was presented which was derived directly from the dynamical equations for the TOMs. It considerably improved the DGA but was a bit cumbersome to use and, more importantly, it was based on the quasi-normal (QN) approximation for the fourth-order moments. Here, we present a new analytic expression for the TOMs which is structurally simpler than the 1994 expression and which avoids the QN approximation. The resulting fit to the LES data is superior to that of the 1994 model.

Canuto, V. M.↗
NASA NTRS2022

Clumped-Isotope Constraint on Upper-Tropospheric Cooling During the Last Glacial Maximum

Ice cores and other paleotemperature proxies, together with general circulation models, have provided information on past surface temperatures and the atmosphere's composition in different climates. Little is known, however, about past temperatures at high altitudes, which play a crucial role in Earth's radiative energy budget. Paleoclimate records at high-altitude sites are sparse, and the few that are available show poor agreement with climate model predictions. These disagreements could be due to insufficient spatial coverage, spatiotemporal biases, or model physics; new records that can mitigate or avoid these uncertainties are needed. Here, we constrain the change in upper-tropospheric temperature at the global scale during the Last Glacial Maximum (LGM) using the clumped-isotope composition of molecular oxygen trapped in polar ice cores. Aided by global three-dimensional chemical transport modeling, we exploit the intrinsic temperature sensitivity of the clumped-isotope composition of atmospheric oxygen to infer that the upper troposphere (effective mean altitude 10–11 km) was 6–9°C cooler during the LGM than during the late preindustrial Holocene. A complementary energy balance approach supports a minor or negligible steepening of atmospheric lapse rates during the LGM, which is consistent with a range of climate model simulations. Proxy-model disagreements with other high-altitude records may stem from inaccuracies in regional hydroclimate simulation, possibly related to land-atmosphere feedbacks.

clumped isotope↗
NASA NTRS1974

Simulated lightning test shuttle .03 scale model

Lightning Attach Point tests were conducted for the space shuttle launch configuration (Orbiter, External Tank and Solid Rocket Boosters). A series of 250 long spark tests (15 to 20 foot sparks) determined that the orbiter may be struck on the nose, windshield brow, tail and wingtips during launch but not on the main engine nozzles which have been shown to be vulnerable to lightning damage. The orbiter main engine and SRB exhaust plumes were simulated electrically with physical models coated with graded resistance paints. The tests showed that the exhaust plumes from the SRB provide additional protection for the main engine nozzles. However, the tests showed that the Orbiter Thermal Protection System (TPS), which has also been shown to be vulnerable to lightning damage, may be struck during launch. Therefore further work is indicated in the areas of swept stroke studies on the model and on TPS panels. Further attach point testing is also indicated on the free-flying orbiter. Photographs of the test setup are shown.

Clifford, D. W.↗
NASA NTRS1988

Numerical Modeling Of Two-Phase, Reactive Flows

Capabilities of several computational schemes combined and enhanced. Advanced Rocket Injector-Combustor Code is computer program for numerical simulation of two-phase flows of chemically reacting or nonreacting components. Applied to sprays of fuel in injectors and combustors, flows in heat exchangers, evaporators, and mass-diffusion systems; and two-dimensional or axisymmetric flows containing liquid or solid suspensions. Used in research for development and verification of various physical models.

Liang, Pak-Yan↗
DOE OSTI2021

MFIX-Exa: A path toward exascale CFD-DEM simulations

MFIX-Exa is a computational fluid dynamics–discrete element model (CFD-DEM) code designed to run efficiently on current and next-generation supercomputing architectures. MFIX-Exa combines the CFD-DEM expertise embodied in the MFIX code—which was developed at NETL and is used widely in academia and industry—with the modern software framework, AMReX, developed at LBNL. The fundamental physics models follow those of the original MFIX, but the combination of new algorithmic approaches and a new software infrastructure will enable MFIX-Exa to leverage future exascale machines to optimize the modeling and design of multiphase chemical reactors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The Steel Equivalency Workbook: An X-Ray Transmission Calculator

The Steel Equivalency Spreadsheet was created as a robust, user-friendly method for the following: 1. Determining if it is possible to image an object. 2. Reducing time needed when selecting equipment and preparing for field imaging activities. 3. Reducing the amount of equipment taken for field imaging activities. 4. Determining a starting point for exposure settings prior to imaging an object. 5. Reducing the amount of dose deposited to complete an imaging operation (in the spirit of As Low As Reasonably Achievable, ALARA). The spreadsheet was developed using various physical models to account for different phenomena. Future implementations aim to expand beyond the Digital Radiography and Computed Tomography Single Munition Scanner (DRCT SMS) in standard configuration to include high energy XRGs (Betatrons) and sub-MeV pulsed XRGs (XRS4). Alternative computational methods that compensate for incoherent scattering when calculating relative transmission are also being pursued. Ultimately the spreadsheet exceeded the developmental goal of having less than 10% average error when comparing calculations to real-world images.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Overcoming the Technical Challenges of Coordinating Distributed Load Resources at Scale (Final Report)

Significant recent research has investigated the potential for loads to provide balancing services to the grid. However, this research has not addressed key issues that may arise when such schemes are applied at scale including: 1. Distribution Network Issues. Coordination of large numbers of loads could result in power flows that violate distribution network constraints; 2. Stability Issues. Certain strategies to control loads can exhibit nonlinearity in the form of period-adding bifurcations and chaos. Other control strategies can potentially synchronize the behavior of large numbers of loads. In both cases, the outcome can be power oscillations and instability; 3. Communication Network Issues. Bidirectional low-latency communication channels between a central controller (or several distributed controllers) and each resource are expensive and likely not necessary for effective coordination. Our research questions were: What network, stability, and communication issues might arise in practice when we coordinate large aggregations of loads? How can we coordinate loads to achieve performance objectives in a cost effective manner while avoiding these issues? The ultimate technical goal of the project was the development of network-aware, communication-constrained, non-disruptive load control strategies with stability guarantees that achieve the performance requirements of typical balancing services at a sufficiently low cost to enable the load aggregator and customer to profit. The overall goal was to establish credibility for load control at scale and contribute to U.S. energy security and environmental goals. The team succeeded in answering these research questions and developing these control strategies. The overall approach was based on the development of three testing environments: a simulation testbed, an experimental testbed (20 physical model houses with window-box air conditioners) coupled with the simulation testbed, and a field testbed (100 actual homes in Austin, TX) coupled with the simulation testbed, which enabled controller testing, identification of issues, controller development, and controller validation. The resulting controller was used to demonstrate fast timescale grid balancing (frequency regulation) by aggregations of physical and virtual air conditioners, with sufficient quality to participate in the electricity market. Cost benefit analysis showed overall benefits to the participating households, load aggregators, and the grid, especially if the control technology was integrated directly into existing programmable communicating thermostats. The project provides a variety of wider benefits. Our technology transfer and outreach activities lead us to choose an open-source licensing commercialization pathway, enabling the project results to be available to researchers, industry, and the public. Furthermore, new grid balancing technologies will increase grid flexibility and will enable higher penetrations of intermittent renewable energy resources, such as wind and solar, to be connected to the grid, reducing its environmental impact, and mitigating climate change to the benefit of society.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Neutron’s dark secret

The existing discrepancy between neutron lifetime measurements in bottle and beam experiments has been interpreted as a sign of the neutron decaying to dark particles. We summarize the current status of this proposal, including a discussion of particle physics models involving such a portal between the Standard Model and a baryonic dark sector. We also review further theoretical developments around this idea and elaborate on the prospects for verifying the neutron dark decay hypothesis in current and upcoming experiments.

Astronomy & Astrophysics↗
NASA NTRS2003

The Electric Propulsion Interactions Code (EPIC): A Member of the NASA Space Environment and Effects Program (SEE) Toolset

Science Applications International Corporation is currently developing the Electric Propulsion Interactions Code, EPIC, as part of a project sponsored by the Space Environments and Effects Program at NASA Marshall Space Flight Center. Now in its second year of development, EPIC is an interactive computer toolset that allows the construction of a 3-D spacecraft model, and the assessment of a variety of interactions between its subsystems and the plume from an electric thruster. This paper reports on the progress of EPZC including the recently added ability to exchange results the NASA Charging Analyzer Program, Nascap-2k. The capability greatly enhances EPIC's range of applicability. Expansion of the toolset's various physics models proceeds in parallel with the overall development of the software. Also presented are recent upgrades of the elastic scattering algorithm in the electric propulsion Plume Tool. These upgrades are motivated by the need to assess the effects of elastically scattered ions on the SIC for ion beam energies that exceed loo0 eV. Such energy levels are expected in future high-power (>10 kW) ion propulsion systems empowered by nuclear sources.

Mikellides, Ioannis G.↗
DOE OSTI2021

Transport model predictions of 225 Ac production cross sections via energetic p, d and α irradiation of 232 Th targets

Monte Carlo transport codes PHITS and MCNP6 were used to calculate the production cross sections of 225,227 Ac, 227,229 Th, 223,225 Ra, and 229,230,231 Pa via the bombardment of a 232 Th target with energetic protons, deuterons, and α-particles. The incident projectile energies ranged between 10 and 800 MeV/nucleon. When possible, the predicted production cross sections were compared with the available experimental data and other predictions. Here, the degree of the codes’ abilities to match the measured data provides a qualitative assessment of the codes’ abilities to predict data from similar, but unmeasured, projectile/target systems. In addition, a comparison between calculated cross sections and data may provide insight into possible improvements in the physics models employed by those transport codes.

223Ra↗
NASA NTRS2017

A Thermal Runaway Failure Model for Low-Voltage BME Ceramic Capacitors with Defects

Reliability of base metal electrode (BME) multilayer ceramic capacitors (MLCCs) that until recently were used mostly in commercial applications, have been improved substantially by using new materials and processes. Currently, the inception of intrinsic wear-out failures in high quality capacitors became much greater than the mission duration in most high-reliability applications. However, in capacitors with defects degradation processes might accelerate substantially and cause infant mortality failures. In this work, a physical model that relates the presence of defects to reduction of breakdown voltages and decreasing times to failure has been suggested. The effect of the defect size has been analyzed using a thermal runaway model of failures. Adequacy of highly accelerated life testing (HALT) to predict reliability at normal operating conditions and limitations of voltage acceleration are considered. The applicability of the model to BME capacitors with cracks is discussed and validated experimentally.

reliability↗