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,405 records · Page 78

Predictive Model for Particle Residence Time Distributions in Riser Reactors. Part 1: Model Development and Validation

Here in this computational study, we model the mixing of biomass pyrolysis vapor with solid catalyst in circulating riser reactors with a focus on the determination of solid catalyst residence time distributions (RTDs). A comprehensive set of 2D and 3D simulations were conducted for a pilot-scale riser using the Eulerian-Eulerian two-fluid modeling framework with and without sub-grid-scale models for the gas-solids interaction. A validation test case was also simulated and compared to experiments, showing agreement in the pressure gradient and RTD mean and spread. For simulation cases, it was found that for accurate RTD prediction, the Johnson and Jackson partial slip solids boundary condition was required for all models and a sub-grid model is useful so that ultra high resolutions grids that are very computationally intensive are not required. Finally, we discovered a 2/3 scaling relation for the RTD mean and spread when comparing resolved 2D simulations to validated unresolved 3D sub-grid-scale model simulations.

09 BIOMASS FUELS↗

Nonlinear simulations of thermo-resistive tearing mode formalism of the density limit

This work is an extension of the previous semi-analytical work (Teng et al 2016 Nucl. Fusion 56 106001), explaining the density limit with a thermo-resistive tearing mode model. The thermo-resistive island growth is calculated with a 3D MHD code M3D-C 1 (Jardin et al 2012 Comput. Sci. Discovery 5 014002). It is shown with nonlinear 3D MHD simulations that impurity radiation stimulates large magnetic islands at the density limit. The impact of thermal perturbations inside and outside of the island is explored. Inside the island, net cooling enhances the island growth and net heating suppresses its growth. Outside the island, thermal perturbations have a much smaller impact on the island growth. When the plasma density is increased towards the density limit, the island changes from being heated to being cooled and grows to a much larger island width. The convergence test of the simulations over the temporal and spatial grid is performed. In conclusion, the numerical model and the semi-analytical model are compared by calculating the Δ' terms, which are defined in the modified Rutherford equation, and good agreement is observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Model Based Predictive Control for Frequency Support in Multi-port Autonomous Reconfigurable Solar Plants

Primary frequency support functionality can be an added advantage for any new grid integrated power electronics resource. With increased penetration of power electronics resources, inertial and primary frequency response capabilities are decreasing. In this research paper, an advanced model-based predictive control (MBPC) algorithm for providing frequency support to a new topology of integrated photovoltaic(PV), batterybased energy storage systems (ESS), and high-voltage direct current (HV dc) systems called multi-port autonomus reconfigurable solar (MARS) plants is based on the synchronverter-based control. The simulation of MARS-HV dc system with proposed control strategy is simulated and validated for MARS connected to low short circuit ratio (SCR) grids in PSCAD/EMTDC simulation environment. The control algorithm proposed showcases better performance in terms of frequency nadir improvement and frequency steady state improvemnt with respect to no frequency control.

Marthi, Phani Ratna Vanamali↗

Data for Grogan et al. "Bringing Hydrologic Realism to Water Markets"

This data set provides model output and post-processing files required to reproduce the results, tables, and figures in the paper "Bringing Hydrologic Realism to Water Markets" by Grogan et al. (in review). Other input data used in this study includes: Lisk, M., Grogan, D., Zuidema, S., Caccese, R., Peklak, D., Zheng, J., Fisher-Vanden, K., Lammers, R., Olmstead, S., & Fowler, L. (2023). Harmonized Database of Western U.S. Water Rights (HarDWR) (Version v1) [Data set]. MSD-LIVE Data Repository. https://doi.org/10.57931/2205619 Two models were used in this study: (1) The University of New Hampshire Water Balance Model WBM, and (2) a Water Market Model. Market model code and model output post-processing code that make use of these data can be found here Model output files are: 1. WBM output files: scenario[x]_wbm_output.zip Where [x] is one of 1, 2, 2a, 3, and 3a Each zipped directory contains 7 gridded NetCDF files, each reporting the 10-year annual average value of a given variable, in units of average mm/day: File Name: wbm_indUseGross_yc.nc; Description: Water withdrawals by industry (part of the urban sector) File Name: wbm_domUseGross_yc.nc; Description: Water withdrawals by the domestic sector (part of the urban sector) File Name: wbm_irrigationGross_yc.nc; Description: Water withdrawals for agriculture File Name: wbm_irrigationExtra_yc.nc; Description: Water withdrawals from unsustainable groundwater for agriculture File Name: wbm_indUseEvap_yc.nc; Description: Consumptive water use by industry File Name: wbm_domUseEvap_yc.nc; Description: Consumptive water use by the domestic sector File Name: wbm_irrigationNet_yc.nc; Description: Consumptive water use by agriculture The file full_cell_area.nc gives the area of each grid cell in km2, which is used for converting water depth to water volume. 2. Water market model output & post processing output Folder: marketTrdSummaries/ Description: Files in this folder are used as input to code 1_WelfareCalculation_actual_trades.R. They summarize historical water right trade transactions in each state. File Name: welfare_gain_by_state_sector.csv; Description: Welfare gains by state and sector, as shown in Figure 3F. Used in code Figure3.R and produced (as a .xlsx file) by code 2_DemandCurves_simulated_trades.R File Name: welfare_data_actual.rdata; Description: welfare gains by WMA from actual historical trades, as shown in Figure 3A. This data is the output of code 1_WelfareCalculation_actual_trades.R File Name: welfare_summary_simulated.xlsx; Description: Welfare gains by state as simulated by the market model in Scenario 1. Produced by code 2_DemandCurves_simulated_trades.R, and used in code 4_WelfareCalculation.R. File Name: welfare_summary_cutoffs.xlsx; Description: Welfare gains by state as simulated by the market model in Scenario 2. Produced by code 3_DemandCurves_simulated_trades_cutoffs.R, and used in code 4_WelfareCalculation.R. File Name: welfare_summary_cutoffs_SGMS.xlsx; Description: Welfare gains by state as simulated by the market model in Scenario 2a. Produced by code 3_DemandCurves_simulated_trades_cutoffs.R, and used in code 4_WelfareCalculation.R. File Name: welfare_data_actual.rdata; Description: Spatial data, actual historical welfare gains by WMA as shown in Figure 3A. Produced by code 4_WelfareCalculation.R and used by code Figure3.R. File Name: welfare_data_simulated.rdata; Description: Spatial data, simulated Scenario 1 welfare gains by WMA as shown in Figure 3B. Produced by code 4_WelfareCalculation.R and used by code Figure3.R. File Name: welfare_data_simulated_cutoffs.rdata; Description: Spatial data, simulated Scenario 2 welfare gains by WMA. Produced by code 4_WelfareCalculation.R and used by code Figure3.R. File Name: welfare_data_simulated_cutoffs_SGMA.rdata; Description: Spatial data, simulated Scenario 2a welfare gains by WMA. Produced by code 4_WelfareCalculation.R and used by code Figure3.R. Additional files are provided for efficient reproduction of tables and figures. These include: File Name: wma_thresold_dates_Scenario2(a).csv; Description: Wet vs. paper right threshold dates for each WMA. Shown in Figure 2A,B. Produced and used by code calculate_thresolds_Figure2.R File Name: WWRTradeBounds (directory); Description: Trade boundary shapefile required to reproduce Figure 3A-D. Used in code Figure3.R File Name: welfare_region_totals.csv; Description: Welfare gains for the entire study region, as shown in Figure 3E. Used in code Figure3.R File Name: Welfare_gain_by_state_sector.csv; Description: Welfare gains by state and sector, as shown in Figure 3F. Used in code Figure3.R and produced (as a .xlsx file) by code 2_DemandCurves_simulated_trades.R File Name: WECC_MERIT_5min_v3b_mask.nc; Description: Gridded file that identified which land grid cells are in the WBM model domain, used for processing in code Figure4.py File Name: Table_1.csv; Description: All data in Table 1, reproducible from WBM output files using code table_1.R

Economics↗

OVERFLOW Analysis of Supersonic Retropropulsion Testing on a Blunt Mars Entry Vehicle Concept

Supersonic Retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle was simulated using the OVERFLOW Computational Fluid Dynamics (CFD) solver. Two nozzle configurations were tested (1E and 1F) as a subset of the seven total configurations in the Descent Systems Study (DSS) testing campaign. A generalized, Adaptive Mesh Refinement (AMR) shock capturing and plume refinement technique was developed and calibrated for producing automatic and optimized grid systems unique to any given testing condition. Simulations were subjected to successive levels of grid refinement in order to determine solution independence from grid resolution. OVERFLOW Reynolds Averaged Navier-Stokes (RANS) solutions produced realistic SRP flow phenomena, including Mach disk normal shocks contained in high-thrust plumes and bow shock triple-points. Solutions for the 1E nozzle configuration were steady. Conversely, a subset of the 1F test conditions were unsteady and exhibited periodic, non-sinsoidal expansions and contractions of the streamwise-oriented plume, resulting in unsteady aerodynamic loads on the vehicle.

Supersonic Retropropulsion↗

Numerical simulations of unsteady, viscous, transonic flow over isolated and cascaded airfoils using a deforming grid

A compressible, unsteady, full Navier-Stokes, finite difference code was developed for modeling transonic flow through two-dimensional, oscillating cascades. The procedure introduces a deforming grid technique to capture the motion of the airfoils. Results using a deforming grid are presented for both isolated and cascaded airfoils. The load histories and unsteady pressure distributions are predicted for the NASA 64A010 isolated airfoil and compared with existing experimental data. Results show that the deforming grid technique can be used to successfully predict the unsteady flow properties around an oscillating airfoil. The deforming grid technique was extended for modeling unsteady flow in a cascade. The use of a deforming grid simplifies the specification of boundary conditions. Unsteady flow solutions similar to the isolated airfoil predictions are found for a NACA 0012 cascade with zero interblade phase angle and zero stagger. Experimental data for these cases are not available for code validation, but computational results are presented to show sample predictions from the code. Applications of the code to typical turbomachinery flow conditions will be presented in future work.

Huff, Dennis L.↗

Numerical simulations of unsteady, viscous, transonic flow over isolated and cascaded airfoils using a deforming grid

A compressible, unsteady, full Navier-Stokes, finite difference code was developed for modeling transonic flow through two-dimensional, oscillating cascades. The procedure introduces a deforming grid technique to capture the motion of the airfoils. Results using a deforming grid are presented for both isolated and cascaded airfoils. The load histories and unsteady pressure distributions are predicted for the NASA 64A010 isolated airfoil and compared with existing experimental data. Results show that the deforming grid technique can be used to successfully predict the unsteady flow properties around an oscillating airfoil. The deforming grid technique was extended for modeling unsteady flow in a cascade. The use of a deforming grid simplifies the specification of boundary conditions. Unsteady flow solutions similar to the isolated airfoil predictions are found for a NACA 0012 cascade with zero interblade phase angle and zero stagger. Experimental data for these cases are not available for code validation, but computational results are presented to show sample predictions from the code. Applications of the code to typical turbomachinery flow conditions will be presented in future work.

Huff, Dennis L.↗

Numerical Simulation of Inlet Bleed with Circular Holes on Plate Under Shock-Wave/Boundary-Layer Interactions

A numerical study was performed to investigate the shock-wave/boundary-layer interactions on a flat plate with bleed through one or more circular holes that vent into a plenum. The bleed-hole patterns considered for the study include in-line multiple holes and staggered multiple-row holes that are configured to simulate the patterns used in inlet bleed systems of high performance aircraft. The focus of the study was to examine how the bleed through multiple holes affect bleed rate and the pressure and Mach number distributions. Since the bleed performance was found sensitive to the change in bleed conditions, a computational procedure was developed to give a good turnaround computational time for parametric studies involving changes in bleed hole geometry and the structure of shock-wave/boundary-layer flowfield. The procedure includes the grid-generation methodology and the flow simulation with solutions from the Navier-Stokes equations. The computational techniques permit analysis of complex bleed systems and make possible the investigation of a broader range of design variables associated with inlet bleed operation.

Chyu, Wei J.↗

Source Term Model for an Array of Vortex Generator Vanes

A source term model was developed for numerical simulations of an array of vortex generators. The source term models the side force created by a vortex generator being modeled. The model is obtained by introducing a side force to the momentum and energy equations that can adjust its strength automatically based on a local flow. The model was tested and calibrated by comparing data from numerical simulations and experiments of a single low-profile vortex generator vane, which is only a fraction of the boundary layer thickness, over a flat plate. The source term model allowed a grid reduction of about seventy percent when compared with the numerical simulations performed on a fully gridded vortex generator without adversely affecting the development and capture of the vortex created. The source term model was able to predict the shape and size of the stream wise vorticity and velocity contours very well when compared with both numerical simulations and experimental data.

Buning, P. G.↗

Numerical Dissipation and Wrong Propagation Speed of Discontinuities for Stiff Source Terms

In compressible turbulent combustion/nonequilibrium flows, the constructions of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and (b) obtaining correct propagation speed of discontinuities for stiff reacting terms on coarse grids share one important ingredient - minimization of numerical dissipation while maintaining numerical stability. Here coarse grids means standard mesh density requirement for accurate simulation of typical non-reacting flows. This dual requirement to achieve both numerical stability and accuracy with zero or minimal use of numerical dissipation is most often conflicting for existing schemes that were designed for non-reacting flows. The goal of this paper is to relate numerical dissipations that are inherited in a selected set of high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities as a function of stiffness of the source term and the grid spacing.

Yee, H. C.↗

Multiphase Large Eddy Simulation of a Pulsed Sieve-Plate Extraction Column - 18125

Presented here is research resulting from a framework of single-phase and multiphase modelling trials. These are intended to accurately simulate the internal flow phenomena arising in pulsed sieve-plate extraction columns representative of those found in industry. Three-dimensional unsteady flow calculations using large eddy simulation coupled with dynamic sub-grid scale modelling have been used to capture the transient effects of sinusoidal flow conditions using the open source CFD code OpenFOAM{sup R}. Multiphase interface interactions are modelled using the volume of fluid method with appropriate heavy and light phase fluidic properties capturing surface tension, density and viscous effects. The system studied consists of a 150 mm diameter column with two perforated plates with a fractional free area of 25 %, operational under 1 Hz pulsing frequency at a 50 mm amplitude and plate spacing of 300 mm. The materials simulated include water (single-phase) and 3 M nitric acid and dodecane / tributyl phosphate mixture at 30 vol% (multiphase). Results are compared with appropriate correlations from literature. The dispersive mixing efficiency of the system was also evaluated via a mixing index value which allows quantification of the ratio of rotational and irrotational flow components. This study has found that the predictive methods employed are effective at capturing the different scales of turbulence present within pulsed sieve-plate extraction columns, and its influence on their operation. As such, the work highlights the ability of large eddy simulation-based modelling to predict the complex flows in such column geometries. Furthermore, the techniques employed provide valuable insight into the hydrodynamic flow characteristics and mixing mechanisms arising within pulsed sieve-plate extraction columns. (authors)

42 ENGINEERING↗

Impacts of Control, Penetration, and Distribution of Embedded Storage Network in Bulk Power System

The current shift in generation mix from fossil fuel plants towards variable and intermittent renewable energy sources is poised to create a future grid with reduced physical inertia and mismatch between generation and demand. Embedded storage, which is a concept of a coordinated network of storage units sited at the interface between the transmission and distribution system, is proposed as a mechanism to provide a buffer between generation and demand. This paper proposes an automated framework to model and integrate embedded storage in large-scale power systems with industry-grade grid-following (GFL) and grid-forming (GFM) control technologies. More importantly, the developed framework is used to explore the impacts of embedded storage control, penetration, location, and capacity in providing fast frequency response to the grid under contingency events such as generator trips and faults. The framework and study are conducted using the transient-stability simulation tool PSS/E and a realistic model of the Puerto Rico grid as a chosen test system. The simulation results show that GFL and GFM embedded storage, distributed throughout the system, with sufficient penetration and capacity, can effectively improve primary frequency response of the system under the studied contingency events.

Battery Energy Storage, embedded storage, grid-for↗

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges

Abstract High-fidelity simulation of transitional and turbulent flows over multi-scale surface roughness presents several challenges. For instance, the complex and irregular geometrical nature of surface roughness makes it impractical to employ conforming structured grids, commonly adopted in large-scale numerical simulations due to their high computational efficiency. One possible solution to overcome this problem is offered by immersed boundary methods, which allow wall boundary conditions to be enforced on grids that do not conform to the geometry of the solid boundary. To this end, a three-dimensional, second-order accurate boundary data immersion method (BDIM) is adopted. A novel mapping algorithm that can be applied to general three-dimensional surfaces is presented, together with a newly developed data-capturing methodology to extract and analyze on-surface flow quantities of interest. A rigorous procedure to compute gradient quantities such as the wall shear stress and the heat flux on complex non-conforming geometries is also introduced. The new framework is validated by performing a direct numerical simulation (DNS) of fully developed turbulent channel flow over sinusoidal egg-carton roughness in a minimal-span domain. For this canonical case, the averaged streamwise velocity profiles are compared against results from the literature obtained with a body-fitted grid. General guidelines on the BDIM resolution requirements for multi-scale roughness simulation are given. Momentum and energy balance methods are used to validate the calculation of the overall skin friction and heat transfer at the wall. The BDIM is then employed to investigate the effect of irregular homogeneous surface roughness on the performance of an LS-89 high-pressure turbine blade at engine-relevant conditions using DNS. This is the first application of the BDIM to realize multi-scale roughness for transitional flow in transonic conditions in the context of high-pressure turbines. The methodology adopted to generate the desired roughness distribution and to apply it to the reference blade geometry is introduced. The results are compared to the case of an equivalent smooth blade.

Engineering↗

Budget-Constrained Sizing of Renewable and Energy Storage Systems for Farm-Scale Ammonia Production Within the Food-Energy-Water Nexus

In the transition toward sustainable agriculture, farms have emerged as eco-friendly pioneers, harnessing clean hybrid wind and solar systems to improve farm performance. A concern in this paradigm is the effective sizing of renewable energy systems to ensure optimal energy use within budget considerations. This research focuses on optimizing renewable energy sizing in small-scale ammonia production to meet specific farm demands and enhance local resilience, emphasizing the interplay between environmental and economic factors. These findings promise increased energy efficiency and sustainability in this innovative agricultural sector. Additionally, our approach considers small-scale ammonia plant needs and the dynamic relationships between ammonia, water, and farm demands. Simulations demonstrate substantial cost savings in farm electricity consumption. Specifically, scenarios with renewable energy integration in the farm can reduce at least 13% electricity cost compared to a grid-dependent system in the 15-year simulation.

ammonia↗

Electron attenuation in free, neutral ethane clusters

The electron effective attenuation length (EAL) in free, neutral ethane clusters has been determined at 40 eV kinetic energy by combining carbon 1s x-ray photoelectron spectroscopy and theoretical lineshape modeling. More specifically, theory is employed to form model spectra on a grid in cluster size (N) and EAL (λ), allowing N and λ to be determined by optimizing the goodness-of-fit χ{sup 2}(N, λ) between model and observed spectra. Experimentally, the clusters were produced in an adiabatic-expansion setup using helium as the driving gas, spanning a range of 100–600 molecules in mean cluster size. The effective attenuation length was determined to be 8.4 ± 1.9 Å, in good agreement with an independent estimate of 10 Å formed on the basis of molecular electron-scattering data and Monte Carlo simulations. The aggregation state of the clusters as well as the cluster temperature and its importance to the derived EAL value are discussed in some depth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A priori mesh quality metrics for three-dimensional hybrid grids

Use of general hybrid grids to attain complex-geometry field simulations poses a challenge on estimation of their quality. Apart from the typical problems of non-uniformity and non-orthogonality, the change in element topology is an extra issue to address. The present work derives and evaluates an a priori mesh quality indicator for structured, unstructured, as well as hybrid grids consisting of hexahedra, prisms, tetrahedra, and pyramids. Emphasis is placed on deriving a direct relation between the quality measure and mesh distortion. The work is based on use of the Finite Volume discretization for evaluation of first order spatial derivatives. The analytic form of the truncation error is derived and applied to elementary types of mesh distortion including typical hybrid grid interfaces. The corresponding analytic expressions provide relations between the truncation error and the degree of stretching, skewness, shearing, torsion, expansion, as well as the type of grid interface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗