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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem

The equilibrium configuration of a plasma in an axially symmetric reactor is described mathematically by a free boundary problem associated with the celebrated Grad-Shafranov equation. The presence of uncertainty in the model parameters introduces the need to quantify the variability in the predictions. This is often done by computing a large number of model solutions on a computational grid for an ensemble of parameter values and then obtaining estimates for the statistical properties of solutions. In this study, we explore the savings that can be obtained using multilevel Monte Carlo methods, which reduce costs by performing the bulk of the computations on a sequence of spatial grids that are coarser than the one that would typically be used for a simple Monte Carlo simulation. We examine this approach using both a set of uniformly refined grids and a set of adaptively refined grids guided by a discrete error estimator. Numerical experiments show that multilevel methods dramatically reduce the cost of simulation, with cost reductions typically on the order of 60 or more and possibly as large as 200. Furthermore, adaptive griding results in more accurate computation of geometric quantities such as x-points associated with the model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2003

Effects of Wing Leading Edge Penetration with Venting and Exhaust Flow from Wheel Well at Mach 24 in Flight

A baseline solution for CFD Point 1 (Mach 24) in the STS-107 accident investigation was modified to include effects of: (1) holes through the leading edge into a vented cavity; and (2) a scarfed, conical nozzle directed toward the centerline of the vehicle from the forward, inboard corner of the landing gear door. The simulations were generated relatively quickly and early in the investigation because simplifications were made to the leading edge cavity geometry and an existing utility to merge scarfed nozzle grid domains with structured baseline external domains was implemented. These simplifications in the breach simulations enabled: (1) a very quick grid generation procedure; and (2) high fidelity corroboration of jet physics with internal surface impingements ensuing from a breach through the leading edge, fully coupled to the external shock layer flow at flight conditions. These simulations provided early evidence that the flow through a two-inch diameter (or larger) breach enters the cavity with significant retention of external flow directionality. A normal jet directed into the cavity was not an appropriate model for these conditions at CFD Point 1 (Mach 24). The breach diameters were of the same order or larger than the local, external boundary-layer thickness. High impingement heating and pressures on the downstream lip of the breach were computed. It is likely that hole shape would evolve as a slot cut in the direction of the external streamlines. In the case of the six-inch diameter breach the boundary layer is fully ingested. The intent of externally directed jet simulations in the second scenario was to approximately model aerodynamic effects of a relatively large internal wing pressure, fueled by combusting aluminum, which deforms the corner of the landing gear door and directs a jet across the windside surface. These jet interactions, in and of themselves, were not sufficiently large to explain observed aerodynamic behavior.

Gnoffo, Peter A.↗
DOE OSTI2014

Greenland Surface Mass Balance as Simulated by the Community Earth System Model. Part II: Twenty-First-Century Changes

This study presents the first twenty-first-century projections of surface mass balance (SMB) changes for the Greenland Ice Sheet (GIS) with the Community Earth System Model (CESM), which includes a new ice sheet component. For glaciated surfaces, CESM includes a sophisticated calculation of energy fluxes, surface albedo, and snowpack hydrology (melt, percolation, refreezing, etc.). To efficiently resolve the high SMB gradients at the ice sheet margins and provide surface forcing at the scale needed by ice sheet models, the SMB is calculated at multiple elevations and interpolated to a finer 5-km ice sheet grid. During a twenty-first-century simulation driven by representative concentration pathway 8.5 (RCP8.5) forcing, the SMB decreases from 372 ± 100 Gt yr -1 in 1980–99 to -78 ± 143 Gt yr -1 in 2080–99. The 2080–99 near-surface temperatures over the GIS increase by 4.7 K (annual mean) with respect to 1980–99, only 1.3 times the global increase (+3.7 K). Snowfall increases by 18%, while surface melt doubles. The ablation area increases from 9% of the GIS in 1980–99 to 28% in 2080–99. Over the ablation areas, summer downward longwave radiation and turbulent fluxes increase, while incoming shortwave radiation decreases owing to increased cloud cover. The reduction in GIS-averaged July albedo from 0.78 in 1980–99 to 0.75 in 2080–99 increases the absorbed solar radiation in this month by 12%. Summer warming is strongest in the north and east of Greenland owing to reduced sea ice cover. In the ablation area, summer temperature increases are smaller due to frequent periods of surface melt.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Real-time Simulation Model of a Vanadium Redox Flow Battery Energy Storage System with Power Electronics Integration

This paper presents a real-time (RT) model of a 20 kW multi-stack vanadium redox flow battery (VRFB) system including power electronics converters in a single real-time framework. Implemented on a Typhoon HIL 604 device, this model is designed to simulate the operation of the flow battery with system-level controls. To provide a realistic overview of a VRFB system, the model incorporates a closed-loop control system for managing the speeds of two centrifugal pumps and an electro-thermal model of the battery stacks, enabling realtime temperature estimation. Furthermore, a two-stage power electronics topology and associated control solution are presented, operating with distinct strategies for interfacing the battery model with the main grid and a local load. Results from simulations demonstrate that the VRFB model is well suited for RT environment, and system can manage power during both charging and discharging cycles, while a SCADA monitor is used to display critical variables for safe, efficient, and reliable operation of the VRFB systems.

Rezende Da Costa Reis Kimpara, Renata [ORNL] (ORCI↗
NASA NTRS2011

Unstructured CFD Aerodynamic Analysis of a Generic UCAV Configuration

Three independent studies from the United States (NASA), Sweden (FOI), and Australia (DSTO) are analyzed to assess the state of current unstructured-grid computational fluid dynamic tools and practices for predicting the complex static and dynamic aerodynamic and stability characteristics of a generic 53-degree swept, round-leading-edge uninhabited combat air vehicle configuration, called SACCON. NASA exercised the USM3D tetrahedral cell-centered flow solver, while FOI and DSTO applied the FOI/EDGE general-cell vertex-based solver. The authors primarily employ the Reynolds Averaged Navier-Stokes (RANS) assumption, with a limited assessment of the EDGE Detached Eddy Simulation (DES) extension, to explore sensitivities to grids and turbulence models. Correlations with experimental data are provided for force and moments, surface pressure, and off-body flow measurements. The vortical flow field over SACCON proved extremely difficult to model adequately. As a general rule, the prospect of obtaining reasonable correlations of SACCON pitching moment characteristics with the RANS formulation is not promising, even for static cases. Yet, dynamic pitch oscillation results seem to produce a promising characterization of shapes for the lift and pitching moment hysteresis curves. Future studies of this configuration should include more investigation with higher-fidelity turbulence models, such as DES.

Frink, Neal T.↗
DOE OSTI2015

Multi-Lab EV Smart Grid Integration Requirements Study. Providing Guidance on Technology Development and Demonstration

The report begins with a discussion of the current state of the energy and transportation systems, followed by a summary of some VGI scenarios and opportunities. The current efforts to create foundational interface standards are detailed, and the requirements for enabling PEVs as a grid resource are presented. Existing technology demonstrations that include vehicle to grid functions are summarized. The report also includes a data-based discussion on the magnitude and variability of PEVs as a grid resource, followed by an overview of existing simulation tools that vi This report is available at no cost from the National Renewable Energy Laboratory (NREL) at www.nrel.gov/publications. can be used to explore the expansion of VGI to larger grid functions that might offer system and customer value. The document concludes with a summary of the requirements and potential action items that would support greater adoption of VGI.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2018

Demeter – A Land Use and Land Cover Change Disaggregation Model

Demeter is an open source Python package that was built to disaggregate projections of future land allocations generated by an integrated assessment model (IAM). Projected land allocation from IAMs is traditionally transferred to Earth System Models (ESMs) in a variety of gridded formats and spatial resolutions as inputs for simulating biophysical and biogeochemical fluxes. Existing tools for performing this translation generally require a number of manual steps which introduces error and is inefficient. Demeter makes this process seamless and repeatable by providing gridded and land cover change (LULCC) products derived directly from an IAM—in this case, the Global Change Assessment Model (GCAM)—in a variety of formats and resolutions commonly used by ESMs. Demeter is publicly available via GitHub and has an extensible output module allowing for future ESM needs to be easily accommodated.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A correction procedure for self-induced velocity of a finite-sized particle in two-way coupled Euler–Lagrange simulations

The importance of incorporating a correction to undo the self-induced perturbation velocity of a particle, when its size becomes comparable to the Eulerian grid, in a two-way coupled Euler–Lagrange (EL) simulation is now well appreciated. The present work improves upon the prior correction procedures in a few important ways. First, the past correction procedures have been scalar-based with the assumption that the quasi-steady force is the source of self-induced velocity perturbation. Here we generalize to a vector correction procedure and thereby the directions of feedback force and relative velocity can be different. This allows the correction procedure to be used even in the presence of added-mass, history, and lift forces. Second, the effect of a nearby wall has been systematically included in the correction procedure. The correction procedure depends on fundamental Oseen solutions of streamwise and transverse regularized feedback forces. We present a Fourier transform-based analytical approach to obtaining these regularized Oseen solutions. We also present a step-by-step numerical procedure for obtaining the Oseen solutions in any EL code. With the analytical or numerical Oseen functions, the correction procedure can be easily implemented in any EL code. Iterations are required in solving the implicit correction equations and it is demonstrated that the correction procedure converges rapidly within three or four iterations. In conclusion, a simple empirical approach is also presented to account for unsteady effects in the correction procedure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

An unstructured C-grid type variational formulation for the sea ice dynamics

Historically, B-grid formulations of sea ice dynamics have been dominant because they have matched the grid type used by ocean models. The reason for the grid match is simple – it facilitates penetration of the curl of ice-ocean stress into the deep ocean with minimal numerical diffusivity because sea ice and ocean velocity are co-located. In recent years, as ocean models have increasingly progressed to C-grids, sea ice models have followed suit on quadrilateral meshes, but the implementation of an unstructured C-grid sea ice models is new. We present an unstructured C-grid discretization of the Elastic Viscous Plastic (EVP) rheology, where the velocity unknowns are discretized at the edges of the mesh cells with n-sides, where typically n is greater than or equal to four, rather than at the vertices, as in the B-grid. Our framework of choice is the Model for Prediction Across Scales (MPAS) within E3SM, the climate model of the U.S. Department of Energy, although our approach is general and could be applied to other models as well. While MPAS-Seaice is currently defined on a B-grid, MPAS-Ocean runs on a C-grid, hence interpolation operators are heavily used when coupled simulations are performed. In this work, we describe a mathematical formulation to transition the dynamics of MPAS-Seaice to a C-grid, in order to ultimately facilitate the coupling with MPAS-Ocean and reduce numerical errors associated with this communication. Numerical results are reported to highlight the features of the method.

58 GEOSCIENCES↗
NASA NTRS2004

Development of a Flexible Framework for Hypersonic Navier-Stoke Space Shuttle Orbiter Meshes

A flexible framework constructing block structured volume grids for hypersonic Navier-Strokes flow simulations was developed for the analysis of the Shuttle Orbiter Columbia. The development of the framework, which was partially basedon the requirements of the primary flow solvers used resulted in an ability to directly correlate solutions contributed by participating groups on a common surface mesh. A foundation was built through the assessment of differences between differnt solvers, which provided confidence for independent assessment of other damage scenarios by team members. The framework draws on the experience of NASA Langley and NASA Ames Research Centers in structured grid generation, and consists of a grid generation, and consist of a grid generation process implemented through a division of responsibilities. The nominal division of labor consisted of NASA Johnson Space Center coordinating the damage scenarios to be analyzed by the Aerothermodynamics Columbia Accident Investigation (ACAI) team, Ames developing the surface grids that described the computational volume about the Orbiter, and Langley improving grid quality of Ames generated data and constructing the final computational volume grids. Distributing the work among the participant in th ACAI team resulted in significantl less time required to construct complete meshes than possible by any individual participant. The approach demonstrated that the One-NASA grid generation team could sustain the demand of for five new meshes to explore new damage scenarios within an aggressive time-line.

Alter, Stephen J.↗
DOE OSTI2021

MIRACL Co‐Simulation platform for control and operation of distributed wind in microgrid

Abstract This paper presents a co‐simulation platform (CSP) developed as a part of the “Microgrids, Infrastructure Resilience and Advanced Controls Launchpad (MIRACL)” project, hereafter called MIRACL‐CSP, to allow simulation‐based testing, demonstration, and evaluations of distributed wind under various grid operating conditions. MIRACL‐CSP provides modular interactions among the power distribution system, distributed wind, and utility decision‐making framework. A hierarchical engine for large‐scale infrastructure co‐simulations (HELICS) is used as the core engine of MIRACL‐CSP to establish time and information coordination among the MIRACL‐CSP modules. In this study, MIRACL‐CSP performance is demonstrated using the IEEE 123‐node test distribution grid modelled in GridLAB‐D and the utility decision support application modelled in Python. The functionality of MIRACL‐CSP is demonstrated through various grid operational scenarios in a microgrid and networked microgrid environment. A co‐simulation study is described that demonstrated MIRACL‐CSP capabilities for the microgrid operational scenario in the IEEE 123‐node test system. MIRACL‐CSP is a generic platform that facilitates distributed wind research for different test systems, applications, and valuation of distributed wind.

17 WIND ENERGY↗
DOE OSTI2016

Effects on image quality of a 2D antiscatter grid in x-ray digital breast tomosynthesis: Initial experience using the dual modality (x-ray and molecular) breast tomosynthesis scanner

Purpose: Radiation scattered from the breast in digital breast tomosynthesis (DBT) causes image degradation, including loss of contrast between cancerous and background tissue. Unlike in 2-dimensional (2D) mammography, an antiscatter grid cannot readily be used in DBT because changing alignment between the tube and detector during the scan would result in unacceptable loss of primary radiation. However, in the dual modality breast tomosynthesis (DMT) scanner, which combines DBT and molecular breast tomosynthesis, the tube and detector rotate around a common axis, thereby maintaining a fixed tube-detector alignment. This C-arm geometry raises the possibility of using a 2D (cellular) focused antiscatter grid. The purpose of this study is to assess change in image quality when using an antiscatter grid in the DBT portion of a DMT scan under conditions of fixed radiation dose. Methods: Two 2D focused prototype grids with 80 cm focal length were tested, one stack-laminated from copper (Cu) and one cast from a tungsten-polymer (W-poly). They were reciprocated using a motion scheme designed to maximize transmission of primary x-ray photons. Grid-in and grid-out scatter-to-primary ratios (SPRs) were measured for rectangular blocks of material simulating 30%, 50%, and 70% glandular tissue compositions. For assessment of changes in image quality through the addition of a grid, the Computerized Imaging Reference Systems, Inc., phantom Model 011A containing a set of 1 cm thick blocks simulating a range of glandular/adipose ratios from 0/100 to 100/0 was used. To simulate 6.5 and 8.5 cm thick compressed breasts, 1 cm thick slices of PMMA were added to the Model 011A phantom. DBT images were obtained with and without the grid, with exposure parameters fixed for a given compressed thickness. Signal-difference-to-noise ratios (SDNRs), contrast, and voxel value-based attenuation coefficients (μ) were measured for all blocks from reconstructed phantom images. Results: For 4, 6, and 8 cm tissue-equivalent block phantom thicknesses, the inclusion of the W-poly grid reduced the SPR by factors of 5, 6, and 5.8, respectively. For the same thicknesses, the copper grid reduced the SPR by factors of 3.9, 4.5, and 4.9. For the 011A phantom, the W-poly grid raised the SDNR of the 70/30 block from 0.8, −0.32, and −0.72 to 0.9, 0.76, and 0.062 for the 4.5, 6.5, and 8.5 cm phantoms, respectively. It raised the SDNR of the 100/0 block from 3.78, 1.95, and 1.0 to 3.79, 3.67, and 3.25 for the 4.5, 6.5, and 8.5 cm phantoms, respectively. Inclusion of the W-poly grid improved the accuracy of image-based μ values for all block compositions. However, smearing of attenuation across slices due to limited angular sampling decreases the sensitivity of voxel values to changing composition compared to theoretical μ values. Conclusions: Under conditions of fixed radiation dose to the breast, use of a 2D focused grid increased contrast, SDNR, and accuracy of estimated attenuation for mass-simulating block compositions in all phantom thicknesses tested, with the degree of improvement depending upon material composition. A 2D antiscatter grid can be usefully incorporated in DBT systems that employ fully isocentric tube-detector rotation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Application of High Performance Computing for Studying Cyclic Variability in Dilute Internal Combustion Engines

Combustion instabilities in dilute internal combustion engines are manifest in cyclic variability (CV) in engine performance measures such as integrated heat release or shaft work. Understanding the factors leading to CV is important in model-based control, especially with high dilution where experimental studies have demonstrated that deterministic effects can become more prominent. Observation of enough consecutive engine cycles for significant statistical analysis is standard in experimental studies but is largely wanting in numerical simulations because of the computational time required to compute hundreds or thousands of consecutive cycles. We have proposed and begun implementation of an alternative approach to allow rapid simulation of long series of engine dynamics based on a low-dimensional mapping of ensembles of single-cycle simulations which map input parameters to output engine performance. This paper details the use Titan at the Oak Ridge Leadership Computing Facility to investigate CV in a gasoline direct-injected spark-ignited engine with a moderately high rate of dilution achieved through external exhaust gas recirculation. The CONVERGE™ CFD software was used to perform single-cycle simulations with imposed variations of operating parameters and boundary conditions selected according to a sparse grid sampling of the parameter space. Using an uncertainty quantification technique, the sampling scheme is chosen similar to a design of experiments grid but uses algorithms designed to minimize the number of samples required to achieve a desired degree of accuracy. The simulations map input parameters to output metrics of engine performance for a single cycle, and by mapping over a large parameter space, results can be interpolated from within that space. This interpolation scheme forms the basis for a low-dimensional ‘metamodel’ (or model of a model) which can be used to mimic the dynamical behavior of corresponding high-dimensional simulations. Simulations of high-EGR spark-ignition combustion cycles within a parametric sampling grid were performed and analyzed statistically, and sensitivities of the physical factors leading to high CV are presented. With these results, the prospect of producing low-dimensional metamodels to describe engine dynamics at any point in the parameter space will be discussed. Additionally, modifications to the methodology to account for nondeterministic effects in the numerical solution environment are proposed.

Finney, Charles E. A.↗
NASA NTRS1991

AFE base flow computations

Hypersonic wake flows behind the Aeroassist Flight Experiment (AFE) geometry are analyzed using two Navier-Stokes flow solvers. Many of the AFE wake features observed in ballistic-range shadowgraphs are simulated using a simple, two-dimensional semicylinder geometry at moderate angles of attack. At free-stream conditions corresponding to a Hypersonic Free Flight Facility (HFFF) AFE experiment, the three-dimensional base flow for the AFE geometry is computed using an ideal-gas, Navier-Stokes solver. The computed results agree reasonably well with the shadowgraphs taken at the HFFF. An ideal-gas and a nonequilibrium Navier-Stokes solver have been coupled and applied to the complete flow around the AFE vehicle at the free-stream conditions corresponding to a nomial trajectory point. Limitations of the coupled ideal-gas and nonequilibrium solution are discussed. The nonequilibrium base flow solution is analyzed for the wake radiation and the radiation profiles along various lines of sight are compared. Finally, the wake unsteadiness is predicted using experimental correlations and the numerical solutions. An adaptive grid code, SAGE, has been used in all the simulations to enhance the solution accuracy. The grid adaptation is found to be necessary in obtaining base flow solutions with accurate flow features.

Venkatapathy, Ethiraj↗
DOE OSTI2024

Improving the estimate of higher-order moments from lidar observations near the top of the convective boundary layer

Abstract. Ground-based lidar data have proven extremely useful for profiling the convective boundary layer (CBL). Many groups have derived higher-order moments (e.g., variance, skewness, fluxes) from high-temporal-resolution lidar data using an autocovariance approach. However, these analyses are highly uncertain near the CBL top when the depth of the CBL (zi) is changing during the analysis period. This is because the autocovariance approach is usually applied to constant height levels and the character of the eddies is changing on either side of the changing CBL top. Here, a new approach is presented wherein the autocovariance analysis is performed on a normalized height grid, with a temporally smoothed zi. Output from a large eddy simulation model demonstrates that deriving higher-order moments from time series on a normalized height grid has better agreement with the slab-averaged quantities than the moments derived from the original height grid.

Rosenberger, Tessa E. (ORCID:0000000333205873)↗
NASA NTRS1992

Experimental investigation of a simulated LOX injector flow field and other nonintrusive measurement efforts

Efforts to improve the characteristics of fuel-oxidizer mixing in liquid rocket combustors have lead to a swirl element design for a liquid oxygen injector which is being considered for use on the STME. For the design which is the subject of this investigation, the oxygen enters the injector element perpendicular to the injector axis and nearly tangent to the circular injector wall. This swirl element is at one end of a tube and the injector exit is at the other. This geometric configuration creates a plume in the shape of a conical sheet. This sheet is either primarily contiguous liquid or droplets depending on the pressure drop in the injector and the distance from the injector exit. Probe-based devices such as two-dimensional grid patternators have been used to investigate simulated LOX injector flow fields (Hulka). The primary work described herein is an effort to use optical techniques to investigate the plume of a swirl injector element. For this investigation, a high pressure (500 psig) cold flow test facility was constructed. Water was used as the LOX simulate and air pressure was used to drive the injector flow field. Laser-induced fluorescence (LIF) from dye seeded into the water was used to obtain quantitative measurements of the time-averaged water concentration distribution in the plume. Scattered laser light and LIF were used for time averaged plume visualization and scattered light from a strobe with a 1 microsecond pulse was used for time-resolved plume visualization. During the Summer Faculty Fellowship for which this report was developed, an additional effort, unrelated to the swirl injector investigation, was made to resolve fluctuations in the combustion product composition in the exhaust of a hybrid rocket motor. A brief description of this effort is included herein.

Hartfield, Roy J., Jr.↗
DOE OSTI2015

Multiscale Modeling of Porous Medium Systems. In: Handbook of Porous Media, Third edition

The use of computers over the last 60 years for simulating subsurface “ow problems has provided an impressive, ever-expanding ability to model processes at a high resolution. Whereas in the 1960s a computational grid with 500 œxed nodes would have pushed the boundaries of computer power [59], simulations today involve millions of spatial grid points and adaptive meshes [40,57].

Dye, Amanda↗
DOE OSTI2023

Mass‐Conserving Downscaling of Climate Model Precipitation Over Mountainous Terrain for Water Resource Applications

Abstract A mass‐conserving method to downscale precipitation from global climate models (GCMs) using sub‐grid‐scale topography and modeled 700‐mb wind direction is presented. Runoff is simulated using a stand‐alone hydrological model, with this and several other methods as inputs, and compared to runoff simulated using historical observations over the western contiguous United States. Results suggest the mitigation of grid‐scale biases is more critical than downscaling for some regions with large wet biases (e.g., the Great Basin and Upper Colorado). In other regions (e.g., the Pacific Northwest) the new method produces more realistic sub‐grid‐scale variability in runoff compared to unadjusted GCM output and a simpler downscaling method. The presented method also brings the runoff centroid timing closer to that simulated with observations for all subregions examined.

Geology↗