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,045 records · Page 58

Model-independent measurement of the Higgs boson associated production with two jets and decaying to a pair of W bosons in proton-proton collisions at $\sqrt{s}=13$ TeV

A model-independent measurement of the differential production cross section of the Higgs boson decaying into a pair of W bosons, with a final state including two jets produced in association, is presented. In the analysis, events are selected in which the decay products of the two W bosons consist of an electron, a muon, and missing transverse momentum. The model independence of the measurement is maximized by employing a discriminating variable, developed through machine learning, that is agnostic to the signal hypothesis. The analysis is based on proton-proton collision data at $\sqrt{s}=13$ TeV collected with the CMS detector from 2016–2018, corresponding to an integrated luminosity of 138 fb −1 . The production cross section is measured as a function of the difference in azimuthal angle between the two jets. The differential cross section measurements are used to constrain Higgs boson couplings within the standard model effective field theory framework.

Hadron-Hadron Scattering↗

Self-similar Reynolds-averaged mechanical–scalar turbulence models for reshocked Richtmyer–Meshkov instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Richtmyer–Meshkov instability-induced mixing in planar geometry derived in the small Atwood number limit are extended to construct models for reshocked Richtmyer–Meshkov mixing. In this study, the models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ modeled either differentially or algebraically. The three- and four-equation models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing. Mixing layer growth parameters and other physical observables were obtained explicitly as functions of the model coefficients and were used to calibrate the model coefficients. Here, the solutions for the singly shocked Richtmyer–Meshkov case for the mixing layer width and the turbulent fields are used to construct piecewise-continuous generalizations of these quantities for times after reshock. For generality, the post-reshock mixing layer width is not assumed to grow with the same power-law as the pre-reshock width, and an impulsive approximation applied to Rayleigh–Taylor instability growth is used to establish the expression for the post-reshock width. A four-equation model is then used to illustrate the spatiotemporal behavior of the mean and turbulent fields and late-time turbulent equation budgets across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for reshocked Richtmyer–Meshkov instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗

Evaluation of current density measurement methods for high-current, low-beta electron beams

Recent current density measurements with a low-beta electron beam (β=0.5-0.75) are strongly affected by electron scatter and Cherenkov limits. These invasive measurements are heavily influenced by the electron energy and the intercepting material properties: density, index of refraction, and atomic number. We evaluate various measurement methods, including x-ray scintillation and Cherenkov emission. To optimize the x-ray scintillation technique for this energy regime, metal foils of differing atomic number (Z) and thickness were placed upstream of a scintillation screen. We have selected Cherenkov emitters, based on a material refractive index, in order to differentiate the contributions from electron scatter and produced fluorescence, Cherenkov emission, and total internal reflection. The observed distributions for both approaches are verified through MCNP6®. Additionally, the simulations provide insight as to what physical processes dominate the measured distribution. We conclude by determining the limits and optimal measurement range for each measurement technique.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exponential Time Differencing Schemes for Fuel Depletion and Transport in Molten Salt Reactors: Theory and Implementation

A numerical framework for modeling depletion and mass transport in liquid-fueled molten salt reactions is presented based on exponential time differencing. The solution method involves using the finite volume method to transform the system of partial differential equations (PDEs) into a much larger system of ordinary differential equations. The key part of this method involves solving for the exponential of a matrix. We explore six different algorithms to compute the exponential in a series of progression problems that explore physical transport phenomena in molten salt reactors. This framework shows good results for solving linear parabolic PDEs with each of the six matrix exponential algorithms. For large problems, the series solvers such as Padé and Taylor have large run times, which can be mitigated by using the Krylov subspace.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The root cause of disruptive NTMs and paths to stable operation in DIII-D ITER baseline scenario plasmas

Analyses of the DIII-D ITER Baseline Scenario database support that the disruptive m,n=2,1 magnetic islands are pressure gradient driven, non-linear instabilities seeded in a sequence of stochastic transient magnetic perturbations, and that the current profile relaxation does not affect the m,n=2,1 island onset rate. At low torque, these Neoclassical Tearing Modes are most commonly seeded by non-linear 3-wave coupling when the differential rotation between the q=1 & q=2 rational surfaces approaches zero. Lack of statistically significant difference between the current profiles of stable and unstable states, as well as lack of correlation between the tearing mode onset rate and the current profile relaxation both reject causality between the current profile evolution and the 2,1 magnetic island onsets in these plasmas. These support that preserving the differential rotation between the q=1 and q=2 rational surfaces is key to long pulse stable operation in the plasma scenario planned for ITER, while optimization of the current profile within the explored parameter space may lead to much weaker improvements than sustaining the differential rotation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Glauber-theory analysis of nuclear reactions on a 12 C target with variational Monte Carlo wave functions

The application of Glauber theory has been playing an increasingly important role with the study of unstable or exotic nuclei. Its adaptation to medium and high-energy nucleus-nucleus collisions is severely limited because one has to evaluate the matrix elements of multiple-scattering operators. The extraction of physical observables has been done using ‘approximate’ Glauber theory whose validity is hard to evaluate. Here, we perform a full calculation of the matrix elements using Monte Carlo integration and analyze the elastic differential cross sections and the total reaction cross sections for p+¹²C, ⁴,⁶He+¹²C, and ¹²C+¹²C collisions. We use the variational Monte Carlo wave functions for ⁴,⁶He and ¹²C obtained by using realistic two- and three-nucleon potentials. We demonstrate the performance of the Glauber-theory calculations by comparing with available experimental data. We further discuss the accuracy of the conventional approximate methods in the light of the cumulant expansion for Glauber’s phase-shift function.

Horiuchi, W. [Osaka Metropolitan University (Japan↗

Optimization of a Welding Procedure for Making Critical Aluminum Welds on the LBNF Absorber Core Block

he LBNF Absorber consists of thirteen 6061-T6 aluminum core blocks. The core blocks are water cooled with de-ionized (DI) water which becomes radioactive during beam operations. The cooling water flows through gun-drilled channels in the core blocks. A weld quality optimization was performed to produce National Aeronautical Standard (NAS) 1514 Class I quality welds on the aluminum core blocks. This was not successful in all cases. An existing Gas Tungsten Arc Welding (GTAW) Welding Procedure Specification (WPS) was fine tuned to minimize, in most cases, and eliminate detectable tungsten inclusions in the welds. All the weld coupons however, passed welding inspection as per the piping code: ASME B31.3 Normal Fluid Service. Tungsten electrode diameter, type, and manufacturer were varied. Some of the samples were pre-heated and others were not. It was observed that larger diameter electrodes, 5/32 in., with pre-heated joints resulted in welds with the least number of tungsten inclusions. It is hypothesized that thinner electrodes breakdown easily and get lodged into the weld pool during the welding process. This breakdown is further enhanced by the large temperature differential between the un-preheated sample and the hot electrode.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Characterization of soil and postlaunch pad debris from Cape Canaveral launch complex and analysis of soil interaction with aqueous HCl

Soil samples were fractionated and analyzed in order to assess the physical and chemical interactions of entrained soil with solid-rocket exhaust clouds. The sandy soil consisted primarily of quartz (silica) particles, 30 to 500 microns in diameter, and also contained seashell fragments. Differential and cumulative soil-mass size distributions are presented along with mineralogy, elemental compositions, and solution pH histories. About 90 percent of the soil mass consisted of particles 165 microns in diameter. Characteristic reaction times in aqueous HC1 slurries varied from a few minutes to several days, and capacities for reaction under acidic conditions varied from 10 to 40 g HCl/kg soil, depending on particle size. Airborne lifetimes of particles 165 microns are conservatively 30 min, and this major grouping is predicted to represent a small short-term chemical sink for up to 5% of the total HC1. The smaller and more minor fractions, below a 165 micron diameter, may act as giant cloud condensation nuclei over much longer airborne lifetimes. Finally, the demonstrated time dependency of neutralization is a complicating factor; it can influence the ability to deduce in-cloud HCl scavenging with reaction and can affect the accuracy of measured chemical compositions of near-field wet deposition.

Pellett, G. L.↗

Automating the parallel processing of fluid and structural dynamics calculations

The NASA Lewis Research Center is actively involved in the development of expert system technology to assist users in applying parallel processing to computational fluid and structural dynamic analysis. The goal of this effort is to eliminate the necessity for the physical scientist to become a computer scientist in order to effectively use the computer as a research tool. Programming and operating software utilities have previously been developed to solve systems of ordinary nonlinear differential equations on parallel scalar processors. Current efforts are aimed at extending these capabilities to systems of partial differential equations, that describe the complex behavior of fluids and structures within aerospace propulsion systems. This paper presents some important considerations in the redesign, in particular, the need for algorithms and software utilities that can automatically identify data flow patterns in the application program and partition and allocate calculations to the parallel processors. A library-oriented multiprocessing concept for integrating the hardware and software functions is described.

Arpasi, Dale J.↗

Automating the parallel processing of fluid and structural dynamics calculations

The NASA Lewis Research Center is actively involved in the development of expert system technology to assist users in applying parallel processing to computational fluid and structural dynamic analysis. The goal of this effort is to eliminate the necessity for the physical scientist to become a computer scientist in order to effectively use the computer as a research tool. Programming and operating software utilities have previously been developed to solve systems of ordinary nonlinear differential equations on parallel scalar processors. Current efforts are aimed at extending these capabilties to systems of partial differential equations, that describe the complex behavior of fluids and structures within aerospace propulsion systems. This paper presents some important considerations in the redesign, in particular, the need for algorithms and software utilities that can automatically identify data flow patterns in the application program and partition and allocate calculations to the parallel processors. A library-oriented multiprocessing concept for integrating the hardware and software functions is described.

Arpasi, Dale J.↗

Spillover, nonlinearity, and flexible structures

Many systems whose evolution in time is governed by Partial Differential Equations (PDEs) are linearized around a known equilibrium before Computer Aided Control Engineering (CACE) is considered. In this case, there are infinitely many independent vibrational modes, and it is intuitively evident on physical grounds that infinitely many actuators would be needed in order to control all modes. A more precise, general formulation of this grave difficulty (spillover problem) is due to A.V. Balakrishnan. A possible route to circumvention of this difficulty lies in leaving the PDE in its original nonlinear form, and adding the essentially finite dimensional control action prior to linearization. One possibly applicable technique is the Liapunov Schmidt rigorous reduction of singular infinite dimensional implicit function problems to finite dimensional implicit function problems. Omitting details of Banach space rigor, the formalities of this approach are given.

Bass, Robert W.↗

Understanding the Variability of Helium Abundance in the Solar Corona Using Three Fluid Modeling and Ultraviolet Observations

The variability of helium abundance in the solar corona and the solar wind is an important signature of solar activity, solar cycle, and solar wind sources, as well as coronal heating processes. Motivated by recently reported remote-sensing UV imaging observations by Helium Resonance Scattering in the Corona and Heliosphere payload sounding rocket of helium abundance in the inner corona on 2009 September 14 near solar minimum, we present the results of the first three-dimensional three-fluid (electrons, protons, and alpha particles) model of tilted coronal streamer belt and slow solar wind that illustrates the various processes leading to helium abundance differentiation and variability. We find good qualitative agreement between the three-fluid model and the coronal helium abundance variability deduced from UV observations of streamers, providing insight on the effects of the physical processes, such as heating, gravitational settling, and interspecies Coulomb friction in the outflowing solar wind that produce the observed features. The study impacts our understanding of the origins of the slow solar wind.

Leon Ofman↗

Stochastic Modeling in a Multimaterial Continuum Mixture Shock Physics Code

Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigation into the Origin of the RF Phase Noise for the LISA Laser Transmitter

The Laser Interferometer Space Antenna (LISA) mission is a collaborative consortium led by the European Space Agency (ESA)in cooperation with the National Aeronautics and Space Administration (NASA). LISA’s goal is to capture gravitational wave events, such as spacetime distortions caused by massive orbiting bodies. Multiple low noise, single frequency master oscillator power amplifier (MOPA) lasers with a polarization-maintaining Yb-fiber amplifier will serve as the light sources for the LISA observatory. The NASA Goddard Space Flight Center (GSFC)is currently developing the laser transmitters for this project. The LISA lasers will carry radio frequency (RF)sidebands for exchanging reference RF clock information between spacecrafts. Excess phase noise, especially carrier-sideband differential phase noise added to the modulation side band by the phase modulator and the fiber amplifier, becomes a concern. In the gravitational wave frequency ranges, leading sources of this RF “differential” phase noise can be both intrinsic and extrinsic. Examples of intrinsic sources include thermal load on the active fiber, thermo-mechanical induced birefringence, and nonlinearities in the various passive and active fibers. Extrinsic sources can include external temperature fluctuations and infrasonic mechanical strain or vibrations. These noise spectra usually include both white and pink (1/f) thermal noise contributions. This work investigates the physical origins and impact of these sources of RF phase noise, and their characteristic spectra and management will be discussed.

Laser Interferometer Space Antenna↗

Explicit block encodings of boundary value problems for many-body elliptic operators

Simulation of physical systems is one of the most promising use cases of future digital quantum computers. In this work we systematically analyze the quantum circuit complexities of block encoding the discretized elliptic operators that arise extensively in numerical simulations for partial differential equations, including high-dimensional instances for many-body simulations. When restricted to rectangular domains with separable boundary conditions, we provide explicit circuits to block encode the many-body Laplacian with separable periodic, Dirichlet, Neumann, and Robin boundary conditions, using standard discretization techniques from low-order finite difference methods. To obtain high-precision, we introduce a scheme based on periodic extensions to solve Dirichlet and Neumann boundary value problems using a high-order finite difference method, with only a constant increase in total circuit depth and subnormalization factor. We then present a scheme to implement block encodings of differential operators acting on more arbitrary domains, inspired by Cartesian immersed boundary methods. We then block encode the many-body convective operator, which describes interacting particles experiencing a force generated by a pair-wise potential given as an inverse power law of the interparticle distance. This work provides concrete recipes that are readily translated into quantum circuits, with depth logarithmic in the total Hilbert space dimension, that block encode operators arising broadly in applications involving the quantum simulation of quantum and classical many-body mechanics.

Kharazi, Tyler [University of California, Berkeley↗

Inner Magnetospheric Physics

The inner magnetosphere extends from just above the topside ionosphere to approximately 8 RE geocentric distance. Magnetospheric physics is a young science that only started to be recognized as a region with the space observations by Explorer 1 in 1958. The region is mostly populated by ionized gas or plasma from Earth’s ionosphere. Plasma populations are differentiated by their energies primarily. From the least energetic to most are the plasmasphere, ring current, and radiation belts, extending from about 1 eV to 10 MeV in energy and from 1,000s cm-3 down to a few particles per cubic centimeter and less, respectively. The solar wind and solar erupted coronal mass ejections (CMEs) arriving and interacting with Earth’s magnetic field creates a dynamo effect that drives million ampere currents along magnetic field lines that close through the ionosphere. The solar wind dynamo also creates a 100s kV electric field across the magnetosphere that drives convective motion of the plasma within it. The solar wind driven currents compress Earth’s magnetic field on the sunward side and greatly extents the field on the nightside to form the magnetotail. The energy stored in the magnetotail is impulsively released when magnetic field lines there merge, releasing energy into the plasma trapped by the magnetic field. Those plasma become the ring current that loses plasma into the atmosphere to produce the aurora and at the same time ring current plasma can be further energized by wave-particle interactions to become the radiation belts. The presentation will review these topics and a few of the underlying physical processes that are involved in this highly coupled planetary system.

inner magnetosphere↗

Grain formation around carbon stars. 1: Stationary outflow models

Asymptotic giant branch (AGB) stars are known to be sites of dust formation and undergo significant mass loss. The outflow is believed to be driven by radiation pressure on grains and momentum coupling between the grains and gas. While the physics of shell dynamics and grain formation are closely coupled, most previous models of circumstellar shells have treated the problem separately. Studies of shell dynamics typically assume the existence of grains needed to drive the outflow, while most grain formation models assume a constant veolcity wind in which grains form. Furthermore, models of grain formation have relied primarily on classical nucleation theory instead of using a more realistic approach based on chemical kinetics. To model grain formation in carbon-rich AGB stars, we have coupled the kinetic equations governing small cluster growth to moment equations which determine the growth of large particles. Phenomenological models assuming stationary outflow are presented to demonstrate the differences between the classical nucleation approach and the kinetic equation method. It is found that classical nucleation theory predicts nucleation at a lower supersaturation ratio than is predicted by the kinetic equations, resulting in significant differences in grain properties. Coagulation of clusters larger than monomers is unimportant for grain formation in high mass-loss models but becomes more important to grain growth in low mass-loss situations. The properties of the dust grains are altered considerably if differential drift velocities are ignored in modeling grain formation. The effect of stellar temperature, stellar luminosity, and different outflow velocities are investigated. The models indicate that changing the stellar temperature while keeping the stellar luminosity constant has little effect on the physical parameters of the dust shell formed. Increasing the stellar luminosity while keeping the stellar temperature constant results in large differences in grain properties. For small outflow velocities, grains form at lower supersaturation ratios and close to the stellar photosphere, resulting in larger but fewer grains. The reverse is true when grains form under high outflow velocities, i.e., they form at higher supersaturation ratios, farther from the star, and are much smaller but at larger quantities.

Egan, Michael P.↗