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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Predicted performance of a tangential viewing hard x-ray camera for the DIII-D high field side lower hybrid current drive experiment

High field side launch of lower hybrid current drive (LHCD) has improved accessibility and penetration over low field side launch on DIII-D. Simulations predict single pass absorption under a wide range of plasma conditions. Hard x-ray (HXR) measurement of LHCD generated fast electron bremsstrahlung (50–250 keV) will validate wave propagation and absorption. Emissivity profiles are recovered from one-dimensional inversion of HXR brightness to determine LH damping location, fast electron slowing down time, and some indication of the fast electron energy. The camera will be implemented by populating 32 tangential sightlines of the existing Gamma Ray Imager with Kromek SPEARTM Cadmium Zinc Telluride (CZT) detectors sensitive to 10–1000 keV photons with 10 keV energy resolution. Expected count rates allow for <0.5 ms time resolution. Pulses are processed using 50 ns shaping time Cremat CR-200 Gaussian shaping modules and are digitized by 25 MHz D-TACQ ACQ216 digitizers. The performance of the HXR camera is evaluated by comparing predicted fast electron density profiles and inverted synthetic brightnesses obtained from the ray-tracing/Fokker–Planck codes GENRAY/CQL3D. Inversions closely matched predicted fast electron profiles for a range of experimental parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Perturbed Ion Temperature and Toroidal Flow Profile Measurements in Rotating Neoclassical Tearing Mode Magnetic Islands

The perturbed ion temperature and toroidal flow were measured in rotating Neoclassical Tearing Modes (NTM) in a tokamak for the first time. These toroidally and radially resolved profiles were obtained by impurity ion spectroscopy in a 2,1 NTM in DIII-D. In agreement with drift-kinetic simulations, the electron temperature profile is flat, while the ion temperature gradient is restored across the magnetic island O-point in the presence of fast ions; the perturbed flow has minima in the O-points and maxima at the X-points. Furthermore, these measurements provide the first confirmation of the theoretically expected ion temperature and flow response to a magnetic island needed to predict the NTM on set threshold scaling for ITER and other future tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

A fast viscous correction method for unsteady transonic flow about airfoils

This study modified the inviscid, small-disturbances, unsteady, transonic code, LTRAN2, so that the viscous effect could be considered for the purpose of improving the code's accuracy without substantially increasing the computational time. By using a conventional integral boundary layer method to determine the displacement thickness and an empirical model of viscous wedge to simulate the suddenly thickened boundary layer behind a shock, it resulted in better agreement with experimental data for both the shock position and the pressure coefficient. Aerodynamic performances were calculated for two supercritical airfoils. When the unsteady procedure was used to obtain steady state results as the initial condition, the unsteady viscous solution for the RAE 2822 airfoil required 100.5% of the computational time of the inviscid solution. When an iterative scheme was used to obtain the steady state results as the initial condition, the unsteady viscous solution for the NLR 7301 airfoil required 88% of the computational time of the inviscid solution.

Lee, S. C.↗
DOE OSTI2020

EMP-Resilient Electric Grid Transformer Analysis

The response of a high - voltage (HV) transformer to fast rise time voltages, such as that from a electromagnetic pulse (EMP) can result in interruption of power distribution and possibly system failure. To help identify these potential occurrences, it is necessary to develop a transformer model that not only captures the input/output response of the transformer but also the internal behavior. The model constructed should cover the frequency band of interest while capturing the internal physical and electrical characteristics. This broad-band, high-fidelity model would enable the prediction of unwanted effects through simulation. A proposed modeling scheme for a HV transformer is described in Part 1 of this report. Part 2 of this report details assessments of internal voltage and electrical field holdoff testing of transformer insulation dielectric breakdown, including comparison of low frequency (DC/60 Hz) holdoff to rise times characteristic of lightning (1 s) and EMP E1 transients (10 - 30 ns). This initial project is a path toward establishing electrical grid transformer failure criteria for EMP voltage transients. We developed modeling methods and measured breakdown electrical field statistical distributions for direct current, 60 Hz, lightning and EMP characteristic voltage rise times. Methods of nanosecond-scale capacitive discharge unit high voltage source development, suggestions for derating of 60 Hz insulation maximum electrical fields for EMP nanosecond pulse voltage withstand rating, and potential methods for increasing transformer resilience to such fast rise time pulses are described.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRSDate not supplied

A Gaussian Process Enhancement to Linear Parameter Varying Models

Simulation and analysis for modern engineering systems now routinely requires the merging of multiple disciplines, physical-domains, time-scales, and data sets — all at ever increasing levels. These capabilities are especially needed in the domain of Advanced Air Mobility, where rapidly emerging vehicle designs are significantly more complex, while having to be both cost-effective and safe. To meet these engineering challenges, machine learning methods are an attractive option for merging models and data across multiple areas while providing uncertainty quantification and maintaining computational efficiency. This paper examines the use of Gaussian process machine learning to generalize and enhance the commonly used class of quasi-Linear Parameter Varying models for fast full-envelope simulation while also supporting control system design and analysis with model uncertainty. Gaussian process machine learning is selected because it: can fuse multiple data sets, enables an easy trade-off between data fitting and smoothing, provides model uncertainty quantification, scales well with increasing complexity, and does not generally require starting from a large training data set. To demonstrate the benefits of the approach, a robust stability analysis with Gaussian process uncertainty is shown for a NASA reference design of an electric quad-rotor air-taxi concept vehicle with motor parameter uncertainty.

Gaussian Process↗
NASA NTRS1992

Relative timing of substorm features in MHD simulations

An investigation of the temporal sequence of substorm phenomena based on three dimensional MHD (magnetohydrodynamic) simulations of magnetic reconnection and plasmoid formation is presented. The investigation utilizes a spatially localized resistivity model which leads to a significantly faster evolution than found in previous investigations. The analysis of the results concentrates on substorm features that have received considerable attention in the past. The formation of magnetic neutral lines, the occurrence of fast flows directed both earthward and tailward, and the magnetic field changes leading to the formation of the substorm current wedge, and to the depolarization of the magnetic field earthward of the reconnection region and its dependence on the spatial distribution of resistivity, are discussed. These phenomena are seen as an integral part of the nonlinear evolution of the three dimensional tearing instability.

Hesse, Michael↗
DOE OSTI2022

Fast flavor instability in hypermassive neutron star disk outflows

Here, we examine the effect of neutrino flavor transformation by the fast flavor instability (FFI) on long-term mass ejection from accretion disks formed after neutron star mergers. Neutrino emission and absorption in the disk set the composition of the disk ejecta, which subsequently undergoes r-process nucleosynthesis upon expansion and cooling. Here we perform 28 time-dependent, axisymmetric, viscous-hydrodynamic simulations of accretion disks around hypermassive neutron stars (HMNSs) of variable lifetime, using a 3-species neutrino leakage scheme for emission and an annular-lightbulb scheme for absorption. We include neutrino flavor transformation due the FFI in a parametric way, by modifying the absorbed neutrino fluxes and temperatures, allowing for flavor mixing at various levels of flavor equilibration, and also in a way that aims to respect the lepton-number preserving symmetry of the neutrino self-interaction Hamiltonian. We find that for a promptly-formed black hole (BH), the FFI lowers the average electron fraction of the disk outflow due to a decrease in neutrino absorption, driven primarily by a drop in electron neutrino/antineutrino flux upon flavor mixing. For a long-lived HMNS, the disk emits more heavy lepton neutrinos and reabsorbs more electron neutrinos than for a BH, with a smaller drop in flux compensated by a higher neutrino temperature upon flavor mixing. The resulting outflow has a broader electron fraction distribution, a more proton-rich peak, and undergoes stronger radiative driving. Disks with intermediate HMNS lifetimes show results that fall in between these two limits. In most cases, the impact of the FFI on the outflow is moderate, with changes in mass ejection, average velocity, and average electron fraction of order ~10%, and changes in the lanthanide/actinide mass fraction of up to a factor ~2.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

Fine pointing of Laser Beams by using Laser Arrays for Applications to CubeSats

A method of fine pointing of laser beams by using laser arrays has been developed. The telescope system combines a lens system and a VCSEL/Photodetector Array. It does not use moving parts. In computer simulations, it is applied to CubeSats that use body pointing. Body pointing was used by the Aerospace Corporation for CubeSats in LEO in NASA’s Optical Communications and Sensors Demonstration (OCSD) program. Computer simulations of this fine pointing capability have been applied previously to CubeSats in the OCSD program. In this paper, computer simulations of laser pointing using this telescope design are applied to CubeSats in LLO, at 100 km. These CubeSats could form part of the LunaNet, the lunar communications and navigation network, part of the NASA ARTEMIS Program. With more accurate pointing, a laser beam with smaller divergence can be used. For the case of the AeroCube-7B vehicle that was used in the OCSD program, computer simulations will show, for example, that the divergence of the output beam can be reduced from approximately 0.06° FWHM to 0.014°. For the proposed electro-optical system, reaction times to pointing changes and vibrations are on a nanosecond time scale, much faster than those for fine pointing systems that use moving parts such as fast steering mirrors, including MEMS, or that use quad-cell photodetectors to improve the body pointing of the CubeSat. Other possible applications are to Optical Multiple Access (OMA)for simultaneously communicating with ground stations at different locations and to Wavelength-Division Multiplexing (WDM)for increasing data rate transmission.

SMD↗
NASA NTRS2021

Fine Pointing of Laser Beams by Using Laser Arrays for Applications to CubeSats

A method of fine pointing of laser beams by using laser arrays has been developed. The telescope system combines a lens system and a VCSEL/Photodetector Array. It does not use moving parts. In computer simulations, it is applied to CubeSats that use body pointing. Body pointing was used by the Aerospace Corporation for CubeSats in LEO in NASA’s Optical Communications and Sensors Demonstration (OCSD) program. Computer simulations of this fine pointing capability have been applied previously to CubeSats in the OCSD program. In this paper, computer simulations of laser pointing using this telescope design are applied to CubeSats in LLO, at 100 km. These CubeSats could form part of the LunaNet, the lunar communications and navigation network, part of the NASA ARTEMIS Program. With more accurate pointing, a laser beam with smaller divergence can be used. For the case of the AeroCube-7B vehicle that was used in the OCSD program, computer simulations will show, for example, that the divergence of the output beam can be reduced from approximately 0.06° FWHM to 0.014°. For the proposed electro-optical system, reaction times to pointing changes and vibrations are on a nanosecond time scale, much faster than those for fine pointing systems that use moving parts such as fast steering mirrors, including MEMS, or that use quad-cell photodetectors to improve the body pointing of the CubeSat. Other possible applications are to Optical Multiple Access (OMA)for simultaneously communicating with ground stations at different locations and to Wavelength-Division Multiplexing (WDM)for increasing data rate transmission.

HEC↗
DOE OSTI2025

Micropolar Elastoplasticity Using a Fast Fourier Transform‐Based Solver

ABSTRACT This work presents a micromechanical spectral formulation for obtaining the full‐field and homogenized response of elastoplastic micropolar composites. A closed‐form radial‐return mapping is derived from thermodynamics‐based micropolar elastoplastic constitutive equations to determine the increment of plastic strain necessary to return the generalized stress state to the yield surface, and the algorithm implementation is verified using the method of numerically manufactured solutions. Then, size‐dependent material response and micro‐plasticity are shown as features that may be efficiently simulated in this micropolar elastoplastic framework. The computational efficiency of the formulation enables the generation of large datasets in reasonable computing times.

42 ENGINEERING↗
DOE OSTI2020

High mass and halo resolution from fast low resolution simulations

Generating mocks for future sky surveys requires large volumes and high resolutions, which is computationally expensive even for fast simulations. Here we try to develop numerical schemes to calibrate various halo and matter statistics in fast low resolution simulations compared to high resolution N-body and hydrodynamic simulations. For the halos, we improve the initial condition resolution and develop a halo finder "relaxed-FoF", where we allow different linking lengths for different halo mass and velocity dispersions. We show that our relaxed-FoF halo finder improves the common statistics, such as halo bias, halo mass function, halo auto power spectrum, cross correlation coefficient with the reference halo catalog, and halo-matter cross power spectrum. We also calibrate small-scale velocities of small halos to improve the power spectrum in redshift space. For the matter statistics, we incorporate the potential gradient descent (PGD) method into fast simulations to improve the matter distribution at nonlinear scales. By building a lightcone output, we show that the PGD method significantly improves the weak lensing convergence tomographic power spectrum. With these improvements FastPM is comparable to the high resolution full N-body simulation of the same mass resolution, with two orders of magnitude fewer time steps. These techniques can be used to improve the halo and matter statistics of FastPM simulations for mock catalogs of future surveys such as DESI and LSST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Development and Applications of an eReaxFF Force Field for Graphitic Anodes of Lithium-Ion Batteries

Graphene is one of the most promising materials for lithium-ion battery anodes due to its superior electronic conductivity, high surface area for lithium intercalation, fast ionic diffusivity and enhanced specific capacity. A reliable description of many battery processes requires an explicit description of electrochemical interactions involving electrons. A detailed atomistic modeling of electronic conduction and non-zero voltage simulations of graphitic materials require the inclusion of an explicit electronic degree of freedom. To enable large length- and time-scale simulations of electron conduction in graphitic anodes, we developed an eReaxFF force field concept describing graphitic materials with an explicit electron. The newly developed force field, verified against quantum chemistry-based data describing, amongst others, electron affinities and equation of states, reproduces the qualitative behavior of electron conductivity in pristine and imperfect graphitic materials at different applied temperatures and voltages. In addition, excess electron localization near a defect site estimated from eReaxFF simulations agree quite well with the corresponding density functional theory calculations. Here, our eReaxFF simulations show the initiation of lithium-metal-plating driven by electron transfer from the graphene surface to the exposed lithium ions demonstrating the method’s potential for studying lithium-graphene interactions with explicit electrons and explain many unresolved electrode and electrode-electrolyte interface processes.

25 ENERGY STORAGE↗
NASA NTRS1983

Experiments using least square lattice filters for the identification of structural dynamics

An approach for identifying the dynamics of large space structures is applied to a free-free beam. In this approach the system's order is determined on-line, along with mode shapes, using recursive lattice filters which provide a least square estimate of the measurement data. The mode shapes determined are orthonormal in the space of the measurements and, hence, are not the natural modes of the structure. To determine the natural modes of the structure, a method based on the fast Fourier transform is used on the outputs of the lattice filter. These natural modes are used to obtain the modal amplitude time series which provides the input data for an output error parameter identification scheme that identifies the ARMA parameters of the difference equation model of the modes. the approach is applied to both simulated and experimental data.

Sundararajan, N.↗
NASA NTRS1992

Inversion of surface parameters using fast learning neural networks

A neural network approach to the inversion of surface scattering parameters is presented. Simulated data sets based on a surface scattering model are used so that the data may be viewed as taken from a completely known randomly rough surface. The fast learning (FL) neural network and a multilayer perceptron (MLP) trained with backpropagation learning (BP network) are tested on the simulated backscattering data. The RMS error of training the FL network is found to be less than one half the error of the BP network while requiring one to two orders of magnitude less CPU time. When applied to inversion of parameters from a statistically rough surface, the FL method is successful at recovering the surface permittivity, the surface correlation length, and the RMS surface height in less time and with less error than the BP network. Further applications of the FL neural network to the inversion of parameters from backscatter measurements of an inhomogeneous layer above a half space are shown.

Dawson, M. S.↗
NASA NTRS1992

Sea ice classification using fast learning neural networks

A first learning neural network approach to the classification of sea ice is presented. The fast learning (FL) neural network and a multilayer perceptron (MLP) trained with backpropagation learning (BP network) were tested on simulated data sets based on the known dominant scattering characteristics of the target class. Four classes were used in the data simulation: open water, thick lossy saline ice, thin saline ice, and multiyear ice. The BP network was unable to consistently converge to less than 25 percent error while the FL method yielded an average error of approximately 1 percent on the first iteration of training. The fast learning method presented can significantly reduce the CPU time necessary to train a neural network as well as consistently yield higher classification accuracy than BP networks.

Dawson, M. S.↗
DOE OSTI2023

A fast Monte Carlo cell-by-cell simulation for radiobiological effects in targeted radionuclide therapy using pre-calculated single-particle track standard DNA damage data

Introduction: We developed a new method that drastically speeds up radiobiological Monte Carlo radiation-track-structure (MC-RTS) calculations on a cell-by-cell basis. Methods: The technique is based on random sampling and superposition of single-particle track (SPT) standard DNA damage (SDD) files from a “pre-calculated” data library, constructed using the RTS code TOPAS-nBio, with “time stamps” manually added to incorporate dose-rate effects. This time-stamped SDD file can then be input into MEDRAS, a mechanistic kinetic model that calculates various radiation-induced biological endpoints, such as DNA double-strand breaks (DSBs), misrepairs and chromosomal aberrations, and cell death. As a benchmark validation of the approach, we calculated the predicted energy-dependent DSB yield and the ratio of direct-to-total DNA damage, both of which agreed with published in vitro experimental data. We subsequently applied the method to perform a superfast cell-by-cell simulation of an experimental in vitro system consisting of neuroendocrine tumor cells uniformly incubated with 177 Lu. Results and discussion: The results for residual DSBs, both at 24 and 48 h post-irradiation, are in line with the published literature values. Our work serves as a proof-of-concept demonstration of the feasibility of a cost-effective “in silico clonogenic cell survival assay” for the computational design and development of radiopharmaceuticals and novel radiotherapy treatments more generally.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2020

Efficient Treatment of Large Active Spaces through Multi-GPU Parallel Implementation of Direct Configuration Interaction

In this study, we have extended our graphical processing unit (GPU)-accelerated direct configuration interaction program to multiple devices, reducing iteration times for configuration spaces of 165 million determinants to only 3 s using NVIDIA P100 GPUs. Similar improvements in the one- and two-particle reduced density matrix formation allow for fast analytical energy gradients and electronic properties. Our parallel algorithm enables the calculation of arbitrarily large configuration spaces (limited only by available system memory), with iteration times of 13 min for an active space of 18 electrons in 18 orbitals (2.4 billion determinants) using six consumer grade NVIDIA 1080Ti GPUs. These advances enable routine molecular dynamics simulations, geometry optimizations, and absorption spectrum calculations for molecules with large configuration spaces, a task that has heretofore required massive computational effort. In this work, we demonstrate the utility of our program by generating the absorption spectrum for diphenyl acetylene at the floating occupation molecular orbital complete active space configuration interaction level of theory. Lastly, several active spaces were investigated to assess the dependence of spectral features on orbital space dimension.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗