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,189 records · Page 66

ORNL Neutron Cross Section Measurements of 90 Zr NCSP ND-1 Task [Slides]

Luiz Leal's idea to integrate the experimental uncertainties, covariances, and sensitivities in this workflow and centralize the process at one location led to the creation of a cohesive tool for data testing. It has been practiced for more than two decades by the NCS and ND groups at ORNL. The development and use of the ORNL codes TSUNAMI and TSURFER within SCALE made it possible. Today, it has been extended to SAMINT which enables the coupling of differential and integral data evaluation in a continuous-energy framework.

43 PARTICLE ACCELERATORS↗

The heavy particle hazard, what physical data are needed?

The physical data required to evaluate the radiation hazard from heavy galactic cosmic rays to astronauts on extended missions are discussed. The spectral characteristics, nuclear interaction parameters, and track structure of particles are emphasized. The data on the lower energy portion of the differential spectrum of the iron group and nuclear fragmentation in tissue and aluminum are tested, and results are shown.

Curtis, S. B.↗

Boundary-layer equations in generalized curvilinear coordinates

A set of higher-order boundary-layer equations is derived valid for three-dimensional compressible flows. The equations are written in a generalized curvilinear coordinate system, in which the surface coordinates are nonorthogonal; the third axis is restricted to be normal to the surface. Also, higher-order viscous terms which are retained depend on the surface curvature of the body. Thus, the equations are suitable for the calculation of the boundary layer about arbitrary vehicles. As a starting point, the Navier-Stokes equations are derived in a tensorian notation. Then by means of an order-of-magnitude analysis, the boundary-layer equations are developed. To provide an interface between the analytical partial differentiation notation and the compact tensor notation, a brief review of the most essential theorems of the tensor analysis related to the equations of the fluid dynamics is given. Many useful quantities, such as the contravariant and the covariant metrics and the physical velocity components, are written in both notations.

Panaras, Argyris G.↗

Synchronized computational architecture for generalized bilateral control of robot arms

A master six degree of freedom Force Reflecting Hand Controller (FRHC) is available at a master site where a received image displays, in essentially real time, a remote robotic manipulator which is being controlled in the corresponding six degree freedom by command signals which are transmitted to the remote site in accordance with the movement of the FRHC at the master site. Software is user-initiated at the master site in order to establish the basic system conditions, and then a physical movement of the FRHC in Cartesean space is reflected at the master site by six absolute numbers that are sensed, translated and computed as a difference signal relative to the earlier position. The change in position is then transmitted in that differential signal form over a high speed synchronized bilateral communication channel which simultaneously returns robot-sensed response information to the master site as forces applied to the FRHC so that the FRHC reflects the feel of what is taking place at the remote site. A system wide clock rate is selected at a sufficiently high rate that the operator at the master site experiences the Force Reflecting operation in real time.

Szakaly, Zoltan F.↗

Effects of the space environment on space-based radar phased-array antenna; status and preliminary observations (LDEF Experiment A0133)

The overall objective of the experiment was to evaluate the effect of the space environment on components considered for a Space-Based Radar (SBR) Phased-Array Antenna. Of primary interest was a study of the degradation of the polyimide film Kapton (DuPont trademark), the material considered for use in the antenna plane. The most striking result of the experiment was the overall good condition of the Kapton antenna planes and Kapton tensile specimens, despite nearly six years of exposure to the space environment. This was largely attributable to the orientation of the Kapton (parallel and flush on the space end) and the stability of the Long Duration Exposure Facility (LDEF) in orbit. However, weathering of exposed Kapton surfaces was not insignificant. Results on elongation and mechanical properties of the plain and the fiberglass-reinforced Kapton are presented. Reduction in strain to failure of flight-exposed Kapton is attributed to surface defects of these specimens. Physical property testing of the materials to date reveals no significant difference between flight-exposed and control material. The second objective was to investigate the interaction between high-voltage electrodes and typical spacecraft contaminants in simulation of discharge triggering across differentially charged dielectric surfaces (spacecraft charging conditions). Electronic data acquisition and memory systems appeared to operate correctly, but very few discharges were recorded. Induced radioactivity, contamination, impacts, and orientation features of atomic oxygen erosion were observed.

J B Whiteside↗

ORNL Neutron Cross Section Measurements of 90 Zr NCSP ND-1 Task [Slides]

Luiz Leal's idea to integrate the experimental uncertainties, covariances, and sensitivities in the ND workflow and to centralize the process at one location created a cohesive tool for data testing. It has been practiced for more than two decades by the NCS and ND groups at ORNL. The development and use of the ORNL codes TSUNAMI and TSURFER within SCALE made it possible. Today, it has been extended to SAMINT, which enables the coupling of differential and integral data evaluation in a continuous-energy framework.

43 PARTICLE ACCELERATORS↗

NESSi : The N on- E quilibrium S ystems S imulation package

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump–probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green’s functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. Here, we present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green’s functions on the L-shaped Kadanoff–Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green’s functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green’s functions. The library employs a discretization of the Kadanoff–Baym contour into time points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Collaborative Research: Enabling multi-scale studies of magnetic reconnection with interpretable data-driven models

The development of accurate reduced descriptions and improved closures for magnetic reconnection is an important and a long‐standing challenge in plasma physics. The four‐fluid approach, and associated closures, that were investigated have the potential to improve the accuracy of plasma fluid models, capturing physical effects which would otherwise require a kinetic description. If successful, this approach could have an important impact for the modeling of laboratory and space plasmas. The major goals of this project were to develop new machine learning (ML) tools based on sparse and symbolic regression techniques, and to extract interpretable and generalizable reduced models (e.g., in the form of partial differential equations - PDEs) from data generated by first principles plasma simulations. Preserving interpretability of such data‐driven models is key to addressing the long‐standing theoretical and numerical challenges. Prior proof‐of‐principle studies have demonstrated the enormous potential of this approach, by recovering the well‐established hierarchy of plasma equations (from Vlasov to MHD) from data produced by particle‐in‐cell (PIC) simulations. Our goal in this project was to extend and apply these new tools to construct better kinetic closures for magnetic reconnection; to derive better models of particle injection and acceleration by this fundamental plasma process; and to use this understanding to accelerate the development of multi‐scale plasma algorithms. While our immediate focus was on the problem of magnetic reconnection, the tools that were will developed are general and applicable to other areas of plasma physics, and more broadly to many‐body phenomena. We anticipate that the development of these multi‐scale models will have a significant impact across different areas of plasma science, from fusion to space and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Iron-Rain Model for Core Formation on Asteroid 4 Vesta

Asteroid 4 Vesta is differentiated into a crust, mantle, and core, as demonstrated by studies of the eucrite and diogenite meteorites and by data from NASA's Dawn spacecraft. Most models for the differentiation and thermal evolution of Vesta assume that the metal phase completely melts within 20 degrees of the eutectic temperature, well before the onset of silicate melting. In such a model, core formation initially happens by Darcy flow, but this is an inefficient process for liquid metal and solid silicate. However, the likely chemical composition of Vesta, similar to H chondrites with perhaps some CM or CV chondrite, has 13-16 weight percent S. For such compositions, metal-sulfide melting will not be complete until a temperature of at least 1350 degrees Centigrade. The silicate solidus for Vesta's composition is between 1100 and 1150 degrees Centigrade, and thus metal and silicate melting must have substantially overlapped in time on Vesta. In this chemically and physically more likely view of Vesta's evolution, metal sulfide drops will sink by Stokes flow through the partially molten silicate magma ocean in a process that can be envisioned as "iron rain". Measurements of eucrites show that moderately siderophile elements such as Ni, Mo, and W reached chemical equilibrium between the metal and silicate phases, which is an important test for any Vesta differentiation model. The equilibration time is a function of the initial metal grain size, which we take to be 25-45 microns based on recent measurements of H6 chondrites. For these sizes and reasonable silicate magma viscosities, equilibration occurs after a fall distance of just a few meters through the magma ocean. Although metal drops may grow in size by merger with other drops, which increases their settling velocities and decreases the total core formation time, the short equilibration distance ensures that the moderately siderophile elements will reach chemical equilibrium between metal and silicate before metal drop merger becomes important. In this model, there must be at least 30 percent melting of the silicate phase when metal melting is complete, corresponding to a crust thickness of at least 30 kilometers on Vesta, consistent with Dawn gravity observations. Greater degrees of silicate melting and a correspondingly thicker crust are possible if Vesta accreted sufficiently rapidly.

Kiefer, Walter S.↗

First Measurement of Inclusive Muon Neutrino Charged Current Triple Differential Cross Section on Argon

The field of accelerator neutrino experiments is entering an era of precision oscillation measurements where the remaining unknown neutrino measurements will be determined. The upcoming DUNE and Hyper-K experiments aim to determine the neutrino mass hierarchy and degree of Charge-Parity (CP) violation in the neutrino sector, providing potential insight on the matter-antimatter imbalance observed in the universe. However, these experiments require highly accurate measurements, and neutrino cross section modeling uncertainties may limit their capabilities. Cross section measurements at current- generation experiments can aid the development of neutrino interaction models to reduce these uncertainties. This is especially true for measurements of neutrino energy, as it drives neutrino oscillations and is of key importance to oscillation experiments. The MicroBooNE experiment uses a Liquid Argon Time Projection Chamber (LArTPC) to produce neutrino-argon cross sections as one of its physics goals. The MicroBooNE detector’s fully active volume, precision reconstruction, and calorimetry information are leveraged in the Wire-Cell analysis to produce a muon neutrino selection that is 92% pure while maintaining 68% efficiency. A reconstruction chain featuring a fully 3D charge reconstruction and a graph-based particle trajectory fit are used to produce accurate measurements of lepton kinematics as well as visible hadronic energy produced in a neutrino interaction. This thesis presents the first neutrino-argon triple-differential cross section measurement, targeting inclusive charged-current final states. Wiener SVD unfolding is used to produce a measurement over neutrino energy, muon momentum, and muon scattering angle. A series of constrained goodness of fit tests are used to demonstrate the validity of MicroBooNE’s model in describing the distribution of reconstructed kinematics seen in data to ensure the accuracy of unfolding. The validated unfolding to neutrino energy represents a step forward in the field of neutrino cross sections, and demonstrates the capabilities of the LArTPC detector.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Accelerating the Structure Exploration of Diverse Bi–Pt Nanoclusters via Physics‐Informed Machine Learning Potential and Particle Swarm Optimization

Bimetallic Bi–Pt nanoclusters exhibit diverse structural motifs, including core-shell, Janus, and mixed alloy configurations, due to the unique bonding characteristics between Bi and Pt atoms. Using density functional theory refinements from ChIMES physically machine-learned potential and CALYPSO particle swarm optimization global searches, 34 Bi20-Pt20 nanoclusters are systematically classified. The results reveal that Bi atoms predominantly occupy surface sites, driven by charge transfer effects. Cohesive energy trends alone prove insufficient for structure differentiation, necessitating a data-driven approach employing principal component analysis and K-means clustering. Furthermore, vibrational, electronic, and infrared spectral analyses provide additional insights into structure-property relationships. The findings offer an original framework for the automated classification and analysis of bimetallic nanoclusters, enhancing the understanding of their stability and functional properties.

bimetallic nanoparticles↗

Threshold photoproduction of η c and η b using holographic QCD

We discuss the possibility that threshold photoproduction of η c , b may be sensitive to the pseudovector 1 + − glueball exchange. We use the holographic construction to identify the pseudovector glueball with the Kalb-Ramond field, minimally coupled to bulk Dirac fermions. We derive the holographic C-odd form factor and its respective charge radius. Using the pertinent Witten diagrams, we derive and analyze the differential photoproduction cross section for η c , b in the threshold regime, including the interference from the dual bulk photon exchange with manifest vector dominance. The possibility of measuring this process at current and future electron facilities is discussed. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Hybrid (PDE+ML) models in the context of land ice modeling

Focal Area(s): Primary focal areas are: predictive modeling through the use of AI-derived model components; advanced methods including network design/optimization/deep learning. Science Challenge: The atmospheric, ocean and ice dynamics components of the Energy Exascale Earth System Model (E3SM) are governed by Partial Differential Equations (PDEs) and significant efforts have been made during the last decades to develop such computational models. Here we propose to fundamentally improve these PDE-based codes by enhancing them with Machine Learning (ML) sub-models for complex, poorly understood physical processes in the context of ice sheet modeling. We propose to train these models with a novel approach that allows the assimilation of the different sources of data available (direct/indirect observations and possibly simulation data), improving on existing simplified models. We also highlight computational challenges originating from the coexistence of PDE-based and ML-based models.

54 ENVIRONMENTAL SCIENCES↗

Observations of charging dynamics

Plasma data from the University of California, San Diego, Auroral Particles experiments on Applied Technology satellites 5 and 6 are used to investigate the dynamics of natural charging events. Both eclipse and daylight charging events are considered, and typical responses are illustrated by data from specific events. Two different physical processes are found to be involved in the charging process. One of these is straightforward: the spacecraft structure potential responds rapidly to changes in the environment, typically changing by hundreds of volts in a few seconds. The other process is more subtle: differential charging and potential barrier formation precede structural charging and determine the time scale; typically, tens of minutes are required for the potential to change by several hundreds of volts. The latter process is found to be predominately responsible for daylight charging on both spacecraft.

Olsen, R. C.↗

Electro-Thermal Simulation Studies of SiC Junction Diodes Containing Screw Dislocations Under High Reverse-Bias Operation

The objective of this work was to conduct a modeling study of SiC P-N junction diodes operating under high reverse biased conditions. Analytical models and numerical simulation capabilities were to be developed for self-consistent electro-thermal analysis of the diode current-voltage (I-V) characteristics. Data from GRC indicate that screw dislocations are unavoidable in large area SiC devices, and lead to changes in the SiC diode electrical response characteristics under high field conditions. For example, device instability and failures linked to internal current filamentation have been observed. The physical origin of these processes is not well understood, and quantitative projections of the electrical behavior under high field and temperature conditions are lacking. Thermal calculations for SiC devices have not been reported in the literature either. So estimates or projections of peak device temperatures and power limitations do not exist. This numerical study and simulation analysis was aimed at resolving some of the above issues. The following tasks were successfully accomplished: (1) Development of physically based models using one- and two-dimensional drift-diffusion theory for the transport behavior and I-V characteristics; (2) One- and two-dimensional heat flow to account for internal device heating. This led to calculations of the internal temperature profiles, which in turn, were used to update the electrical transport parameters for a self-consistent analysis. The temperature profiles and the peak values were thus obtainable for a given device operating condition; (3) Inclusion of traps assumed to model the presence of internal screw dislocations running along the longitudinal direction; (4) Predictions of the operating characteristics with and without heating as a function of applied bias with and without traps. Both one and two-dimensional cases were implemented; (5) Assessment of device stability based on the operating characteristics. The presence of internal non-uniformities, particularly filamentary structures, was probed and demonstrated; (6) Cause and physical origins of filamentary behavior and unstable I-V characteristics were made transparent; (7) It was demonstrated that diodes containing defects would be more prone to thermal breakdown associated with the temperature dependent decrease in the thermal conductivity; and (8) Finally, negative differential resistance (S-shaped NDR) which can potential lead to device instability and filamentary behavior was shown to occur for diodes containing a line of defects such as could be associated with a screw dislocation line.

Joshi, R. P.↗

Physics-based Tests to Identify the Accuracy of Solar Wind Ion Measurements: A Case Study with the Wind Faraday Cups

We present techniques for comparing measurements of velocity, temperature, and density with constraints imposed by the plasma physics of magnetized bi-Maxwellian ions. Deviations from these physics-based constraints are interpreted as arising from measurement errors. Two million ion spectra from the Solar Wind Experiment Faraday Cup instruments on the Wind spacecraft are used as a case study. The accuracy of velocity measurements is determined by the fact that differential flow between hydrogen and helium should be aligned with the ambient magnetic field. Modeling the breakdown of field alignment suggests velocity uncertainties are less than 0.16% in magnitude and 3deg in direction. Temperature uncertainty is found by examining the distribution of observed temperature anisotropies in high-beta solar wind intervals where the firehose, mirror, and cyclotron microinstabilities should drive the distribution to isotropy. The presence of a finite anisotropy at high beta suggests overall temperature uncertainties of 8%. Hydrogen and helium number densities are compared with the electron density inferred from observations of the local electron plasma frequency as a function of solar wind speed and year. We find that after accounting for the contribution of minor ions, the results are consistent with a systematic offset between the two instruments of 34%. The temperature and density methods are sensitive to non-Maxwellian features such as heat flux and proton beams and as a result are more suited to slow solar wind where these features are rare. These procedures are of general use in identifying the accuracy of observations from any solar wind ion instrument.

Kasper, J. C.↗

Redox Variations in Early Solar System Materials and Implications for Late Stage Planetary Accretion and Planet Formation

Oxygen fugacity plays an important role in determining the detailed physical and chemical aspects of planets and their building blocks. Basic chemical properties such as the amount of oxidized Fe in a mantle (as FeO), the nature of alloying elements in the core (S, C, H, O, Si), and the solubility of various volatile elements in the silicate and metallic portions of embryos and planets can influence physical properties such as the size of the core, the liquidus and solidus of the mantle and core, and the speciation of volatile compounds contributing to atmospheres. This paper will provide an overview of the range of fO2 variation observed in primitive and differentiated materials that may have participated in accretion (cosmic dust, Star-dust and meteorites), a comparison to observations of planetary fO2 (Mercury, Mars and Earth), and a discus-sion of timing of variation of fO2 within both early and later accreted materials. This overview is meant to promote discussion and interaction between students of these two stages of planet formation to identify areas where more work is needed.

Righter, K.↗

Physics-driven learning of Wasserstein GAN for density reconstruction in dynamic tomography

Object density reconstruction from projections containing scattered radiation and noise is of critical importance in many applications. Existing scatter correction and density reconstruction methods may not provide the high accuracy needed in many applications and can break down in the presence of unmodeled or anomalous scatter and other experimental artifacts. Incorporating machine-learning models could prove beneficial for accurate density reconstruction, particularly in dynamic imaging, where the time evolution of the density fields could be captured by partial differential equations or by learning from hydrodynamics simulations. In this work, we demonstrate the ability of learned deep neural networks to perform artifact removal in noisy density reconstructions, where the noise is imperfectly characterized. Here, we use a Wasserstein generative adversarial network (WGAN), where the generator serves as a denoiser that removes artifacts in densities obtained from traditional reconstruction algorithms. We train the networks from large density time-series datasets, with noise simulated according to parametric random distributions that may mimic noise in experiments. The WGAN is trained with noisy density frames as generator inputs, to match the generator outputs to the distribution of clean densities (time series) from simulations. A supervised loss is also included in the training, which leads to an improved density restoration performance. In addition, we employ physics-based constraints such as mass conservation during the network training and application to further enable highly accurate density reconstructions. Our preliminary numerical results show that the models trained in our frameworks can remove significant portions of unknown noise in density time-series data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗