Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational 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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗
DOE OSTI2024

Exploring the impact of multivalent substitution on high-temperature electrical conductivity in doped lanthanum chromite perovskites: An experimental and ab-initio study

Doped LaCrO 3 perovskites hold promise as robust materials for electrical interconnects and sensor applications in harsh environments. In this work, we investigated the high-temperature behavior of La 1–x Sr x CrO 3 , La 1–x Ca x CrO 3 and La 0.8 Sr 0.2 Cr 1-x Mn x O3 (0.1 ≤ x ≤ 0.4) through a combination of computational modeling and physical characterization up to 1500 °C. Crystalline structural properties were determined and compared with ab-initio calculations, which demonstrated excellent agreement with experimental findings. High-temperature electrical conductivity measurements were performed under different atmospheres. Calcium and strontium/manganese co-doped lanthanum chromites exhibited typical semiconductor exponential trends and conductivity showed proportional correlation with substitutions levels up to 30 %. The DFT modelling was completed up to 1500 °C, including low and high temperature chromite phases and oxygen vacancies insertion. Calculations were correlated with experimental electrical properties. This work expands the understanding of doped lanthanum chromites and paves the way for the development of materials suitable for demanding high-temperature applications.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Elastic nucleon-pion scattering at m π = 200 MeV from lattice QCD

Elastic nucleon-pion scattering amplitudes are computed using lattice QCD on a single ensemble of gauge field configurations with N f = 2 + 1 dynamical quark flavors and m π = 200 MeV. The s-wave scattering lengths with both total isospins I = 1/2 and I = 3/2 are inferred from the finite-volume spectrum below the inelastic threshold together with the I = 3/2 p-wave containing the Δ(1232) resonance. The amplitudes are well-described by the effective range expansion with parameters constrained by fits to the finite-volume energy levels, enabling a determination of the I = 3/2 scattering length with statistical errors below 5%, while the I = 1/2 scattering length is somewhat less precisely evaluated. Systematic errors due to excited states and the influence of higher partial waves are controlled, providing a step toward future computations down to physical light quark masses with multiple lattice spacings and volumes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Surrogate modeling for efficiently, accurately and conservatively estimating measures of risk

We present a surrogate modeling framework for conservatively estimating measures of risk from limited realizations of an expensive physical experiment or computational simulation. Risk measures combine objective probabilities with the subjective values of a decision maker to quantify anticipated outcomes. Given a set of samples, we construct a surrogate model that produces estimates of risk measures that are always greater than their empirical approximations obtained from the training data. These surrogate models limit over-confidence in reliability and safety assessments and produce estimates of risk measures that converge much faster to the true value than purely sample-based estimates. We first detail the construction of conservative surrogate models that can be tailored to a stakeholder’s risk preferences and then present an approach, based on stochastic orders, for constructing surrogate models that are conservative with respect to families of risk measures. Our surrogate models include biases that permit them to conservatively estimate the target risk measures. We provide theoretical results that show that these biases decay at the same rate as the L 2 error in the surrogate model. Numerical demonstrations confirm that risk-adapted surrogate models do indeed overestimate the target risk measures while converging at the expected rate.

42 ENGINEERING↗
DOE OSTI2021

The Rise of Neural Networks for Materials and Chemical Dynamics

Machine learning (ML) is quickly becoming a premier tool for modeling chemical processes and materials. ML-based force fields, trained on large data sets of high-quality electron structure calculations, are particularly attractive due their unique combination of computational efficiency and physical accuracy. This Perspective summarizes some recent advances in the development of neural network-based interatomic potentials. Designing high-quality training data sets is crucial to overall model accuracy. One strategy is active learning, in which new data are automatically collected for atomic configurations that produce large ML uncertainties. Another strategy is to use the highest levels of quantum theory possible. Transfer learning allows training to a data set of mixed fidelity. A model initially trained to a large data set of density functional theory calculations can be significantly improved by retraining to a relatively small data set of expensive coupled cluster theory calculations. These advances are exemplified by applications to molecules and materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Material Fracturing and Failure Simulation Datasets

Fracturing is a fundamental physics phenomena with broad relevance across multiple domains, ranging from infrastructure integrity, aerospace durability, reservoir production, and seismic events. We present a diverse dataset of simulated fracture evolution and material failure generated from two numerical solvers: the phase-field method and the combined finite-discrete element method (FDEM). These solvers differ in formulation, physical fidelity, and computational efficiency. The dataset includes five materials: PBX, anisotropic shale, tungsten, aluminum, and steel. For each, phase-field simulations span 400,000 cases: 200,000 under uniaxial tension and 200,000 under biaxial tension. The computationally expensive FDEM simulations include 90,000 split evenly among PBX, shale, and tungsten under uniaxial loading. All simulations begin with randomized initial fracture patterns. Each entry includes temporal data capturing fracture propagation dynamics. This comprehensive dataset is designed to support the development of foundational or surrogate machine learning approaches for predicting material failure. While no such models are introduced here, the dataset lays a robust foundation for advancing future research and innovation in these areas.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Porting ATLAS Fast Calorimeter Simulation to GPUs with Performance Portable Programming Models

FastCaloSim is a parameterized simulation of the particle energy response and of the energy distribution in the ATLAS calorimeter. It is a relatively small and self-contained package with massive inherent parallelism and captures the essence of GPU offloading via important operations like data transfer, memory initialization, floating point operations, and reduction. It was identified by the High Energy Physics Center for Computational Excellence project as a good testbed for evaluating the performance and ease of portability of programming models. In this paper, we will discuss the results of our evaluation of the porting process to Kokkos, SYCL, Alpaka, OpenMP and std::par (nvc++), and compare performance on NVIDIA, AMD and Intel GPUs, as well as multicore CPUs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Operator-level quantum acceleration of non-logconcave sampling

Sampling from probability distributions of the form 𝝈 ∝ e −𝜷V , where V is a continuous potential, is a fundamental task across physics, chemistry, biology, computer science, and statistics. However, when V is nonconvex, the resulting distribution becomes non-logconcave, and classical methods such as Langevin dynamics often exhibit poor performance. We introduce a quantum algorithm that provably accelerates a broad class of continuous-time sampling dynamics. For Langevin dynamics, our method encodes the target Gibbs measure into the amplitudes of aquantum state, identified as the kernel of a block matrix derived from a factorization of the Witten Laplacian operator. This connection enables Gibbs sampling via singular value thresholding and yields up to a quartic quantum speedup over best-knownclassical Langevin-based methods in the non-logconcave setting. Building on this framework, we further develop the first quantum algorithm that accelerates replica exchange Langevin diffusion, a widely used method for sampling from complex, rugged energy landscapes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Preliminary Introduction of a Free Vortex Wake Method Into OpenFAST

Increasing wind turbine size can have substantial economic benefits. However, increasing wind turbine turbine size requires the use of models that accurately capture the aerodynamic response to large deflections and nonstraight blade geometries. Additionally, the installation of offshore floating turbines requires turbine models that can handle large rotor motion. Free vortex wake (FVW) methods are capable of modeling this type of complex physics while remaining computationally tractable to perform the many simulations necessary for turbine design. In this work, a FVW model is added to the National Renewable Energy Laboratory engineering tool OpenFAST to allow for the aerodynamic modeling of highly flexible turbines along with the aero-hydro-servo-elastic response capabilities of OpenFAST. Code implementation is not complete and thus no results are presented. However, areas of potential modeling improvement are discussed, such as modifications to the circulation calculation and the need for different FVW input model specifications when high-thrust cases are being modeled.

17 WIND ENERGY↗
DOE OSTI2023

Microscopic activated dynamics theory of the shear rheology and stress overshoot in ultradense glass-forming fluids and colloidal suspensions

We formulate a particle and force level, activated dynamics-based statistical mechanical theory for the continuous startup nonlinear shear rheology of ultradense glass-forming hard sphere fluids and colloidal suspensions in the context of the elastically collective nonlinear Langevin equation approach and a generalized Maxwell model constitutive equation. Activated structural relaxation is described as a coupled local-nonlocal event involving caging and longer range collective elasticity which controls the characteristic stress relaxation time. Theoretical predictions for the deformation-induced enhancement of mobility, the onset of relaxation acceleration at remarkably low values of stress, strain, or shear rate, apparent power law thinning of the steady-state structural relaxation time and viscosity, a nonvanishing activation barrier in the shear thinning regime, an apparent Herschel–Buckley form of the shear rate dependence of the steady-state shear stress, exponential growth of different measures of a yield or flow stress with packing fraction, and reduced fragility and dynamic heterogeneity under deformation were previously shown to be in good agreement with experiments. The central new question we address here is the defining feature of the transient response—the stress overshoot. In contrast to the steady-state flow regime, understanding the transient response requires an explicit treatment of the coupled nonequilibrium evolution of structure, elastic modulus, and stress relaxation time. We formulate a new quantitative model for this aspect in a physically motivated and computationally tractable manner. Theoretical predictions for the stress overshoot are shown to be in good agreement with experimental observations in the metastable ultradense regime of hard sphere colloidal suspensions as a function of shear rate and packing fraction, and accounting for deformation-assisted activated motion appears to be crucial for both the transient and steady-state responses.

Mechanics↗
DOE OSTI2021

Energy Exascale Earth System Model v2.0

First release of version 2 of the Energy Exascale Earth System Model. The atmosphere component remains EAM. Major changes since version 1 include: all column-physics parameterizations are computed on a separate grid that has approximately half the number of points of the dynamics grid, a new nonhydrostatic dynamical core (running in hydrostatic mode) with semi-Lagrangian tracer transport, CLUBB updated from v1 to v2, a new convective trigger (dCAPE/ULL) based on the dynamic Convective Available Potential Energy (CAPE) (dCAPE) and the Unrestricted Launch Level (ULL) concepts is used in ZM. minimum cloud droplet number changed, gravity wave drag energy conservation fixed and new tunings used, dust emission size distribution changed to emit more coarse dust particles The land component is still ELM. Major changes since version 1 include: using SNICAR-AD for radiation in snow to match the sea-ice model and fixing bugs in snow compaction and water state calculation. The ocean component remains MPAS-ocean. Major change since version 1 include: Redi isopycnal mixing has been updated, tested, and tuned in combination with the Gent-McWilliams parameterization, a sign error was fixed in the 3rd-order flux routines, the mesh used in low-resolution coupled cases was modified and the time steps adjusted, new regionally refined meshes were created, one focused on North America and another on the Southern Ocean, including ice shelf cavities. The sea-ice component remains MPAS-seaice. Major changes since version 1 include: A new heat- and freshwater-conserving coupling of frazil ice, turning off of SSH filtering, addition of SNICAR-AD and snow grain aging. The land-ice component remains MPAS-Albany-landIce (MALI) and is a static ice sheet. There is more out-of-the-box support for cryosphere configurations including the new regionally refined configuration around Antarctica with ice shelf cavities. The river model is MOSART. The half-degree river mesh was redone so it no longer treats Black and Caspian seas as ocean. The coupler remains cpl7/MCT. Major changes since version 1 include handling of ice shelf melt fluxes (heat / freshwater exchange with the ocean), and data icebergs. All components allow regional refinement of their meshes and two separate example of refinement, one in and around North America and one around Antarctica and the Southern Ocean, are provided.

E3SM Project, DOE↗
DOE OSTI2021

Diffraction properties of cylindrically bent KAP crystals in energy range of 2.3–7.5 keV using synchrotron radiation

Verification of physics models and computer simulations are heavily reliant upon the accuracy of experimental measurements. Calibration of instrument responses becomes an important step to achieve this goal. This paper presents systematic studies of bent potassium acid phthalate (KAP) crystals using Lawrence Berkeley National Laboratories, Advanced Light Source, beamline 9.3.1 in the energy range of 2.3 to 7.5 keV. A set of KAP crystals, gradually bent from flat up to a 50.8 mm cylindrical curvature. The measured integrated reflectivity for this set of KAP crystals shows good agreement with the X-ray Oriented Program (XOP) calculations when adjusting the Debye–Waller temperature factor and using the multilamellar model in the calculations. Significant differences in rocking curve profiles were observed between experimental measurements and theory. A forward convolution model and software code were developed to include experimental parameters, allowing the investigation of the difference between measurements and calculations. After considering the experimental parameters, good agreements were obtained for the rocking curve profiles for all bending radii with a unique set of parameters. Our results show that XOP can be a useful and reliable tool to predict performance of cylindrically bent KAP crystals in this energy range.

Wu, Ming↗
DOE OSTI1993

Nuclear Safety [Vol. 34, No. 1, January-March 1993]

Nuclear Safety is a review journal that covers significant developments in the field of nuclear safety. Its scope includes the analysis and control of hazards associated with nuclear energy, operations involving fissionable materials, and the products of nuclear fission and their effects on the environment. Primary emphasis is on safety in reactor design, construction, and operation; however, the safety aspects of the entire fuel cycle, including fuel fabrication, spent-fuel processing, nuclear waste disposal, handling of radioisotopes, and environmental effects of these operations, are also treated. Table of Contents for this issue follows. GENERAL SAFETY CONSIDERATIONS: 1 An Integrated Regional Approach To Risk Management of Industrial Systems, S. Chakraborty and R. Stratton; 9 Annual Technical Meeting of the NRC Cooperative Severe Accident Research Program, E. G. Silver; ACCIDENT ANALYSIS: 13 Deterministic Severe Accident Criteria as Severe Accident Design Criteria and Policy for the New Production Reactor-Heavy Water Reactor, P. T. Rhoads; 20 Proposed Deterministic Severe Accident Criteria for the Heavy Water Reactor-New Production Reactor Containment, K. D. Bergeron, S. E. Slezak, and C. E. Leach; 33 Quasi-Static Core Liquid Level Depression and Long-Term Core Uncovery During a PWR LOCA, Y. Kukita, R. R. Schultz, H. Nakamura, and J. Katayama; 49 LIRA: An Advanced Containment System to Minimize the Accidental Radioactivity Releases, A. Turricchia; 63 Errata to “A Review of Hydrogen Production During Melt/Water Interaction in LWRs,” Vol. 33, No. 4; CONTROL AND INSTRUMENTATION: 64 The Nuplex 80+™ Advanced Control Complex from ABB Combustion Engineering, F. Ridolfo, D. Harmon, and K. Scarola; DESIGN FEATURES: 76 A Scheme for Passive Isolation of the Containment of a Reactor, A. K. Ghosh, V. V. Raj, and K. Kakodkar; 84 Utility Requirements for Safety in the Passive Advanced Light-Water Reactor, T. U. Marston, W. H. Layman, and G. Bockhold, Jr.; WASTE AND SPENT FUEL MANAGEMENT: 91 Activities Related to Waste and Spent Fuel Management, Compiled by M. D. Muhlheim and E. G. Silver; OPERATING EXPERIENCES: 103 Testing Deficiencies in Auxiliary Systems Feedwater Systems, J. D. Kueck; 110 Reactor Shutdown Experience, Compiled by J. W. Cletcher; 113 Selected Safety-Related Events, Compiled by G. A. Murphy; 115 Operating U.S. Power Reactors, Compiled by M. D. Muhlheim and E. G. Silver; RECENT DEVELOPMENTS: 133 General Administrative Activities, Compiled by M. D. Muhlheim and E. G. Silver; 145 Reports, Standards, and Safety Guides, D. S. Queener; 150 Proposed Rule Changes as of Sept. 30, 1992; ANNOUNCEMENTS: 75 Harvard Short Course on Risk: Science, Assessment, and Management; 102 International Conference on Reactor Physics and Reactor Computations; 102 The Fourth International Conference on Nuclear Waste Reprocessing and Waste Management [RECOD ’94]; 132 International Symposium on the Radiological Exposure of the Population of the European Community from Radioactivity in the Mediterranean Sea [MARINA-MED]; 162 International High-Level Radioactive Waste Conference ‘Technology Integration—Above and Below”; 154 The Authors; 158 Indexes to Nuclear Safety, Volume 33.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Baylor University Campus-Wide Deep Dive

In January 2020, staff members from the Engagement and Performance Operations Center (EPOC) and the Lonestar Education And Research Network (LEARN) met with researchers and staff at Baylor University for the purpose of a Campus-Wide Deep Dive into research drivers. The goal of this meeting was to help characterize the requirements for five campus research use cases and to enable cyberinfrastructure support staff to better understand the needs of the researchers they support. Profiled scientific use cases included: - Experimental High Energy Physics (HEP) - Proton Computed Tomography (pCT) - Nutrition and Relation to Digestive Microbiome - Baylor University Core Research Facilities - Molecular Quantum-dot Cellular Automata (QCA), and Material Science of Quantum Computing - Modeling and Simulation of Low-Dimensional and Nano-Structured Materials - Computational Fluid Dynamics Material for this event included the written documentation from each of the research areas at Baylor University, documentation about the current state of technology support, and a write-up of the discussion that took place in person. The Case Studies highlighted the ongoing challenges that Baylor University has in supporting a cross-section of established and emerging research use cases. Each Case Study mentioned unique challenges which were summarized into common needs. These included: - Tradeoffs for network/software security, and usability of the resulting infrastructure. Better communication to set expectations and understand realities is required. - Computation use on campus is widespread and healthy. While no major problems were uncovered, upgrades to maintain current usage patterns and encourage growth will be required. - Storage is a critical need for enterprise use cases and research. In particular, a campus wide ‘storage architecture’ to support research use cases (e.g. instruments, data sharing) is required in the 2-5 year time window. - Instrumentation on campus is healthy and expanding. Technology must scale with this in the form of computation and storage. - Working with LEARN to upgrade network capacity (in multiples of 10G, or upgrades to 100G) will be required in the 1-3 year time frame. - Network monitoring and visibility will help to establish external science use cases. - Data sharing via portal systems is not currently a critical need, but growing in scope. EPOC can assist Baylor with options.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

Surrogate Modeling For Efficiently Accurately and Conservatively Estimating Measures of Risk

We present a surrogate modeling framework for conservatively estimating measures of risk from limited realizations of an expensive physical experiment or computational simulation. We adopt a probabilistic description of risk that assigns probabilities to consequences associated with an event and use risk measures, which combine objective evidence with the subjective values of decision makers, to quantify anticipated outcomes. Given a set of samples, we construct a surrogate model that produces estimates of risk measures that are always greater than their empirical estimates obtained from the training data. These surrogate models not only limit over-confidence in reliability and safety assessments, but produce estimates of risk measures that converge much faster to the true value than purely sample-based estimates. We first detail the construction of conservative surrogate models that can be tailored to the specific risk preferences of the stakeholder and then present an approach, based upon stochastic orders, for constructing surrogate models that are conservative with respect to families of risk measures. The surrogate models introduce a bias that allows them to conservatively estimate the target risk measures. We provide theoretical results that show that this bias decays at the same rate as the L 2 error in the surrogate model. Our numerical examples confirm that risk-aware surrogate models do indeed over-estimate the target risk measures while converging at the expected rate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Biokinetic and Dosimetric Models [Slides]

Presentation describes internal dosimetry as an intersection of disciplines involving Physiology, Anatomy, Physics, Mathematics and Computer Science, which some people call "The art of Internal Dosimetry". The basic components for internal dose calculations are: Biokinetic (metabolic) models describing the intake, distribution, retention and excretion of radionuclides in the body, Dosimetric models describing the interaction of radiation within the several body organs and tissues, and System of dose limitation.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2022

Resilience Enhancements through Deep Learning Yields

This report documents the Resilience Enhancements through Deep Learning Yields (REDLY) project, a three-year effort to improve electrical grid resilience by developing scalable methods for system operators to protect the grid against threats leading to interrupted service or physical damage. The computational complexity and uncertain nature of current real-world contingency analysis presents significant barriers to automated, real-time monitoring. While there has been a significant push to explore the use of accurate, high-performance machine learning (ML) model surrogates to address this gap, their reliability is unclear when deployed in high-consequence applications such as power grid systems. Contemporary optimization techniques used to validate surrogate performance can exploit ML model prediction errors, which necessitates the verification of worst-case performance for the models.

24 POWER TRANSMISSION AND DISTRIBUTION↗