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

Constraints on Decreases in Eta Carinae's Mass-loss from 3D Hydrodynamic Simulations of Its Binary Colliding Winds

Recent work suggests that the mass-loss rate of the primary star Eta-A in the massive colliding wind binary Eta Carinae dropped by a factor of 2-3 between 1999 and 2010. We present result from large- (+/- 1545 au) and small- (+/- 155 au) domain, 3D smoothed particle hydrodynamics (SPH) simulations of Eta Car's colliding winds for three Eta-A mass-loss rates ( (dot-M(sub Eta-A) = 2.4, 4.8 and 8.5 × 10(exp −4) M(solar)/ yr), investigating the effects on the dynamics of the binary wind-wind collision (WWC). These simulations include orbital motion, optically thin radiative cooling and radiative forces. We find that dot-M Eta-A greatly affects the time-dependent hydrodynamics at all spatial scales investigated. The simulations also show that the post-shock wind of the companion star Eta-B switches from the adiabatic to the radiative-cooling regime during periastron passage (Phi approx.= 0.985-1.02). This switchover starts later and ends earlier the lower the value of dot-M Eta-A and is caused by the encroachment of the wind of Eta-A into the acceleration zone of Eta-B's wind, plus radiative inhibition of Eta-B's wind by Eta-A. The SPH simulations together with 1D radiative transfer models of Eta-A's spectra reveal that a factor of 2 or more drop in dot-M EtaA should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in dot- M Eta-A was likely by a factor of approx. < 2 and occurred after 2004. We speculate that most of the recent observed changes in Eta Car are due to a small increase in the WWC opening angle that produces significant effects because our line of sight to the system lies close to the dense walls of the WWC zone. A modest decrease in dot-M Eta-A may be responsible, but changes in the wind/stellar parameter of Eta-B, while less likely, cannot yet be fully ruled out. We suggest observations during Eta-Car's next periastron in 2014 to further test for decreases in dot-M Eta-A. If dot-M Eta-A is declining and continues to do so, the 2014 X-ray minimum should be even shorter than that of 2009.

Eta Carinae↗

The dimension of attractors underlying periodic turbulent Poiseuille flow

A lower bound on the Liapunov dimenison, D-lambda, of the attractor underlying turbulent, periodic Poiseuille flow at a pressure-gradient Reynolds number of 3200 is calculated, on the basis of a coarse-grained (16x33x8) numerical solution, to be approximately 352. Comparison of Liapunov exponent spectra from this and a higher-resolution (16x33x16) simulation on the same spatial domain shows these spectra to have a universal shape when properly scaled. On the basis of these scaling properties, and a partial exponent spectrum from a still higher-resolution (32x33x32) simulation, it is argued that the actual dimension of the attractor underlying motion of the given computational domain is approximately 780. It is suggested that this periodic turbulent shear flow is deterministic chaos, and that a strange attractor does underly solutions to the Navier-Stokes equations in such flows.

Keefe, Laurence↗

In Silico Human Mobility Data Science: Leveraging Massive Simulated Mobility Data (Vision Paper)

Human mobility data science using trajectories or check-ins of individuals has many applications. Recently, we have seen a plethora of research efforts that tackle these applications. However, research progress in this field is limited by a lack of large and representative datasets. The largest and most commonly used dataset of individual human trajectories captures fewer than 200 individuals, while datasets of individual human check-ins capture fewer than 100 check-ins per city per day. Thus, it is not clear if findings from the human mobility data science community would generalize to large populations. Since obtaining massive, representative, and individual-level human mobility data is hard to come by due to privacy considerations, the vision of this work is to embrace the use of data generated by large-scale socially realistic microsimulations. Informed by both real data and leveraging social and behavioral theories, massive spatially explicit microsimulations may allow us to simulate entire megacities at the person level. The simulated worlds, which do not capture any identifiable personal information, allow us to perform “in silico” experiments using the simulated world as a sandbox in which we have perfect information and perfect control without jeopardizing the privacy of any actual individual. In silico experiments have become commonplace in other scientific domains such as chemistry and biology, permitting experiments that foster the understanding of concepts without any harm to individuals. This work describes challenges and opportunities for leveraging massive and realistic simulated alternate worlds for in silico human mobility data science.

97 MATHEMATICS AND COMPUTING↗

Phase decorrelation, streamwise vortices and acoustic radiation in mixing layers

Several direct numerical simulations were performed and analyzed to study various aspects of the early development of mixing layers. Included are the phase jitter of the large-scale eddies, which was studied using a 2-D spatially-evolving mixing layer simulation; the response of a time developing mixing layer to various spanwise disturbances; and the sound radiation from a 2-D compressible time developing mixing layer.

Ho, C. M.↗

A spatial model of wind shear and turbulence for flight simulation

A three dimensional model which combines measurements of wind shear in the real atmosphere with three dimensional Monte Carlo simulated turbulence was developed. The wind field over the body of an aircraft can be simulated and all aerodynamic loads and moments calculated.

Campbell, C. W.↗

New Approaches to Quantifying Transport Model Error in Atmospheric CO2 Simulations

In recent years, much progress has been made in observing CO2 distributions from space. However, the use of these observations to infer source/sink distributions in inversion studies continues to be complicated by difficulty in quantifying atmospheric transport model errors. We will present results from several different experiments designed to quantify different aspects of transport error using the Goddard Earth Observing System, Version 5 (GEOS-5) Atmospheric General Circulation Model (AGCM). In the first set of experiments, an ensemble of simulations is constructed using perturbations to parameters in the model s moist physics and turbulence parameterizations that control sub-grid scale transport of trace gases. Analysis of the ensemble spread and scales of temporal and spatial variability among the simulations allows insight into how parameterized, small-scale transport processes influence simulated CO2 distributions. In the second set of experiments, atmospheric tracers representing model error are constructed using observation minus analysis statistics from NASA's Modern-Era Retrospective Analysis for Research and Applications (MERRA). The goal of these simulations is to understand how errors in large scale dynamics are distributed, and how they propagate in space and time, affecting trace gas distributions. These simulations will also be compared to results from NASA's Carbon Monitoring System Flux Pilot Project that quantified the impact of uncertainty in satellite constrained CO2 flux estimates on atmospheric mixing ratios to assess the major factors governing uncertainty in global and regional trace gas distributions.

Ott, L.↗

Groundwater Recharge Estimated by Land Surface Models: An Evaluation in the Conterminous U.S.

Estimating diffuse recharge of precipitation is fundamental to assessing groundwater sustainability. Diffuse recharge is also the process through which climate and climate change directly affect groundwater. In this study, we evaluated diffuse recharge over the conterminous U.S. simulated by a suite of land surface models (LSMs) that were forced using a common set of meteorological input data. Simulated annual recharge exhibited spatial patterns that were similar among the LSMs, with the highest values in the eastern U.S. and Pacific Northwest. However, the magnitudes of annual recharge varied significantly among the models and were associated with differences in simulated ET, runoff and snow. Evaluation against two independent datasets did not answer the question of whether the ensemble mean performs the best, due to inconsistency between those datasets. The amplitude and timing of seasonal maximum recharge differed among the models, influenced strongly by model physics governing deep soil moisture drainage rates and, in cold regions, snowmelt. Evaluation using in situ soil moisture observations suggested that true recharge peaks 1-3 months later than simulated recharge, indicating systematic biases in simulating deep soil moisture. However, recharge from lateral flows and through preferential flows cannot be inferred from soil moisture data, and the seasonal cycle of simulated groundwater storage actually compared well with in situ groundwater observations. Long-term trends in recharge were not consistently correlated with either precipitation trends or temperature trends. This study highlights the need to employ dynamic flow models in LSMs, among other improvements, to enable more accurate simulation of recharge.

snowmelt↗

High-Fidelity Simulation of Turbulent Flow Past a Gaussian Bump

A spanwise-periodic computation of a turbulent flow past a Gaussian bump is performed in the form of a hybrid direct numerical simulation and wall-resolved large-eddy simulation. A fourth-order spatially-accurate flow solver is employed to perform the simulation, using 10.2 billion grid points for a Reynolds number of 170000 based on the bump height. The key findings from the simulation are reported in the acceleration and deceleration flow regions associated with the bump shape. Significant anisotropy in the normal Reynolds stresses, along both the wall-normal and streamwise directions, is observed within the acceleration region. The ratio between the Reynolds shear stress and turbulent kinetic energy in that region also experiences significant deviations from the norms of a zero pressure gradient turbulent boundary layer. The chosen Reynolds number generates strong flow separation in the adverse pressure gradient region, which is in contrast with a previous simulation at half the Reynolds number that only indicated incipient separation. An internal layer generated in the acceleration region evolves into a free shear layer that develops in the deceleration region and separates. Proper modeling of this inner layer appears crucial to predict the flow separation. Surface curvature effects on the attached flow development are also discussed.

Turbulent Boundary Layer↗

An Integrated Framework for Memory-Centric Analysis: From Trace Collection to Co-Design

The memory wall phenomenon—where advances in processor performance significantly outpace those in memory subsystems—poses a fundamental challenge for contemporary computing systems. In memory-bound applications, memory subsystem behavior dominates performance, yet existing analysis approaches present significant limitations: detailed microarchitectural simulators require days to weeks to simulate modest workloads; hardware performance counters provide only aggregate statistics that obscure temporal and spatial access patterns; and scaled simulation approaches face challenges in capturing certain behaviors that emerge at larger scales. These limitations reflect a processor-centric design philosophy increasingly misaligned with memory-bound workloads where detailed understanding of memory access patterns, cache hierarchy interactions, and contention is critical for effective optimization. This paper presents an integrated framework for memory-centric analysis that enables effective hardware-software co-design. We describe practical trace collection techniques, including hardware-assisted processor tracing with minimal overhead and portable software-based instrumentation with statistical sampling. We present multi-perspective analysis methods that examine memory behavior from temporal, sequential, spatial, and relational viewpoints, revealing distinct optimization opportunities invisible in aggregate metrics. We detail an architectural modeling framework that uses sampled traces with temporal interpolation and confidence-based filtering to evaluate cache and memory configurations. Evaluation on representative benchmarks demonstrates that this framework achieves practical accuracy (L2 cache errors of 2.64\%, confidence-filtered L3 errors of 9.92\%, bandwidth errors of 7.33\%) while providing substantial speedup (26.8×) over cycle-accurate simulation, enabling rapid design space exploration. We demonstrate how this integrated framework enables systematic identification of both hardware optimizations (memory controller tuning, bank partitioning, NUMA configuration) and software optimizations (data layout restructuring, prefetching strategies, memory-aware scheduling). Through this comprehensive treatment of the memory-centric analysis pipeline—from trace collection through architectural modeling to co-design application—we provide researchers and practitioners with practical techniques for addressing memory bottlenecks in contemporary computing systems.

Gajaria, Dhruv Mayur↗

Evaluation of a Head-Worn Display with Ambient Vision Cues for Unusual Attitude Recovery

A Commercial Aviation Safety Team (CAST) study of 18 loss-of-control events determined that a lack of external visual references was a contributing factor in 17 of these events. CAST recommended that manufacturers should develop and implement virtual day-VMC display systems, such as synthetic vision (SV) or equivalent systems (CAST Safety Enhancement, SE-200). In support of this recommended action, CAST has requested studies to define minimum requirements for virtual day-visual meteorological conditions (VMC) displays to improve flight crew awareness of airplane attitude. NASA’s research in Virtual day-VMC displays, known as synthetic vision systems, are intended to support intuitive flight crew attitude awareness similar to a day-VMC-like environment, especially if they could be designed to create visual dominance. A study was conducted to evaluate the utility of ambient vision (AV) cues paired with virtual Head-Up Display (HUD) symbology on a prototype head-worn display (HWD) during recovery from unusual attitudes in a simulated environment. The virtual-HUD component meets the requirement that the HWD may be used as an equivalent display to the HUD. The presence of AV cueing leverages the potential that a HWD has over the HUD for spatial disorientation prevention. The simulation study was conducted as a single-pilot operation, under realistic flight scenarios, with off-nominal events occurring that were capable of inducing unusual attitudes. Independent variables of the experiment included: 1) AV capability (on vs off) 2) AV display opaqueness (transparent vs opaque) and display location (HWD vs traditional head- down displays); AV cues were only present when the HWD was being worn by the subject pilot.

Nicholas, Stephanie N.↗

AVIRIS and TIMS data processing and distribution at the land processes distributed active archive center

The U.S. Government has initiated the Global Change Research program, a systematic study of the Earth as a complete system. NASA's contribution of the Global Change Research Program is the Earth Observing System (EOS), a series of orbital sensor platforms and an associated data processing and distribution system. The EOS Data and Information System (EOSDIS) is the archiving, production, and distribution system for data collected by the EOS space segment and uses a multilayer architecture for processing, archiving, and distributing EOS data. The first layer consists of the spacecraft ground stations and processing facilities that receive the raw data from the orbiting platforms and then separate the data by individual sensors. The second layer consists of Distributed Active Archive Centers (DAAC) that process, distribute, and archive the sensor data. The third layer consists of a user science processing network. The EOSDIS is being developed in a phased implementation. The initial phase, Version 0, is a prototype of the operational system. Version 0 activities are based upon existing systems and are designed to provide an EOSDIS-like capability for information management and distribution. An important science support task is the creation of simulated data sets for EOS instruments from precursor aircraft or satellite data. The Land Processes DAAC, at the EROS Data Center (EDC), is responsible for archiving and processing EOS precursor data from airborne instruments such as the Thermal Infrared Multispectral Scanner (TIMS), the Thematic Mapper Simulator (TMS), and Airborne Visible and Infrared Imaging Spectrometer (AVIRIS). AVIRIS, TIMS, and TMS are flown by the NASA-Ames Research Center ARC) on an ER-2. The ER-2 flies at 65000 feet and can carry up to three sensors simultaneously. Most jointly collected data sets are somewhat boresighted and roughly registered. The instrument data are being used to construct data sets that simulate the spectral and spatial characteristics of the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) instrument scheduled to be flown on the first EOS-AM spacecraft. The ASTER is designed to acquire 14 channels of land science data in the visible and near-IR (VNIR), shortwave-IR (SWIR), and thermal-IR (TIR) regions from 0.52 micron to 11.65 micron at high spatial resolutions of 15 m to 90 m. Stereo data will also be acquired in the VNIR region in a single band. The AVIRIS and TMS cover the ASTER VNIR and SWIR bands, and the TIMS covers the TIR bands. Simulated ASTER data sets have been generated over Death Valley, California, Cuprite, Nevada, and the Drum Mountains, Utah using a combination of AVIRIS, TIMS, amd TMS data, and existing digital elevation models (DEM) for the topographic information.

Mah, G. R.↗

A time-parallel method for scalable heat transfer simulations of additive manufacturing

Here, a major challenge in simulating the thermal behavior in additive manufacturing processes is the disparate length and time scales between transport phenomena occurring in the melt pool and the component. A common simulation approach relies on spatial decomposition for parallel computing, but due to the nature of heat transfer in AM, where most of the computational expenditure is localized near the melt pool, the computational speedup from spatial parallelization saturates quickly. Therefore, additional parallelism by means of time-domain decomposition is needed to fully take advantage of high-performance computing (HPC) resources. This work introduces a time-parallel method to improve the computational scalability of additive manufacturing simulations on HPC systems, while maintaining high temporal resolution of heat transfer near the melt pool. The method, inspired by the nonlinear paraexp formalism, performs an iterative superposition of nonlinear solutions to the initial value problem, integrating the heat equation across overlapping time-parallel intervals. For a single layer of the NIST AMB2018–01 L7 benchmark problem, the method achieves a 38.51x speedup in wall-clock time with a maximum error in the global temperature solution of 0.99%. This reduces the total solution time from 196.72 min to 5.11 min on 128 nodes of the ORNL Frontier supercomputer. The tradeoff between accuracy and total wall-clock time is investigated and recommendations for time-parallel deployment for AM problems are made.

Additive manufacturing↗

Has Reducing Ship Emissions Brought Forward Global Warming?

Abstract Ships brighten low marine clouds from emissions of sulfur and aerosols, resulting in visible “ship tracks”. In 2020, new shipping regulations mandated an ∼80% reduction in the allowed fuel sulfur content. Recent observations indicate that visible ship tracks have decreased. Model simulations indicate that since 2020 shipping regulations have induced a net radiative forcing of +0.12 Wm −2 . Analysis of recent temperature anomalies indicates Northern Hemisphere surface temperature anomalies in 2022–2023 are correlated with observed cloud radiative forcing and the cloud radiative forcing is spatially correlated with the simulated radiative forcing from the 2020 shipping emission changes. Shipping emissions changes could be accelerating global warming. To better constrain these estimates, better access to ship position data and understanding of ship aerosol emissions are needed. Understanding the risks and benefits of emissions reductions and the difficultly in robust attribution highlights the large uncertainty in attributing proposed deliberate climate intervention.

54 ENVIRONMENTAL SCIENCES↗

Kinetic study of shock formation and particle acceleration in laser-driven quasi-parallel magnetized collisionless shocks

Quasi-parallel magnetized collisionless shocks are believed to be one of the most efficient accelerators in the universe. Compared to quasi-perpendicular shocks, quasi-parallel shocks are more difficult to form in the laboratory and to simulate because of their large spatial scales and long formation times. Our two-dimensional particle-in-cell simulations show that the early stages of quasi-parallel shock formation are achievable in experiments planned for the National Ignition Facility and that particles accelerated by diffusive shock acceleration (DSA) are expected to be observable in the experiment. Repetitive ion acceleration by crossings of the shock front, a key feature of DSA, is seen in the simulations. Other characteristic features of quasi-parallel shocks such as upstream wave excitation by energetic ions are also observed, and energy partition between the ions and the electrons in the downstream of the shock is briefly discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Selecting the spatial resolution of satellite sensors required for global monitoring of land transformations

The paper provides preliminary evidence for the spatial resolutions required to monitor land transformations at broad scales. This is obtained from simulations of imagery at various spatial resolutions between 125 and 4000 m derived from Landsat MSS imagery. Consideration is given to the various types of spatial images detectable by remotely-sensed systems, as well as to the difficulties associated in disentangling permanent land transformations from shorter term changes such as phenological and interannual changes.

Townshend, J. R. G.↗

Modeling of high speed chemically reacting flow-fields

The SPARK3D and SPARK3D-PNS computer programs were developed to model 3-D supersonic, chemically reacting flow-fields. The SPARK3D code is a full Navier-Stokes solver, and is suitable for use in scramjet combustors and other regions where recirculation may be present. The SPARK3D-PNS is a parabolized Navier-Stokes solver and provides an efficient means of calculating steady-state combustor far-fields and nozzles. Each code has a generalized chemistry package, making modeling of any chemically reacting flow possible. Research activities by the Langley group range from addressing fundamental theoretical issues to simulating problems of practical importance. Algorithmic development includes work on higher order and upwind spatial difference schemes. Direct numerical simulations employ these algorithms to address the fundamental issues of flow stability and transition, and the chemical reaction of supersonic mixing layers and jets. It is believed that this work will lend greater insight into phenomenological model development for simulating supersonic chemically reacting flows in practical combustors. Currently, the SPARK3D and SPARK3D-PNS codes are used to study problems of engineering interest, including various injector designs and 3-D combustor-nozzle configurations. Examples, which demonstrate the capabilities of each code are presented.

Drummond, J. P.↗

Crack Imaging and Quantification in Aluminum Plates with Guided Wave Wavenumber Analysis Methods

Guided wavefield analysis methods for detection and quantification of crack damage in an aluminum plate are presented in this paper. New wavenumber components created by abrupt wave changes at the structural discontinuity are identified in the frequency-wavenumber spectra. It is shown that the new wavenumbers can be used to detect and characterize the crack dimensions. Two imaging based approaches, filter reconstructed imaging and spatial wavenumber imaging, are used to demonstrate how the cracks can be evaluated with wavenumber analysis. The filter reconstructed imaging is shown to be a rapid method to map the plate and any existing damage, but with less precision in estimating crack dimensions; while the spatial wavenumber imaging provides an intensity image of spatial wavenumber values with enhanced resolution of crack dimensions. These techniques are applied to simulated wavefield data, and the simulation based studies show that spatial wavenumber imaging method is able to distinguish cracks of different severities. Laboratory experimental validation is performed for a single crack case to confirm the methods' capabilities for imaging cracks in plates.

Yu, Lingyu↗

Estimation of discontinuous coefficients in parabolic systems - Applications to reservoir simulation

Spline-based techniques for estimating spatially varying parameters that appear in parabolic distributed systems (typical of those found in reservoir simulation problems) are presented. In particular, the problem of determining discontinuous coefficients is discussed, estimating both the functional shape and points of discontinuity for such parameters. In addition, the ideas may also be applied to problems with unknown initial conditions and unknown parameters appearing in terms representing external forces. Convergence results and a summary of numerical performance of the resulting algorithms are given.

Lamm, Patricia K.↗