Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Pore architecture controls on mineral reactivity

Mineral dissolution rates measured in natural environments are much slower than those measured in laboratory settings. This project tested the hypothesis that the way fluid flows through rocks in natural systems creates areas where mineral dissolution is fast and areas where mineral dissolution is slow. This hypothesis was tested with a combination of laboratory experiments and numerical simulation. We demonstrated a separation of fluid flow pathways and rates of mineral dissolution in laboratory experiments for the first time using an experimental approach where we created synthetic rocks that have different ratios of connected and dead-end pathways for fluid flow. In laboratory experiments with higher proportions of dead-end pathways, the mineral dissolution rates were slower. We found that where fluid flows through connected pathways the continuous refreshing of fluid at the mineral surface creates conditions where dissolution is fast. In contrast, where fluid either flows slowly through poorly connected pathways or is stagnant in dead-end pathways, mineral dissolution is slow. The results from this project suggest that the overall slowing of rates of mineral dissolution is important when the proportion of dead-end pathways is greater than ~40%. This project informs our understanding of the way that fluids react with rocks in carbon dioxide sequestration and enhanced geothermal projects where fluids are purposefully injected into rocks for energy applications.

58 GEOSCIENCES↗
DOE OSTI2025

Relativistic gas accretion onto supermassive black hole binaries from inspiral through merger

Accreting supermassive black hole binaries are powerful multimessenger sources emitting both gravitational and electromagnetic (EM) radiation. Understanding the accretion dynamics of these systems and predicting their distinctive EM signals is crucial to informing and guiding upcoming efforts aimed at detecting gravitational waves produced by these binaries. To this end, accurate numerical modeling is required to describe both the spacetime and the magnetized gas around the black holes. In this paper, we present two key advances in this field of research. First, we have developed a novel 3D general relativistic magnetohydrodynamics (GRMHD) framework that combines multiple numerical codes to simulate the inspiral and merger of supermassive black hole binaries starting from realistic initial data and running all the way through merger. Throughout the evolution, we adopt a simple but functional prescription to account for gas cooling through photon emission. Next, we have applied our new computational method to follow the time evolution of a circular, equal-mass, nonspinning black hole binary for ∼200 orbits, starting from a separation of 20⁢𝑟 𝑔 and reaching the postmerger evolutionary stage of the system. We have shown how mass continues to flow toward the binary even after the binary “decouples” from its surrounding disk, but the accretion rate onto the black holes diminishes. We have identified how the minidisks orbiting each black hole are slowly drained and eventually dissolve as the binary compresses. We confirm previous findings that the system’s luminosity decreases by a factor of a few during inspiral; however, we observe an abrupt increase by ∼50% in this quantity at the time of merger, likely accompanied by an equally abrupt change in spectrum. Lastly, we have demonstrated that during the inspiral, fluid ram pressure regulates the fraction of the magnetic flux transported to the binary that attaches to the black holes’ horizons.

Accretion disk & black-hole plasma↗
DOE OSTI2020

Optimal economic operation of liquid petroleum products pipeline systems

The majority of overland transport needs for crude petroleum and refined petroleum products are met using pipelines. Numerous studies have developed optimization methods for design of these systems in order to minimize construction costs while meeting capacity requirements. Here, we formulate problems to optimize the operations of existing single liquid commodity pipeline systems subject to physical flow and pump engineering constraints. The objectives are to maximize the economic value created for users of the system and to minimize operating costs. We present a general computational method for this class of continuous, non-convex nonlinear programs, and examine the use of pump operating settings and flow allocations as decision variables. Furthermore, the approach is applied to compute optimal operating regimes and perform engineering economic sensitivity analyses for a case study of a crude oil pipeline developed using publicly available data.

02 PETROLEUM↗
DOE OSTI2020

Scattering using real-time path integrals

Background: Path integrals are a powerful tool for solving problems in quantum theory that are not amenable to a treatment by perturbation theory. Most path integral computations require an analytic continuation to imaginary time. While imaginary time treatments of scattering are possible, imaginary time is not a natural framework for treating scattering problems. More importantly, quantum algorithms for calculating path integrals require real-time evolution. Purpose: Here, we test a recently introduced method for performing direct calculations of scattering observables using real-time path integrals in order to understand the challenges facing real-time path integral calculations of scattering observables. Method: The computations are based on a new interpretation of the path integral as the expectation value of a potential functional on cylinder sets of continuous paths with respect to a complex probability distribution. The method can in principle be applied to arbitrary short-range potentials. Results: The method is applied to compute matrix elements of Møller wave operators applied to narrow wave packets. These are used to calculate half-shell sharp-momentum transition matrix elements for one-dimensional potential scattering. The calculations for half-shell transition operator matrix elements converge to the numerical solution of the Lippmann-Schwinger equation. Conclusions: This work presents a proof in principle that scattering observables can be computed using real-time Feynman path integrals. While the computational method is not efficient, it can be improved. It provides a laboratory for studying quantum computational algorithms that are applicable to scattering problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Initial OpenStudio Profiling Results

OpenStudio’s performance has not historically been an area of much work, but as it has successfully replaced ad hoc model generation solutions, the performance of the software is more and more central to continuing success. This report describes an initial effort to profile OpenStudio, describes the problems encountered, the solutions to those problems, and some early recommendations for further work should funding become available. The approach taken here is to use special software, referred to as profilers, to assess the code and how it executes. This approach is more appropriate for this kind of software than the checkpoint-style timing that is often done with numerical codes. Profiling was most successful on the MacOS platform, where Apple’s Instruments software was able to decipher the complexities of OpenStudio’s command line execution of a workflow. Even with the limited exploration of performance done here, the team quickly ran into limitations imposed on the code by the stateless architecture, and the team recommends an evaluation of this architecture as a good next step to improve performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

PHEV Distribution Grid Integration and Smart Systems Testing: Cooperative Research and Development Final Report

To continue NREL's mission, which includes preparing technologies and markets with "speed and scale," NREL researchers will collaborate with Toyota on plug-in electric vehicle technology. Toyota will provide a collection of plug-in hybrid electric vehicles ("PHEVs" [and ultimately including PEVs], collectively, the "Vehicles") to conduct infrastructure interaction testing at NREL facilities. NREL will use available resources in the Vehicle Testing and Integration Facility (VTIF), its staff parking garage with numerous commercial Electric Vehicle Supply Equipment ("EVSEs"), and its Energy Systems Integration Facility (ESIF) to safely execute research experiments that highlight potential distribution system power quality challenges related to plug-in electric vehicles. The expected outcome of this project is a better understanding of the correlation of individual vehicle power quality attributes to a system of vehicles on a distribution network. The data collected and the simulations conducted will be used to guide future experiments and project development that aid in the deployment of plug-in electric vehicles ("PEV").

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Acceptance criteria for in situ surveillance of MSR materials based on thermally-loaded mechanical test articles

This report describes practices and acceptance test procedures for designing, running, and maintaining a material surveillance program in a future operating molten salt reactor. The programs described here rely on test data from passively actuated mechanical test articles inserted into critical regions of the reactor and periodically removed for out-of-reactor testing. The report defines definite acceptance procedures, based on the results of these tests, to determine whether a component can continue to operate accounting for the accumulation of environmentally-assisted mechanical damage in the component materials to date, and extrapolated out through the next inspection period. Additionally, the report describes work on a software tool implementing many of the surveillance methods and procedures described here and progress on simplified methods for inferring damage accumulation in the test articles, based on out-of-reactor thermal cycling, that do not rely on sophisticated numerical analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Scientific challenges to characterizing the wind resource in the marine atmospheric boundary layer

Abstract. With the increasing level of offshore wind energy investment, it is correspondingly important to be able to accurately characterize the wind resource in terms of energy potential as well as operating conditions affecting wind plant performance, maintenance, and lifespan. Accurate resource assessment at a particular site supports investment decisions. Following construction, accurate wind forecasts are needed to support efficient power markets and integration of wind power with the electrical grid. To optimize the design of wind turbines, it is necessary to accurately describe the environmental characteristics, such as precipitation and waves, that erode turbine surfaces and generate structural loads as a complicated response to the combined impact of shear, atmospheric turbulence, and wave stresses. Despite recent considerable progress both in improvements to numerical weather prediction models and in coupling these models to turbulent flows within wind plants, major challenges remain, especially in the offshore environment. Accurately simulating the interactions among winds, waves, wakes, and their structural interactions with offshore wind turbines requires accounting for spatial (and associated temporal) scales from O(1 m) to O(100 km). Computing capabilities for the foreseeable future will not be able to resolve all of these scales simultaneously, necessitating continuing improvement in subgrid-scale parameterizations within highly nonlinear models. In addition, observations to constrain and validate these models, especially in the rotor-swept area of turbines over the ocean, remains largely absent. Thus, gaining sufficient understanding of the physics of atmospheric flow within and around wind plants remains one of the grand challenges of wind energy, particularly in the offshore environment. This paper provides a review of prominent scientific challenges to characterizing the offshore wind resource using as examples phenomena that occur in the rapidly developing wind energy areas off the United States. Such phenomena include horizontal temperature gradients that lead to strong vertical stratification; consequent features such as low-level jets and internal boundary layers; highly nonstationary conditions, which occur with both extratropical storms (e.g., nor'easters) and tropical storms; air–sea interaction, including deformation of conventional wind profiles by the wave boundary layer; and precipitation with its contributions to leading-edge erosion of wind turbine blades. The paper also describes the current state of modeling and observations in the marine atmospheric boundary layer and provides specific recommendations for filling key current knowledge gaps.

17 WIND ENERGY↗
DOE OSTI2026

Harnessing the Second-Order Metal−Insulator Transition for Neuromorphic Computing

Vanadium oxides are widely studied phase change materials for brain-inspired computing architectures. Systems like VO 2 and V 2 O 3 exhibit first-order metal−insulator transitions (MITs) with hysteresis and percolative switching, increasing stochasticity and device variability. Here, we focus on the less explored Magnéli phase V 4 O 7 , which undergoes a continuous, non-hysteretic, second-order MIT. This surprisingly enables highly reproducible volatile resistive switching in spiking-neuron-type devices. We synthesize V 4 O 7 films, characterize their structural and transport properties, and demonstrate voltage and current-driven threshold switching with electrothermal feedback. In a Pearson–Anson oscillator, V 4 O 7 devices produce stable, tunable spiking across 20–200 kHz, with consistent operation among multiple devices. We introduce a numerical analog leaky-integrate-and-fire (aLIF) model that captures waveform shapes and their dependence on load resistance, temperature, and voltage. Furthermore, these findings suggest that second-order MIT materials like V 4 O 7 are promising for deterministic, scalable spiking neuron arrays for neuromorphic computing.

V4O7↗
DOE OSTI2023

On fast simulation of dynamical system with neural vector enhanced numerical solver

The large-scale simulation of dynamical systems is critical in numerous scientific and engineering disciplines. However, traditional numerical solvers are limited by the choice of step sizes when estimating integration, resulting in a trade-off between accuracy and computational efficiency. To address this challenge, we introduce a deep learning-based corrector called Neural Vector (NeurVec), which can compensate for integration errors and enable larger time step sizes in simulations. Our extensive experiments on a variety of complex dynamical system benchmarks demonstrate that NeurVec exhibits remarkable generalization capability on a continuous phase space, even when trained using limited and discrete data. NeurVec significantly accelerates traditional solvers, achieving speeds tens to hundreds of times faster while maintaining high levels of accuracy and stability. Moreover, NeurVec’s simple-yet-effective design, combined with its ease of implementation, has the potential to establish a new paradigm for fast-solving differential equations based on deep learning.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

An experimental database of cell performance for vanadium redox flow battery

The continual growth in energy demand has resulted in the deployment of renewable energy generators to reduce the impact of fossil fuel dependence. However, these generators often suffer from intermittency and require energy storage when there is over-generation and the subsequent release of this stored energy at high demand. One promising energy storage technology which can provide a solution to improve energy management and grid stability, is the redox flow battery. Among the numerous flow battery systems, vanadium redox flow battery is the most iconic solution to large scale energy storage, giving a more efficient link between energy production, especially from renewables, and energy demand. The aim of the current database is to characterize the performance of the cell design and to provide training/validation data for physical model or data-driven model. The database includes hundreds of experimental cell performance data of vanadium redox flow battery with various current densities for multiple charge-discharge cycles. All the cell parameters, chemical parameters, material parameters, operation parameters and thermodynamic parameters of the cell system are listed. Coulomb, voltaic and energy efficiencies are also provided. The database will be helpful for researchers in the field of redox flow batteries.

Gao, Peiyuan↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2025

A new tubular ironless linear induction motor using non-interrupted winding method

The ironless linear induction motor device can be used for any application which requires a large motion span or rotation restrictedly by the ironless feature due to the strong magnetic field environment. One of the possible applications is ITER diagnostics shutters. The tokamak reactor systems have a severe radiation issue during the operation. The traditional Kapton-insulated wires have significantly lower radiation resistance than an expected large radiation dose at ITER Motion Stark Effect (MSE) shutter location. Furthermore, the traditional single-sided linear induction motor winds the wire to make the N-turn pitch coils independently and then connect them in series to make a one-phase winding. Therefore, it is a big challenge to make so many radiation-protection connectors in such a limited space as inside the vacuum vessel of fusion reactors. The approved mineral-insulated kinds of cables need a large bend radius for the end turns, so the traditional single-sided linear motor with two-layer three-phase windings need to be upgraded to the tubular linear motors. The new proposed tubular ironless linear induction motor can not only overcome the large bend radius issue with the mineral-insulated cables but also allow the motor coils wound continuously around the tubular-shaped stator, forward and backward N times to make a N-turn winding for each phase, in order to resolve the radiation protection issue with the winding connection. The 3D transient Maxwell electromagnetic model with a tubular three-phase linear induction motor has been analyzed. Furthermore, this paper will present the detailed design and winding method of ironless linear induction motor actuators, and results of the 3D transient numerical simulation of the actuator performance.

Ironless↗
DOE OSTI2023

Finite-difference time-domain methods

The finite-difference time-domain (FDTD) method is a widespread numerical tool for full-wave analysis of electromagnetic fields in complex media and for detailed geometries. Applications of the FDTD method cover a range of time and spatial scales, extending from subatomic to galactic lengths and from classical to quantum physics. Technology areas that benefit from the FDTD method include biomedicine — bioimaging, biophotonics, bioelectronics and biosensors; geophysics — remote sensing, communications, space weather hazards and geolocation; metamaterials — sub-wavelength focusing lenses, electromagnetic cloaks and continuously scanning leaky-wave antennas; optics — diffractive optical elements, photonic bandgap structures, photonic crystal waveguides and ring-resonator devices; plasmonics — plasmonic waveguides and antennas; and quantum applications — quantum devices and quantum radar. This Primer summarizes the main features of the FDTD method, along with key extensions that enable accurate solutions to be obtained for different research questions. Additionally, hardware considerations are discussed, plus examples of how to extract magnitude and phase data, Brillouin diagrams and scattering parameters from the output of an FDTD model. Furthermore, the Primer ends with a discussion of ongoing challenges and opportunities to further enhance the FDTD method for current and future applications.

42 ENGINEERING↗
DOE OSTI2024

Hot-Roll Fabrication of Anisotropic Nanograin Nd-Fe-B Magnet

Nd-Fe-B based magnets have the highest energy product among all permanent magnets, which is required for numerous clean energy technologies. For higher temperature applications (T > 150°C), additions of heavy rare earth elements (HREEs) such as Dy are required to maintain sufficient coercivity during operation. Additions of Dy are expensive. Thus, it is desirable to reduce the need for HREEs by reducing the grain size to the nanoscale, which increases the coercivity and decreases its temperature dependence. Here, we report a novel nanograin Nd-Fe-B magnet fabrication method that is continuous and inexpensive. The process uses mechanically milled Nd-Fe-B melt-spun flakes as feedstock powder that is packed into a metal vessel and then hot rolled to form a fully dense and highly textured strip magnet with tailored thicknesses, down to 800 µm. Finally, using this process, fully dense nanograin bulk magnets can be synthesized in minutes compared to the traditional multi-step processes that are typically low throughput.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Comparative analysis of atmospheric parameters from high-resolution spectroscopic sky surveys: APOGEE, GALAH, Gaia -ESO

SDSS-IV APOGEE-2, GALAH, and Gaia-ESO are high-resolution, ground-based, multi-object spectroscopic surveys providing fundamental stellar atmospheric parameters and multiple elemental abundance ratios for hundreds of thousands of stars of the Milky Way. Data from these and other surveys contribute to investigations of the history and evolution of the Galaxy. We undertake a comparison between the most recent data releases of these surveys to investigate the accuracy and precision of derived parameters by placing the abundances on an absolute scale. We also discuss the correlations in parameter and abundance differences as a function of main parameters. Uncovering the variants provides a basis to continue the efforts of future sky surveys. Quality samples from the APOGEE-GALAH (15 537 stars), APOGEE-GES (804 stars), and GALAH-GES (441 stars) overlapping catalogs were collected. We investigated the mean variants between the surveys, and linear trends were also investigated. We compared the slope of correlations and mean differences with the reported uncertainties. The average and scatter of v rad , T eff , log g, [M/H], and v micro , along with numerous species of elemental abundances in the combined catalogs, show that in general there is a good agreement between the surveys. We find large radial velocity scatters ranging from 1.3 km s -1 to 4.4 km s -1 when comparing the three surveys. We observe some weak trends (e.g., in ΔT eff vs. Δlog g for the APOGEE-GES stars) and a clear correlation in the v micro - Δv micro planes in the APOGEE-GALAH common sample. For [α/H], [Ti/H] (APOGEE-GALAH giants), and [Al/H] (APOGEE-GALAH dwarfs) potential strong correlations are discovered as a function of the differences in the main atmospheric parameters, and we also find weak trends for other elements. In general we find good agreement between the three surveys within their respective uncertainties. However, there are certain regimes in which strong variants exist, which we discuss. There are still offsets larger than 0.1 dex in the absolute abundance scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Stochastic Learning Approach for Binary Optimization: Application to Bayesian Optimal Design of Experiments

Here, we present a novel stochastic approach to binary optimization suited for optimal experimental design (OED) for Bayesian inverse problems governed by mathematical models such as partial differential equations. The OED utility function, namely, the regularized optimality criterion, is cast into a stochastic objective function in the form of an expectation over a multivariate Bernoulli distribution. The probabilistic objective is then solved by using a stochastic optimization routine to find an optimal observational policy. This formulation (a) is generally applicable to binary optimization problems with soft constraints and is ideal for OED and sensor placement problems; (b) does not require differentiability of the original objective function (e.g., a utility function in OED applications) with respect to the design variable, and thus it enables direct employment of sparsity-enforcing penalty functions such as $\ell_0$, without needing to utilize a continuation procedure or apply a rounding technique; (c) exhibits much lower computational cost than traditional gradient-based relaxation approaches; and (d) can be applied to both linear and nonlinear OED problems with proper choice of the utility function. The proposed approach is analyzed from an optimization perspective with detailed convergence analysis of the optimization approach and is also analyzed from a machine learning perspective with correspondence to policy gradient reinforcement learning. The approach is demonstrated numerically by using an idealized two-dimensional Bayesian linear inverse problem and validated by extensive numerical experiments carried out for sensor placement in a parameter identification setup.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Asynchronous distributed-memory task-parallel algorithm for compressible flows on unstructured 3D Eulerian grids

Here, we discuss the implementation of a finite element method, used to numerically solve the Euler equations of compressible flows, using an asynchronous runtime system (RTS). The algorithm is implemented for distributed-memory machines, using stationary unstructured 3D meshes, combining data-, and task-parallelism on top of the Charm++ RTS. Charm++’s execution model is asynchronous by default, allowing arbitrary overlap of computation and communication. Task-parallelism allows scheduling parts of an algorithm independently of, or dependent on, each other. Built-in automatic load balancing enables continuous redistribution of computational load by migration of work units based on real-time CPU load measurement. The RTS also features automatic checkpointing, fault tolerance, resilience against hardware failure, and supports power-, and energy-aware computation. We demonstrate scalability up to 25 x 10 9 cells at $\mathscr{O}$10 4 compute cores and the benefits of automatic load balancing for irregular workloads. The full source code with documentation is available at https://quinoacomputing.org.

42 ENGINEERING↗