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,441 records · Page 80

LAVA Simulations for the 3rd AIAA CFD High Lift Prediction Workshop with Body Fitted Grids

In response to the 3rd AIAA CFD High Lift Prediction Workshop, the workshop cases were analyzed using Reynolds-averaged Navier-Stokes flow solvers within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework. For the workshop cases the advantages and limitations of both overset-structured an unstructured polyhedral meshes were assessed. The workshop included 3 cases: a 2D airfoil validation case, a mesh convergence study using the High Lift Common Research Model, and a nacelle/pylon integration study using the JAXA (Japan Aerospace Exploration Agency) Standard Model. The 2D airfoil case from the workshop is used to verify the implementation of the Spalart-Allmaras turbulence model along with some of its variants within the solver. The High Lift Common Research Model case is used to assess solver performance and accuracy at varying mesh resolutions, as well as identify the minimum mesh fidelity required for LAVA on this class of problem. The JAXA Standard Model case is used to assess the solver's sensitivity to the turbulence model and to compare the structured and unstructured mesh paradigms. These workshop cases have helped establish best practices for high lift flow configurations for the LAVA solver.

LAVA↗

Uncertain Representations of Sub-Grid Pollutant Transport in Chemistry-Transport Models and Impacts on Long-Range Transport and Global Composition

Sub-grid transport, by convection and turbulence, is known to play an important role in lofting pollutants from their source regions. Consequently, the long-range transport and climatology of simulated atmospheric composition are impacted. This study uses the Goddard Earth Observing System, Version 5 (GEOS-5) atmospheric model to study pollutant transport. The baseline model uses a Relaxed Arakawa-Schubert (RAS) scheme that represents convection through a sequence of linearly entraining cloud plumes characterized by unique detrainment levels. Thermodynamics, moisture and trace gases are transported in the same manner. Various approximate forms of trace-gas transport are implemented, in which the box-averaged cloud mass fluxes from RAS are used with different numerical approaches. Substantial impacts on forward-model simulations of CO (using a linearized chemistry) are evident. In particular, some aspects of simulations using a diffusive form of sub-grid transport bear more resemblance to space-biased CO observations than do the baseline simulations with RAS transport. Implications for transport in the real atmosphere will be discussed. Another issue of importance is that many adjoint/inversion computations use simplified representations of sub-grid transport that may be inconsistent with the forward models: implications will be discussed. Finally, simulations using a complex chemistry model in GEOS-5 (in place of the linearized CO model) are underway: noteworthy results from this simulation will be mentioned.

Pawson, Steven↗

Simulated rarefied aerodynamics of the Magellan spacecraft during aerobraking

Aerodynamic loads upon the Magellan spacecraft during aerobraking through the atmosphere of Venus are computed at off-design attitudes with a direct simulation Monte Carlo (DSMC) particle method. Simulated rarefied flows at nominal altitudes near 140 km and an entry speed of 8.6 km/s were compared to simulated and analytic free molecular results. Aerodynamic moments, forces, and heating for rarefied entry at all attitudes were 7-10 percent below free molecular results. All moments acted to restore the vehicle to its nominal zero-pitch, zero-yaw attitude. Suggested canting of the solar panels is an innovative configuration to assess gas-surface interaction during aerobraking. The resulting roll torques about the central body-axis as predicted in rarefied flow simulations were nearly twice that predicted for free molecular flow, although differences became less distinct for thermal accommodation coefficients well below unity. Roll torques increased dramatically with reduced accommodation coefficients employed in the simulation. In the DSMC code, periodic free-molecule boundary conditions and a coarse computational grid and body resolution served to minimize the simulation size and cost while retaining solution validity.

Haas, Brian L.↗

Vlasov simulations of electron-ion collision effects on damping of electron plasma waves

Collisional effects can play an essential role in the dynamics of plasma waves by setting a minimum damping rate and by interfering with wave-particle resonances. Kinetic simulations of the effects of electron-ion pitch angle scattering on Electron Plasma Waves (EPWs) are presented here. In particular, the effects of such collisions on the frequency and damping of small-amplitude EPWs for a range of collision rates and wave phase velocities are computed and compared with theory. Both the Vlasov simulations and linear kinetic theory find the direct contribution of electron-ion collisions to wave damping significantly reduced from that obtained through linearized fluid theory. To our knowledge, this simple result has not been published before. Simulations have been carried out using a grid-based (Vlasov) approach, based on a high-order conservative finite difference method for discretizing the Fokker-Planck equation describing the evolution of the electron distribution function. Details of the implementation of the collision operator within this framework are presented. Such a grid-based approach, which is not subject to numerical noise, is of particular interest for the accurate measurements of the wave damping rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Potential of Data Center Controls in Grid Services

The rapid proliferation of large data centers brings both challenges and opportunities for grid reliability. The data center resources and their potential flexibility have the potential to contribute resources to grid operations. Through capabilities like energy shifting and resource coordination, data centers can help reduce their net demand on the transmission network, as well as provide additional grid services to support reliable operation on the grid. While transient and long-term grid planning and operations are the scenarios that draw most attention, the quasi-steady state timeseries (QSTS) operation of data centers and grid bring interesting scenarios that can help evaluate the data center controls to aid grid services. This work is focused on modeling data centers for QSTS applications – incorporating the AI data center load profiles and building on the PNNL digital twin model for the thermal management loads to enable simulation studies to reveal the impact of data center controls on grid performance. This includes the integration of a QSTS battery and natural gas generator model to incorporate local resource impacts to the system. The simulation study is performed with a modified IEEE 24-Bus transmission system. Scenarios are focused on evaluating the data center load impacts on the transmission system and leveraging both data center and local generation controls to mitigate those impacts and provide additional grid services. The data center controls revealed the ability to contribute to two main kinds of grid services: preventing congestion on a weak grid by coordinating the data center resources with the collocated BESS and onsite generation; and the ability to help the grid operations during stressed times of operation like during a contingency. Leveraging these and other capabilities has the potential to help data centers become grid responsive assets, aiding in both their integration into the power system and grid reliability.

power grid simulation↗

Plasma particle simulation of electrostatic ion thrusters

Charge exchange collisons between beam ions and neutral propellant gas can result in erosion of the accelerator grid surfaces of an ion engine. A particle in cell (PIC) is developed along with a Monte Carlo method to simulate the ion dynamics and charge exchange processes in the grid region of an ion thruster. The simulation is two-dimensional axisymmetric and uses three velocity components (2d3v) to investigate the influence of charge exchange collisions on the ion sputtering of the accelerator grid surfaces. An example calculation has been performed for an ion thruster operated on xenon propellant. The simulation shows that the greatest sputtering occurs on the downstream surface of the grid, but some sputtering can also occur on the upstream surface as well as on the interior of the grid aperture.

Peng, Xiaohang↗

Impedance-Based Stability Analysis in Grid Interconnection Impact Study Owing to the Increased Adoption of Converter-Interfaced Generators

This study investigates the emerging harmonic stability concerns to be addressed by grid planners in generation interconnection studies, owing to the increased adoption of renewable energy resources connected to the grid via power electronic converters. The wideband and high-frequency electromagnetic transient (EMT) characteristics of these converter-interfaced generators (CIGs) and their interaction with the grid impedance are not accurately captured in the typical dynamic studies conducted by grid planners. This paper thus identifies the desired components to be studied and subsequently develops a practical process for integrating a new CIG into a grid with the existing CIGs. The steps of this process are as follows: the impedance equation of a CIG using its control dynamics and an interface filter to the grid, for example, an LCL filter (inductor-capacitor-inductor type), is developed; an equivalent impedance model including the existing CIGs nearby and the grid observed from the point of common coupling are derived; the system stability for credible operating scenarios is assessed. Detailed EMT simulations validate the accuracy of the impedance models and stability assessment for various connection scenarios. Here, by complementing the conventional EMT simulation studies, the proposed analytical approach enables grid planners to identify critical design parameters for seamlessly integrating a new CIG and ensuring the reliability of the grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A regularized deconvolution method for turbulent closure modeling in implicitly filtered large-eddy simulation

Turbulence closure models that are based on a regularized deconvolution method (RDM) are proposed for application to implicitly filtered Large-Eddy Simulation (LES). This method reconstructs sub-grid scale contributions by approximately recovering filtered quantities through an optimization approach. Physical principles on the boundedness and conservation of reconstructed quantities are preserved by constraining the optimization problem. Further, to account for effects of scales that are not resolved in implicitly filtered LES, a reconstruction method is proposed that combines sub-grid projection and energy-similarity to enrich unresolved scales. With this method, closure models are developed for turbulent scalar fluxes and turbulence–chemistry interaction that are intrinsically consistent with the LES formulation. These closure models are evaluated in application to a piloted turbulent jet flame with inhomogeneous inlet composition, in which combustion is represented by a two-scalar manifold model. The RDM-formulation is benchmarked against simulations obtained using a presumed PDF method and eddy diffusivity closure model. In partially-premixed flame regions, it is observed that RDM improves the prediction of turbulent mixing, and faster grid convergence is observed for results with RDM compared to simulations employing a presumed PDF method.

42 ENGINEERING↗

Simulation of the elementary evolution operator with the motional states of an ion in an anharmonic trap

Following a recent proposal of L. Wang and D. Babikov [J. Chem. Phys. 137, 064301 (2012)], we theoretically illustrate the possibility of using the motional states of a Cd{sup +} ion trapped in a slightly anharmonic potential to simulate the single-particle time-dependent Schrödinger equation. The simulated wave packet is discretized on a spatial grid and the grid points are mapped on the ion motional states which define the qubit network. The localization probability at each grid point is obtained from the population in the corresponding motional state. The quantum gate is the elementary evolution operator corresponding to the time-dependent Schrödinger equation of the simulated system. The corresponding matrix can be estimated by any numerical algorithm. The radio-frequency field which is able to drive this unitary transformation among the qubit states of the ion is obtained by multi-target optimal control theory. The ion is assumed to be cooled in the ground motional state, and the preliminary step consists in initializing the qubits with the amplitudes of the initial simulated wave packet. The time evolution of the localization probability at the grids points is then obtained by successive applications of the gate and reading out the motional state population. The gate field is always identical for a given simulated potential, only the field preparing the initial wave packet has to be optimized for different simulations. We check the stability of the simulation against decoherence due to fluctuating electric fields in the trap electrodes by applying dissipative Lindblad dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Performance Evaluation of a Two-Dimensional Flood Model on Heterogeneous High-Performance Computing Architectures

This paper describes the implementation of a two-dimensional hydrodynamic flood model with two different numerical schemes on heterogeneous high-performance computing architectures. Both schemes were able to solve the nonlinear hyperbolic shallow water equations using an explicit upwind first-order approach on finite differences and finite volumes, respectively, and were conducted using MPI and CUDA. Four different test cases were simulated on the Summit supercomputer at Oak Ridge National Laboratory. Both numerical schemes scaled up to 128 nodes (768 GPUs) with a maximum 98.2x speedup of over 1 GPU. The lowest run time for the 10 day Hurricane Harvey event simulation at 5 meter resolution (272 million grid cells) was 50 minutes. GPUDirect communication proved to be more convenient than the standard communication strategy. Both strong and weak scaling are shown.

Sharif, Md Bulbul↗

Feasibility Analysis of Converter-Interfaced Combined Heat and Power System

As a promising new design concept, the converter-interfaced combined heat and power (CHP) system is coupled to the bulk grid through a rectifier and a grid-ready tied inverter. Compared to the traditional directly-coupled CHP system, it removes the requirement for oversizing the CHP generator, limits the short-circuit contribution of the generator and simplifies the grid integration process of CHP system. This paper evaluates the economic benefits of this concept by calculating the annualized Return-on-Investment (ROI) and comparing it to the directly-coupled system. The economic analysis includes timeseries simulations to compute energy transactions with the bulk grid as well as sizing the equipment to calculate the capital and operational costs. Obtained results indicate that in majority of user cases evaluated, the converter-interfaced CHP systems can provide better ROI than directly-coupled systems. Given the additional technical benefits provided by inverter-based distributed energy resources (DERs), the proposed concept is proved to be technically viable and economically feasible.

03 NATURAL GAS↗

Binary black hole accretion from a circumbinary disk: Gas dynamics inside the central cavity

We present the results of two-dimensional (2D) hydrodynamical simulations of circumbinary disk accretion using the finite-volume code DISCO. This code solves the 2D viscous Navier-Stokes equations on a high-resolution moving mesh which shears with the fluid flow, greatly reducing advection errors in comparison with a fixed grid. We perform a series of simulations for binary mass ratios in the range 0.026 ≤ q ≤ 1.0, each lasting longer than a viscous time so that we reach a quasi-steady accretion state. In each case, we find that gas is efficiently stripped from the inner edge of the circumbinary disk and enters the cavity along accretion streams, which feed persistent 'mini disks' surrounding each black hole. We find that for q ≳ 0.1, the binary excites eccentricity in the inner region of the circumbinary disk, creating an overdense lump which gives rise to enhanced periodicity in the accretion rate. The dependence of the periodicity on mass ratio may provide a method for observationally inferring mass ratios from measurements of the accretion rate. We also find that for all mass ratios studied, the magnitude of the accretion onto the secondary is sufficient to drive the binary toward larger mass ratio. This suggests a mechanism for biasing mass-ratio distributions toward equal mass.

79 ASTRONOMY AND ASTROPHYSICS↗

Transonic flow simulation of prop-fan configurations

This paper presents the results of a study that used a three-dimensional transonic analysis method to analyze flow around a wing-mounted prop-fan configuration. The transonic analysis method employs a cell-oriented finite volume approach to solve the full potential equation in conservative form over a numerically generated surface-fitted grid. The propeller slipstream effects are simulated using a linearized transpiration boundary condition. The results of analyses of flow over clean wing/body, wing/body/nacelle, and wing/body nacelle with slipstream are compared to test data.

Samant, S. S.↗

High-Resolution Regional Wave Hindcast for the U.S. Alaska Coast

Regional long-term wave hindcast with state-of-art third-generation models is essential for characterization of regional wave resources. This report summarizes modeling efforts for the simulation of the wave climate along southern Alaska, which used an unstructured, nested-grid modeling approach that incorporates a global-regional nested grid using WAVEWATCH III ® and the high-resolution SWAN (Simulating Waves Nearshore) model. Wave resource and wave bulk parameters were simulated for a 32-year period from 1979 to 2010 and were subsequently validated with wave buoy data in the model domain. Error statistics of model skills were calculated at all buoy stations. Overall, model results, both resource and bulk parameters, match observations well. The model hindcast was able to reproduce the seasonal variation of the sea state with large waves that occur in the winter and early spring months when wind forcing is strong and the calm sea state during the summer when wind forcing is weak. The hindcast also considers the effect of sea ice on wave growth, which affects the southern Bering Sea. The hindcast results, including resource and bulk parameters, were archived at every model grid point for the entire model domain at three-hour intervals. In addition, spectral results were saved at National Data Buoy Center buoy stations and selected virtual buoy locations at hourly intervals. The nested-grid modeling framework employed in this study provides a powerful and efficient modeling approach for accurately simulating wave climate at regional and long-term temporal scales with sufficiently fine resolutions in the nearshore region.

54 ENVIRONMENTAL SCIENCES↗

Application of Rotor Disk Model in the Overflow CFD Code

To reduce the complexity and to reduce the cost of turnaround time of simulating unsteady flows of rotating blades with moving overset grids, amid-fidelity option is provided in the OVERFLOW computational fluid dynamics analysis to model rotors and propellers as infinitely thin disks. The effect of the flow through the disk is represented by the momentum source terms in the governing Navier-Stokes equations. The model is loosely coupled with the rotorcraft comprehensive analysis CAMRADII. A brief description of the coupling algorithm is provided. Simulations of a single rotor and multiple rotors are shown, including comparisons with measurements and with the full rotor where all blades are modeled, and unsteady simulation performed with moving grids for rotating geometry,

ARMD↗

Modeling aerosols and their interactions with shallow cumuli during the 2007 CHAPS field study

The Weather Research and Forecasting model coupled with chemistry (WRF-Chem) is used to simulate relationships between aerosols and clouds in the vicinity of Oklahoma City during the June 2007 Cumulus Humilis Aerosol Processing Study (CHAPS). The regional scale simulation completed using 2 km horizontal grid spacing evaluates four important relationships between aerosols and shallow cumulus clouds observed during CHAPS. First, the model reproduces the trends of higher nitrate volume fractions in cloud droplet residuals compared to interstitial non-activated aerosols, as measured using the Aerosol Mass Spectrometer. Comparing simulations with cloud chemistry turned on and off, we show that nitric acid vapor uptake by cloud droplets explains the higher nitrate content of cloud droplet residuals. Second, as documented using an offline code, both aerosol water and other inorganics (OIN), which are related to dust and crustal emissions, significantly affect predicted aerosol optical properties. Reducing the OIN content of wet aerosols by 50% significantly improves agreement of model predictions with measurements of aerosol optical properties. Third, the simulated hygroscopicity of aerosols is too high as compared to their hygroscopicity derived from cloud condensation nuclei and particle size distribution measurements, indicating uncertainties associated with simulating size-dependent chemical composition and treatment of aerosol mixing state within the model. Fourth, the model reasonably represents the observations of the first aerosol indirect effect where pollutants in the vicinity of Oklahoma City increase cloud droplet number concentrations and decrease the droplet effective radius. While previous studies have often focused on cloud-aerosol interactions in stratiform and deep convective clouds, this study highlights the ability of regional-scale models to represent some of the important aspects of cloud-aerosol interactions associated with fields of short-lived shallow cumuli.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Grid-Integrated Electric Mobility Model (GEM) v1.0

Transportation is the fastest-growing source of greenhouse gas (GHG) emissions and energy consumption globally. The convergence of shared mobility, vehicle automation, and electrification has the potential to drastically reduce transportation impacts, but requires careful integration with rapidly evolving electricity systems. We have developed the GEM Model (Grid-Integrated Electric Mobility) to examine these interactions with a U.S.-wide simulation framework encompassing private electric vehicles (EVs); shared automated EVs (SAEVs); charging infrastructure; controlled EV charging; and a grid economic dispatch model to simulate mobility futures exclusively using EVs. We find that an SAEV fleet 9% the size of today's active vehicles can satisfy trip demand with only 2.6 million chargers (0.2 per EV). Controlled EV charging can also reduce electricity demand variability, significantly reducing GHG emissions and decreasing solar curtailment by about one-third. While private EVs with uncontrolled charging would reduce GHG emissions by 53% compared to gasoline vehicles, SAEVs could achieve a 70% reduction.

Sheppard, Colin↗