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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

SINDA '85/FLUINT, Version 2.2

SINDA, Systems Improved Numerical Differencing Analyzer, software system developed for solving physical problems governed by diffusion-type equations modeled by lumped-parameter representations. Used as general thermal analyzer with resistor-and-capacitor-network representations, adapted to wide range of problems represented by such differential equations as Fourier, Poisson, or LaPlace. Software system updated. FLUINT, FLUid INTegrator, is advanced, one-dimensional fluid-analysis program solving equations of arbitrary fluid-flow networks. Written in FORTRAN and ASSEMBLER.

Cullimore, Brent↗
NASA NTRS1995

Numerical investigations of low-density nozzle flow by solving the Boltzmann equation

A two-dimensional finite-difference code to solve the BGK-Boltzmann equation has been developed. The solution procedure consists of three steps: (1) transforming the BGK-Boltzmann equation into two simultaneous partial differential equations by taking moments of the distribution function with respect to the molecular velocity u(sub z), with weighting factors 1 and u(sub z)(sup 2); (2) solving the transformed equations in the physical space based on the time-marching technique and the four-stage Runge-Kutta time integration, for a given discrete-ordinate. The Roe's second-order upwind difference scheme is used to discretize the convective terms and the collision terms are treated as source terms; and (3) using the newly calculated distribution functions at each point in the physical space to calculate the macroscopic flow parameters by the modified Gaussian quadrature formula. Repeating steps 2 and 3, the time-marching procedure stops when the convergent criteria is reached. A low-density nozzle flow field has been calculated by this newly developed code. The BGK Boltzmann solution and experimental data show excellent agreement. It demonstrated that numerical solutions of the BGK-Boltzmann equation are ready to be experimentally validated.

Deng, Zheng-Tao↗
NASA NTRS2004

A Model for the Oxidation of Carbon Silicon Carbide Composite Structures

A mathematical theory and an accompanying numerical scheme have been developed for predicting the oxidation behavior of carbon silicon carbide (C/SiC) composite structures. The theory is derived from the mechanics of the flow of ideal gases through a porous solid. The result of the theoretical formulation is a set of two coupled nonlinear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations are solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of the spatial location and time. The local rate of carbon oxidation is determined using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The nonlinear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual method, allowing for the solution of the differential equations numerically. The numerical method is demonstrated by utilizing the method to model the carbon oxidation and weight loss behavior of C/SiC specimens during thermogravimetric experiments. The numerical method is used to study the physics of carbon oxidation in carbon silicon carbide composites.

Sullivan, Roy M.↗
NASA NTRS2018

Picometer Differential Wavefront Metrology by Nonlinear Zernike Wavefront Sensing for LUVOIR

We propose the Nonlinear Zernike wavefront sensor (NLZWFS) for out-of-band differential wavefront sensing to augment primary mirror stability on LUVOIR and similar mission concepts during exoplanet coronagraphy. This new data analysis paradigm involving a full polychromatic scalar physical optics model for the phase-shifting Zernike wavefront sensor removes the linearity assumptions which would otherwise prevent accurate sensing. We show Monte-Carlo simulations of NLZWFS and focus-diverse phase retrieval to understand the exposure times necessary to achieve picometer-level stability in the telescope wavefront.

Redding, David C.↗
DOE OSTI2022

Publishing unbinned differential cross section results

Machine learning tools have empowered a qualitatively new way to perform differential cross section measurements whereby the data are unbinned, possibly in many dimensions. Unbinned measurements can enable, improve, or at least simplify comparisons between experiments and with theoretical predictions. Furthermore, many-dimensional measurements can be used to define observables after the measurement instead of before. There is currently no community standard for publishing unbinned data. While there are also essentially no measurements of this type public, unbinned measurements are expected in the near future given recent methodological advances. The purpose of this paper is to propose a scheme for presenting and using unbinned results, which can hopefully form the basis for a community standard to allow for integration into analysis workflows. This is foreseen to be the start of an evolving community dialogue, in order to accommodate future developments in this field that is rapidly evolving.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

On single-crystal total scattering data reduction and correction protocols for analysis in direct space

Data reduction and correction steps and processed data reproducibility in the emerging single-crystal total-scattering-based technique of three-dimensional differential atomic pair distribution function (3D-ΔPDF) analysis are explored. All steps from sample measurement to data processing are outlined using a crystal of CuIr 2 S 4 as an example, studied in a setup equipped with a high-energy X-ray beam and a flat-panel area detector. Computational overhead as pertains to data sampling and the associated data-processing steps is also discussed. Various aspects of the final 3D-ΔPDF reproducibility are explicitly tested by varying the data-processing order and included steps, and by carrying out a crystal-to-crystal data comparison. Situations in which the 3D-ΔPDF is robust are identified, and caution against a few particular cases which can lead to inconsistent 3D-ΔPDFs is noted. Although not all the approaches applied herein will be valid across all systems, and a more in-depth analysis of some of the effects of the data-processing steps may still needed, the methods collected herein represent the start of a more systematic discussion about data processing and corrections in this field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Observable CMB B -modes from cosmological phase transitions

A B -mode polarization signal in the cosmic microwave background (CMB) is widely regarded as smoking gun evidence for gravitational waves produced during inflation. Here, we demonstrate that tensor perturbations sourced during noninflationary epochs can yield non-negligible B -mode signals, which can in principle complicate the interpretation of future observational data. As a case study, we consider tensor perturbations sourced in the bubble collision stage of a first-order cosmological phase transition occurring in a secluded dark sector. Although phase transitions arise from causal subhorizon physics, they nevertheless exhibit a white noise power spectrum on superhorizon scales. Power is suppressed on the large scales relevant for CMB B -mode polarization, but it is not necessarily negligible. We show that for appropriately chosen phase transition parameters, the maximal B -mode amplitude can compete with inflationary predictions that can be tested with current and future experiments. These scenarios can be differentiated by performing measurements on multiple angular scales, since the phase transition signal predicts peak power on smaller scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A measurement of proton-carbon forward scattering in a proof-of-principle test of the EMPHATIC spectrometer

The next generation of long-baseline neutrino experiments will be capable of precision measurements of neutrino oscillation parameters, precision neutrino-nucleus scattering, and unprecedented sensitivity to physics beyond the Standard Model. Reduced uncertainties in neutrino fluxes are necessary to achieve high precision and sensitivity in these future precise neutrino measurements. New measurements of hadron-nucleus interaction cross sections are needed to reduce uncertainties of neutrino fluxes. We report measurements of the differential cross-section as a function of scattering angle for proton-carbon interactions with a single charged particle in the final state at beam momenta of 20, 30, and 120 GeV/c. These measurements are the result of a beam test for EMPHATIC, a hadron-scattering and hadron-production experiment. The total, elastic and inelastic cross-sections are also extracted from the data and compared to previous measurements. These results can be used in current and future long-baseline neutrino experiments, and demonstrate the feasibility of future measurements by an upgraded EMPHATIC spectrometer.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Atomically sharp domain walls in an antiferromagnet

The interest in understanding scaling limits of magnetic textures such as domain walls spans the entire field of magnetism from its physical fundamentals to applications in information technologies. Here, we explore antiferromagnetic CuMnAs in which imaging by x-ray photoemission reveals the presence of magnetic textures down to nanoscale, reaching the detection limit of this established microscopy in antiferromagnets. We achieve atomic resolution by using differential phase-contrast imaging within aberration-corrected scanning transmission electron microscopy. We identify abrupt domain walls in the antiferromagnetic film corresponding to the Néel order reversal between two neighboring atomic planes. Our work stimulates research of magnetic textures at the ultimate atomic scale and sheds light on electrical and ultrafast optical antiferromagnetic devices with magnetic field–insensitive neuromorphic functionalities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Molecular remodeling in Populus PdKOR RNAi roots profiled using LC-MS/MS proteomics

Plant endo-β-1,4-glucanases belonging to the Glycoside Hydrolase Family 9 have functional roles in cell wall biosynthesis and remodeling via endohydrolysis of (1→4)-β-D-glucosidic linkages. Modification of cell wall chemistry via RNAi-mediated downregulation of Populus deltoides KOR1 (PdKOR), a endo-β-1,4-glucanase gene, in Populus deltoides has been shown to have functional consequences for the composition of secondary metabolome and the ability of modified roots to interact with beneficial microbes. The molecular remodeling that underlies the observed differences at metabolic, physiological, and morphological levels in roots is not well understood. Here we used a LC-MS/MS-based proteome profiling approach to survey the molecular remodeling in root tissues of PdKOR and control plants. A total of 14316 peptides were identified and these mapped to 7139 P. deltoides proteins. Based on 90% sequence identity, the measured protein accessions represent 1187 functional protein groups. Analysis of GO categories and specific individual proteins showed differential expression of proteins relevant to plant-microbe interactions, cell wall chemistry, and metabolism. The new proteome dataset serves as a useful resource for deriving new hypotheses and empirical testing pertaining to functional roles of proteins and pathways in differential priming of plant roots to interactions with microbes.

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

Isochronous and period-doubling diagrams for symplectic maps of the plane

Symplectic mappings of the plane serve as key models for exploring the fundamental nature of complex behavior in nonlinear systems. Central to this exploration is the effective visualization of stability regimes, which enables the interpretation of how systems evolve under varying conditions. While the area-preserving quadratic Hénon map has received significant theoretical attention, a comprehensive description of its mixed parameter-space dynamics remain lacking. This limitation arises from early attempts to reduce the full two-dimensional phase space to a one-dimensional projection, a simplification that resulted in the loss of important dynamical features. Consequently, there is a clear need for a more thorough understanding of the underlying qualitative aspects. This paper aims to address this gap by revisiting the foundational concepts of reversibility and associated symmetries, first explored in the early works of G.D. Birkhoff. We extend the original framework proposed by Hénon by adding a period-doubling diagram to his isochronous diagram, which allows to represents the system’s bifurcations and the groups of symmetric periodic orbits that emerge in typical bifurcations of the fixed point. A qualitative and quantitative explanation of the main features of the region of parameters with bounded motion is provided, along with the application of this technique to other symplectic mappings, including cases of multiple reversibility. Modern chaos indicators, such as the Reversibility Error Method (REM) and the Generalized Alignment Index (GALI), are employed to distinguish between various dynamical regimes in the mixed space of variables and parameters. These tools prove effective in differentiating regular and chaotic dynamics, as well as in identifying twistless orbits and their associated bifurcations. Additionally, we discuss the application of these methods to real-world problems, such as visualizing dynamic aperture in accelerator physics, where our findings have direct relevance.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Spatial coupling of gyrokinetic simulations, a generalized scheme based on first-principles

In this paper, we present a scheme that spatially couples two gyrokinetic codes using first-principles. Coupled equations are presented and a necessary and sufficient condition for ensuring accuracy is derived. This new scheme couples both the field and the particle distribution function. The coupling of the distribution function is only performed once every few time-steps, using a five-dimensional (5D) grid to communicate the distribution function between the two codes. This 5D grid interface enables the coupling of different types of codes and models, such as particle and continuum codes, or delta-f and total-f models. Transferring information from the 5D grid to the marker particle weights is achieved using a new resampling technique. Demonstration of the coupling scheme is shown using two XGC gyrokinetic simulations for both the core and edge. We also apply the coupling scheme to two continuum simulations for a one-dimensional advection–diffusion problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Periodic Korteweg–de Vries soliton potentials generate quasisymmetric magnetic field strength in a finite plasma- β equilibrium

Quasisymmetry (QS) is a hidden symmetry of the magnetic field strength, B , that enables effective confinement of charged particles in a fully three-dimensional (3D) toroidal plasma equilibrium. Such equilibria are typically modeled by the ideal magnetohydrostatic (MHS) equations. The nonlinear, overdetermined nature of the QS MHS equations severely complicates our understanding of the interplay between 3D shaping, equilibrium properties such as pressure and rotational transform, and B . Progress has been made through expansions near the magnetic axis; however, a more comprehensive theory is desirable. Using a combination of analysis and regression on a large dataset of numerically optimized quasisymmetric stellarators, we demonstrate that there is a hidden lower dimensionality of B on a magnetic flux surface with connections to the theory of periodic solitons. We show that B on a flux surface is determined by three or at most four flux functions, each of which determines a critical value of the derivative of B along the field line. While being consistent with the near-axis models, our results are global and hold even on the last closed flux surface.

Differential geometry↗
DOE OSTI2023

Black and gray box learning of amplitude equations: Application to phase field systems

Here, we present a data-driven approach to learning surrogate models for amplitude equations and illustrate its application to interfacial dynamics of phase field systems. In particular, we demonstrate learning effective partial differential equations describing the evolution of phase field interfaces from full phase field data. We illustrate this on a model phase field system, where analytical approximate equations for the dynamics of the phase field interface (a higher-order eikonal equation and its approximation, the Kardar-Parisi-Zhang equation) are known. For this system, we discuss data-driven approaches for the identification of equations that accurately describe the front interface dynamics. When the analytical approximate models mentioned above become inaccurate, as we move beyond the region of validity of the underlying assumptions, the data-driven equations outperform them. In these regimes, going beyond black box identification, we explore different approaches to learning data-driven corrections to the analytically approximate models, leading to effective gray box partial differential equations.

42 ENGINEERING↗
DOE OSTI2022

Real-Time Interactive 4D-STEM Phase-Contrast Imaging From Electron Event Representation Data: Less computation with the right representation

The arrival of direct electron detectors (DED) with high frame-rates in the field of scanning transmission electron microscopy has enabled many experimental techniques that require collection of a full diffraction pattern at each scan position, a field which is subsumed under the name four dimensional-scanning transmission electron microscopy (4D-STEM). DED frame rates approaching 100 kHz require data transmission rates and data storage capabilities that exceed commonly available computing infrastructure. Current commercial DEDs allow the user to make compromises in pixel bit depth, detector binning or windowing to reduce the per-frame file size and allow higher frame rates. This change in detector specifications requires decisions to be made before data acquisition that may reduce or lose information that could have been advantageous during data analysis. The 4D Camera, a DED with 87 kHz frame-rate developed at Lawrence Berkeley National Laboratory, reduces the raw data to a linear-index encoded electron event representation (EER). Here we show with experimental data from the 4D Camera that linear-index encoded EER and its direct use in 4D-STEM phase contrast imaging methods enables real-time, interactive phase-contrast from large-area 4D-STEM datasets. Furthermore, we detail the computational complexity advantages of the EER and the necessary computational steps to achieve real-time interactive ptychography and center-of-mass differential phase contrast using commonly available hardware accelerators.

4D-STEM↗
DOE OSTI2024

Study of electro-physical properties of CsPbBr3 perovskite single crystals

Methods have been developed for the synthesis and growth of lead halide perovskites CsPbBr3, which can be used as detectors of optical, X-ray and gamma-radiation. Single crystals of these compounds were grown by the Bridgman method. The melt temperature was 580 °C, the gradient at the crystallization front was 5 K/cm, and the ampoule lowering speed was 3 mm/hour. In this work, the electro-physical properties of the inorganic perovskite CsPbBr3 were investigated. Crystals for research were cut from the beginning and end of the ingot. The optical absorption spectra were measured depending on the crystal thickness and temperature, the differential resistivity dif and the product  were determined from the I-V characteristics. It was established that dif practically does not change for samples made from the beginning and the end of the ingot, and is dif  21010 Ωcm at 100 V and 293 K.

Sklyarchuk, V.↗
DOE OSTI2022

Computational Fluid Dynamics and Heat Transfer Modeling of a Dimpled Heat Exchanger

Multiphysics and analytical calculations were conducted for a heat exchanger with passive, natural circulation flow. A glycol/water working fluid convects the heat to a dimpled heat exchanger shell, which subsequently transfers the heat to the soil, which acts as the ultimate heat sink. Because the system is fully-passive, it is not subject to the expenses, maintenance, and mechanical breakdowns associated with moving parts. Density, heat capacity, and thermal conductivity material properties were measured for various soil samples, and subsequently included as input for the soil heat conduction model. The soil model was coupled to a computational fluid dynamics (CFD) heat exchanger model that included the dynamic Smagorinsky large eddy simulation and k- omega turbulence models. The analysis showed that the fluid dynamics and heat transfer models worked properly, albeit at a slow pace. Nevertheless, the coupled CFD/heat conduction simulation ran long enough to determine a key parameter—the amount of heat conducted from the heat exchanger to the ground. This unique performance value, along with experimental data, was used as input for stand-alone, fast-running CFD models, as well as boundaries to obtain solutions to partial differential equations for soil heat conduction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

First experimental study of multiple orientation muon tomography, with image optimization in sparse data environments

Due to the high penetrating power of cosmic ray muons, they can be used to probe very thick and dense objects. As charged particles, they can be tracked by ionization detectors, determining the position and direction of the muons. With detectors on either side of an object, particle direction changes can be used to extract scattering information within an object. This can be used to produce a scattering intensity image within the object related to density and atomic number. Such imaging is typically performed with a single detector-object orientation, taking advantage of the more intense downward flux of muons, producing planar imaging with some depth-of-field information in the third dimension. Several simulation studies have been published with multi-orientation tomography, which can form a three-dimensional representation faster than a single orientation view. In this work we present the first experimental multiple orientation muon tomography study. Experimental muon-scatter based tomography was performed using a concrete filled steel drum with several different metal wedges inside, between detector planes. Data was collected from different detector-object orientations by rotating the steel drum. The data collected from each orientation were then combined using two different tomographic methods. Results showed that using a combination of multiple depth-of-field reconstructions, rather than a traditional inverse Radon transform approach used for CT, resulted in more useful images for sparser data. As cosmic ray muon flux imaging is rate limited, the imaging techniques were compared for sparse data. Using the combined depth-of-field reconstruction technique, fewer detector-object orientations were needed to reconstruct images that could be used to differentiate the metal wedge compositions.

Applied Physics (physics.app-ph)↗