Search NASA⌕ Search

SEARCH · Search NASA

Results for “differentiable simulations”

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 271 records · Page 15

Optimizing multigrid reduction-in-time and Parareal coarse-grid operators for linear advection

Parallel-in-time methods, such as multigrid reduction-in-time (MGRIT) and Parareal, provide an attractive option for increasing concurrency when simulating time-dependent partial differential equations (PDEs) in modern high-performance computing environments. While these techniques have been very successful for parabolic equations, it has often been observed that their performance suffers dramatically when applied to advection-dominated problems or purely hyperbolic PDEs using standard rediscretization approaches on coarse grids. In this paper, we apply MGRIT or Parareal to the constant-coefficient linear advection equation, appealing to existing convergence theory to provide insight into the typically nonscalable or even divergent behavior of these solvers for this problem. To overcome these failings, we replace rediscretization on coarse grids with improved coarse-grid operators that are computed by applying optimization techniques to approximately minimize error estimates from the convergence theory. Therefore, one of our main findings is that, in order to obtain fast convergence as for parabolic problems, coarse-grid operators should take into account the behavior of the hyperbolic problem by tracking the characteristic curves. Our approach is tested for schemes of various orders using explicit or implicit Runge–Kutta methods combined with upwind-finite-difference spatial discretizations. In all cases, we obtain scalable convergence in just a handful of iterations, with parallel tests also showing significant speed-ups over sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗

Strength-conductivity synergy via lean alloy design: A study on novel solid stir extrusion of an Al-2Cu-0.1 Nb-0.15Zr (~wt.%) alloy

Aluminum alloys find wide applications in power and transmission cables, and electrical conductors due to their advantageous properties of high electrical conductivity, and lower cost per unit. However, the lower tensile strength compared to copper poses a challenge for the mechanical stability of aluminum-based cables. In this study, we addressed this issue by adopting a process-based alloy design approach, creating a lean alloy with a dilute composition of Al-2Cu-0.1 Nb-0.15Zr (~wt.%). Here, the alloy was processed using a friction stir based novel SolidStir® (SSE) technique, followed by low-temperature aging. The resulting aluminum wire exhibited improved strength (240 MPa) while simultaneously enhancing electrical conductivity (64%IACS). The study employed computational simulations (Thermocalc) and experimental techniques (differential scanning calorimetry and transmission electron microscopy) to thoroughly investigate the microstructure. Microstructural examinations revealed that the precipitation kinetics in SSE and SSE+aged conditions played a significant role in enhancing strength and electrical conductivity. These findings demonstrate the successful achievement of high strength and high electrical conductivity in aluminum wire through a lean alloy design approach.

36 MATERIALS SCIENCE↗

Multidimensional quantitative phenotypic and molecular analysis reveals neomorphic behaviors of p53 missense mutants

Abstract Mutations in the TP53 tumor suppressor gene occur in >80% of the triple-negative or basal-like breast cancer. To test whether neomorphic functions of specific TP53 missense mutations contribute to phenotypic heterogeneity, we characterized phenotypes of non-transformed MCF10A-derived cell lines expressing the ten most common missense mutant p53 proteins and observed a wide spectrum of phenotypic changes in cell survival, resistance to apoptosis and anoikis, cell migration, invasion and 3D mammosphere architecture. The p53 mutants R248W, R273C, R248Q, and Y220C are the most aggressive while G245S and Y234C are the least, which correlates with survival rates of basal-like breast cancer patients. Interestingly, a crucial amino acid difference at one position—R273C vs. R273H—has drastic changes on cellular phenotype. RNA-Seq and ChIP-Seq analyses show distinct DNA binding properties of different p53 mutants, yielding heterogeneous transcriptomics profiles, and MD simulation provided structural basis of differential DNA binding of different p53 mutants. Integrative statistical and machine-learning-based pathway analysis on gene expression profiles with phenotype vectors across the mutant cell lines identifies quantitative association of multiple pathways including the Hippo/YAP/TAZ pathway with phenotypic aggressiveness. Further, comparative analyses of large transcriptomics datasets on breast cancer cell lines and tumors suggest that dysregulation of the Hippo/YAP/TAZ pathway plays a key role in driving the cellular phenotypes towards basal-like in the presence of more aggressive p53 mutants. Overall, our study describes distinct gain-of-function impacts on protein functions, transcriptional profiles, and cellular behaviors of different p53 missense mutants, which contribute to clinical phenotypic heterogeneity of triple-negative breast tumors.

60 APPLIED LIFE SCIENCES↗

Geometry optimization speedup through a geodesic approach to internal coordinates

We present a new geodesic-based method for geometry optimization in a basis set of redundant internal coordinates. Overall, our method updates the molecular geometry by following the geodesic generated by a displacement vector on the internal coordinate manifold, which dramatically reduces the number of steps required to converge to a minimum. Our method can be implemented in any existing optimization code, requiring only implementation of derivatives of the Wilson B-matrix and the ability to numerically solve an ordinary differential equation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Neuromorphic scaling advantages for energy-efficient random walk computations

Computing stands to be radically improved by neuromorphic computing (NMC) approaches inspired by the brain's incredible efficiency and capabilities. Most NMC research, which aims to replicate the brain's computational structure and architecture in man-made hardware, has focused on artificial intelligence; however, less explored is whether this brain-inspired hardware can provide value beyond cognitive tasks. We demonstrate that high-degree parallelism and configurability of spiking neuromorphic architectures makes them well-suited to implement random walks via discrete time Markov chains. Such random walks are useful in Monte Carlo methods, which represent a fundamental computational tool for solving a wide range of numerical computing tasks. Additionally, we show how the mathematical basis for a probabilistic solution involving a class of stochastic differential equations can leverage those simulations to provide solutions for a range of broadly applicable computational tasks. Despite being in an early development stage, we find that NMC platforms, at a sufficient scale, can drastically reduce the energy demands of high-performance computing platforms.

59 BASIC BIOLOGICAL SCIENCES↗

Janssen

Differentiable and Modular Optical Microscopy Simulations and Inversions

Mukherjee, Debangshu [Oak Ridge National Laborator↗

Optimum spectral bands for rock discrimination

Using stepwise discriminant analysis on spectral reflectance and spectral emissivity data collected by a Multispectral Scanner and Data System, mounted in an NC-130B aircraft and flown at an altitude of approximately 3 km, spectral bands were ranked as to their usefulness in separating specific rock types and rock alteration products in seven geologically diverse Utah sites. The optimum band for rock discrimination included the 1.18 to 1.30 micron interval, and the optimum combination of bands comprised the 1.18 to 1.30, 4.50 to 4.75, 0.46 to 0.50, 1.52 to 1.73, and 2.10 to 2.36 micron intervals. It is concluded that the spectral interval combination was more successful in differentiating geologic materials than either simulated Multispectral Scanner bands or simulated Thematic Mapper bands.

Siegrist, A. W.↗

Numerical observations of dynamic behavior in two-dimensional compressible convection

Using numerical simulations of the anelastic partial differential equations, convective flow in a two-dimensional polytropic atmosphere is investigated. Fixing the aspect ratio at unity and choosing the viscosity as the order parameter, the numerical solutions have indicated a variety of dynamic behavior including (1) oscillatory motion with a period of order the eddy turnover time, (2) quasi-periodic motion with a long period much longer than an eddy turnover time, and (3) at least two branches of solutions in certain regimes.

Ginet, G. P.↗

Effects of the space environment on space-based radar phased-array antenna; status and preliminary observations (LDEF Experiment A0133)

The overall objective of the experiment was to evaluate the effect of the space environment on components considered for a Space-Based Radar (SBR) Phased-Array Antenna. Of primary interest was a study of the degradation of the polyimide film Kapton (DuPont trademark), the material considered for use in the antenna plane. The most striking result of the experiment was the overall good condition of the Kapton antenna planes and Kapton tensile specimens, despite nearly six years of exposure to the space environment. This was largely attributable to the orientation of the Kapton (parallel and flush on the space end) and the stability of the Long Duration Exposure Facility (LDEF) in orbit. However, weathering of exposed Kapton surfaces was not insignificant. Results on elongation and mechanical properties of the plain and the fiberglass-reinforced Kapton are presented. Reduction in strain to failure of flight-exposed Kapton is attributed to surface defects of these specimens. Physical property testing of the materials to date reveals no significant difference between flight-exposed and control material. The second objective was to investigate the interaction between high-voltage electrodes and typical spacecraft contaminants in simulation of discharge triggering across differentially charged dielectric surfaces (spacecraft charging conditions). Electronic data acquisition and memory systems appeared to operate correctly, but very few discharges were recorded. Induced radioactivity, contamination, impacts, and orientation features of atomic oxygen erosion were observed.

J B Whiteside↗

Low-Power SOI CMOS Transceiver

The work aims at developing a low-power Silicon on Insulator Complementary Metal Oxide Semiconductor (SOI CMOS) Transceiver for deep-space communications. RF Receiver must accomplish the following tasks: (a) Select the desired radio channel and reject other radio signals, (b) Amplify the desired radio signal and translate them back to baseband, and (c) Detect and decode the information with Low BER. In order to minimize cost and achieve high level of integration, receiver architecture should use least number of external filters and passive components. It should also consume least amount of power to minimize battery cost, size, and weight. One of the most stringent requirements for deep-space communication is the low-power operation. Our study identified that two candidate architectures listed in the following meet these requirements: (1) Low-IF receiver, (2) Sub-sampling receiver. The low-IF receiver uses minimum number of external components. Compared to Zero-IF (Direct conversion) architecture, it has less severe offset and flicker noise problems. The Sub-sampling receiver amplifies the RF signal and samples it using track-and-hold Subsampling mixer. These architectures provide low-power solution for the short- range communications missions on Mars. Accomplishments to date include: (1) System-level design and simulation of a Double-Differential PSK receiver, (2) Implementation of Honeywell SOI CMOS process design kit (PDK) in Cadence design tools, (3) Design of test circuits to investigate relationships between layout techniques, geometry, and low-frequency noise in SOI CMOS, (4) Model development and verification of on-chip spiral inductors in SOI CMOS process, (5) Design/implementation of low-power low-noise amplifier (LNA) and mixer for low-IF receiver, and (6) Design/implementation of high-gain LNA for sub-sampling receiver. Our initial results show that substantial improvement in power consumption is achieved using SOI CMOS as compared to standard CMOS process. Potential advantages of SOI CMOS for deep-space communication electronics include: (1) Radiation hardness, (2) Low-power operation, and (3) System-on-Chip (SOC) solutions.

Fujikawa, Gene↗

Aeroelastic Analysis of Helicopter Rotor Blades Incorporating Anisotropic Piezoelectric Twist Actuation

A simple aeroelastic analysis of a helicopter rotor blade incorporating embedded piezoelectric fiber composite, interdigitated electrode blade twist actuators is described. The analysis consists of a linear torsion and flapwise bending model coupled with a nonlinear ONERA based unsteady aerodynamics model. A modified Galerkin procedure is performed upon the rotor blade partial differential equations of motion to develop a system of ordinary differential equations suitable for dynamics simulation using numerical integration. The twist actuation responses for three conceptual fullscale blade designs with realistic constraints on blade mass are numerically evaluated using the analysis. Numerical results indicate that useful amplitudes of nonresonant elastic twist, on the order of one to two degrees, are achievable under one-g hovering flight conditions for interdigitated electrode poling configurations. Twist actuation for the interdigitated electrode blades is also compared with the twist actuation of a conventionally poled piezoelectric fiber composite blade. Elastic twist produced using the interdigitated electrode actuators was found to be four to five times larger than that obtained with the conventionally poled actuators.

Wilkie, W. Keats↗

Template Matching Approach to Signal Prediction

A new approach to signal prediction and prognostic assessment of spacecraft health resolves an inherent difficulty in fusing sensor data with simulated data. This technique builds upon previous work that demonstrated the importance of physics-based transient models to accurate prediction of signal dynamics and system performance. While models can greatly improve predictive accuracy, they are difficult to apply in general because of variations in model type, accuracy, or intended purpose. However, virtually any flight project will have at least some modeling capability at its disposal, whether a full-blown simulation, partial physics models, dynamic look-up tables, a brassboard analogue system, or simple hand-driven calculation by a team of experts. Many models can be used to develop a predict, or an estimate of the next day s or next cycle s behavior, which is typically used for planning purposes. The fidelity of a predict varies from one project to another, depending on the complexity of the simulation (i.e. linearized or full differential equations) and the level of detail in anticipated system operation, but typically any predict cannot be adapted to changing conditions or adjusted spacecraft command execution. Applying a predict blindly, without adapting the predict to current conditions, produces mixed results at best, primarily due to mismatches between assumed execution of spacecraft activities and actual times of execution. This results in the predict becoming useless during periods of complicated behavior, exactly when the predict would be most valuable. Each spacecraft operation tends to show up as a transient in the data, and if the transients are misaligned, using the predict can actually harm forecasting performance. To address this problem, the approach here expresses the predict in terms of a baseline function superposed with one or more transient functions. These transients serve as signal templates, which can be relocated in time and space against the signal background. One then has the ability to reconstruct a signal regardless of the precise timing of the transients. During operation, one applies the actual start times of spacecraft activities as they occur, and produces a reconstructed, accurate predict in real-time. This general approach is valid under two important conditions. First, the transients themselves must be time-invariant. Second, the transients must be reasonably consistent with respect to different operating points. Both of these assumptions are generally valid, but in the case of a complicated system with numerous types of overlapping transients, this approach may not be effective. Fortunately, there tend to be few transients in spacecraft telemetry of sensor quantities or low-level health and status information because these signals rarely reflect multiple different types of operation. Furthermore, if the predict is at least reasonably close to actual operation, the shift either in time, or in the operating point when the transient occurs - is likely to be small. The proposed approach considers three ways to recognize a transient. The first and most reliable is to use a different signal that identifies operating mode - often the transient will be correlated to a change in operating mode, and this change can usually be detected positively from discrete signals in spacecraft telemetry. If there is no useful mode signal, the second option is to identify the transient template by hand, and detect the onset of the transient using a curve-fitting approach. Finally, if one elects not to choose by hand, there is an option for automatic selection. This technique has been applied to sensor data from several JPL missions and industrial applications, demonstrating an improvement in accurate prediction of future behavior and early detection of system problems. This technique is applicable to practically any time-varying, quantitative sensor measurement.

Mackey, Ryan↗

An Observational Approach to Low Lunar Frozen Orbit Design

Low lunar frozen orbits can be used to enable long-duration missions at the Moon while also minimizing the fuel usage necessary for orbit maintenance. There is no analytical expression for determining these frozen orbits, so a high-fidelity simulation approach is required. A differential correction targeting approach was used to find initial conditions that zeroed out argument of periapsis and eccentricity growth. This approach was used during the Lunar Reconnaissance Orbiter extended mission period for four different frozen orbit instances. Observations of the evolution of key orbit elements (e.g. periapsis altitude and the eccentricity vector) led to an understanding of the orbit dynamics and helped to formulate a simple observational approach to finding frozen or quasi-frozen low lunar orbits that looks at minimizing the spread in the periselene altitude evolution.

Moon Frozen Orbit LRO↗

An Observational Approach to Low Lunar Frozen Orbit Design

Low lunar frozen orbits can be used to enable long-duration missions at the Moon while also minimizing the fuel usage necessary for orbit maintenance. There is no analytical expression for determining these frozen orbits, so a high-fidelity simulation approach is required. A differential correction targeting approach was used to find initial conditions that zeroed out argument of periapsis and eccentricity growth. This approach was used during the Lunar Reconnaissance Orbiter extended mission period for four different frozen orbit instances. Observations of the evolution of key orbit elements (e.g. periapsis altitude and the eccentricity vector) led to an understanding of the orbit dynamics and helped to formulate a simple observational approach to finding frozen or quasi-frozen low lunar orbits that looks at minimizing the spread in the periselene altitude evolution.

Moon Frozen Orbit LRO↗

pnnl/JAX-CanVeg

Differentiable land surface model reimplementing an existing simulator, CANOAK, in JAX—a Google-developed Python package for high-performance machine learning research using automatic differentiation. The model's purpose is to perform hybrid land surface modeling that seamlessly couples process-based components with deep neural networks

Jiang, Peishi↗

Development of a Single-Phase, Transient, Subchannel Code, within the MOOSE Multi-Physics Computational Framework

Subchannel codes have been widely used for thermal-hydraulics analyses in nuclear reactors. This paper details the development of a novel subchannel code within the Idaho National Laboratory’s (INL) Multi-physics Object Oriented Simulation Environment (MOOSE). MOOSE is a parallel computational framework targeted at the solution of systems of coupled, nonlinear partial differential equations, that often arise in the simulation of nuclear processes. As such, it includes codes/modules able to solve the multiple linear and nonlinear physics that describe a nuclear reactor, under normal operation conditions or accidents. This includes thermal-hydraulics, fuel performance, and neutronics codes, between others. A MOOSE-based subchannel code is a new addition to the fleet of INL-developed codes, based on the MOOSE framework. In this work, we present the derivation of the subchannel equations for a single-phase fluid, we proceed with the description of the algorithm that is used to solve these equations and describe how this algorithm was implemented within MOOSE. We also present how this code can be coupled to the BISON fuel performance code. Next, we verify the friction model and the turbulent mixing model. We calibrate the turbulent modeling parameters for momentum mixing and enthalpy mixing, C T , β. We validate the code using experimental results and last demonstrate the coupling capabilities using a simple example.

42 ENGINEERING↗

Geothermal reservoir simulation

The prediction of long-term geothermal reservoir performance and the environmental impact of exploiting this resource are two important problems associated with the utilization of geothermal energy for power production. Our research effort addresses these problems through numerical simulation. Computer codes based on the solution of partial-differential equations using finite-element techniques are being prepared to simulate multiphase energy transport, energy transport in fractured porous reservoirs, well bore phenomena, and subsidence.

Mercer, J. W., Jr.↗