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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Structure and Dynamics of Coronal Plasma

Brief summaries of the four published papers produced within the present performance period of NASA Grant NAGW-4081 are presented. The full text of the papers are appended to the report. The first paper titled "Coronal Structures Observed in X-rays and H-alpa Structures" was published in the Kofu Symposium proceedings. The study analyzes cool and hot behavior of two x-ray events, a small flare and a surge. It was found that a large H-alpha surge appears in x-rays as a very weak event, while a weak H-alpha feature corresponds to the brightest x-ray emission on the disk at the time of the observation. Calculations of the heating necessary to produce these signatures, and implications for the driving and heating mechanisms of flares vs. surges are presented. The second paper "Differential Magnetic Field Shear in an Active Region" has been published in The Astrophysical Journal. The study compared the three dimensional extrapolation of magnetic fields with the observed coronal structure in an active region. Based on the fit between observed coronal structure throughout the volume of the region and the calculated magnetic field configurations, the authors propose a differential magnetic field shear model for this active region. The decreasing field shear in the outer portions of the AR may indicate a continual relaxation of the magnetic field with time, corresponding to a net transport of helicity outward. The third paper "Difficulties in Observing Coronal Structure" has been published in the journal Solar Physics. This paper discusses the evidence that the temperature and density structure of the corona are far more complicated than had previously been thought. The discussion is based on five studies carried out by the group on coronal plasma properties, showing that any one x-ray instrument does see all of the plasma present in the corona, that hot and cool material may appear to be co-spatial at a given location in the corona, and that simple magnetic field extrapolations provide only a poor fit to the observed structure. The fourth paper "Analysis and Comparison of Loop Structures Imaged with NIXT and Yohkoh/SXT" has been published in Astronomy and Astrophysics. This paper analyzes and compares a variety of coronal loops, deriving loop pressure and emission measure from loop models. They are able to determine the volume filling factor in the corona, which is found to be in the range 0.001 to 0.01 for compact loops, and of order 1 for large structures. The small values suggest highly filamented structures, especially at lower temperatures.

Golub, Leon↗
NASA NTRS2005

Mathematical Metaphors: Problem Reformulation and Analysis Strategies

This paper addresses the critical need for the development of intelligent or assisting software tools for the scientist who is working in the initial problem formulation and mathematical model representation stage of research. In particular, examples of that representation in fluid dynamics and instability theory are discussed. The creation of a mathematical model that is ready for application of certain solution strategies requires extensive symbolic manipulation of the original mathematical model. These manipulations can be as simple as term reordering or as complicated as discovery of various symmetry groups embodied in the equations, whereby Backlund-type transformations create new determining equations and integrability conditions or create differential Grobner bases that are then solved in place of the original nonlinear PDEs. Several examples are presented of the kinds of problem formulations and transforms that can be frequently encountered in model representation for fluids problems. The capability of intelligently automating these types of transforms, available prior to actual mathematical solution, is advocated. Physical meaning and assumption-understanding can then be propagated through the mathematical transformations, allowing for explicit strategy development.

Thompson, David E.↗
NASA NTRS2017

Measuring N2 Pressure Using Cyanobacteria Discipline: Geomicrobiology

The evolution of Earth's atmosphere has been governed by biological evolution. Dinitrogen (N2) has been a major constituent of Earth's atmosphere throughout the planet's history, yet only a few constraints exist for the partial pressure of N2 (pN2). In this study we evaluate two new potential proxies for pN2: the physical spacing between heterocysts and the isotopic signature of nitrogen fixation in filamentous cyanobacteria. Heterocyst-forming filamentous cyanobacteria are some of the oldest photosynthetic microorganisms on Earth, and debated fossilized specimens have been found in sedimentary rocks as old as 2 Ga. These organisms overcome nitrogen limitation in their aqueous environment through cellular differentiation along their filaments. The specialized cells that develop, known as heterocysts, fix the nitrogen and laterally distribute it to neighboring cells along the filaments. Because the concentration of the dissolved N2 available to the filaments correlates directly with pN2, any preservable physiological response of the organism to the changed N2 availability constitutes a potential proxy for pN2. In the laboratory, we have examined how pN2 is reflected in the heterocyst spacing pattern and in the isotopic signature of nitrogen fixation by subjecting the representative species Anabaena cylindrica and Anabaena variabilis to different N2 partial pressures during growth at constant temperature and lighting (in media free of combined nitrogen). We show experimentally that the distance between heterocysts and the nitrogen isotope fractionation measured in bulk biomass reflect the pN2 experienced by Anabaena cylindrica. Current work is investigating these responses in Anabaena variabilis. When heterocystous cyanobacteria fossilize, these morphological and isotopic signatures should preserve information about pN2 at that time. Application of this relationship to the rock record may provide a paleoproxy to complement the two existing geobarometers.

nitrogen↗
DOE OSTI2025

Constrained or unconstrained? Neural-network-based equation discovery from data

Throughout many fields, practitioners often rely on differential equations to model systems. Yet, for many applications, the theoretical derivation of such equations and/or the accurate resolution of their solutions may be intractable. Instead, recently developed methods, including those based on parameter estimation, operator subset selection, and neural networks, allow for the data-driven discovery of both ordinary and partial differential equations (PDEs), on a spectrum of interpretability. The success of these strategies is often contingent upon the correct identification of representative equations from noisy observations of state variables and, as importantly and intertwined with that, the mathematical strategies utilized to enforce those equations. Specifically, the latter has been commonly addressed via unconstrained optimization strategies. Representing the PDE as a neural network, we propose to discover the PDE (or the associated operator) by solving a constrained optimization problem and using an intermediate state representation similar to a physics-informed neural network (PINN). The objective function of this constrained optimization problem promotes matching the data, while the constraints require that the discovered PDE is satisfied at a number of spatial collocation points. We present a penalty method and a widely used trust-region barrier method to solve this constrained optimization problem, and we compare these methods on numerical examples. Our results on several example problems demonstrate that the latter constrained method outperforms the penalty method, particularly for higher noise levels or fewer collocation points. This work motivates further exploration into using sophisticated constrained optimization methods in scientific machine learning, as opposed to their commonly used, penalty-method or unconstrained counterparts. For both of these methods, we solve these discovered neural network PDEs with classical methods, such as finite difference methods, as opposed to PINNs-type methods relying on automatic differentiation. Here, we briefly highlight how simultaneously fitting the data while discovering the PDE improves the robustness to noise and other small, yet crucial, implementation details.

Data-driven discovery↗
DOE OSTI2021

Designing optimal experiments: an application to proton Compton scattering

Interpreting measurements requires a physical theory, but the theory’s accuracy may vary across the experimental domain. To optimize experimental design, and so to ensure that the substantial resources necessary for modern experiments are focused on acquiring the most valuable data, both the theory uncertainty and the expected pattern of experimental errors must be considered. We develop a Bayesian approach to this problem, and apply it to the example of proton Compton scattering. Chiral Effective Field Theory (χEFT) predicts the functional form of the scattering amplitude for this reaction, so that the electromagnetic polarizabilities of the nucleon can be inferred from data. With increasing photon energy, both experimental rates and sensitivities to polarizabilities increase, but the accuracy of χEFT decreases. Our physics-based model of χEFT truncation errors is combined with present knowledge of the polarizabilities and reasonable assumptions about experimental capabilities at HI γ S and MAMI to assess the information gain from measuring specific observables at specific kinematics, i.e. to determine the relative amount by which new data are apt to shrink uncertainties. The strongest gains would likely come from new data on the spin observables Σ 2x and Σ 2x' at ω≃140 to 200 MeV and 40° to 120°. These would tightly constrain γ E1E1 –γ E1M2 . New data on the differential cross section between 100 and 200 MeV and over a wide angle range will substantially improve constraints on α E1 –β M1 , γ π and γ M1M1 –γ M1E2 . Good signals also exist around 160 MeV for Σ 3 and Σ 2z' . As a result, such data will be pivotal in the continuing quest to pin down the scalar polarizabilities and refine understanding of the spin polarizabilities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Complete quantification of parametric uncertainties in ( d , p ) transfer reactions

Deuteron-induced transfer reactions are a popular probe in nuclear structure and nuclear astro physics studies. The interpretation of these transfer measurements relies on reaction theory that takes as input effective interactions between the nucleons and the target nucleus. Previous work quantified the uncertainty associated with the optical potentials between the nucleons and the target. In this study, we extend that work by also including the parameters of the mean field associated with the overlap function of the final bound state, thus obtaining the full parametric uncertainty on transfer observables. We use Bayesian Markov Chain Monte Carlo simulations to obtain parameter posterior distributions. We use elastic-scattering cross sections to constrain the optical potential parameters and use the asymptotic normalization coefficient of the final state to constrain the bound state interaction. We then propagate these posteriors to the transfer angular distributions and obtain confidence intervals for this observable. We study (d,p) reactions on 14 C, 16 O, and 48 Ca at energies in the range $E$ $d$ = 7 - 24 MeV. Our results show a strong reduction in uncertainty by using the asymptotic normalization coefficient as a constraint, particularly for those reactions most sensitive to ambiguities in the mean field. For those reactions, the importance of constraining the bound state interaction is equal to that of constraining the optical potentials. Furthermore, the case of 14 C is an outlier because the cross section is less sensitive to the nuclear interior. When minimal constraints are used on the parameters of the nucleon-target interaction, the 1σ uncertainties on the differential cross sections are large (~ 500-100%). However, if elastic-scattering data and the asymptotic normalization coefficient are used in the analysis, with an error of 10% (5%), this uncertainty reduces to ~ 30% (~ 15%).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Resolving the Discontinuity Between Thermal Scattering Data and Fast Data

The task of resolving discontinuities between thermal cross sections and fast cross sections in AMPX processing is part of a broader project to provide guidance on the potential effect of unknown thermal neutron scattering law (TNSL) data, primarily from the standpoint of criticality safety. The primary goal of this work is to deliver code fixes in AMPX that increase confidence that SCALE libraries provide physically correct cross sections across the full range of neutron energies encountered in transport applications. This report discusses the resolution of two sources of discontinuity. The first issue occurs when the Bragg edges in an evaluated nuclear data file do not extend to 5 eV. The second issue concerns parameters chosen during the processing of the nuclear data library—specifically, the parameter that defines the energy at which the thermal data ends and the fast nuclide data begins. This report also discusses a third issue, which is related to the calculation of probability distributions from a densely gridded double-differential cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

In-situ sensor monitoring of multi-class gas porosity formation in laser powder bed fusion using convolutional neural network

In-situ monitoring of defect formation remains a significant challenge in the laser powder bed fusion (LPBF) process. Recent advances have enabled real-time defect detection with machine learning and in-situ sensing technologies; however, most studies focus on binary classification of keyhole pores, limiting nuanced multi-class pore differentiation and formation mechanisms. This work introduces a multi-class pore detection framework (no pore, small pores < 15 µm, and large pores > 15 µm) by leveraging photodiode sensor data alongside high-fidelity synchrotron X-ray imaging. The 15 µm threshold is selected to distinguish between two fundamentally different defect mechanisms, following the physical size-mechanism boundary established by prior high-resolution synchrotron X-ray characterization of Al6061 LPBF. Distinguishing these classes is critical because large keyhole pores are structurally detrimental, whereas small gas pores are often benign, requiring different process control strategies. Thermal emission monitoring data collected simultaneously with high-speed X-ray imaging at the Stanford Synchrotron Radiation Lightsource (SSRL), are correlated with subsurface melt pool dynamics to establish ground truth. Continuous Wavelet Transform (CWT) with optimized parameters converts the photodiode time-series signals into time–frequency images, facilitating feature extraction. Convolutional Neural Networks (CNN) are then applied for real-time multi-class pore classification in an average inference time of 1 ms per signal window. It achieves 79% accuracy and an Area Under the Receiver Operating Characteristic curve (AUC ROC) score of 0.89 with five-fold cross-validation. The results demonstrate that coupling CWT-based feature engineering with CNN architecture enables reliable multi-class pore detection in Al6061 builds using affordable in-situ sensors. This approach advances scalable and affordable quality assurance in additive manufacturing by moving beyond binary defect detection toward more nuanced classification of porosity mechanisms with in-situ sensors and machine learning.

Laser powder bed fusion, Multi-class pores, In-sit↗
DOE OSTI2020

Quantum nondemolition photon counting with a hybrid electromechanical probe

Quantum nondemolition (QND) measurements of photons is a much pursued endeavor in the field of quantum optics and quantum information processing. Here, we propose a hybrid optoelectromechanical platform that integrates a cavity system with a hybrid electromechanical probe for QND photon counting. Building upon a mechanical-mode-mediated nonperturbative electro-optical dispersive coupling, our protocol performs the QND photon counting measurement by means of the current-voltage characteristics of the probe. In particular, we show that the peak voltage shift of the differential conductance is linearly dependent on the photon occupation number, thus providing a sensitive measure of the photon number, especially in the strong optomechanical coupling regime. Here, given that our proposed hybrid system is compatible with state-of-the-art experimental techniques, we discuss its implementations and anticipate applications in quantum optics and polariton physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Photoproduction of η mesons off the proton for 1.2 < E γ < 4.7GeV using CLAS at Jefferson Laboratory

Photoproduction cross sections are reported for the reaction $\gamma p\to p\eta$ using energy-tagged photons and the CLAS spectrometer at Jefferson Laboratory. The $\eta$ mesons are detected in their dominant charged decay mode, $\eta\to \pi^+\pi^-\pi^0$, and results on differential cross sections are presented for incident photon energies between 1.2 and 4.7 GeV. These new $\eta$ photoproduction data are consistent with earlier CLAS results but extend the energy range beyond the nucleon resonance region into the Regge regime. The normalized angular distributions are also compared with the experimental results from several other experiments, and with predictions of $\eta$ MAID\,2018 and the latest solution of the Bonn-Gatchina coupled-channel analysis. Differential cross sections $d\sigma/dt$ are presented for incident photon energies $E_\gamma > 2.9$ GeV ($W > 2.5$ GeV), and compared with predictions which are based on Regge trajectories exchange in the $t$-channel (Regge models). Overall, the data confirm the expected dominance of $\rho$, $\omega$ vector-meson exchange in an analysis by the Joint Physics Analysis Center.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Realizing Kagome Band Structure in Two-Dimensional Kagome Surface States of RV 6 Sn 6 (R = Gd, Ho)

Here, we report angle resolved photoemission experiments on a newly discovered family of kagome metals RV 6 Sn 6 (R = Gd, Ho). Intrinsic bulk states and surface states of the vanadium kagome layer are differentiated from those of other atomic sublattices by the real-space resolution of the measurements with a small beam spot. Characteristic Dirac cone, saddle point, and flat bands of the kagome lattice are observed. Our results establish the two-dimensional (2D) kagome surface states as a new platform to investigate the intrinsic kagome physics.

2-dimensional systems↗
NASA NTRS1975

In situ spectroradiometric calibration of EREP imagery and estuarine and coastal oceanography of Block Island sound and adjacent New York coastal waters

The author has identified the following significant results. The first part of the study resulted in photographic procedures for making multispectral positive images which greatly enhance the color differences in land detail using an additive color viewer. An additive color analysis of the geologic features near Willcox, Arizona using enhanced black and white multispectral positives allowed compilation of a significant number of unmapped geologic units which do not appear on geologic maps of the area. The second part demonstrated the feasibility of utilizing Skylab remote sensor data to monitor and manage the coastal environment by relating physical, chemical, and biological ship sampled data to S190A, S190B, and S192 image characteristics. Photographic reprocessing techniques were developed which greatly enhanced subtle low brightness water detail. Using these photographic contrast-stretch techniques, two water masses having an extinction coefficient difference of only 0.07 measured simultaneously with the acquisition of S190A data were readily differentiated.

Yost, E. F.↗
NASA NTRS1992

Modeling of the WSTF frictional heating apparatus in high pressure systems

In order to develop a computer program able to model the frictional heating of metals in high pressure oxygen or nitrogen a number of additions have been made to the frictional heating model originally developed for tests in low pressure helium. These additions include: (1) a physical property package for the gases to account for departures from the ideal gas state; (2) two methods for spatial discretization (finite differences with quadratic interpolation or orthogonal collocation on finite elements) which substantially reduce the computer time required to solve the transient heat balance; (3) more efficient programs for the integration of the ordinary differential equations resulting from the discretization of the partial differential equations; and (4) two methods for determining the best-fit parameters via minimization of the mean square error (either a direct search multivariable simplex method or a modified Levenburg-Marquardt algorithm). The resulting computer program has been shown to be accurate, efficient and robust for determining the heat flux or friction coefficient vs. time at the interface of the stationary and rotating samples.

Skowlund, Christopher T.↗
NASA NTRS1991

Numerical simulations of turbulent thermal convection with differential rotation

The solar nebula, from which the planets in our solar system formed, featured a disk of gas and dust grains in rapid, differential rotation, and at some stage was likely to have been unstable to thermal convection. This situation is suspected by many to lead to significant turbulent Reynolds stress production and angular momentum transport in such systems, and estimates of transport rates have been attempted from unsubstantiated phenomenological models. In order to determine the circumstances and physical conditions under which our own planetary system formed and to explain recent observations of young stellar systems, it is necessary to develop realistic models of heat and angular momentum transport for such flows. Developing an understanding of complicated flows featuring thermal convection, rotation, and shear is also of wide interest in stellar astrophysics and in planetary and terrestrial atmospheric studies. The ultimate objective is to develop workable models based on the numerical simulations for constructing global solar nebula models; viz., relatively simple prescriptions for heat and angular momentum fluxes from given system parameters (e.g., ratios of rotation, shear, and convective lapse rates) are characterized, quantified, and developed. Toward this end, our program has been to attempt to understand the behavior of the direct numerical simulations of Boussinesq convection, which, despite the complexity of the results, is still an overly simplified approximation to the real system and should be more amenable to analysis. These results are also intended to be tested against turbulence models, especially those designed for atmospheric boundary layers, and may provide a basis for subgrid-scale models. In order to make the numerical simulations more realistic with regard to the solar nebula problem, a fully compressible code that will allow incorporation of large density stratifications and realistic thermodynamic and radiative properties is developed. In order to explore the properties of these flows at the very high values of Re found in natural systems and the very low values of Pr found in most astrophysical contexts, we will need to employ large-eddy simulations for which we want to determine the most appropriate subgrid-scale model to incorporate.

Cabot, William H.↗
NASA NTRS1997

Structure and Dynamics of Coronal Plasmas

During the past year this grant has funded research in the interaction between magnetic fields and the hot plasma in the solar outer atmosphere. The following is a brief summary of the published papers, abstracts and talks which have been supported. The paper 'Coronal Structures Observed in X-rays and H-alpha Structures' was published in the Kofu Symposium proceedings. The study analyzes cool and hot behavior of two x-ray events, a small flare and a surge. We find that a large H-alpha surge appears in x-rays as a very weak event, while a weak H-alpha feature corresponds to the brightest x-ray emission on the disk at the time of the observation. Calculations of the heating necessary to produce these signatures, and implications for the driving and heating mechanisms of flares vs. surges are presented. A copy of the paper is appended to this report. The paper 'Differential Magnetic Field Shear in an Active Region' has been published in The Astrophysical Journal. We have compared the 3D extrapolation of magnetic fields with the observed coronal structure in an active region. Based on the fit between observed coronal structure throughout the volume of the region and the calculated magnetic field configurations, we propose a differential magnetic field shear model for this active region. The decreasing field shear in the outer portions of the AR may indicate a continual relaxation of the magnetic field with time, corresponding to a net transport of helicity outward. The paper 'Difficulties in Observing Coronal Structure' has been accepted for publication in the journal Solar Physics. In this paper we discuss the evidence that the temperature and density structure of the corona are far more complicated than had previously been thought. The discussion is based on five studies carried out by our group on coronal plasma properties, showing that any one x-ray instrument does see all of the plasma present in the corona, that hot and cool material may appear to be co-spatial at a given location in the corona, and that simple magnetic field extrapolations provide only a poor fit to the observed structure. A copy of the paper is appended to this report.

Golub, Leon↗
DOE OSTI2020

Comparison between PARFUME and Bison Using the AGR-2 Irradiation Experiment

This report documents comparisons between Fuel Model (PARFUME) model predictions versus Bison for selected compacts from the second irradiation test of the Advanced Gas Reactor (AGR) program that occurred from June 2010 to October 2013 in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). PARFUME is a fuel performance analysis and modeling code, used for evaluating gas-reactor tristructural isotropic (TRISO) coated particle fuel for prismatic, pebble bed, plate, and cylindrical type fuel geometries. PARFUME is an integrated mechanistic computer code that evaluates the thermal, mechanical, and physico-chemical behavior of TRISO coated-fuel particles and the probability for fuel failure given the particle-to-particle statistical variations in physical dimensions and material properties that arise during the fuel fabrication process. Bison is a nuclear fuel performance application built using the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element library (Permann 2020). Bison is capable of modeling multiple fuel forms in a wide variety of dimensions and geometries. It solves coupled nonlinear partial differential equations, including heat conduction, mechanics, fission product species transport etc., in a fully implicit manner. Comparisons between PARFUME and Bison were performed using four compacts from the AGR-2 experiment. Selected outputs were chosen based on their impact on the probability of the SiC layer failing. In general there was good agreement between PARFUME and Bison with the exception of predicting the gap formed between the buffer and IPyC. Further comparisons between PARFUME and Bison are planned to further develop Bison’s capabilities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Salt Hydrate Eutectic Mixtures for Near-ambient Thermal Energy Storage Applications

In this work, we synthesized and characterized Na2SO4.10H2O – Na2CO3.10H2O eutectic salt hydrate phase change material (PCM) using Differential Scanning Calorimetry (DSC), Attenuated total reflectance infrared spectroscopy (ATR-IR), and X-ray Diffraction (XRD) techniques. The eutectic point of the PCM mixture was predicted based on binary eutectic effect and ion exchange theory. The XRD and ATR-IR results indicate that the eutectic PCM mixture was formed based on only physical mixing without chemical interaction. Moreover, DSC results showed that the eutectic concentration is 50wt% Na2SO4.10H2O – 50wt% Na2CO3.10H2O with melting temperature of 25.18°C and melting enthalpy of 187.46 J/g. Also, 78% of melting enthalpy was recoverable after 10 melt/freeze cycles. This eutectic mixture has the potential to be used in building thermal comfort applications.

Gluesenkamp, Kyle↗