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 OSTI2012

Chance Constrained Optimal Power Flow: Risk-Aware Network Control under Uncertainty

When uncontrollable resources fluctuate, Optimum Power Flow (OPF), routinely used by the electric power industry to re-dispatch hourly controllable generation (coal, gas and hydro plants) over control areas of transmission networks, can result in grid instability, and, potentially, cascading outages. This risk arises be cause OPF dispatch is computed without awareness of major uncertainty, in particular fluctuations in renewable output. As a result, grid operation under the OPF with renewable variability can lead to frequent conditions where power line flow ratings are significantly exceeded. Such a condition, which is borne by simulations of real grids, would likely resulting in automatic line tripping to protect lines from thermal stress, a risky and undesirable outcome which compromises stability. Smart grid goals include a commitment to large penetration of highly fluctuating renewables, thus calling to reconsider current practices, in particular the use of standard OPF. Our Chance Constrained (CC) OPF corrects the problem and mitigates dangerous renewable fluctuations with minimal changes in the current operational procedure. Assuming availability of a reliable wind forecast parameterizing the distribution function of the uncertain generation, our CC-OPF satisfies all the constraints with high probability while simultaneously minimizing the cost of economic re-dispatch. CC-OPF allows efficient implementation, e.g. solving a typical instance over the 2746-bus Polish network in 20s on a standard laptop.

Energy Planning, Policy, & Economy(29)↗
NASA NTRS2001

Impact of Lidar Wind Sounding on Mesoscale Forecast

An Observing System Simulation Experiment (OSSE) was conducted to study the impact of airborne lidar wind sounding on mesoscale weather forecast. A wind retrieval scheme, which interpolates wind data from a grid data system, simulates the retrieval of wind profile from a satellite lidar system. A mesoscale forecast system based on the PSU/NCAR MM5 model is developed and incorporated the assimilation of the retrieved line-of-sight wind. To avoid the "identical twin" problem, the NCEP reanalysis data is used as our reference "nature" atmosphere. The simulated space-based lidar wind observations were retrieved by interpolating the NCEP values to the observation locations. A modified dataset obtained by smoothing the NCEP dataset was used as the initial state whose forecast was sought to be improved by assimilating the retrieved lidar observations. Forecasts using wind profiles with various lidar instrument parameters has been conducted. The results show that to significantly improve the mesoscale forecast the satellite should fly near the storm center with large scanning radius. Increasing lidar firing rate also improves the forecast. Cloud cover and lack of aerosol degrade the quality of the lidar wind data and, subsequently, the forecast.

Miller, Timothy L.↗
NASA NTRS2000

Geometry Modeling and Grid Generation for "Icing Effects" and "Ice Accretion" Simulations on Airfoils

There are two distinct icing-related problems for airfoils that can be simulated. One is predicting the effects of ice on the aerodynamic performance of airfoils when ice geometry is known ("icing effects" study). The other is simulating ice accretion under specified icing conditions ("ice accretion" simulation). This paper will address development of two different software packages for two-dimensional geometry preparation and grid generation for both "icing effects" and "ice accretion" studies.

Choo, Yung↗
DOE OSTI2014

A Sparse Grid Method for Bayesian Uncertainty Quantification with Application to Large Eddy Simulation Turbulence Models

There is wide agreement that the accuracy of turbulence models suffer from their sensitivity with respect to physical input data, the uncertainties of user-elected parameters, as well as the model inadequacy. However, the application of Bayesian inference to systematically quantify the uncertainties in parameters, by means of exploring posterior probability density functions (PPDFs), has been hindered by the prohibitively daunting computational cost associated with the large number of model executions, in addition to daunting computation time per one turbulence simulation. In this effort, we perform in this paper an adaptive hierarchical sparse grid surrogate modeling approach to Bayesian inference of large eddy simulation (LES). First, an adaptive hierarchical sparse grid surrogate for the output of forward models is constructed using a relatively small number of model executions. Using such surrogate, the likelihood function can be rapidly evaluated at any point in the parameter space without simulating the computationally expensive LES model. This method is essentially similar to those developed in Zhang et al. (Water Resour Res 49:6871–6892, 2013) for geophysical and groundwater models, but is adjusted and applied here for a much more challenging problem of uncertainty quantification of turbulence models. Through a numerical demonstration of the Smagorinsky model of two-dimensional flow around a cylinder at sub-critical Reynolds number, our approach is proven to significantly reduce the number of costly LES executions without losing much accuracy in the posterior probability estimation. Here, the model parameters are calibrated against synthetic data related to the mean flow velocity and Reynolds stresses at different locations in the flow wake. The influence of the user-elected LES parameters on the quality of output data will be discussed.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2004

Propulsion Simulations Using Advanced Turbulence Models with the Unstructured Grid CFD Tool, TetrUSS

A computational investigation has been completed to assess the capability of TetrUSS for exhaust nozzle flows. Three configurations were chosen for this study (1) an axisymmetric supersonic jet, (2) a transonic axisymmetric boattail with solid sting operated at different Reynolds number and Mach number, and (3) an isolated non-axisymmetric nacelle with a supersonic cruise nozzle. These configurations were chosen because existing experimental data provided a means for measuring the ability of TetrUSS for simulating complex nozzle flows. The main objective of this paper is to validate the implementation of advanced two-equation turbulence models in the unstructured-grid CFD code USM3D for propulsion flow cases. USM3D is the flow solver of the TetrUSS system. Three different turbulence models, namely, Menter Shear Stress Transport (SST), basic k epsilon, and the Spalart-Allmaras (SA) are used in the present study. The results are generally in agreement with other implementations of these models in structured-grid CFD codes. Results indicate that USM3D provides accurate simulations for complex aerodynamic configurations with propulsion integration.

Abdol-Hamid, Khaled S.↗
DOE OSTI2012

Effect of Terrestrial and Marine Organic Aerosol on Regional and Global Climate: Model Development, Application, and Verification with Satellite Data

In this DOE project the improvements to parameterization of marine primary organic matter (POM) emissions, hygroscopic properties of marine POM, marine isoprene derived secondary organic aerosol (SOA) emissions, surfactant effects, new cloud droplet activation parameterization have been implemented into Community Atmosphere Model (CAM 5.0), with a seven mode aerosol module from the Pacific Northwest National Laboratory (PNNL)'s Modal Aerosol Model (MAM7). The effects of marine aerosols derived from sea spray and ocean emitted biogenic volatile organic compounds (BVOCs) on microphysical properties of clouds were explored by conducting 10 year CAM5.0-MAM7 model simulations at a grid resolution 1.9° by 2.5° with 30 vertical layers. Model-predicted relationship between ocean physical and biological systems and the abundance of CCN in remote marine atmosphere was compared to data from the A-Train satellites (MODIS, CALIPSO, AMSR-E). Model simulations show that on average, primary and secondary organic aerosol emissions from the ocean can yield up to 20% increase in Cloud Condensation Nuclei (CCN) at 0.2% Supersaturation, and up to 5% increases in droplet number concentration of global maritime shallow clouds. Marine organics were treated as internally or externally mixed with sea salt. Changes associated with cloud properties reduced (absolute value) the model-predicted short wave cloud forcing from -1.35 Wm -2 to -0.25 Wm -2 . By using different emission scenarios, and droplet activation parameterizations, this study suggests that addition of marine primary aerosols and biologically generated reactive gases makes an important difference in radiative forcing assessments. All baseline and sensitivity simulations for 2001 and 2050 using global-through-urban WRF/Chem (GU-WRF) were completed. The main objective of these simulations was to evaluate the capability of GU-WRF for an accurate representation of the global atmosphere by exploring the most accurate configuration of physics options in GWRF for global scale modeling in 2001 at a horizontal grid resolution of 1° x 1°. GU-WRF model output was evaluated using observational datasets from a variety of sources including surface based observations (NCDC and BSRN), model reanalysis (NCEP/ NCAR Reanalysis and CMAP), and remotely-sensed data (TRMM) to evaluate the ability of GU-WRF to simulate atmospheric variables at the surface as well as aloft. Explicit treatment of nanoparticles produced from new particle formation in GU-WRF/Chem-MADRID was achieved by expanding particle size sections from 8 to 12 to cover particles with the size range of 1.16 nm to 11.6m. Simulations with two different nucleation parameterizations were conducted for August 2002 over a global domain at a 4º by 5º horizontal resolution. The results are evaluated against field measurement data from the 2002 Aerosol Nucleation and Real Time Characterization Experiment (ANARChE) in Atlanta, Georgia, as well as satellite and reanalysis data. We have also explored the relationship between clean marine aerosol optical properties and ocean surface wind speed using remotely sensed data from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) on board the CALIPSO satellite and the Advanced Microwave Scanning Radiometer (AMSR-E) on board the AQUA satellite. Detailed data analyses were carried out over 15 regions selected to be representative of different areas of the global ocean for the time period from June 2006 to April 2011. We show that for very low (less than 4 m s -1 ) and very high (more than 12 m s -1 ) wind speed conditions the mean CALIPSO-derived aerosol optical depth (AOD) has little dependency on the surface wind speed. For an intermediate (between 4 and 12 m s -1 ) marine AOD was linearly correlated with the surface wind speed values, with a slope of 0.0062 s m -1 . Results of our study suggest that considerable improvements to both optical properties of marine aerosols and their production mechanisms can be achieved by discriminating clean marine aerosols (or sea salt particles) from all other types of aerosols present over the ocean.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

Simulation of wave propagation in three-dimensional random media

Quantitative error analysis for simulation of wave propagation in three dimensional random media assuming narrow angular scattering are presented for the plane wave and spherical wave geometry. This includes the errors resulting from finite grid size, finite simulation dimensions, and the separation of the two-dimensional screens along the propagation direction. Simple error scalings are determined for power-law spectra of the random refractive index of the media. The effects of a finite inner scale are also considered. The spatial spectra of the intensity errors are calculated and compared to the spatial spectra of intensity. The numerical requirements for a simulation of given accuracy are determined for realizations of the field. The numerical requirements for accurate estimation of higher moments of the field are less stringent.

Coles, William A.↗
NASA NTRS2011

Properties of a Variable-Delay Polarization Modulator

We investigate the polarization modulation properties of a variable-delay polarization modulator (VPM). The VPM modulates polarization via a variable separation between a polarizing grid and a parallel mirror. We find that in the limit where the wavelength is much larger than the diameter of the metal wires that comprise the grid, the phase delay derived from the geometric separation between the mirror and the grid is sufficient to characterize the device. However, outside of this range, additional parameters describing the polarizing grid geometry must be included to fully characterize the modulator response. In this paper, we report test results of a VPM at wavelengths of 350 micron and 3 mm. Electromagnetic simulations of wire grid polarizers were performed and are summarized using a simple circuit model that incorporates the loss and polarization properties of the device.

Chuss, David T.↗
NASA NTRS2010

Validation of a Three-Dimensional Ablation and Thermal Response Simulation Code

The 3dFIAT code simulates pyrolysis, ablation, and shape change of thermal protection materials and systems in three dimensions. The governing equations, which include energy conservation, a three-component decomposition model, and a surface energy balance, are solved with a moving grid system to simulate the shape change due to surface recession. This work is the first part of a code validation study for new capabilities that were added to 3dFIAT. These expanded capabilities include a multi-block moving grid system and an orthotropic thermal conductivity model. This paper focuses on conditions with minimal shape change in which the fluid/solid coupling is not necessary. Two groups of test cases of 3dFIAT analyses of Phenolic Impregnated Carbon Ablator in an arc-jet are presented. In the first group, axisymmetric iso-q shaped models are studied to check the accuracy of three-dimensional multi-block grid system. In the second group, similar models with various through-the-thickness conductivity directions are examined. In this group, the material thermal response is three-dimensional, because of the carbon fiber orientation. Predictions from 3dFIAT are presented and compared with arcjet test data. The 3dFIAT predictions agree very well with thermocouple data for both groups of test cases.

Chen, Yih-Kanq↗
DOE OSTI2022

Development of a Configurable DERMS Test System in GridAPPS-D

Distributed Energy Resources and their increased penetration within the power grid are a suitable method to provide essential reliability services. In order to do so, they must operate in aggregate. Aggregation of a wide variety of systems with varying standards and protocols is a complex task, which could benefit from rules and standards to support interoperability. To this end, the Portland State University Power Lab has developed the Energy Grid of Things project, which includes a Distributed Energy Resource Management System to provide aggregation and demonstrate the usage of an Energy Services Interface: a set of rules governing the interface between resource owners and aggregator businesses. This thesis presents a Modeling Environment that provides a simulation and test platform reflecting the complex needs of an aggregator. A grid simulator is required that is able to reflect the effects of thousands of Distributed Energy Resources operating in tandem; furthermore, these effects must be communicated to external actors such as the Grid Operator, which provides grid service requests to the aggregator and requires feedback. The Modeling Environment provides all of this functionality by leveraging the GridAPPS-D platform to develop an interoperable class-based system that provides configurable simulations with modular input and output APIs. The Modeling Environment uses abstractions to allow any feasible type of input to be reflected within the model as a Distributed Energy Resource, and serves as a proof-of-concept for the development of aggregation simulation and testing within GridAPPS-D.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Effect of Secondary Circulations on the Surface–Atmosphere Exchange of Energy at an Isolated Semi-arid Forest

Afforestation in semi-arid regions can potentially enhance the global carbon sink by increasing the terrestrial biomass. However, the survival of planted forests under such extreme environmental conditions is not guaranteed a priori, and critically depends on the surface–atmosphere exchange of energy. Here, we investigate the pine forest Yatir in Israel, an example of a man-made semi-arid ecosystem, by means of large-eddy simulations. We focus on the interaction between surface–atmosphere exchange and secondary circulations that couple the isolated forest to the surrounding shrubland. The large-eddy simulations feature a grid resolution that resolves the forest canopy in several layers, and are initialized by satellite data and Doppler lidar, eddy-covariance and radiosonde measurements. We perform three large-eddy simulations with different geostrophic wind speeds to investigate the influence of those wind speeds on the surface–atmosphere exchange. We reproduce the measured mean updrafts above the forest and mean downdrafts above the shrubland, which increase in strength with decreasing geostrophic wind speed. The largest updrafts emerge above the older, denser part of the forest, triggering secondary circulations. The spatial extent of these circulations does not cover the entire forest area, although we observe a reduced aerodynamic resistance in the regions of updraft. Our simulations indicate that the enhanced surface–atmosphere exchange of the Yatir forest is not sufficient to compensate for the increased net radiation, due to the lower albedo of the forest with respect to the surroundings, resulting in higher air temperatures inside the forest. However, the difference between the forest and shrubland temperatures decreases with increasing geostrophic wind speed due to reduction in the aerodynamic resistance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simulation of gas mixture dynamics in a pipeline network using explicit staggered-grid discretization

Here we develop an explicit staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge boundary values that are more complex than for a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Meso-Beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

The Control-B simulation experiment was designed to improve on the Control-A experiment performed with the GMASS model. This experiment addressed several inadequacies with the first smooth terrain numerical simulation by including: (1) increased nested-grid resolution to better define the simulated gravity waves, (2) increased horizontal diffusion to remove outflow boundary condition noise, and (3) the use of reanalyzed rawinsonde data and surface observations in the initial state to increase the definition of the observed jet streak as well as other low-level features. A smoothed-terrain dry simulation with the nested-grid GMASS model has revealed many important aspects of the processes which resulted in the generation of gravity waves in the region and time when and where they were observed. However, the vertical structure, number, and characteristics of the waves are still quite different from observed waves as diagnosed thus necessitating future improved simulations. However, this control simulation has produced substantial insight into processes which occur on many spatial scales over a 30 hour time period thus allowing one to draw promising inferences as to the mechanisms for the complex process which occurred in nature during the CCOPE case study. The theoretical aspects of the project have focus on understanding the nature of the ageostrophic circulations which are produced in idealized models of the atmosphere in which the troposphere is modeled in one of two ways. The first model assumes that the lower atmosphere can be represented as a single layer of homogeneous fluid whose upper surface is free to exhibit vertical displacement. Two-dimensional internal convergence (divergence) occurring during the adjustment to an asymptotic equilibrium state from an ageostrophic initial state whose momentum structure is representative of a midlatitude localized zonal wind anomaly will cause the free upper surface of the homogeneous atmosphere to rise (fall), and therefore the response can be viewed as being physically three-dimensional. The second model assumes that the troposphere can be represented by an unbounded continuously stratified Boussinesq fluid of constant Brunt-Vaisala frequency N, where the vertical gradient of the basic state potential temperature profile allows for the existence of vertically propagating internal inertia-gravity waves.

Lin, Yuh-Lang↗
DOE OSTI2017

Adaptive RAS/SPS System Settings for Improving Grid Reliability and Asset Utilization through Predictive Simulation and Controls

Remedial Action Schemes (RAS) or Special Protection Schemes (SPS) are used throughout the bulk transmission system as a non-wires method of increasing transmission transfer capability. However, as the population grows within the Western Interconnection, more power needs to be transferred with less ability to build lines. This leads to a greater number of RAS, and the subsequent interaction between these disparate schemes can have unintended consequences. While the large number of RAS within the Western Electricity Coordinating Council (WECC) protect for different things, which can range from the violation of thermal limitations to transient stability issues, they can be described by common methods of operation. In order to understand these methods, a survey was sent out to members of the WECC Modeling & Validation Working Group (MVWG). The response received provides a necessary base in which to understand and build models that can be used within WECC. This report provides information on the survey results regarding RAS/SPS modeling practices for planning and operational studies within the WECC footprint, which is foundational for on-going and future work on RAS. Pacific Northwest National Laboratory (PNNL) and its industry partners are developing innovative methods for adaptively setting RAS parameters based on realistic and near real-time operational conditions to improve power grid reliability and grid asset utilization.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

Features in simulation of crystal growth using the hyperbolic PFC equation and the dependence of the numerical solution on the parameters of the computational grid

We investigate the three-dimensional mathematical model of crystal growth called PFC (Phase Field Crystal) in a hyperbolic modification. This model is also called the modified model PFC (originally PFC model is formulated in parabolic form) and allows to describe both slow and rapid crystallization processes on atomic length scales and on diffusive time scales. Modified PFC model is described by the differential equation in partial derivatives of the sixth order in space and second order in time. The solution of this equation is possible only by numerical methods. Previously, authors created the software package for the solution of the Phase Field Crystal problem, based on the method of isogeometric analysis (IGA) and PetIGA program library. During further investigation it was found that the quality of the solution can strongly depends on the discretization parameters of a numerical method. In this report, we show the features that should be taken into account during constructing the computational grid for the numerical simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2004

Particle Acceleration and Magnetic Field Generation in Electron-Positron Relativistic Shocks

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., Buneman, Weibel and other two-stream instabilities) created in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic electron-positron jet front propagating into an ambient electron-positron plasma with and without initial magnetic fields. We find small differences in the results for no ambient and modest ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. The non-linear fluctuation amplitudes of densities, currents, electric, and magnetic fields in the electron-positron shock are larger than those found in the electron-ion shock studied in a previous paper. This comes from the fact that both electrons and positrons contribute to generation of the Weibel instability. Additionally, we have performed simulations with different electron skin depths. We find that growth times scale inversely with the plasma frequency, and the sizes of structures created by the Weibel instability scale proportional to the electron skin depth. This is the expected result and indicates that the simulations have sufficient grid resolution. While some Fermi acceleration may occur at the jet front, the majority of electron and positron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying nonuniform, small-scale magnetic fields which contribute to the electron's (positron's) transverse deflection behind the jet head. This small scale magnetic field structure is appropriate to the generation of "jitter" radiation from deflected electrons (positrons) as opposed to synchrotron radiation. The jitter radiation has different properties than synchrotron radiation calculated assuming a a uniform magnetic field. The jitter radiation resulting from small scale magnetic field structures may be important for understanding the complex time structure and spectral evolution observed in gamma-ray bursts or other astrophysical sources containing relativistic jets and relativistic collisionless shocks.

Nishikawa, K.-I.↗