Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid 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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Relating Convection to GCM Grid-Scale Fields Using Cloud-Resolving Model Simulation of a Squall Line Observed during MC3E Field Experiment

In this study, a WRF (Weather Research and Forecasting) model is used as a cloud-resolving model to simulate a squall line observed on 20 May 2011 in the Southern Great Plains (SGP) of the United States. The model output is then used to examine the relationships between convective precipitation and coarse-grained variables averaged over a range of subdomain sizes equivalent to various global climate model horizontal resolutions. The objective is to determine to what extent convection within the subdomains can be related to these “large-scale” variables, thus that they can potentially serve as closure in convective parameterization. Results show that convective precipitation is well correlated with the vertical velocity at 500 hPa, column integrated moisture convergence and CAPE change due to large-scale advective forcing (dCAPE) for various subdomain sizes, but the correlation decreases with decreasing subdomain size. dCAPE leads convective precipitation for all subdomain sizes examined; however, the lead time decreases with decreasing subdomain size. Moisture convergence leads convective precipitation for subdomain sizes greater than 32 km but has no lead time for smaller subdomain sizes. Mid-tropospheric vertical velocity has no lead time or slightly lags convective precipitation. The lead/lag composite analysis with respect to maximum precipitation time indicates that peaks of large-scale variables increase with decreasing subdomain size. The peaks of 500 hPa vertical velocity and column integrated moisture convergence occur at the same time as maximum precipitation, but maximum dCAPE leads maximum precipitation by twelve minutes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Counterrotating prop-fan simulations which feature a relative-motion multiblock grid decomposition enabling arbitrary time-steps

Improvements are presented of a computer algorithm developed for the time-accurate flow analysis of rotating machines. The flow model is a finite volume method utilizing a high-resolution approximate Riemann solver for interface flux definitions. The numerical scheme is a block LU implicit iterative-refinement method which possesses apparent unconditional stability. Multiblock composite gridding is used to orderly partition the field into a specified arrangement of blocks exhibiting varying degrees of similarity. Block-block relative motion is achieved using local grid distortion to reduce grid skewness and accommodate arbitrary time step selection. A general high-order numerical scheme is applied to satisfy the geometric conservation law. An even-blade-count counterrotating unducted fan configuration is chosen for a computational study comparing solutions resulting from altering parameters such as time step size and iteration count. The solutions are compared with measured data.

Janus, J. Mark↗
NASA NTRSDate not supplied

An Overset Grid Motion Driver for Resolving Entry Vehicle CFD Simulations with Large Changes in Wake Orientation

Time-accurate computational fluid dynamics (CFD) simulations using overset meshes can enable reduced computational cost compared to a rigid, monolithic mesh for blunt-body atmospheric entry vehicles. If the vehicle has a significant change in attitude, the high-density region of the overset mesh designed for resolving the subsonic wake can be kept in position as the wake orientation changes relative to the body. Current simulations generally keep the wake mesh fixed in attitude while the body rotates, assuming the wake orientation remains relatively constant. However, changes in the relative velocity over a long-duration simulation or effects of the body geometry can invalidate this assumption. This work presents a method for calculating the desired mesh orientation to resolve a blunt body wake using the atmosphere-relative velocity of the vehicle. A motion driver using these calculations is implemented in a CFD-in-the-loop flight model bridging the FUN3D flow solver and POST2 trajectory propagator. Overset meshing is performed using the Yoga extension for FUN3D. Results are presented for simulations with both overset and monolithic meshes under large changes in relative velocity (and therefore wake orientation). The results demonstrate that the motion driver successfully orients an overset wake mesh to capture the subsonic wake with sufficient resolution despite variable body rotation and body-relative freestream velocity.

CFD↗
NASA NTRSDate not supplied

An Overset Grid Motion Driver for Resolving Entry Vehicle CFD Simulations with Large Changes in Wake Orientation

Time-accurate computational fluid dynamics (CFD) simulations using overset meshes can enable reduced computational cost compared to a rigid, monolithic mesh for blunt-body atmospheric entry vehicles. If the vehicle has a significant change in attitude, the high-density region of the overset mesh designed for resolving the subsonic wake can be kept in position as the wake orientation changes relative to the body. Current simulations generally keep the wake mesh fixed in attitude while the body rotates, assuming the wake orientation remains relatively constant. However, changes in the relative velocity over a long-duration simulation or effects of the body geometry can invalidate this assumption. This work presents a method for calculating the desired mesh orientation to resolve a blunt body wake using the atmosphere-relative velocity of the vehicle. A motion driver using these calculations is implemented in a CFD-in-the-loop flight model bridging the FUN3D flow solver and POST2 trajectory propagator. Overset meshing is performed using the Yoga extension for FUN3D. Results are presented for simulations with both overset and monolithic meshes under large changes in relative velocity (and therefore wake orientation). The results demonstrate that the motion driver successfully orients an overset wake mesh to capture the subsonic wake with sufficient resolution despite variable body rotation and body-relative freestream velocity.

Computational Fluid Dynamics↗
DOE OSTI2020

Impacts of Stochastic Mixing in Idealized Convection-Permitting Simulations of Squall Lines

This study investigates impacts of altering subgrid-scale mixing in “convection-permitting” kilometer-scale horizontal-grid-spacing (Δ h ) simulations by applying either constant or stochastic multiplicative factors to the horizontal mixing coefficients within the Weather Research and Forecasting Model. In quasi-idealized 1-km Δ h simulations of two observationally based squall-line cases, constant enhanced mixing produces larger updraft cores that are more dilute at upper levels, weakens the cold pool, rear-inflow jet, and front-to-rear flow of the squall line, and degrades the model’s effective resolution. Reducing mixing by a constant multiplicative factor has the opposite effect on all metrics. Completely turning off parameterized horizontal mixing produces bulk updraft statistics and squall-line mesoscale structure closest to an LES “benchmark” among all 1-km simulations, although the updraft cores are too undilute. The stochastic mixing scheme, which applies a multiplicative factor to the mixing coefficients that varies stochastically in time and space, is employed at 0.5-, 1-, and 2-km Δ h . It generally reduces midlevel vertical velocities and enhances upper-level vertical velocities compared to simulations using the standard mixing scheme, with more substantial impacts at 1- and 2-km Δ h compared to 0.5-km Δ h . Further, the stochastic scheme also increases updraft dilution to better agree with the LES for one case, but has less impact on the other case. Stochastic mixing acts to weaken the cold pool but without a significant impact on squall-line propagation. It also does not affect the model’s overall effective resolution unlike applying constant multiplicative factors to the mixing coefficients.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The June 2012 North American Derecho: A Testbed for Evaluating Regional and Global Climate Modeling Systems at Cloud‐Resolving Scales

Abstract In this paper, we introduce a testbed for evaluating and comparing climate modeling systems at cloud resolving scales using hindcasts of the June 2012 North American derecho. To demonstrate its utility for model intercomparison, the testbed is applied to two models: the regionally‐refined Simple Cloud‐Resolving E3SM Atmosphere Model (SCREAM) at 6.5, 3.25 and 1.625 km grid spacing and the Weather Research and Forecasting (WRF) model with 3.2 and 1.6 km grid spacing. We find the simulation results to be highly sensitive to the initial conditions (ICs), initialization time, and model configurations, with ICs from the Rapid Refresh producing the best simulation. Significant improvement is identified in both models as horizontal grid spacing is refined. While a propagation delay of approximately 2 hr is found in both models, SCREAM at 1.625 km simulates the observed bow echo structure of the derecho well and predicts strong surface gusts that exceed 30 m/s. In comparison, WRF has difficulty producing surface wind over 25 m/s, with wind gusts in WRF 42%–46% lower than in SCREAM. However, WRF has a lower bias in simulating cloud top temperature and extent, but overestimates precipitation intensity. Both models reproduce the observed outgoing longwave radiation spatial patterns well (Pearson correlation >0.88), but, compared with NEXRAD observations, simulate generally larger areas of composite radar reflectivity >40 dBZ and underestimate the precipitating area by ∼47%.

2012 North American derecho↗
NASA NTRS1995

Navier-Stokes Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP

The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time, load balancing, and scalability will be provided.

Vadyak, Joseph↗
DOE OSTI2026

Recent progress on coarse graining simulations

We focus on coarse graining simulations based on the primary conservation equations, effectively codesigned physics and algorithms, and low-Mach-number corrected (LMC) hydrodynamics. Simulation methods involve LANL’s x-Radiation-Adaptive-Grid-Eulerian Large-Eddy Simulation, Besnard-Harlow-Rauenzahn (BHR) Reynolds-Averaged Navier-Stokes (RANS) approach, and Dynamic BHR – a paradigm bridging RANS and LES. A relevant question addressed relates to whether 3D RANS and RANS/LES hybrids – the industry standards for aerospace and automotive research, are presently relevant for practical variable-density applications involving shocked and accelerated interface instabilities. Furthermore, recent simulations of the GaTECH inclined mixing-layer shock-tube and NIF ICF-capsule experiments are used to demonstrate issues, challenges, and potential for 3D coarse grained LMC simulation strategies for robustly simulating complex transitional and coupled hydrodynamics-multiphysics with coarser resolution. Present LES readiness to provide accurate predictions at scale is demonstrated – whereas 3D RANS and RANS/LES bridging do not appear impactful in this context.

42 ENGINEERING↗
NASA NTRSDate not supplied

Simulating Transonic Buffet Aerodynamics for the Boeing Transonic Truss-Braced Wing Aircraft

The flow past the Boeing Transonic Truss-Braced Wing (TTBW) aircraft is simulated with a Hybrid Reynolds-Averaged Navier Stokes Large-Eddy Simulations (HRLES) turbulence modeling approach to study transonic buffet onset and evaluate the predictive capability of the numerical approach. All simulations used structured overset curvilinear grids for the free-air configuration, with all simulations performed using the Launch, Ascent, and Vehicle Aerodynamics (LAVA) computational framework. Aerodynamic loads and surface pressure obtained from the HRLES solutions are compared to results from unsteady RANS (URANS) and experiments conducted at NASA Ames 11-by 11-Foot Transonic Wind Tunnel. Overall good agreement is obtained in the predicted loads and surface pressure. Unsteady pressure data from HRLES and URANS simulations is utilized to compute power spectral density (PSD) to predict transonic buffet onset and compared to the Kulite data obtained from the experiment. Comparisons of the PSD spectrum also show reasonable agreement with the wind tunnel data, with HRLES and URANS predicting buffet onset at a slightly earlier angle-of-attack than the experiment.

AATT↗
NASA NTRS1991

Simulation Of Flow In A Turbine Cascade

Report discusses application of lower-upper (LU) implicit numerical-integration scheme, together with interactive grid-generation method, to simulation of inviscid, compressible flow in two-dimensional model of turbine cascade. Study part of continuing effort to develop techniques for computation of complicated three-dimensional flows in turbomachinery.

Choo, Yung K.↗
NASA NTRS2023

Optimizing the Resolution of Hydrodynamic Simulations for MCRaT Radiative Transfer Calculations

Despite their discovery about half a century ago, the Gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and techniques. One such tool is the PLUTO hydro dynamics code, which is used to numerically simulate GRB outflows. PLUTO uses Adaptive Mesh Refinement to focus computational efforts on the portion of the grid that contains the simulated jet. Another tool is the Monte Carlo Radiation Transfer (MCRaT) code, which predicts electromagnetic signatures of GRBs by conducting photon scatterings within a jet using PLUTO. The effects of the underlying resolution of a PLUTO simulation with respect to MCRaT post-processing radiative transfer results have not yet been quantified. We analyze an analytic spherical outflow and a hydrodynamically simulated GRB jet with MCRaT at varying spatial and temporal resolutions and quantify how decreasing both resolutions affect the resulting mock observations. We find that changing the spatial resolution changes the hydrodynamic properties of the jet, which directly affect the MCRaT mock observable peak energies. We also find that decreasing the temporal resolution artificially decreases the high energy slope of the mock observed spectrum, which increases both the spectral peak energy and the luminosity. We show that the effects are additive when both spatial and temporal resolutions are modified. Our results allow us to understand how decreased hydrodynamic temporal and spatial resolutions affect the results of post-processing radiative transfer calculations, allowing for the optimization of hydrodynamic simulations for radiative transfer codes.

Jose Arita-Escalante↗
NASA NTRS1987

Analysis of subgrid-scale eddy viscosity with use of results from direct numerical simulations

Without resort to any modeling, subgrid-scale eddy viscosity is computed from the results of high-resolution (64 cubed and 128 cubed grid points) direct numerical simulations of three-dimensional homogeneous isotropic decaying turbulence. In these simulations the eddy viscosity peaks sharply at the cutoff wavenumber, in rough agreement with the results of Kraichnan (1976). In addition, in the low-wavenumber range the eddy viscosity may be negative, contrary to the generally accepted concept of subgrid-scale eddy viscosity. Some possible explanations of this behavior are discussed.

Domaradzki, J. Andrzej↗
DOE OSTI2023

Analyzing at-scale distribution grid response to extreme temperatures

Threats against power grids continue to increase, as extreme weather conditions and natural disasters (extreme events) become more frequent. Hence, there is a need for the simulation and modeling of power grids to reflect realistic conditions during extreme events conditions, especially distribution systems. Herein, this paper presents a modeling and simulation platform for electric distribution grids which can estimate overall power demand during extreme weather conditions. The presented platform's efficacy is shown by demonstrating estimation of electrical demand for 1) Electricity Reliability Council of Texas (ERCOT) during winter storm Uri in 2021, and 2) alternative hypothetical scenarios of integrating Distributed Energy Resources (DERs), weatherization, and load electrification. In comparing to the actual demand served by ERCOT during the winter storm Uri of 2021, the proposed platform estimates approximately 34 GW of peak capacity deficit. These numbers are consistent with state-of-the-art prediction results published in the literature. For the case of the future electrification of heating loads, peak capacity of 78 GW (124% increase) is estimated, which would be reduced to 47 GW (38% increase) with the adoption of efficient heating appliances and improved thermal insulation. Integrating distributed solar PV and storage into the grid causes improvement in the local energy utilization and hence reduces the potential unmet energy by 31% and 40%, respectively.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Evaluating Large-Eddy Simulations of Boundary-Layer Cloud Response to Climate Change using MAGIC Observations

Over the cool waters off California, a persistent and extensive deck of low-lying clouds forms atop a shallow, turbulent layer of moist air extending 500-2000 meters above the ocean, above which the air is much warmer and drier. Its west edge varies with weather and season and is typically half way to Hawaii. This cloud deck and similar cloud decks elsewhere reflect sunlight back to space, cooling the underlying ocean and indeed the entire planet. They must be accurately represented by a climate model to correctly simulate our current climate or future global warming. Most climate models suggest these cloud decks will thin and break up in a warmer CO 2 -enriched climate, so extra sunlight would reach the ocean surface and amplify warming worldwide. But these models have difficulty accurately simulating the cloud thickness and extent because the turbulent cloud-forming air currents are much too small to be simulated by their coarse computational grid. In 2012-2013, the DOE Atmospheric Radiation Measurement MAGIC program instrumented a container ship transiting between California and Hawaii to observe clouds and rain, small particles (aerosols), and weather conditions. Data were obtained from 20 4-6 day cruises in each direction. This study tested whether a sophisticated ‘large-eddy simulation’ (LES) of the turbulent air motions that create the clouds in this region accurately predicts the variations in cloud properties seen across the different cruises. If so, that model is more believable for predicting cloud changes in this region in a future warmer climate. Our LES used a numerical mesh over a 6.4x6.4 km region with 50 m horizontal and 5 m vertical grid spacing to simulate the cloud-forming eddies and the sharp overlying transition between moist and dry air. To better compare with MAGIC observations, we configured our model to simulate a region that moved along with the ship. We used a weather forecast model to specify effects of the large-scale air flow. We compared our results with balloon ascents, a radar and other instruments pointed upward at the clouds, measurements of downward sunlight and infrared radiation, aerosols, and weather variables for all suitable cruise periods. Our LES was remarkably skillful in predicting the cloudiness variations measured by the ship. MAGIC also sampled mesoscale cells, regions of cloud thickening and thinning over distances of 10-50 km. An LES run over a larger region 100 km on a side shows these cells form spontaneously, enhanced by strong radiative heat loss at the cloud tops.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Quantifying the uncertainty introduced by discretization and time-averaging in two-fluid model predictions

The two-fluid model (TFM) has become a tool for the design and troubleshooting of industrial fluidized bed reactors. To use TFM for scale up with confidence, the uncertainty in its predictions must be quantified. Here, we study two sources of uncertainty: discretization and time-averaging. First, we show that successive grid refinement may not yield grid-independent transient quantities, including cross-section–averaged quantities. Successive grid refinement would yield grid-independent time-averaged quantities on sufficiently fine grids. A Richardson extrapolation can then be used to estimate the discretization error, and the grid convergence index gives an estimate of the uncertainty. Richardson extrapolation may not work for industrial-scale simulations that use coarse grids. We present an alternative method for coarse grids and assess its ability to estimate the discretization error. Second, we assess two methods (autocorrelation and binning) and find that the autocorrelation method is more reliable for estimating the uncertainty introduced by time-averaging TFM data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2021

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↗
NASA NTRS2016

Similar Estimates of Temperature Impacts on Global Wheat Yield by Three Independent Methods

The potential impact of global temperature change on global crop yield has recently been assessed with different methods. Here we show that grid-based and point-based simulations and statistical regressions (from historic records), without deliberate adaptation or CO2 fertilization effects, produce similar estimates of temperature impact on wheat yields at global and national scales. With a 1 C global temperature increase, global wheat yield is projected to decline between 4.1% and 6.4%. Projected relative temperature impacts from different methods were similar for major wheat-producing countries China, India, USA and France, but less so for Russia. Point-based and grid-based simulations, and to some extent the statistical regressions, were consistent in projecting that warmer regions are likely to suffer more yield loss with increasing temperature than cooler regions. By forming a multi-method ensemble, it was possible to quantify 'method uncertainty' in addition to model uncertainty. This significantly improves confidence in estimates of climate impacts on global food security.

Climate impacts↗