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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
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 at the comparable simulation time. 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.-L.↗
NASA NTRS2005

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

Shock acceleration is a 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 three-dimensional 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. New 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. Furthermore, the nonlinear fluctuation amplitudes of densities, currents, and electric and magnetic fields in the electron-positron shock are larger than those found in the electron-ion shock studied in a previous paper at a comparable simulation time. This comes from the fact that both electrons and positrons contribute to generation of the Weibel instability. In addition, 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 tine Weibel instability scale proportionally 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 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.↗
DOE OSTI2019

Photovoltaic and Behind-the-Meter Battery Storage: Advanced Smart Inverter Controls and Field Demonstration

Electric utilities have little visibility of the electrical distribution system, and consequently, limited diagnostic capabilities. The distribution grid was designed for a unidirectional power flow, where energy is supplied by few large centralized power plants; however, this is set to change to meet California’s aggressive decarbonization goals. The large-scale deployment of distributed renewable generation, such as photovoltaics (PV), can have many negative effects on an unprepared grid. As illustrated in the California “duck curve,” PV generation modifies load profiles during the day, which causes steep ramping in the evening. This effect indicates the need to revise the electrical grid’s design and operation. This project sought to (a) create a centralized resource in California to test and validate distribution technology controls with industry standard PV, storage, and high-fidelity sensors, (b) support strategic and operational decisions for new grid architectures by providing simulation models, (c) promote new ways to control clusters of PV smart inverters, in accordance with California Rule 21, and (d) innovate, develop, and field test a predictive controller to maximize economic profit for the customer while supporting the grid. The controller was built using the state-of-the-art model predictive control methodology to optimally control behind-the-meter PV and battery storage. In consideration of the duck curve, the controller optimally controls the battery dispatch by charging during excess generation periods and discharging during the critical afternoon demand hours. The controller was evaluated in annual simulations and revealed the potential cost-effectiveness of behind-the-meter battery storage. The simulations showed that the annual electricity bill could be reduced by as much as 35 percent, with a payback period of the investment in battery storage in about 6 years – significantly shorter than the manufacturer’s 10-year warranty. All developed simulation models, the grid event library, and the MPC controller are open-source and available online.

Model Predictive control↗
NASA NTRS1996

A Mixed Finite Volume Element Method for Flow Calculations in Porous Media

A key ingredient in the simulation of flow in porous media is the accurate determination of the velocities that drive the flow. The large scale irregularities of the geology, such as faults, fractures, and layers suggest the use of irregular grids in the simulation. Work has been done in applying the finite volume element (FVE) methodology as developed by McCormick in conjunction with mixed methods which were developed by Raviart and Thomas. The resulting mixed finite volume element discretization scheme has the potential to generate more accurate solutions than standard approaches. The focus of this paper is on a multilevel algorithm for solving the discrete mixed FVE equations. The algorithm uses a standard cell centered finite difference scheme as the 'coarse' level and the more accurate mixed FVE scheme as the 'fine' level. The algorithm appears to have potential as a fast solver for large size simulations of flow in porous media.

Jim E Jones↗
NASA NTRS2018

Analysis of the Planetary Boundary Layer Height During DISCOVER-AQ Baltimore–Washington, D.C., with Lidar and High-Resolution WRF Modeling

The daytime planetary boundary layer (PBL) was examined for the Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality (DISCOVER-AQ) Baltimore (Maryland)–Washington, D.C., campaign of July 2011 using PBL height (PBLH) retrievals from aerosol backscatter measurements from ground-based micropulse lidar (MPL), the NASA Langley Research Center airborne High Spectral Resolution Lidar-1 (HSRL-1), and the Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP) on the Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite. High-resolution Weather Research and Forecasting (WRF) Model simulations with horizontal grid spacing of 1 km and different combinations of PBL schemes, urban parameterization, and sea surface temperature inputs were evaluated against PBLHs derived from lidars, ozonesondes, and radiosondes. MPL and WRF PBLHs depicted a growing PBL in the morning that reached a peak height by midafternoon. WRF PBLHs calculated from gridded output profiles generally showed more rapid growth and higher peak heights than did the MPLs, and all WRF–lidar differences were dependent on model configuration, PBLH calculation method, and synoptic conditions. At inland locations, WRF simulated an earlier descent of the PBL top in the afternoon relative to the MPL retrievals and radiosonde PBLHs. At Edgewood, Maryland, the influence of the Chesapeake Bay breeze on the PBLH was captured by both the ozonesonde and WRF data but generally not by the MPL PBLH retrievals because of generally weaker gradients in the aerosol backscatter profile and limited normalized relative backscatter data near the top height of the marine layer.

Hegarty, Jennifer D.↗
DOE OSTI2022

Enhancing Local Grid Resilience with Small Hydropower Hybrids: Proving the concept through demonstration, simulation, and analysis with Idaho Falls Power

Large hydropower, connected to the transmission system and typically possessing significant ability to balance grid frequency, has long been central to black start plans in regions where it is present. Small hydropower possesses most of the attributes required for black start but is often connected to distribution or sub-transmission systems and has less ability to balance grid frequency. Integrating energy storage such as batteries or ultracapacitors increases the combined asset’s ability to balance frequency. This asset integration enables a bottom-up grid restoration paradigm in which critical electric loads on the local distribution system can be powered even when the regional transmission system is down. This report presents a field demonstration conducted with Idaho Falls Power that proved this concept. The contribution of energy storage in restoring the small hydro-dominated distribution grid and its operational sensitivity across different control settings are further analyzed using high-fidelity simulations. Readers will learn about technical details on the field demonstration setup, energy storage contribution to black start, detailed simulation steps of islanded distribution grid restoration, and usage of field demonstration measurements for transient model refinements. Collectively this report points to a great opportunity for small hydropower to enhance resilience of local electric grids.

13 HYDRO ENERGY↗
DOE OSTI2017

Assessment of Vertical Mesh Refinement in Concurrently Nested Large-Eddy Simulations Using the Weather Research and Forecasting Model

To facilitate multiscale simulation using the Weather Research and Forecasting Model, vertical mesh refinement for one-way concurrent nested simulation was recently introduced. Grid refinement in the vertical dimension removes the requirement of different grid aspect ratios on the bounding versus the nested domain, such that results from refinement are in the horizontal directions only, and thereby can also reduce numerical errors on the bounding domain arising from large aspect ratios in the presence of complex terrain. Herein, the impacts of vertical grid refinement on the evolving downstream flow in nested large-eddy simulations are evaluated in relation to other model configuration choices, including turbulence subfilter-scale (SFS) stress models, mesh configuration, and alternative methods for calculating several near-surface flow parameters. Although vertical nesting requires coarsening of the vertical grid on the bounding domain, leading to a smaller range of resolved turbulence scales in the nest’s lateral boundary conditions, parameter values within the nested domains are generally only minimally impacted, relative to nesting using the same vertical grid on each domain. Two dynamic SFS models examined herein generally improved the simulated mean wind speed, turbulence kinetic energy, stresses and spectra, on both domains, and accelerated equilibration rates within nested domains, relative to two constant coefficient models. A new method of extrapolating horizontal velocity components to near-surface locations at nested domain lateral boundaries, and a correction to the calculation of deformation elements near the surface, are each shown to slightly alter the mean parameter values, yet only minimally impact equilibration rates within the nested domain.

14 SOLAR ENERGY↗
NASA NTRS2014

Uncertainties in the Modelled CO2 Threshold for Antarctic Glaciation

frequently cited atmospheric CO2 threshold for the onset of Antarctic glaciation of approximately780 parts per million by volume is based on the study of DeConto and Pollard (2003) using an ice sheet model and the GENESIS climate model. Proxy records suggest that atmospheric CO2 concentrations passed through this threshold across the Eocene-Oligocene transition approximately 34 million years. However, atmospheric CO2 concentrations may have been close to this threshold earlier than this transition, which is used by some to suggest the possibility of Antarctic ice sheets during the Eocene. Here we investigate the climate model dependency of the threshold for Antarctic glaciation by performing offline ice sheet model simulations using the climate from 7 different climate models with Eocene boundary conditions (HadCM3L, CCSM3, CESM1.0, GENESIS, FAMOUS, ECHAM5 and GISS_ER). These climate simulations are sourced from a number of independent studies, and as such the boundary conditions, which are poorly constrained during the Eocene, are not identical between simulations. The results of this study suggest that the atmospheric CO2 threshold for Antarctic glaciation is highly dependent on the climate model used and the climate model configuration. A large discrepancy between the climate model and ice sheet model grids for some simulations leads to a strong sensitivity to the lapse rate parameter.

simulation↗
DOE OSTI2012

CENTRAL REGIONS OF BARRED GALAXIES: TWO-DIMENSIONAL NON-SELF-GRAVITATING HYDRODYNAMIC SIMULATIONS

The inner regions of barred galaxies contain substructures such as off-axis shocks, nuclear rings, and nuclear spirals. These substructures may affect star formation, and control the activity of a central black hole (BH) by determining the mass inflow rate. We investigate the formation and properties of such substructures using high-resolution, grid-based hydrodynamic simulations. The gaseous medium is assumed to be infinitesimally thin, isothermal, and non-self-gravitating. The stars and dark matter are represented by a static gravitational potential with four components: a stellar disk, a bulge, a central BH, and a bar. To investigate various galactic environments, we vary the gas sound speed, c{sub s} , as well as the mass of the central BH, M{sub BH}. Once the flow has reached a quasi-steady state, off-axis shocks tend to move closer to the bar major axis as c{sub s} increases. Nuclear rings shrink in size with increasing c{sub s} , but are independent of M{sub BH}, suggesting that the ring position is not determined by the Lindblad resonances. Rings in low-c{sub s} models are narrow since they are occupied largely by gas on x{sub 2}-orbits and well decoupled from nuclear spirals, while they become broad because of large thermal perturbations in high-c{sub s} models. Nuclear spirals persist only when either c{sub s} is small or M{sub BH} is large; they would otherwise be destroyed completely by the ring material on eccentric orbits. The shape and strength of nuclear spirals depend sensitively on c{sub s} and M{sub BH} such that they are leading if both c{sub s} and M{sub BH} are small, weak trailing if c{sub s} is small and M{sub BH} is large, and strong trailing if both c{sub s} and M{sub BH} are large. While the mass inflow rate toward the nucleus is quite small in low-c{sub s} models because of the presence of a narrow nuclear ring, it becomes larger than 0.01 M{sub Sun} yr{sup -1} when c{sub s} is large, providing a potential explanation of nuclear activity in Seyfert galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Phase transitions of macromolecular microsphere composite hydrogels based on the stochastic Cahn–Hilliard equation

We use the stochastic Cahn–Hilliard equation to simulate the phase transitions of the macromolecular microsphere composite (MMC) hydrogels under a random disturbance. Based on the Flory–Huggins lattice model and the Boltzmann entropy theorem, we develop a reticular free energy suit for the network structure of MMC hydrogels. Taking the random factor into account, with the time-dependent Ginzburg-Landau (TDGL) mesoscopic simulation method, we set up a stochastic Cahn–Hilliard equation, designated herein as the MMC-TDGL equation. The stochastic term in the equation is constructed appropriately to satisfy the fluctuation-dissipation theorem and is discretized on a spatial grid for the simulation. A semi-implicit difference scheme is adopted to numerically solve the MMC-TDGL equation. Some numerical experiments are performed with different parameters. The results are consistent with the physical phenomenon, which verifies the good simulation of the stochastic term.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Incorporation of aerosol into the COSPv2 satellite lidar simulator for climate model evaluation

Abstract. Atmospheric aerosol has substantial impacts on climate, air quality and biogeochemical cycles, and its concentrations are highly variable in space and time. A key variability to evaluate within models that simulate aerosol is the vertical distribution, which influences atmospheric heating profiles and aerosol–cloud interactions, to help constrain aerosol residence time and to better represent the magnitude of simulated impacts. To ensure a consistent comparison between modeled and observed vertical distribution of aerosol, we implemented an aerosol lidar simulator within the Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package version 2 (COSPv2). We assessed the attenuated total backscattered (ATB) signal and the backscatter ratios (SRs) at 532 nm in the U.S. Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1). The simulator performs the computations at the same vertical resolution as the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP), making use of aerosol optics from the E3SMv1 model as inputs and assuming that aerosol is uniformly distributed horizontally within each model grid box. The simulator applies a cloud masking and an aerosol detection threshold to obtain the ATB and SR profiles that would be observed above clouds by CALIOP with its aerosol detection capability. Our analysis shows that the aerosol distribution simulated at a seasonal timescale is generally in good agreement with observations. Over the Southern Ocean, however, the model does not produce the SR maximum as observed in the real world. Comparison between clear-sky and all-sky SRs shows little differences, indicating that the cloud screening by potentially incorrect model clouds does not affect the mean aerosol signal averaged over a season. This indicates that the differences between observed and simulated SR values are due not to sampling errors, but to deficiencies in the representation of aerosol in models. Finally, we highlight the need for future applications of lidar observations at multiple wavelengths to provide insights into aerosol properties and distribution and their representation in Earth system models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2003

An Aerodynamic Analysis of a Spinning Missile with Dithering Canards

A generic spinning missile with dithering canards is used to demonstrate the utility of an overset structured grid approach for simulating the aerodynamics of rolling airframe missile systems. The approach is used to generate a modest aerodynamic database for the generic missile. The database is populated with solutions to the Euler and Navier-Stokes equations. It is used to evaluate grid resolution requirements for accurate prediction of instantaneous missile loads and the relative aerodynamic significance of angle-of-attack, canard pitching sequence, viscous effects, and roll-rate effects. A novel analytical method for inter- and extrapolation of database results is also given.

Meakin, Robert L.↗
NASA NTRS1994

Towards a Comprehensive Computational Simulation System for Turbomachinery

The objective of this work is to develop algorithms associated with a comprehensive computational simulation system for turbomachinery flow fields. This development is accomplished in a modular fashion. These modules includes grid generation, visualization, network, simulation, toolbox, and flow modules. An interactive grid generation module is customized to facilitate the grid generation process associated with complicated turbomachinery configurations. With its user-friendly graphical user interface, the user may interactively manipulate the default settings to obtain a quality grid within a fraction of time that is usually required for building a grid about the same geometry with a general-purpose grid generation code. Non-Uniform Rational B-Spline formulations are utilized in the algorithm to maintain geometry fidelity while redistributing grid points on the solid surfaces. Bezier curve formulation is used to allow interactive construction of inner boundaries. It is also utilized to allow interactive point distribution. Cascade surfaces are transformed from three-dimensional surfaces of revolution into two-dimensional parametric planes for easy manipulation. Such a transformation allows these manipulated plane grids to be mapped to surfaces of revolution by any generatrix definition. A sophisticated visualization module is developed to al-low visualization for both grid and flow solution, steady or unsteady. A network module is built to allow data transferring in the heterogeneous environment. A flow module is integrated into this system, using an existing turbomachinery flow code. A simulation module is developed to combine the network, flow, and visualization module to achieve near real-time flow simulation about turbomachinery geometries. A toolbox module is developed to support the overall task. A batch version of the grid generation module is developed to allow portability and has been extended to allow dynamic grid generation for pitch changing turbomachinery configurations. Various applications with different characteristics are presented to demonstrate the success of this system.

Shih, Ming-Hsin↗
DOE OSTI2022

Methods for estimating X-ray machine output through measurement and simulation

We report ball grid arrays are increasingly being applied in the electronics industry and may require X-ray inspection to ensure the integrity and correct placement of solder pins. However, as the architecture of integrated circuits continues to narrow while simultaneously growing more complex, the risk of electronic failure due to radiation damage increases. While medical X-ray devices have been held to high standards and are repeatedly shown to be well characterized, devices used for electronic inspection are often lacking detailed characterization. This study presents unique methods to solve for important properties in X-ray inspection devices such as source to object distance and energy spectrum. This information can then be applied to Monte Carlo models to achieve better overall dose estimates to electronics, which will lead to superior manufactured products. Since X-ray devices can vary greatly in source characteristics, this work investigates spectral measurement and Monte Carlo representation of three X-ray devices. For a Philips SRO 33 100 medical diagnostic device, the spectral output followed expected trends given by the prediction software SpekCalc and Spektr. For the Dage XD7500NT, direct measurement showed a spectral artifact that through the use of Gafchromic films, was shown to be a contributing effect in the dose output. For the Rad Source RS1800, a high powered irradiation device, direct spectral measurement was not achieved. However, a Monte Carlo model using an assumed spectra was found to match ion chamber measurements to a high degree.

47 OTHER INSTRUMENTATION↗
NASA NTRS1980

An electrostatic analog for generating cascade grids

Accurate and efficient numerical simulation of flows through turbomachinery blade rows depends on the topology of the computational grids. These grids must reflect the periodic nature of turbomachinery blade row geometries and conform to the blade shapes. Three types of grids can be generated that meet these minimal requirements: through-flow grids, O-type grids, and C-type grids. A procedure which can be used to generate all three types of grids is presented. The resulting grids are orthogonal and can be stretched to capture the essential physics of the flow. A discussion is also presented detailing the extension of the generation procedure to three dimensional geometries.

Adamczyk, J. J.↗