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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Numerical simulations of unsteady, viscous, transonic flow over isolated and cascaded airfoils using a deforming grid

A compressible, unsteady, full Navier-Stokes, finite difference code was developed for modeling transonic flow through two-dimensional, oscillating cascades. The procedure introduces a deforming grid technique to capture the motion of the airfoils. Results using a deforming grid are presented for both isolated and cascaded airfoils. The load histories and unsteady pressure distributions are predicted for the NASA 64A010 isolated airfoil and compared with existing experimental data. Results show that the deforming grid technique can be used to successfully predict the unsteady flow properties around an oscillating airfoil. The deforming grid technique was extended for modeling unsteady flow in a cascade. The use of a deforming grid simplifies the specification of boundary conditions. Unsteady flow solutions similar to the isolated airfoil predictions are found for a NACA 0012 cascade with zero interblade phase angle and zero stagger. Experimental data for these cases are not available for code validation, but computational results are presented to show sample predictions from the code. Applications of the code to typical turbomachinery flow conditions will be presented in future work.

Huff, Dennis L.↗
NASA NTRS1987

Numerical simulations of unsteady, viscous, transonic flow over isolated and cascaded airfoils using a deforming grid

A compressible, unsteady, full Navier-Stokes, finite difference code was developed for modeling transonic flow through two-dimensional, oscillating cascades. The procedure introduces a deforming grid technique to capture the motion of the airfoils. Results using a deforming grid are presented for both isolated and cascaded airfoils. The load histories and unsteady pressure distributions are predicted for the NASA 64A010 isolated airfoil and compared with existing experimental data. Results show that the deforming grid technique can be used to successfully predict the unsteady flow properties around an oscillating airfoil. The deforming grid technique was extended for modeling unsteady flow in a cascade. The use of a deforming grid simplifies the specification of boundary conditions. Unsteady flow solutions similar to the isolated airfoil predictions are found for a NACA 0012 cascade with zero interblade phase angle and zero stagger. Experimental data for these cases are not available for code validation, but computational results are presented to show sample predictions from the code. Applications of the code to typical turbomachinery flow conditions will be presented in future work.

Huff, Dennis L.↗
NASA NTRS1994

Numerical Simulation of Inlet Bleed with Circular Holes on Plate Under Shock-Wave/Boundary-Layer Interactions

A numerical study was performed to investigate the shock-wave/boundary-layer interactions on a flat plate with bleed through one or more circular holes that vent into a plenum. The bleed-hole patterns considered for the study include in-line multiple holes and staggered multiple-row holes that are configured to simulate the patterns used in inlet bleed systems of high performance aircraft. The focus of the study was to examine how the bleed through multiple holes affect bleed rate and the pressure and Mach number distributions. Since the bleed performance was found sensitive to the change in bleed conditions, a computational procedure was developed to give a good turnaround computational time for parametric studies involving changes in bleed hole geometry and the structure of shock-wave/boundary-layer flowfield. The procedure includes the grid-generation methodology and the flow simulation with solutions from the Navier-Stokes equations. The computational techniques permit analysis of complex bleed systems and make possible the investigation of a broader range of design variables associated with inlet bleed operation.

Chyu, Wei J.↗
NASA NTRS2003

Source Term Model for an Array of Vortex Generator Vanes

A source term model was developed for numerical simulations of an array of vortex generators. The source term models the side force created by a vortex generator being modeled. The model is obtained by introducing a side force to the momentum and energy equations that can adjust its strength automatically based on a local flow. The model was tested and calibrated by comparing data from numerical simulations and experiments of a single low-profile vortex generator vane, which is only a fraction of the boundary layer thickness, over a flat plate. The source term model allowed a grid reduction of about seventy percent when compared with the numerical simulations performed on a fully gridded vortex generator without adversely affecting the development and capture of the vortex created. The source term model was able to predict the shape and size of the stream wise vorticity and velocity contours very well when compared with both numerical simulations and experimental data.

Buning, P. G.↗
NASA NTRS2012

Numerical Dissipation and Wrong Propagation Speed of Discontinuities for Stiff Source Terms

In compressible turbulent combustion/nonequilibrium flows, the constructions of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and (b) obtaining correct propagation speed of discontinuities for stiff reacting terms on coarse grids share one important ingredient - minimization of numerical dissipation while maintaining numerical stability. Here coarse grids means standard mesh density requirement for accurate simulation of typical non-reacting flows. This dual requirement to achieve both numerical stability and accuracy with zero or minimal use of numerical dissipation is most often conflicting for existing schemes that were designed for non-reacting flows. The goal of this paper is to relate numerical dissipations that are inherited in a selected set of high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities as a function of stiffness of the source term and the grid spacing.

Yee, H. C.↗
DOE OSTI2018

Multiphase Large Eddy Simulation of a Pulsed Sieve-Plate Extraction Column - 18125

Presented here is research resulting from a framework of single-phase and multiphase modelling trials. These are intended to accurately simulate the internal flow phenomena arising in pulsed sieve-plate extraction columns representative of those found in industry. Three-dimensional unsteady flow calculations using large eddy simulation coupled with dynamic sub-grid scale modelling have been used to capture the transient effects of sinusoidal flow conditions using the open source CFD code OpenFOAM{sup R}. Multiphase interface interactions are modelled using the volume of fluid method with appropriate heavy and light phase fluidic properties capturing surface tension, density and viscous effects. The system studied consists of a 150 mm diameter column with two perforated plates with a fractional free area of 25 %, operational under 1 Hz pulsing frequency at a 50 mm amplitude and plate spacing of 300 mm. The materials simulated include water (single-phase) and 3 M nitric acid and dodecane / tributyl phosphate mixture at 30 vol% (multiphase). Results are compared with appropriate correlations from literature. The dispersive mixing efficiency of the system was also evaluated via a mixing index value which allows quantification of the ratio of rotational and irrotational flow components. This study has found that the predictive methods employed are effective at capturing the different scales of turbulence present within pulsed sieve-plate extraction columns, and its influence on their operation. As such, the work highlights the ability of large eddy simulation-based modelling to predict the complex flows in such column geometries. Furthermore, the techniques employed provide valuable insight into the hydrodynamic flow characteristics and mixing mechanisms arising within pulsed sieve-plate extraction columns. (authors)

42 ENGINEERING↗
DOE OSTI2024

Impacts of Control, Penetration, and Distribution of Embedded Storage Network in Bulk Power System

The current shift in generation mix from fossil fuel plants towards variable and intermittent renewable energy sources is poised to create a future grid with reduced physical inertia and mismatch between generation and demand. Embedded storage, which is a concept of a coordinated network of storage units sited at the interface between the transmission and distribution system, is proposed as a mechanism to provide a buffer between generation and demand. This paper proposes an automated framework to model and integrate embedded storage in large-scale power systems with industry-grade grid-following (GFL) and grid-forming (GFM) control technologies. More importantly, the developed framework is used to explore the impacts of embedded storage control, penetration, location, and capacity in providing fast frequency response to the grid under contingency events such as generator trips and faults. The framework and study are conducted using the transient-stability simulation tool PSS/E and a realistic model of the Puerto Rico grid as a chosen test system. The simulation results show that GFL and GFM embedded storage, distributed throughout the system, with sufficient penetration and capacity, can effectively improve primary frequency response of the system under the studied contingency events.

Battery Energy Storage, embedded storage, grid-for↗
DOE OSTI2023

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges

Abstract High-fidelity simulation of transitional and turbulent flows over multi-scale surface roughness presents several challenges. For instance, the complex and irregular geometrical nature of surface roughness makes it impractical to employ conforming structured grids, commonly adopted in large-scale numerical simulations due to their high computational efficiency. One possible solution to overcome this problem is offered by immersed boundary methods, which allow wall boundary conditions to be enforced on grids that do not conform to the geometry of the solid boundary. To this end, a three-dimensional, second-order accurate boundary data immersion method (BDIM) is adopted. A novel mapping algorithm that can be applied to general three-dimensional surfaces is presented, together with a newly developed data-capturing methodology to extract and analyze on-surface flow quantities of interest. A rigorous procedure to compute gradient quantities such as the wall shear stress and the heat flux on complex non-conforming geometries is also introduced. The new framework is validated by performing a direct numerical simulation (DNS) of fully developed turbulent channel flow over sinusoidal egg-carton roughness in a minimal-span domain. For this canonical case, the averaged streamwise velocity profiles are compared against results from the literature obtained with a body-fitted grid. General guidelines on the BDIM resolution requirements for multi-scale roughness simulation are given. Momentum and energy balance methods are used to validate the calculation of the overall skin friction and heat transfer at the wall. The BDIM is then employed to investigate the effect of irregular homogeneous surface roughness on the performance of an LS-89 high-pressure turbine blade at engine-relevant conditions using DNS. This is the first application of the BDIM to realize multi-scale roughness for transitional flow in transonic conditions in the context of high-pressure turbines. The methodology adopted to generate the desired roughness distribution and to apply it to the reference blade geometry is introduced. The results are compared to the case of an equivalent smooth blade.

Engineering↗
DOE OSTI2024

Budget-Constrained Sizing of Renewable and Energy Storage Systems for Farm-Scale Ammonia Production Within the Food-Energy-Water Nexus

In the transition toward sustainable agriculture, farms have emerged as eco-friendly pioneers, harnessing clean hybrid wind and solar systems to improve farm performance. A concern in this paradigm is the effective sizing of renewable energy systems to ensure optimal energy use within budget considerations. This research focuses on optimizing renewable energy sizing in small-scale ammonia production to meet specific farm demands and enhance local resilience, emphasizing the interplay between environmental and economic factors. These findings promise increased energy efficiency and sustainability in this innovative agricultural sector. Additionally, our approach considers small-scale ammonia plant needs and the dynamic relationships between ammonia, water, and farm demands. Simulations demonstrate substantial cost savings in farm electricity consumption. Specifically, scenarios with renewable energy integration in the farm can reduce at least 13% electricity cost compared to a grid-dependent system in the 15-year simulation.

ammonia↗
DOE OSTI2014

Electron attenuation in free, neutral ethane clusters

The electron effective attenuation length (EAL) in free, neutral ethane clusters has been determined at 40 eV kinetic energy by combining carbon 1s x-ray photoelectron spectroscopy and theoretical lineshape modeling. More specifically, theory is employed to form model spectra on a grid in cluster size (N) and EAL (λ), allowing N and λ to be determined by optimizing the goodness-of-fit χ{sup 2}(N, λ) between model and observed spectra. Experimentally, the clusters were produced in an adiabatic-expansion setup using helium as the driving gas, spanning a range of 100–600 molecules in mean cluster size. The effective attenuation length was determined to be 8.4 ± 1.9 Å, in good agreement with an independent estimate of 10 Å formed on the basis of molecular electron-scattering data and Monte Carlo simulations. The aggregation state of the clusters as well as the cluster temperature and its importance to the derived EAL value are discussed in some depth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

A priori mesh quality metrics for three-dimensional hybrid grids

Use of general hybrid grids to attain complex-geometry field simulations poses a challenge on estimation of their quality. Apart from the typical problems of non-uniformity and non-orthogonality, the change in element topology is an extra issue to address. The present work derives and evaluates an a priori mesh quality indicator for structured, unstructured, as well as hybrid grids consisting of hexahedra, prisms, tetrahedra, and pyramids. Emphasis is placed on deriving a direct relation between the quality measure and mesh distortion. The work is based on use of the Finite Volume discretization for evaluation of first order spatial derivatives. The analytic form of the truncation error is derived and applied to elementary types of mesh distortion including typical hybrid grid interfaces. The corresponding analytic expressions provide relations between the truncation error and the degree of stretching, skewness, shearing, torsion, expansion, as well as the type of grid interface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Characteristics of Precipitation and Mesoscale Convective Systems Over the Peruvian Central Andes in Multi 5‐Year Convection‐Permitting Simulations

Abstract Using the Weather Research and Forecasting model with two planetary boundary layer schemes, ACM2 and MYNN, convection‐permitting model (CPM) regional climate simulations were conducted for a 6‐year period, including a one‐year spin‐up period, at a 15‐km grid spacing covering entire South America and a nested convection‐permitting 3‐km grid spacing covering the Peruvian central Andes region. These two CPM simulations along with a 4‐km simulation covering South America produced by National Center for Atmospheric Research (NCAR), three gridded precipitation products, and rain gauge data in Peru and Brazil, are used to document the characteristics of precipitation and mesoscale convective systems (MCSs) in the Peruvian central Andes region. Results show that all km‐scale simulations generally capture the spatiotemporal patterns of precipitation and MCSs at both seasonal and diurnal scales, although biases exist in aspects such as precipitation intensity and MCS frequency, size, propagation speed, and associated precipitation intensity. The 3‐km simulation using MYNN scheme generally outperforms the other simulations in capturing seasonal and diurnal precipitation over the mountain, while both it and the 4‐km simulation demonstrate superior performance in the western Amazon Basin, based on the comparison to the gridded precipitation products and gauge data. Dynamic factors, primarily low‐level jet and terrain‐induced uplift, are the key drivers for precipitation and MCS genesis along the east slope of the Andes, while thermodynamic factors control the precipitation and MCS activity in the western Amazon Basin and over elevated mountainous regions. The study suggests model improvements and better model configurations for future regional climate projections.

Meteorology & Atmospheric Sciences↗
DOE OSTI2014

THE SIGNATURE OF INITIAL CONDITIONS ON MAGNETOHYDRODYNAMIC TURBULENCE

We demonstrate that the initial correlation between velocity and current density fluctuations can lead to the formation of enormous current sheets in freely evolving magnetohydrodynamic (MHD) turbulence. These coherent structures are observed at the peak of the energy dissipation rate and are the carriers of long-range correlations despite all of the nonlinear interactions during the formation of turbulence. The size of these structures spans our computational domain, dominating the scaling of the energy spectrum, which follows a E∝k {sup –2} power law. As the Reynolds number increases, the curling of the current sheets due to Kelvin-Helmholtz-type instabilities and reconnection modifies the scaling of the energy spectrum from k {sup –2} toward k {sup –5/3}. This transition occurs due to the decorrelation of the velocity and the current density which is proportional to Re{sub λ}{sup −3/2}. Finite Reynolds number behavior is observed without reaching a finite asymptote for the energy dissipation rate even for a simulation of Re{sub λ} ≅ 440 with 2048{sup 3} grid points. This behavior demonstrates that even state-of-the-art numerical simulations of the highest Reynolds numbers can be influenced by the choice of initial conditions and consequently they are inadequate to deduce unequivocally the fate of universality in MHD turbulence. Implications for astrophysical observations are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1997

On the Numerical Prediction of Flow Transition

Knowledge from dynamical systems theory is used to study numerical uncertainties in direct numerical simulation of transition from laminar to turbulent flows. Spurious behavior due to underresolved grids and temporal discretizations is illustrated with 2-D CFD simulations.

Yee, H. C.↗
DOE OSTI2023

Faster-than-real-time Simulation with Demonstration for Resilient DER Integration

The US electric grid is facing operational, stability, and security challenges. Transmission system operators need some measure of visibility into distribution system renewable generation. Distribution system generation needs to support transmission system voltage. The grid is experiencing an expansion in measurement systems. How to take full advantage of this expansion and defend against attacks, both cyber and physical, poses additional challenges. The Faster-than-real-time Simulation with demonstration for Resilient DER Integration project set out to do the following: a. Flatten the voltage profile through the feeders and system for cost saving and voltage stabilization needs. b. Increase the amount of intermittent distributed energy resources (IDERs) that could be deployed on a utility feeder and provide 100% or more energy needed for the demands on that feeder, and c. based on an accurate model (Digital Twin) of the utilities system, be able to detect any abnormalities on the utilities distribution system. To manage the voltage and increase IDER penetration (a,b), Graph Trace Analysis is employed in a time-series, optimal power flow to coordinate the time-varying feedback control setpoints of a distribution feeder’s utility control devices. Under the coordinated control are a Load Tap Changing Transformer, a voltage regulator, and five switched capacitor banks. The feeder serves over 2000 customers, the feeder secondaries are modeled, and the feeder has 2.3 MW of PV generation, corresponding to a 17.4% penetration of PV generation. The feeder model has over 12,000 components, where every customer load bus and PV generator are modeled. The accuracy of the power flow solution is compared against historical meter voltage measurements, the improvement in conservation voltage reduction energy savings as a function of the coordinated control desired voltage profile is investigated, and the increase in PV penetration of the coordinated control over the existing control is presented. To achieve improved control performance while observing system operation constraints, bellwether Advanced Metering Infrastructure (AMI) voltage measurements are used to adjust the desired voltage profile used by the optimal power flow analysis. To detect and alleviate or negate attacks or failures on the distribution and transmission utility grids (c) the grid needs to be resilient and self-healing. In this project software was designed to do just that. At the center of the software is an Integrated System Model (ISM) that spans from transmission to secondary distribution. The ISM is employed in real-time abnormality detection, voltage stability forecasting, and multi-mode control. Testing results are presented for: 1—attacks on utility infrastructure; 2—energy savings from optimal control; 3—distribution system control response during a low voltage transmission system event; 4—cyber-attacks on PV inverters, where physical inverters are used in hard-ware-in-the-simulation-loop studies. Contributions of this work include real-time analysis that spans from three-phase transmission through secondary distribution; an approach for detecting abnormalities that employs measurements from three independent measurement systems; and a multi-mode distribution system control that responds to cyber-attacks, physical attacks, equipment failures, and transmission system needs.

Integrated System Model, Graph Trace Analysis, Adv↗
NASA NTRS1995

Interstellar clouds in high-speed, supersonic flows: Two-dimensional simulations

We present a series of gasdynamical simulations of the interaction of a dense, cool interstellar cloud with a high-speed, supersonic wind that confines and accelerates the embedded cloud. Our goal is to attempt to determine if such clouds can survive various potentially disruptive instabilities, that occur at their peripheries, long enough to be accelerated to speeds which are comparable to the wind velocity. These simulations are performed using two-dimensional, Eulerian gas dynamics on both an axisymmetric (about the cloud axis) and 'slab' geometric grid. The spatial and temporal resolutions of the simulations are varied over a wide range to investigate the effects of small-scale instabilities on the overall acceleration of clouds and the development of large-scale, disruptive instabilities. Also, we study the effects of wind/cloud Mach number variations by changing the wind speed constant at about 12 km/s (which corresponds to a cloud temperature of 10,000 K). The current simulations track the evolution of clouds as they are accelerated to speeds approximately 4-5 times greater than their internal sound speeds. Furthermore, the models with the highest resolution were extended far beyond quasi-linear Rayleigh-Taylor growth times reaching 6-7 Rayleigh-Taylor growth times for the largest scale instabilities before being terminated because of the accumulation of errors at the rear grid boundary.

Schiano, A. V. R.↗
DOE OSTI2023

Architecture for Co-Simulation of Transportation and Distribution Systems with Electric Vehicle Charging at Scale in the San Francisco Bay Area

This work describes the Grid-Enhanced, Mobility-Integrated Network Infrastructures for Extreme Fast Charging (GEMINI) architecture for the co-simulation of distribution and transportation systems to evaluate EV charging impacts on electric distribution systems of a large metropolitan area and the surrounding rural regions with high fidelity. The current co-simulation is applied to Oakland and Alameda, California, and in future work will be extended to the full San Francisco Bay Area. It uses the HELICS co-simulation framework to enable parallel instances of vetted grid and transportation software programs to interact at every model timestep, allowing high-fidelity simulations at a large scale. This enables not only the impacts of electrified transportation systems across a larger interconnected collection of distribution feeders to be evaluated, but also the feedbacks between the two systems, such as through control systems, to be captured and compared. The findings are that with moderate passenger EV adoption rates, inverter controls combined with some distribution system hardware upgrades can maintain grid voltages within ANSI C.84 range A limits of 0.95 to 1.05 p.u. without smart charging. However, EV charging control may be required for higher levels of charging or to reduce grid upgrades, and this will be explored in future work.

24 POWER TRANSMISSION AND DISTRIBUTION↗