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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Maximum likelihood identification of aircraft parameters with unsteady aerodynamic modelling

A simplified aerodynamic force model based on the physical principle of Prandtl's lifting line theory and trailing vortex concept has been developed to account for unsteady aerodynamic effects in aircraft dynamics. Longitudinal equations of motion have been modified to include these effects. The presence of convolution integrals in the modified equations of motion led to a frequency domain analysis utilizing Fourier transforms. This reduces the integro-differential equations to relatively simple algebraic equations, thereby reducing computation time significantly. A parameter extraction program based on the maximum likelihood estimation technique is developed in the frequency domain. The extraction algorithm contains a new scheme for obtaining sensitivity functions by using numerical differentiation. The paper concludes with examples using computer generated and real flight data

Keskar, D. A.↗
NASA NTRS2005

Coupled Gravity and Elevation Measurements of Ice Sheet Mass Change

We measured surface gravity and position at ten locations about two glaciological measurement networks located on the South-central Greenland Ice during June 2004. Six of the individual sites of the first network were occupied the previous year. At the repeat sites we were able to measure annual accumulation rate and surface displacement by referencing measurements to aluminum poles left in the firn the previous year. We occupied 4 additional sites at a second measurement network for the first time since initial observations were last made at the network in 1981. At each individual site, we operated a GPS unit for 90 minutes - the unit was operated simultaneously with a base station unit in Sondrestrom Fjord so as to enable differential, post-processing of the data. We installed an aluminum, accumulation-rate-pole at each site. The base section of the pole also served as the mount for the GPS antenna. A new, Scintrex gravimeter was used at each site and relative gravity measurements were tied to the network of absolute gravity stations in Sondrestrom. We measured snow physical properties in two shallow pits. This report summarizes our observations and data analysis.

Jezek, K. C.↗
NASA NTRS2020

Extending Science from Lunar Laser Ranging

The Lunar Laser Ranging (LLR) experiment has accumulated 50 years of range data of improving accuracy from ground stations to the laser retroreflector arrays (LRAs) on the lunar surface. The upcoming decade offers several opportunities to break new ground in data precision through the deployment of the next generation of single corner-cube lunar retroreflectors and active laser transponders. This is likely to expand the LLR station network. Lunar dynamical models and analysis tools have the potential to improve and fully exploit the long temporal baseline and precision allowed by millimetric LLR data. Some of the model limitations are outlined for future efforts. Differential observation techniques will help mitigate some of the primary limiting factors and reach unprecedented accuracy. Such observations and techniques may enable the detection of several subtle signatures required to understand the dynamics of the Earth- Moon system and the deep lunar interior. LLR model improvements would impact multi- disciplinary fields that include lunar and planetary science, Earth science, fundamental physics, celestial mechanics and ephemerides.

Vishnu Viswanathan↗
DOE OSTI2025

$\mathrm{SageNet}$: Fast Neural Network Emulation of the Stiff-amplified Gravitational Waves from Inflation

Accurate modeling of the inflationary gravitational waves (GWs) requires time-consuming, iterative numerical integrations of differential equations to take into account their backreaction on the expansion history. To improve computational efficiency while preserving accuracy, we present the Stiff-amplified Gravitational-wave Emulator Network (SageNet), a deep learning framework designed to replace conventional numerical solvers (code available at https://github.com/YifangLuo/SageNet). SageNet employs a long short-term memory architecture to emulate the present-day energy density spectrum of the inflationary GWs with possible stiff amplification, Ω GW (f). Trained on a data set of 25,689 numerically generated solutions, SageNet allows accurate reconstructions of Ω GW (f) and generalizes well to a wide range of cosmological parameters; 90.9% of the test emulations with randomly distributed parameters exhibit errors of under 4%. In addition, SageNet demonstrates its ability to learn and reproduce the artificial, adaptive sampling patterns in numerical calculations, which implement denser sampling of frequencies around changes in spectral indices in Ω GW (f). The dual capability of learning both physical and artificial features of the numerical GW spectra establishes SageNet as a robust alternative to exact numerical methods. Finally, our benchmark tests show that SageNet reduces the computation time from tens of seconds to milliseconds, achieving a speedup of ∼10 4 times over standard CPU-based numerical solvers with the potential for further acceleration on GPU hardware. These capabilities make SageNet a powerful tool for accelerating Bayesian inference procedures for extended cosmological models. In a broad sense, the SageNet framework offers a fast, accurate, and generalizable solution to modeling cosmological observables whose theoretical predictions demand costly differential equation solvers.

Astronomy data modeling↗
DOE OSTI2025

NEML2: An efficient and modular multiphysics constitutive modeling library for hybrid computing environments

This paper presents NEML2, an open-source, high-performance library developed for constitutive material modeling, designed to support the flexible and modular development of models for complex material behavior. Building on the foundational structure of its predecessor, NEML, the NEML2 library introduces significant improvements, including enhanced vectorization, automatic differentiation, and seamless integration with PyTorch, facilitating the application of machine learning techniques in material simulations. NEML2 provides a C++ backend with Python bindings, enabling users to create custom material models that can be executed efficiently on both CPU and GPU platforms. The library also supports coupling with Multiphysics simulation frameworks like MOOSE, making it suitable for realistic simulations involving coupled physical processes. Rigorous quality assurance through unit and regression testing ensures the reliability of results, while the extensible, user-friendly design encourages collaboration and reproducibility across the scientific community. This paper provides an overview of NEML2’s architecture, core features, and applications, highlighting its impact on accelerating material qualification and advancing computational methods in materials science.

GPU↗
DOE OSTI2022

Enhancing high-fidelity nonlinear solver with reduced order model

Abstract We propose the use of reduced order modeling (ROM) to reduce the computational cost and improve the convergence rate of nonlinear solvers of full order models (FOM) for solving partial differential equations. In this study, a novel ROM-assisted approach is developed to improve the computational efficiency of FOM nonlinear solvers by using ROM’s prediction as an initial guess. We hypothesize that the nonlinear solver will take fewer steps to the converged solutions with an initial guess that is closer to the real solutions. To evaluate our approach, four physical problems with varying degrees of nonlinearity in flow and mechanics have been tested: Richards’ equation of water flow in heterogeneous porous media, a contact problem in a hyperelastic material, two-phase flow in layered porous media, and fracture propagation in a homogeneous material. Overall, our approach maintains the FOM’s accuracy while speeding up nonlinear solver by 18–73% (through suitable ROM-assisted FOMs). More importantly, the proximity of ROM’s prediction to the solution space leads to the improved convergence of FOMs that would have otherwise diverged with default initial guesses. We demonstrate that the ROM’s accuracy can impact the computational efficiency with more accurate ROM solutions, resulting in a better cost reduction. We also illustrate that this approach could be used in many FOM discretizations (e.g., finite volume, finite element, or a combination of those). Since our ROMs are data-driven and non-intrusive, the proposed procedure can easily lend itself to any nonlinear physics-based problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

W$^+$W$^-$ boson pair production in proton-proton collisions at $\sqrt{s} =$ 13 TeV

A measurement of the W+W- boson pair production cross section in proton-proton collisions at s=13 TeV is presented. The data used in this study are collected with the CMS detector at the CERN LHC and correspond to an integrated luminosity of 35.9 fb-1. The W+W- candidate events are selected by requiring two oppositely charged leptons (electrons or muons). Two methods for reducing background contributions are employed. In the first one, a sequence of requirements on kinematic quantities is applied allowing a measurement of the total production cross section, 117.6±6.8 pb, which agrees well with the theoretical prediction. Fiducial cross sections are also reported for events with zero or one jet, and the change in the zero-jet fiducial cross section with the jet transverse momentum threshold is measured. Normalized differential cross sections are reported within the fiducial region. A second method for suppressing background contributions employs two random forest classifiers. The analysis based on this method includes a measurement of the total production cross section and also a measurement of the normalized jet multiplicity distribution in W+W- events. Finally, a dilepton invariant mass distribution is used to probe for physics beyond the standard model in the context of an effective field theory, and constraints on the presence of dimension-6 operators are derived.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Power Flow Geometry and Approximation

Here, the power flow equations are important in numerous power systems problems of practical interest which consider alternating current power flow (ACPF) physics. Perhaps the most well studied being the alternating current optimal power flow problem (ACOPF), seeking to optimize the operation of an electric power system. Due to their non-linearity, problems which include the power flow equations are typically challenging, particularly in optimization. Interestingly, the set of solutions to the power flow equations forms a smooth manifold. As a result, differential geometry can be used to describe and analyze this set of equations. This approach has proven effective in several engineering applications (e.g., solving ACOPF and analyzing the solution space boundary). Central to the success of this approach is an understanding of the power flow manifold's geometry. In this work, we develop the geometric and topological properties of this manifold using concepts from differential geometry. After demonstrating the convenience of this manifold's representation as a function's graph, computational methods are emphasized: we develop retractions, error bounds for linear approximation, and formulas for evaluating the Riemannian metric (including associated objects such as geodesics and the curvature tensor). Scalar curvature and the second fundamental form play a new role in quantifying the quality of linear approximations, like the popular direct current approximation. All functions are implemented in Julia and available in an online repository. Proofs are included for completeness.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Inverse Analysis with Variational Autoencoders: A Comparison of Shallow and Deep Networks

Inverse problems are applied to determine unknown properties by matching observational data with a physical model that often takes many parameters as input. To overcome the underconstrained nature of inverse problems and achieve good performance, an approach is presented involving regularization with a technique known as a variational autoencoder (VAE), which is trained to map a high-dimensional parameter space with a complex structure to a low-dimensional latent space with a simple structure. We apply this approach to unconditioned realizations of the parameters (heterogeneous hydraulic fields) for a hydrogeological inverse problem. Two types of hydraulic conductivity fields are used to evaluate the characterization for the different levels of heterogeneity complexity of the physical inputs. This approach keeps the computational cost of generating the training data low. The reason is unconditioned realizations neither rely on the observational data used to perform the inverse analysis nor require any groundwater flow model forward runs. In addition, this approach applies regularization on a low-dimensional latent space from the VAE and increases optimization efficiency through automatic differentiation. Furthermore, two different neural network (NN) structures are tested for their utility in using the VAE for inverse analysis. The performance of a deep, convolutional neural network strongly depends on the dimensionality of the latent space, which requires tuning. In contrast, a shallow, dense neural network provides consistently accurate characterization without tuning. Furthermore, our approach evaluates the advantages of the shallow, dense neural network over the deep, convolutional one and enables future application to a wide range of inverse problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Recycling of Pretreated Polyolefin-Based Ocean-Bound Plastic Waste by Incorporating Clay and Rubber

Plastic waste found in oceans has become a major concern because of its impact on marine organisms and human health. There is significant global interest in recycling these materials, but their reclamation, sorting, cleaning, and reprocessing, along with the degradation that occurs in the natural environment, all make it difficult to achieve high quality recycled resins from ocean plastic waste. To mitigate these limitations, various additives including clay and rubber were explored. In this study, we compounded different types of ocean-bound (o-HDPE and o-PP) and virgin polymers (v-LDPE and v-PS) with various additives including a functionalized clay, styrene-multi-block-copolymer (SMB), and ethylene-propylene-based rubber (EPR). Physical observation showed that all blends containing PS were brittle due to the weak interfaces between the polyolefin regions and the PS domains within the polymer blend matrix. Blends containing clay showed rough surfaces and brittleness because of the non-uniform distribution of clay particles in the polymer matrix. To evaluate the properties and compatibility of the blends, characterizations using differential scanning calorimetry (DSC), scanning electron microscopy (SEM), and small-amplitude oscillatory shear (SAOS) rheology were carried out. The polymer blend (v-LDPE, o-HDPE, o-PP) containing EPR showed improved elasticity. Incorporating additives such as rubber could improve the mechanical properties of polymer blends for recycling purposes.

42 ENGINEERING↗
NASA NTRS1973

Marshall system for aerospace system simulation (MARSYAS), user's manual

The capabilities of the Marshall system for aerospace system simulation (MARSYAS) and how to use it are described. MARSYAS is a software system that allows easy setup and control of the simulation of the dynamics of large physical systems on a digital computer. The physical systems are modeled in the form of block diagrams or equations. The blocks can have multiple inputs and multiple outputs, and they can be nested to form hierarchies. The block diagrams can contain transfer functions, nonlinear and logical functions, equations, analog computer elements and FORTRAN programs. The input format of the equations can be combinations of nonlinear, time-varying differential equations and algebraic equations in their original format. MARSYAS could also serve as a storage and retrieval system for models as a basis for a model configuration control system on a central time-shared computer. The outputs of the simulation system can be not only time-responses but also other analysis data such as frequency response, power spectrum and stability parameters. The MARSYAS translator is written in FORTRAN running on the Univac 1108 computer under the EXEC 8 operating system.

Ventre, A.↗
NASA NTRS1981

Design and analysis of a torsion braid pendulum displacement transducer

The dynamic properties at various temperatures of braids impregnated with polymer can be measured by using the braid as the suspension of a torsion pendulum. This report describes the electronic and mechanical design of a torsional braid pendulum displacement transducer which is an advance in the state of the art. The transducer uses a unique optical design consisting of refracting quartz windows used in conjunction with a differential photocell to produce a null signal. The release mechanism for initiating free torsional oscillation of the pendulum has also been improved. Analysis of the precision and accuracy of the transducer indicated that the maximum relative error in measuring torsional amplitude was approximately 0. A serious problem inherent in all instruments which use a torsional suspension was analyzed: misalignment of the physical and torsional axes of the torsional member which results in modulation of the amplitude of the free oscillation.

Rind, E.↗
NASA NTRS1987

Distortions of perceived visual directions out of pictures

Two experiments are reported examining judgments from 16 subjects who indicated the apparent direction of a photographed pointer that was rotated to different physical positions while being photographed. The photographs themselves were rotated about a vertical axis to several positions with respect to the subjects' central viewing axis. The results replicate the well-known distortion in apparent direction associated with photographed pointers positioned to project directly out of the plane of the photograph. This effect has been described by Goldstein (1979) as the 'differential rotation effect' because its magnitude is reduced as the depicted angle of the pointer becomes less orthogonal to the photograph. Analysis of the two-dimensional properties of the projected images shows that this differential rotation is related to projected angles on the surface of the photograph. This analysis may explain why circular objects often do not appear to be correctly drawn in the periphery of geometrically correct projections.

Ellis, Stephen R.↗
NASA NTRS2013

Modelling of Equilibrium Between Mantle and Core: Refractory, Volatile, and Highly Siderophile Elements

Siderophile elements have been used to constrain conditions of core formation and differentiation for the Earth, Mars and other differentiated bodies [1]. Recent models for the Earth have concluded that the mantle and core did not fully equilibrate and the siderophile element contents of the mantle can only be explained under conditions where the oxygen fugacity changes from low to high during accretion and the mantle and core do not fully equilibrate [2,3]. However these conclusions go against several physical and chemical constraints. First, calculations suggest that even with the composition of accreting material changing from reduced to oxidized over time, the fO2 defined by metal-silicate equilibrium does not change substantially, only by approximately 1 logfO2 unit [4]. An increase of more than 2 logfO2 units in mantle oxidation are required in models of [2,3]. Secondly, calculations also show that metallic impacting material will become deformed and sheared during accretion to a large body, such that it becomes emulsified to a fine scale that allows equilibrium at nearly all conditions except for possibly the length scale for giant impacts [5] (contrary to conclusions of [6]). Using new data for D(Mo) metal/silicate at high pressures, together with updated partitioning expressions for many other elements, we will show that metal-silicate equilibrium across a long span of Earth s accretion history may explain the concentrations of many siderophile elements in Earth's mantle. The modeling includes refractory elements Ni, Co, Mo, and W, as well as highly siderophile elements Au, Pd and Pt, and volatile elements Cd, In, Bi, Sb, Ge and As.

Righter, K.↗
DOE OSTI2022

How Well Does Kohn–Sham Regularizer Work for Weakly Correlated Systems?

Kohn–Sham regularizer (KSR) is a differentiable machine learning approach to finding the exchange-correlation functional in Kohn–Sham density functional theory that works for strongly correlated systems. Here we test KSR for a weak correlation. We propose spin-adapted KSR (sKSR) with trainable local, semilocal, and nonlocal approximations found by minimizing density and total energy loss. We assess the atoms-to-molecules generalizability by training on one-dimensional (1D) H, He, Li, Be, and Be 2+ and testing on 1D hydrogen chains, LiH, BeH 2 , and helium hydride complexes. The generalization error from our semilocal approximation is comparable to other differentiable approaches, but our nonlocal functional outperforms any existing machine learning functionals, predicting ground-state energies of test systems with a mean absolute error of 2.7 mH.

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

Ultrasensitive single-cell proteomics workflow identifies >1000 protein groups per mammalian cell

Here, we report on the combination of nanodroplet sample preparation, ultra-low-flow nanoLC, high-field asymmetric ion mobility spectrometry (FAIMS), and the latest-generation Orbitrap Eclipse Tribrid mass spectrometer for greatly improved single-cell proteome profiling. FAIMS effectively filtered out singly charged ions for more effective MS analysis of multiply charged peptides, resulting in an average of 1056 protein groups identified from single HeLa cells without MS1-level feature matching. This is 2.3 times more identifications than without FAIMS and a far greater level of proteome coverage for single mammalian cells than has been previously reported for a label-free study. Differential analysis of single microdissected motor neurons and interneurons from human spinal tissue indicated a similar level of proteome coverage, and the two subpopulations of cells were readily differentiated based on single-cell label-free quantification.

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

Layered supramolecular hydrogels from thioglycosides

Low molecular weight hydrogels are made of small molecules that aggregate via noncovalent interactions. Here, comprehensive characterization of the physical and chemical properties of hydrogels made from thioglycolipids of the disaccharides lactose and cellobiose with simple alkyl chains is reported. While thiolactoside hydrogels are robust, thiocellobioside gels are metastable, precipitating over time into fibrous crystals that can be entangled to create pseudo-hydrogels. Rheology confirms the viscoelastic solid nature of these hydrogels with storage moduli ranging from 10–600 kPa. Additionally, thiolactoside hydrogels are thixotropic which is a desirable property for many potential applications. Freeze-fracture electron microscopy of xerogels shows layers of stacked sheets that are entangled into networks. These structures are unique compared to the fibers or ribbons typically reported for hydrogels. Differential scanning calorimetry provides gel-to-liquid phase transition temperatures ranging from 30 to 80 °C. Prodan fluorescence spectroscopy allows assignment of phase transitions in the gels and other lyotropic phases of high concentration samples. Phase diagrams are estimated for all hydrogels at 1–10 wt% from 5 to ≥ 80 °C. In conclusion, these hydrogels represent a series of interesting materials with unique properties that make them attractive for numerous potential applications.

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

A Tailored Convolutional Neural Network for Nonlinear Manifold Learning of Computational Physics Data Using Unstructured Spatial Discretizations

In this work, we propose a nonlinear manifold learning technique based on deep convolutional autoencoders that is appropriate for model order reduction of physical systems in complex geometries. Convolutional neural networks have proven to be highly advantageous for compressing data arising from systems demonstrating a slow-decaying Kolmogorov n-width. However, these networks are restricted to data on structured meshes. Unstructured meshes are often required for performing analyses of real systems with complex geometry. Our custom graph convolution operators based on the available differential operators for a given spatial discretization effectively extend the application space of deep convolutional autoencoders to systems with arbitrarily complex geometry that are typically discretized using unstructured meshes. We propose sets of convolution operators based on the spatial derivative operators for the underlying spatial discretization, making the method particularly well suited to data arising from the solution of partial differential equations. We demonstrate the method using examples from heat transfer and fluid mechanics and show better than an order of magnitude improvement in accuracy over linear methods.

97 MATHEMATICS AND COMPUTING↗