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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Toward the large-eddy simulation of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96(exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.

Erlebacher, G.↗
NASA NTRS1992

Toward the large-eddy simulation of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 (exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine-grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.

Erlebacher, G.↗
DOE OSTI2022

Optimizing detachment control using the magnetic configuration of divertors

Abstract As tokamak research moves to reactor conditions, the control of a stable, optimally-detached divertor plasma has become increasingly relevant. Simple predictions of such detachment control have been performed previously using the detachment location sensitivity (DLS) model (Lipschultz et al 2016 Nucl. Fusion 56 056007). In this study the DLS model is extended and combined with SOLPS-ITER simulations of isolated divertor grids to study the effects of alternate divertor magnetic field properties on detachment control. The DLS model predicts that divertors can achieve easier access to detachment through a long connection length, a high total flux expansion, and a high average magnetic field in the divertor compared with that at the x -point. SOLPS-ITER simulations suggest an even stronger impact of total flux expansion and connection length on detachment access than the DLS model predicts. In terms of detachment evolution, both simulation and modelling show a high gradient in the total magnetic field and low B p o l B are able to more easily keep a detachment front at a desired poloidal location. Regions with no stable detachment front locations can arise for high magnetic field gradients directed towards the target. Significant differences have been found between impurity scan simulations and the DLS model. These differences may be attributed to sources and sinks of power and electron pressure, and should be explored further in future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

CFD Analysis of RLUOB Zone 1 HEPA Filter Plenum and Testing Manifolds

The purpose of this computational fluid dynamics (CFD) analysis is to ensure that Camfil Farr’s upstream and downstream injection and sampling manifolds can meet or exceed the requirements outlined in the ASME AG-1 1997 a(2000) Code on Nuclear Air and Gas Treatment for the testing of HEPA and adsorbent filters. This paper will present a numerical simulation of airflow in the Radiological Laboratory Utility Office Building (RLUOB) zone 1 HEPA filter plenum and testing manifolds using the commercial CFD software ANSYS FLUENT 2020R1. The CFD analysis focuses on the investigation of the air flow distribution and air-aerosol mixing uniformity. The evaluation was done for all steps of the modeling process: grid generation, physics setup, simulation, and post-processing. The mass flow rate in each section of the zone 1 injection and sampling manifolds is also reported.

42 ENGINEERING↗
DOE OSTI2015

CHOLLA: A NEW MASSIVELY PARALLEL HYDRODYNAMICS CODE FOR ASTROPHYSICAL SIMULATION

We present Computational Hydrodynamics On ParaLLel Architectures (Cholla ), a new three-dimensional hydrodynamics code that harnesses the power of graphics processing units (GPUs) to accelerate astrophysical simulations. Cholla models the Euler equations on a static mesh using state-of-the-art techniques, including the unsplit Corner Transport Upwind algorithm, a variety of exact and approximate Riemann solvers, and multiple spatial reconstruction techniques including the piecewise parabolic method (PPM). Using GPUs, Cholla evolves the fluid properties of thousands of cells simultaneously and can update over 10 million cells per GPU-second while using an exact Riemann solver and PPM reconstruction. Owing to the massively parallel architecture of GPUs and the design of the Cholla code, astrophysical simulations with physically interesting grid resolutions (≳256{sup 3}) can easily be computed on a single device. We use the Message Passing Interface library to extend calculations onto multiple devices and demonstrate nearly ideal scaling beyond 64 GPUs. A suite of test problems highlights the physical accuracy of our modeling and provides a useful comparison to other codes. We then use Cholla to simulate the interaction of a shock wave with a gas cloud in the interstellar medium, showing that the evolution of the cloud is highly dependent on its density structure. We reconcile the computed mixing time of a turbulent cloud with a realistic density distribution destroyed by a strong shock with the existing analytic theory for spherical cloud destruction by describing the system in terms of its median gas density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

From plot to scale: ex-ante assessment of conservation agriculture in Zambia

This study combined bottom-up and top-down approaches to assess the ex-ante effects of conservation agriculture (CA)-based systems in Zambia considering both biophysical and economic factors and prevailing farm systems characteristics. For continuous maize cropping we compared a CA-based system of no-tillage with crop residue retention to a control system of conventional tillage with crop residue removal. First, we simulated yield effects that were calibrated and evaluated against multiple datasets, including on-farm agronomic trials from two seasons and six sites. Next, we extrapolated our simulations to all maize-growing areas in Zambia using gridded climate and soil datasets. Then simulated yields (in kg ha( -1 )) were combined with economic data from a nationally-representative household survey to construct economic indicators including benefit-cost ratios (based on gross benefits and variable costs both in $ ha( -1 )) that captured the implicit value of crop residues and labor demands. The field scale (per ha) indicators were scaled out using harvested areas as an expansion factor. All indicators were calculated over 3-, 10-, and 20-year simulation periods using an interpolated sequence of historical climate data. Finally, we conducted a spatial farm typology analysis to help understand the spatial variation in our field-scale indicators and provide insights into trade-offs and the suitability of CA-based systems for farmers. Average changes in yield from using CA-based systems (compared with the control) at the district scale ranged from -37% to 70% (average 33%), with a similar range of changes in benefit-cost ratios once economic factors were included, in addition to intra-district yield variability. Combining the changes in benefit cost ratios with maize harvested area resulted in an average annual change in district-scale net benefit ranging from US $ - 3.9 to US $9.9 million (with an average of US $1.1 million). The heterogeneity in biophysical and economic factors gave a ranking of provinces different according to biophysical or economic indicators, reinforcing the importance of coupling biophysical and economic approaches. The spatial farm typology analysis highlighted the specific contexts of farmers relevant to the suitability of CA, such as their mineral fertilizer applications rates, ownership of livestock, and prevailing soil texture and rainfall.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2004

HiPEP Ion Optics System Evaluation Using Gridlets

Experimental measurements are presented for sub-scale ion optics systems comprised of 7 and 19 aperture pairs with geometrical features that are similar to the HiPEP ion optics system. Effects of hole diameter and grid-to-grid spacing are presented as functions of applied voltage and beamlet current. Recommendations are made for the beamlet current range where the ion optics system can be safely operated without experiencing direct impingement of high energy ions on the accelerator grid surface. Measurements are also presented of the accelerator grid voltage where beam plasma electrons backstream through the ion optics system. Results of numerical simulations obtained with the ffx code are compared to both the impingement limit and backstreaming measurements. An emphasis is placed on identifying differences between measurements and simulation predictions to highlight areas where more research is needed. Relatively large effects are observed in simulations when the discharge chamber plasma properties and ion optics geometry are varied. Parameters investigated using simulations include the applied voltages, grid spacing, hole-to-hole spacing, doubles-to-singles ratio, plasma potential, and electron temperature; and estimates are provided for the sensitivity of impingement limits on these parameters.

Willliams, John D.↗
DOE OSTI2020

A Model-Predictive Hierarchical-Control Framework for Aggregating Residential DERs to Provide Grid Regulation Services

This paper develops a hierarchical control frame-work to aggregate and to manage behind-the-meter distributed energy resources (DERs), which will be ubiquitous in future distribution systems. In the proposed framework, firstly, each controller in the hierarchy determines the flexibility of the DERs such that the obtained flexibility is feasible with respect to its operational purview. For example, the operational purview of a home energy management system may only consider consumer comfort preferences, while that for an aggregator or a grid controller may consider network voltage management as well. Based on the feasible flexibility, optimal setpoints for the DERs is then determined by the hierarchical controllers to help the distribution power network in voltage regulation, coordination issues with existing transmission-level conventional generators, etc. Therefore, the proposed strategy, which is based on model-predictive control, can be effectively utilized by the distribution network to coordinate several DERs to provide grid regulation services. Numerical simulations performed on the IEEE 37-bus test system demonstrate the efficacy of the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Projecting Future Energy Production from Operating Wind Farms in North America. Part I: Dynamical Downscaling

Abstract New simulations at 12-km grid spacing with the Weather and Research Forecasting (WRF) Model nested in the MPI Earth System Model (ESM) are used to quantify possible changes in wind power generation potential as a result of global warming. Annual capacity factors (CF; measures of electrical power production) computed by applying a power curve to hourly wind speeds at wind turbine hub height from this simulation are also used to illustrate the pitfalls in seeking to infer changes in wind power generation directly from low-spatial-resolution and time-averaged ESM output. WRF-derived CF are evaluated using observed daily CF from operating wind farms. The spatial correlation coefficient between modeled and observed mean CF is 0.65, and the root-mean-square error is 5.4 percentage points. Output from the MPI-WRF Model chain also captures some of the seasonal variability and the probability distribution of daily CF at operating wind farms. Projections of mean annual CF (CF A ) indicate no change to 2050 in the southern Great Plains and Northeast. Interannual variability of CF A increases in the Midwest, and CF A declines by up to 2 percentage points in the northern Great Plains. The probability of wind droughts (extended periods with anomalously low production) and wind bonus periods (high production) remains unchanged over most of the eastern United States. The probability of wind bonus periods exhibits some evidence of higher values over the Midwest in the 2040s, whereas the converse is true over the northern Great Plains. Significance Statement Wind energy is playing an increasingly important role in low-carbon-emission electricity generation. It is a “weather dependent” renewable energy source, and thus changes in the global atmosphere may cause changes in regional wind power production (PP) potential. We use PP data from operating wind farms to demonstrate that regional simulations exhibit skill in capturing actual power production. Projections to the middle of this century indicate that over most of North America east of the Rocky Mountains annual expected PP is largely unchanged, as is the probability of extended periods of anomalously high or low production. Any small declines in annual PP are of much smaller magnitude than changes due to technological innovation over the last two decades.

Meteorology & Atmospheric Sciences↗
NASA NTRS2024

Extending the Utility of Space-Borne Snow Water Equivalent Observations Over Vegetated Areas With Data Assimilation

Snow is a vital component of the earth system, yet no snow-focused satellite remote sensing platform currently exists. In this study, we investigate how synthetic observations of snow water equivalent (SWE) representative of a synthetic aperture radar remote sensing platform could improve spatiotemporal estimates of snowpack. We use a fraternal twin observing system simulation experiment, specifically investigating how much snow simulated using widely used models and forcing data could be improved by assimilating synthetic observations of SWE. We focus this study across a 24° x 37° domain in the western USA and Canada, simulating snow at 250 m resolution and hourly time steps in water year 2019. We perform two data assimilation experiments, including (1) a simulation excluding synthetic observations in forests where canopies obstruct remote sensing retrievals and (2) a simulation inferring snow distribution in forested grid cells using synthetic observations from nearby canopy-free grid cells. Results found that, relative to a nature run, or assumed true simulation of snow evolution, assimilating synthetic SWE observations improved average SWE biases at maximum snowpack timing in shrub, grass, crop, bare-ground, and wetland land cover types from 14 %, to within 1 %. However, forested grid cells contained a disproportionate amount of SWE volume. In forests, SWE mean absolute errors at the time of maximum snow volume were 111 mm and average SWE biases were on the order of 150 %. Here the data assimilation approach that estimated forest SWE using observations from the nearest canopy-free grid cells substantially improved these SWE biases (18 %) and the SWE mean absolute error (27 mm). Simulations employing data assimilation also improved estimates of the temporal evolution of both SWE and runoff, even in spring snowmelt periods when melting snow and high snow liquid water content prevented synthetic SWE retrievals. In fact, in the Upper Colorado River region, melt-season SWE biases were improved from 63 % to within 1 %, and the Nash–Sutcliffe efficiency of runoff improved from −2.59 to 0.22. These results demonstrate the value of data assimilation and a snow-focused globally relevant remote sensing platform for improving the characterization of SWE and associated water availability.

Justin Pflug↗
DOE OSTI2018

A conservative scheme for electromagnetic simulation of magnetized plasmas with kinetic electrons

A conservative scheme has been formulated and confirmed for gyrokinetic particle simulations of electromagnetic waves and instabilities in magnetized plasmas. An electron continuity equation derived from the drift kinetic equation is used to time advance the electron density perturbation by using the perturbed mechanical flow calculated from the parallel vector potential, and the parallel vector potential is solved by using the perturbed canonical flow from the perturbed distribution function. In gyrokinetic particle simulations using this new scheme, the shear Alfvén wave dispersion relation in the shearless slab and continuum damping in the sheared cylinder have been recovered. The new scheme overcomes the stringent requirement in the conventional perturbative simulation method that perpendicular grid size needs to be as small as electron collisionless skin depth even for the long wavelength Alfvén waves. The new scheme also circumvents the problem in the typical approach that an unphysically large parallel electric field arises due to the inconsistency between electrostatic potential calculated from the perturbed density and vector potential calculated from the perturbed canonical flow. Lastly, the gyrokinetic particle simulations of the Alfvén waves in sheared cylinder have superior numerical properties compared with the fluid simulations, which suffer from numerical difficulties associated with singular mode structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Quantum and Classical Dynamics Simulations of ATP Hydrolysis in Solution

ATP hydrolysis is a key reaction in living cells that drives many cellular processes. The reaction, which involves gamma phosphate cleavage from ATP, converting it to ADP, has been suggested to occur via an associative or dissociative mechanism dependent upon the surrounding environment. Prior quantum chemical studies suffered from short simulation time scales failing to capture free energy contributions due to relaxation of the surrounding aqueous environment. We have developed a highly parallelized QM/MM implementation in the NAMD and OpenAtom simulation packages, using the dual grid, dual length scale method for combined plane-wave and Euler exponential spline-based QM/MM simulations. This approach, using message-driven parallel quantum and classical dynamics, permits sufficient time scale simulations for quantum chemical events such as ATP hydrolysis and is found to accurately and reliably include the free energy contributions of solvent relaxation to hydrolysis. In this paper, we describe the application of the dual grid, dual length plane-wave-based QM/MM method to study both the associative and dissociative mechanisms of ATP hydrolysis, accounting for the free energy contribution from solvent relaxation, as well as for the key role of Mg2+ in the reaction.

Chemistry↗
DOE OSTI2024

Updraft Width Modulates Ambient Atmospheric Controls on Convective Cloud Depth

Abstract The depth of convective clouds affects vertical transport of atmospheric constituents, influencing downstream weather and climate. Atmospheric controls on the maximum depth reached by moist convection are investigated with radar‐tracked convective cells tagged with sounding‐derived atmospheric parameters from a field campaign in central Argentina. Regression analyses show that narrow (<12‐km diameter) and wide (>16‐km diameter) cell depths respond to disparate factors, where cell areas are defined using composite reflectivity signatures. Undiluted lifted parcel indices including convective available potential energy (CAPE) and level of neutral buoyancy (LNB) are top predictors of wide cell maximum depth while mid‐tropospheric relative humidity is the top predictor of narrow cell maximum depth. Because narrow cells are more numerous than wide cells, the overall outcome of the full cell population does not strongly correlate with CAPE and LNB conditions. Tracked cells and atmospheric conditions in a simulation with 3‐km grid spacing covering the field campaign produce similar results to those observed. Narrow cells that are relatively deep have a cooler and moister mid‐troposphere with weaker free tropospheric subsidence, while relatively deep wide cells have much warmer and moister lower tropospheric conditions. These atmospheric differences are present 1 hr before cell initiation at both a fixed observing site and variable cell initiation locations. Simulated narrow cell maximum equivalent potential temperature decreases with height at a rate similar to the ambient vertical gradient, causing these cells to fall short of their LNB and supporting the view that entrainment‐driven dilution is a dominant control on their depth.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

Regional Climate Simulation and Data Assimilation with Variable-Resolution GCMs

Variable resolution GCMs using a global stretched grid (SG) with enhanced regional resolution over one or multiple areas of interest represents a viable new approach to regional climateklimate change and data assimilation studies and applications. The multiple areas of interest, at least one within each global quadrant, include the major global mountains and major global monsoonal circulations over North America, South America, India-China, and Australia. They also can include the polar domains, and the European and African regions. The SG-approach provides an efficient regional downscaling to mesoscales, and it is an ideal tool for representing consistent interactions of globaYlarge- and regionallmeso- scales while preserving the high quality of global circulation. Basically, the SG-GCM simulations are no different from those of the traditional uniform-grid GCM simulations besides using a variable-resolution grid. Several existing SG-GCMs developed by major centers and groups are briefly described. The major discussion is based on the GEOS (Goddard Earth Observing System) SG-GCM regional climate simulations.

Fox-Rabinovitz, Michael S.↗
NASA NTRS1987

Flow solution on a dual-block grid around an airplane

The compressible flow around a complex fighter-aircraft configuration (fuselage, cranked delta wing, canard, and inlet) is simulated numerically using a novel grid scheme and a finite-volume Euler solver. The patched dual-block grid is generated by an algebraic procedure based on transfinite interpolation, and the explicit Runge-Kutta time-stepping Euler solver is implemented with a high degree of vectorization on a Cyber 205 processor. Results are presented in extensive graphs and diagrams and characterized in detail. The concentration of grid points near the wing apex in the present scheme is shown to facilitate capture of the vortex generated by the leading edge at high angles of attack and modeling of its interaction with the canard wake.

Eriksson, Lars-Erik↗
NASA NTRS2018

Generalised Primal-Dual Grids for Unstructured Co-Volume Schemes

The generation of high-quality staggered unstructured grids is considered, leading to the development of a new optimisation-based strategy designed to construct weighted 'Regular-Power' tessellations appropriate for co-volume type numerical discretisation techniques. This new framework aims to extend the conventional Delaunay-Voronoi primal-dual structure; seeking to assemble generalised orthogonal tessellations with enhanced geometric quality. The construction of these grids is motivated by the desire to improve the performance and accuracy of numerical methods based on unstructured co-volume type schemes, including various staggered grid techniques for the simulation of fluid dynamics and hyperbolic transport. In this study, a new hybrid optimisation strategy is proposed; seeking to optimise the geometry, topology and weights associated with general, two-dimensional Regular-Power tessellations using a combination of gradient-ascent and energy-based techniques. The performance of this new method is tested experimentally, with a range of complex, multi-resolution primal-dual grids generated for various coastal and regional ocean modelling applications.

power diagrams↗
NASA NTRS2006

Aerothermodynamic Analyses of Towed Ballutes

A ballute (balloon-parachute) is an inflatable, aerodynamic drag device for application to planetary entry vehicles. Two challenging aspects of aerothermal simulation of towed ballutes are considered. The first challenge, simulation of a complete system including inflatable tethers and a trailing toroidal ballute, is addressed using the unstructured-grid, Navier-Stokes solver FUN3D. Auxiliary simulations of a semi-infinite cylinder using the rarefied flow, Direct Simulation Monte Carlo solver, DSV2, provide additional insight into limiting behavior of the aerothermal environment around tethers directly exposed to the free stream. Simulations reveal pressures higher than stagnation and corresponding large heating rates on the tether as it emerges from the spacecraft base flow and passes through the spacecraft bow shock. The footprint of the tether shock on the toroidal ballute is also subject to heating amplification. Design options to accommodate or reduce these environments are discussed. The second challenge addresses time-accurate simulation to detect the onset of unsteady flow interactions as a function of geometry and Reynolds number. Video of unsteady interactions measured in the Langley Aerothermodynamic Laboratory 20-Inch Mach 6 Air Tunnel and CFD simulations using the structured grid, Navier-Stokes solver LAURA are compared for flow over a rigid spacecraft-sting-toroid system. The experimental data provides qualitative information on the amplitude and onset of unsteady motion which is captured in the numerical simulations. The presence of severe unsteady fluid - structure interactions is undesirable and numerical simulation must be able to predict the onset of such motion.

Gnoffo, Peter A.↗
NASA NTRS2000

NPSS on NASA's Information Power Grid: Using CORBA and Globus to Coordinate Multidisciplinary Aeroscience Applications

This paper describes a project to evaluate the feasibility of combining Grid and Numerical Propulsion System Simulation (NPSS) technologies, with a view to leveraging the numerous advantages of commodity technologies in a high-performance Grid environment. A team from the NASA Glenn Research Center and Argonne National Laboratory has been studying three problems: a desktop-controlled parameter study using Excel (Microsoft Corporation); a multicomponent application using ADPAC, NPSS, and a controller program-, and an aviation safety application running about 100 jobs in near real time. The team has successfully demonstrated (1) a Common-Object- Request-Broker-Architecture- (CORBA-) to-Globus resource manager gateway that allows CORBA remote procedure calls to be used to control the submission and execution of programs on workstations and massively parallel computers, (2) a gateway from the CORBA Trader service to the Grid information service, and (3) a preliminary integration of CORBA and Grid security mechanisms. We have applied these technologies to two applications related to NPSS, namely a parameter study and a multicomponent simulation.

Lopez, Isaac↗