Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial 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 289 records · Page 16

Current-driven resistive ballooning modes in axially bounded solar flare plasmas

The most unstable current-driven resistive modes of an axially bounded coronal loop are found in computer simulations to exhibit the spatial structure of ballooning modes. The observed modes are not confined to mode rational surfaces, but instead have broad radial extent. A theory assuming ballooning mode spatial structure predicts that a minimum current should be required for linear instability, and that, when the mode is unstable, the linear growth rate scales linearly with the resistivity eta below a critical resistivity, and scales as cu root of eta for larger resistivities. Both predictions are borne out by simulation results. Both theory and simulation analyses of the mode suggest that the strong radial structure of the mode near the ends of the system is the primary contributing factor to the instability of the mode. A helical current sheet is formed in the nonlinear evolution of the mode near the edge of the current channel and is accompanied by a strong radial gradient in the current and partial current reversal.

Otani, Niels F.↗

Sensitivity of a Cloud-Resolving Model to Bulk and Explicit Bin Microphysical Schemes: Cloud Microphysics and Storm Dynamics Interactions - Part 2

Part I of this paper compares two simulations, one using a bulk and the other a detailed bin microphysical scheme, of a long-lasting, continental mesoscale convective system with leading convection and trailing stratiform region. Diagnostic studies and sensitivity tests are carried out in Part II to explain the simulated contrasts in the spatial and temporal variations by the two microphysical schemes and to understand the interactions between cloud microphysics and storm dynamics. It is found that the fixed raindrop size distribution in the bulk scheme artificially enhances rain evaporation rate and produces a stronger near surface cool pool compared with the bin simulation. In the bulk simulation, cool pool circulation dominates the near-surface environmental wind shear in contrast to the near-balance between cool pool and wind shear in the bin simulation. This is the main reason for the contrasting quasi-steady states simulated in Part I. Sensitivity tests also show that large amounts of fast-falling hail produced in the original bulk scheme not only result in a narrow trailing stratiform region but also act to further exacerbate the strong cool pool simulated in the bulk parameterization. An empirical formula for a correction factor, r(q(sub r)) = 0.11q(sub r)(exp -1.27) + 0.98, is developed to correct the overestimation of rain evaporation in the bulk model, where r is the ratio of the rain evaporation rate between the bulk and bin simulations and q(sub r)(g per kilogram) is the rain mixing ratio. This formula offers a practical fix for the simple bulk scheme in rain evaporation parameterization.

Li, Xiaowen↗

Improved methods for large-eddy simulations of turbulence

By using Fourier transforms for evaluating spatial derivatives, the accuracy of the large-eddy simulation of homogeneous isotropic turbulence is improved. Especially, the treatment of certain terms that arise in filtering the equations is considerably improved in both speed and accuracy. Use of vorticity as the principal variable is shown to be a viable and potentially useful alternative to the primitive variables. A method of deriving conservation properties of numerical schemes is given which is much simpler than previous methods and is widely applicable. The methods are applied to the computation of homogeneous isotropic turbulence, and it is found that the subgrid scale model is improved by using finite differences in place of 'exact' derivatives.

Mansour, N. N.↗

A three-dimensional laser ray-tracing methodology for radiation-hydrodynamics simulations

We report on a methodology for performing laser ray-tracing in three spatial dimensions for radiation-hydrodynamics simulation codes. Our method, which is an extension of that developed in Haines et al., Comput. Fluids 201, 104478 (2020), utilizes an automatically generated separate mesh for the laser ray-tracing from the radiation-hydrodynamics mesh. This enables the laser mesh to be tailored to minimize ray noise with significantly fewer rays than would be required when the ray-tracing is performed on the radiation-hydrodynamics mesh, primarily by allowing the use of high-aspect-ratio cells that are not suitable for hydrodynamics solvers. For a planar target, we show that our method provides a ≈ 100× reduction in computational expense to achieve a fixed level of ray noise relative to ray-tracing directly on the radiation-hydrodynamics mesh. The relatively low ray requirement also enables efficient computation of cross-beam energy transfer. Each cell in the logically cubic laser mesh is a non-convex dodecahedron with triangular sides, and numerical integration of the ray trajectories and inverse bremsstrahlung is performed by mapping each cell to the unit cube. We will describe our methodology in detail as well as its implementation in the xRAGE radiation-hydrodynamics code, discuss performance, and present the results from applying the methodology to test problems with analytic solutions for laser ray-tracing through a quadratic density gradient with an analytic solution as well as for a laser-driven heat front. In 3D radiation-hydrodynamics simulations of laser-driven experiments performed on the National Ignition Facility, laser ray-tracing with our methodology uses less than 1% of total computational time while introducing acceptably low levels of ray noise.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SPARTA: A flux adjustment methodology to interpret complex experiments

For the accurate determination of reactivity from a detector count rate, correction of spatial effects is of prime importance. This spatial correction is often provided using simulation methodologies, but this may introduce a bias if the result of the experiment is also used as input data for the simulation. Here, this work presents a flux adjustment methodology able to infer experimental reactivity and correction of spatial effects without the need for a simulation. It can process the signal from a complex experiment such as a heat balance measurement in the TREAT reactor, where control rods are continuously adjusted to maintain a constant power. In the present work, this methodology successfully computed the reactivity and the local spatial variation of the flux of a generated signal. It also proved to be robust against noise and errors on kinetic parameters and provides a credible interpretation of a heat balance experiment in TREAT. Efficiency of flux adjustment methods for complex experiment enable a better experiment interpretation less reliant on nuclear data evaluation.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The sensitivity of terrestrial carbon storage to climate change

The role of the terrestrial biosphere in controlling atmospheric CO2 levels during climate perturbations is estimated. Simulations are used to calculate the local geographical distribution of vegetation during the last glacial maximum. The known changes in sea level at this time, together with the simulated climate-driven spatial arrangement of vegetation, result in a mass transfer of carbon from the terrestrial biosphere to the atmosphere ranging from 30 Gt (corresponding to 15 ppm CO2) to -50 Gt (25 ppm). Thus, although the biosphere may have contributed to the decrease in atmospheric CO2 at 80 ppm, known to have occurred at 18 kyr, it does not seem to have been a dominant factor. For simulations run with twice the present-day CO2 levels, strong negative feedbacks appear which remove 235 Gt of carbon from the atmosphere.

Prentice, Katharine C.↗

Direct Evidence of Multispecies Hydrodynamics in Ignition-Scale Hohlraums

A targeted experiment at the National Ignition Facility (NIF) confirms the presence of multispecies hydrodynamics in inertial confinement fusion hohlraums relevant to ignition. Here, the effects are identified by filling the gold hohlraum with a deuterium-tritium (DT) gas mixture instead of helium. As the hohlraum is heated by the NIF lasers, it implodes inward, compressing and heating the DT, which leads to fusion. The resulting DT-fusion neutrons are measured in space, time, yield, angle, and energy. A distinct, peaked, triangular shape in the radial neutron emission profile provides evidence of a “leaky piston” effect caused by the interpenetration of DT into the expanding gold. This process reduces the reversibility of DT compression and decreases neutron generation on the temporal trailing edge compared to the leading edge. These results are well described by multispecies hydrodynamics simulations, which reproduce the observed spatial and temporal features, as well as the total neutron yield and angularly resolved energy spectra. In contrast, conventional simulations that use only single-species physics fail to match the experimental data. They overpredict the neutron yield, produce a flatter-than-observed spatial profile, and show excessive emission at late times compared to the experimental measurements.

Higginson, Drew P. [Lawrence Livermore National La↗

Vortex simulation of three-dimensional mixing layers

Spanwise structures of 'nominally' two-dimensional (2D), spatially growing, turbulent mixing layers are simulated numerically by a three-dimensional (3D) vortex method. Small-amplitude 3D disturbances are introduced into an otherwise 2D flow field. Results show that a large-scale spanwise variation of the flow field is produced because of amplification of initial disturbances. Pairs of counterrotating streamwise vortices are formed as a result of stretching of primary spanwise vortices. The streamwise vortices are formed at a fixed spanwise location which depends on initial disturbances. The calculated magnitude of streamwise vorticity is close to that of the spanwise vorticity. The results also suggest that the presence of background disturbances may be essential for the formation of streamwise vortices. The results are in good qualitative agreement with experiments.

Inoue, Osamu↗

An evolving MHD vortex street model for quasi-periodic solar wind fluctuations

Magnetohydrodynamic (MHD) simulation are used to provide a dynamical basis for the 'vortex street' model of the quasi-periodic meridional flow observed by Voyager 2 in the outer heliosphere. Various observations suggest the existence near the current sheet at solar minimum, of a vorticity distribution of two opposite shear layers with an antisymmetric staggered velocity pattern due to structured high-speed wind surrounding low-speed equatorial flow. It is shown that this flow pattern leads to the formation of a highly stable vortex street through the nonlinear interaction of the two shear layers. Spatial profiles of various simulated parameters (velocity, density, meridional flow angle and the location of magnetic sector boundaries) and their relative locations in the quasi-steady vortex street are generally in good agreement with the observations.

Siregar, Edouard↗

Prediction of Broadband Blade-Wake Back-Scatter Noise from a Hovering Ideally Twisted Rotor using OVERFLOW2-ANOPP2

This work illustrates the use of a hybrid RANS/LES framework for the prediction of broadband noise via integral formulations from an ideally twisted hovering rotor. Four simulations were conducted, which investigated the effect of different turbulence modeling and spatiotemporal resolutions. The aerodynamic performance and aeroacoustics predicted by these four computational cases were compared against data acquired in the Small Hover Anechoic Chamber at the NASA Langley Research Center and against previous simulations conducted using the lattice-Boltzmann solver, PowerFLOW. The simulation case with the finest spatial resolution had the most favorable aerodynamic performance and broadband noise comparison with measured data; however, the simulation case utilizing the SST-DDES turbulence model predicted the most accurate tonal noise signature. Aerodynamic flowfield visualizations identified near-wake trailing edge vortex shedding, perpendicular blade-vortex interactions, and tip vortex formation, which are known to cause broadband noise. The fluid dynamics associated with a new broadband noise source, blade-wake back-scatter, were also identified. Lastly, a broadband noise source visualization technique was utilized to confirm the presence of blade-wake back-scatter as a dominant broadband noise source.

blade-wake back-scatter↗

Impact of Spatial Sampling on Continuity of MODIS-VIIRS Land Surface Reflectance Products: A Simulation Approach

With the increasing need to construct long-term climate-quality data records to understand, monitor, and predict climate variability and change, it is vital to continue systematic satellite measurements along with the development of new technology for more quantitative and accurate observations. The Suomi National Polar-orbiting Partnership mission provides continuity in monitoring the Earths surface and its atmosphere in a similar fashion as the heritage MODIS instruments onboard the National Aeronautics and Space Administrations Terra and Aqua satellites. In this paper, we aim at quantifying the consistency of Aqua MODIS and Suomi-NPP Visible Infrared Imaging Radiometer Suite (VIIRS) Land Surface Reflectance (LSR) and NDVI products as related to their inherent spatial sampling characteristics. To avoid interferences from sources of measurement and/or processing errors other than spatial sampling, including calibration, atmospheric correction, and the effects of the bidirectional reflectance distribution function, the MODIS and VIIRSLSR products were simulated using the Landsat-8s Operational Land Imager (OLI) LSR products. The simulations were performed using the instruments point spread functions on a daily basis for various OLI scenes over a 16-day orbit cycle. It was found that the daily mean differences due to discrepancies in spatial sampling remain below 0.0015 (1) in absolute surface reflectance at subgranule scale (i.e., OLI scene size).We also found that the MODISVIIRS product intercomparisons appear to be minimally impacted when differences in the corresponding view zenith angles (VZAs) are within the range of -15deg to -35deg (VZA(sub v) - VZA(sub m)), where VIIRS and MODIS footprints resemble in size. In general, depending on the spatial heterogeneity of the OLI scene contents, per-grid-cell differences can reach up to 20.Further spatial analysis of the simulated NDVI and LSR products revealed that, depending on the user accuracy requirements for product intercomparisons, spatial aggregations may be used. It was found that if per-grid-cell differences on the order of 10(in LSR or NDVI) are tolerated, the product intercomparisons are expected to be immune from differences in spatial sampling.

VIIRS↗

Simulation of electrostatic ion instabilities in the presence of parallel currents and transverse electric fields

A spatially two-dimensional electrostatic PIC simulation code was used to study the stability of a plasma equilibrium characterized by a localized transverse dc electric field and a field-aligned drift for L is much less than Lx, where Lx is the simulation length in the x direction and L is the scale length associated with the dc electric field. It is found that the dc electric field and the field-aligned current can together play a synergistic role to enable the excitation of electrostatic waves even when the threshold values of the field aligned drift and the E x B drift are individually subcritical. The simulation results show that the growing ion waves are associated with small vortices in the linear stage, which evolve to the nonlinear stage dominated by larger vortices with lower frequencies.

Nishikawa, K.-I.↗

Simulation of electrostatic turbulence due to sheared flows parallel and transverse to the magnetic field

A spatially two-dimensional electrostatic particle simulation code is used to examine the stability of a plasma equilibrium characterized by a localized transverse dc electric field and a magnetic-field-aligned electron drift for L much less than Lx, where Lx is the simulation length in the x direction and L is the scale length associated with the dc electric field. It is found that the dc electric field and the field-aligned current can together play a synergistic role to enable the excitation of electrostatic waves even when the threshold values of the field-aligned drift and the E x B drift are individually subcritical. The simulation results indicate that a broadband turbulence is associated with such an equilibrium.

Nishikawa, K.-I.↗

Hairpin Vortices and the Final Stages of Transition

A spatially developing direct numerical simulation has been performed for flow over a flat plate that is subjected to a one-time fluid injection through an elongated slit in the wall. The flow parameters have been chosen to closely approximate the experimental conditions of Haidari, Taylow, and Smith (AIAA-89-0964). A hairpin vortex quickly develops near the upstream end of the slit, and a pair of necklace vortices form around the slow-moving injection fluid. As seen in the experiments and reported in Haidari and Smith (in review, JFM), the hairpin vortex spawns both in-line and sidelobe secondary vortices. However, no subdsidiary vortices (those formed by the inviscid deformation of a vortex-line bundle) are observed. At later times, a set of three different types of vortices are identified: hairpin vortex structures with heads that rise away from the wall horseshoe-shaped vortices that do not rise out of the boundary layer, and quasi-streamwise vortices. These structures interact with each other and with the wall layer to generate new vortices that are similar in structure to those mentioned above, although a particular parent vortex may have an offspring that more nearly resembles another member of the set. Perturbation velocity and vertical vorticity contours reveal an arrowhead shape of the highly disturbed region that is reminiscent of a turbulent spot. Spatially averaged velocity profiles in the highly disturbed area are nonlaminar, but as yet do not show typical low-of-the-wall behavior.

Singer, Bart A.↗

Recalibration of the NASA Exterior Effects Room

The Exterior Effects Room (EER) is a psychoacoustic testing facility at NASA Langley Research Center which primarily focuses on testing human response to aircraft noise. The 39-seat auditorium houses a real-time spatial audio system and flyover simulation environment. The audio server utilizes an implementation of three-dimensional vector base amplitude panning (VBAP), a perceptual spatial audio technique exploiting loudspeaker triplets to place a virtual sound source at an arbitrary spatial position. Due to the irregular room geometry and non-uniform loudspeaker setup, the audio server applies equalization filtering to compensate for spectral coloration attributed to loudspeaker installation, crossover filtering, and delay/gain offsets. These filters can incorporate measurements taken at multiple listening points, allowing a more extensive listening area than filters derived from measurements taken at a single listening point. This work endeavors to update the existing equalization filter generation process by producing a new set of filters for various human subject testing scenarios.

calibration↗

Recalibration of the NASA Exterior Effects Room

The Exterior Effects Room (EER) is a psychoacoustic testing facility at NASA Langley Research Center which primarily focuses on testing human response to aircraft noise. The 39-seat auditorium houses a real-time spatial audio system and flyover simulation environment. The audio server utilizes an implementation of three-dimensional vector base amplitude panning (VBAP), a perceptual spatial audio technique exploiting loudspeaker triplets to place a virtual sound source at an arbitrary spatial position. Due to the irregular room geometry and nonuniform loudspeaker setup, the audio server applies equalization filtering to compensate for spectral coloration attributed to loudspeaker installation, crossover filtering, and delay/gain offsets. These filters can incorporate measurements taken at multiple listening points, allowing a more extensive listening area than filters derived from measurements taken at a single listening point. This work endeavors to update the existing equalization filter generation process by producing a new set of filters for various human subject testing scenarios.

Signal Processing↗

An integrated modeling framework with open architecture for phase field simulation of multi-component alloys

An integrated modeling framework (PanPhaseField) has been developed, which enables a direct and fast coupling between CALPHAD calculations and large-scale phase field simulations for multi-component alloys. Further, it adopts an open architecture allowing for integration of user-defined phase field models in a plug-and-play manner by taking full advantage of the user-friendly graphical interface of Pandat software. The developed modeling platform becomes an enabling tool that can be used to simulate the evolution of spatially varying microstructures of industrial complex alloys for various engineering applications.

36 MATERIALS SCIENCE↗

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part I: Lagrangian and indirect Euler AMR algorithms

Many applications of physics and engineering involve wide ranges of time and spatial scales. The numerical simulation of localized small scales such as shock waves and material interfaces requires a large number of computational cells in these regions. For these applications, Lagrangian and Arbitrary-Lagrangian-Eulerian (ALE) related methods are engaging since the moving mesh feature naturally brings mesh cells on shock discontinuities and material interfaces are carefully captured. In addition, Adaptive-Mesh-Refinement (AMR) strategies aim to optimize computational resources by concentrating finer mesh cells only in areas of interest while using coarser cells elsewhere. A key but challenging AMR requirement consists in efficiently distributing the computational effort to achieve high accuracy without the prohibitive computational costs associated with uniformly fine grids. Here, in this document, the coupling of the p4est AMR library with a cell-centered Lagrangian scheme is presented with the goal to perform reliable 3D Lagrangian-AMR and indirect Euler-AMR multi-material simulations. In particular, it is shown that starting from a 3D indirect ALE code, the memory management and load balancing requirements can be delegated to an external library (here the p4est library) to unlock ALE-AMR capabilities. First, we present a strategy to transcribe the octant-based connectivity of the 3D AMR framework with that of an unstructured mesh of polygonal cells used in Lagrangian hydrodynamics. Then, we show how refinement and coarsening operations must be adapted to the particular Lagrangian framework to ensure the conservation of volume during those steps. Finally, several numerical test cases are presented that demonstrate the capabilities of the Lagrangian-AMR and indirect Euler-AMR algorithms.

3D cell-centered Lagrangian numerical scheme↗