Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical integration”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A hereditary integral, transient network approach to modeling permanent set and viscoelastic response in polymers

An efficient numerical framework is presented for modeling viscoelasticity and permanent set of polymers. It is based on the hereditary integral form of transient network theory, in which polymer chains belong to distinct networks each with different natural equilibrium states. Chains continually detach from previously formed networks and reattach to new networks in a state of zero stress. The free energy of these networks is given in terms of the deformation gradient relative to the configuration at which the network was born . A decomposition of the kernel for various free energies allows for a recurrence relationship to be established, bypassing the need to integrate over all time history. The technique is established for both highly compressible and nearly incompressible materials through the use of neo-Hookean, Blatz–Ko, Yeoh, and Ogden-Hill material models. Multiple examples are presented showing the ability to handle rate-dependent response and residual strains under complex loading histories.

Finite element method↗
DOE OSTI2026

WellPINN: Accurate Well Representation for Transient Fluid Pressure Diffusion in Subsurface Reservoirs With Physics‐Informed Neural Networks

Accurate representation of pumping wells is essential for reliable reservoir characterization and simulation of operational scenarios in subsurface flow models. Physics-informed neural networks (PINNs) are emerging as a promising alternative to numerical models for reservoir modeling, offering seamless integration of monitoring data and governing physical equations. However, existing PINN-based studies face major challenges in capturing fluid pressure near wells when using a source/sink term, particularly during the early stages after pumping begins. We address this problem by introducing WellPINN, a workflow in which an initially trained PINN infers fluid pressure across the entire reservoir domain using a large equivalent well radius. This initial PINN solution is then locally refined around the well by a set of subdomain PINNs that are trained for smaller equivalent well radii. Continuity across these subdomain interfaces as well as at the initial condition is ensured by hard-constraining each PINN on its subdomain boundary. Our results demonstrate WellPINN as the first workflow of its kind to focus on accurate inference of fluid pressure from pumping rates throughout the entire injection period, significantly advancing the potential of PINNs for inverse modeling and operational scenario simulations. All data and code for this paper are openly available at https://doi.org/10.20350/DIGITALCSIC/17260.

58 GEOSCIENCES↗
DOE OSTI2025

Flame Characterization of a Hydrogen Micromix Burner using Experimental and Numerical Methods - Extended Abstract

Micromix combustion presents a promising solution for integrating hydrogen into gas turbine combustors, offering rapid and compact air-fuel mixing that effectively shortens the combustion zone and significantly reduces nitrogen oxides (NOX) emissions. A major challenge in hydrogen combustion systems is flashback, which is mitigated through short mixing lengths and high flow velocities. This study investigates three fuel injector designs – zero swirl, low swirl, and high swirl – for a 6-mm inner diameter nozzle, where hydrogen is injected through six 0.3 mm diameter holes in a jet-in-crossflow configuration relative to the airstream. Numerical simulations indicate stable flame formation, moderate peak flame temperatures, and subsequent low NOx potential for all designs. These injectors have been fabricated using additive manufacturing, and upcoming experiments will explore flame structure and NOx behavior across varying fuel injection angles, equivalence ratios and pre-mixer velocities. These efforts aim to advance the development of low-emission hydrogen micromix combustion.

hydrogen↗
DOE OSTI2024

Surfaceology for colored Yukawa theory

Arkani-Hamed and collaborators have recently shown that scattering amplitudes for colored theories can be expressed as integrals over combinatorial objects simply constructed from surfaces decorated by kinematic data. In this paper we extend the curve integral formalism to theories with colored fermionic matter and present a compact formula for the all-loop, all-genus, all-multiplicity amplitude integrand of a colored Yukawa theory. The curve integral formalism makes certain properties of the amplitudes manifest and repackages non-trivial numerators into a single combinatorial object. We also present an efficient formula for L-loop integrated amplitudes in terms of a sum over 2 L combinatorial determinants.

1/N Expansion↗
DOE OSTI2025

Using "AI Poincare" to analyze non-linear integrable optics

This study dives into the applicability of using automated discovery of conserved quantities in dynamical systems relevant to accelerator physics. Specifically, we explore the performance of AI Poincaré in analyzing numerical trajectory data obtained using the McMillan system of non-linear integrable optics. A comprehensive evaluation of the algorithm's performance is conducted through diverse methodologies. These include the analysis of the estimated number of conserved quantities embedded in a dataset and the deviation of interpolated points on the inferred manifold with respect to points in actually in the dataset. the investigation identifies an optimal range of perturbation distances where the underlying manifold extraction algorithm inside AI Poincaré exhibits optimal performance. Additionally, an improved neural network architecture is proposed based on the observed results. Finally, we apply the algorithm to preliminary experimental data from the Integrable Optics Test Accelerator at Fermilab to successfully infer the number of conserved quantities even in the presence of fast decoherence of the measured signal.

Osmanov, Lazare [Free U. Tbilisi]↗
DOE OSTI2025

On nonlocal problems with Neumann boundary conditions: scaling and convergence for nonlocal operators and solutions

Formulations of Neumann-type boundary conditions for boundary value problems in the nonlocal framework are beset with difficulties, some related to the choice of a proper scaling. Here we identify a space-dependent scaling for a nonlocal Neumann operator, for which we prove linear in δ (δ being the radius for the support for the kernel) convergence of the Neumann operator and $\mathcal{O}$(δ 2 ) convergence of solutions to their classical counterparts. The pointwise-like convergence of the nonlocal normal operator is cast as a new type of two-scale operator-point convergence, which we call condensated convergence . The results hold for general integrable kernels, a setting which is favored in numerical simulations. We support this analysis with numerical convergence studies using a piecewise linear discontinuous Galerkin discretization and show an $\mathcal{O}$(δ 2 ) rate of convergence of solutions, also exhibiting an $\mathcal{O}$(h 2 ) convergence, where h is the mesh size.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

CSAC observations from 2025 CHEDAR review [Slides]

Overall, the CSAC is impressed by the strong start of CHEDAR, a multi-institutional center with broader scope than its predecessor. The CSAC was impressed and energized by discussions with center students and staff; the enthusiasm for science is universal and the sense of the contribution in context to the NNSA mission is mostly understood. The Center is well positioned to build on a strong start and (1) enhance integration of code-development targeted to, and numerical design in on-going experimental campaigns, (2) focus and better tell the story of the astrophysical connection for on-going experimental campaigns, and (3) expand and deepen ties between Center students (and faculty) and the NNSA National Laboratories. Overall, the CSAC was encouraged and energized by the materials presented. Now, the Center has to recognize the need to continue growing strategically. There is abundant access to experimental facilities – is that a challenge to focused effort? The PIs will need to balance the pace of experimental campaigns and analysis and simulation so that informed progress and focused investigation can develop in each campaign. There was a clear need to connect simulations to the goals of experimental campaigns – what is the best way to cultivate that connection? Are there adequate resources in the near term to develop code capability? What is the balance between close personal mentorship and collaboration across the separate institutions that compose the Center? The CSAC feels confident that the team leading CHEDAR will find creative solution to these challenges.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Considerations for a Co-Located Aquaculture and Wave Energy Deployment: Quantifying Demand and Assessing Integration

The co-location of marine energy and aquaculture is a concept of increasing value and interest in the United States, as the desire for sustainably produced seafood and renewable energy continues to grow. With both industries being fairly nascent in their national development and facing challenges, the deployment of a wave-powered aquaculture farm requires numerous considerations. This report analyzes the energy demands of an integrated multi-trophic aquaculture system at various farm scales and discusses the preliminary technical, regulatory, and logistical considerations required for a co-located deployment. A pilot-scale deployment is outlined with potentially ideal characteristics, and an optimized planning and permitting approach is introduced.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Nyström type exponential integrators for strongly magnetized charged particle dynamics

Solving for charged particle motion in electromagnetic fields (i.e. the particle pushing problem) is a computationally intensive component of particle-in-cell (PIC) methods for plasma physics simulations. This task is especially challenging when the plasma is strongly magnetized due numerical stiffness arising from the wide range of time scales between highly oscillatory gyromotion and long term macroscopic behavior. A promising approach to solve these problems is by a class of methods known as exponential integrators that can solve linear problems exactly and are A-stable. This work extends the standard exponential integration framework to derive Nyström-type exponential integrators that integrates the Newtonian equations of motion as a second-order differential equation directly. In particular, we derive second-order and third-order Nyström-type exponential integrators for strongly magnetized particle pushing problems. Numerical experiments show that the Nyström-type exponential integrators exhibit significant improvement in computation speed over the standard exponential integrators.

general physics↗
DOE OSTI2025

RE-INTEGRATE EMT Simulation Software: DAE Solvers and Automation

Existing electromagnetic transient (EMT) simulation tools face challenges in accelerating EMT simulations, especially for very large-scale power networks. To tackle this issue, next generation EMT simulation tools such as RE-INTEGRATE EMT are being researched upon. Such tools should be equipped with automation capabilities and advanced numerical differential-algebraic equation (DAE) solvers. In this paper, the DAE solvers incorporated within the RE-INTEGRATE EMT simulation tool are discussed. In particular, a modified ODEINT-based DAE solver and the ARKODE solver from SUN-DIALS are leveraged within RE-INTEGRATE EMT. In addition, the automation implemented within RE-INTEGRATE EMT to automate the DAE generation (replacing the need of manual discretization and assembling DAEs) is discussed. Different use cases were implemented using the RE-INTEGRATE EMT tool and were validated with respect to baseline simulations.

Marthi, Phani Ratna Vanamali [ORNL] (ORCID:0000000↗
DOE OSTI2025

Stress accommodation in nanoscale dolan bridges designed for superconducting qubits

Josephson junctions are the principal circuit element in numerous superconducting quantum information devices and can be readily integrated into large-scale electronics. However, device integration at the wafer scale necessarily depends on having a reliable, high-fidelity, and high-yield fabrication method for creating Josephson junctions. When creating Al/AlO x based superconducting qubits, the standard Josephson junction fabrication method relies on a sub-micron suspended resist bridge, known as a Dolan bridge, which tends to be particularly fragile and can often times fracture during the resist development process, ultimately resulting in device failure. In this work, we demonstrate a unique Josephson junction lithography mask design that incorporates stress-relief channels. Our simulation results show that the addition of stress-relief channels reduces the lateral stress in the Dolan bridge by more than 70% for all the bridge geometries investigated. In practice, our novel mask design significantly increased the survivability of the bridge during device processing, resulting in 100% yield for over 100 Josephson junctions fabricated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $\mathscr{O}$($v$/$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $\mathscr{O}$($v$/$c$) Eulerian-frame number density and is consistent, to $\mathscr{O}$($v$/$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Integrated photonic structures for photon-mediated entanglement of trapped ions

Trapped atomic ions are natural candidates for quantum information processing and have the potential to realize or improve quantum computing, sensing, and networking. These applications often require the collection of individual photons emitted from ions into guided optical modes, in some cases for the production of entanglement between separated ions. Proof-of-principle demonstrations of such photon collection from trapped ions have been performed using high-numerical-aperture lenses or cavities and single-mode fibers, but integrated photonic elements in ion-trap structures offer advantages in scalability and manufacturability over traditional optics. In this paper we analyze structures monolithically fabricated with an ion trap for collecting ion-emitted photons, coupling them into waveguides, and manipulating them via interference. We calculate geometric limitations on collection efficiency for this scheme, simulate a single-layer grating that shows performance comparable to demonstrated free-space optics, and discuss practical fabrication and fidelity considerations. Based on this analysis, we conclude that integrated photonics can support scalable systems of trapped ions that can distribute quantum information via photon-mediated entanglement.

Knollmann, Felix W.↗
DOE OSTI2026

Hybrid Symbolic-Numerical Modeling and Parametric Stability Analysis of DC–AC Power Systems

Hybrid DC-AC power systems integrating diverse inverter-based resources (IBRs) and multi-terminal high-voltage direct current (MTDC) networks represent a promising architecture for future power grids, while introducing challenges for modeling, stability analysis, and control design. This paper develops a hybrid symbolic-numerical modeling framework and tool to characterize the parametric small-signal stability of DC-AC coupled power systems. The proposed approach constructs parametric state-space models to enable efficient representation of system dynamics under varying control parameters and network configurations, with target parameters retained as symbolic variables and the remainder treated numerically. The stability analysis framework covers eigenvalue, sensitivity, and stability boundary and region characterization. Enhanced linear matrix inequality (LMI) techniques are proposed to directly certify small-signal stability over regions of parameter space while also reducing the conservativeness and computational burden. The resulting tools and frameworks enable rapid parametric model construction across diverse grid conditions, thereby facilitating stability-informed control and operation in complex DC–AC power systems.

DC–AC power systems↗
DOE OSTI2025

Boosting efficiency and reducing graph reliance: Basis adaptation integration in Bayesian multi-fidelity networks

The computational cost of high-fidelity numerical models makes outer-loop analysis, which requires repeated interrogation of the model such as uncertainty quantification, computationally demanding. Multi-fidelity methods, which construct a surrogate model using data from an ensemble of models of varying cost and accuracy, can substantially reduce the cost of outer-loop analysis. However, these methods can be difficult to apply when the model ensemble does not admit a clear hierarchy a priori and the correlations between models are low. Consequently, in this paper, we present a multi-fidelity method that leverages dimension reduction to enhance the correlation between models, thereby reducing the amount of data needed to train a surrogate from an unordered ensemble of models. Our method utilizes basis adaptation to build low-dimensional polynomial chaos expansions of each model and employs Multi-fidelity Networks to encode the relationships among models. We show that the resulting method exhibit two notable advantages over its counterpart: (1) enhanced accuracy (both reduced bias and variance); and (2) reduced dependency on the graph structure encoding relationships among models. We demonstrate the approach on an analytical test problem and a challenging finite element model for a spent nuclear fuel. Our method produces a surrogate model that is significantly more accurate than either a single-fidelity surrogate or a multi-fidelity surrogate constructed without basis adaptation.

42 ENGINEERING↗
DOE OSTI2026

Real-time plasma monitoring framework for advanced plasma control and ML-research in DIII-D

Real-time and adaptive plasma control is crucial for robust tokamak operation, requiring sensitivity and tolerance measurements of the plasma state. This paper presents the implementation of an integrated real-time plasma monitoring framework on the DIII-D tokamak to support advanced control approaches, including machine-learning (ML) methods. The system is built on the SHIELD framework, a high-performance modular architecture that provides a unified pipeline for integrating diverse diagnostics. The framework leverages high-bandwidth digitizers, fast numerical processing, and deterministic, low-latency interconnects to stream high-fidelity data from diagnostics such as electron cyclotron emission (ECE), beam emission spectroscopy (BES), CO interferometers, and a visible tangential divertor camera (TangTV). The system’s validity is demonstrated through direct comparisons of real-time and offline data. Furthermore, we present two key applications of the developed plasma monitoring system with ML-based plasma control strategies, including real-time divertor detachment and active Alfvén Eigenmode control. As a result, this work presents a robust and scalable approach for integrating high-frequency, multidimensional diagnostics into advanced control algorithms for future fusion devices.

AI/ML↗
DOE OSTI2026

Experimental characterization and analysis of phase change material-based thermal energy storage system for refrigerated display case

Refrigerated display cases that are used to store and exhibit food products in supermarkets and retail spaces consume a significant portion of the buildings’ total electricity. Importantly, the refrigeration-related energy cost and demand charges are greatly affected by the time-of-use electricity pricing and demand rates, which are at their maximum during peak hours, typically when refrigeration energy consumption is also high. Using energy storage to shift the refrigeration load from peak to off-peak hours can greatly reduce the operational costs in the supermarket. This study demonstrates a phase change material-based thermal energy storage (TES) system, specifically designed in stackable units, that can be integrated with an open vertical refrigerated display case. We perform numerical and experimental characterization that includes finite-difference modeling for the TES, prototype fabrication, and laboratory evaluation, followed by a preliminary system-level analysis to predict the impact of TES on the refrigerated case performance, energy use, and energy cost. The results show that the dedicated latent TES for refrigerated cases can provide a specific energy of 50.4 Wh/kg and a specific power of 15.5 W/kg. The TES can be charged during 12 h of the off-peak period and discharged at various rates during 4 to 6 h of the peak period, thereby shifting the refrigeration load from the peak to the off-peak period. Consequently, annual cost savings up to 19% can be achieved, depending on the thermal load, the summer/winter peak electricity pricing, and the transition temperature of the phase change material used.

25 ENERGY STORAGE↗
DOE OSTI2026

On the roles of welding residual stresses in determination of fracture toughness in austenitic stainless steel SUS 304 pipeline girth welds

Welding residual stresses especially the high tensile stresses are proved to have negative impacts on the fatigue and fracture behaviors of welded structures. In this study, a virtual fabrication of test specimens from welding process to specimen preparation was carried out by numerical simulation. An austenitic stainless steel multi-pass pipe welding was simulated by transient thermal–mechanical finite element analysis, the residual stresses were then mapped into the test specimen to evaluate fracture toughness. The findings in this study confirmed that, residual stress can be high in a sub-sized compact tensile specimen, which may accelerate or hinder the crack propagation during actual fatigue and fracture tests as reported in recent years. The influence of the cutting location and orientation of the specimen on fracture performance was investigated systematically to provide a fundamental understanding of welding residual stress and necessary insights into the specimen preparation procedure. Considering the limitation of measuring techniques and the complexity of the stress distribution, the developed numerical model can be a very useful tool to elucidate the stress evolution and quantify the effect of remaining welding stress on fracture toughness.

Fracture behavior↗