Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GRID MODELING”

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 919 records · Page 51

Numerical Simulation of Two-grid Ion Optics Using a 3D Code

A three-dimensional ion optics code has been developed under NASA's Project Prometheus to model two grid ion optics systems. The code computes the flow of positive ions from the discharge chamber through the ion optics and into the beam downstream of the thruster. The rate at which beam ions interact with background neutral gas to form charge exchange ions is also computed. Charge exchange ion trajectories are computed to determine where they strike the ion optics grid surfaces and to determine the extent of sputter erosion they cause. The code has been used to compute predictions of the erosion pattern and wear rate on the NSTAR ion optics system; the code predicts the shape of the eroded pattern but overestimates the initial wear rate by about 50%. An example of use of the code to estimate the NEXIS thruster accelerator grid life is also presented.

ion optics↗

On the Synergy Between Numerics and Subgrid Scale Modeling in LES of Stratified Flows: Grid Convergence of a Stratocumulus-Topped Boundary Layer

The effectiveness of a linear upwinding scalar advection scheme to suppress numerical dispersion errors near sharp inversions in large-eddy simulations of a nocturnal stratocumulus-topped boundary layer is assessed. Linear upwinding is a trade-off between non-dissipative and non-linear positive definite advection schemes. It is shown that linear upwinding does not negatively impact the model's grid convergence properties and a sharp inversion free of numerical artifacts is maintained. Even though mean profiles and turbulence fluxes show good grid convergence characteristics the liquid water amount varies significantly with grid resolution. The entrainment rate is identical for all resolutions and independent of the liquid water amount. For the present stratocumulus case, the impact of cloud-top radiative cooling is negligible and turbulence is largely driven by convection emanating from the surface.

Matheou, Georgios↗

A Grid Intensity Map-Based Source Model for Monte Carlo Modeling of IMRT Field Irradiation - Paper 5

An accurate source model is important for TPS (Treatment Planning System) dose calculation and verification, especially for IMRT (Intensity Modulated Radiotherapy) irregular and inhomogeneous fields. Monte Carlo simulation is the most accurate way to calculate dose distributions in the heterogeneous media like human bodies. A full radiotherapy modeling generally includes Linac modeling and phantom dose calculation. Full Linac modeling is very time-consuming and requires highly detailed knowledge of the accelerator composition and structure, however, which is very difficult to obtain for commercial interests. Often all relevant parameters are not known with sufficient accuracy, and thus a time-consuming tuning of e.g. the electron beam spot size and energy must be performed. A multiple source model (MSM) regards particles coming from different parts of an accelerator as if they are coming from different virtual sources. MSM uses a set of virtual sources representing the treatment beam (e.g. primary photons, head scatter, electron contamination) by a set of analytic functions or histograms with parameters derived from dosimetric measurements or a full MC simulation. Tillikainen et al built a three sub-source model by matching with measurement data using an automatic process to optimize the free model parameters for a Varian 2100 C/D accelerator. Lin et al have built a convenient five-photon-sub-source model for Varian 2300C 6 MV-X beam based on Monte Carlo simulated Linac for ARTS (Accurate Radiotherapy System). The particle generation of MSM is very fast compared to full Linac simulation, and its sampling efficiency is not limited by the size of the phase space file (PSF). However, generally MSM is generated based on regular field PSF and needs an irregular component to form the beam shape meanwhile apt to bringing error Linac modeling and MSM producing from PSF analysis actually both are multi-freedom optimization process by adjusting several parameters (e.g. electron energy and angular divergence). Linac modeling adjusts its parameters from the beginning of electron hitting target through all components to fit the simulation to the measurement. MSM tunes its PSF-derived parameters by optimizing to minimize the deviations between the simulation and the measurement. Since the full Linac modeling and the MSM model both need to fine tune their parameters to match the simulation to the measurement, and sometimes maybe offsets the real accelerator for its multi-parameter optimization process. So, why not fine tune directly all source model parameters just behind all Linac components to fit the simulation to the measurement? A grid intensity-based surface source model was suggested here, which assumed the particle emitting from the bottom surface of Linac last beam trimmer MLC (Multi Leaf Collimator). The particle emitting position is sampled on the grid intensity map optimized by TPS, and the particle weight is adjusted by its original grid intensity. This source model not only takes advantage of the Monte Carlo method in the heterogeneous media, but also avoids effectively the transportation speed shortcoming of Monte Carlo simulation through Linac components. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Hybrid-RL-MPC4CLR (Hybird-Reinforcement-Learning-Model-Predictive-Control-for-Reserve-Policy-Assisted-Critical-Load-Restoration-in-Distribution-Grids)

Hybrid-RL-MPC4CLR was developed as a hybrid controller for active distribution grid critical load restoration, combining deep reinforcement learning (RL) and model predictive control (MPC) aiming at maximizing total restored load following an extreme event. The RL determines a policy for quantifying operating reserve requirements, thereby hedging against uncertainty, while the MPC models grid operations incorporating the RL policy actions (i.e., reserve requirements), renewable (wind and solar) power predictions, and load demand forecasts. The developers formulated the reserve requirement determination problem as a sequential decision-making problem based on the Markov Decision Process (MDP) and design an RL learning environment based on the OpenAI Gym framework and MPC simulation. The RL agent reward and MPC objective function aim to maximize and monotonically increase total restored load and minimize load shedding and renewable power curtailment. The software is developed using various software packages in Python. The MPC's optimal power flow (OPF) model is implemented using the Pyomo package, the RL simulation environment is implemented using the MPC simulation with various scenarios of renewable energy and load demand profiles and power outage beginning times, based on the OpenAI Gym framework. The RL agent training is performed using the RLlib Ray package. The RL algorithm is trained offline using historical forecasts of renewable generation and load demand profiles. Simulation analysis and performance tests are conducted using a modified IEEE 13-bus distribution test feeder containing wind turbine, photovoltaic, microturbine, and battery.

Eseye, Abinet Tesfaye↗

Robust trajectory-constrained frequency control for microgrids considering model linearization error

Grid supportive modes integrated within inverter-based resources can improve the frequency response of renewable-rich microgrids. The synthesis of grid supportive modes to guarantee frequency trajectory constraints under a predefined disturbance set is challenging but essential. To tackle this challenge, a numerical optimal control (NOC)-based control synthesis methodology is proposed. Without loss of generality, a wind-diesel fed microgrid is studied, where we aim to design grid supportive functions in the wind turbine. In the control design, linearized models are used, and the linearization-induced errors are quantitatively analyzed by reachability and interval arithmetics and represented in the form of interval uncertainties. Then, the NOC problem can be formulated into a robust mixed-integer linear program. The control structure is strategically configured into two levels to realize online deployment. The proposed control is verified on the modified 33-node microgrid with a full-order three-phase nonlinear model in Simulink. In conclusion, the simulation results show the effectiveness of the proposed control paradigm and the necessity of considering linearization-induced uncertainty.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Influence of SMAP Soil Moisture Retrieval Assimilation on Runoff Estimation Across South Asia

This study was designed to characterize and quantify the influence of surface soil moisture assimilation on estimated runoff (surface flow and baseflow) and hydraulically-routed streamflow across three large river basins in South Asia that are at risk of impending water stress. Soil Moisture Active Passive (SMAP) surface soil moisture retrievals were assimilated into the Noah-MP land surface model. The gridded runoff was hydraulically routed to obtain volumetric streamflow values for a river network using a runoff routing module (Hydrological Modeling and Analysis Platform, HyMAP). The open loop (OL, model-only) and data assimilated (DA, includes soil moisture retrieval assimilation) Noah-MP runoff estimates highlighted the improvements in estimated total runoff across irrigated areas. Soil moisture assimilation impacted baseflow more relative to surface runoff. The HyMAP-based OL and DA streamflow generally underestimated the streamflow at upstream stations and overestimated the streamflow at downstream stations within the Indus basin due to missing physics related to reservoir operations. For stations located in the Ganges–Brahmaputra basins, the OL and DA estimation performance varied. The OL and DA results showed that the assimilation of soil moisture retrievals improves gridded runoff and volumetric streamflow across irrigated areas. Considerable relative change in streamflow (>70% increase in magnitude relative to OL) is noted after assimilation across the highly irrigated lower Indus basin and high precipitation regions in Bangladesh. Improving the modeling system’s representativeness of ground conditions via inclusion of water management information could improve large-scale streamflow modeling across South Asia.

soil moisture↗

Atmospheric model development in support of SEASAT. Volume 4: Forecast model sensitivity study

Atmospheric analysis and prediction models of varying (grid) resolution were developed. The models were tested using real observational data for the purpose of assessing the impact of grid resolution on short range numerical weather prediction. The work statement was amended to include the performance of sensitivity tests using a coarse mesh (63 x 63 x 5 level) prediction model in order to identify and order factors which might mask or impair the utility of SEASAT data on short range weather prediction. Such factors included: initial conditions; topography; surface friction; latent heating; diffusion of momentum and temperature; and computational devices such as tendency truncators, pressure smoothers, and temporal filters.

Lewit, H. L.↗

Fidelity Analysis of Complex Large-Scale Simulation Models of PEs (MARS) in Future Power Grids

With the proliferation of power electronic (PE)-based resources in the electric grid, there have been numerous different kinds of simulation models developed and studied to understand their dynamic behavior. Study of PE-based resources using each simulation model has its own distinct advantages and disadvantages. Therefore it is very important to carefully select the appropriate model for the evaluation of the PE-based resources. In this paper, one such effort is made to analyze the different types of models in the case of a multi-port autonomous reconfigurable solar power plant (MARS). The MARS integrates a hybrid photovoltaic (PV)-energy storage system (ESS) plant to high-voltage direct current (HVdc) system and transmission alternating current (ac) grid. The dynamic models are evaluated under different test cases varying the short-circuit ratio (SCR) and control algorithms. Recommendations are provided on the equipment study needs based on the analysis completed.

Marthi, Phani Ratna Vanamali↗

Hybrid Modeling of Three-Phase Grid-Supporting Inverters for Dynamic Studies

Grid technologies connected by power electronic converter (PEC) interfaces continually implement grid support functions mandated by grid codes and standards. The transition to converter-based generation demands precise PEC models to assess system dynamics, which have been previously overlooked in conventional power systems. This study proposes a hybrid method for analyzing grid-connected three-phase PEC dynamics with the IEEE standard 1547-2018 Volt-VAr mode that combines physics and data-driven techniques. The physics model reflects the PEC’s internal behavior, whereas the data-driven modeling technique evaluates the grid-supporting capabilities of the smart PEC. The system identification approach is used to generate dynamic PEC models based on changing grid voltage and measured current injected into the grid by the PEC. In the Volt-VAr support mode, a detailed topological model including switches is utilized to compare the goodness-of-fit of the extracted hybrid dynamic model. The results demonstrate that the hybrid PEC model in the Volt-VAr mode accurately matches the dynamics with the topological model.

Subedi, Sunil↗

Theory for transport in magnetized plasmas

We present a two fluid model for anomalous thermal, particle and momentum transport which has been derived from kinetic theory using a small parameter ε~10-2. This small parameter can come from either the ratio of wave frequency to cyclotron frequency or the ratio of potential to thermal energies. This fluid model is strongly nonlinear and has an exact closure within the accuracy ε. It can be used to simulate the full radius of a tokamak including internal and edge transport barriers at the same time without introducing any local assumptions for physics model or grid size. It involves a full transport matrix, thus giving possibility for both density and heat pinches. The model can be derived from the Vlasov equation, adding also close collisions and includes particle pinches of the type found in the magnetosphere and in the Massachusetts Institute of Technology elevated dipole machine. Pinches are also in agreement for the heat pinch experiments on Doublet III D including simultaneous particle and heat pinches. For tokamaks it gives the right scaling of confinement time versus heating power. Finally, we have also recovered the Turbulent Equipartition result that the density scales as inverse safety factor q in typical simulations. Thus this model is both fundamental and of high practical use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transonic analysis of the F-16A with under-wing fuel tanks - An application of the TranAir full-potential code

Initial results obtained using the TranAir transonic full-potential code are presented for a fully configured F-16A geometry, including two 370-gal under-wing fuel tanks, at freestream Mach numbers of 0.6 and 0.9 and angle of attack of 4 deg. The geometry is modeled using surface panels, and the flow field is defined by a rectangular array of flow-field grid points. The paneled geometry is embedded within this rectangular flow field. By avoiding the use of surface-conforming flow-field grids, the modeling generality afforded by surface panels can be utilized to analyze very complex configurations in the transonic flow regime.

Madson, Michael D.↗

Automated rectification and geocoding of SAR imagery

An automated post-processing system has been developed for rectification and geocoding of SAR (Synthetic Aperture Radar) imagery. The system uses as input a raw uncorrected image from the operational SAR correlator, and produces as a standard output a rectified and geocoded product. The accurate geolocation of SAR image pixels is provided by a spatial transformation model which maps the slant range-azimuth SAR image pixels into their location on a prespecified map grid. This model predicts the geodetic location of each pixel by utilizing: the sensor platform position; a geoid model; the parameters of the data collection system and the processing parameters used in the SAR correlator. Based on their geodetic locations, the pixels are mapped by using the desired cartographic projection equations. This rectification and geocoding technique has been tested with Seasat and SIR-B images. The test results demonstrate absolute location uncertainty of less than 50 m and relative distortion (scale factor and skew) of less than 0.1 percent relative to local variations from the assumed geoid.

Kwok, R.↗

Bgridi - Interactive three-dimensional turbomachinery grid generation system with applications

An interactive, graphics-based grid generation system for 3D turbomachinery geometries has been developed. The system consists of separate modules for geometry modeling and grid generation. The grid generation method generates a series of 2D grids in the blade-to-blade passage to build up the 3D grid. A Poisson equation with forcing functions selected to control orthogonality and spacing on all boundaries is solved to generate the grid. Interactive definition of Bezier curves and surfaces as internal boundaries is used to improve control of grid quality. A multi-block data structure simpifies the creation of structured H-grids about complex turbomachinery geometries. The multi-block strategy facilitates the creation of a grid in the tip flow region, which is an important contributor to losses within the passage. The grid generation system is applied to several complex geometries including a simple radial impeller, a conventional turbine rotor, the SSME LOX turbine rotor with and without a tip flow grid, a multi-passage inducer-impeller, and a tip-flow cascade. In addition, a calculation of the tip-flow cascade flow field is shown.

Shoemaker, J. M.↗

Positive-Sequence Modeling of Droop-Controlled Grid-Forming Inverters for Transient Stability Simulation of Transmission Systems

Here, this paper describes a positive-sequence model to represent two widely reported droop-controlled grid-forming inverters for bulk power systems simulation study. Methods of how to develop the equivalent voltage source behind impedance to represent inverters with and without inner control loops, modeling of P-f and Q-V droop controls, active and reactive power limiting modeling, and modeling of fault current limiting controls are described in detail. The model is implemented in commercially available positive-sequence simulation tools and compared with detailed electromagnetic transient models in both a single-grid-forming-inverter infinite-bus system and a modified IEEE 39-bus system. Finally, the model is tested on the full US Western Electricity Coordinating Council (WECC) system. Study results show that the model has a good level of accuracy compared to the detailed EMT models, and at the same time also achieves a high computational efficiency, which is suitable for the bulk power system simulation study.

42 ENGINEERING↗

Expanding the computational frontier of multi-scale atmospheric simulation to advance understanding of low cloud / climate feedbacks: Final Technical Report

Marine boundary layer clouds are a challenge for climate models. They cover much of the oceans, but they are driven by small-scale turbulent eddies only a few hundred meters across, which in turn respond to the cloud formation. Global climate models have a grid spacing that is far too large to simulate such fluid motions, so they parameterize these cloud formation processes. Different climate modeling groups have designed different parameterizations in which the clouds turn out to be differently sensitive to a warming climate. This drives uncertainties in our best guess at the sensitivity of global warming to greenhouse gas increases. If these cloud-forming eddies could be directly simulated using the well-known equations of fluid motion, they would no longer need to be parameterized, removing a major source of climate modeling uncertainty. In this project, we overcame software engineering challenges to successfully implemented ‘ultraparameterization’ (UP), the first global model that does this, and we tested how well it works. We simulated five-year periods with present-day temperatures and with a warmer climate, and we investigated how the UP-simulated clouds responded to climate – the ‘cloud feedback’ problem. We found little response of clouds at all latitudes to the imposed climate change, which is within the range of predictions of conventional global climate models. UP is a variation on superparameterization, in which small cloud-resolving models (CRMs) are embedded in each column of the global model. In UP, the CRM grid is fine enough (250 m horizontal × 20 m vertical) to explicitly capture boundary-layer turbulent eddies and associated clouds. Because only one small columnar patch is simulated within each climate model grid cell, this is a million-fold more efficient than simulating the entire globe on this same CRM grid., but achieves much of the same effect for the cloud properties. It doesn’t work perfectly. For instance, like conventional climate models, UP simulates too little subtropical stratocumulus cloud, a bias that we are continuing to work to reduce. However, because it directly simulates the turbulent cloud-forming processes, UP is inherently more plausible for simulating how clouds will change in a perturbed climate. In future, we hope to apply UP to another key climate modeling issue: cloud-aerosol interaction and the effect of human-produced aerosols on the climate change we have already experienced and that which is likely to come.

54 ENVIRONMENTAL SCIENCES↗

Computational Analysis of Cryogenic Flow Through a Control Valve

The initial efforts to develop the capability to model valves used in rocket engine component testing at Stennis Space Center are documented. An axisymmetric model of a control valve with LN2 as the working fluid was developed. The goal was to predict the effect of change in the plug/sear region of the valve prior to testing. The valve flow coefficient was predicted for a range of plug positions. Verification of the calculations was carried out to quantify the uncertainty in the numerical answer. The modeled results compared well qualitatively to experimental trends. Additionally, insights into the flow processes in the valve were obtained. Benefits from the verification process included the ability to use coarser grids and insight into ways to reduce computational time by using double precision accuracy and non-integer grid ratios. Future valve modeling activities will include shape optimization of the valve/seat region and dynamic grid modeling.

Danes, Russell↗

A Model Calibration Method for Grid-Forming Inverters Using Iterative Bayesian Optimization

As inverter-based resources (IBRs) are rapidly deployed, especially at the distribution and microgrid levels, the need to include their accurate representations in power systems models increases. With well-calibrated models, IBR-interactions can be examined, preventing any stability and operational issues, especially in islanded systems. This paper proposes a method to calibrate generic inverter models using Bayesian optimization, but with a parameter-grouping approach to improve speed and accuracy. This approach is illustrated for a grid-forming inverter using synthetic data and field measurements from a simple IBRbased microgrid. Tests show that the calibration algorithm has modest computation cost and is robust to noisy measurements..

Biswas, Shuchismita↗

Continuous Dynamic Grid Adaptation In A Global Atmospheric Model: Application And Refinement (Scientific Final Report)

This project had goals of advancing the performance capabilities of the numerical general circulation model EULAG and using it to produce a fully operational atmospheric global climate model (AGCM) that can employ either static or dynamic grid stretching for targeted phenomena. The resulting AGCM combined EULAG's advanced dynamics core with the 'physics' of the NCAR Community Atmospheric Model (CAM). Effort discussed below shows how we improved model performance and tested both EULAG and the coupled CAM-EULAG in several ways to demonstrate the grid stretching and ability to simulate very well a wide range of scales, that is, multi-scale capability. We leveraged our effort through interaction with an international EULAG community that has collectively developed new features and applications of EULAG, which we exploited for our own work summarized here. Overall, the work contributed to over 40 peer-reviewed publications and over 70 conference/workshop/seminar presentations, many of them invited.

54 ENVIRONMENTAL SCIENCES↗