Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differentiable physics”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Surface Equilibration Mechanism Controls the Stability of a Model Codeposited Glass Mixture of Organic Semiconductors

While previous work has identified the conditions for preparing ultrastable single-component organic glasses by physical vapor deposition (PVD), little is known about the stability of codeposited mixtures. Here, we prepared binary PVD glasses of organic semiconductors, TPD (N,N'-Bis(3-methylphenyl)-N,N'-diphenylbenzidine) and m-MTDATA (4,4',4"-Tris[phenyl(m-tolyl)amino]triphenylamine), with a 50:50 mass concentration over a wide range of substrate temperatures (T sub ). The enthalpy and kinetic stability are evaluated with differential scanning calorimetry and spectroscopic ellipsometry. Binary organic semiconductor glasses with exceptional thermodynamic and kinetic stability comparable to the most stable single-component organic glasses are obtained when deposited at T sub = 0.78–0.90T g (where T g is the conventional glass transition temperature). When deposited at 0.94T g , the enthalpy of the m-MTDATA/TPD glass equals that expected for the equilibrium liquid at that temperature. Thus, the surface equilibration mechanism previously advanced for single-component PVD glasses is also applicable for these codeposited glasses. Furthermore, these results provide an avenue for designing high-performance organic electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Uncertainty quantification in scientific machine learning: Methods, metrics, and comparisons

Neural networks (NNs) are currently changing the computational paradigm on how to combine data with mathematical laws in physics and engineering in a profound way, tackling challenging inverse and ill-posed problems not solvable with traditional methods. However, quantifying errors and uncertainties in NN-based inference is more complicated than in traditional methods. This is because in addition to aleatoric uncertainty associated with noisy data, there is also uncertainty due to limited data, but also due to NN hyperparameters, overparametrization, optimization and sampling errors as well as model misspecification. Although there are some recent works on uncertainty quantification (UQ) in NNs, there is no systematic investigation of suitable methods towards quantifying the total uncertainty effectively and efficiently even for function approximation, and there is even less work on solving partial differential equations and learning operator mappings between infinite-dimensional function spaces using NNs. In this work, we present a comprehensive framework that includes uncertainty modeling, new and existing solution methods, as well as evaluation metrics and post-hoc improvement approaches. Further, to demonstrate the applicability and reliability of our framework, we present an extensive comparative study in which various methods are tested on prototype problems, including problems with mixed input-output data, and stochastic problems in high dimensions. In the Appendix, we include a comprehensive description of all the UQ methods employed. Further, to help facilitate the deployment of UQ in Scientific Machine Learning research and practice, we present and develop in [1] an open-source Python library (github.com/Crunch-UQ4MI/neuraluq), termed NeuralUQ, that is accompanied by an educational tutorial and additional computational experiments.

11 physics-informed neural networks↗
DOE OSTI2021

Modeling atmospheric data and identifying dynamics Temporal data-driven modeling of air pollutants

Atmospheric modelling has recently experienced a surge with the advent of deep learning. Most of these models, however, predict concentrations of pollutants following a data-driven approach in which the physical laws that govern their behaviors and relationships remain hidden. Furthermore, with the aid of real-world air quality data collected hourly in different stations throughout Madrid, we present a case study using a series of data-driven techniques with the following goals: (1) Find systems of ordinary differential equations that model the concentration of pollutants and their changes over time; (2) assess the performance and limitations of our model using stability analysis; (3) reconstruct the time series of chemical pollutants not measured in certain stations using delay coordinate embedding results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

An Open‐Source, Physics‐Based, Tropical Cyclone Downscaling Model With Intensity‐Dependent Steering

Abstract An open‐source, physics‐based tropical cyclone (TC) downscaling model is developed, in order to generate a large climatology of TCs. The model is composed of three primary components: (a) a random seeding process that determines genesis, (b) an intensity‐dependent beta‐advection model that determines the track, and (c) a non‐linear differential equation set that determines the intensification rate. The model is entirely forced by the large‐scale environment. Downscaling ERA5 reanalysis data shows that the model is generally able to reproduce observed TC climatology, such as the global seasonal cycle, genesis locations, track density, and lifetime maximum intensity distributions. Inter‐annual variability in TC count and power‐dissipation is also well captured, on both basin‐wide and global scales. Regional TC hazard estimated by this model is also analyzed using return period maps and curves. In particular, the model is able to reasonably capture the observed return period curves of landfall intensity in various sub‐basins around the globe. The incorporation of an intensity‐dependent steering flow is shown to lead to regionally dependent changes in power dissipation and return periods. Advantages and disadvantages of this model, compared to other downscaling models, are also discussed.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Severe, short-term sleep restriction reduces gut microbiota community richness but does not alter intestinal permeability in healthy young men

Sleep restriction alters gut microbiota composition and intestinal barrier function in rodents, but whether similar effects occur in humans is unclear. This study aimed to determine the effects of severe, short-term sleep restriction on gut microbiota composition and intestinal permeability in healthy adults. Fecal microbiota composition, measured by 16S rRNA sequencing, and intestinal permeability were measured in 19 healthy men (mean ± SD; BMI 24.4 ± 2.3 kg/m 2 , 20 ± 2 years) undergoing three consecutive nights of adequate sleep (AS; 7–9 h sleep/night) and restricted sleep (SR; 2 h sleep/night) in random order with controlled diet and physical activity. α-diversity measured by amplicon sequencing variant (ASV) richness was 21% lower during SR compared to AS (P = 0.03), but α-diversity measured by Shannon and Simpson indexes did not differ between conditions. Relative abundance of a single ASV within the family Ruminococcaceae was the only differentially abundant taxon (q = 0.20). No between-condition differences in intestinal permeability or β-diversity were observed. Findings indicated that severe, short-term sleep restriction reduced richness of the gut microbiota but otherwise minimally impacted community composition and did not affect intestinal permeability in healthy young men.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Voronoi cell analysis: The shapes of particle systems

Many physical systems can be studied as collections of particles embedded in space, often evolving in time. Natural questions arise concerning how to characterize these arrangements—are they ordered or disordered? If they are ordered, how are they ordered and what kinds of defects do they possess? Voronoi tessellations, originally introduced to study problems in pure mathematics, have become a powerful and versatile tool for analyzing countless problems in pure and applied physics. Here we explain the basics of Voronoi tessellations and the shapes that they produce and describe how they can be used to characterize many physical systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
NASA NTRS1992

Wave equations and computational models for sonic boom propagation through a turbulent atmosphere

The improved simulation of sonic boom propagation through the real atmosphere requires greater understanding of how the transient acoustic pulses popularly termed sonic booms are affected by atmospheric turbulence. A nonlinear partial differential equation that can be used to simulate the effects of smaller-scale atmospheric turbulence on sonic boom waveforms is described. The equation is first order in the time derivative and involves an extension of geometrical acoustics to include diffraction phenomena. Various terms in the equation are explained in physical terms. Such terms include those representing convection at the wave speed, diffraction, molecular relaxation, classical dissipation, and nonlinear steepening. The atmospheric turbulence enters through an effective sound speed, which varies with all three spatial coordinates, and which is the sum of the local sound speed and the component of the turbulent flow velocity projected along a central ray that connects the aircraft trajectory with the listener.

Pierce, Allan D.↗
DOE OSTI2020

The Locus Algorithm III: A Grid Computing system to generate catalogues of optimised pointings for Differential Photometry

This paper discusses the hardware and software components of the Grid Computing system used to implement the Locus Algorithm to identify optimum pointings for differential photometry of 61,662,376 stars and 23,799 quasars. The scale of the data, together with initial operational assessments demanded a High Performance Computing (HPC) system to complete the data analysis. Grid computing was chosen as the HPC solution as the optimum choice available within this project. The physical and logical structure of the National Grid computing Infrastructure informed the approach that was taken. That approach was one of layered separation of the different project components to enable maximum flexibility and extensibility.

Creaner, Oisń↗
DOE OSTI2023

Learning heterogeneous reaction kinetics from X-ray videos pixel by pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries and electrocatalysts. Experimental characterizations of such materials by operando microscopy produce rich image datasets, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

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

Lattice Disorder and Oxygen Migration Pathways in Pyrochlore and Defect-Fluorite Oxides

Atomic-scale disorder plays an important role in the chemical and physical properties of oxide materials. The structural flexibility of pyrochlore-type oxides allows for crystal-chemical engineering of these properties. Compositional modification can push pyrochlore oxides toward a disordered defect-fluorite structure with anion Frenkel pair defects that facilitate oxygen migration. The local structure of the long-range average cubic defect-fluorite was recently claimed to consist of randomly arranged orthorhombic weberite-type domains. Here, we show, using low-temperature neutron total-scattering experiments, that this is not the case for Zr-rich defect-fluorites. By analyzing data from the pyrochlore/defect-fluorite Y 2 Sn 2–x Zr x O 7 series using a combination of neutron pair distribution function and big-box modelling, we have differentiated and quantified the relationship between anion sub-lattice disorder and Frenkel defects. These details directly influence the energy landscape for oxygen migration and are crucial for simulations and design of new materials with improved properties.

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

Unfolding quantum computer readout noise

Abstract In the current era of noisy intermediate-scale quantum computers, noisy qubits can result in biased results for early quantum algorithm applications. This is a significant challenge for interpreting results from quantum computer simulations for quantum chemistry, nuclear physics, high energy physics (HEP), and other emerging scientific applications. An important class of qubit errors are readout errors. The most basic method to correct readout errors is matrix inversion, using a response matrix built from simple operations to probe the rate of transitions from known initial quantum states to readout outcomes. One challenge with inverting matrices with large off-diagonal components is that the results are sensitive to statistical fluctuations. This challenge is familiar to HEP, where prior-independent regularized matrix inversion techniques (“unfolding”) have been developed for years to correct for acceptance and detector effects, when performing differential cross section measurements. We study one such method, known as iterative Bayesian unfolding, as a potential tool for correcting readout errors from universal gate-based quantum computers. This method is shown to avoid pathologies from commonly used matrix inversion and least squares methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Kinetic mixing, dark photons and extra dimensions. Part III. Brane localized dark matter

Extra dimensions have proven to be a very useful tool in constructing new physics models. In earlier work, we began investigating toy models for the 5-D analog of the kinetic mixing/vector portal scenario where the interactions of dark matter, taken to be, e.g., a complex scalar, with the brane-localized fields of the Standard Model (SM) are mediated by a massive U(1) D dark photon living in the bulk. These models were shown to have many novel features differentiating them from their 4-D analogs and which, in several cases, avoided some well-known 4-D model building constraints. However, these gains were obtained at the cost of the introduction of a fair amount of model complexity, e.g., dark matter Kaluza-Klein excitations. In the present paper, we consider an alternative setup wherein the dark matter and the dark Higgs, responsible for U(1) D breaking, are both localized to the ‘dark’ brane at the opposite end of the 5-D interval from where the SM fields are located with only the dark photon now being a 5-D field. The phenomenology of such a setup is explored for both flat and warped extra dimensions and compared to the previous more complex models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Why Boron Nitride is such a Selective Catalyst for the Oxidative Dehydrogenation of Propane

Boron–containing materials, and in particular boron nitride, have recently been identified as highly selective catalysts for the oxidative dehydrogenation of alkanes such as propane. To date, no mechanism exists that can explain both the unprecedented selectivity, the observed surface oxyfunctionalization, and the peculiar kinetic features of this reaction. Here, we combine catalytic activity measurements with quantum chemical calculations to put forward a bold new hypothesis. We argue that the remarkable product distribution can be rationalized by a combination of surface–mediated formation of radicals over metastable sites, and their sequential propagation in the gas phase. Based on known radical propagation steps, we quantitatively describe the oxygen pressure–dependent relative formation of the main product propylene and by–product ethylene. Free radical intermediates most likely differentiate this catalytic system from less selective vanadium–based catalysts.

03 NATURAL GAS↗
DOE OSTI2021

Asymmetric Fraunhofer spectra in a topological insulator-based Josephson junction

Josephson junctions with topological insulators as their weak link (S–TI–S junctions) are predicted to host Majorana fermions, which are key to creating qubits for topologically protected quantum computing. But the details of the S–TI–S current-phase relation and its interplay with magnetic fields are not well understood. We fabricate a Bi 2 Se 3 junction with NbTi leads and measure the Fraunhofer patterns of the junction with applied in-plane fields. We observe that asymmetric Fraunhofer patterns appear in the resistance maps of B z vs B x,y , with aperiodic node spacings. These asymmetric patterns appear even at zero parallel field and for temperatures up to 1 K. The anomalous features are compared to asymmetric Fraunhofer patterns expected for finite Cooper pair momentum shifts as well as geometric effects. We show that the geometric effects can dominate, independent of in-plane field magnitude. Furthermore, these results are important for differentiating geometrical phase shifts from those caused by Cooper pair momentum shifting, Majorana mode signatures, or other unconventional superconducting behavior.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Ab initio nucleon-nucleus elastic scattering with chiral effective field theory uncertainties

Effective interactions for nucleon-nucleus elastic scattering from first principles require the use of the same nucleon-nucleon interaction in the structure and reaction calculations and a consistent treatment of the relevant operators at each order. Systematic investigations of the effect of truncation uncertainties of chiral nucleon-nucleon ( N N ) forces have been carried out for scattering observables in the two- and three-nucleons system as well as for bound-state properties of light nuclei. Here we extend this type of study to proton and neutron elastic scattering for 16 O and 12 C. Using the frameworks of the spectator expansion of multiple scattering theory as well as the no-core shell model, we employ one specific chiral interaction from the LENPIC collaboration and consistently calculate the leading-order effective nucleon-nucleus interaction up to the third chiral order, from which we extract elastic scattering observables. Then we apply pointwise as well as correlated uncertainty quantification for the estimation of the chiral truncation error. We calculate and analyze proton elastic scattering observables for 16 O and 12 C as well as neutron elastic scattering observables for 12 C between 65 and 185 MeV projectile kinetic energy. We find in all cases qualitatively similar results for the chiral truncation uncertainties as in few-body systems and assess them using similar diagnostic tools. The order-by-order convergence of the scattering observables for 16 O and 12 C is very reasonable around 100 MeV, while for higher energies the chiral expansion parameter becomes too large for convergence. Comparing proton and neutron differential cross sections reveals that their description is comparable up to around 100 MeV. Here, we also find a nearly perfect correlation between the differential cross section for neutron scattering and the N N Wolfenstein amplitudes for small momentum transfers. The diagnostic tools for studying order-by-order convergence in observables in few-body systems can be employed for observables in nucleon-nucleus scattering with only minor modifications provided the momentum scale in the problems is not too large. We also find that the chiral N N interaction on which our study is based gives a very good description of differential cross sections for 16 O and 12 C as low as 65-MeV projectile energy, particularly in the forward direction. In addition, the very forward direction of the neutron differential cross section mirrors the behavior of the N N interaction amazingly well.

6 ≤ A ≤ 19↗
DOE OSTI2025

Study of the light scalar a 0 ( 980 ) through the decay D 0 → a 0 ( 980 ) − e + ν e with a 0 ( 980 ) − → η π −

Using 7.93 fb − 1 of e + e − collision data collected at a center-of-mass energy of 3.773 GeV with the BESIII detector, we present an analysis of the decay D 0 → η π − e + ν e . The branching fraction of the decay D 0 → a 0 ( 980 ) − e + ν e with a 0 ( 980 ) − → η π − is measured to be ( 0.86 ± 0.1 7 stat ± 0.0 5 syst ) × 10 − 4 . The decay dynamics of this process is studied with a single-pole parametrization of the hadronic form factor and the Flatté formula describing the a 0 ( 980 ) line shape in the differential decay rate. The product of the form factor f + a 0 ( 0 ) and the Cabibbo-Kobayashi-Maskawa matrix element | V c d | is determined for the first time with the result f + a 0 ( 0 ) | V c d | = 0.126 ± 0.01 3 stat ± 0.00 3 syst . Published by the American Physical Society 2025

Ablikim, M.↗
NASA NTRS1975

Special analysis of community annoyance with aircraft noise reported by residents in the vicinity of JFK Airport, 1972

During the summer of 1972, about 1500 residents were interviewed twice in 11 communities near JFK airport. Detailed aircraft operations reports were also collected for this period, and an effort has been made to analyze recorded human response data in relation to a number of physical exposure parameters. A series of exposure indexes, based on an arithmetic integration of aircraft operations, were correlated with summated aircraft noise annoyance responses. None of these correlations were as good as the CNR index which assumes a logrithmetic integration of numbers of aircraft exposures and includes a day-night differential weighting of 10:1. There were substantial variations in average annoyance responses among communities with similar CNR exposures, substantiating previous findings that attitudinal and other personal variables also play an important role in determining annoyance differences.

Borsky, P. N.↗