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

Automated Detection of Qubit Structures on Quantum Chips

In this study, YOLO (You Only Look Once), a well-known object detection and image segmentation model, is used to detect qubits on a quantum chip. The model is trained by exposing it to various images of objects of interest and by tweaking various training parameters to allow the model to learn from a limited and highly specialized dataset of SEM imagery. The model was provided with numerous images of qubits and various prominent features of note located on a qubit. These detections were then integrated into a simple GUI, allowing for real-time feedback alongside SEM usage. This will allow us to create various functions enabling automated functions like chip-wide scanning, auto-focusing, and image correction.

Perjuste, Ruth [Wellesley Coll.]↗
DOE OSTI2020

Acquisition and Implementation of a Comprehensive Environmental Permits Linking Tool at Savannah River Site - 20234

Historically, Savannah River Nuclear Solutions LLC (SRNS) tracked environmental regulatory commitments and the requirements from hundreds of permits at the Savannah River Site (SRS) using several separate methods, making integrated compliance assurance cumbersome and labor-intensive. When SRNS experienced an increase in environmental issues in 2017, SRNS management and U.S. Department of Energy - Savannah River (DOE-SR) management decided a single, proactive approach was needed to capture environmental permit information (including regulations, Consent Orders, DOE Orders, and any other state or federally issued statement of requirements), track the tasks necessary to ensure compliance with these requirements, and thereby mitigate the risk of noncompliance. SRNS developed a list of mandatory objectives that the tool must meet to function as a Comprehensive Environmental Permits Linking Tool (CEPLT). A key requirement was the ability to map Site permits to their governed locations and display the associated requirements at the compliance point (e.g., outfall, stack, waste unit, etc.). Several options included modifying existing onsite resources, building a custom onsite solution, purchasing an off-the-shelf solution, and contracting an offsite developer to build a custom solution. SRNS concluded an off-the-shelf solution with configuration and customization options would provide the flexibility to fit the unique needs of Savannah River Site (SRS) while taking advantage of industry-tested software and providing a reduced deployment timeline. SRNS chose Gensuite{sup R} a, a cloud-based solution that offers numerous a la carte applications in the environmental, health, and safety arenas, as the best candidate. SRNS selected three (3) integrated applications (Compliance Calendar, Permit Manager, and Mapper) to function as the CEPLT. The Compliance Calendar module allows for creation and tracking of regulatory commitment tasks assigned to responsible environmental professionals. Permit Manager organizes permits and other requirement documents, linking the commitments in each to Compliance Calendar tasks and/or implementing procedures. Mapper provides GIS capability for mapping the data from the other two modules to their physical onsite locations. The CEPLT was configured to allow for other SRS Site Tenants to eventually utilize the applications. Additionally, DOE-SR uses the CEPLT to provide an overview of contractor environmental compliance activities and to organize DOE-specific documents and tasks. The CEPLT fulfilled the requirement of meeting current compliance needs as well as providing the ability to grow as new organizations are incorporated and functionality is expanded. SRNS now uses the CEPLT to more effectively manage permit requirements, improve knowledge transfer, and increase the Site's overall protection of the environment, the Site worker, and the public. Possible future uses of the system include integration of mobile applications for timely communication of potential non-compliant conditions and DOE complex-wide deployment allowing for enhanced DOE site and Head Quarters oversight. Implementation of CEPLT will result in cost savings, both in terms of dollars and man-hours. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Short-Term Probabilistic Solar Forecasting via Reinforcement Learning over ECMWF

In this paper, we present an innovative reinforcement learning approach for short-term solar forecasting, leveraging data from the European Centre for Medium-Range Weather Forecasts (ECMWF). The methodology begins with the application of the System Advisor Model (SAM) to transform various ECMWF numerical weather prediction members into predictive photovoltaic power generation. To enhance the precision of deterministic forecasting, we introduce a dynamic model selection algorithm based on Q-learning. This algorithm dynamically identifies and utilizes the most accurate ensemble member for forecasting purposes. Furthermore, we employ a support vector regression surrogate model with a Gaussian distribution to generate probabilistic forecasts, providing a holistic view of solar energy generation uncertainty. To expedite the training process and make it more practical for real-world applications, we integrate a rolling update workflow. This innovative workflow reduces the training period from months to a mere 19 days, making our method highly efficient. Numerical results of the case study show that in comparison to benchmark models, the proposed method improves the deterministic and probabilistic solar forecasting accuracy by up to 40.84% and 48.42%, respectively.

ensemble forecasting↗
DOE OSTI2020

A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing is developed using a general buoyancy–shear–drag model. Analytic solutions to the simplified model equations corresponding to each instability are derived separately, which are then used to calibrate the model coefficients to predict specific values of the mixing layer growth parameters and exponents. The buoyancy-shear–drag equations for the bubble and spike mixing layer widths are then solved numerically, and a turbulent diffusivity (or viscosity) is constructed dimensionally as the product of the mixing layer width, h, and its time-derivative, dh/dt. Surrogate turbulent fields are then constructed by multiplying a presumed, approximate self-similar spatial profile by appropriate functions of h and dh/dt. Using several simplifying approximations, the turbulent diffusion equations satisfied by the mean mass fraction and mean shear velocity are solved analytically in a reference frame moving with the mean advection velocity, and it is shown that the mean fields evolve in space and time as expected. The explicit modeling and solution of turbulent transport equations (e.g., for the turbulent kinetic energy and its dissipation rate or a lengthscale) are not required. By using separate equations for the bubble and spike evolution, the model naturally captures the asymmetry induced by an increasing density contrast between the heavy and light fluids. The model also includes molecular dissipation and diffusion, and therefore, can describe transition to fully-developed turbulence if the initial turbulent diffusivity/viscosity is much smaller than the molecular diffusivity/viscosity. This model is applied to constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and variable-density Kelvin–Helmholtz mixing layers to demonstrate its utility. It is shown that the numerical solutions of the model calibrated using specific values of the instability growth parameters and exponents: (1) produces mixing layer widths in agreement with the expected self-similar growth power laws; (2) gives spatiotemporal profiles of turbulent fields that are expected and consistent with previous results; (3) predicts the expected power-law growths and decays of the spatially-integrated turbulent fields, and; (4) gives spatiotemporal profiles of the mean fields that are expected and consistent with previous results. This relatively simple model can be implemented numerically, and is expected to be useful for fundamental as well as astrophysical and inertial confinement fusion applications.

42 ENGINEERING↗
DOE OSTI2021

Electromagnetic Transient Simulation Algorithms for Evaluation of Large-Scale Extreme Fast Charging Systems (Distribution Grid Models)

The distribution and transmission grids are observing an increased penetration of power electronics in loads and generations. For example, there is increasing interest in integrating in extreme fast charging (XFC) systems for fast charging of electrical vehicles. As these systems are integrated, developing high-fidelity electromagnetic transient model of XFC systems in distribution grids and evaluating their interactions with the power grid would be of significant interest. This model will be utilized for design of XFC systems, to identify upgrades in distribution and/or transmission grids, for planning purposes by transmission planners or operators or owners, among others. It can also be utilized in operations for improved reliable performance of the grid and/or XFC station. The challenge with simulating these models is the high computational complexity introduced by the large number of states present in the system and the time-step needed to simulate the system. In this paper, advanced simulations algorithms are applied to reduce the computational complexity of simulating large-scale XFC systems. The algorithms include numerical stiffness-based segregation, time constant-based segregation, clustering and aggregation on differential algebraic equations (DAEs), and multi-order integration approaches. While the first three algorithms split the matrix that needs to be inverted from a large matrix to much smaller matrices, the final algorithm reduces the computational burden of applying higher-order integration approaches in the complete system. The comparison made in the previous sentence is with respect to use of homogeneous integration approaches used in conventional electromagnetic transient simulators like power systems computer aided design (PSCAD). The approaches mentioned here have resulted in speed-up of 36x in the simulation of a single distribution system with 15 XFCs.

Debnath, Suman↗
DOE OSTI2025

Modular Hydronic Room Conditioning System (CRADA NFE-24-10120 Final Report)

This project presents the development and validation of high-fidelity numerical models for fin-tube heat exchangers to enable accurate performance prediction and informed design optimization. The modeling framework integrates detailed geometric specifications, thermophysical property data, and system-level constraints to simulate the heat exchanger’s behavior under a range of operating conditions. The model is calibrated using real-world product specifications and validated against experimental data collected from controlled cooling and heating tests. In cooling mode, the model captures the overall trends in capacity and outlet air temperature but tends to underpredict latent effects, especially at lower air flow rates. In heating mode, the simulation consistently overestimates both the thermal capacity and outlet air temperature, indicating the need for refinement in air-side heat transfer assumptions. Despite these deviations, the model provides a solid foundation for optimization, allowing key design variables to be tuned within physical and performance-based constraints. This research advances the ability to simulate, validate, and optimize fin-tube heat exchanger designs with greater confidence and efficiency.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Ultra-narrow-linewidth hybrid-integrated self-injection locked laser at 780 nm

Narrow-linewidth lasers are essential across a wide range of applications, including classical and quantum sensing, trapped ion systems, position/navigation/timing systems, optical clocks, and microwave frequency synthesizers. In the visible and near-visible spectrum, low-noise lasers are particularly important for laser trapping and cooling techniques, which are vital for trapped ion quantum computing, sensing, and atomic clocks. In this context, our work showcases a hybrid-integrated narrow-linewidth laser that operates at 780 nm, achieving a self-heterodyne linewidth of 105 Hz. To validate the experimental results, we performed a numerical analysis that combines insights from a many-body theory applied to the gain region with a travelling-wave model to capture the laser dynamics. Our investigation further delves into how the linewidth of the self-injection locked lasers is influenced by the parameters of micro-ring resonators, aiming to assess the potential for achieving Hz-level integrated laser linewidths at 780 nm. This work not only demonstrates the technical feasibility of Hz-level narrow-linewidth lasers but also lays the groundwork for future explorations in the field.

Prokoshin, Artem (ORCID:0009000699760238)↗
DOE OSTI2025

Time-domain all-frequency stable formulation for low-frequency electromagnetic simulation with Newmark-β time integration

An implicitly Coulomb-gauged A-ϕ formulation has previously been proposed and validated for finite ele- ment simulations of low-frequency and multiscale electromag- netic problems in the frequency domain. This formulation has demonstrated numerical stability across all frequencies, with its accuracy, efficiency, and iterative convergence established in various frequency-domain scenarios. However, direct time- domain computation is often preferable for wideband electro- magnetic problems and is typically indispensable in nonlinear and multiphysics simulations. In this work, the A-ϕ formulation is extended to the time domain. By incorporating the well-known Newmark-β time integration scheme, the proposed formulation is validated through capacitive and inductive test cases. The results confirm the solution’s accuracy and demonstrate the formulation’s stability in the time domain.

Mekonnen, Minyechil↗
DOE OSTI2024

Integration of 5G and Time Sensitive Networks in Fossil Energy Generation Systems: A Case Study

Precise timing and data transmission within stringent time constraints are critical for numerous applications, such as robotics, virtual and augmented reality, industrial automation, energy and medical business, and various other sectors. Time-sensitive networks (TSN) and fifth-generation wireless communications (5G) are crucial for industrial communications, enabling convergent communication for various services using a common network core. Applications that are time-sensitive and require deterministic communications with low latency fall into this category, such as generation plant systems. This paper presents a simulated model of 5G-TSN for a fossil power plant based on the wired network parameters implemented in situ. Metrics analysis, comparison and future work are presented. © 2024 IEEE.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Powering the Blue Economy: Progress Exploring Marine Renewable Energy Integration With Ocean Observations

The blue economy is a dynamic and rapidly growing movement that captures the interplay between economic, social, and ecological sustainability of the ocean and encompasses numerous maritime sectors and activities (e.g., commerce and trade; living resources; renewable energy; minerals, materials, and freshwater; and ocean health and data). The demand for ocean data to inform scientific, risk reduction, and national security needs is leading to a large increase in the number of deployed ocean observation and monitoring systems, most of which require increased power. Because ocean observation systems are often placed in remote locations, they primarily rely on energy storage (or in some cases in situ energy generation) to power instruments and equipment, which imposes limits on sampling rates, deployment times, and spatiotemporal resolution of data. The U.S. Department of Energy Water Power Technologies Office is exploring the potential for marine renewable energy (MRE) devices (largely wave and tidal energy converters) to provide power to support multiple blue economy opportunities. A portion of these opportunities focus on power at sea markets for providing power in off-grid and offshore locations to support a variety of ocean-based activities, including ocean observation and navigation, underwater vehicle charging, marine aquaculture, marine algae farming, and seawater mining. Initially, research has focused on better understanding how and where MRE can provide a consistent source of reliable power to extend ocean observing missions, including operation of autonomous underwater vehicles. Online surveys as well as phone and inperson interviews were conducted with experts in the field of ocean observing systems and observatories to gather end-user requirements, determine energy needs, identify opportunities for codevelopment, and pinpoint constraints for MRE to meet those needs. The surveys and interviews provided feedback on the potential for powering devices and vehicles using MRE, including identifying common themes and challenges that will inform foundational research and development steps needed to advance the integration of MRE with ocean observing systems. In most cases, additional power generation on the order of watts was identified as significantly beneficial to enhancing ocean observations capabilities.

marine renewable energy, powering the blue economy↗
DOE OSTI2023

The reflected entanglement spectrum for free fermions

We consider the reflected entropy and the associated entanglement spectrum for free fermions reduced to two intervals in 1 + 1 dimensions. Working directly in the continuum theory the reflected entropy can be extracted from the spectrum of a singular integral equation whose kernel is determined by the known free fermion modular evolved correlation function. We find the spectrum numerically and analytically in certain limits. For intervals that almost touch the reflected entanglement spectrum approaches the spectrum of the thermal density matrix. This suggests that the reflected entanglement spectrum is well suited to the task of extracting physical data of the theory directly from the ground state wave function.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Contributions of biofilm-induced flow heterogeneities to solute retention and anomalous transport features in porous media

We report that microbial biofilms are ubiquitous within porous media and the dynamics of their growth influence surface and subsurface flow patterns which impacts the physical properties of porous media and large-scale transport of solutes. Here, a two-dimensional pore-scale numerical model was used to evaluate the impact of biofilm-induced flow heterogeneities on conservative transport. Our study integrates experimental biofilm images of Paenibacillus 300A strain in a microfluidic device packed with cylindrical grains in a hexagonal distribution, with mathematical modeling. Biofilm is represented as a synthetic porous structure with locally varying physical properties that honors the impact of biofilm on the porous medium. We find that biofilm plays a major role in shaping the observed conservative transport dynamics by enhancing anomalous characteristics. More specifically, when biofilm is present, the pore structure in our geometry becomes more spatially correlated. We observe intermittent behavior in the Lagrangian velocities that switches between fast transport periods and long trapping events. Our results suggest that intermittency enhances solute spreading in breakthrough curves which exhibit extreme anomalous slope at intermediate times and very marked late solute arrival due to solute retention. The efficiency of solute retention by the biofilm is controlled by a transport regime which can extend the tailing in the concentration breakthrough curves. These results indicate that solute retention by the biofilm exerts a strong control on conservative solute transport at pore-scale, a role that to date has not received enough attention.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

PICSAR-QED: a Monte Carlo module to simulate strong-field quantum electrodynamics in particle-in-cell codes for exascale architectures

Abstract Physical scenarios where the electromagnetic fields are so strong that quantum electrodynamics (QED) plays a substantial role are one of the frontiers of contemporary plasma physics research. Investigating those scenarios requires state-of-the-art particle-in-cell (PIC) codes able to run on top high-performance computing (HPC) machines and, at the same time, able to simulate strong-field QED processes. This work presents the PICSAR-QED library, an open-source, portable implementation of a Monte Carlo module designed to provide modern PIC codes with the capability to simulate such processes, and optimized for HPC. Detailed tests and benchmarks are carried out to validate the physical models in PICSAR-QED, to study how numerical parameters affect such models, and to demonstrate its capability to run on different architectures (CPUs and GPUs). Its integration with WarpX, a state-of-the-art PIC code designed to deliver scalable performance on upcoming exascale supercomputers, is also discussed and validated against results from the existing literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Learning PDFs through interpretable latent representations in Mellin space

Representing the parton distribution functions (PDFs) of the proton and other hadrons through flexible, high-fidelity parametrizations has been a long-standing goal of particle physics phenomenology. This is particularly true since the chosen parametrization methodology can play an influential role in the ultimate PDF uncertainties as extracted in QCD global analyses; these, in turn, are often determinative of the reach of experiments at the LHC and other facilities to nonstandard physics, including at large 𝑥, where parametrization effects can be significant. In this study, we explore a series of encoder-decoder machine-learning (ML) models with various neural-network topologies as efficient means of reconstructing PDFs from meaningful information stored in an interpretable latent space. Given recent effort to pioneer synergies between QCD analyses and lattice-gauge calculations, we formulate a latent representation based on the behavior of PDFs in Mellin space, i.e., their integrated moments, and test the ability of various models to decode PDFs from this information faithfully. We introduce a numerical package, PDFdecoder, which implements several encoder-decoder models to reconstruct PDFs with high fidelity and use this end-to-end tool to explore how such neural-network-based models might connect PDF parametrizations to underlying properties like their Mellin moments. We additionally dissect patterns of learned correlations between encoded Mellin moments and reconstructed PDFs that suggest opportunities for further improvements to ML-based approaches to PDF parametrizations and uncertainty quantification.

Machine learning↗
DOE OSTI2023

Continuum Damage Mechanics Modeling of High-Temperature Flaw Propagation: Application to Creep Crack Growth In 316H Standardized Specimens and Nuclear Reactor Components

Predicting creep crack growth (CCG) of flaws found during operation in high-temperature alloy components is essential for assessing the remaining lifetime of those components. While defect assessment procedures are available for this purpose in design codes, these are limited in their range of applicability. This study assesses the application of a local damage-based finite-element methodology as a more general technique for the prediction of CCG at high temperatures on a variety of structural configurations. Numerical results for stainless steel 316H, which are validated against experimental data, show the promise of this approach. This integration of continuum damage mechanics (CDM) based methodologies, together with adequate inelastic models, into assessment procedures can therefore inform the characterization of CCG under complex operating conditions, while avoiding excessive conservatism. This article shows that such modeling frameworks can be calibrated to experimental data and used to demonstrate that the degree of tri-axiality ahead of a growing creep crack affects its rate of growth. The framework is also successfully employed in characterizing CCG in a realistic reactor pressure vessel geometry under an arbitrary loading condition. These results are particularly relevant to the nuclear power industry for defect assessment and inspections as part of codified practices of structural components with flaws in high-temperature reactors.

42 ENGINEERING↗
DOE OSTI2021

Learning Forecasts of Rare Stratospheric Transitions from Short Simulations

Abstract Rare events arising in nonlinear atmospheric dynamics remain hard to predict and attribute. We address the problem of forecasting rare events in a prototypical example, sudden stratospheric warmings (SSWs). Approximately once every other winter, the boreal stratospheric polar vortex rapidly breaks down, shifting midlatitude surface weather patterns for months. We focus on two key quantities of interest: the probability of an SSW occurring, and the expected lead time if it does occur, as functions of initial condition. These optimal forecasts concretely measure the event’s progress. Direct numerical simulation can estimate them in principle but is prohibitively expensive in practice: each rare event requires a long integration to observe, and the cost of each integration grows with model complexity. We describe an alternative approach using integrations that are short compared to the time scale of the warming event. We compute the probability and lead time efficiently by solving equations involving the transition operator, which encodes all information about the dynamics. We relate these optimal forecasts to a small number of interpretable physical variables, suggesting optimal measurements for forecasting. We illustrate the methodology on a prototype SSW model developed by Holton and Mass and modified by stochastic forcing. While highly idealized, this model captures the essential nonlinear dynamics of SSWs and exhibits the key forecasting challenge: the dramatic separation in time scales between a single event and the return time between successive events. Our methodology is designed to fully exploit high-dimensional data from models and observations, and has the potential to identify detailed predictors of many complex rare events in meteorology.

Finkel, Justin↗
DOE OSTI2020

Diversified Therapeutic Phage Cocktails from Close Relatives of the Target Bacterium

This project tackles the antibiotic resistance crisis, developing a new method for discovering numerous efficacious bacteriophages for therapeutic cocktails against bacterial pathogens. The phage therapy approach to infectious disease, recently rekindled in U.S. medicine, requires numerous phages for each bacterial pathogen. Our approach 1) uses Sandia-unique software to identify dormant phages (prophages) integrated into bacterial chromosomes, 2) identifies prophage-laden bacteria that are close relatives of the target pathogenic strain to be killed, and 3) engineers away properties of these phages that are undesirable for therapy. We have perfected our phage-finding software, implemented our phage therapy strategy by targeting the pathogen Pseudomonas aeruginosa, and prepared new software to assist the phage engineering. We then turned toward Burkholderia pathogens, aiming to overcome the difficulty to transform these bacteria with a novel phage conjugation approach. Our work demonstrates the validity of a new approach to phage therapy for killing antibiotic resistant pathogens.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Computing the Instantaneous Collision Probability between Satellites using Characteristic Function Inversion

The probability that two satellites overlap in space at a specified instant of time is called their instantaneous collision probability. Assuming Gaussian uncertainties and spherical satellites, this probability is the integral of a Gaussian distribution over a sphere. This paper shows how to compute the probability using an established numerical procedure called characteristic function inversion. The collision probability in the short-term encounter scenario is also evaluated with this approach, where the instant at which the probability is computed is the time of closest approach between the objects. Python and R code is provided to evaluate the probability in practice. Overall, the approach has been established for over fifty years, is implemented in existing software, does not rely on analytical approximations, and can be used to evaluate two and three dimensional collision probabilities.

79 ASTRONOMY AND ASTROPHYSICS↗