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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Essentially nonoscillatory postprocessing filtering methods

High order accurate centered flux approximations used in the computation of numerical solutions to nonlinear partial differential equations produce large oscillations in regions of sharp transitions. Here, we present a new class of filtering methods denoted by Essentially Nonoscillatory Least Squares (ENOLS), which constructs an upgraded filtered solution that is close to the physically correct weak solution of the original evolution equation. Our method relies on the evaluation of a least squares polynomial approximation to oscillatory data using a set of points which is determined via the ENO network. Numerical results are given in one and two space dimensions for both scalar and systems of hyperbolic conservation laws. Computational running time, efficiency, and robustness of method are illustrated in various examples such as Riemann initial data for both Burgers' and Euler's equations of gas dynamics. In all standard cases, the filtered solution appears to converge numerically to the correct solution of the original problem. Some interesting results based on nonstandard central difference schemes, which exactly preserve entropy, and have been recently shown generally not to be weakly convergent to a solution of the conservation law, are also obtained using our filters.

Lafon, F.↗
NASA NTRS1994

An Approach for Dynamic Grids

An approach is presented for the generation of two-dimensional, structured, dynamic grids. The grid motion may be due to the motion of the boundaries of the computational domain or to the adaptation of the grid to the transient, physical solution. A time-dependent grid is computed through the time integration of the grid speeds which are computed from a system of grid speed equations. The grid speed equations are derived from the time-differentiation of the grid equations so as to ensure that the dynamic grid maintains the desired qualities of the static grid. The grid equations are the Euler-Lagrange equations derived from a variational statement for the grid. The dynamic grid method is demonstrated for a model problem involving boundary motion, an inviscid flow in a converging-diverging nozzle during startup, and a viscous flow over a flat plate with an impinging shock wave. It is shown that the approach is more accurate for transient flows than an approach in which the grid speeds are computed using a finite difference with respect to time of the grid. However, the approach requires significantly more computational effort.

Slater, John W.↗
NASA NTRS2010

Analysis of Instabilities in Non-Axisymmetric Hypersonic Boundary Layers Over Cones

Hypersonic flows over circular cones constitute one of the most important generic configurations for fundamental aerodynamic and aerothermodynamic studies. In this paper, numerical computations are carried out for Mach 6 flows over a 7-degree half-angle cone with two different flow incidence angles and a compression cone with a large concave curvature. Instability wave and transition-related flow physics are investigated using a series of advanced stability methods ranging from conventional linear stability theory (LST) and a higher-fidelity linear and nonlinear parabolized stability equations (PSE), to the 2D eigenvalue analysis based on partial differential equations. Computed N factor distribution pertinent to various instability mechanisms over the cone surface provides initial assessments of possible transition fronts and a guide to corresponding disturbance characteristics such as frequency and azimuthal wave numbers. It is also shown that strong secondary instability that eventually leads to transition to turbulence can be simulated very efficiently using a combination of advanced stability methods described above.

Li, Fei↗
DOE OSTI2021

Differentiable programming for online training of a neural artificial viscosity function within a staggered grid Lagrangian hydrodynamics scheme

Lagrangian methods to solve the inviscid Euler equations produce numerical oscillations near shock waves. A common approach to reducing these oscillations is to add artificial viscosity (AV) to the discrete equations. The AV term acts as a dissipative mechanism that attenuates oscillations by smearing the shock across a finite number of computational cells. However, AV introduces several control parameters that are not determined by the underlying physical model, and hence, in practice are tuned to the characteristics of a given problem. We seek to improve the standard quadratic-linear AV form by replacing it with a learned neural function that reduces oscillations relative to exact solutions of the Euler equations, resulting in a hybrid numerical-neural hydrodynamic solver. Because AV is an artificial construct that exists solely to improve the numerical properties of a hydrodynamic code, there is no offline ‘viscosity data’ against which a neural network can be trained before inserting into a numerical simulation, thus requiring online training. We achieve this via differentiable programming, i.e. end-to-end backpropagation or adjoint solution through both the neural and differential equation code, using automatic differentiation of the hybrid code in the Julia programming language to calculate the necessary loss function gradients. A novel offline pre-training step accelerates training by initializing the neural network to the default numerical AV scheme, which can be learned rapidly by space-filling sampling over the AV input space. We find that online training over early time steps of simulation is sufficient to learn a neural AV function that reduces numerical oscillations in long-term hydrodynamic shock simulations. These results offer an early proof-of-principle that online differentiable training of hybrid numerical schemes with novel neural network components can improve certain performance aspects existing in purely numerical schemes.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

The physical properties of double degenerate common proper motion binaries

Spectral types and spectrophotometry are presented for 21 double degenerate (DD) common proper motion binaries, along with estimates of their colors, absolute visual and bolometric magnitudes, and cooling ages. The oldest pairs in the sample are 9 x 10 to the 9th yr; the differential cooling ages range from 0.01 to 0.84. The median and mean separations of the DD pairs are 426 and 407 Au, respectively, both apparently smaller than the WD+MS values. The average UVW motions and velocity dispersions are significantly larger than the average velocities and dispersions associated with selected samples of single white dwarfs and MS+WD binaries when the latter are restricted to the same color/Mv range as the DD systems. This may be a result of the dynamical inflation of the velocity dispersion of DD systems due to their extremely ancient total stellar ages.

Sion, Edward M.↗
DOE OSTI2022

Generalized Bayesian Monte Carlo Evaluation [Slides]

Bayesian Monte Carlo is a tool to address imperfect data & models, non-linear models, non-normal PDFs, and integration of differential and integral data. Simple propagation to criticality shows effects of new $\mathcal{L}$($β|z,γ$).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Measurement of Charmonium Production in $p + p$ and $p + d$ Interactions in the Fermilab SeaQuest Experiment

The Fermilab SeaQuest experiment has measured dimuon events from the interactions of 120GeV proton beam on liquid hydrogen and deuterium targets with dimuon mass between 2 and 9GeV. These dimuon events contain both the Drell-Yan process and the charmonium ($J/\psi$ and $\psi^\prime$) production. Unlike the Drell-Yan process which probes the antiquark distributions in the nucleons, the charmonium production is sensitive to both quark and gluon distributions. SeaQuest has extracted the $\sigma^{pd}/2\sigma^{pp}$ ratio as well as the differential cross sections for charmonium production in the kinematic region of $0.4 < x_F < 0.9$. The $\sigma^{pd}/2\sigma^{pp}$ ratio for charmonium production are found to be significantly different from that of the Drell-Yan process. The measured differential cross sections for charmonium production are compared with theoretical calculations using Color Evaporation Model and Non-Relativistic QCD model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Selective Detection of Diglycolamide Diastereomers By Metal-Ligand Cluster Formation

Introduction (120 words): Differentiation of isomers using mass spectrometry is a challenge, as these species tend to produce ions of identical mass-to-charge ratio. Therefore, separation with mass spectrometry must rely on chemical differences between isomeric species. Separation of stereoisomers is especially challenging, as these species often have very similar chemistry. However, recent investigations of the solvent extraction behavior of diastereomers of diglycolamide molecules have shown that different stereoisomers can have significantly different lanthanide extraction efficiencies. These studies showed that the R,S stereoisomer of 2,2'-oxybis(N,N-dioctylpropanamide efficiently extracts lanthanides, while the S,S stereoisomer does not. Thus, gas-phase metal complexation may provide a pathway for differentiating these species using mass spectrometry. Methods (120 words): Europium, samarium, and holmium clusters with the R,S, and the S,S diastereomers of 2,2'-oxybis(N,N-didecylpropanamide (trivially tetradecyl diglycolamide, or mTDDGA) were generated in gas-phase using the electrospray ionization source of a Bruker (Billerica, MA, USA) micrOTOF-Q II quadrupole time-of-flight mass spectrometer. Solutions of aqueous lanthanide nitrate and diglycolamide ligands in 2-propanol were added to 50%:50% by volume acetonitrile:2-propanol to generate the spray solutions. Ligand concentration was maintained at approximately 6 µM, while lanthanide concentration was maintained at either 6 µM or 30 µM . High resolution, high mass accuracy spectra were recorded using the time-of-flight. Mass accuracy was ensured using external calibration with Agilent (Santa Clara, CA, USA) ESI-L Low Concentration tuning mix. Preliminary data (300 words): Electrospray of solutions containing only R,S-mTDDGA or S,S-m-TDDGA yield abundant peaks at m/z= 721.7, corresponding to [mTDDGA+H]+1, m/z = 743.7, corresponding to [mTDDGA+Na]+1, and m/z = 1464.4, corresponding to the sodium-bound dimer, [Na(mTDDGA)2]+1. Spray of solutions containing 6 µM europium (Eu) with 3 µM R,S-mTDDGA produced an abundant cluster ion with an isotopic envelope centered at m/z=772.0, corresponding to [Eu(R,S-mTDDGA)3]3+. Europium-containing cluster ions with isotopic envelopes apexing at m/z = 1012.3 and 1189.0 correspond to the [Eu(R,S-mTDDGA)4]3+ and [Eu(R,S-mTDDGA)5]3+ cluster ions, respectively. Protonated and sodiated R,S-mTDDGA are present as well, along with the sodium-bound dimer. In contrast, spray of solutions containing 6 µM europium (Eu) with 3 µM S,S-mTDDGA diastereomer yielded only protonated and sodiated ligand ions; no Eu- S,S-mTDDGA ions were present. This mirrors the solution-phase behavior, where the S,S diestereomer extracts lanthanides with significantly less efficiency than the R,S diestereomer. Experiments with holmium and samarium yielded the same results. Increasing the lanthanide concentration to 30 µM resulted in limited formation of lanthanide-S,S-mTDDGA clusters. However, the abundance of the [Ln(S,S-mTDDGA)3]3+ clusters remained low compared to the abundance of the protonated and sodiated ligand. Additionally, the abundance of the [Ln(R,S-mTDDGA)3]3+ cluster ion was significantly higher than that of the protonated and sodiated ligand in experiments with spray solutions containing 30 µM lanthanide with R,S-mTDDA, supporting the hypothesis that R,S-mTDDGA has a higher intrinsic affinity for the lanthanides than S,S-mTDDGA. Thus, the intrinsic complexation with lanthanides appears to offer a route to differentiate diglycolamide diestereomers in the gas-phase. Novelty (20): Stereoisomers can be differentiated in the gas phase by exploitation of differences in metal complexation chemistry.

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

Achieving Cyber-Resilience for Power Systems using a Learning, Model-Assisted Blockchain Framework

The secure integration and management of distributed energy resources (DER) and power aggregators in the electric grid requires secure communications and a physics-aware Command and Control (C2) strategy. A Blockchain (BC)-based overlay network was developed to provide a security layer for the existing power grid network that mitigates risks in current and legacy network and C2 protocols. By integrating a Model-Assisted Machine Learning (MAML) framework with a Secure Blockchain Overlay Network (SBON) a defense-in-depth strategy was achieved. In our approach, the MAML framework leveraged a smart contract framework to gather network data and learn the dynamics of DER to develop detection strategies for attacks targeting sensors and actuators used by DER. The MAML framework learned dynamical systems models for individual DERs to detect sensor attacks. For DER we utilized a Digital Twin (DT) to accelerate the learning process for a model resistant to stealthy attacks. The project created DT for PV inverters and BESS. The DTs were coupled with a model-assisted, data-driven learning of DER behavior. Specifically, we evaluated architectures for model-based learning with model-free fine-tuning. Additionally, differential privacy techniques were used to obfuscate data, while still allowing the computation of attack detection results based on obfuscated data. The SBON developed leverages a private permissioned blockchain network orchestrated with the Hyperledger Fabric framework. To connect the cyber world, which orchestrates the blockchain fabric, and the physical world where the power network resides, we developed a system implementation to enable the secure interaction of the physical world and the abstracted blockchain.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

A differential game solution to the Coplanar tail-chase aerial combat problem

Numerical results obtained in a simplified version of the one on one aerial combat problem are presented. The primary aim of the data is to specify the roles of pursuer and evader as functions of the relative geometry and of the significant physical parameters of the problem. Numerical results are given in a case in which the slower aircraft is more maneuverable than the faster aircraft. A third order dynamic model of the relative motion is described, for which the state variables are relative range, bearing, and heading. The ranges at termination are arbitary in the present version of the problem, so the weapon systems of both aircraft can be visualized as forward firing high velocity weapons, which must be aimed at the tail pipe of the evader. It was found that, for the great majority of the ralative geometries, each aircraft can evade the weapon system of the other.

Merz, A. W.↗
NASA NTRS2017

Implicit Formulation of Muscle Dynamics in OpenSim

Astronauts lose bone and muscle mass during spaceflight. Exercise countermeasure is the primary method for counteracting bone and muscle mass loss in space. New spacecraft exercise device concepts are currently being developed for the NASAs new crew exploration vehicle. The NASA Digital Astronaut Project (DAP) uses computational modeling to help determine if the new exercise devices will be effective as countermeasures. The NASA Digital Astronaut Project is developing the ability to utilize predictive simulation to provide insight into the change in kinematics and kinetics with a change in device and gravitational environment (1-g versus 0-g). For example, in space exercise the subject's body weight is applied in addition to the loads prescribed for musculoskeletal maintenance. How and where these loads are applied obviously directly impacts bone and tissue loads. Additionally, due to space vehicle structural requirements, exercise devices are often placed on vibration isolation systems. This changes the apparent impedance or stiffness of the device as seen by the user. Data collection under these conditions is often impractical and limited. Predictive modeling provides a means to have a virtual subject to test hypotheses. Predictive simulation provides a virtual subject for which we are able to perform studies such as sensitivity to device loading and vibration isolation without the need for laboratory kinematic or kinetic test data.Direct Collocation optimization provides an efficient means to perform task based optimization and predictive modeling. It is relatively straight forward to structure a physical exercise task in a Direct Collocation mathematical formulation: perform a motion such that you start at an initial pose, achieve a given amount of deflection i.e a squat, return to the initial pose, and minimize muscle activation cost. Direct Collocation is advantageous in that it does not require numerical integration to evaluate the objective function. Instead, the system dynamics are transformed to discrete time and the optimizer is constrained such that the solution is not considered to be a valid unless the dynamic equations are satisfied at all time points. The simulation and optimization are effectively done simultaneously. Due to the implicit integration, time steps can be more coarse than in a differential equation solver. In a gait scenario this means that that the model constraints and cost function are evaluated at 100 nodes in the gait cycle versus 10,000 integration steps in a variable-step forward dynamic simulation. Furthermore, no time is wasted on accurate simulations of movements that are far from the optimum. Constrained optimization algorithms require a Jacobian matrix that contains the partial derivatives of each of the dynamic constraints with respect to of each of the state and control variables at all time points. This is a large but sparse matrix. An implicit dynamics formulation requires computation of the dynamic residuals f as a function of the states x and their derivatives, and controls u:f(x, dxdt, u) 0If the dynamics of musculoskeletal system are formulated implicitly, the Jacobian elements are often available analytically, eliminating the need for numerical differentiation; this is obviously computationally advantageous. Additionally, implicit formulation of musculoskeletal dynamics do not suffer from singularities from low mass bodies, zero muscle activation, or other stiff system or

physical exercise↗
DOE OSTI2025

Machine Learning Techniques for Data Reduction of Climate Applications

Scientists conduct large-scale simulations to compute derived quantities-of-interest (QoI) from primary data. Often, QoI are linked to specific features, regions, or time intervals, such that data can be adaptively reduced without compromising the integrity of QoI. For many spatiotemporal applications, these QoI are binary in nature and represent presence or absence of a physical phenomenon. We present a pipelined compression approach that first uses neural-network-based techniques to derive regions where QoI are highly likely to be present. Then, we employ a Guaranteed Autoencoder (GAE) to compress data with differential error bounds. GAE uses QoI information to apply low-error compression to only these regions. This results in overall high compression ratios while still achieving downstream goals of simulation or data collections. Experimental results are presented for climate data generated from the E3SM Simulation model for downstream quantities such as tropical cyclone and atmospheric river detection and tracking. These results show that our approach is superior to comparable methods in the literature.

Li, Xiao [University of Florida]↗
DOE OSTI2021

Poisson_CCD: A dedicated simulator for modeling CCDs

A dedicated simulator, Poisson_CCD, has been constructed which models astronomical CCDs by solving Poisson’s equation numerically and simulating charge transport within the CCD. The potentials and free carrier densities within the CCD are self-consistently solved for, giving realistic results for the charge distribution within the CCD storage wells. The simulator has been used to model the CCDs that are being used to construct the LSST digital camera. The simulator output has been validated by comparing its predictions with several different types of CCD measurements, including astrometric shifts, brighter-fatter induced pixel–pixel covariances, saturation effects, and diffusion spreading. The code is open source and freely available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

On the force exerted on a non-spherical asymmetric dust grain from homogeneous, stationary, isotropic, non-magnetized plasma

Whereas the conventional wisdom suggests that the force between non-magnetized homogeneous, stationary, isotropic plasma, and the dust grain is only possible for the case of relative plasma–grain velocity, it is shown that stationary non-spherical asymmetric dust grain immersed in stationary, non-magnetized, isotropic plasma can experience a force caused by the grain–plasma interactions. The component of the force due to scattering of plasma particles in the limit of infinite Debye length is considered analytically. Both the particle scattering and absorption force components are modeled numerically in the limits of infinite and finite Debye length using a newly developed 2D3V Aspherical Particle-in-Cell code. The code simulates interactions of dust grain of selected non-spherical asymmetric shape with plasmas using dust shape conforming coordinates. Here, the simulations confirm the existence of the force on non-spherical asymmetric grain in stationary non-magnetized plasma and show that the plasma screening effects can lead to reversal of the force direction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Field-level reconstruction from foreground-contaminated 21-cm maps

Current and upcoming 21-cm experiments will soon be able to map 21-cm spatial fluctuations in three dimensions for a wide range of redshifts. However, bright foreground contamination and the nature of radio interferometry create significant challenges, making it difficult to access rich cosmological information from the Fourier modes that lie within the “foreground wedge”. Here, in this work, we introduce two approaches aiming to reconstruct the full 21-cm density field, including the missing modes in the wedge: (a) a field-level inference under an effective field theory (EFT) framework; (b) a diffusion-based deep generative model trained on simulations. Under the EFT framework, we implement a fully differentiable forward model that maps the initial conditions of matter fluctuations to the observed, foreground-filtered 21-cm maps. This enables a gradient-based sampler to simultaneously sample the initial conditions and bias parameters, allowing a physically motivated mode reconstruction. Alternatively, we apply a variational diffusion model to perform 21-cm density reconstruction at the map level. Our model is trained on semi-numerical simulations over a wide range of astrophysical parameters. Our results from both approaches should provide improved cosmological constraints from the field level and also enable cross-correlation between experiments that have little or no overlapping modes.

cosmological perturbation theory↗
DOE OSTI2024

Fiducial and differential cross-section measurements of electroweak $W\gamma jj$ production in $pp$ collisions at $\sqrt{s} = 13$ TeV with the ATLAS detector

The observation of the electroweak production of a W boson and a photon in association with two jets, using pp collision data at the Large Hadron Collider at a centre of mass energy of $\sqrt{s} = 13$ TeV, is reported. The data were recorded by the ATLAS experiment from 2015 to 2018 and correspond to an integrated luminosity of 140 fb -1 . This process is sensitive to the quartic gauge boson couplings via the vector boson scattering mechanism and provides a stringent test of the electroweak sector of the Standard Model. Events are selected if they contain one electron or muon, missing transverse momentum, at least one photon, and two jets. Multivariate techniques are used to distinguish the electroweak $W\gamma jj$ process from irreducible background processes. The observed significance of the electroweak $W\gamma jj$ process is well above six standard deviations, compared to an expected significance of 6.3 standard deviations. Fiducial and differential cross sections are measured in a fiducial phase space close to the detector acceptance, which are in reasonable agreement with leading order Standard Model predictions from MADGRAPH5+PYTHIA8 and SHERPA. The results are used to constrain new physics effects in the context of an effective field theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Update on NNL TSL Evaluations and Validation [Slides]

12 new/revised TSL evaluations were contributed to ENDF/B-VIII.1. The thermal diffusion length $L$ is an integral property of a single material’s absorption and scattering cross sections (both differential and integral). No other neutron reactions or materials are involved. The MC21-calculated $L$ for water is consistent with the spread of experimental data and is sufficiently sensitive to different H-H 2 O TSL physics models to use the method as a TSL integral performance benchmark. Modern high-quality diffusion experiments at elevated $T$ would allow direct low-cost physics benchmarking of water TSLs when public elevated-$T$ critical benchmarks are limited. NNL is working with LLNL and RPI to develop experimental PNDA capability. An ICSBEP Fundamental Physics benchmark of the Nassar and Murphy PNDA experiment is being developed to provide an example of how to evaluate these type of benchmarks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗