Search NASA⌕ Search

SEARCH · Search NASA

Results for “fast time simulation”

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 739 records · Page 41

Dynamics of a fractal set of first-order magnetic phase transitions in frustrated Lu 2 CoMnO 6

The axial next-nearest-neighbor Ising model predicts a fractal (infinite) set of phases with incommensurate wave vectors that are separated by first-order phase boundaries. This complexity results from a simple frustration condition between nearest- and next-nearest-neighbor interactions along a chain of Ising spins. Using x-ray photon correlation spectroscopy (XPCS), we investigate the surprising antiferromagnetic dynamics that emerge from such a complex phase diagram over a wide range of temperatures. In this work, we present XPCS measurements of the frustrated magnetic chain compound Lu 2 CoMnO 6 and Monte Carlo simulations. Incommensurate magnetic Bragg peaks slide towards commensurate “up-up-down-down” spin order with decreasing temperature and increasing time. Both simulation and experiment support a counterintuitive “upside-down” temperature dependence of the magnetic dynamics: at higher temperatures in the region of first-order phase boundaries, slower dynamics are observed where the speckle maintains its coherence. At the lowest temperatures, where part of the sample adopts commensurate order, the dynamics speed up and result in fast decoherence.

36 MATERIALS SCIENCE↗

Rocket-Based Combined Cycle Engine Concept Development

The development of rocket-based combined cycle (RBCC) propulsion systems is part of a 12 year effort under both company funding and contract work. The concept is a fixed geometry integrated rocket, ramjet, scramjet, which is hydrogen fueled and uses hydrogen regenerative cooling. The baseline engine structural configuration uses an integral structure that eliminates panel seals, seal purge gas, and closeout side attachments. Engine A5 is the current configuration for NASA Marshall Space Flight Center (MSFC) for the ART program. Engine A5 models the complete flight engine flowpath of inlet, isolator, airbreathing combustor, and nozzle. High-performance rocket thrusters are integrated into the engine enabling both low speed air-augmented rocket (AAR) and high speed pure rocket operation. Engine A5 was tested in GASL's new Flight Acceleration Simulation Test (FAST) facility in all four operating modes, AAR, RAM, SCRAM, and Rocket. Additionally, transition from AAR to RAM and RAM to SCRAM was also demonstrated. Measured performance demonstrated vision vehicle performance levels for Mach 3 AAR operation and ramjet operation from Mach 3 to 4. SCRAM and rocket mode performance was above predictions. For the first time, testing also demonstrated transition between operating modes.

Ratekin, G.↗

Optimal Control Modification Adaptive Law for Time-Scale Separated Systems

Recently a new optimal control modification has been introduced that can achieve robust adaptation with a large adaptive gain without incurring high-frequency oscillations as with the standard model-reference adaptive control. This modification is based on an optimal control formulation to minimize the L2 norm of the tracking error. The optimal control modification adaptive law results in a stable adaptation in the presence of a large adaptive gain. This study examines the optimal control modification adaptive law in the context of a system with a time scale separation resulting from a fast plant with a slow actuator. A singular perturbation analysis is performed to derive a modification to the adaptive law by transforming the original system into a reduced-order system in slow time. A model matching conditions in the transformed time coordinate results in an increase in the actuator command that effectively compensate for the slow actuator dynamics. Simulations demonstrate effectiveness of the method.

Nguyen, Nhan T.↗

In-Situ Strain Analysis of Potential Habitat Composites Exposed to a Simulated Long-Term Lunar Radiation Exposure

NASA is studying the effects of long-term space radiation on potential multifunctional composite materials for habitats to better determine their characteristics in the harsh space environment. Two composite materials were selected for the study and were placed in a test stand that simulated the stresses of a pressure vessel wall on the material. The samples in the test stand were exposed to radiation at either a fast dose rate or a slow dose rate, and their strain and temperature was recorded during the exposure. It was found that during a fast dose rate exposure the materials saw a decreased strain with time, or a shrinking of the materials. Given previous radiation studies of polymers, this is believed to be a result of crosslinking occurring in the matrix material. However, with a slow dose rate, the materials saw an increase in strain with time, or a stretching of the materials. This result is consistent with scission or degradation of the matrix occurring, possibly due to oxidative degradation.

Rojdev, Kristina↗

Software Aids Visualization of Computed Unsteady Flow

Unsteady Flow Analysis Toolkit (UFAT) is a computer program that synthesizes motions of time-dependent flows represented by very large sets of data generated in computational fluid dynamics simulations. Prior to the development of UFAT, it was necessary to rely on static, single-snapshot depictions of time-dependent flows generated by flow-visualization software designed for steady flows. Whereas it typically takes weeks to analyze the results of a largescale unsteady-flow simulation by use of steady-flow visualization software, the analysis time is reduced to hours when UFAT is used. UFAT can be used to generate graphical objects of flow visualization results using multi-block curvilinear grids in the format of a previously developed NASA data-visualization program, PLOT3D. These graphical objects can be rendered using FAST, another popular flow visualization software developed at NASA. Flow-visualization techniques that can be exploited by use of UFAT include time-dependent tracking of particles, detection of vortex cores, extractions of stream ribbons and surfaces, and tetrahedral decomposition for optimal particle tracking. Unique computational features of UFAT include capabilities for automatic (batch) processing, restart, memory mapping, and parallel processing. These capabilities significantly reduce analysis time and storage requirements, relative to those of prior flow-visualization software. UFAT can be executed on a variety of supercomputers.

Kao, David↗

Optimal Control Modification for Time-Scale Separated Systems

Recently a new optimal control modification has been introduced that can achieve robust adaptation with a large adaptive gain without incurring high-frequency oscillations as with the standard model-reference adaptive control. This modification is based on an optimal control formulation to minimize the L2 norm of the tracking error. The optimal control modification adaptive law results in a stable adaptation in the presence of a large adaptive gain. This study examines the optimal control modification adaptive law in the context of a system with a time scale separation resulting from a fast plant with a slow actuator. A singular perturbation analysis is performed to derive a modification to the adaptive law by transforming the original system into a reduced-order system in slow time. A model matching conditions in the transformed time coordinate results in an increase in the actuator command that effectively compensate for the slow actuator dynamics. Simulations demonstrate effectiveness of the method.

Adaptive Control↗

Controls-structures-interaction dynamics during RCS control of the Orbiter/SRMS/SSF configuration

During the assembly flights of the Space Station Freedom (SSF), the Orbiter will either dock with the SSF and retract to the final berthed position, or will grapple the SSF using the Shuttle Remote Manipulator System (SRMS) and maneuver the SRMS coupled vehicles to their final berthed position. The SRMS method is expected to take approximately one to one and a half hours to complete and require periodic attitude corrections by either the Orbiter or the SSF reaction control system (RCS) or continuous control by a control moment gyro (CMG) system with RCS desaturation as required. Free drift of the attached vehicles is not currently thought to be acceptable because the desired system attitude will quickly deteriorate due to unbalanced gravity gradient and aerodynamic torques resulting in power generation problems, thermodynamic control problems, and communications problems. This paper deals with the simulation and control of the SRMS during trunnion/latch interaction dynamics and during RCS maneuvers. The SRMS servo drive joints have highly non-linear elastic characteristics which tend to degrade sensitive control strategies. In addition the system natural frequencies are extremely low and depend on the drive joint deflections and SRMS geometric position. The lowest mean period of oscillation for the Orbiter/SRMS/SSF(MB6) system in brakes hold mode positioned near the final berthed position is approximately 120 seconds. A detailed finite element model of the SRMS has been developed and used in a newly developed SRMS systems dynamics simulation to investigate the non-linear transient response dynamics of the Orbiter/SRMS/SSF systems. The present SRMS control strategy of brakes only recommended by the Charles Draper Labs is contrasted with a robust controller developed by the authors. The robust controller uses an optimal inear quadratic regulator (LQR) to optimally place the closed-loop poles of a multivariable continuous-time system within the common region of an open sector with the sector angle plus or minus 45 degrees from the negative real axis, and the left-hand side of a parallel to the imaginary axis in the complex s-plane. This guarantees that the critical damping ratio for the desired control modes is equal to or in excess of 0.707. The matrix sign function is used for solving the Riccati equations which appear in the controller design procedure. Fast and stable algorithms have recently been developed for the computation of the matrix sign function. Simulation results are given which demonstrate the potential CSI involvement for the current SRMS control system and the proposed control system.

Schliesing, J. A.↗

Comparing the Performance of Julia on CPUs versus GPUs and Julia-MPI versus Fortran-MPI: a case study with MPAS-Ocean (Version 7.1)

Abstract. Some programming languages are easy to develop at the cost of slow execution, while others are fast at runtime but much more difficult to write. Julia is a programming language that aims to be the best of both worlds – a development and production language at the same time. To test Julia's utility in scientific high-performance computing (HPC), we built an unstructured-mesh shallow water model in Julia and compared it against an established Fortran-MPI ocean model, the Model for Prediction Across Scales–Ocean (MPAS-Ocean), as well as a Python shallow water code. Three versions of the Julia shallow water code were created: for single-core CPU, graphics processing unit (GPU), and Message Passing Interface (MPI) CPU clusters. Comparing identical simulations revealed that our first version of the Julia model was 13 times faster than Python using NumPy, where both used an unthreaded single-core CPU. Further Julia optimizations, including static typing and removing implicit memory allocations, provided an additional 10–20× speed-up of the single-core CPU Julia model. The GPU-accelerated Julia code was almost identical in terms of performance to the MPI parallelized code on 64 processes, an unexpected result for such different architectures. Parallelized Julia-MPI performance was identical to Fortran-MPI MPAS-Ocean for low processor counts and ranges from 2× faster to 2× slower for higher processor counts. Our experience is that Julia development is fast and convenient for prototyping but that Julia requires further investment and expertise to be competitive with compiled codes. We provide advice on Julia code optimization for HPC systems.

54 ENVIRONMENTAL SCIENCES↗

CCUS 2024, Interpreting the strain tensor Larry Murdoch Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring during CCUS. We have demonstrated this method by deploying strainmeters at shallow depths (30 to 40m) and then conducting injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. We have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS. Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and our objective was to evaluate opportunities for strain monitoring during characterization and monitoring for CCUS. Our approach was to deploy strainmeters at shallow depths (30 to 40m) and then conduct injection well tests in an underlying reservoir at 530m depth. The results indicate that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. In conclusion, we have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS.

Murdoch, Larry↗

Status Report on Fast Flux Test Facility Mechanistic Fuel Failure Experiment Analysis with BISON for Post Irradiation Examination Support

The renewed interest in metallic U-Zr nuclear fuel alloy has led to a drive for deeper understanding of the mechanisms driving the phenomena observed under irradiation conditions. The Department of Energy Advanced Fuel Campaign has developed infrastructure to support metallic fuel development, including Post Irradiation Examination (PIE) of legacy Fast Flux Test Facility (FFTF) Mechanistic Fuel Failure (MFF) experiments. The PIE performed on legacy FFTF MFF experiments gives insight on metallic fuel performance and can address the lack of knowledge and scarcity of reliable data identified in several studies over recent years. Unfortunately, PIE efforts can cost significant time and resources which can impede the progress of metallic U-Zr fuel development. Metallic U-Zr fuel performance modeling can be used to inform PIE efforts on regions of interest for relevant investigations and can help understand phenomena observed in PIE. This report demonstrates the current progress of FFTF MFF fuel performance simulations using the BISON fuel performance code and discusses the support provided by simulation to PIE efforts. Progress in temperature, profilometry, fission gas release, plenum pressure, and zirconium redistribution simulation results have been demonstrated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SREMI: Super-resolution electromagnetic imaging with single-channel ground-penetrating radar

High-resolution near-surface imaging has important applications in civil engineering, infrastructure inspection, military threat detection, geological characterization, and lunar and planetary exploration. Zero-offset, singlechannel ground penetrating radar (GPR) imaging is an established technique for near-surface target imaging and sensing but often suffers from low spatial resolution and imaging artifacts, especially of deep structures. In response, we formulate the GPR imaging as a dual-sparsity optimization problem, and develop a super-resolution electromagnetic imaging method based on a fast iterative shrinkage-thresholding algorithm. We develop our GPR imaging method in the framework of electromagnetic exploding-reflectors simulation theory, therefore the imaging method is computationally efficient. In this work, we demonstrate through synthetic and field data examples that our method can produce sharper, more reliable images with fewer artifacts compared with single-pass reverse-time migration GPR method, thus leading to improved near-surface interpretation and object identification.

58 GEOSCIENCES↗

Efficient mapping between void shapes and stress fields using Deep Convolutional Neural Networks with sparse data

Establishing fast and accurate structure-to-property relationships is an important component in the design and discovery of advanced materials. Physics-based simulation models like the finite element method (FEM) are often used to predict deformation, stress, and strain fields as a function of material microstructure in material and structural systems. Such models may be computationally expensive and time intensive if the underlying physics of the system is complex. This limits their application to solve inverse design problems and identify structures that maximize performance. In such scenarios, surrogate models are employed to make the forward mapping computationally efficient to evaluate. However, the high dimensionality of the input microstructure and the output field of interest often renders such surrogate models inefficient, especially when dealing with sparse data. Deep convolutional neural network (CNN) based surrogate models have shown great promise in handling such high-dimensional problems. In this paper, a single ellipsoidal void structure under a uniaxial tensile load represented by a linear elastic, high-dimensional and expensive-to-query, FEM model. We consider two deep CNN architectures, a modified convolutional autoencoder framework with a fully connected bottleneck and a UNet CNN, and compare their accuracy in predicting the von Mises stress field for any given input void shape in the FEM model. Additionally, a sensitivity analysis study is performed using the two approaches, where the variation in the prediction accuracy on unseen test data is studied through numerical experiments by varying the number of training samples from 20 to 100.

surrogate modeling; convolutional neural networks;↗

Normalizing flows for high-dimensional detector simulations

Whenever invertible generative networks are needed for LHC physics, normalizing flows show excellent performance. In this work, we investigate their performance for fast calorimeter shower simulations with increasing phase space dimension. We use fast and expressive coupling spline transformations applied to the CaloChallenge datasets. In addition to the base flow architecture we also employ a VAE to compress the dimensionality and train a generative network in the latent space. We evaluate our networks on several metrics, including high-level features, classifiers, and generation timing. Our findings demonstrate that invertible neural networks have competitive performance when compared to autoregressive flows, while being substantially faster during generation.

Ernst, Florian↗

Theoretical prediction of the impact of Auger recombination on charge collection from an ion track

A recombination mechanism that significantly reduces charge collection from very dense ion tracks in silicon devices was postulated by Zoutendyk et al. The theoretical analysis presented here concludes that Auger recombination is such a mechanism and is of marginal importance for higher density tracks produced by 270-MeV krypton, but of major importance for higher density tracks. The analysis shows that recombination loss is profoundly affected by track diffusion. As the track diffuses, the density and recombination rate decrease so fast that the linear density (number of electron-hole pairs per unit length) approaches a non-zero limiting value as t yields infinity. Furthermore, the linear density is very nearly equal to this limiting value in a few picoseconds or less. When Auger recombination accompanies charge transport processes that have much longer time scales, it can be simulated by assigning a reduced linear energy transfer to the ion.

Edmonds, Larry D.↗

Development of Comprehensive Reduced Kinetic Models for Supersonic Reacting Shear Layer Simulations

Large-scale simulations of multi-dimensional unsteady turbulent reacting flows with detailed chemistry and transport can be computationally extremely intensive even on distributed computing architectures. With the development of suitable reduced chemical kinetic models, the number of scalar variables to be integrated can be decreased, leading to a significant reduction in the computational time required for the simulation with limited loss of accuracy in the results. A general MATLAB-based automated mechanism reduction procedure is presented to reduce any complex starting mechanism (detailed or skeletal) with minimal human intervention. Based on the application of the quasi steady-state (QSS) approximation for certain chemical species and on the elimination of the fast reaction rates in the mechanism, several comprehensive reduced models, capable of handling different fuels such as C2H4, CH4 and H2, have been developed and thoroughly tested for several combustion problems (ignition, propagation and extinction) and physical conditions (reactant compositions, temperatures, and pressures). A key feature of the present reduction procedure is the explicit solution of the concentrations of the QSS species, needed for the evaluation of the elementary reaction rates. In contrast, previous approaches relied on an implicit solution due to the strong coupling between QSS species, requiring computationally expensive inner iterations. A novel algorithm, based on the definition of a QSS species coupling matrix, is presented to (i) introduce appropriate truncations to the QSS algebraic relations and (ii) identify the optimal sequence for the explicit solution of the concentration of the QSS species. With the automatic generation of the relevant source code, the resulting reduced models can be readily implemented into numerical codes.

Zambon, A. C.↗

North Pacific Decadal Variability in the GEOS-5 Atmosphere-Ocean Model

This study examines the mechanisms of the Pacific decadal oscillation (PDO) in the GEOS-5 general circulation model. The model simulates a realistic PDO pattern that is resolved as the first empirical orthogonal function (EOF) of winter sea surface temperature (SST). The simulated PDO is primarily forced by Aleutian low through Ekman transport and surface fluxes, and shows a red spectrum without any preferred periodicity. This differs from the observations, which indicate a greater role of El Nino-Southern Oscillation (ENSO) forcing, and likely reflects the too short time scale of the simulated ENSO. The geostrophic transport in response to the Aleutian low is limited to the Kuroshio-Oyashio Extension, and is unlikely the main controlling factor in this model, although it reinforces the Ekman-induced SST anomalies. The delay between the Aleutian low and the PDO is relatively short (1 year) suggesting that the fast Ekman response (rather than Rossby wave propagation) sets the SST pattern immediately following an Aleutian low fluctuation. The atmospheric feedback (response to the SST) is only about 25 of the forcing and never evolves into an Aleutian low completely, instead projecting onto the North Pacific Oscillation (NPO), a meridional dipole in sea level pressure (SLP). The lack of preferred periodicity and weak atmospheric response bothindicate a coupled oscillation is an unlikely mechanism for the PDO in this model. In agreement with recent studies, the NPO is correlated with the North Pacific Gyre Oscillation (NPGO), which is another leading EOF of the North Pacific SST. A possible connection between the PDO and the NPGO is discussed.

Pacific Decadal Variability↗

Capabilities, Calibration, and Impact of the ISS-RAD Fast Neutron Detector

In the current NASA crew radiation health risk assessment framework, estimates for the neutron contributions to crew radiation exposure largely rely on simulated data with sizeable uncertainties due to the lack of experimental measurements inside the ISS. Integrated in the ISS-RAD instrument, the ISS-RAD Fast Neutron Detector (FND) will deploy to the ISS on one of the next cargo supply missions. Together with the ISS-RAD Charged Particle Detector, the FND will perform, for the first time, routine and precise direct neutron measurements inside the ISS between 0.5 and 80 MeV. The measurements will close the NASA Medical Operations Requirement to monitor neutrons inside the ISS and impact crew radiation health risk assessments by reducing uncertainties on the neutron contribution to crew exposure, enabling more efficient mission planning. The presentation will focus on the FND detection mechanism, calibration results and expectations about the FND's interaction with the mixed radiation field inside the ISS.

Leitgab, Martin↗

Dense 3D-Reconstruction from Monocular Image Sequences for Computationally Constrained UAS

The ability to find safe landing sites over complex 3D terrain is an essential safety feature for fully autonomous small unmanned aerial systems (UAS), which requires on-board perception for 3D reconstruction and terrain analysis if the overflown terrain is unknown. This is a challenge for UAS that are limited in size, weight and computational power, such as small rotorcrafts executing autonomous missions on Earth, or in planetary applications such as the Mars Helicopter. For such a computationally constraint system, we propose a structure from motion approach that uses inputs from a single downward facing camera to produce dense point clouds of the overflown terrain in real time. In contrast to existing approaches, our method uses metric pose information from a visual-inertial odometry algorithm as camera pose priors, which allows deploying a fast pose refinement step to align camera frames such that a conventional stereo algorithm can be used for dense 3D reconstruction. We validate the performance of our approach with extensive evaluations in simulation, and demonstrate the feasibility with data from UAS flights.

Brockers, Roland↗