Results for “grid simulator” — Page 32

NASA NTRS1999

Regional Climate Simulation with a Variable Resolution Stretched Grid GCM: The Regional Down-Scaling Effects

The results obtained with the variable resolution stretched grid (SG) GEOS GCM (Goddard Earth Observing System General Circulation Models) are discussed, with the emphasis on the regional down-scaling effects and their dependence on the stretched grid design and parameters. A variable resolution SG-GCM and SG-DAS using a global stretched grid with fine resolution over an area of interest, is a viable new approach to REGIONAL and subregional CLIMATE studies and applications. The stretched grid approach is an ideal tool for representing regional to global scale interactions. It is an alternative to the widely used nested grid approach introduced a decade ago as a pioneering step in regional climate modeling. The GEOS SG-GCM is used for simulations of the anomalous U.S. climate events of 1988 drought and 1993 flood, with enhanced regional resolution. The height low level jet, precipitation and other diagnostic patterns are successfully simulated and show the efficient down-scaling over the area of interest the U.S. An imitation of the nested grid approach is performed using the developed SG-DAS (Data Assimilation System) that incorporates the SG-GCM. The SG-DAS is run with withholding data over the area of interest. The design immitates the nested grid framework with boundary conditions provided from analyses. No boundary condition buffer is needed for the case due to the global domain of integration used for the SG-GCM and SG-DAS. The experiments based on the newly developed versions of the GEOS SG-GCM and SG-DAS, with finer 0.5 degree (and higher) regional resolution, are briefly discussed. The major aspects of parallelization of the SG-GCM code are outlined. The KEY OBJECTIVES of the study are: 1) obtaining an efficient DOWN-SCALING over the area of interest with fine and very fine resolution; 2) providing CONSISTENT interactions between regional and global scales including the consistent representation of regional ENERGY and WATER BALANCES; 3) providing a high computational efficiency for future SG-GCM and SG-DAS versions using PARALLEL codes.

Fox-Rabinovitz, Michael S.↗

DOE OSTI2012

Is there a single spot size and grid for intensity modulated proton therapy? Simulation of head and neck, prostate and mesothelioma cases

Purpose: To assess the quality of dose distributions in real clinical cases for different dimensions of scanned proton pencil beams. The distance between spots (i.e., the grid of delivery) is optimized for each dimension of the pencil beam. Methods: The authors vary the {sigma} of the initial Gaussian size of the spot, from {sigma}{sub x} = {sigma}{sub y} = 3 mm to {sigma}{sub x} = {sigma}{sub y} = 8 mm, to evaluate the impact of the proton beam size on the quality of intensity modulated proton therapy (IMPT) plans. The distance between spots, {Delta}x and {Delta}y, is optimized on the spot plane, ranging from 4 to 12 mm (i.e., each spot size is coupled with the best spot grid resolution). In our Hyperion treatment planning system (TPS), constrained optimization is applied with respect to the organs at risk (OARs), i.e., the optimization tries to satisfy the dose objectives in the planning target volume (PTV) as long as all planning objectives for the OARs are met. Three-field plans for a nasopharynx case, two-field plans for a prostate case, and two-field plans for a malignant pleural mesothelioma case are considered in our analysis. Results: For the head and neck tumor, the best grids (i.e., distance between spots) are 5, 4, 6, 6, and 8 mm for {sigma} = 3, 4, 5, 6, and 8 mm, respectively. {sigma} {<=} 5 mm is required for tumor volumes with low dose and {sigma}{<=} 4 mm for tumor volumes with high dose. For the prostate patient, the best grid is 4, 4, 5, 5, and 5 mm for {sigma} = 3, 4, 5, 6, and 8 mm, respectively. Beams with {sigma} > 3 mm did not satisfy our first clinical requirement that 95% of the prescribed dose is delivered to more than 95% of prostate and proximal seminal vesicles PTV. Our second clinical requirement, to cover the distal seminal vesicles PTV, is satisfied for beams as wide as {sigma} = 6 mm. For the mesothelioma case, the low dose PTV prescription is well respected for all values of {sigma}, while there is loss of high dose PTV coverage for {sigma} > 5 mm. The best grids have a spacing of 6, 7, 8, 9, and 12 mm for {sigma} = 3, 4, 5, 6, and 8 mm, respectively. Conclusions: The maximum acceptable proton pencil beam {sigma} depends on the volume treated, the protocol of delivery, and optimization of the plan. For the clinical cases, protocol and optimization used in this analysis, acceptable {sigma}s are {<=} 4 mm for the head and neck tumor, {<=} 3 mm for the prostate tumor and {<=} 6 mm for the malignant pleural mesothelioma. One can apply the same procedure used in this analysis when given a ''class'' of patients, a {sigma} and a clinical protocol to determine the optimal grid spacing.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

NASA NTRSDate not supplied

The Effect of Trailing Edge Grid Resolution in Large-Eddy Simulations

Ceramic Matrix Composite (CMC) with environmental barrier coating can tolerate significantly higher temperatures compared with the conventional metal alloy blades and thus can help achieve higher thermal efficiency in gas turbine engines. Due to the complex fabrication process of CMC blades, larger blade trailing edge thicknesses and larger leading-edge radii of curvature are expected. These features significantly alter aerodynamic performance especially near the trailing edge and thus the loss profile. In previous studies, the authors used the large eddy simulation (LES) as well as the Reynolds Averaged Navier-Stokes (RANS) coupled with an intermittency function-based transition model to simulate an aerodynamics performance of three different CMC blades and compared the predictions against the data acquired at NASA Glenn Transonic Turbine Blade Cascade Rig. Although we were able to successfully predict different characteristics of aerodynamic performance among three different blades (e.g., pressure loading and integrated losses), LES consistently predicted a small recirculation bubble at the end of the suction side. This was not inconsistent with the physics of the flow, but the pressure loading, locally, near the trailing edge, and as a result, the measured wake profiles did not match the experimental data. To delve further into these issues, a new LES grid with a much finer mesh near the trailing edge was generated where the spanwise extent of the grid was increased by 50%. An LES simulation for the CMC blade with the largest trailing edge thicknesses at the high-freestream turbulence (Tu=13%) condition showed that using such a mesh is necessary for LES to accurately capture the transition phenomena on the suction side. The resulting total pressure coefficient profile and the shift of the wake profile were greatly improved.

CMC↗

NASA NTRSDate not supplied

The Effect of Trailing Edge Grid Resolution in Large-Eddy Simulations

Ceramic Matrix Composite (CMC) with environmental barrier coating can tolerate significantly higher temperatures compared with the conventional metal alloy blades and thus can help achieve higher thermal efficiency in gas turbine engines. Due to the complex fabrication process of CMC blades, larger blade trailing edge thicknesses and larger leading-edge radii of curvature are expected. These features significantly alter aerodynamic performance especially near the trailing edge and thus the loss profile. In previous studies, the authors used the large eddy simulation (LES) as well as the Reynolds Averaged Navier-Stokes (RANS) coupled with an intermittency function-based transition model to simulate an aerodynamics performance of three different CMC blades and compared the predictions against the data acquired at NASA Glenn Transonic Turbine Blade Cascade Rig. Although we were able to successfully predict different characteristics of aerodynamic performance among three different blades (e.g., pressure loading and integrated losses), LES consistently predicted a small recirculation bubble at the end of the suction side. This was not inconsistent with the physics of the flow, but the pressure loading, locally, near the trailing edge, and as a result, the measured wake profiles did not match the experimental data. To delve further into these issues, a new LES grid with a much finer mesh near the trailing edge was generated where the spanwise extent of the grid was increased by 50%. An LES simulation for the CMC blade with the largest trailing edge thicknesses at the high-freestream turbulence (Tu=13%) condition showed that using such a mesh is necessary for LES to accurately capture the transition phenomena on the suction side. The resulting total pressure coefficient profile and the shift of the wake profile were greatly improved.

LES↗

DOE OSTI2017

A Hubble Space Telescope Survey for Novae in M87. II. Snuffing out the Maximum Magnitude–Rate of Decline Relation for Novae as a Non-standard Candle, and a Prediction of the Existence of Ultrafast Novae

The extensive grid of numerical simulations of nova eruptions from the work of Yaron et al. first predicted that some classical novae might significantly deviate from the Maximum Magnitude–Rate of Decline (MMRD) relation, which purports to characterize novae as standard candles. Kasliwal et al. have announced the observational detection of a new class of faint, fast classical novae in the Andromeda galaxy. These objects deviate strongly from the MMRD relationship, as predicted by Yaron et al. Recently, Shara et al. reported the first detections of faint, fast novae in M87. These previously overlooked objects are as common in the giant elliptical galaxy M87 as they are in the giant spiral M31; they comprise about 40% of all classical nova eruptions and greatly increase the observational scatter in the MMRD relation. We use the extensive grid of the nova simulations of Yaron et al. to identify the underlying causes of the existence of faint, fast novae. These are systems that have accreted, and can thus eject, only very low-mass envelopes, of the order of 10{sup −7}–10{sup −8} M {sub ⊙}, on massive white dwarfs. Such binaries include, but are not limited to, the recurrent novae. These same models predict the existence of ultrafast novae that display decline times, t {sub 2,} to be as short as five hours. We outline a strategy for their future detection.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Voronoi Meshing to Accurately Capture Geological Structure in Subsurface Simulations

Mesh generation lies at the interface of geological modeling and reservoir simulation. Highly skewed or very small grid cells may be necessary to accurately capture the geometry of geological features, but the resulting poorly scaled or small grid cells can have a substantial negative impact on simulator accuracy and speed. One way to minimize numerical errors caused by gridding complex structures is to simulate on high-quality Voronoi meshes, which reduce grid orientation effects in fluid flow. This work presents a complete methodology to create Voronoi simulation grids, model fluid flow in complex geological systems, and visualize the results. A recently developed Voronoi meshing method that can automatically generate provably good unstructured meshes that conform to input surfaces creating closed volumes is used. Initially an analytical benchmark simulation is presented to validate the quality of the meshes and simulation results and demonstrate the superiority of simulation results using Voronoi meshes over flexed-hexahedral meshes on a domain with internal features. Next, meshes are created for test structures representing four of the most common geological features in the subsurface: layering, pinch-out, an interior lens that tapers to zero thickness on all sides and a fault with offset. Two benchmark flow simulations are run for each test structure. Finally, a realistic geological example for CO 2 injection into an anticline is simulated. Three realizations of the Voronoi mesh at the same resolution are generated for the simulations. Each mesh is highly refined near the injection wells and coarse in areas of less interest. These three meshes are used to model the CO 2 plume in the subsurface as it migrates to the top of the structure and then fills downward. Simulations on the meshes with randomly generated elements inside the input volumes each give slightly different fingering patterns for the viscous-unstable buoyant gas flow. The results presented in this work show a promising step towards utilizing fully automated Voronoi meshing for subsurface flow simulations in complex geology.

58 GEOSCIENCES↗

DOE OSTI2023

Residual estimation for grid modification in wall-modeled large eddy simulation using unstructured high-order methods

Here, the accuracy and computational cost of a large eddy simulation are highly dependent on the computational grid. Building optimal grids manually from a priori knowledge is not feasible in most practical use cases; instead, solution-adaptive strategies can provide a robust and cost-efficient method to generate a grid with the desired accuracy. We adapt the residual estimation algorithm developed by Toosi and Larsson for Discontinuous Galerkin Spectral Elements Methods (DGSEM) to guide the grid-adaptation process. The core of the method is the computation of the estimated modeling residual using the polynomial basis functions used in DGSEM, and the averaging of the estimated residual over each element. The final method is assessed in multiple channel flow test cases and for the transonic flow over an airfoil, in both cases making use of mortar interfaces between elements with hanging nodes. The method is found to be robust and reliable, and to provide solutions on grids with significantly fewer elements at comparable accuracy compared to when using human-generated grids.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2019

Kestrel Results at Liftoff Conditions for a Space Launch System Configuration in Proximity to the Launch Tower

Aerodynamic data books for Space Launch System vehicles require databases for the integrated forces and moments and section loads during liftoff and transition to the ascent phase of flight. While the force and moment database can be generated from wind tunnel results, computational analyses are necessary to provide the extensive surface information required to generate proper lineloads. Of the two flight regimes, the liftoff problem is the more costly and complex situation to simulate, as it requires modeling of the vehicle in proximity to the launch tower. The effects of massive separation on the leeward pressure fields of both the tower and vehicle are not well captured with RANS methods, necessitating the use of more advanced methods, such as Delayed Detached Eddy Simulation, in conjunction with computational grids sufficiently refined to resolve the wakes. Details on the computational setup for employing the Kestrel flow solver to address the liftoff problem are presented. The methodology involves the use of independent unstructured near-body grids for the vehicle and the tower, overset by a solution adaptive Cartesian off-body grid. Results from the simulations are compared to experimental results from a test in the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel.

Computational fluid dynamics↗

NASA NTRS2001

Time-Dependent Simulation of Incompressible Flow in a Turbopump Using Overset Grid Approach

This paper reports the progress being made towards complete unsteady turbopump simulation capability by using overset grid systems. A computational model of a turbo-pump impeller is used as a test case for the performance evaluation of the MPI, hybrid MPI/Open-MP, and MLP versions of the INS3D code. Relative motion of the grid system for rotor-stator interaction was obtained by employing overset grid techniques. Unsteady computations for a turbo-pump, which contains 114 zones with 34.3 Million grid points, are performed on Origin 2000 systems at NASA Ames Research Center. The approach taken for these simulations, and the performance of the parallel versions of the code are presented.

Kiris, Cetin↗

DOE OSTI2018

On the role of finite grid extent in SOLPS-ITER edge plasma simulations for JET H-mode discharges with metallic wall

The impact of the finite grid size in SOLPS-ITER edge plasma simulations is assessed for JET H-mode discharges with a metal wall. For a semi-horizontal divertor configuration it is shown that the separatrix density is at least 30% higher when a narrow scrape-off layer (SOL) grid width is chosen in SOLPS-ITER compared to the case for which the SOL grid width is maximised. The density increase is caused by kinetic neutrals being not confined inside the divertor region because of the reduced extent of the plasma grid. In this case, an enhanced level of reflections of energetic neutrals at the low-field side (LFS) metal divertor wall is observed. This leads to a shift of the ionisation source further upstream which must be accounted for as a numerical artefact. An overestimate in the cooling at the divertor entrance is observed in this case, identified by a reduced heat flux decay parameters λ q div . Otherwise and further upstream the mid-plane heat decay length λ q parameter is not affected by any change in divertor dissipation. This confirms the assumptions made for the ITER divertor design studies, i.e. that λ q upstream is essentially set by the assumptions for the ratio radial to parallel heat conductivity. It is also shown that even for attached conditions the decay length relations λ ne > λ Te > λ q hold in the near-SOL upstream. Thus for interpretative edge plasma simulations one must take the (experimental) value of λ ne into account, rather than λ q , as the former actually defines the required minimum upstream SOL grid extent.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

A Machine‐Learning‐Assisted Stochastic Cloud Population Model as a Parameterization of Cumulus Convection

Abstract A machine‐learning‐assisted stochastic cloud population model is coupled with the Advanced Research Weather Research and Forecasting (WRF) model to represent fluctuations in the cloud‐base mass flux associated with the life cycles and interactions among cumulus convection cells. In this cloud population model, the size distribution and the associated cloud‐base mass flux of the convective cells are related to their previous state and to the change in the total convective area via a transition function. The convective area tendency in turn is assumed to depend on the cloud‐base mass flux that is resolved by the host WRF model. The transition function is represented by a single hidden‐layer neural network trained by the evolution of convective cell size distributions in a 1‐km grid‐spacing WRF simulation run over the Australian Monsoon region. At every grid point of the host model, the cloud population model predicts the cell size and cloud‐base mass flux distributions from which a random sample of cells is fed to an entraining parcel model that calculates precipitation as well as the associated liquid water potential temperature and total moisture tendencies. These tendencies are averaged over the cells and provided to the host model. Several regional simulations are performed over tropical and midlatitude domains to test this as a potential approach to scale‐aware parameterization. It is shown that such an approach could be a new promising path to simulating realistic precipitation statistics and propagation of precipitation associated with the Madden‐Julian Oscillation while maintaining realistic depictions of the diurnal cycle over both land and ocean.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Examining System-Wide Impacts of Solar PV Control Systems with a Power Hardware-in-the-Loop Platform

High penetration levels of distributed solar PV power generation can lead to adverse power quality impacts such as excessive voltage rise, voltage flicker, and reactive power values that result in unacceptable voltage levels. Advanced inverter control schemes have been proposed that have the potential to mitigate many power quality concerns. However, closed-loop control may lead to unintended behavior in deployed systems as complex interactions can occur between numerous operating devices. In order to enable the study of the performance of advanced control schemes in a detailed distribution system environment, a Hardware-in-the-Loop (HIL) platform has been developed. In the HIL system, GridLAB-D, a distribution system simulation tool, runs in real-time mode at the Pacific Northwest National Laboratory (PNNL) and supplies power system parameters at a point of common coupling to hardware located at the National Renewable Energy Laboratory (NREL). Hardware inverters interact with grid and PV simulators emulating an operational distribution system and power output from the inverters is measured and sent to PNNL to update the real-time distribution system simulation. The platform is described and initial test cases are presented. The platform is used to study the system-wide impacts and the interactions of controls applied to inverters that are integrated into a simulation of the IEEE 8500-node test feeder, with inverters in either constant power factor control or active volt/VAR control. We demonstrate that this HIL platform is well-suited to the study of advanced inverter controls and their impacts on the power quality of a distribution feeder. Additionally, the results from HIL are used to validate GridLAB-D simulations of advanced inverter controls.

14 SOLAR ENERGY↗

DOE OSTI2023

Recent Hydrodynamic Modeling Enhancements in OpenFAST

To improve the hydrodynamic modeling capabilities of the offshore wind design and modeling tool OpenFAST, a restructuring of the HydroDyn hydrodynamics module was undertaken with several new features implemented. The generation of the wave field is now separated from HydroDyn into a new module called SeaState. Unlike previous versions of HydroDyn, which precomputed the wave kinematics at the undisplaced positions of the hydrodynamic nodes, the new SeaState module computes the complete time history of the wave field at the vertices of a user-defined wave grid. During the simulation, the wave kinematics at any point within the grid can be efficiently interpolated, allowing the strip-theory wave loads to be evaluated based on the wave kinematics at the displaced structure position. A phase correction was also implemented for the potential-flow wave-exciting loads to account for large structure displacements. Several wave-stretching methods were implemented for the strip-theory solution along with a wave-load redistribution method that ensures smooth variation of the nodal loads as the nodes of a discrete hydrodynamic mesh enter and exit the water. The load smoothing is included to avoid the excitation of unphysical high-frequency structural vibration during hydroelastic simulations. The MacCamy-Fuchs diffraction correction for the strip-theory solution and the capability to insert a constrained NewWave into extreme stochastic sea states with directional spreading were also added. The present article documents the formulations of the new improvements to HydroDyn with example applications and numerical results.

constrained wave↗

DOE OSTI2022

Recent Hydrodynamic Modeling Enhancements in OpenFAST: Preprint

To improve the hydrodynamic modeling capabilities of the offshore wind design and modeling tool OpenFAST, a restructuring of the HydroDyn hydrodynamics module was undertaken with several new features implemented. The generation of the wave field is now separated from HydroDyn into a new module called SeaState. Unlike previous versions of HydroDyn, which precomputed the wave kinematics at the undisplaced positions of the hydrodynamic nodes, the new SeaState module computes the complete time history of the wave field at the vertices of a user-defined wave grid. During the simulation, the wave kinematics at any point within the grid can be efficiently interpolated, allowing the strip-theory wave loads to be evaluated based on the wave kinematics at the displaced structure position. A phase correction was also implemented for the potential-flow wave-exciting loads to account for large structure displacements. Several wave-stretching methods were implemented for the strip-theory solution along with a wave-load redistribution method that ensures smooth variation of the nodal loads as the nodes of a discrete hydrodynamic mesh enter and exit the water. The load smoothing is included to avoid the excitation of unphysical high-frequency structural vibration during hydroelastic simulations. The MacCamy-Fuchs diffraction correction for the strip-theory solution and the capability to insert a constrained NewWave into extreme stochastic sea states with directional spreading were also added. The present article documents the formulations of the new improvements to HydroDyn with example applications and numerical results.

constrained wave↗

DOE OSTI2025

Control Design and EMT Simulation Toward the Restoration of Faulted Wind-Dominant Grids

Present power engineers have been challenged to employ electromagnetic transient (EMT) simulations to study the restoration of wind-dominant grids. This paper addresses this hurdle by developing strategies to: (i) reliably control the dc-link voltage of grid-forming wind turbines with dc-coupled batteries, (ii) tune controllers to stably interconnect wind power plants into a grid under recovery, and (iii) autonomously make restoration decisions if transmission faults occur during restoration. These contributions are assessed via EMT simulations of wind-dominant versions of the WSCC 9-bus grid and a 22-bus power system. Furthermore, this paper is significant to address recommendations by the North American Electric Reliability Corporation.

17 WIND ENERGY↗

DOE OSTI2012

Simulation and Analysis of Vehicle-to-Grid Operations in Microgrid

Abstract not provided

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2016

Temperature and precipitation extremes in century-long gridded observations, reanalyses, and atmospheric model simulations

Knowledge about long-term changes in climate extremes is vital to better understand multidecadal climate variability and long-term changes and to place today's extreme events in a historical context. While global changes in temperature and precipitation extremes since the midtwentieth century are well studied, knowledge about century-scale changes is limited. This paper analyses a range of largely independent observations-based data sets covering 1901–2010 for long-term changes and interannual variability in daily scale temperature and precipitation extremes. We compare across data sets for consistency to ascertain our confidence in century-scale changes in extremes. We find consistent warming trends in temperature extremes globally and in most land areas over the past century. For precipitation extremes we find global tendencies toward more intense rainfall throughout much of the twentieth century; however, local changes are spatially more variable. While global time series of the different data sets agree well after about 1950, they often show different changes during the first half of the twentieth century. In regions with good observational coverage, gridded observations and reanalyses agree well throughout the entire past century. Finally, simulations with an atmospheric model suggest that ocean temperatures and sea ice may explain up to about 50% of interannual variability in the global average of temperature extremes, and about 15% in the global average of moderate precipitation extremes, but local correlations are mostly significant only in low latitudes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016

Large-eddy simulation, fuel rod vibration and grid-to-rod fretting in pressurized water reactors

Grid-to-rod fretting (GTRF) in pressurized water reactors is a flow-induced vibration phenomenon that results in wear and fretting of the cladding material on fuel rods. GTRF is responsible for over 70% of the fuel failures in pressurized water reactors in the United States. Predicting the GTRF wear and concomitant interval between failures is important because of the large costs associated with reactor shutdown and replacement of fuel rod assemblies. The GTRF-induced wear process involves turbulent flow, mechanical vibration, tribology, and time-varying irradiated material properties in complex fuel assembly geometries. This paper presents a new approach for predicting GTRF induced fuel rod wear that uses high-resolution implicit large-eddy simulation to drive nonlinear transient dynamics computations. The GTRF fluid–structure problem is separated into the simulation of the turbulent flow field in the complex-geometry fuel-rod bundles using implicit large-eddy simulation, the calculation of statistics of the resulting fluctuating structural forces, and the nonlinear transient dynamics analysis of the fuel rod. Ultimately, the methods developed here, can be used, in conjunction with operational management, to improve reactor core designs in which fuel rod failures are minimized or potentially eliminated. Furthermore, robustness of the behavior of both the structural forces computed from the turbulent flow simulations and the results from the transient dynamics analyses highlight the progress made towards achieving a predictive simulation capability for the GTRF problem.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗