Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical physics methods”

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 1,279 records · Page 71

Learning functional priors and posteriors from data and physics

In this work, we develop a new Bayesian framework based on deep neural networks to be able to extrapolate in space-time using historical data and to quantify uncertainties arising from both noisy and gappy data in physical problems. Specifically, the proposed approach has two stages: (1) prior learning and (2) posterior estimation. At the first stage, we employ the physics-informed Generative Adversarial Networks (PI-GAN) to learn a functional prior either from a prescribed function distribution, e.g., Gaussian process, or from historical data and physics. At the second stage, we employ the Hamiltonian Monte Carlo (HMC) method to estimate the posterior in the latent space of PI-GANs. In addition, we use two different approaches to encode the physics: (1) automatic differentiation, used in the physicsinformed neural networks (PINNs) for scenarios with explicitly known partial differential equations (PDEs), and (2) operator regression using the deep operator network (DeepONet) for PDE-agnostic scenarios. We then test the proposed method for (1) meta-learning for one-dimensional regression, and forward/inverse PDE problems (combined with PINNs); (2) PDE-agnostic physical problems (combined with DeepONet), e.g., fractional diffusion as well as saturated stochastic (100-dimensional) flows in heterogeneous porous media; and (3) spatial-temporal regression problems, i.e., inference of a marine riser displacement field using experimental data from the Norwegian Deepwater Programme (NDP). The results demonstrate that the proposed approach can provide accurate predictions as well as uncertainty quantification given very limited scattered and noisy data, since historical data could be available to provide informative priors. In summary, the proposed method is capable of learning flexible functional priors, e.g., both Gaussian and non-Gaussian process, and can be readily extended to big data problems by enabling mini-batch training using stochastic HMC or normalizing flows since the latent space is generally characterized as low dimensional.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Phase Demodulation by Frequency Chirping in Coherence Microwave Photonic Interferometry

This paper presents a signal processing method to demodulate the optical interference phase of cascaded individual optical fiber intrinsic Fabry-Perot interferometric (IFPI) sensors in a coherent microwave-photonic interferometry (CMPI) distributed sensing system. The new method utilizes the chirp effect of electro-optic modulator (EOM) to create a quasi-quadrature optical interference phase shift between two adjacent pulses which correspond to two adjacent reflection points in the time domain. The phase shift can be controlled by adjusting the bias voltage that is applied to the EOM. The interference phase is calculated by elliptically fitting the phase shift. The interference phase change is proportional to the optical path difference (OPD) change of the interferometer, and the sign can be used to differentiate the increase or decrease of the OPD. The method is demonstrated for distributed strain sensing, showing good linearity, high resolution and large dynamic range.

42 ENGINEERING↗

Understanding Material Strength Variabilities and Uncertainties for Component Qualification (2020 LDRD Final Report)

At its core, material specification and materials qualification require an understanding of the uncertainties in relation to a given performance metric. When simulations are used to inform the design and evaluation of these metrics, it is essential to understand the uncertainties in the underlying models and parameters that feed into those simulations. In this work, we have created a statistical framework to address uncertainties in the materials strength modeling used in integrated simulations. The framework is based on the Bayesian methodology, which allows for uncertainties to be updated as new data become available. The results of this work show how model parameters and their uncertainties can be updated as new and different types of strength data are added, as well as methods for determining which future experiment has the greatest potential to reduce the uncertainty. The statistical framework developed in this project will be utilized to assess and propagate the uncertainties in the strength models used in the programmatic simulations and to create a way to formally update those uncertainties as new data become available.

36 MATERIALS SCIENCE↗

Cyber-Physical System Implementation for Manufacturing With Analytics in the Cloud Layer

Effective and efficient modern manufacturing operations require the acceptance and incorporation of the fourth industrial revolution, also known as Industry 4.0. Traditional shop floors are evolving their production into smart factories. To continue this trend, a specific architecture for the cyber-physical system is required, as well as a systematic approach to automate the application of algorithms and transform the acquired data into useful information. This work makes use of an approach that distinguishes three layers that are part of the existing Industry 4.0 paradigm: edge, fog, and cloud. Each of the layers performs computational operations, transforming the data produced in the smart factory into useful information. Trained or untrained methods for data analytics can be incorporated into the architecture. A case study is presented in which a real-time statistical control process algorithm based on control charts was implemented. The algorithm automatically detects changes in the material being processed in a computerized numerical control (CNC) machine. The algorithm implemented in the proposed architecture yielded short response times. The performance was effective since it automatically adapted to the machining of aluminum and then detected when the material was switched to steel. The data were backed up in a database that would allow traceability to the line of g-code that performed the machining.

97 MATHEMATICS AND COMPUTING↗

Three-Receiver Quantum Broadcast Channels: Classical Communication with Quantum Non-unique Decoding

In network communication, it is common in broadcasting scenarios for there to exist a hierarchy among receivers based on information they decode due, for example, to different physical conditions or premium subscriptions. This hierarchy may result in varied information quality, such as higher-quality video for certain receivers. This is modeled mathematically as a degraded message set, indicating a hierarchy between messages to be decoded by different receivers, where the default quality corresponds to a common message intended for all receivers, a higher quality is represented by a message for a smaller subset of receivers, and so forth. We extend these considerations to quantum communication, exploring three-receiver quantum broadcast channels with two- and three-degraded message sets. Our technical tool involves employing quantum non-unique decoding, a technique we develop by utilizing the simultaneous pinching method. Here, we construct one-shot codes for various scenarios and find achievable rate regions relying on various quantum Rényi mutual information error exponents. Our investigation includes a comprehensive study of pinching across tensor product spaces, presenting our findings as the asymptotic counterpart to our one-shot codes. By employing the non-unique decoding, we also establish a simpler proof to Marton’s inner bound for two-receiver quantum broadcast channels without the need for more involved techniques. Additionally, we derive no-go results and demonstrate their tightness in special cases.

Salek, Farzin [Technical University of Munich (Ger↗

Towards a Quantum Algorithm for the Incompressible Nonlinear Navier-Stokes Equations

In this work, we present novel concepts for quantum algorithms to solve transient, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the incompressible Navier-Stokes equations as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts solving nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. We propose a new framework based on matrix product states (MPSs) and matrix product operators (MPOs), in addition to the Krylov subspace methods. For example, the solution variables of the Navier-Stokes equations are represented by MPSs, and the linear and nonlinear terms are processed by MPOs. The time evolution of the operators is attained by a fast-forwarding algorithm using Krylov subspace methods. Furthermore, we discuss various techniques for efficient encoding of MPSs, measurement reduction for MPOs, and use of tensor operations to treat multi-variate, multi-physics characteristics of Navier-Stokes.

Gopalakrishnan Meena, Murali [ORNL] (ORCID:0000000↗

The importance of phonon spectrum measurements for accurate prediction of thermal scattering laws [Slides]

Although approximate phonon spectra (e.g. Debye model) are sufficient to reproduce integral quantities (e.g. total INS cross section) in reasonable agreement with the measured data, more accurate phonon spectra are needed to achieve agreement with the differential data. SANS component is essential for reproducing the total scattering cross section; this is the case for nanodiamonds and other materials where SANS components are present (such as porous materials like nuclear graphite). A format proposal for SANS data in ENDF or GNDS is needed. Therefore, differential (e.g. phonon spectrum) and integral (e.g. total cross section) measurements are necessary in the process of creating an ENDF evaluation. We have successfully developed a method by which ENDF evaluations can be directly compared to the relevant experimentally measured quantities. INS experiments at SNS, ILL, ISIS, and, in the future, ESS, as well as total cross section measurements at the VESUVIO instrument at ISIS, are accessible through user programs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

System and method for deriving a process-based specification

A system and method for deriving a process-based specification for a system is disclosed. The process-based specification is mathematically inferred from a trace-based specification. The trace-based specification is derived from a non-empty set of traces or natural language scenarios. The process-based specification is mathematically equivalent to the trace-based specification. Code is generated, if applicable, from the process-based specification. A process, or phases of a process, using the features disclosed can be reversed and repeated to allow for an interactive development and modification of legacy systems. The process is applicable to any class of system, including, but not limited to, biological and physical systems, electrical and electro-mechanical systems in addition to software, hardware and hybrid hardware-software systems.

Hinchey, Michael Gerard↗

Review of recent activities with MOOSE, an open-source finite element & finite volume multi-fidelity simulation framework

Modeling and simulation are an increasing part of engineering. This is undoubtedly driven by the high costs of constructing experimental facilities, but also enabled by the exponential increase in computing powers over the last decades, which allows computational models to be closer than ever to reality. One of the main drivers for the development of MOOSE is supporting advanced nuclear reactor simulations. A challenging aspect of modeling advanced nuclear reactors is the plurality of physics involved, including neutronics, thermal hydraulics and fuel performance. These physics are all coupled to some extent and are generally solved in a sequential but iterative fashion. The United States (U.S.) national laboratories have been developing MOOSE, an open source multiphysics framework since its inception at the Idaho National Laboratory (INL) in 2008. This framework enables seamless coupling of multiphysics simulations and facilitates the implementation of new physics and material governing laws. It is continuously expanded with novel numerical methods and new pre-implemented physics module. Numerous applications, developed within the Department of Energy (DOE) laboratories, academia, and industry, including outside of nuclear engineering, have been developed to study specialized physics problems. International collaborations are welcome on this open-source modeling and simulation project.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Interface PINNs (I-PINNs): A physics-informed neural networks framework for interface problems

Here, we present a novel physics-informed neural networks (PINNs) framework for modeling interface problems, termed Interface PINNs (I-PINNs). I-PINNs uses different neural networks for any two subdomains separated by a sharp interface such that the neural networks differ only through their activation functions while the other parameters remain identical. The performance of I-PINNs, conventional PINNs, and other existing domain-decomposition PINNs methods such as extended PINNs (XPINNs) and multi-domain PINN (M-PINN) is compared through several one-dimensional, two-dimensional, and three-dimensional benchmark elliptic interface problems. The results demonstrate that I-PINNs provides a root-mean-square-error accuracy, at least two orders of magnitude better than conventional PINNs and XPINNs at approximately one-tenth of the computational cost of conventional PINNs and half the cost of XPINNs. Additionally, while I-PINNs and M-PINN provide comparable accuracies, M-PINN is found to be approximately 50% more expensive.

42 ENGINEERING↗

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗

Compact representation and long-time extrapolation of real-time data for quantum systems using the ESPRIT algorithm

Representing real-time data as a sum of complex exponentials provides a compact form that enables both denoising and extrapolation. As a fully data-driven method, the Estimation of Signal Parameters via Rotational Invariance Techniques (ESPRIT) algorithm is agnostic to the underlying physical equations, making it broadly applicable to various observables and experimental or numerical setups. In this work, we consider applications of the ESPRIT algorithm primarily to extend real-time dynamical data from simulations of quantum systems. We evaluate ESPRIT's performance in the presence of noise and compare it to other extrapolation methods. We demonstrate its ability to extract information from short-time dynamics to reliably predict long-time behavior and determine the minimum time interval required for accurate results. We discuss how this insight can be leveraged in numerical methods that propagate quantum systems in time, and we show how ESPRIT can predict infinite-time values of dynamical observables, offering a purely data-driven approach to characterizing quantum phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

(Invited) Continuum Mathematical Modeling of Water Electrolysis: A Tutorial

Widespread use of hydrogen energy is contingent on the development of reliable and economical sources of hydrogen. Electrolysis from renewably-derived low-carbon electricity is a potentially viable method of hydrogen generation. Prime among the electrolysis technologies are those utilizing ion-conducting polymers (ionomers) including proton-exchange-membrane water electrolyzer (PEMWE). However, these technologies need to exhibit increased efficiency, performance, and durability to become commercially viable. Like most electrochemical devices, PEMWEs involve multiple components (e.g., catalyst, ionomer, transport layers, membrane, plates) and multiple phases, with phenomena occurring across different time and length scales. Furthermore, it is difficult to experimentally probe many of the species and phenomena during operation. Thus, mathematical modeling at the continuum level has been an invaluable aid in exploring, understanding, and optimizing PEMWE cell and components. Furthermore, this is especially true in the highly coupled and complex physics and chemistries that occur with the membrane-electrode assembly (MEA). The physics in a typical volume-averaged non-isothermal model include multiphase transport in porous media, concentrated-solution theory, Ohm's Law, and Butler-Volmer kinetics, and ion, gas, and water transport in the ionomer.

Dizon, Arthur↗

Thermochemical nonequilibrium in atomic hydrogen at elevated temperatures

A numerical study of the nonequilibrium flow of atomic hydrogen in a cascade arc was performed to obtain insight into the physics of the hydrogen cascade arc. A rigorous mathematical model of the flow problem was formulated, incorporating the important nonequilibrium transport phenomena and atomic processes which occur in atomic hydrogen. Realistic boundary conditions, including consideration of the wall electrostatic sheath phenomenon, were included in the model. The governing equations of the asymptotic region of the cascade arc were obtained by writing conservation of mass and energy equations for the electron subgas, an energy conservation equation for heavy particles and an equation of state. Finite-difference operators for variable grid spacing were applied to the governing equations and the resulting system of strongly coupled, stiff equations were solved numerically by the Newton-Raphson method.

Scott, R. K.↗

Assessments of epistemic uncertainty using Gaussian stochastic weight averaging for fluid-flow regression

Here, we use Gaussian stochastic weight averaging (SWAG) to assess the epistemic uncertainty associated with neural-network-based function approximation relevant to fluid flows. SWAG approximates a posterior Gaussian distribution of each weight, given training data, and a constant learning rate. Having access to this distribution, it is able to create multiple models with various combinations of sampled weights, which can be used to obtain ensemble predictions. The average of such an ensemble can be regarded as the 'mean estimation', whereas its standard deviation can be used to construct 'confidence intervals', which enable us to perform uncertainty quantification (UQ) with regard to the training process of neural networks. We utilize representative neural-network-based function approximation tasks for the following cases: (i) a two-dimensional circular-cylinder wake; (ii) the DayMET dataset (maximum daily temperature in North America); (iii) a three-dimensional square-cylinder wake; and (iv) urban flow, to assess the generalizability of the present idea for a wide range of complex datasets. SWAG-based UQ can be applied regardless of the network architecture, and therefore, we demonstrate the applicability of the method for two types of neural networks: (i) global field reconstruction from sparse sensors by combining convolutional neural network (CNN) and multi-layer perceptron (MLP); and (ii) far-field state estimation from sectional data with two-dimensional CNN. We find that SWAG can obtain physically-interpretable confidence-interval estimates from the perspective of epistemic uncertainty. This capability supports its use for a wide range of problems in science and engineering.

97 MATHEMATICS AND COMPUTING↗

An Uncertainty Management Framework for Integrated Gas-Electric Energy Systems

In many parts of the world, electric power systems have seen a significant shift toward generation from renewable energy and natural gas. Because of their ability to flexibly adjust power generation in real time, gas-fired power plants are frequently seen as the perfect partner for variable renewable generation. However, this reliance on gas generation increases interdependence and propagates uncertainty between power grids and gas pipelines and brings coordination and uncertainty management challenges. To address these issues, we propose an uncertainty management framework for uncertain, but bounded gas consumption by gas-fired power plants. The admissible ranges are computed based on a joint optimization problem for the combined gas and electricity networks, which involves chance-constrained scheduling for the electric grid and a novel robust optimization formulation for the natural-gas network. This formulation ensures feasibility of the integrated system with a high probability, while providing a tractable numerical formulation. A key advance with respect to existing methods is that our method is based on a physically accurate, validated model for transient gas pipeline flows. Our case study benchmarks our proposed formulation against methods that ignore how reserve activation impacts the fuel use of gas power plants and only consider predetermined gas consumption. Here, the results demonstrate the importance of considering uncertainty to avoid operating constraint violations and curtailment of gas to the generators.

24 POWER TRANSMISSION AND DISTRIBUTION↗

AI-Constrained Bottom-Up Ecohydrology and Improved Prediction of Seasonal, Interannual, and Decadal Flood and Drought Risks

Focal Areas: (2) Predictive modeling through the use of AI techniques and AI-derived model components; the use of AI and other tools to design a prediction system composed of a hierarchy of models (3)Insight gleaned from complex data (both observed & simulated) using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge-guided AI

54 ENVIRONMENTAL SCIENCES↗

Predicting Catalyst Surface Stability Under Reaction Conditions Using Deep Reinforcement Learning and Machine Learning Potentials

Catalysts are critical for most large-scale energy intensive chemical transformation processes, such as energy storage, liquid fuel production, and the formation of chemical building blocks. The catalyst composition, structure, and morphology impact the performance under reaction conditions and influences properties like activity and selectivity. Many industrial catalysts offer imperfect activity/selectivity, or contain expensive metals. Further, catalysts can deactivate over time as the morphology changes or harsh reactive environments alter the surface structure. Methods to automatically model and predict the kinetics of how catalyst surfaces will restructure would enable engineers to design around these challenges, improve performance, increase longevity. This project investigated a specific type of machine learning model, deep reinforcement learning, machine learning models to act as surrogates for the physical system, and compared the results of these approaches with a specialized high-throughput experimental synthesis and measurement platform.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗