Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative limiters”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Physics-Driven Deep-Learning Network for Solving Nonlinear Inverse Problems

Geosteering inversion, which can be viewed as a nonlinear inverse problem, is an important technique used by directional drilling. Traditional methods that rely on iterative procedures and regularization are sensitive to the selection of initial values and can be slow due to convergence issues. In industrial applications, a lookup table is used to produce fast predictions. However, performance is not guaranteed by this approach due to the limitation of the hardware. In this paper, we propose a novel physics-driven deep-learning framework for providing a fast, accurate surrogate to solve the inverse problem. Particularly, leveraged by the forward physical model and 1D convolutional neural network (1D-CNN), the proposed method provides more reliable solutions to the inverse problem with improved performance. In addition, a new physics-driven loss function is introduced to accommodate both the model misfit and the data misfit. In conclusion, our experiments demonstrate the effectiveness of our method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A New Shutdown Dose Rate Benchmark Problem for Representative Fusion Applications

Here, this work introduces a new benchmark problem for calculating shutdown dose rates (SDDRs) aimed at fusion reactor applications. The model is designed to represent a simplified version of a typical ITER port plug. The responses of interest include neutron flux, gamma flux, and gamma SDDR at 12 different locations scattered throughout the port. This article outlines the geometry specifications of the problem, provides material definitions for the components, specifies the required responses to be calculated, and presents the source definition information. The need for this benchmark arises from the limited availability of publicly accessible references, with only one benchmark representing the typical dimensions and materials found in fusion systems. This existing benchmark has been cited extensively, reflecting the demand within the scientific community to test both established and novel workflows for SDDR calculations. However, since its presentation at a conference in 2011, the results have become increasingly well known. Moreover, the absence of formal publication and peer review has led to the details of this benchmark being extracted from secondary sources, such as subsequent studies that reference it. As a result, analysts are left with significant flexibility in interpreting the key parameters, which can be adjusted to account for unknown systematic errors, ultimately reproducing the already well-known responses. This new benchmark serves as an updated version of that earlier work, with the aim of providing a more reliable description of the materials and their impurities, which is crucial for assessing activation and subsequent gamma emission. Additionally, it seeks to provide a geometry that more closely represents an ITER port plug. The improvements in the problem definition will lead to a more reproducible benchmark problem, while also presenting the radiation transport community with a completely new challenge. The results will be published in a future article to allow analysts adequate time to analyze this problem independently.

Benchmark↗
DOE OSTI2020

Scalable Asynchronous Domain Decomposition Solvers

We discuss how parallel implementations of linear iterative solvers generally alternate between phases of data exchange and phases of local computation. Increasingly large problem sizes and more heterogeneous compute architectures make load balancing and the design of low latency network interconnects that are able to satisfy the communication requirements of linear solvers very challenging tasks. In particular, global communication patterns such as inner products become increasingly limiting at scale. We explore the use of asynchronous communication based on one-sided Message Passing Interface primitives in the context of domain decomposition solvers. In particular, a scalable asynchronous two-level Schwarz method is presented. We discuss practical issues encountered in the development of a scalable solver and show experimental results obtained on a state-of-the-art supercomputer system that illustrate the benefits of asynchronous solvers in load balanced as well as load imbalanced scenarios. Using the novel method, we can observe speedups of up to four times over its classical synchronous equivalent.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Accelerating nuclear fuel development and qualification: Modeling and simulation integrated with separate-effects testing

In this work, an approach to transform and accelerate nuclear fuel development and qualification is outlined. The approach exploits advanced modeling and simulation at the outset to inform constituent and system selection and to enable integral fuel performance analyses. Analyses using these tools identify and prioritize the most important fuel performance parameters and phenomena for subsequent targeted characterization with separate-effects tests. Separate-effects testing spans out-of-pile and in-pile tests and is meant to iterate with and inform engineering-scale integral fuel performance analyses throughout the development process. Exercising this cycle in an agile fashion will increase confidence in the integral fuel performance predictions while reducing uncertainties. This process sets the stage for executing a much more limited set of well-defined integral irradiation tests designed to validate engineering-scale fuel performance codes and to confirm the performance and safety of the fuel system under prototypic conditions. This approach will reduce the time for development and qualification of a new fuel system, and it will also reduce associated costs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Comparison Calculations of Subcritical Limits in ANSI/ANS-8.1 [Slides]

The calculations presented in this paper were performed with ENDF/B-VIII.0 and SCALE 6.3.1 to provide users with perspective on what the current k eff values are for the SCLs listed in ANS-8.1. There are multiple SCLs that exceed the reference value of 0.98. This may indicate that these values need to be adjusted by the WG for the next iteration of the standard.

ENDF↗
DOE OSTI2025

Copper-dependent halogenase catalyses unactivated C−H bond functionalization

Carbon–hydrogen (C–H) bonds are the foundation of essentially every organic molecule, making them an ideal place to do chemical synthesis. The key challenge is achieving selectivity for one particular C(sp 3 )−H bond. In recent years, metalloenzymes have been found to perform C(sp 3 )−H bond functionalization. Despite substantial progresses in the past two decades, enzymatic halogenation and pseudohalogenation of unactivated C(sp 3 )−H—providing a functional handle for further modification—have been achieved with only non-haem iron/α-ketoglutarate-dependent halogenases, and are therefore limited by the chemistry possible with these enzymes. Here, in this work, we report the discovery and characterization of a previously unknown halogenase ApnU, part of a protein family containing domain of unknown function 3328 (DUF3328). ApnU uses copper in its active site to catalyse iterative chlorinations on multiple unactivated C(sp 3 )−H bonds. By taking advantage of the softer copper centre, we demonstrate that ApnU can catalyse unprecedented enzymatic C(sp 3 )−H bond functionalization such as iodination and thiocyanation. Using biochemical characterization and proteomics analysis, we identified the functional oligomeric state of ApnU as a covalently linked homodimer, which contains three essential pairs—one interchain and two intrachain—of disulfide bonds. The metal-coordination active site in ApnU consists of binuclear type II copper centres, as revealed by electron paramagnetic resonance spectroscopy. This discovery expands the enzymatic capability of C(sp 3 )−H halogenases and provides a foundational understanding of this family of binuclear copper-dependent oxidative enzymes.

biocatalysis↗
DOE OSTI2023

High-pressure and temperature neural network reactive force field for energetic materials

Reactive force fields for molecular dynamics have enabled a wide range of studies in numerous material classes. These force fields are computationally inexpensive compared with electronic structure calculations and allow for simulations of millions of atoms. However, the accuracy of traditional force fields is limited by their functional forms, preventing continual refinement and improvement. Therefore, we develop a neural network-based reactive interatomic potential for the prediction of the mechanical, thermal, and chemical responses of energetic materials at extreme conditions. The training set is expanded in an automatic iterative approach and consists of various CHNO materials and their reactions under ambient and shock-loading conditions. Further, this new potential shows improved accuracy over the current state-of-the-art force fields for a wide range of properties such as detonation performance, decomposition product formation, and vibrational spectra under ambient and shock-loading conditions.

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

A Field Guide to Corralling the Chaos: A Conceptual Framework for Using Models to Guide Opportunistic Field Studies of Natural Disturbances

Watersheds regulate biogeochemical processes and provide ecosystem services to human societies, but disturbances can fundamentally alter these processes across space and time. Determining when and where to sample to capture disturbance impacts in watersheds remains a central challenge. Manipulation studies and long-term monitoring are often constrained by scope, and opportunistic studies often lack pre-disturbance data needed to statistically determine disturbance impacts. We identify a persistent knowledge gap: the absence of a clear, transferable framework to guide opportunistic disturbance research where pre-disturbance data collection is not a feasible option. To address this gap, we present a conceptual framework that intentionally integrates modeling and empirical observation in an iterative, stepwise model–experiment workflow. We demonstrate its application through two contrasting case studies: wildfire impacts on headwater streams using a pre-disturbance preparedness approach, and saltwater flooding impacts on coastal forests using an ‘ex-post-facto’ approach. From these applications, we assess strengths, limitations, and the critical role of team science for transferability across disturbance types and study designs. Broadly, this framework offers a scalable path towards more rigorous, timely, and actionable disturbance science that can inform watershed management, hazard risk reduction, and ecosystem resilience.

Coastal Biogeochemistry↗
DOE OSTI2024

Quantum subspace expansion in the presence of hardware noise

Finding ground state energies on current quantum processing units (QPUs) using algorithms such as the variational quantum eigensolver (VQE) continues to pose challenges. Hardware noise severely affects both the expressivity and trainability of parameterized quantum circuits, limiting them to shallow depths in practice. Here, we demonstrate that both issues can be addressed by synergistically integrating VQE with a quantum subspace expansion, allowing for an optimal balance between quantum and classical computing capabilities and costs. We perform a systematic benchmark analysis of the iterative quantum-assisted eigensolver in the presence of hardware noise. We determine ground state energies of 1D and 2D mixed-field Ising spin models on noisy simulators and the IBM QPUs ibmq_quito (5 qubits) and ibmq_guadalupe (16 qubits). To maximize accuracy, we propose a suitable criterion to select the subspace basis vectors according to the trace of the noisy overlap matrix. Finally, we show how to systematically approach the exact solution by performing controlled quantum error mitigation based on probabilistic error reduction on the noisy backend fake_guadalupe.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Tula: Optimizing Time, Cost, and Generalization in Distributed Large-Batch Training

Distributed training increases the number of batches processed per iteration either by scaling-out (adding more nodes) or scaling-up (increasing the batch-size). However, the largest configuration does not necessarily yield the best performance. Horizontal scaling introduces additional communication overhead, while vertical scaling is constrained by computation cost and device memory limits. Thus, simply increasing the batch-size leads to diminishing returns: training time and cost decrease initially but eventually plateaus, creating a knee-point in the time/cost vs. batch-size pareto curve. The optimal batch-size therefore depends on the underlying model, data and available compute resources. Large batches also suffer from worse model quality due to the well-known “generalization gap”. In this paper, we present Tula, an online service that automatically optimizes time, cost, and convergence quality for large-batch training of convolutional models. It combines parallel-systems modeling with statistical performance prediction to identify the optimal batchsize. Tula predicts training time and cost within 7.5−14% error across multiple models, and achieves up to 20× overall speedup and improves test accuracy by ≈9% on average over standard large-batch training on various vision tasks, thus successfully mitigating the generalization gap and accelerating training at the same time.

Tyagi, Sahil [ORNL] (ORCID:0009000783144745)↗
DOE OSTI2025

Lightweight Metal Stamping Optimization Enabled by Artificial Intelligence

Successfully manufacturing an automotive body structure made via the sheet metal stamping process depends upon simultaneous consideration of component design, tooling design, stamping process control, and material properties. In many cases, introducing lightweight sheet materials (e.g., aluminum alloys, magnesium alloys, advanced high strength steels) holds the potential to significantly reduce vehicle weight, but challenges the stamping process by introducing materials with inherently less ductility. Successful and repeatable applications require co-developing the stamping process controls with the varying material properties, including formability. During the stamping process, as soon as the forming limit of the sheet is exceeded, the material shows localized necking which quickly leads to splits. Controlling process variability to avoid these material splits will enable deployment of less formable, lighter, and stronger materials for stamped automotive components. A typical optimization procedure for manufacturing requires an iterative process involving parameter setting, execution of computational simulations, and modifying the parameters. The entire process demands substantial computational time, making it impractical for real-time feedback towards rapid corrective actions required for in-line control for running production processes. To overcome this challenge, artificial intelligence (AI) can be leveraged to determine optimal manufacturing parameters within a single manufacturing cycle time. This research proposes an in-line optimization framework incorporating a trained AI model to predict kidney-shaped die forming. Preliminary results indicate that the AI framework can accurately predict draw-in values based on a given parameter set, a process referred to as forward prediction. Furthermore, the AI framework can also predict the optimal parameter set that leads to the desired draw-in values, referred to as inverse optimization (or backward prediction). This research has been performed in collaborations with USCAR (US Council for Automotive Research) and AutoForm. The members of USCAR are Ford, GM, and Stellantis.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Achievement of a high-density, high-confinement, and high-beta tokamak plasma regime in DIII-D, and implications for a lower-current path for ITER and FPP

Experiments on DIII-D have demonstrated a density-confinement synergy that enables sustainment of high performance in a previously unattained parameter regime of simultaneous very high energy confinement quality (H 98y2 ≥ 1.5), very high line-average density Greenwald fraction (ƒ Gr = πa 2 < n >/I P ≥ 1.4), and high toroidal beta (β T ≥ 3%). Tokamak operation in this regime is essential for a compact steady-state FPP, as well as for Q=10 with 500 MW of fusion power in ITER at I P << 15 MA. These experiments leveraged the knowledge that, in the high-poloidal-beta (β P ) regime, impurity and density gradients can enhance turbulence stabilization caused by high α MHD (α MHD ~(dβ P )⁄dr). This was described by theoretical predictions and gyrokinetic transport simulations [M.T. Kotschenreuther et al, 2024 Nucl. Fusion, 64 076033], and later confirmed by experiments on DIII-D [S. Ding et al, 2024 Nature 629 555]. To increase both β P and β T , the new experiments increased the ideal-wall stability β N -limit by using a smaller plasma-outer wall distance and higher triangularity in the plasma cross section (top/bottom average δ~0.9), enabled by the recent “shape & volume rise” (SVR) modification to the DIII-D divertor. The higher triangularity also contributed to achieving higher ƒ Gr by enabling higher pedestal density. At high density, the pedestal is ballooning limited and exhibits small and frequent ELMs, while the divertor is near detachment even without any impurity seeding. High plasma performance was attained and sustained reproducibly, with the eventual terminations brought about by an MHD mode destabilized as the current profile slowly continued to evolve. A path to stationary fully noninductive operation might include ECH injection to reduce both core impurity accumulation and the electron collisionality, thus increasing the bootstrap current. These experiments provide the first experimental demonstration of the ƒ Gr , H 98y2 , and β T values required simultaneously for ITER Q = 10 at I P < 10 MA, pointing to practical ways to improve the energy confinement in a fusion reactor.

Garofalo, Andrea M. [General Atomics, San Diego, C↗
DOE OSTI2025

Iterative Reconstruction for Multimodal Neutron Tomography

Here, we describe a unified framework for model-based iterative 3-D reconstruction of multimodal neutron transmission, hydrogen-scatter, and induced-fission images from low resolution data recorded using 14.1-MeV neutrons and the associated-particle imaging (API) technique. The framework, which was developed to facilitate use in challenging field-deployment scenarios, is centered around physics-based system models and a total variation (TV) constrained implementation of the simultaneous iterative reconstruction technique (SIRT). Modified to solve a statistically weighted least squares (WLS) problem, the SIRT algorithm is accelerated using ordered subsets and Nesterov’s momentum for which we derive a near-optimal value of the governing Lipschitz constant. The approach enables the reconstruction of images that are high resolution compared to the acquired data and is robust to both limited statistics and a limited number of projection angles. Moreover, the framework is fast enough to be practical. Example images are provided that demonstrate both the ability to perform fast-neutron imaging of high-atomic-number materials with low radiation dose and the benefit of multimodal neutron imaging to identify key materials.

Hydrogen scatter↗
DOE OSTI2024

Full time-dependent SOLPS-ITER simulation of the SPARC tokamak: actuator design for particle and divertor condition control *

Abstract This paper presents the application of full time-dependent SOLPS-ITER simulations for actuator design in the SPARC tokamak. This study employs both the EIRENE module, a neutral solver, and the B2.5 plasma module in a time-dependent mode. This is in contrast to most SOLPS simulations, which focus on steady-state solutions, where the neutral distribution is evolved without any time limit or for a time step of 1 ⋅ 10 − 3 second, which is several orders of magnitude larger than the fluid plasma time step. The time-dependent EIRENE was tested with a fixed B2.5 background and compared with a simple conductance based model in a simplified pump chamber geometry. This comparison aimed to verify the reliability of the neutral relaxation timescale derived from the time-dependent EIRENE. Subsequently, a full time-dependent simulation was performed in a realistic geometry, with the Monte-Carlo neutral time step synchronized with the plasma fluid time step. The numerical setup of the code, including relative time steps and the size of the census data used to store Monte-Carlo particle information is considered. The full-time dependent simulations are then applied to inform the design of the SPARC louver structure, which affects divertor plasma parameters by regulating the neutral conductance from the divertor to the pump. The response of the plasma and neutral parameters was captured on a timescale that enables the design of the actuator to consider time-dependent control capability. It was found that changing the louver opacity has an equivalent effect as varying the gas throughput via puff actuation. Therefore, equivalent divertor plasma conditions can be obtained from both actuators, while the neutral pressure distribution in the pump and divertor differs for each actuator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Overview of the ASDEX Upgrade results

After a 26-month vent ASDEX Upgrade (AUG) went back in operation with a newly designed upper W-divertor suitable for alternative divertor configurations (featuring in-vessel coils and cryo-pump). Parameter scans and an extensive set of measurements were obtained and their interpretation is ongoing. Prompted by the ITER wall change, dedicated experiments on non-boronized plasma startup were contrasted to that employing asymmetric and more symmetric boronizations. The asymmetric boronization proved to be as beneficial as the more symmetric one, which is in contrast to previous model calculations assuming perfect sticking of boron (measurements suggest sticking ≈ 0.3). In the startup phase also the impurity influxes at the outboard limiters were investigated contrasting the unboronized case featuring cold edges (low- Z radiation) to the boronized case, in which the lifetime of the boron layers could be estimated. Pedestal stability investigations revealed that the quasi-continuous exhaust (QCE) regime is obtained when ballooning modes are active in the vicinity of the separatrix and the global peeling-ballooning stability is high enough. Thus, at high enough shaping and high gas flux both can be achieved and QCE is a consequence. The closely related enhanced D-alpha (EDA) mode is not clearly distinguishable from QCE, e.g. the quasi-coherent mode characteristic for EDA also shows up in QCE. In QCE the impurity transport is behaving benign as could be measured for Ne with a novel analysis method making use of a comprehensive set of CXRS measurements. For high radiative fractions the regime of the X-point radiator (XPR) is accessible at AUG and the understanding of its access conditions and behaviour is further developed. Due to the localized radiative cooling at the X-point the XPR can be well diagnosed and thus controlled. For negative triangularity shapes, further experiments at increased shaping resulted in strongly heated L-mode plasmas avoiding ELMs. Two integrated modelling approaches towards ITER suggested that core W-accumulation will be no issue for ITER and that the fusion yield in ITER may be Q = 12 (i.e. ITPA20-IL scaling is too pessimistic). Further, investigations of the ITER ramp-down in AUG provide insights into maintaining position control. Various aspects of shattered pellet injection were investigated in AUG and one of the results show that with increasing Ne fraction the radiation during the current quench increases and the current decay becomes faster.

ASDEX Upgrade↗
DOE OSTI2023

Learning local equivariant representations for large-scale atomistic dynamics

Abstract A simultaneously accurate and computationally efficient parametrization of the potential energy surface of molecules and materials is a long-standing goal in the natural sciences. While atom-centered message passing neural networks (MPNNs) have shown remarkable accuracy, their information propagation has limited the accessible length-scales. Local methods, conversely, scale to large simulations but have suffered from inferior accuracy. This work introduces Allegro, a strictly local equivariant deep neural network interatomic potential architecture that simultaneously exhibits excellent accuracy and scalability. Allegro represents a many-body potential using iterated tensor products of learned equivariant representations without atom-centered message passing. Allegro obtains improvements over state-of-the-art methods on QM9 and revMD17. A single tensor product layer outperforms existing deep MPNNs and transformers on QM9. Furthermore, Allegro displays remarkable generalization to out-of-distribution data. Molecular simulations using Allegro recover structural and kinetic properties of an amorphous electrolyte in excellent agreement with ab-initio simulations. Finally, we demonstrate parallelization with a simulation of 100 million atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Applying Constrained Bayesian Optimization to the Design of Critical Experiments

Often when planning a criticality experiment, many design configurations are iteratively investigated with a Monte Carlo transport code. The goal is that the experiment will be optimal with respect to some variable, like the fraction of fissions occurring at a certain energy range, while simultaneously being critical. Unfortunately, the Monte Carlo transport simulations are expensive, which can ultimately limit the number of configurations that can be explored. In this work, we present how Gaussian processes (GPs) can be used as a reduced-order model in a constrained Bayesian optimization (CBO) algorithm to design a criticality experiment. The GPs replace the Monte Carlo transport simulations that explore the design space. The CBO algorithm efficiently identifies new points in the design space to run the Monte Carlo transport code while respecting the criticality constraint. It does so in a manner that both improves the accuracy of the GP and finds the approximate global optimum. We demonstrate the performance of CBO with the design of a Thermal Epithermal eXperiment (TEX) for the criticality safety validation of nuclear waste models of the Hanford Tank Farm.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

SYS-5620-Project - Achieving Robust Laser Performance utilizing Historical Shot Experiments - A Systems Engineering Proposal

The National Ignition Facility (NIF) at Lawrence Livermore National Laboratory precisely guides, amplifies, reflects, and focuses 192 powerful laser beams into a target about the size of a pencil eraser in a few billionths of a second, delivering more than 2 million joules of ultraviolet energy and 500 trillion watts of peak power. A crucial goal of the system is to trigger precise implosions of fuel capsules. This is achieved by delivering all 192 beams at user-specified times and locations on the target, minimizing any deviation from the requested performance. Power requirements can vary substantially on each experiment and the facility supports numerous amplifier pumping configurations and their attendant nonlinear effects. To achieve the tight performance required across such a broad array of configurations, constant comparison of measured and requested power delivery are tracked and long-term trends analyzed as a guide to understanding future performance. A very common question asked is given the current state of the laser today – how well would a similar experiment from the past perform today? Additionally, if one were to specify an alternate amplifier configuration using today’s model would and damage limits be exceeded and would there be an increase in performance. The physics model used to make these predictions and equipment protection checks is called the Virtual Beamline (VBL). VBL is used with an incoming desired pulse shape and energy to be delivered on target, and then does an iterative solve to predict the needed injected pulse in the front-end of the system to achieve this result. To effectively guide and predict future NIF experiment performance, laser scientists explore current amplifier configurations and compare them with historical data utilizing a tool called the Reverify Toolbox.

42 ENGINEERING↗