SEARCH · Search NASA

Results for “predictive 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS2006

Applying Parallel Adaptive Methods with GeoFEST/PYRAMID to Simulate Earth Surface Crustal Dynamics

This viewgraph presentation reviews the use Adaptive Mesh Refinement (AMR) in simulating the Crustal Dynamics of Earth's Surface. AMR simultaneously improves solution quality, time to solution, and computer memory requirements when compared to generating/running on a globally fine mesh. The use of AMR in simulating the dynamics of the Earth's Surface is spurred by future proposed NASA missions, such as InSAR for Earth surface deformation and other measurements. These missions will require support for large-scale adaptive numerical methods using AMR to model observations. AMR was chosen because it has been successful in computation fluid dynamics for predictive simulation of complex flows around complex structures.

GeoFEST↗

DOE OSTI2024

Validation of SOLPS-ITER and EDGE2D-EIRENE simulations for H, D, and T JET ITER-like wall low-confinement mode plasmas

Both experiments and simulations with SOLPS-ITER and EDGE2D-EIRENE show that the onset of detachment for the low-field side (LFS) divertor – defined here as the line-averaged upstream density (〈n e 〉 edge ) at which the plasma flux to the LFS target (I LFS-plate ) starts to decrease with increasing 〈n e 〉 edge – is independent of the isotope mass. However, there are three major simulation-experiment discrepancies: (i) the absolute values of I LFS-plate and the electron density (n e ) in the LFS divertor at the onset of detachment are significantly lower in simulations, i.e., approximately a factor of 2 for I LFS-plate and a factor of 3-4 for n e ; (ii) the degree of detachment – defined here as the difference between I LFS-plate at the onset of detachment and at an 〈n e 〉 edge value close to the density limit – is smaller in simulations compared to experiments; and (iii) the experimentally observed larger degree of detachment for D and T plasmas compared to H plasmas cannot be clearly distinguished from the simulation results. There are strong indications that discrepancy (i) is to a large extent caused by neglecting Lyman-opacity effects in our simulations. The simulations predict a similar net volumetric recombination source for all isotopes due to the fact that molecule-activated recombination (MAR) compensates for the reduced electron–ion recombination (EIR) for H, whereas MAR is negligible for D and T. This similar net volumetric recombination source for all isotopes leads to an isotope-independent degree of detachment in simulations. An analysis of the Balmer-α and Balmer-$γ$ emission confirms the underestimate of MAR in simulations (especially for D and T) for the JET metallic wall, which was previously observed for devices with a carbon wall. The underestimate of MAR is an important cause for discrepancy (ii) and the fact that there is a stronger underestimate of MAR for D and T than for H explains discrepancy (iii). Extending the plasma grid to the vessel wall increases I LFS-plate and n e at the onset of detachment by 25%, and the EIR source increases by 80% in detached conditions. Hence, while the extended grid results are closer to the experimental observations, the previously described qualitative discrepancies still persist.

Balmer emission↗

DOE OSTI2024

Effects of drive pulse shape on graded metal pushered single shell capsule implosions on the National Ignition Facility

Graded metal pushered single shells (PSS) are a viable alternative to low-Z capsules (Z is the atomic number) for indirect drive inertial confinement fusion implosions due to enhanced core tamping and radiation trapping, but they can be compromised by the pusher mixing with the fuel. We compare 2-shock and 3-shock laser pulses for Be/Cr PSS capsules filled with deuterium–tritium gas fuel at 6 mg/cc density. 1D radiation-hydrodynamic simulations predict higher core compression and, hence, ∼2× higher fusion yield for the 3-shock drive than for 2-shock. Nevertheless, we observe similar core ion temperatures and fusion yields for both drives. The implosion burn duration is 25% shorter and the core volume is ∼2.5× smaller for the 3-shock drive than for 2-shock, consistent with a higher compression. 1D LASNEX mix simulations using a buoyancy-drag model matching the measured yields also agree with the observed core sizes and burn durations and suggest ∼40% and ∼70% yield degradations for 2-shock and 3-shock drives due to hydrodynamic instabilities and atomic mix at the pusher–fuel interface. At the same time, 2D HYDRA simulations show that mid-mode (2–250) instability degradations are negligible for the 2-shock implosion (9%) and significant (45%) for 3-shock. Subtracting these from the 1D mix simulations, we infer similar degradations from high-mode instabilities and atomic mix for both drives. Due to its robustness to mid-mode instabilities, future pusher–gas mix studies will use the 2-shock drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2020

Anomalous Ion Heating and Pole Erosion in the 12.5-kW Hall Effect Rocket with Magnetic Shielding (HERMeS)

Wear measurements of the inner pole cover in HERMeS have shown as much as ten times higher erosion rates than past numerical simulations predicted. Also, measurements of ion velocity distribution functions (IVDF) suggest that ion heating occurs near the chamfered regions of the acceleration channel and around the inner pole surface. In our companion paper we present theoretical evidence that the plasma conditions in this region can excite instabilities in the lower hybrid frequency. Here, we provide further evidence of their presence through comparisons between numerical simulations and measurements. After including in the 2-D axisymmetric hydrodynamic-PIC code Hall2De an idealized model of the ion temperature based on a scaling with the lower hybrid frequency, we find the agreement between the computed and measured IVDF is significantly improved. In addition, the computed erosion rates at the inner pole surface are now in close agreement with wear test measurements.

Frieman, Jason↗

DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗

NASA NTRS2012

Validation of Aquarius Measurements Using Radiative Transfer Models at L-Band

Aquarius/SAC-D was launched in June 2011 by NASA and CONAE (Argentine space agency). Aquarius includes three L-band (1.4 GHz) radiometers dedicated to measuring sea surface salinity. We report detailed comparisons of Aquarius measurements with radiative transfer model predictions. These comparisons were used as part ofthe initial assessment of Aquarius data. In particular, they were used successfully to estimate the radiometer calibration bias and stability. Further comparisons are being performed to assess the performance of models in the retrieval algorithm for correcting the effect of sources of geophysical "noise" (e.g. the galactic background, atmospheric attenuation and reflected signal from the Sun). Such corrections are critical in bringing the error in retrieved salinity down to the required 0.2 practical salinity unit (psu) on monthly global maps at 150 km by 150 km resolution. The forward models making up the Aquarius simulator have been very useful for preparatory studies in the years leading to Aquarius' launch. The simulator includes various components to compute effects ofthe following processes on the measured signal: 1) emission from Earth surfaces (ocean, land, ice), 2) atmospheric emission and absorption, 3) emission from the Sun, Moon and celestial Sky (directly through the antenna sidelobes or after reflection/scattering at the Earth surface), 4) Faraday rotation, and 5) convolution of the scene by the antenna gain patterns. Since the Aquarius radiometers tum-on in late July 2011, the simulator has been used to perform a first order validation of the data. This included checking the order of magnitude ofthe signal over ocean, land and ice surfaces, checking the relative amplitude of signal at different polarizations, and checking the variation with incidence angle. The comparisons were also used to assess calibration bias and monitor instruments calibration drift. The simulator is also being used in the salinity retrieval. For example, initial assessments of the salinity retrieved from Aquarius data showed degradation in accuracy at locations where glint from the galactic sky background was important. This was traced to an inaccurate correction for the Sky glint. We present comparisons of the simulator prediction to the Aquarius data in order to assess the performances of the models of various physical processes impacting the measurements, such as the effect of sea surface roughness, the impact of the celestial Sky and the Sun emission scattered at the rough ocean surface. We discuss what components of the simulator appear reliable and which ones need improvements. Improved knowledge on the radiative transfer models at L-band will not only lead to better salinity retrieved from Aquarius data, it will also allow be beneficial for SMOS or the upcoming SMAP mission.

Dinnat, E.↗

NASA NTRS2023

Hybrid CFD Engineering Model of Plume Induced Erosion and Crater Formation During Descent of Lunar Landers

With rapidly increased worldwide interest in landing on the moon, the issue of Plume Surface Interactions (PSI) is gaining attention. Hazards posed by lander plume induced dust and debris, as well as landing site deformation can be mitigated when better understood through predictive simulations. As simulation enabling computational power continues to increase, hybrid Computational Fluid Dynamics (CFD)/Engineering models provide the immediate ability to conduct parametric/trade studies driving design decisions for landers. Reduced order erosion engineering models apply correlations of the surface erosion rate to the plume induced surface forces with the correlations anchored to flight observations from Apollo LM landings. The MSFC propulsion fluid dynamics branch has developed such a hybrid model for predicting the plume induced viscous erosive regression of the Lunar surface beneath a landing vehicle. The (Descent Interpolated Gas Granular Erosion Model) DIGGEM was originally implemented as a post processing tool to calculate induced erosion rates through vehicle descent using CFD solutions of the vehicle plume at several fixed altitudes. This capability has since been advanced in the Loci/Chem-DIGGEM model to allow transient, moving vehicle, fully coupled viscous erosion modeling of vehicle descent PSI. This model has recently been used to make preflight predictions of the erosive regression of the ground beneath a Commercial Lunar Payload Services (CLPS) vehicle in support of measurements to be made by the Stereo Cameras for Lunar Plume Surface Studies (SCALPSS) instrument.

Plume Surface Interaction↗

NASA NTRS2023

Hybrid CFD Engineering Model of Plume Induced Erosion and Crater Formation During Descent of Lunar Landers

Plume Surface Interaction↗

NASA NTRS2001

Osmotic forces and gap junctions in spreading depression: a computational model

In a computational model of spreading depression (SD), ionic movement through a neuronal syncytium of cells connected by gap junctions is described electrodiffusively. Simulations predict that SD will not occur unless cells are allowed to expand in response to osmotic pressure gradients and K+ is allowed to move through gap junctions. SD waves of [K+]out approximately 25 to approximately 60 mM moving at approximately 2 to approximately 18 mm/min are predicted over the range of parametric values reported in gray matter, with extracellular space decreasing up to approximately 50%. Predicted waveform shape is qualitatively similar to laboratory reports. The delayed-rectifier, NMDA, BK, and Na+ currents are predicted to facilitate SD, while SK and A-type K+ currents and glial activity impede SD. These predictions are consonant with recent findings that gap junction poisons block SD and support the theories that cytosolic diffusion via gap junctions and osmotic forces are important mechanisms underlying SD.

Spreading Cortical Depression/physiology↗

NASA NTRS1992

Spatial orientation of the vestibular system

1. A simplified three-dimensional state space model of visual vestibular interaction was formulated. Matrix and dynamical system operators representing coupling from the semicircular canals and the visual system to the velocity storage integrator were incorporated into the model. 2. It was postulated that the system matrix for a tilted position was a composition of two linear transformations of the system matrix for the upright position. One transformation modifies the eigenvalues of the system matrix while another rotates the pitch and roll eigenvectors with the head, while maintaining the yaw axis eigenvector approximately spatially invariant. Using this representation, the response characteristics of the pitch, roll, and yaw eye velocity were obtained in terms of the eigenvalues and associated eigenvectors. 3. Using OKAN data obtained from monkeys and comparing to the model predictions, the eigenvalues and eigenvectors of the system matrix were identified as a function of tilt to the side or of tilt to the prone positions, using a modification of the Marquardt algorithm. The yaw eigenvector for right-side-down tilt and for downward pitch cross-coupling was approximately 30 degrees from the spatial vertical. For the prone position, the eigenvector was computed to be approximately 20 degrees relative to the spatial vertical. For both side-down and prone positions, oblique OKN induced along eigenvector directions generated OKAN which decayed to zero along a straight line with approximately a single time constant. This was verified by a spectral analysis of the residual sequence about the straight line fit to the decaying data. The residual sequence was associated with a narrow autocorrelation function and a wide power spectrum. 4. Parameters found using the Marquardt algorithm were incorporated into the model. Diagonal matrices in a head coordinate frame were introduced to represent the direct pathway and the coupling of the visual system to the integrator. Model simulations predicted the behavior of yaw and pitch OKN and OKAN when the animal was upright, as well as the cross-coupling in the tilted position. The trajectories in velocity space were also accurately simulated. 5. There were similarities between the monkey eigenvectors and human perception of the spatial vertical. For side-down tilts and downward eye velocity cross-coupling, there was only an Aubert (A) effect. For upward eye velocity cross-coupling there were both Muller (E) and Aubert (A) effects. The mean of the eigenvectors for upward and downward eye velocities overlay human 1 x g perceptual data.(ABSTRACT TRUNCATED AT 400 WORDS).

NASA Discipline Neuroscience↗

NASA NTRS1989

Optical fiber sensors and signal processing for intelligent structure monitoring

Few mode optical fibers have been shown to produce predictable interference patterns when placed under strain. The use is described of a modal domain sensor in a vibration control experiment. An optical fiber is bonded along the length of a flexible beam. Output from the modal domain sensor is used to suppress vibrations induced in the beam. A distributed effect model for the modal domain sensor is developed. This model is combined with the beam and actuator dynamics to produce a system suitable for control design. Computer simulations predict open and closed loop dynamic responses. An experimental apparatus is described and experimental results are presented.

Thomas, Daniel↗

DOE OSTI2024

First-principles study of the magneto-Raman effect in van der Waals layered magnets

Magneto-Raman spectroscopy has been used to study spin-phonon coupling in two-dimensional (2D) magnets. Raman spectra of CrI 3 show a strong dependence on the magnetic order within a layer and between the layers. Here we carry out the first systematic theoretical investigation of the magneto-Raman effect in 2D magnets by performing density functional theory calculations and developing a generalized polarizability model. Our first-principles simulations well reproduce experimental Raman spectra of CrI 3 with different magnetic states. The model reveals how the change of spin orientation in each layer is coupled to the layer’s vibration to induce or eliminate the spin-dependent anti-symmetric off-diagonal terms in the Raman tensor for altering the selection rules. We also uncover that the correlation between phonon modes and magnetic orders is a universal phenomenon, which should exist in other phonon modes and 2D magnets. Our predictive simulations and modeling are expected to guide the research in 2D magnets.

36 MATERIALS SCIENCE↗

NASA NTRS1991

Large-eddy simulation of compressible turbulence

The increase in the range of length scales with increasing Reynolds number limits the direct simulation of turbulent flows to relatively simple geometries and low Reynolds numbers. However, since most flows of engineering interest occur at much higher Reynolds number than is currently within the capabilities of full simulation, prediction of these flow fields can only be obtained by solving some suitably-averaged set of governing equations. In the traditional Reynolds-averaged approach, the Navier-Stokes equations are averaged over time. This in turn yields correlations between various turbulence fluctuations. It is these terms, e.g. the Reynolds stresses, for which a turbulence model must be derived. Turbulence modeling of incompressible flows has received a great amount of attention in the literature. An area of research that has received comparatively less attention is the modeling of compressible turbulent flows. An approach to simulating compressible turbulence at high Reynolds numbers is through the use of Large-Eddy Simulation (LES). In LES the dependent variables are decomposed into a large-scale (resolved) component and a sub-grid scale component. It is the small-scale components of the velocity field which are presumably more homogeneous than the large scales and, therefore, more easily modeled. Thus, it seems plausible that simpler models, which should be more universal in character than those employed in second-order closure schemes, may be developed for LES of compressible turbulence. The objective of the present research, therefore, is to explore models for the Large-Eddy Simulation of compressible turbulent flows. Given the recent successes of Zeman in second order closure modeling of compressible turbulence, model development was guided by principals employed in second-order closures.

Squires, Kyle D.↗

NASA NTRS2002

The Relationship Between Fidelity and Learning in Aviation Training and Assessment

Flight simulators can be designed to train pilots or assess their flight performance. Low-Fidelity simulators maximize the initial learning rate of novice pilots and minimize initial costs; whereas, expensive, high-fidelity simulators predict the realworld in-flight performance of expert pilots (Fink & Shriver, 1978 Hays & Singer 1989; Kinkade & Wheaton. 1972). Although intuitively appealing and intellectually convenient to generalize concepts of learning and assessment, what holds true for the role of fidelity in assessment may not always hold true for learning, and vice versa. To bring clarity to this issue, the author distinguishes the role of fidelity in learning from its role in assessment as a function of skill level by applying the hypothesis of Alessi (1988) and reviewing the Laughery, Ditzian, and Houtman (1982) study on simulator validity. Alessi hypothesized that there is it point beyond which one additional unit of flight-simulator fidelity results in a diminished rate of learning. The author of this current paper also suggests the existence of an optimal point beyond which one additional unit of flight-simulator fidelity results in a diminished rate of practical assessment of nonexpert pilot performance.

Noble, Cliff↗

NASA NTRS2022

Acoustic Modeling of Novel Over-The-Rotor Acoustic Liner Concepts Using COMSOL Multiphysics Simulation Software

This paper presents simulation predictions of novel over-the-rotor acoustic liner configurations compared with experimental data. This is a follow on to the investigation on the effect of grazing flow on grooved over-the-rotor casing treatments. The purpose of this modeling investigation stems from the determination that traditional liner modeling techniques are inadequate for nontraditional internal liner cavity geometries. A series of these expansion chamber type liner design concepts were previously developed by Glenn Research Center and tested using the Langley Research Center Normal Incidence Tube. The measured normalized impedance and absorption coefficient is compared with that obtained from the COMSOL Multiphysics® simulation. Additionally, using a parametric case study, another series of simulations were conducted to determine what additional benefits, if any, this genre of liner has to offer the acoustics community. A total of 53 permutations of this concept was investigated. The results indicate that with longer internal fins, the resonance is sifted to lower frequencies.

Acoustics↗

DOE OSTI2025

Three-dimensional reconstruction of inertial confinement fusion hot-spot plasma from x-ray and nuclear diagnostics on OMEGA

Multidimensional effects degrade the neutron yield and the compressed areal density of laser-direct-drive inertial confinement fusion implosions of layered deuterium–tritium cryogenic targets on the OMEGA Laser System with respect to 1D radiation-hydrodynamic simulation predictions. A comprehensive physics-informed 3D reconstruction effort is under way to infer hot-spot and shell conditions at stagnation from four x-ray and seven neutron detectors distributed around the OMEGA target chamber. Neutron diagnostics, providing measurements of the neutron yield, hot-spot flow velocity, and apparent ion-temperature distribution, are used to infer the mode-1 perturbation at stagnation. The x-ray imagers record the shape of the hot-spot plasma to diagnose mode-1 and mode-2 perturbations. A deep-learning convolutional neural network trained on an extensive set of 3D radiation-hydrodynamic simulations is used to interpret the x-ray and nuclear measurements to infer the 3D profiles of the hot-spot plasma conditions and the amount of laser energy coupled to the hot-spot plasma. A 3D simulation database shows that larger mode-1 asymmetries are correlated with higher hot-spot flow velocities and reduced laser-energy coupling and neutron yield. Three-dimensional hot-spot reconstructions from x-ray measurements indicate that higher amounts of residual kinetic energy are correlated with higher measured hot-spot flow velocities, consistent with 3D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Effect of Illumination Area on the Ultrafast Temporal Response of MSM GaN Photodiodes

We investigate the influence of spatial illumination profiles on the temporal response of metal–semiconductor–metal GaN photodiodes. Using both simulation and experimental measurements, we compare the response curves under two scenarios: illumination confined to the active area between electrodes and extended illumination beyond the contacts. The results show that limiting the beam illumination to the active region significantly sharpens the response, reducing the long decay tails associated with slow carrier drift from peripheral regions. As a result, the experimental data closely match simulation predictions, confirming that illumination geometry plays a critical role in optimizing photodiode performance for ultrafast detection applications.

Carrier drift time↗

NASA NTRS2019

TPSAS-NF1676L-32920-DND

In June 2019, a full-scale crash test of a Fokker F28 Fellowship aircraft will be conducted at NASA Langley Research Center?s Landing and Impact Research (LandIR) Facility. The F28 is a high-performance twin-turbo fan narrow-body aircraft with seating in a 3+2 configuration. The MK4000 variant, used in this test, is capable of carrying up to 85 passengers on medium range routes. The F28 was first type certified by the Federal Aviation Administration (FAA) in 1969 and now the majority of the F28 fleet has retired from service in the United States. In 2016, the FAA and NASA Langley Research Center (LaRC) signed an interagency agreement for conducting a research program to obtain test data that will support the development of airframe level crash requirements for transport category airplanes [1]. The objectives of the full-scale crash test can be divided into six categories: (1) To compare and contrast responses in identical aircraft undergoing vertical only to combined vertical and horizontal loading conditions, (2) To examine the effects of horizontal loading on aircraft structure during a crash event, (3) To generate data for the use in calibration of computer modelling efforts, (4) To generate data from onboard Anthropomorphic Test Devices (ATDs) for the evaluation of injury, (5) To obtain data from experimental seats, and (6) To obtain data from new and novel ATD designs including Warrior Injury Assessment MANikin (WIAMan) [2], Test device for Human Occupant Response (THOR), and other newly developed child ATDs. The focus of this presentation will be to document finite element model development of the full-scale F28 aircraft and to present preliminary test-analysis predictions. NASA obtained the full-scale F28 aircraft, plus three fuselage sections and two sets of wings with funding through the NASA Aviation Safety Program in 1998. In addition to the hardware, NASA purchased a full NASTRAN model of the airframe that had been developed by the Dutch manufacturer, Fokker. Beginning in 2016, the NASTRAN model was converted to LS-DYNA? [4, 5] format and modified by combining parts, adding missing parts of internal structure, and including ballast for loading weights, etc. It contains: 89,223 nodes; 24,065 beam elements; 55,404 shell elements; 29,044 solid elements; 740 parts; 80 material definitions; and, 46 Constrained Nodal Rigid Bodies (CNRBs). The aircraft model will be executed in LSDYNA to simulate the test article impact onto a 2-ft.-high bed of soil under combined velocity conditions of 70-ft/s forward and 30-ft/s vertical velocity. Pre-test simulation predictions will be generated and correlated with test data.

Karen E Jackson↗