Search NASA⌕ Search

SEARCH · Search NASA

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Deuteration Effects on the Physical and Optoelectronic Properties of Donor–Acceptor Conjugated Polymers

The significant differences in scattering cross sections between deuterium and protium are unique to neutron scattering techniques and have been a long-standing area of interest within the neutron scattering community. Researchers have explored selective deuteration to manipulate scattering contrast in soft matter systems, leading to the widespread use of deuterium labeling in materials development. As deuteration changes the atomic mass, it alters physical properties such as molecular volume, polarizability, and polarity, which in turn may affect noncovalent interactions and crystal ordering. Despite previous studies, there remains a limited understanding of how deuteration impacts donor–acceptor (DA) conjugated polymers. To address this, we synthesized deuterated DPP polymers and systematically investigated the effects of side-chain deuteration on their thermal stability, crystal packing, morphology, and optoelectronic properties. We found that deuteration increased the melting and crystallization temperatures of DPP polymers, although it did not significantly alter their morphology, molecular packing, or charge mobility. These properties were assessed by using atomic force microscopy (AFM), X-ray scattering, and thin-film transistor device measurements, respectively, for DPP polymers. Our work shows that deuterium labeling could be a powerful method for controlling scattering length density, enabling neutrons to study the structure and dynamics of conjugated polymers without impacting their electronic performance.

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

Learning Constitutive Relations From Soil Moisture Data via Physically Constrained Neural Networks

Abstract The constitutive relations of the Richardson‐Richards equation encode the macroscopic properties of soil water retention and conductivity. These soil hydraulic functions are commonly represented by models with a handful of parameters. The limited degrees of freedom of such soil hydraulic models constrain our ability to extract soil hydraulic properties from soil moisture data via inverse modeling. We present a new free‐form approach to learning the constitutive relations using physically constrained neural networks. We implemented the inverse modeling framework in a differentiable modeling framework, JAX, to ensure scalability and extensibility. For efficient gradient computations, we implemented implicit differentiation through a nonlinear solver for the Richardson‐Richards equation. We tested the framework against synthetic noisy data and demonstrated its robustness against varying magnitudes of noise and degrees of freedom of the neural networks. We applied the framework to soil moisture data from an upward infiltration experiment and demonstrated that the neural network‐based approach was better fitted to the experimental data than a parametric model and that the framework can learn the constitutive relations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Data-Driven Invertible Neural Surrogates of Atmospheric Transmission

We present Data-Driven Invertible Neural Surrogates of Atmospheric transmission, or DINSAT. DINSAT is a novel framework for inferring an atmospheric transmission profile from a spectral scene. This framework leverages a lightweight, physics-based simulator that is automatically tuned -- by virtue of autodifferentiation and differentiable programming -- to construct a surrogate atmospheric profile to model the observed data. The framework has utility in (i) performing atmospheric correction, (ii) recasting spectral data between various modalities (e.g. radiance and reflectance at the surface and at the sensor), and (iii) inferring atmospheric transmission profiles, such as absorbing bands and their relative magnitudes. We demonstrate the utility of these methods by performing a canonical atmospheric correction task for the purposes of further analysis - in this case, target detection within a scene.

Koch, James V.↗
DOE OSTI2024

Machine-learning based model reduction for partial differential equations

We develop a novel synergistic approach between model reduction and machine learning. The specific goal of this project is to aid in the construction of reduced order models for basis functions that are custom-made to represent the solution of partial differential equations. Partial differential equations (PDEs) are one of the main mathematical tools for describing physical phenomena. However, due to either efficiency or necessity, for many real-world problems, we are interested in constructing reduced order models (ROMs) which focus only on the explicit computation of subsets of the active spatio-temporal scales in the problem, while treating the interaction with the rest of the scales approximately. The task of accurate representation of such interactions (usually called memory terms) constitutes a vast area of research known as model reduction. PI Stinis has significant expertise in the construction of ROMs for complex systems. In addition, in recent work with the project key participant Qadeer, they have utilized machine learning to acquire custom-made basis functions (CBFs) to expand the solutions of PDEs. In the proposed work, we will merge the two concepts by constructing ROMs for subsets of the CBFs needed to represent the solution of a PDE. Specifically, we will use the Mori-Zwanzig model reduction formalism to construct ROMs for subsets of CBFs for nonlinear PDEs of various complexity, as well as investigate the usage of CBFs in the spectral vanishing viscosity method for problems that can form shocks in finite time. The outcome of the research is aimed to be proof-of-concept about a novel synergistic approach between model reduction and machine learning, thus advancing the field of scientific machine learning. Such a capability will benefit the efficient modeling of physical systems appearing in various areas of interest to the DOE.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1979

Origin of cosmic magnetic fields

Reasons for describing cosmic hydromagnetic phenomena by introducing electric currents are summarized demonstrating that change from the traditional magnetic field description to a current description provides new aspects of cosmic electrodynamics. It is shown that the kink-instability of electric currents is the basic physical phenomenon responsible for magnetic flux generation in cosmic physics. A model which starts with a poloidal magnetic field is discussed; differential motions change the field configuration so that kinetic energy is transferred into magnetic energy of induced toroidal fields; the currents associated with these become unstable when the toroidal energy exceeds the poloidal energy resulting in an amplification of the original poloidal field.

Alfven, H.↗
NASA NTRS1984

The influence of a wind tunnel on helicopter rotational noise: Formulation of analysis

An analytical model is discussed that can be used to examine the effects of wind tunnel walls on helicopter rotational noise. A complete physical model of an acoustic source in a wind tunnel is described and a simplified version is then developed. This simplified model retains the important physical processes involved, yet it is more amenable to analysis. The simplified physical model is then modeled as a mathematical problem. An inhomogeneous partial differential equation with mixed boundary conditions is set up and then transformed into an integral equation. Details of generating a suitable Green's function and integral equation are included and the equation is discussed and also given for a two-dimensional case.

Mosher, M.↗
NASA NTRS1986

Alexandrite laser source for atmospheric lidar measurements

During the past years, there has been a marked increase in interest in the applications of vibronic solid state lasers to meteorology and atmospheric physics. Two airborne lidar programs are now under development in France. The differential absorption lidar (DIAL) method with vibronic solid state lasers is very attractive for water vapor, temperature and pressure measurements. Alexandrite laser and titanium-sapphire are both suitable for these applications. However, only alexandrite rods are commercially available. The requirements on the laser source for airborne dial applications are two fold: (1) a restriction on laser linewidth and a requirement on stability and tunability with a good spectral purity; and (2) a requirement on the time separation between the two pulses. These constraints are summarized.

Pelon, J.↗
NASA NTRS1986

An analysis of morphologic variation in simple lunar craters

Previous photogeologic classification of lunar craters have differentiated them into various morphologic types. In an attempt to identify the physical factors responsible for these morphologic variations, this study examined 44 fresh, simple lunar craters between 1.1 and 14 km in diameter on the maria. Photogeologic analysis was performed using Apollo stereo panoramic photography, while a quantitative analysis of cross-sectional shape was done using Fourier analysis. These analyses showed that the sample could not be split into two distinct types (flat- and round-floored) as had been done previously. The two crater types at best represent end-members of a distribution of morphologies. Some correlation was found between the degree of relative degradation and floor roundness, implying that flat floors are primary and round floors are degradational features. No difference was found between the average depth-to-diameter ratio of flat- and round-floored craters, indicating that the factor of two variation in depth-to-diameter ratio in this sample of fresh craters is a primary morphologic feature and that minor degradation of fresh craters has no significant effect on average crater depth.

Ravine, M. A.↗
NASA NTRS1998

On the Circulation Manifold for Two Adjacent Lifting Sections

The circulation functional relative to two adjacent lifting sections is studied for two cases. In the first case we consider two adjacent circles. The circulation is computed as a function of the displacement of the secondary circle along the axis joining the two centers and of the angle of attack of the secondary circle, The gradient of such functional is computed by deriving a set of elliptic functions with respect both to their argument and to their Period. In the second case studied, we considered a wing-flap configuration. The circulation is computed by some implicit mappings, whose differentials with respect to the variation of the geometrical configuration in the physical space are found by divided differences. Configurations giving rise to local maxima and minima in the circulation manifold are presented.

Zannetti, Luca↗
NASA NTRS1998

A System Concept for the Advanced Post-TRMM Rainfall Profiling Radars

Atmospheric latent heating field is fundamental to all modes of atmospheric circulation and upper mixed layer circulations of the ocean. The key to understanding the atmospheric heating process is understanding how and where precipitation occurs. The principal atmospheric processes which link precipitation to atmospheric circulation include: (1) convective mass fluxes in the form of updrafts and downdrafts; (2) microphysical. nucleation and growth of hydrometeors; and (3) latent heating through dynamical controls on the gravitation-driven vertical mass flux of precipitation. It is well-known that surface and near-surface rainfall are two of the key forcing functions on a number of geophysical parameters at the surface-air interface. Over ocean, rainfall variation contributes to the redistribution of water salinity, sea surface temperature, fresh water supply, and marine biology and eco-system. Over land, rainfall plays a significant role in rainforest ecology and chemistry, land hydrology and surface runoff. Precipitation has also been closely linked to a number of atmospheric anomalies and natural hazards that occur at various time scales, including hurricanes, cyclones, tropical depressions, flash floods, droughts, and most noticeable of all, the El Ninos. From this point of view, the significance of global atmospheric precipitation has gone far beyond the science arena - it has a far-reaching impact on human's socio-economic well-being and sustenance. These and many other science applications require the knowledge of, in a global basis, the vertical rain structures, including vertical motion, rain intensity, differentiation of the precipitating hydrometeors' phase state, and the classification of mesoscale physical structure of the rain systems. The only direct means to obtain such information is the use of a spaceborne profiling radar. It is important to mention that the Tropical Rainfall Measuring Mission (TRMM) have made a great stride forward towards this ultimate goal. The Precipitation Radar (PR) aboard the TRMM satellite is the first ever spaceborne radar dedicated to three-dimensional, global precipitation measurements over the tropics and the subtropics, as well as the detailed synopsis of a wide range of tropical rain storm systems. In only twelve months since launch, the PR, together with other science instruments abroad the satellite have already provided unprecedented insights into the rainfall systems. It is anticipated the a lot more exciting and important rain observations would be made by TRMM throughout its mission duration. While TRMM has provided invaluable data to the user community, it is only the first step towards advancing our knowledge on rain processes and its contributions to climate variability. It is envisioned that a TRMM follow-on mission is needed in such a way to capitalize on the pioneering information provided by TRMM, and its instrument capability must be extended beyond TRMM in such a way to fully address the key science questions from microphysical to climatic time scale. In fact, a number of new and innovative mission concepts have recently put forth for this purpose. Almost all of these new concepts have suggested the utility of a more advanced, high-resolution, Doppler-enabled, vertical profiling radar that can provide multi-parameter observations of precipitation. In this paper, a system concept for a second- gene ration precipitation radar (PR-2) which addresses the above requirements will be described.

Im, Eastwood↗
NASA NTRS2004

Estimation of Length-Scales in Soils by MRI

Soil can be best described as an unconsolidated granular media that forms porous structure. The present macroscopic theory of water transport in porous media rests upon the continuum hypothesis that the physical properties of porous media can be associated with continuous, twice-differentiable field variables whose spatial domain is a set of centroids of Representative Elementary Volume (REV) elements. MRI is an ideal technique to estimate various length-scales in porous media. A 0.267 T permanent magnet at NASA GRC was used for this study. A 2D or 3D spatially-resolved porosity distribution were obtained from the NMR signal strength from each voxel and the spin-lattice relaxation time. A classical spin-warp imaging with Multiple Spin Echos (MSE) was used to evaluate proton density in each voxel. Initial resolution of 256 x 256 was subsequently reduced by averaging neighboring voxels and the porosity convergence was observed. A number of engineered "space candidate" soils such as Isolite(trademark), Zeoponics(trademark), Turface(trademark), and Profile(trademark) were used. Glass beads in the size range between 50 microns to 2 mm were used as well. Initial results with saturated porous samples have shown a good estimate of the average porosity consistent with the gravimetric porosity measurement results. For Profile(trademark) samples with particle sizes ranging between 0.25 to 1 mm and characteristic interparticle pore size of 100 microns the characteristic Darcy scale was estimated to be about delta(sub REV) = 10 mm. Glass beads porosity show clear convergence toward a definite REV which stays constant throughout homogeneous sample. Additional information is included in the original extended abstract.

N E Daidzic↗
NASA NTRS2010

Turbulence and the Stabilization Principle

Further results of research, reported in several previous NASA Tech Briefs articles, were obtained on a mathematical formalism for postinstability motions of a dynamical system characterized by exponential divergences of trajectories leading to chaos (including turbulence). To recapitulate: Fictitious control forces are introduced to couple the dynamical equations with a Liouville equation that describes the evolution of the probability density of errors in initial conditions. These forces create a powerful terminal attractor in probability space that corresponds to occurrence of a target trajectory with probability one. The effect in ordinary perceived three-dimensional space is to suppress exponential divergences of neighboring trajectories without affecting the target trajectory. Con sequently, the postinstability motion is represented by a set of functions describing the evolution of such statistical quantities as expectations and higher moments, and this representation is stable. The previously reported findings are analyzed from the perspective of the authors Stabilization Principle, according to which (1) stability is recognized as an attribute of mathematical formalism rather than of underlying physics and (2) a dynamical system that appears unstable when modeled by differentiable functions only can be rendered stable by modifying the dynamical equations to incorporate intrinsic stochasticity.

Zak, Michail↗
NASA NTRS2017

Investigation of Global Particulate Nitrate from the AeroCom Phase III Experiment

An assessment of global particulate nitrate and ammonium aerosol based on simulations from nine models participating in the Aerosol Comparisons between Observations and Models (AeroCom) phase III study is presented. A budget analysis was conducted to understand the typical magnitude, distribution, and diversity of the aerosols and their precursors among the models. To gain confidence regarding model performance, the model results were evaluated with various observations globally, including ground station measurements over North America, Europe, and east Asia for tracer concentrations and dry and wet depositions, as well as with aircraft measurements in the Northern Hemisphere mid-to-high latitudes for tracer vertical distributions. Given the unique chemical and physical features of the nitrate occurrence, we further investigated the similarity and differentiation among the models by examining (1) the pH-dependent NH3 wet deposition; (2) the nitrate formation via heterogeneous chemistry on the surface of dust and sea salt particles or thermodynamic equilibrium calculation including dust and sea salt ions; and (3) the nitrate coarse-mode fraction (i.e., coarse/total). It is found that HNO3, which is simulated explicitly based on full O3-HOx-NOx-aerosol chemistry by all models, differs by up to a factor of 9 among the models in its global tropospheric burden. This partially contributes to a large difference in NO3(-), whose atmospheric burden differs by up to a factor of 13. The atmospheric burdens of NH3 and NHC 4 differ by 17 and 4, respectively. Analyses at the process level show that the large diversity in atmospheric burdens of NO3(-), NH3, and NHC4(+) is also related to deposition processes. Wet deposition seems to be the dominant process in determining the diversity in NH3 and NHC 4 lifetimes. It is critical to correctly account for contributions of heterogeneous chemical production of nitrate on dust and sea salt, because this process overwhelmingly controls atmospheric nitrate production (typically greater than 80 %) and determines the coarse- and fine-mode distribution of nitrate aerosol.

Ammonium aerosol based on simulations↗
NASA NTRS1980

An analytical transformation technique for generating uniformly spaced computational mesh

An analytical transformation method which can map arbitrary physical coordinate grid distribution into desired computational coordinate with uniform grid distribution is derived. The transformation function and its higher derivatives are differentiable. Salient features include; (1) precise control of grid sizes; (2) more than one location of clustered grids; (3) exact positioning of particular computational nodes in the physical plane; and (4) ensuring several patches of uniformly spaced grids in the physical plane for the higher accuracies (such as at the boundaries), etc., while keeping the variation of grid spacing continuous to avoid numerical instability.

Oh, Y. H.↗
NASA NTRS1987

On determining important aspects of mathematical models: Application to problems in physics and chemistry

The use of parametric and functional gradient sensitivity analysis techniques is considered for models described by partial differential equations. By interchanging appropriate dependent and independent variables, questions of inverse sensitivity may be addressed to gain insight into the inversion of observational data for parameter and function identification in mathematical models. It may be argued that the presence of a subset of dominantly strong coupled dependent variables will result in the overall system sensitivity behavior collapsing into a simple set of scaling and self similarity relations amongst elements of the entire matrix of sensitivity coefficients. These general tools are generic in nature, but herein their application to problems arising in selected areas of physics and chemistry is presented.

Rabitz, Herschel↗
DOE OSTI2026

Combination and interpretation of differential Higgs boson production cross sections in proton-proton collisions at $ \sqrt{s}=13 $ TeV

Precision measurements of Higgs boson differential production cross sections are a key tool to probe the properties of the Higgs boson and test the standard model. New physics can affect both Higgs boson production and decay, leading to deviations from the distributions that are expected in the standard model. In this paper, combined measurements of differential spectra in a fiducial region matching the experimental selections are performed, based on analyses of four Higgs boson decay channels (γγ, ZZ$^{(*)}$, WW$^{(*)}$, and ττ) using proton-proton collision data recorded with the CMS detector at $ \sqrt{s}=13 $ TeV, corresponding to an integrated luminosity of 138 fb$^{−1}$. The differential measurements are extrapolated to the full phase space and combined to provide the differential spectra. A measurement of the total Higgs boson production cross section is also performed using the γγ and ZZ decay channels, with a result of $ {53.4}_{-2.9}^{+2.9}{\left(\textrm{stat}\right)}_{-1.8}^{+1.9}\left(\textrm{syst}\right) $ pb, consistent with the standard model prediction of 55.6 ± 2.5 pb. The fiducial measurements are used to compute limits on Higgs boson couplings using the κ-framework and the SM effective field theory.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Spatiotemporal control of laser intensity using differentiable programming

Optical techniques for spatiotemporal control can produce laser pulses with custom amplitude, phase, or polarization structure. In nonlinear optics and plasma physics, the use of structured pulses typically follows a forward design approach, in which the efficacy of a known structure is analyzed for a particular application. Inverse approaches, in contrast, enable the discovery of new structures with the potential for superior performance. Here, an implementation of the unidirectional pulse propagation equation that supports automatic differentiation is combined with gradient-based optimization to design structured pulses with features that are advantageous for a range of nonlinear optical and plasma-based applications: (1) a longitudinally uniform intensity over an extended region, (2) a superluminal intensity peak that travels many Rayleigh ranges with constant duration, spot size, and amplitude, and (3) a laser pulse that ionizes a gas to form a uniform column of plasma. In the final case, optimizing the full spatiotemporal structure improves the performance by a factor of 15 compared to optimizing only spatial or only temporal structure, highlighting the advantage of spatiotemporal control.

automatic differentiation↗
NASA NTRS1998

Simulation of Ultra-Small MOSFETs Using a 2-D Quantum-Corrected Drift-Diffusion Model

The continued down-scaling of electronic devices, in particular the commercially dominant MOSFET, will force a fundamental change in the process of new electronics technology development in the next five to ten years. The cost of developing new technology generations is soaring along with the price of new fabrication facilities, even as competitive pressure intensifies to bring this new technology to market faster than ever before. To reduce cost and time to market, device simulation must become a more fundamental, indeed dominant, part of the technology development cycle. In order to produce these benefits, simulation accuracy must improve markedly. At the same time, device physics will become more complex, with the rapid increase in various small-geometry and quantum effects. This work describes both an approach to device simulator development and a physical model which advance the effort to meet the tremendous electronic device simulation challenge described above. The device simulation approach is to specify the physical model at a high level to a general-purpose (but highly efficient) partial differential equation solver (in this case PROPHET, developed by Lucent Technologies), which then simulates the model in 1-D, 2-D, or 3-D for a specified device and test regime. This approach allows for the rapid investigation of a wide range of device models and effects, which is certainly essential for device simulation to catch up with, and then stay ahead of, electronic device technology of the present and future. The physical device model used in this work is the density-gradient (DG) quantum correction to the drift-diffusion model [Ancona, Phys. Rev. B 35(5), 7959 (1987)]. This model adds tunneling and quantum smoothing of carrier density profiles to the drift-diffusion model. We used the DG model in 1-D and 2-D (for the first time) to simulate both bipolar and unipolar devices. Simulations of heavily-doped, short-base diodes indicated that the DG quantum corrections do not have a large effect on the IN characteristics of electronic devices without heteroj unction s. On the other hand, ultra-small MOSFETs certainly exhibit important quantum effects that the DG model will include: quantum repulsion of the inversion and gate charges from the oxide interfaces, and quantum tunneling through thin gate oxides. We present initial results of 2-D DG simulations of ultra-small MOSFETs. Subtle but important issues involving the specification of the model, boundary conditions, and interface constraints for DG simulation of MOSFETs will also be illuminated.

Biegal, Bryan A.↗