Search NASA⌕ Search

SEARCH · Search NASA

Results for “space grid”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Flexible User-Defined Domain Decomposition in Kilometer-Scale E3SM Land Model Simulation

The Energy Exascale Earth System Model (E3SM) Land Model (ELM) has been extended to kilometer-scale (km-ELM) resolutions, enabling high-fidelity simulations of terrestrial processes at 1 km x 1 km grid spacing. In ELM, domain decomposition partitions the computational domain across processors, ensuring efficient parallel execution. Currently, round-robin decomposition is applied, providing a straightforward way to distribute computational workload. As ELM continues evolving at the kilometer-scale (km-scale), particularly with integrating lateral flow modeling, decomposition strategies must also account for the increased workload and data movement. This paper introduces a flexible user-defined domain decomposition framework, allowing users to customize domain partitioning based on application requirements. The impact of different decomposition strategies is evaluated across various applications concerning computation, communication, and I/O. Results demonstrate that while 1D partitioning yields superior I/O performance, k-nearest neighbors (KNN) clustering effectively reduces inter-process communication overhead. This study lays the groundwork for scalable partitioning in large-scale land surface simulations, enhancing next-generation Earth system modeling.

Wang, Dali [ORNL] (ORCID:0000000168065108)↗
DOE OSTI2020

Robust late twenty-first century shift in the regional monsoons in RegCM-CORDEX simulations

Abstract We use an unprecedented ensemble of regional climate model (RCM) projections over seven regional CORDEX domains to provide, for the first time, an RCM-based global view of monsoon changes at various levels of increased greenhouse gas (GHG) forcing. All regional simulations are conducted using RegCM4 at a 25 km horizontal grid spacing using lateral and lower boundary forcing from three General Circulation Models (GCMs), which are part of the fifth phase of the Coupled Model Inter-comparison Project (CMIP5). Each simulation covers the period from 1970 through 2100 under two Representative Concentration Pathways (RCP2.6 and RCP8.5). Regional climate simulations exhibit high fidelity in capturing key characteristics of precipitation and atmospheric dynamics across monsoon regions in the historical period. In the future period, regional monsoons exhibit a spatially robust delay in the monsoon onset, an increase in seasonality, and a reduction in the rainy season length at higher levels of radiative forcing. All regions with substantial delays in the monsoon onset exhibit a decrease in pre-monsoon precipitation, indicating a strong connection between pre-monsoon drying and a shift in the monsoon onset. The weakening of latent heat driven atmospheric warming during the pre-monsoon period delays the overturning of atmospheric subsidence in the monsoon regions, which defers their transitioning into deep convective states. Monsoon changes under the RCP2.6 scenario are mostly within the baseline variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Future projections in the climatology of global low-level jets from CORDEX-CORE simulations

The potential changes in the strength and location of five low-level jets (LLJs) located within four Coordinated Regional Climate Downscaling Experiment (CORDEX) domains are examined for present and future climate conditions using an ensemble of simulations conducted with the RegCM4 regional model at a 25 km horizontal grid spacing. Lateral and lower boundary forcing fields are from three General Circulation Models (GCMs), and we analyse a historical period (1995–2014) along with two future periods (2041–2060 and 2080–2099) under the Representative Concentration Pathways 2.6 and 8.5. The RegCM4, as driven by the GCMs, is capable of capturing most of the observed climatological features of the LLJs, both in terms of spatial location and seasonal evolution. Analysis of the influence of global warming on the LLJs shows a consistent strengthening of the jets and a shift in their location under both warming scenarios. The Monsoon and West African westerly LLJs exhibit a northward shift, while the Caribbean and South American LLJs present a westward expansion. The use of an ensemble of high-resolution simulations is found to provide a key element for a robust assessment of changes in LLJs associated with future global warming scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Future projections in tropical cyclone activity over multiple CORDEX domains from RegCM4 CORDEX-CORE simulations

The characteristics of tropical cyclone (TC) activity over five TC basins lying within four Coordinated Regional Downscaling Experiment (CORDEX) domains are examined for present and future climate conditions using a new ensemble of simulations completed as part of the CORDEX-CORE initiative with the regional climate model RegCM4. The simulations are conducted at a 25 km horizontal grid spacing and are driven by three CMIP5 general circulation models (GCMs) under two Representative Concentration Pathways (RCP2.6 and RCP8.5). The RegCM4 captures most features of the observed TC climatology, except for the TC intensity, which is thus statistically adjusted using a bias correction procedure to account for the effect of the coarse model resolution. The RegCM4 exhibits an improved simulation of several TC statistics compared to the driving GCMs, over most basins analyzed. In future climate conditions we find significant increases in TC frequency over the North Indian Ocean, the Northwest Pacific and Eastern Pacific regions, which are consistent with an increase in mid-tropospheric relative humidity. The North Atlantic and Australasia regions show a decrease in TC frequency, mostly associated with an increase in wind shear. We also find a consistent increase in future storm rainfall rates associated with TCs and in the frequency of the most intense TCs over most domains. Finally, our study shows robust responses often, but not always, in line with previous studies, still implying the presence of significant uncertainties in the projection of TC characteristics, which need to be addressed using large ensembles of simulations with high-resolution models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Anticipating how rain-on-snow events will change through the 21st century: lessons from the 1997 new year’s flood event

The California-Nevada 1997 New Year’s flood was an atmospheric river (AR)-driven rain-on-snow (RoS) event and remains the costliest in their history. The joint occurrence of saturated soils, rainfall, and snowmelt generated inundation throughout northern California-Nevada. Although AR RoS events are projected to occur more frequently with climate change, the warming sensitivity of their flood drivers across scales remains understudied. We leverage the regionally refined mesh capabilities of the Energy Exascale Earth System Model (RRM-E3SM) to recreate the 1997 New Year’s flood with horizontal grid spacings of 3.5 km across California, with forecast lead times of up to 4 days, and across six warming levels ranging from pre-industrial conditions to +3.5° C. We describe the sensitivity of the flood drivers to warming including AR duration and intensity, precipitation phase, intensity and efficiency, snowpack mass and energy changes, and runoff efficiency. Our findings indicate current levels of climate change negligibly influence the flood drivers. At warming levels ≥ 1.7° C, AR hazard potential increases, snowpack nonlinearly decreases, antecedent soil moisture decreases (except where the snowline retreats), and runoff decreases (except in the southern Sierra Nevada where antecedent snowpack persists). Storm total precipitation increases, but at rates below warming-induced increases in saturation-specific humidity. Warming intensifies short-duration, high-intensity rainfall, particularly where snowfall-to-rainfall transitions occur. This study highlights the nonlinear tradeoffs in 21st-century RoS flood hazards with warming and provides water management and infrastructure investment adaptation considerations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Comparison of Observations and Predictions of Daytime Planetary-Boundary-Layer Heights and Surface Meteorological Variables in the Columbia River Gorge and Basin During the Second Wind Forecast Improvement Project

The second Wind Forecast Improvement Project (WFIP2) is an 18-month field campaign in the Pacific Northwest U.S.A., whose goal is to improve the accuracy of numerical-weather-prediction forecasts in complex terrain. The WFIP2 campaign involved the deployment of a large suite of in situ and remote sensing instrumentation, including eight 915-MHz wind-profiling radars, and surface meteorological stations. The evolution and annual variability of the daytime convective planetary-boundary-layer (PBL) height is investigated using the wind-profiling radars. Three models with different horizontal grid spacing are evaluated: the Rapid Refresh, the High-Resolution Rapid Refresh, and its nested version. The results are used to assess errors in the prediction of PBL height within the experimental and control versions of the models, with the experimental versions including changes and additions to the model parametrizations developed during the field campaign, and the control version using the parametrizations present in the National Oceanic and Atmospheric Administration/National Centers for Environmental Prediction operational version of the models at the start of the project. Results show that the high-resolution models outperform the low-resolution versions, the experimental versions perform better compared with the control versions, model PBL height estimations are more accurate on cloud-free days, and model estimates of the PBL height growth rate are more accurate than model estimates of the rate of decay. Finally, using surface sensors, we assess surface meteorological variables, finding improved surface irradiance and, to a lesser extent, improved 2-m temperature in the experimental version of the model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Climate Change Impacts on the South American Monsoon System and Its Surface–Atmosphere Processes Through RegCM4 CORDEX-CORE Projections

Here, we evaluate projected changes in surface water and energy balances and surface–atmosphere coupling in the South American Monsoon System (SAMS) for the end of the century (2080–2099). The analyses are based on two ensemble datasets, which follow Representative Concentration Pathway 8.5 in the future period, and cover four subdomains (Northern and Southern—NAMZ and SAMZ—Amazon, La Plata Basin—LPB, and Southern Southeast Brazil—SSB). One ensemble consists of three Global Climate Models (HadGEM2-ES, MPI-ESM-MR and NorESM1-M), while the other consists of their dynamically downscaled version at 25 km horizontal grid spacing using Regional Climate Model version 4 (RegCM4). As both ensembles are able in reproducing the annual cycle of the components of the surface water and energy balances in the present climate, they can be used in the study of future climate. During the wet season (November–March), both ensembles project a decrease in precipitation over NAMZ and SAMZ (an exception is RegCM4 that projects a slight increase in SAMZ), and an increase across the LPB and SSB. These changes do not cause retreat or expansion of the monsoon area over the continent, which is similar to the present climate (1995–2014). For the wet season, the ensembles are in line with the presence of a strong surface–atmosphere coupling in LPB and SSB, weak coupling in SAMZ and very weak coupling in NAMZ. For future climate, the coupling is even weaker in NAMZ, which may be a driver for the negative changes in precipitation. For the other subdomains, while the ensembles project similar signals of precipitation changes, they disagree with the surface-atmosphere coupling highlighting the uncertainties in future climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Multi-species collisions for delta-f gyrokinetic simulations: Implementation and verification with GENE

we report a multi-species linearized collision operator based on the model developed by Sugama et al. has been implemented in the nonlinear gyrokinetic code, GENE. Such a model conserves particles, momentum, and energy to machine precision, and is shown to have negative definite free energy dissipation characteristics, satisfying Boltzmann’s H-theorem, including for realistic mass ratio. Finite Larmor Radius (FLR) effects have also been implemented into the local version of the code. For the global version of the code, the collision operator has been developed to allow for block-structured velocity space grids, allowing for computationally tractable collisional global simulations. The validity of the collision operator has been demonstrated by relaxation and conservation tests, as well as appropriate benchmarks. The newly implemented operator shall be used in future simulations to study magnetically confined fusion plasma turbulence and transport in more extreme regions with higher collisionality.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A projection method for particle resampling

Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase-space due to their better scalability than continuum approaches with respect to dimension. Complex processes collectively referred to as particle noise hamper long time simulations with particle methods. One approach to address this problem is particle mesh adaptivity, or remapping, known as particle resampling and remeshing. Here, this work introduces a resampling method that projects particles to and from a (finite element) function space. The method is simple, using standard sparse linear algebra and finite element techniques, and it preserves all moments up to the order of a polynomial represented exactly by the continuum function space. It is distinguished from most other mesh-based methods in that new particle positions and number are decoupled from the mesh, allowing particle and continuum meshes to be adapted relatively independently. While this work is developed with structured particle and continuum phase-space grids on 1X + 1V Vlasov-Poisson models of Landau damping and two-stream instability, the method is well-suited to unstructured grids. Stable long time dynamics are demonstrated up to time T = 500. Reproducibility artifacts and data are publicly available.

Kinetic methods↗
DOE OSTI2020

Three-dimensional seismic characterization and imaging of the Soda Lake geothermal field

Accurate characterization of subsurface geophysical properties and detection of the fault system are essential for geothermal energy exploration and production. The Soda Lake geothermal field is in western Nevada with a complex fault system. Previous seismic characterization only produced a low-resolution, smooth velocity model along with a simple, conceptual fault model. Using optimized correlation-based full-waveform inversion, wavefield-separation-based reverse-time migration, and automatic fault detection techniques, we present 3D seismic characterization for the Soda Lake geothermal field using 3D surface seismic data acquired with Vibroseis sources. Here, we obtain 3D high-resolution velocity, density, and acoustic impedance models, 3D seismic images with different grid spacings, and a high-resolution fault system. Consistency check between the constructed faults and currently active injection and production geothermal wells verifies that our seismic inversion and imaging results and detected faults are reliable. These results can provide valuable information for optimizing well placement and geothermal energy production at the Soda Lake geothermal field.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗
DOE OSTI2024

Influence of alloy solidification path on melt pool behavior in additive manufacturing

Numerical models used to study transport phenomena in laser-powder bed fusion processes often rely on assumptions and simplifications to reduce their computational expense. One common simplification is in the description of latent heat evolution during the solid-liquid phase change (i.e., the solidification pathway), justified by the fact that the mushy zone thickness is similar to the numerical grid spacing used for continuum transport models. The lack of resolution of transport phenomena in the mushy zone motivates the use of computationally convenient solidification paths such as linear or sigmoidal relationships over pathways derived from fundamental solidification theory such as equilibrium or Scheil models. In the present work, an uncertainty quantification (UQ) framework is used to analyze the influence of solidification pathway selection on the solidification dynamics and melt pool geometries in laser based additive manufacturing (AM) of IN625. Results show the solidification pathway has a quantifiable influence on the cooling rate at the liquidus isotherm, mushy zone thickness, and solidification time. Due to similarities in the latent heat evolution at the beginning of solidification, the equilibrium and Scheil models predict similar cooling rates near the liquidus isotherm, however the wider freezing range of Scheil leads to a wider mushy zone compared to equilibrium. The non-physical latent heat release profiles of sigmoidal and linear paths lead to significant overpredictions of cooling rates at the liquidus isotherm compared to equilibrium and Scheil. Finally, these results indicate that careful consideration should be given to the choice of solidification pathway to ensure reliable model predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system

The development of an implicit, unconditionally stable, numerical method for solving the Vlasov–Poisson system in one dimension using a phase-space grid is presented. The algorithm uses the Crank–Nicolson discretization scheme and operator splitting allowing for direct solution of the finite difference equations. This method exactly conserves particle number, enstrophy and momentum. A variant of the algorithm which does not use splitting also exactly conserves energy but requires the use of iterative solvers. This algorithm has no dissipation and thus fine-scale variations can lead to oscillations and the production of negative values of the distribution function. We find that overall, the effects of negative values of the distribution function are relatively benign. We consider a variety of test cases that have been used extensively in the literature where numerical results can be compared with analytical solutions or growth rates. We examine higher-order differencing and construct higher-order temporal updates using standard composition methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Changes in Precipitation From North Atlantic Tropical Cyclones Under RCP Scenarios in the Variable-Resolution Community Atmosphere Model

Decreasing climate models' grid spacing improves the representation of tropical cyclones at decadal time scales. In this study, a variable-resolution (VR) version of the Community Atmosphere Model 5 (CAM5-VR) is utilized to study North Atlantic tropical cyclone climatology in ensemble historical climate simulations and under two Representative Concentration Pathway (RCP) projections (RCP4.5 and RCP8.5). Basin-wide tropical cyclone counts decrease in the RCP simulations, although landfalling storm counts do not show as straightforward of a pattern, especially when focusing on regional changes. Lifetime maximum intensity metrics suggest that tropical cyclones increase in strength in the RCP ensembles. However, despite increases in tropical cyclone-related precipitation rates and the amount of precipitation produced per storm with warming, the annual average Rx5day from tropical cyclones over the eastern United States decreases due to less landfalling storms. This work is part of a continued effort to quantify how tropical cyclone-induced hazards may change in future climates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Characterizing Subsiding Shells in Shallow Cumulus Using Doppler Lidar and Large‐Eddy Simulation

Abstract The existence of subsiding shells on the periphery of shallow cumulus clouds has major implications concerning the parameterization of shallow convection, with the mass exchange between the shell and cloudy air representing a significant deviation from the commonly used bulk‐plume parameterization. We examine the structure and frequency of subsiding shells in shallow cumulus convection using Doppler lidars at the Atmospheric Radiation Measurement Southern Great Plains facility in the central United States and at the Jülich ObservatorY for Cloud Evolution in western Germany. Doppler lidar indicates that the vertical subsiding shell extent is asymmetric, while shell width is typically ~100 m. Large‐eddy simulation can reasonably simulate the observed shell structure using a grid spacing of 10 m and suggests that much of the observed asymmetry is not a result of transient cloud evolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Characterizing uncertainties in ground “truth” of precipitation over complex terrain through high-resolution numerical modeling

Ground observation of precipitation over complex terrain is subject to large uncertainties due to inadequate sampling. This study explores a method that combines limited gauge data and a high-resolution numerical simulation to quantify the precipitation uncertainties in central Himalaya. Specifically, the Coefficient of spatial Variability (CV) of precipitation and the minimum Number of Required Stations (NRS) to obtain areal-mean precipitation ground truth values within a 0.25-degree area are investigated using fine-scale meteorological simulation at 1.5 km grid spacing. Evaluation over a densely-gauged area demonstrates comparable CV and NRS values between station observations and simulations. The simulated CV and NRS values in a larger area show a strong and positive dependence on each other and an expected positive (negative) correlation with topographic complexity (temporal scale). The proposed method sheds lights on evaluating precipitation products and holds promise for informing the layout of rain gauge networks in complex terrain.

Ouyang, Lin↗
DOE OSTI2021

Impact of Lateral Flow on Surface Water and Energy Budgets over the Southern Great Plains – A Modeling Study

As we see the horizontal grid spacing decrease, treatment of hydrologic processes in land surface models, such as the lateral flow of surface and subsurface flow, need to be explicitly represented. Unlike previous studies that mainly focused on the mountainous regions, in this study the offline WRF-Hydro model is employed to study the impact of lateral flow on soil moisture and energy fluxes over the relatively flat southern Great Plains (SGP). The vast amount of measurements over the SGP provide an unique opportunity to assess the model behavior. In addition, newly developed land surface properties and input forcing are ingested into the model, in an attempt to reduce uncertainties associated with the initial and boundary forcing and help to identify model deficiencies. Our results show that the more realistic inputs (parameters, soil types, forcing) lead to larger underestimation of latent heat flux and dry bias, indicating the existence of model structural uncertainty (embedded errors) in WRF-Hydro that need to be characterized to inform future model development efforts. Including lateral flow processes partly mitigates the model deficiencies in representing hydrologic processes and alleviates the dry bias. In particular, both surface and subsurface lateral flow increase soil moisture mainly over the lower elevations, except that subsurface flow also affects soil moisture over steeper terrains. Additional simulations are performed to assess the effect of routing resolution on model results. When LSM resolution is high, noticeable differences in soil moisture are produced between different routing resolutions especially over steep terrain. Whereas when LSM resolution is coarse, differences between routing resolutions become negligible, especially over flat terrain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Convective-To-Total Precipitation Ratio and the “Drizzling” Bias in Climate Models

Overestimation of precipitation frequency and duration while underestimating intensity, that is, the “drizzling” bias, has been a long-standing problem of global climate models. Here we explore this issue from the perspective of precipitation partitioning. We found that most models in the Climate Model Intercomparison Project Phase 5 (CMIP5) have high convective-to-total precipitation (PC/PR) ratios in low latitudes. Convective precipitation has higher frequency and longer duration but lower intensity than non-convective precipitation in many models. As a result, the high PC/PR ratio contributes to the “drizzling” bias over low latitudes. The PC/PR ratio and associated “drizzling” bias increase as model resolution coarsens from 0.5° to 2.0°, but the resolution's effect weakens as the grid spacing increases from 2.0° to 3.0°. Some of the CMIP6 models show reduced “drizzling” bias associated with decreased PC/PR ratio. Thus, more reasonable precipitation partitioning, along with finer model resolution should alleviate the “drizzling” bias within current climate models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗