Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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

Multi-objective optimization of peel and shear strengths in ultrasonic metal welding using machine learning-based response surface methodology

Ultrasonic metal welding (UMW) is a solid-state joining technique with varied industrial applications. Despite of its numerous advantages, UMW has a relative narrow operating window and is sensitive to variations in process conditions. As such, it is imperative to quantitatively characterize the influence of welding parameters on the resulting joint quality. The quantification model can be subsequently used to optimize the parameters. Conventional response surface methodology (RSM) usually employs linear or polynomial models, which may not be able to capture the intricate, nonlinear input-output relationships in UMW. Furthermore, some UMW applications call for simultaneous optimization of multiple quality indices such as peel strength, shear strength, electrical conductivity, and thermal conductivity. To address these challenges, this paper develops a machine learning (ML)-based RSM to model the input-output relationships in UMW and jointly optimize two quality indices, namely, peel and shear strengths. The performance of various ML methods including spline regression, Gaussian process regression (GPR), support vector regression (SVR), and conventional polynomial regression models with different orders is compared. A case study using experimental data shows that GPR with radial basis function (RBF) kernel and SVR with RBF kernel achieve the best prediction accuracy. The obtained response surface models are then used to optimize a compound joint strength indicator that is defined as the average of normalized shear and peel strengths. In addition, the case study reveals different patterns in the response surfaces of shear and peel strengths, which has not been systematically studied in the literature. While developed for the UMW application, the method can be extended to other manufacturing processes.

42 ENGINEERING↗
DOE OSTI2023

Parameterization and Explicit Modeling of Cloud Microphysics: Approaches, Challenges, and Future Directions

Cloud microphysical processes occur at the smallest end of scales among cloud-related processes and thus must be parameterized not only in large-scale global circulation models (GCMs) but also in various higher-resolution limited-area models such as cloud-resolving models (CRMs) and large-eddy simulation (LES) models. Instead of giving a comprehensive review of existing microphysical parameterizations that have been developed over the years, this study concentrates purposely on several topics that we believe are understudied but hold great potential for further advancing bulk microphysics parameterizations: multi-moment bulk microphysics parameterizations and the role of the spectral shape of hydrometeor size distributions; discrete vs “continuous” representation of hydrometeor types; turbulence-microphysics interactions including turbulent entrainment-mixing processes and stochastic condensation; theoretical foundations for the mathematical expressions used to describe hydrometeor size distributions and hydrometeor morphology; and approaches for developing bulk microphysics parameterizations. Also presented are the spectral bin scheme and particle-based scheme (especially, super-droplet method) for representing explicit microphysics. Their advantages and disadvantages are elucidated for constructing cloud models with detailed microphysics that are essential to developing processes understanding and bulk microphysics parameterizations. Particle-resolved direct numerical simulation (DNS) models are described as an emerging technique to investigate turbulence-microphysics interactions at the most fundamental level by tracking individual particles and resolving the smallest turbulent eddies in turbulent clouds. Outstanding challenges and future research directions are explored as well.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The EGS Collab Project: Learnings from Experiment 1

The primary objective of the EGS Collab Project sponsored by DOE is to increase the understanding needed to efficiently implement enhanced geothermal systems (EGS). One goal of the EGS Collab project is to create a collaborative research environment in which to study stimulation of crystalline rock at the 10 meter scale. Key to this effort is the collection of high quality data to allow comparison to numerical coupled process models in an effort to build confidence in the codes and modeling techniques used. In response to this, the EGS Collab team has created an underground test bed at the Sanford Underground Research Facility (SURF) in Lead SD at a depth of approximately 1.5 km to examine hydraulic fracturing (Experiment 1). We are currently designing a second test bed aimed at investigating shear stimulation (Experiment 2). At the Experiment 1 location, we have characterized our host rock using laboratory testing and numerous field-based geophysical and geological techniques, and created a well-instrumented test bed to allow us to carefully monitor stimulation events and flow tests. In addition to the installed geophysical sensors, we have used tracer tests, differences in the ambient microbial communities at flow collection locations, and cold water injection to inform us about dynamic flow pathways. In Experiment 1, we have hydraulically stimulated the host rock in a number times at several locations in one well, creating new fractures that connect to existing fractures between the injection and production boreholes. We have performed long-term ambient and chilled water injection tests as an analog to EGS, and have monitored system changes resulting from these water injections through geophysical monitoring, flow and pressure measurements, tracer tests, and microbiology. Here, we summarize the tests performed, issues identified including poroelastic and thermoelastic effects, Joule-Thomson effects, restarting effects, indications of flow channeling, and the primary learnings from Experiment 1.

Enhanced Geothermal Systems, EGS Collab, stimulati↗
DOE OSTI2021

Fast multiscale contrast independent preconditioners for linear elastic topology optimization problems

The goal of this work is to present a fast and viable approach for the numerical solution of the high-contrast state problems arising in topology optimization. The optimization process is iterative, and the gradients are obtained by an adjoint analysis, which requires the numerical solution of large high-contrast linear elastic problems with features spanning several length scales. The size of the discretized problems forces the utilization of iterative linear solvers with solution time dependent on the quality of the preconditioner. The lack of clear separation between the scales, as well as the high-contrast, imposes severe challenges on the standard preconditioning techniques. Thus, here we propose new methods for the high-contrast elasticity equation with performance independent of the high-contrast and the multi-scale structure of the elasticity problem. The solvers are based on two-levels domain decomposition techniques with a carefully constructed coarse level to deal with the high-contrast and multi-scale nature of the problem. The construction utilizes spectral equivalence between scalar diffusion and each displacement block of the elasticity problems and, in contrast to previous solutions proposed in the literature, is able to select the appropriate dimension of the coarse space automatically. The new methods inherit the advantages of domain decomposition techniques, such as easy parallelization and scalability. Finally, the presented numerical experiments demonstrate the excellent performance of the proposed methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗
DOE OSTI2022

Photoinduced evolution of lattice orthorhombicity and conceivably enhanced ferromagnetism in LaMnO 3 membranes

Ultrashort laser pulses have been utilized to dynamically drive phase transitions in correlated quantum materials. Of particular interest is whether phases not achievable in thermal equilibrium can be induced in complex oxides with intricately coupled lattice, electron and spin degrees of freedom. Here, we tracked atomic motions in LaMnO 3 following photoexcitation with MeV ultrafast electron diffraction (MeV-UED) technique. We found that the light excited state exhibits numerous signatures different from thermal equilibrium ones, including nearly conserved Bragg intensities, strongly suppressed La cation and oxygen anion displacements, and the long-range lattice orthorhombicity evolution. Furthermore, using first-principles calculations, we predict that the ferromagnetic ordering and conductivity are both enhanced upon laser excitation due to the reduction of the lattice orthorhombicity. This work benefits from recent advance in fabrication of membrane films with high epitaxial quality and in MeV-UED with large momentum space access and high temporal resolution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Characterizing the Sample Selection for Supernova Cosmology

Type Ia supernovae (SNe Ia) are used as distance indicators to infer the cosmological parameters that specify the expansion history of the universe. Parameter inference depends on the criteria by which the analysis SN sample is selected. Only for the simplest selection criteria and population models can the likelihood be calculated analytically, otherwise it needs to be determined numerically, a process that inherently has error. Numerical errors in the likelihood lead to errors in parameter inference. This article presents toy examples where the distance modulus is inferred given a set of SNe at a single redshift. Parameter estimators and their uncertainties are calculated using Monte Carlo techniques. The relationship between the number of Monte Carlo realizations and numerical errors is presented. The procedure can be applied to more realistic models and used to determine the computational and data management requirements of the transient analysis pipeline.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Surface Heater Fabrication Using Micro-Lithography for Transpiration Cooling Heat Transfer Coefficient Measurements

Abstract Transpiration cooling is able to provide more uniform coolant coverage than film cooling to effectively protect the component surface from contacting the hot gas. Due to numerous coolant ejection outlets within a small area at the target surface, the experimental thermo-fluid investigation on transpiration cooling becomes a significant challenge. Two classic methods to investigate film cooling, the steady-state foil heater method and the transient thermography technique, fail for transpiration cooling because the foil heater would block numerous coolant outlets, and the semi-infinite solid conduction model no longer holds for porous plates. In this study, a micro-lithography method to fabricate a silver coil pattern on top of the additively manufactured polymer porous media as the surface heater was proposed. The circuit was deliberately designed to cover the solid surface in a combination of series connection, and parallel connection to ensure the power in each unit cell area at the target surface was identical. With uniform heat flux generation, the steady-state tests were conducted to obtain distributions of a pair of parameters, adiabatic cooling effectiveness, and heat transfer coefficient (HTC). The results showed that the adiabatic cooling effectiveness could reach 0.65 with a blowing ratio lower than 0.5. Meanwhile, the heat transfer coefficient ratio (hf/h0) of transpiration cooling was close to 1 with a small blowing ratio at 0.125. A higher HTC ratio was observed for smaller pitch-to-diameter cases due to more turbulence intensity generated at the target surface.

Engineering↗
DOE OSTI2021

Injection data analysis using material balance time for CO 2 storage capacity estimation in deep closed saline aquifers

Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.

02 PETROLEUM↗
DOE OSTI2026

Identifying stochastic dynamics via finite expression methods

Modeling stochastic differential equations (SDEs) is crucial for understanding complex dynamical systems in various scientific fields. Recent methods often employ neural network-based models, which typically represent SDEs through a combination of deterministic and stochastic terms. However, these models usually lack interpretability and have difficulty in generalizing beyond their training domain. Here, this paper introduces the Finite Expression Method (FEX), a symbolic learning approach designed to derive interpretable mathematical representations of the deterministic component of SDEs. For the stochastic component, we integrate FEX with advanced generative modeling techniques to provide a comprehensive representation of SDEs. The numerical experiments on linear, nonlinear, and multidimensional SDEs demonstrate that FEX generalizes well beyond the training domain and delivers more accurate long-term predictions compared to neural network-based methods. The symbolic expressions identified by FEX not only improve prediction accuracy but also offer valuable scientific insights into the underlying dynamics of the systems.

Complex dynamical systems↗
DOE OSTI2021

The electronic structure and deexcitation pathways of an isolated metalloporphyrin ion resolved by metal L-edge spectroscopy

The local electronic structure of the metal-active site and the deexcitation pathways of metalloporphyrins are crucial for numerous applications but difficult to access by commonly employed techniques. Here, we applied near-edge X-ray absorption mass spectrometry and quantum-mechanical restricted active space calculations to investigate the electronic structure of the metal-active site of the isolated cobalt(III) protoporphyrin IX cation (CoPPIX+) and its deexcitation pathways upon resonant absorption at the cobalt L-edge. The experiments were carried out in the gas phase, thus allowing for control over the chemical state and molecular environment of the metalloporphyrin. The obtained mass spectra reveal that resonant excitations of CoPPIX+ at the cobalt L3-edge lead predominantly to the formation of the intact radical dication and doubly charged fragments through losses of charged and neutral side chains from the macrocycle. The comparison between experiment and theory shows that CoPPIX + is in a 3 A 2g triplet ground state and that competing excitations to metal-centred non-bonding and antibonding σ* molecular orbitals lead to distinct deexcitation pathways.

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

The Continuous Electron Beam Accelerator Facility at 12 GeV

This review paper describes the energy-upgraded Continuous Electron Beam Accelerator Facility (CEBAF) accelerator. This superconducting linac has achieved 12 GeV beam energy by adding 11 new high-performance cryomodules containing 88 superconducting cavities that have operated cw at an average accelerating gradient of 20 MV / m . After reviewing the attributes and performance of the previous 6 GeV CEBAF accelerator, we discuss the upgraded CEBAF accelerator system in detail with particular attention paid to the new beam acceleration systems. In addition to doubling the acceleration in each linac, the upgrade included improving the beam recirculation magnets, adding more helium cooling capacity to allow the newly installed modules to run cold, adding a new experimental hall, and improving numerous other accelerator components. We review several of the techniques deployed to operate and analyze the accelerator performance and document system operating experience and performance. In the final portion of the document, we present much of the current planning regarding projects to improve accelerator performance and enhance operating margins, and our plans for ensuring CEBAF operates reliably into the future. For the benefit of potential users of CEBAF, the performance and quality measures for the beam delivered to each of the experimental halls are summarized in the Appendix. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Deep reinforcement learning for complex evaluation of one-loop diagrams in quantum field theory

In this paper we present a technique based on deep reinforcement learning that allows for numerical analytic continuation of integrals that are often encountered in one-loop diagrams in quantum field theory. To extract certain quantities of two-point functions, such as spectral densities, mass poles or multiparticle thresholds, it is necessary to perform an analytic continuation of the correlator in question. At one-loop level in Euclidean space, this results in the necessity to deform the integration contour of the loop integral in the complex plane of the square of the loop momentum, to avoid nonanalyticities in the integration plane. Using a toy model for which an exact solution is known, we train a reinforcement learning agent to perform the required contour deformations. We report our study shows great promise for an agent to be deployed in iterative numerical approaches used to compute nonperturbative two-point functions, such as the quark propagator Dyson-Schwinger equation, or more generally, Fredholm equations of the second kind, in the complex domain.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Stochastic minibatch approach to the ptychographic iterative engine

The ptychographic iterative engine (PIE) is a widely used algorithm that enables phase retrieval at nanometer-scale resolution over a wide range of imaging experiment configurations. By analyzing diffraction intensities from multiple scanning locations where a probing wavefield interacts with a sample, the algorithm solves a difficult optimization problem with constraints derived from the experimental geometry as well as sample properties. The effectiveness at which this optimization problem is solved is highly dependent on the ordering in which we use the measured diffraction intensities in the algorithm, and random ordering is widely used due to the limited ability to escape from stagnation in poor-quality local solutions. In this study, we introduce an extension to the PIE algorithm that uses ideas popularized in recent machine learning training methods, in this case minibatch stochastic gradient descent. Our results demonstrate that these new techniques significantly improve the convergence properties of the PIE numerical optimization problem.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

High-Temperature and High-Pressure In Situ Magic Angle Spinning Nuclear Magnetic Resonance Spectroscopy

Nuclear magnetic resonance (NMR) spectroscopy represents an important technique to understand the structure and bonding environment of molecules. There exists a drive to conduct characterization experiment under conditions relevant to the chemical process of interest. To address this, in situ high-temperature, high-pressure MAS NMR methods have been developed to enable the observation of chemical interactions over a range of pressures (vacuum to several hundred bar) and temperatures (well below 0°C to 250°C). Further, the chemical identity of the samples can be comprised of solids, liquids, and gases or mixtures of the three. The method incorporates all-zirconia NMR rotors (sample holder for MAS NMR) which can be sealed using a threaded cap to compress an O-ring. This rotor exhibits great chemical resistance, temperature compatibility, low NMR background, and can withstand high pressures. These factors combined enable its use in a wide range of system combinations which permit its use in diverse fields as carbon sequestration, catalysis, material science, geochemistry, and biology. The flexibility of this technique makes it an attractive option for scientists from numerous disciplines.

Nuclear Magnetic Resonance, NMR, in situ, spectros↗
DOE OSTI2025

Strategies for simulating the time evolution of Hamiltonian lattice field theories

Simulating the time evolution of quantum field theories given some Hamiltonian H requires developing algorithms for implementing the unitary operator e -iHt . A variety of techniques exist that accomplish this task, with the most common technique used so far being Trotterization, which is a special case of the application of a product formula. However, other techniques exist that promise better asymptotic scaling in certain parameters of the theory being simulated, the most efficient of which are based on the concept of block encoding. In this work we study the performance of such algorithms in simulating lattice field theories. We derive and compare the asymptotic gate complexities of several commonly used simulation techniques in application to Hamiltonian lattice field theories. Using the scalar $\hat{φ}$ 4 theory as a test, we also perform numerical studies and compare the gate costs required by product formulas and signal-processing-based techniques to simulate time evolution. For the latter, we use the linear combination of unitaries (LCU) construction augmented with the quantum Fourier transform circuit to switch between the field and momentum eigenbases, which leads to immediate order-of-magnitude improvement in the cost of preparing the block encoding. Further, this paper also includes a pedagogical review of the techniques used, in particular product formulas, LCU, qubitization, quantum signal processing, as well as the technique for simulating geometrically-local Hamiltonians developed by Haah, Hastings, Kothari, and Low.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗
DOE OSTI2024

Artificial Intelligence/Machine Learning Technology in Power System Applications

The primary purpose of this report is to provide an overview of the advancement in artificial intelligence and machine learning (AI/ML) technologies and their applications in power systems. It offers a foundation for understanding the transformative role of AI/ML in power systems and aims to stimulate further research and development in this area. This report begins with a historical perspective of AI/ML technologies, then explores their advancement to today’s prominence. The document highlights key contributors to the success of AI/ML technologies, including increased computational power, greater data availability, innovative algorithms, and advanced tools. It further introduces various AI/ML techniques, including supervised, unsupervised and reinforcement learning, graph neural networks, and generative AI. It also emphasizes the critical importance of ensuring the safety, security, and trustworthiness of these AI/ML techniques within this sector. The report reviews the recent representative advancements in various power system applications enhanced by AI/ML techniques, underscoring key developments and their transformative impact as evidenced by numerous studies. It also explores both the opportunities and challenges associated with the application of AI/ML technologies to improve power system applications. While the report extensively covers AI/ML applications in power systems, focusing primarily on the technical and operational aspects, it may not thoroughly explore the sociopolitical, economic, and broader regulatory implications of AI/ML integration in power systems. AI/ML techniques hold significant potential for enhancing power system applications; however, they are not omnipotent. It is crucial to acknowledge their limitations and understand that they may not be able to address all challenges in the power system domain. Various factors must be considered that influence the implementation, adoption, and effectiveness of AI/ML solutions, including but not limited to safety, security, transparency, and trustworthiness. Additionally, the incorporation of advanced human–machine interfaces is essential, as it enables humans to validate the effectiveness of AI/ML solutions while remaining actively engaged, fostering trust in AI/ML deployment. Finally, the report summarizes AI/ML research activities supported by the Department of Energy (DOE) Office of Electricity (OE) through the Advanced Grid Modeling (AGM) program. The work aligns with the interests and mission of DOE-OE AGM, with the report serving as a resource for identifying existing progress and for pinpointing future applications within AI/ML that need further exploration and support.

24 POWER TRANSMISSION AND DISTRIBUTION↗