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 181 records · Page 10

Solar seeing and the spatial properties of the five-minute oscillations.

Numerical simulation of observations of the spatial properties of the five-minute oscillations, assuming that the oscillations are internal gravity waves excited by granular convection according to the theory of Thomas et al. (1971). The simulation includes the effects of seeing and finite aperture. The details of the simulation are chosen to model the observational method of Frazier (1968). The results show that the peak in the observed power spectrum of the oscillations can occur at a wavelength considerably longer than the true wavelength of the oscillations. In particular, the peak in Frazier's observed power spectra at wavelength of about 5000 km is consistent with the considerably shorter true wavelength of about 1500 km predicted by the gravity-wave theory.

Thomas, J. H.↗

Global 1 km land surface parameters for kilometer-scale Earth system modeling

Earth system models (ESMs) are progressively advancing towards the kilometer scale (“k-scale”). However, the surface parameters for land surface models (LSMs) within ESMs running at the k-scale are typically derived from coarse-resolution and outdated datasets. This study aims to develop a new set of global land surface parameters with a resolution of 1 km for multiple years from 2001 to 2020, utilizing the latest and most accurate available datasets. Specifically, the datasets consist of parameters related to land use and land cover, vegetation, soil, and topography. Differences between the newly developed 1 km land surface parameters and conventional parameters emphasize their potential for higher accuracy due to the incorporation of the most advanced and latest data sources. To demonstrate the capability of these new parameters, we conducted 1 km resolution simulations using the E3SM Land Model version 2 (ELM2) over the contiguous United States. Our results demonstrate that land surface parameters contribute to significant spatial heterogeneity in ELM2 simulations of soil moisture, latent heat, emitted longwave radiation, and absorbed shortwave radiation. On average, about 31 % to 54 % of spatial information is lost by upscaling the 1 km ELM2 simulations to a 12 km resolution. Using eXplainable Machine Learning (XML) methods, the influential factors driving the spatial variability and spatial information loss of ELM2 simulations were identified, highlighting the substantial impact of the spatial variability and information loss of various land surface parameters, as well as the mean climate conditions. The comparison against four benchmark datasets indicates that ELM generally performs well in simulating soil moisture and surface energy fluxes. The new land surface parameters are tailored to meet the emerging needs of k-scale LSM and ESM modeling with significant implications for advancing our understanding of water, carbon, and energy cycles under global change. The 1 km land surface parameters are publicly available at https://doi.org/10.5281/zenodo.10815170 (Li et al., 2024).

54 ENVIRONMENTAL SCIENCES↗

Simulations of the gross solar magnetic field during sunspot cycle 21

Regarding new bipolar magnetic regions as sources of flux, the evolution of the radial component of the solar photospheric magnetic field during 1976-1984 has been simulated with a spatial resolution of about 34,000 km, and the corresponding evolution of its absolute value averaged over the visible disk is derived. For nominal values of the transport parameters, this simulated gross field is in close, though imperfect, agreement with the observed gross field and its associated indices of solar activity. By analyzing the response of the simulated gross field to variations in the transport parameters and the source properties, it is found that the simulated field originates in newly erupted bipolar regions. The lifetimes of these regions are almost always less than three mo. Consequently, the strength of the simulated gross field is a measure of the current level of solar activity, and any recurrent patterns with lifetimes in excess of six mo must reflect the continuing eruption of new flux at 'active longitudes' rather than the persistence of old flux in long-lived magnetic structures.

Sheeley, Neil R., Jr.↗

Navier-Stokes Simulation of UH-60A Rotor/Wake Interaction Using Adaptive Mesh Refinement

Time-dependent Navier-Stokes simulations have been carried out for a flexible UH-60A rotor in forward flight, where the rotor wake interacts with the rotor blades. These flow conditions involved blade vortex interaction and dynamic stall, two common conditions that occur as modern helicopter designs strive to achieve greater flight speeds and payload capacity. These numerical simulations utilized high-order spatial accuracy and delayed detached eddy simulation. Emphasis was placed on understanding how improved rotor wake resolution affects the prediction of the normal force, pitching moment, and chord force of the rotor. Adaptive mesh refinement was used to highly resolve the turbulent rotor wake in a computationally efficient manner. Moreover, blade vortex interaction was found to trigger dynamic stall. Time-dependent flow visualization was utilized to provide an improved understanding of the numerical and physical mechanisms involved with three-dimensional dynamic stall.

UH-60A RotorWake↗

Large-Eddy Simulation of Boundary Layer Transition on Swept Wings

The large-eddy simulation of the spatial evolution of a stationary crossflow vortex packet in a three-dimensional boundary layer was performed. Although a coarse grid was used (compared to that required by a direct numerical simulation) the essential features of the disturbance evolution, such as the spanwise disturbance spreading and the vortex rollover, were captured accurately. The eddy viscosity became significant only in the late nonlinear stages of the simulation.

Huai, Xiaoli↗

Onset of Reconnection in the near Magnetotail: PIC Simulations

Using 2.5-dimensional particle-in-cell (PIC) simulations of magnetotail dynamics, we investigate the onset of reconnection in two-dimensional tail configurations with finite Bz. Reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. We found a clear distinction between stable and unstable cases, dependent on deformation amplitude and ion/electron mass ratio. The threshold appears consistent with electron tearing. The evolution prior to onset, as well as the evolution of stable cases, are largely independent of the mass ratio, governed by integral flux tube entropy conservation as imposed in MHD (magnetohydrodynamics). This suggests that ballooning instability in the tail should not be expected prior to the onset of tearing and reconnection. The onset time and other onset properties depend on the mass ratio, consistent with expectations for electron tearing. At onset,we found electron anisotropies T⊥⁄ T∥ (bottom tail divided by parallel tail) equals 1.1-1.3, raising growth rates and wavenumbers. Our simulations have provided a quantitative onset criterion that is easily evaluated in MHD simulations, provided the spatial resolution is sufficient. The evolution prior to onset and after the formation of a neutral line does not depend on the electron physics, which should permit an approximation by MHD simulations with appropriate dissipation terms.

PIC↗

Effect of Spatial Filtering in Implicit Large-Eddy Simulations of Separated Flows

The relatively high Reynolds number of turbulent flows encountered in various applications puts these problems well beyond the reach of direct numerical simulation (DNS) at present. Meanwhile, lower-fidelity Reynolds-averaged Navier-Stokes (RANS) calculations are known to be not accurate enough in complex problems, such as smooth-body flow separation and other flows involving highly-unsteady phenomena. Hence, given the current infeasibility of DNS and the unsatisfactory performance of RANS, intermediate techniques such as large-eddy simulation (LES) and hybrid RANS-LES, whose fidelity lie between RANS and DNS, have received much attention for application to various problems of practical importance. Modeling of the effect of missing scales on resolved scales, also known as subgrid-scale (SGS) modeling, is an important subject for LES. SGS models can be broadly categorized as explicit or implicit approaches. The explicit approach is based on an SGS model that explicitly appears in the governing equations expressed in the form of so-called “filtered Navier-Stokes equations”, which describe the evolution of the turbulence scales resolved by the LES grid. The effect of the scales unresolved by the grid is represented by the SGS model. The implicit modeling approach, on the other hand, does not employ an explicit model but instead treats the intrinsic dissipation of the numerical discretization scheme as an implicit SGS model. An LES without an explicit SGS model is commonly termed as an implicit LES (ILES). The relative merits of one SGS modeling approach over another is a subject of ongoing debate. We have opted to employ an ILES methodology, based on high-order compact finite-difference and spatial filtering schemes, in our recent investigations of separated flow problems [1, 2]. Further discussion of our preference of ILES over explicit LES is provided in Uzun and Malik [2]. The spatial filtering operation, described in the next section, is treated as an implicit SGS model for the ILES. Some observations made during the course of our recent investigations, which pointed out to excessive numerical dissipation in certain parts of the flowfield, prompted us to take a closer look at the potential effect of the spatial filter on ILES predictions. This technical note is therefore devoted to spatial filter effects in the context of a high Reynolds number, transonic shock-induced separated flow.

Ali Uzun↗

Development of Computational Materials Workflows for Additively Manufactured Metallic Materials to Enable Accelerated Prediction of Fatigue Performance

The maturation of computational materials approaches for fatigue performance prediction in a qualification and certification process is stifled by the ability to validate complex, microstructure-based simulations. Such a validation strategy bears immediate challenges including generating accurate virtual microstructures, efficiently solving physics-based mechanical simulations over relevant spatial and temporal scales, and acquiring high-fidelity calibration and validation data at the appropriate length scale. This presentation will overview these common challenges and present a case study to demonstrate a computational materials workflow for additively manufactured metallic materials. In this study, process-specific defects are characterized using segmented X-Ray micro-computed tomography measurements and overlaid on virtual microstructures. Accelerated crystal plasticity-based fatigue simulations are performed to demonstrate cyclic evolution and localization of mechanical fields in the vicinity of defects in response to their precise spatial configuration. An example of how this computational materials workflow may support next-generation qualification is discussed.

computational materials↗

Multizone Modeling of Black Hole Accretion and Feedback in 3D GRMHD: Bridging Vast Spatial and Temporal Scales

Simulating accretion and feedback from the horizon scale of supermassive black holes (SMBHs) out to galactic scales is challenging because of the vast range of scales involved. Elaborating on H. Cho et al., we describe and test a "multizone" technique, which is designed to tackle this difficult problem in three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. While short-timescale variability should be interpreted with caution, the method is demonstrated to be well-suited for finding dynamical steady states over a wide dynamic range. We simulate accretion on a nonspinning SMBH ($a\ast$ = 0) using initial conditions and the external galactic potential from a large-scale galaxy simulation and achieve a steady state over eight decades in radius. As found in H. Cho et al., the density scales with radius as ρ ∝ r –1 inside the Bondi radius R B , which is located at R B = 2 × 10 5 r g (≈60 pc for M87), where r g is the gravitational radius of the SMBH; the plasma-β is ~ unity, indicating an extended magnetically arrested state; the mass accretion rate $\dot{M}$ is ≈1% of the analytical Bondi accretion rate ${\dot{M}}_{{\rm{B}}};$ and there is continuous energy feedback out to ≈100R B (or beyond > kpc) at a rate $\approx 0.02\dot{M}{c}^{2}$. Surprisingly, no ordered rotation in the external medium survives as the magnetized gas flows to smaller radii, and the final steady solution is very similar to when the exterior has no rotation. Using the multizone method, we simulate GRMHD accretion over a wide range of Bondi radii, R B ~ 10 2 –10 7 r g , and find that $\dot{M}/{\dot{M}}_{{B}}\approx {({R}_{{B}}/6\,{r}_{g})}^{-0.5}$.

79 ASTRONOMY AND ASTROPHYSICS↗

Violation of Field Line Conservation and Associated Spatial Scales in Particle-In-Cell Simulations and MMS Data

For the first time, space flight technology exists to detect, in situ, violation of magnetic field line conservation. The violation of magnetic line conservation on scales smaller than the system size is a necessary and sufficient condition for finite magnetic field reconnection. We demonstrate that violation of line conservation produces a detectable, structured signature both in particle-in-cell simulations of reconnection and in data from the Magnetospheric Multi-Scale mission. In particle-in-cell simulations of asymmetric reconnection, the quantity-which we call M-that identifies this violation achieves significant values in electron skin depth-scale layers that extend away from the separator, with higher values emerging on the low-density, high-magnetic-field side of the current sheet. At the separator, M owes largely to perpendicular gradients in the parallel electric field, while it attains its highest values away from the separator in dispersed, layered structures associated with gradients in the perpendicular nonideal electric field and electron transport. Sub-ion scale bipolar forms of the quantity also appear further from the separator. In two MMS burst data intervals detecting the electron diffusion region, we find that M exceeds measurement uncertainties both at the separator and near the separatrices. One interval has highly sheared reconnecting fields and the other a stronger guide field. For one event, we determine the location and scale of M and the inner electron diffusion region relative to electron outflows and the magnetic separatrices. The measure can therefore serve as a potent diagnostic for magnetic reconnection in space measurements.

Wendel, Deirdre E.↗

Numerical simulation of transition in boundary layers

Different implications of employment of the temporal or spatial model in direct numerical simulations of transition based on the complete Navier-Stokes equations are outlined, and it is noted that in situations where approximations are not justified, the temporal model cannot be employed. The discussion is then concentrated on developing numerical methods for solving the Navier-Stokes equations based on the spatial model that are applicable for realistic situations for which the temporal model would not be valid. It is concluded that in the future numerical simulations based on the spatial approach will replace, to a large extent, the temporal simulations.

Fasel, H.↗

Impacts of aerosol size-dependent below-cloud scavenging on tropospheric aerosol in the NASA GEOS model

Cloud scavenging is the major sink for a suite of tropospheric aerosols, and thus largely affects the global distribution and lifetime of aerosols as well as the formation of clouds. Simplified model parameterizations of cloud scavenging of aerosols substantially contribute to large uncertainties in the simulated aerosol burdens, spatial distributions, lifetimes, as well as aerosol radiative forcing. Here we evaluate the wet scavenging scheme in the current NASA GEOS model, develop more physically-based parameterizations of below-cloud scavenging, and assess the impacts of the new scavenging parameterizations on the GEOS simulations of aerosols. The current below-cloud scavenging (BCS; washout) scheme uses a first-order removal with a constant scavenging coefficient. We replaced the current scheme with a size-dependent parameterization of Croft et al.(Atmos. Chem. Phys., 2009; hereafter referred to as C09). The C09 BCS parameterizations calculate aerosol below-cloud scavenging coefficients based on a size-dependent collision efficiency between aerosol particles and raindrops (or snow crystals). A look-up table of compiled collision efficiencies provides values for sixty aerosol sizes and nine rain flux rates. Snow-aerosol collision efficiencies are a function of the sixty aerosol sizes only. We assume a log-normal size distribution for aerosols (bulk) in GOCART. The collision efficiencies are then calculated online using a bilinear interpolation from the look-up table for the given aerosol sizes and precipitation flux rates in the model.

aerosols↗

Three-Dimensional Grid Visualization for Planning Activities: A Dubai Case Study

National Laboratory of the Rockies (NLR), in collaboration with the Dubai Electricity and Water Authority (DEWA) and Infra-X, has undertaken the Energy Visualization Analysis Project. The aim of this project is to enhance analytical and 3D visualization capabilities for distribution network planning and renewable energy integration. As modern grid continues to evolve with large-scale solar PV deployment and emerging distributed energy resources (DERs), the ability to effectively analyze, visualize, and communicate complex grid behaviors has become increasingly critical. The project focuses on developing empirical use cases based on real distribution feeder data and engineering workflows, ensuring the outcomes are directly aligned with operational environment. Through time-series power flow simulations and nodal hosting capacity analysis, the study quantifies the impacts of high PV penetration on voltage and thermal limits within representative 11 kV feeders. These analyses identify specific nodes and conditions where DER integration challenges arise. Furthermore, a Battery Energy Storage System (BESS) optimization algorithm was applied to determine the optimal size and placement of storage systems that can mitigate network constraints and enhance hosting capacity. The comparative results between base-case and BESS-augmented scenarios clearly demonstrate improvements in network stability and load management efficiency. In parallel, the NLR team developed an immersive 3D visualization framework, enabling interactive exploration of grid simulations using commodity head-mounted display (HMD) systems. This framework transforms conventional 2D simulation data into spatially intuitive visual environments - allowing engineers to analyze feeder conditions, PV hosting potential, and BESS effects in real time. This report represents the first foundational phase in establishing a visualization-driven analytical ecosystem. It provides a methodological foundation for data integration, visualization architecture, and simulation-based decision support, paving the way for large-scale adoption of immersive visualization across DEWA's Smart Grid Initiative, R&D activities, and future network resilience studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Evolving the SpaceFOM: Lessons Learned and Future Development

For over fifty years, simulation has been a cornerstone technology for space missions. In many cases, one single simulator, developed by one team, is not enough to meet the simulation objectives. Astronaut crew training, docking, vehicle testing and complex mission planning all require several simulations to be connected and to interoperate, in order to achieve the goal of the simulation. To facilitate the development of interoperable simulations, an open standard, the SISO SpaceFOM was developed in collaboration between government, industry and academia. The standard was released in 2020. The initial release of SpaceFOM focuses on the core interoperability topics for space simulation: (i) specifying spatial data using several, well defined reference frames, (ii) managing time (including real-time, hard real-time, and as-fast-as-possible execution), (iii) composition of systems with subsystems such as space vehicles, and (iv) execution control, such as initializing, running, pausing, and terminating simulations. A number of requirements for the next version of SpaceFOM have been identified and proposed solutions are already under development. One key area is the object classes used for data exchange, where a richer set of classes have been proposed, based on practical use cases. For time management, the use of variable time resolutions and variable time steps is also proposed. For execution control, save and restore of simulations is considered as well as improved robustness, to handle failing components. A number of clarifications are also under consideration. This paper provides both a forward-looking view of the technical developments underway and a retrospective on the standard's evolution. The authors emphasize that prioritizing core challenges and delivering a practical first version were crucial steps. While not every idea was implemented in the initial release, the timely delivery of a functional standard enabled real-world applications, the growth of a user community, and the collection of valuable insights to shape future versions. The authors look forward to working with the next version of SpaceFOM together with the growing user community.

FOMs↗

Rethinking materials simulations: Blending direct numerical simulations with neural operators

Abstract Materials simulations based on direct numerical solvers are accurate but computationally expensive for predicting materials evolution across length- and time-scales, due to the complexity of the underlying evolution equations, the nature of multiscale spatiotemporal interactions, and the need to reach long-time integration. We develop a method that blends direct numerical solvers with neural operators to accelerate such simulations. This methodology is based on the integration of a community numerical solver with a U-Net neural operator, enhanced by a temporal-conditioning mechanism to enable accurate extrapolation and efficient time-to-solution predictions of the dynamics. We demonstrate the effectiveness of this hybrid framework on simulations of microstructure evolution via the phase-field method. Such simulations exhibit high spatial gradients and the co-evolution of different material phases with simultaneous slow and fast materials dynamics. We establish accurate extrapolation of the coupled solver with large speed-up compared to DNS depending on the hybrid strategy utilized. This methodology is generalizable to a broad range of materials simulations, from solid mechanics to fluid dynamics, geophysics, climate, and more.

36 MATERIALS SCIENCE↗

Modeling Macro- and Micro-Scale Turbulent Mixing and Chemistry in Engine Exhaust Plumes

Simulation of turbulent mixing and chemical processes in the near-field plume and plume-vortex regimes has been successfully carried out recently using a reduced gas phase kinetics mechanism which substantially decreased the computational cost. A detailed mechanism including gas phase HOx, NOx, and SOx chemistry between the aircraft exhaust and the ambient air in near-field aircraft plumes is compiled. A reduced mechanism capturing the major chemical pathways is developed. Predictions by the reduced mechanism are found to be in good agreement with those by the detailed mechanism. With the reduced chemistry, the computer CPU time is saved by a factor of more than 3.5 for the near-field plume modeling. Distributions of major chemical species are obtained and analyzed. The computed sensitivities of major species with respect to reaction step are deduced for identification of the dominant gas phase kinetic reaction pathways in the jet plume. Both the near field plume and the plume-vortex regimes were investigated using advanced mixing models. In the near field, a stand-alone mixing model was used to investigate the impact of turbulent mixing on the micro- and macro-scale mixing processes using a reduced reaction kinetics model. The plume-vortex regime was simulated using a large-eddy simulation model. Vortex plume behind Boeing 737 and 747 aircraft was simulated along with relevant kinetics. Many features of the computed flow field show reasonable agreement with data. The entrainment of the engine plumes into the wing tip vortices and also the partial detrainment of the plume were numerically captured. The impact of fluid mechanics on the chemical processes was also studied. Results show that there are significant differences between spatial and temporal simulations especially in the predicted SO3 concentrations. This has important implications for the prediction of sulfuric acid aerosols in the wake and may partly explain the discrepancy between past numerical studies (that employed parabolic or temporal approximations) and the measured data. Finally to address the major uncertainty in the near-field plume modeling related to the plume processing of sulfur compounds and advanced model was developed to evaluate its impact on the chemical processes in the near wake. A comprehensive aerosol model is developed and it is coupled with chemical kinetics and the axisymmetric turbulent jet flow models. The integrated model is used to simulate microphysical processes in the near-field jet plume, including sulfuric acid and water binary homogeneous nucleation, coagulation, non-equilibrium heteromolecular condensation, and sulfur-induced soot activation. The formation and evolution of aerosols are computed and analyzed. The computed results show that a large number of ultra-fine (0.3--0.6 nm in radius) volatile HSO4 - HO embryos are generated in the near-field plume. These embryos further grow in size by self coagulation and condensation. Soot particles can be activated by both heterogeneous nucleation and scavenging of H2SO4-H2O aerosols. These activated soot particles can serve as water condensation nuclei for contrail formation. Conditions under which ice contrails can form behind aircrafts are studied. The sensitivities of the threshold temperature for contrail formation with respect to aircraft propulsion efficiency, relative humidity, and ambient pressure are evaluated. The computed aerosol properties for different extent of fuel sulfur conversion to S(VI) (SO3 and H2SO4) in engine are examined and the results are found to be sensitive to this conversion fraction.

Menon, Suresh↗

Modeling Being "Lost": Imperfect Situation Awareness

Being "lost" is an exemplar of imperfect Situation Awareness/Situation Understanding (SA/SU) -- information/knowledge that is uncertain, incomplete, and/or just wrong. Being "lost" may be a geo-spatial condition - not knowing/being wrong about where to go or how to get there. More broadly, being "lost" can serve as a metaphor for uncertainty and/or inaccuracy - not knowing/being wrong about how one fits into a larger world view, what one wants to do, or how to do it. This paper discusses using agent based modeling (ABM) to explore imperfect SA/SU, simulating geo-spatially "lost" intelligent agents trying to navigate in a virtual world. Each agent has a unique "mental map" -- its idiosyncratic view of its geo-spatial environment. Its decisions are based on this idiosyncratic view, but behavior outcomes are based on ground truth. Consequently, the rate and degree to which an agent's expectations diverge from ground truth provide measures of that agent's SA/SU.

Middleton, Victor E.↗