Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical solution”

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 757 records · Page 42

Numerical investigation of the AP1000 response following loss-of-coolant accident using PCTRAN and CFD to support R&D of SMRs

Small Modular Reactors (SMRs) present a promising solution for the future of sustainable energy, offering advantages such as reduced waste generation, advanced passive safety features, and potential cost efficiencies. Ongoing research focuses on the design and development of SMRs, addressing challenges through numerical simulations and experimental test data. Given that many next-generation reactors are first-of-a-kind (FOK), current operational reactors serve as valuable benchmarks for understanding various thermal-hydraulic phenomena during postulated design basis accidents (DBAs), including loss-of-coolant accidents (LOCA), main steam line breaks (MSLB), and steam generator tube ruptures (SGTR). This study employs the AP1000 model Personal Computer Transient Analyzer (PCTRAN) to simulate reactor responses to small-break LOCA scenarios, specifically analyzing break sizes of 2 and 10 inches. Additionally, ANSYS FLUENT software is utilized to assess containment responses to large-break LOCAs, focusing on the quantification of decay heat removal via natural convection. Verification of the PCTRAN results is achieved using Westinghouse data, with CFD results corroborating the findings. The outcomes demonstrate strong agreement with Westinghouse data, confirming the accuracy of the simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Novel Co‐Catalytic Activities of Solid and Liquid Phase Catalysts in High‐Rate Li‐Air Batteries

Abstract Li‐air batteries are considered strong candidates for the next‐generation energy storage systems designed for electrical transportation. However, low cyclability and current rates are two major drawbacks that hinder them from further realization. These issues necessitate the discovery of novel materials to significantly enhance the redox process of discharge products. In this study, a novel catalytic system comprised of tin sulfide (SnS) nanoflakes as a solid catalyst and tin iodide (SnI 2 ) as a dual‐functional electrolyte additive is discovered. This system enables operating the battery at high current rates up to 10 000 mA g −1 (corresponding to 1 mA cm −2 ). The SnS catalyst shows outstanding catalytic activity for both oxygen reduction and evolution reactions compared to carbon, noble metals, and other transition metal dichalcogenides. It also exhibits good structural integrity at high rates. The computations indicate numerous possible oxygen reduction sites without oxygen dissociations on the SnS surface through solution mechanism that is likely responsible for the formation of Li 2 O 2 . The calculations also indicate that the role of the SnI 2 is not only reacting with the lithium anode to provide protection but reducing the charge potential by promoting catalytic decomposition of the Li 2 O 2 . This work provides new novel additives for designing high‐rate Li‐air batteries.

25 ENERGY STORAGE↗

Stochastic projective splitting

Here, we present a new, stochastic variant of the projective splitting (PS) family of algorithms for inclusion problems involving the sum of any finite number of maximal monotone operators. This new variant uses a stochastic oracle to evaluate one of the operators, which is assumed to be Lipschitz continuous, and (deterministic) resolvents to process the remaining operators. Our proposal is the first version of PS with such stochastic capabilities. We envision the primary application being machine learning (ML) problems, with the method’s stochastic features facilitating “mini-batch” sampling of datasets. Since it uses a monotone operator formulation, the method can handle not only Lipschitz-smooth loss minimization, but also min–max and noncooperative game formulations, with better convergence properties than the gradient descent-ascent methods commonly applied in such settings. The proposed method can handle any number of constraints and nonsmooth regularizers via projection and proximal operators. We prove almost-sure convergence of the iterates to a solution and a convergence rate result for the expected residual, and close with numerical experiments on a distributionally robust sparse logistic regression problem.

97 MATHEMATICS AND COMPUTING↗

Joint optimization of electric bus charging infrastructure, vehicle scheduling, and charging management

High upfront costs of vehicles and charging infrastructure as well as the lack of knowledge related to infrastructure planning and electric bus system operation are major obstacles to the implementation of battery electric buses (BEBs). To tackle the obstacles and promote BEB adoption, a comprehensive optimization framework was developed to address the combined charging infrastructure planning, vehicle scheduling, and charging management problem for BEB systems, with the goal to minimize the total cost of ownership. The problem was formulated as a mixed-integer non-linear problem. A genetic algorithm-based approach was then proposed to solve the problem. Last, three alternative scenarios based on a sub-transit network in Salt Lake City, Utah, were analyzed and compared with the optimal scenario results in the numerical experiments. Our comparison results demonstrate the effectiveness of the proposed model and solution algorithm in determining a cost-efficient planning strategy for BEB systems.

33 ADVANCED PROPULSION SYSTEMS↗

In Situ Velocity‐Strain Sensitivity Near the San Jacinto Fault Zone Analyzed Through Train Tremors

Abstract We utilize train tremors as P‐wave seismic sources to investigate velocity‐strain sensitivity near the San Jacinto Fault Zone. A dense nodal array deployed at the Piñon Flat Observatory is used to detect and identify repeating train energy emitted from a railway in the Coachella valley. We construct P‐wave correlation functions across the fault zone and estimate the spatially averaged dt / t versus strain sensitivity to be 6.25 × 10 4 . Through numerical simulations, we explore how the sensitivity decays exponentially with depth. The optimal solution reveals a subsurface sensitivity of 1.2 × 10 5 and a depth decay rate of 0.05 km −1 . This sensitivity aligns with previous findings but is toward the higher end, likely due to the fractured fault‐zone rocks. The depth decay rate, previously unreported, is notably smaller than assumed in empirical models. This raises the necessity of further investigations of this parameter, which is crucial to study stress and velocity variations at seismogenic depth.

Geology↗

Gluon helicity in the nucleon from lattice QCD and machine learning

We present the first lattice QCD determination of the light cone gluon helicity correlation parton distribution function (PDF) with numerical evidence toward disfavoring negative gluon polarization in the nucleon. We present a solution for eliminating an inevitable contamination term that dominates the Euclidean correlations and makes determining gluon helicity PDF unfeasible. The proposed synergy between lattice QCD and artificial intelligence offers a superior platform to alleviate the defining challenge of extracting quark and gluon PDFs from the lattice data that are available in a limited domain due to a finite range of accessible hadron momenta. We suggest a systematically improvable method to extract PDFs from the lattice data, independent of inadequate parametrizations. The result of the gluon helicity will improve our understanding of the role of spin in the strong interaction and the nucleon-spin structure.

Lattice QCD↗

Modifying the Asynchronous Jacobi Method for Data Corruption Resilience

Moving scientific computation from high-performance computing (HPC) and cloud computing (CC) environments to devices on the edge, i.e., physically near instruments of interest, has received tremendous interest in recent years. Such edge computing environments can operate on data in situ, offering enticing benefits over data aggregation to HPC and CC facilities that include avoiding costs of transmission, increased data privacy, and real-time data analysis. Because of the inherent unreliability of edge computing environments, new fault-tolerant approaches must be developed before the benefits of edge computing can be realized. Motivated by algorithm-based fault tolerance, a variant of the asynchronous Jacobi (ASJ) method is developed that achieves resilience to data corruption by rejecting solution approximations from neighbor devices according to a bound derived from convergence theory. Numerical results on a two-dimensional Poisson problem show that the new rejection criterion, along with a novel approximation to the shortest path length on which the criterion depends, restores convergence for the ASJ variant in the presence of certain types data corruption. Numerical results are obtained for when the singular values in the analytic bound are approximated. Additional linear systems are also explored, one with a more dense sparsity pattern and one that includes advection. All results indicate that successful resilience to data corruption depends on whether the bound tightens fast enough to reject corrupted data before the iteration evolution deviates significantly from that predicted by the convergence theory defining the bound. This observation generalizes to future work on algorithm-based fault tolerance for other asynchronous algorithms, including upcoming approaches that leverage Krylov subspaces.

97 MATHEMATICS AND COMPUTING↗

Numerical methods for fractional Fokker–Planck equation with multiplicative Marcus Lévy noises

The Fokker–Planck equation (FPE) is an important deterministic tool for investigating stochastic dynamical systems. In this paper, we consider the space-time fractional FPE driven by multiplicative Marcus Lévy noises. Efficient numerical schemes are presented to solve the equations. Stability and convergence of the methods are also discussed. We give some numerical experiments to validate our schemes, and examine the effects of parameters on solutions. Additionally, we analyze the maximal likely trajectories and the critical time for the change of the most probability location.

Mathematics↗

Bayesian inference of Stochastic reaction networks using Multifidelity Sequential Tempered Markov Chain Monte Carlo

Stochastic reaction network models are often used to explain and predict the dynamics of gene regulation in single cells.These models usually involve several parameters, such as the kinetic rates of chemical reactions, that are not directly measurable and must be inferred from experimental data. Bayesian inference provides a rigorous probabilistic frame-work for identifying these parameters by finding a posterior parameter distribution that captures their uncertainty.Traditional computational methods for solving inference problems such as Markov Chain Monte Carlo methods based on classical Metropolis-Hastings algorithm involve numerous serial evaluations of the likelihood function, which in turn requires expensive forward solutions of the chemical master equation (CME). We propose an alternate approach based on a multifidelity extension of the Sequential Tempered Markov Chain Monte Carlo (ST-MCMC) sampler. This algorithm is built upon Sequential Monte Carlo and solves the Bayesian inference problem by decomposing it into a sequence of efficiently solved subproblems that gradually increase both model fidelity and the influence of the observed data. We reformulate the finite state projection (FSP) algorithm, a well-known method for solving the CME, to produce a hierarchy of surrogate master equations to be used in this multifidelity scheme. To determine the appropriate fidelity, we introduce a novel information-theoretic criteria that seeks to extract the most information about the ultimate Bayesian posterior from each model in the hierarchy without inducing significant bias. This novel sampling scheme is tested with high performance computing resources using biologically relevant problems.

97 MATHEMATICS AND COMPUTING↗

Neuromorphic scaling advantages for energy-efficient random walk computations

Computing stands to be radically improved by neuromorphic computing (NMC) approaches inspired by the brain's incredible efficiency and capabilities. Most NMC research, which aims to replicate the brain's computational structure and architecture in man-made hardware, has focused on artificial intelligence; however, less explored is whether this brain-inspired hardware can provide value beyond cognitive tasks. We demonstrate that high-degree parallelism and configurability of spiking neuromorphic architectures makes them well-suited to implement random walks via discrete time Markov chains. Such random walks are useful in Monte Carlo methods, which represent a fundamental computational tool for solving a wide range of numerical computing tasks. Additionally, we show how the mathematical basis for a probabilistic solution involving a class of stochastic differential equations can leverage those simulations to provide solutions for a range of broadly applicable computational tasks. Despite being in an early development stage, we find that NMC platforms, at a sufficient scale, can drastically reduce the energy demands of high-performance computing platforms.

59 BASIC BIOLOGICAL SCIENCES↗

The Technical and Economic Potential of the H2@Scale Hydrogen Concept within the United States

The U.S. energy system is evolving as society and technologies change. Renewable electricity generation - especially from wind and solar - is growing rapidly, and alternative energy sources are being developed and implemented across the residential, commercial, transportation, and industrial sectors to take advantage of their cost, security, and health benefits. Systemic changes present numerous challenges to grid resiliency and energy affordability, creating a need for synergistic solutions that satisfy multiple applications while yielding system-wide cost and emissions benefits. One such solution is an integrated hydrogen energy system. This is the focus of H2@Scale - a U.S. Department of Energy (DOE) initiative led by the Office of Energy Efficiency and Renewable Energy’s Hydrogen and Fuel Technologies Office. H2@Scale brings together stakeholders to advance affordable hydrogen production, transport, storage, and utilization in multiple energy sectors. The H2@Scale concept involves hydrogen as an energy intermediate. Hydrogen can be produced from various conventional and renewable energy sources including as a responsive load on the electric grid. Hydrogen has many current applications and many more potential applications, such as energy for transportation - used directly in fuel cell electric vehicles (FCEVs), as a feedstock for synthetic fuels, and to upgrade oil and biomass - feedstock for industry (e.g., for ammonia production, metals refining, and other end uses), heat for industry and buildings, and electricity storage. Owing to its flexibility and fungibility, a hydrogen intermediate could link energy sources that have surplus availability to markets that require energy or chemical feedstocks, benefiting both. This document builds upon a growing body of analyses of hydrogen as an energy intermediate by reporting the results from our initial analysis of the potential impacts of the H2@Scale vision by the mid-21 st century for the 48 contiguous U.S. states. Previous estimates have been based on expert elicitation and focused on hydrogen demands. We build upon them, first, by estimating hydrogen’s serviceable consumption potential for possible hydrogen applications and the technical potential for producing hydrogen from various resources. We define the serviceable consumption potential as the quantity of hydrogen that would be consumed to serve the portion of the market that could be captured without considering economics (i.e., if the price of hydrogen were $0/kg over an extended period); thus, it can be considered an upper bound for the size of the market. We define the technical potential as the resource potential constrained by real-world geography and system performance, but not by economics. We then compare the cumulative serviceable consumption potential with the technical potential of a number of possible sources. Second, we estimate economic potential: the quantity of hydrogen at an equilibrium price at which suppliers are willing to sell and consumers are willing to buy the same quantity of hydrogen. We believe this method provides a deeper understanding than was available in the previous analyses. We develop economic potentials for multiple scenarios across various market and technology-advancement assumptions.

08 HYDROGEN↗

Welded and Hydrogen Charged Zircaloy-4 and Welded Stainless Steel 347 Property Data and Microstructures for the Target Solution Vessel and Support Lines of SHINE

SHINE Medical Technologies is developing a facility that will produce molybdenum-99 for medical isotope procedures. As part of this facility, a metal tank, called the target solution vessel (TSV), will hold a uranyl sulfite solution. The TSV must withstand neutron irradiation near room temperature and potential corrosion from the aqueous solution. The large cylindrical TSV will be fabricated by welding and will have numerous pipes and connections also attached via welding. The vessel and those pipes will have varying thicknesses so both gas tungsten arc welding (GTAW) and flux-cored arc welding (FCAW) will be necessary to fabricate different sections. The original material under consideration for the TSV was Zircaloy-4 (Zry-4), but it is known to form a hydride which can degrade its mechanical properties. Thus, some investigation here is focused on the effects of hydrogen uptake in the Zry-4. Additionally, the alternative material being considered is AISI 347, a stainless steel. Both materials have little existing data for their neutron irradiation behavior below 100°C, and both have open questions on the weld behavior under neutron irradiation. Testing has focused on characterizing their weld properties with tensile tests and performing neutron irradiation of samples in the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL).

07 ISOTOPE AND RADIATION SOURCES↗

First Ever Field Pilot on Alaska's North Slope to Validate the Use of Polymer Floods for Heavy Oil EOR a.k.a Alaska North Slope Field Laboratory (ANSFL)

Alaska’s high viscosity oil resources that range between 20–30+ billion barrels represent about a third of known North Slope original oil in place (OOIP). These resources are primarily concentrated in the Schrader Bluff formation (also called West Sak on the Western North Slope) and Ugnu reservoirs and are categorized as “viscous oils” and “heavy oils” owing to their in-situ viscosities between 5–10,000 cP and up to a million+ cP respectively. The viscous oil deposits are relatively deeper (2,000 – 5,000 ft), whereas the heavy oils are somewhat shallower (2,000 – 4,000 ft). The typically shallow depths and the proximity to the continuous permafrost results in relatively lower formation temperatures and pressures, and consequently higher viscosities. The vertical depth vs. viscosity delineated in Paskvan et al. (2016) differentiates the viscous and heavy oils. As depicted in Paskvan et al. (2016), currently the main focus (referred to as “developing”) is on the viscous oils in the Schrader Bluff formation in the Milne Point Unit (MPU). Notwithstanding this Alaska North Slope (ANS) specific categorization, we use the industry adopted, all-inclusive term “heavy oil” for all high viscosity oils. Resource characterization and additional details can be found in topical publications of Paskvan et al. (2016) and Targac et al. (2005). Despite the vast resource base, the development pace, vis-à-vis the production of heavy oils has been very slow and limited due to multiple factors such as cost, logistics, challenging arctic environment, poor waterflood sweep efficiency due to mobility contrasts, and significantly high minimum miscibility pressures (MMP). Most importantly, typical or standard thermal methods that are commonplace elsewhere (Canada, California) are inapplicable due to the continuous permafrost. As a consequence, cumulative production of heavy and viscous oils is a little over 1% of OOIP slope wide and currently, there is hardly any production from Ugnu. However, on a broader level, these unfavorable factors are outweighed by the fact that (1) these resources, within the established infrastructure, are too large to ignore because of their strategic importance to the Nation and the State of Alaska and (2) Prudhoe Bay type diluent crude oil is still available for heavy oil transport through the Trans Alaska Pipeline System (TAPS). Similarly, from a reservoir standpoint, the following factors also are important offsets: (1) favorable rock characteristics of Schrader Bluff; (2) the promise demonstrated by the initial scoping studies (Seright 2010, 2011) suggesting significant increase of heavy oil recovery using polymer flooding; (3) successful field implementation in Canada, China and elsewhere in the world, and (4) availability of the existing pairs of horizontal injector-producer in Schrader Bluff The foregoing was recognized as the best readily available opportunity for significant investment by the US Department of Energy and the field operator Hilcorp Alaska LLC to conduct the first ever field scale experiment to test the polymer flooding technology to unlock the vast heavy oil resources on ANS. With this primary goal in mind, the research team embarked on a ~4.5 years long project that focused on the field polymer pilot complemented by supporting laboratory and simulation studies. As documented in this final report, over the course of the project, many lessons have been learned and valuable field and supporting laboratory data has been collected, which also is complemented by numerical reservoir simulations. We have been able to establish the injectivity of polymer solution, evidence of significant reduction in the water cut of previously waterflooded pattern, effective propagation of a hydrolyzed polyacrylamide (HPAM), benefits of low salinity water, provide practical guidance on handling of produced fluids containing breakthrough polymer, fit-for-purpose forecast-worthy history matched simulation model, polymer EOR benefit of 700-1000 bopd over waterflood, and most importantly a low polymer utilization factor of ~1.7 lb/stb. In summary this project is deemed as a scientific, technical and economic success, having met all objectives, fulfilled deliverables and within budget, providing impetus to apply polymer EOR throughout the Milne Point Field paving the way for even heavier viscosity oils in the Ugnu area, eventually extending the economic life of TAPS.

02 PETROLEUM↗

Binary Control Pulse Optimization for Quantum Systems

Quantum control aims to manipulate quantum systems toward specific quantum states or desired operations. Designing highly accurate and effective control steps is vitally important to various quantum applications, including energy minimization and circuit compilation. In this paper we focus on discrete binary quantum control problems and apply different optimization algorithms and techniques to improve computational efficiency and solution quality. Specifically, we develop a generic model and extend it in several ways. We introduce a squared L 2 -penalty function to handle additional side constraints, to model requirements such as allowing at most one control to be active. We introduce a total variation (TV) regularizer to reduce the number of switches in the control. We modify the popular gradient ascent pulse engineering (GRAPE) algorithm, develop a new alternating direction method of multipliers (ADMM) algorithm to solve the continuous relaxation of the penalized model, and then apply rounding techniques to obtain binary control solutions. We propose a modified trust-region method to further improve the solutions. Our algorithms can obtain high-quality control results, as demonstrated by numerical studies on diverse quantum control examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Two-Scale Neural Networks for Partial Differential Equations with Small Parameters

We propose a two-scale neural network method for solving partial differential equations (PDEs) with small parameters using physics-informed neural networks (PINNs). We directly incorporate the small parameters into the architecture of neural networks. The proposed method enables solving PDEs with small parameters in a simple fashion, without adding Fourier features or other computationally taxing searches of truncation parameters. Various numerical examples demonstrate reasonable accuracy in capturing features of large derivatives in the solutions caused by small parameters.

97 MATHEMATICS AND COMPUTING↗

Radioactive Cleaning Robot

Here at Fermi National Accelerator Laboratory, particle accelerators are constantly running, giving off excess radioactive dust. This dust can be extremely harmful to people and the environment, posing numerous issues. The goal of this project is to conceptualize an automated cleaning solution. With this main goal in mind, the project started with a brainstorming session. The main concept relied on the use of a powder distributing system paired with a vacuum to easily attract and trap the dust. From this main concept, multiple sketches were produced, showing a robot with pieced-together systems. With the initial idea down, each subsystem could be brought further with the inner workings being sketched out. After sketches and proof of concept were polished, the next step was to begin 3D modeling, using NX, to bring life into the sketches. With every day came new ideas, modifying the previous sketches and models as necessary. Soon after the bigger structures were fleshed out, the focus could be shifted onto the less pressing subsystems. Adding a conveyor belt and spinning brush attachment to the powder hopper allowed for more even powder distribution. Along with that, the addition of spring-loaded brushes on the vacuum nozzle ensures constant dust agitation, improving overall functionality. A few other quality of life improvements, such as a removable dust reservoir with a ramp gave a push towards real life functionality. After considering the logistics of this robot functioning in real-life, pressure plates were provided to stop the robot from crashing into its surroundings. Having all of the individual models lead to an overall assembly. Making this robot programmable allowed for this project to meet the design requirements, providing a solution to the previous dilemma.

43 PARTICLE ACCELERATORS↗

Automated Adversary Emulation for Cyber-Physical Systems via Reinforcement Learning

Adversary emulation is an offensive exercise that provides a comprehensive assessment of a system’s resilience against cyber attacks. However, adversary emulation is typically a manual process, making it costly and hard to deploy in cyber-physical systems (CPS) with complex dynamics, vulnerabilities, and operational uncertainties. In this paper, we develop an automated, domain-aware approach to adversary emulation for CPS. We formulate a Markov Decision Process (MDP) model to determine an optimal attack sequence over a hybrid attack graph with cyber (discrete) and physical (continuous) components and related physical dynamics. We apply model-based and model-free reinforcement learning (RL) methods to solve the discrete-continuous MDP in a tractable fashion. As a baseline, we also develop a greedy attack algorithm and compare it with the RL procedures. We summarize our findings through a numerical study on sensor deception attacks in buildings to compare the performance and solution quality of the proposed algorithms.

Bhattacharya, Arnab↗

An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system

The development of an implicit, unconditionally stable, numerical method for solving the Vlasov–Poisson system in one dimension using a phase-space grid is presented. The algorithm uses the Crank–Nicolson discretization scheme and operator splitting allowing for direct solution of the finite difference equations. This method exactly conserves particle number, enstrophy and momentum. A variant of the algorithm which does not use splitting also exactly conserves energy but requires the use of iterative solvers. This algorithm has no dissipation and thus fine-scale variations can lead to oscillations and the production of negative values of the distribution function. We find that overall, the effects of negative values of the distribution function are relatively benign. We consider a variety of test cases that have been used extensively in the literature where numerical results can be compared with analytical solutions or growth rates. We examine higher-order differencing and construct higher-order temporal updates using standard composition methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗