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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Modeling bristle lift-off in idealized brush seal configurations

We attempt in this paper to develop a model for the flow through brush seals and determine their elastic behavior in order to predict the dependence of brush/journal clearance on geometry and operating conditions. Several idealizations regarding brush seal configuration, flow conditions, and elastic behavior are made in the analysis in order to determine closed form parametric dependence. This formulation assumes that there is no initial interference between the bristle tip and the rotor. Also, interbristle, bristle-backing plate, and bristle-rotor friction is neglected. The bristle bundle or the brush seal as it is alternately called is assumed homogeneous and isotropic on a macroscopic scale so that a physical property like permeability is uniform. The fluid is assumed to be homogeneous, incompressible, viscous, and flowing under steady conditions. A schematic of a brush seal is shown. If the nominal bristle-shaft interference is absent then under static conditions the bristles may deflect axially due to the imposed pressure differential. This axial deflection may create a clearance permitting leakage flow in excess of that which occurs through the porous matrix formed by the bristle bundles. Under dynamic conditions the Couette flow created by shaft motion could be strong enough to cause bristle deflection and once again a clearance may develop. The paper proposes a means to determine this clearance (or at least describe its parametric dependence on geometry and operating conditions) under static as well as dynamic conditions.

Modi, Vijay↗
NASA NTRS1998

A Novel Sensor for Attitude Determination Using Global Positioning System Signals

An entirely new sensor approach for attitude determination using Global Positioning System (GPS) signals is developed. The concept involves the use of multiple GPS antenna elements arrayed on a single sensor head to provide maximum GPS space vehicle availability. A number of sensor element configurations are discussed. In addition to the navigation function, the array is used to find which GPS space vehicles are within the field-of-view of each antenna element. Attitude determination is performed by considering the sightline vectors of the found GPS space vehicles together with the fixed boresight vectors of the individual antenna elements. This approach has clear advantages over the standard differential carrier-phase approach. First, errors induced by multipath effects can be significantly reduced or eliminated altogether. Also, integer ambiguity resolution is not required, nor do line biases need to be determined through costly and cumbersome self-surveys. Furthermore, the new sensor does not require individual antennas to be physically separated to form interferometric baselines to determine attitude. Finally, development potential of the new sensor is limited only by antenna and receiver technology development unlike the physical limitations of the current interferometric attitude determination scheme. Simulation results indicate that accuracies of about 1 degree (3 omega) are possible.

Crassidis, John L.↗
NASA NTRS2003

Development of Algorithms for Control of Humidity in Plant Growth Chambers

Algorithms were developed to control humidity in plant growth chambers used for research on bioregenerative life support at Kennedy Space Center. The algorithms used the computed water vapor pressure (based on measured air temperature and relative humidity) as the process variable, with time-proportioned outputs to operate the humidifier and de-humidifier. Algorithms were based upon proportional-integral-differential (PID) and Fuzzy Logic schemes and were implemented using I/O Control software (OPTO-22) to define and download the control logic to an autonomous programmable logic controller (PLC, ultimate ethernet brain and assorted input-output modules, OPTO-22), which performed the monitoring and control logic processing, as well the physical control of the devices that effected the targeted environment in the chamber. During limited testing, the PLC's successfully implemented the intended control schemes and attained a control resolution for humidity of less than 1%. The algorithms have potential to be used not only with autonomous PLC's but could also be implemented within network-based supervisory control programs. This report documents unique control features that were implemented within the OPTO-22 framework and makes recommendations regarding future uses of the hardware and software for biological research by NASA.

Costello, Thomas A.↗
NASA NTRS2010

Fabrication of Lanthanum Telluride 14-1-11 Zintl High-Temperature Thermoelectric Couple

The development of more efficient thermoelectric couple technology capable of operating with high-grade heat sources up to 1,275 K is key to improving the performance of radioisotope thermoelectric generators. Lanthanum telluride La3-xTe4 and 14-1-11 Zintls (Yb14MnSb11) have been identified as very promising materials. The fabrication of advanced high-temperature thermoelectric couples requires the joining of several dissimilar materials, typically including a number of diffusion bonding and brazing steps, to achieve a device capable of operating at elevated temperatures across a large temperature differential (up to 900 K). A thermoelectric couple typically comprises a heat collector/ exchanger, metallic interconnects on both hot and cold sides, n-type and ptype conductivity thermoelectric elements, and cold-side hardware to connect to the cold-side heat rejection and provide electrical connections. Differences in the physical, mechanical, and chemical properties of the materials that make up the thermoelectric couple, especially differences in the coefficients of thermal expansion (CTE), result in undesirable interfacial stresses that can lead to mechanical failure of the device. The problem is further complicated by the fact that the thermoelectric materials under consideration have large CTE values, are brittle, and cracks can propagate through them with minimal resistance. The inherent challenge of bonding brittle, high-thermal-expansion thermoelectric materials to a hot shoe material that is thick enough to carry the requisite electrical current was overcome. A critical advantage over prior art is that this device was constructed using all diffusion bonds and a minimum number of assembly steps. The fabrication process and the materials used are described in the following steps: (1) Applying a thin refractory metal foil to both sides of lanthanum telluride. To fabricate the n-type leg of the advanced thermoelectric couple, the pre-synthesized lanthanum telluride coupon was diffusion bonded to the metal foil using a thin adhesion layer. (2) Repeating a similar process for the 14-1-11 Zintl p-type leg of the advanced thermoelectric couple. (3) Bonding thick CTE-matched metal plates on the metallized lanthanum telluride and Yb14MnSb11 to form the hot and cold sides of the thermoelectric couple. The calculated conversion efficiency of such an advanced couple would be about 10.5 percent, about 35 percent better than heritage radioisotope thermoelectric technology that relies on Si-Ge alloys. In addition, unlike Si-Ge alloys, these materials can be combined with many other thermoelectric materials optimized for operation at lower temperatures to achieve conversion efficiency in excess of 15 percent (a factor of 2 increase over heritage technology).

Ravi, Vilupanur A.↗
NASA NTRS2009

Optimization of Layer Densities for Spacecraft Multilayered Insulation Systems

Numerous tests of various multilayer insulation systems have indicated that there are optimal densities for these systems. However, the only method of calculating this optimal density was by a complex physics based algorithm developed by McIntosh. In the 1970's much data were collected on the performance of these insulation systems with many different variables analyzed. All formulas generated included number of layers and layer density as geometric variables in solving for the heat flux, none of them was in a differentiable form for a single geometric variable. It was recently discovered that by converting the equations from heat flux to thermal conductivity using Fourier's Law, the equations became functions of layer density, temperatures, and material properties only. The thickness and number of layers of the blanket were merged into a layer density. These equations were then differentiated with respect to layer density. By setting the first derivative equal to zero, and solving for the layer density, the critical layer density was determined. Taking a second derivative showed that the critical layer density is a minimum in the function and thus the optimum density for minimal heat leak, this is confirmed by plotting the original function. This method was checked and validated using test data from the Multipurpose Hydrogen Testbed which was designed using McIntosh's algorithm.

Johnson, W. L.↗
DOE OSTI2025

Nyström type exponential integrators for strongly magnetized charged particle dynamics

Solving for charged particle motion in electromagnetic fields (i.e. the particle pushing problem) is a computationally intensive component of particle-in-cell (PIC) methods for plasma physics simulations. This task is especially challenging when the plasma is strongly magnetized due numerical stiffness arising from the wide range of time scales between highly oscillatory gyromotion and long term macroscopic behavior. A promising approach to solve these problems is by a class of methods known as exponential integrators that can solve linear problems exactly and are A-stable. This work extends the standard exponential integration framework to derive Nyström-type exponential integrators that integrates the Newtonian equations of motion as a second-order differential equation directly. In particular, we derive second-order and third-order Nyström-type exponential integrators for strongly magnetized particle pushing problems. Numerical experiments show that the Nyström-type exponential integrators exhibit significant improvement in computation speed over the standard exponential integrators.

general physics↗
DOE OSTI2022

First-passage time statistics on surfaces of general shape: Surface PDE solvers using Generalized Moving Least Squares (GMLS)

Here, we develop numerical methods for computing statistics of stochastic processes on surfaces of general shape with drift-diffusion dynamics d X t = a (X t ) dt + b(X t ) d W t . We formulate descriptions of Brownian motion and general drift-diffusion processes on surfaces. We consider statistics of the form u (x) = E x [$∫^{τ}_{0}$ g (X t ) dt ] + E x [ f (X τ )] for a domain Ω and the exit stopping time τ = inf t { t >0 | X i Ω}, where f , g are general smooth functions. For computing these statistics, we develop high-order Generalized Moving Least Squares (GMLS) solvers for associated surface PDE boundary-value problems based on Backward- Kolmogorov equations. We focus particularly on the mean First Passage Times (FPTs) given by the case f = 0, g = 1 where u (x) = E x [τ]. We perform studies for a variety of shapes showing our methods converge with high-order accuracy both in capturing the geometry and the surface PDE solutions. We then perform studies showing how statistics are influenced by the surface geometry, drift dynamics, and spatially dependent diffusivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

An implementation of a high-order generalized finite difference method for solving the time-harmonic cold plasma wave equation in toroidal geometry

A high-order physics-informed meshless finite difference numerical technique is introduced for solving the time-harmonic cold plasma wave equation in toroidal geometries, presenting a novel application of the generalized finite difference (GFD) method to plasma wave simulations. The algorithm employs an irregular distribution of computational points, with local point density informed by the shortest wavelength derived from the cold plasma dispersion relation. Numerical stability and robustness are addressed using regularization techniques. The algorithm, implemented for two spatial dimensions, solves for the wave electric field and is demonstrated to achieve convergence rates of $\mathcal{O}$($\mathcal{h}$ $\mathcal{P}$ )⁠. Verification tests reproduce plane wave solutions, and example simulations of ion cyclotron resonance heating and electron cyclotron resonance heating demonstrate its capability, approaching realistic tokamak plasma scenarios. This work contributes to laying a foundation for the GFD method to be used in more sophisticated, optimized, and physically realistic full-wave simulations in time-harmonic plasma wave research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

New insights into the electronic structure of α-U and δ-Pu

Here, this work presents the results of a theoretical study of the electronic structure of two actinide metals, α -U and δ -Pu. We compare our ab-initio results obtained with the recently developed self-consistent Vertex corrected GW approach with previously published experimental measurements such as photo-electron spectroscopy, for the occupied density of states, and bremsstralung isochromat spectroscopy (BIS) and inverse photo-electron spectroscopy (IPES), for the unoccupied density of states. Our ab-initio approach includes all important relativistic effects (it is based on Dirac’s equation) and it represents the first application of the Vertex corrected GW approach in the physics of actinides. Overall, our theoretical results are in good agreement with the experimental data, which supports the level of approximations which our theoretical method is based upon. By comparing our vertex corrected GW results with our results obtained with less sophisticated approaches (local density approximation and self-consistent GW) we differentiate the strength of correlation effects in Uranium and Plutonium. Also, our theoretical results allow us to elucidate the subtle differences between the previously published experimental BIS and IPES data on the unoccupied density of states in α -U.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Towards a Quantum Algorithm for the Incompressible Nonlinear Navier-Stokes Equations

In this work, we present novel concepts for quantum algorithms to solve transient, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the incompressible Navier-Stokes equations as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts solving nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. We propose a new framework based on matrix product states (MPSs) and matrix product operators (MPOs), in addition to the Krylov subspace methods. For example, the solution variables of the Navier-Stokes equations are represented by MPSs, and the linear and nonlinear terms are processed by MPOs. The time evolution of the operators is attained by a fast-forwarding algorithm using Krylov subspace methods. Furthermore, we discuss various techniques for efficient encoding of MPSs, measurement reduction for MPOs, and use of tensor operations to treat multi-variate, multi-physics characteristics of Navier-Stokes.

Gopalakrishnan Meena, Murali [ORNL] (ORCID:0000000↗
NASA NTRS1998

Mars Volatiles and Climate Surveyor (MVACS) Integrated Payload for the Mars Polar Lander Mission

The Mars Volatiles and Climate Surveyor (MVACS) integrated payload for the Mars Polar Lander will be launched in January 1999, with a scheduled landing on Mars' south-polar layered deposits in December 1999. Over the course of its 90-day nominal mission during the martian southern spring and summer seasons, it will make in situ measurements that will provide new insights into the behavior and distribution of martian volatiles, MVACS consists of four major instrument systems: a surface stereo imager (SSI), which will acquire multispectral stereo images of the surface and atmosphere; a 2-m robotic arm (RA), which will dig a O.5-m deep trench and acquire surface and subsurface samples that will be imaged by a focusable robotic arm camera (RAC), which will take close-up images of surface and subsurface samples at a spatial resolution of 21 micron; a meteorology package (MET), which will make the first measurements of surface pressure, temperature, and winds in Mars' southern hemisphere and employ a tunable diode laser (TDL) spectrometer to measure the water-vapor concentration and isotopic composition of CO2 in the martian atmosphere; and a thermal and evolved gas analyzer (TEGA), which will use differential scanning calorimetry and TDL-evolved gas analysis to determine the concentrations of ices, adsorbed volatiles, and volatile-bearing minerals in surface and subsurface soil samples. The unique in situ measurements made by MVACS at its high-latitude landing site will define a number of important aspects of the physical, isotopic, and chemical nature of the martian near-surface and subsurface environment that will be valuable in better understanding Mars meteorites and returned samples, as well as in the search for martian resources that could be utilized by humans.

Paige, D. A.↗
DOE OSTI2022

Trace Impurities Identified as Forensic Signatures in CMX-5 Fuel Pellets Using X-ray Spectroscopic Techniques

We report that small-particle analysis is a highly promising emerging forensic tool for analysis of interdicted special nuclear materials. Integration of microstructural, morphological, compositional, and molecular impurity signatures could provide significant advancements in forensic capabilities. We have applied rapid, high-sensitivity, hard X-ray synchrotron chemical imaging to analyze impurity signatures in two differently fabricated fuel pellets from the 5th Collaborative Materials Exercise (CMX5) of the IAEA Nuclear Forensics International Working Group. The spatial distributions, chemical compositions, and morphological and molecular characteristics of impurities were evaluated using X-ray absorption near-edge structure (XANES) and X-ray fluorescence chemical imaging to discover principal impurities, their granularity, particle sizes, modes of occurrence (distinct grains vs incorporation in the UO 2 lattice), and sources and mechanisms of incorporation. Differences in UO 2+x stoichiometry were detected at the microscale in nominally identical UO 2 ceramics (CMX5-A and CMX5-B), implying the presence of multiple UO 2 host phases with characteristic microstructures and feedstock compositions. Al, Fe, Ni, W, and Zr impurities and integrated impurity signature analysis identified distinctly different pellet synthesis and processing methods. For example, two different Al, W, and Zr populations in the CMX5-B sample indicated a more complex processing history than the CMX5-A sample. K-edge XANES measurements reveal both metallic and oxide forms of Fe and Ni but with different proportions between each sample. Altogether, these observations suggest multiple sources of impurities, including fabrication (e.g., force-sieving) and feedstock (mineral oxides). This study demonstrates the potential of synchrotron techniques to integrate different signatures across length scales (angstrom to micrometer) to detect and differentiate between contrasting UO 2 fuel fabrication techniques.

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

Peroxide-Selective Reduction of O 2 at Redox-Inactive Rare-Earth(III) Triflates Generates an Ambiphilic Peroxide

Metal peroxides are key species involved in a range of critical biological and synthetic processes. Rare-earth (group III and the lanthanides; Sc, Y, La–Lu) peroxides have been implicated as reactive intermediates in catalysis; however, reactivity studies of isolated, structurally characterized rare-earth peroxides have been limited. Herein, we report the peroxide-selective (93–99% O22-) reduction of dioxygen (O 2 ) at redox-inactive rare-earth triflates in methanol using a mild metallocene reductant, decamethylferrocene (Fc*). The first molecular praseodymium peroxide ([Pr III 2( O 2 2- )(18C6) 2 (EG)2][OTf]4; 18C6 = 18-crown-6, EG = ethylene glycol, - OTf = - O 3 SCF 3 ; 2-Pr) was isolated and characterized by single-crystal X-ray diffraction, Raman spectroscopy, and NMR spectroscopy. 2-Pr displays high thermal stability (120 °C, 50 mTorr), is protonated by mild organic acids [pK a1 (MeOH) = 5.09 ± 0.23], and engages in electrophilic (e.g., oxygen atom transfer) and nucleophilic (e.g., phosphate-ester cleavage) reactivity. Our mechanistic studies reveal that the rate of oxygen reduction is dictated by metal-ion accessibility, rather than Lewis acidity, and suggest new opportunities for differentiated reactivity of redox-inactive metal ions by leveraging weak metal–ligand binding events preceding electron transfer.

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

Status of the WPEC subgroup 46 - Efficient and effective use of integral experiments for nuclear data validation

The present paper summarizes the current status of the activities of the on-going WPEC subgroup 46 which was the last significant initiative of Massimo before he passed away. The goal of WPEC/SG46 is to define, test and document a methodology to provide unambiguous feedback to the nuclear data evaluator community, based on the joint use of integral experiments and data assimilation techniques. Part of this effort resulted in a renewed analysis of the original Target Accuracy Requirement (TAR) exercise of WPEC/SG26, by adding more diverse nuclear systems and parameters, including reaction channels correlations and a coarser energy group structure; and making use of the progress made in the most recent evaluations in terms of both nuclear data and covariance matrices. The preliminary outcomes of the updated TAR exercise, documented by various groups worldwide are clear already: uncertainty reductions are required for many nuclide-reaction pairs in a variety of energy range if the target uncertainty requirements set by the industry are to be met, especially for k{sub eff}. The inclusion of correlations between the various reaction channels had a major impact on the magnitude of the required uncertainty reduction. Those uncertainty reductions are unlikely to be met by differential measurements alone. However, the selection of relevant integral experiment and a subsequent adjustment procedure may help meet these requirements. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

ENDF/B-VIII.0 Augmented Covariance Data The first iteration [Slides]

Nuclear data are necessary for reliable modeling and simulation of the next generation of nuclear reactors. However, the variation in the ratio of the computed to experimental values (C/E) for certain types of nuclear systems is much less than predicted by evaluated nuclear data file (ENDF)/B covariances. Figure 1 provides an example of a set of metal-plutonium-fueled, fast-spectrum (PU-MET-FAST) integral experiments from the International Criticality Safety Benchmark Evaluation Project (ICSBEP) in the Oak Ridge National Laboratory (ORNL) VALID database. The variation in the C/E values, shown with one standard deviation error bars, is 100% covered by both the SCALE and ENDB/VIII.0 covariance data. This is due to the comparisons to integral data that are essential during the evaluation process. However, the ENDF evaluations represent uncertainties and correlations in differential data only; they do not reflect the impact of the comparison to integral data in the covariance evaluations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

APOLLO: a facility-scale differentiable virtual accelerator for Fermilab

As the design complexity of modern accelerators grows, there is more interest in using advanced simulations that have fast execution time or yield additional insights like gradients. The FAST/IOTA facility has been working on implementing and experimentally validating an end-to-end digital twin that is both fast and gradient-aware, allowing for rapid prototyping of new software and experiments with minimal beam time costs. Our framework integrates physics and ML codes for linac and ring simulation through a set of generic interfaces between surrogate and physics-based sections. To reproduce device inputs and outputs, system state is exposed as a deterministic discrete event simulator. Because Fermilab is undergoing control system transition, both EPICS and ACNET frontends are supported. Recently, we have begun transitioning to a new community lattice standard, PALS, as well as developing standardized infrastructure for data ingest and normalization to prepare for model calibration during FAST proton injector commissioning. We discuss implementation details as well as challenges, and future plans to extend modelling to main complex proton accelerators like PIPII and Booster.

Kuklev, Nikita [Fermilab]↗
NASA NTRS1998

Differential Absorption Lidar (DIAL) Measurements of Atmospheric Water Vapor Utilizing Robotic Aircraft

A new unpiloted air vehicle (UAV) based water vapor DIAL system will be described. This system is expected to offer lower operating costs, longer test duration and severe weather capabilities. A new high-efficiency, compact, light weight, diode-pumped, tunable Cr:LiSAF laser will be developed to meet the UAV payload weight and size limitations and its constraints in cooling capacity, physical size and payload. Similarly, a new receiver system using a single mirror telescope and an avalanche photo diode (APD) will be developed. Projected UAV parameters are expected to allow operation at altitudes up to 20 km, endurance of 24 hrs and speed of 400 km/hr. At these conditions measurements of water vapor at an uncertainty of 2-10% with a vertical resolution of 200 m and horizontal resolution of 10 km will be possible.

Hoang, Ngoc↗