Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid code”

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,747 records · Page 97

The development of a CFD potential method for the analysis of tilt-rotors

A new CFD potential code, FPX (eXtended Full-Potential), has been developed for application to both helicopters and tilt-rotors. The code solves the unsteady, three-dimensional full potential equation and is an extension of the rotor code, FPR. Both entropy and viscosity corrections are included to enhance the physical modeling capabilities. A number of efficiency related modifications have yielded a factor of two speed-up in the code. An axial flow capability has been added to treat tilt-rotor in forward flight (cruise mode). In order to employ streamwise periodicity and accurately solve for the propagation of acoustic signals in the tip region, an H-H topology has been added to the basic O-H grid system. Computations are performed for the XV-15 Standard and ATB blades at high-speed conditions. Comparisons are made for the blade aerodynamics and the induced fuselage cabin pressure for a range of Mach numbers. Grid generation, wake treatment, and far-field wall treatment are identified as problem areas with recommendations for future research.

Bridgeman, John O.↗

Progress In Incompressible Pump Flow Calculations

Steady and unsteady flows for propulsion systems are efficiently simulated by solving the incompressible Navier-Stokes equations. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Current computations use one equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard kappa - epsilon model equations. The resulting computer code is applied to the flow analysis inside an advanced rocket pump impeller in steadily rotating reference frames. Numerical results are compared with experimental measurements. The effects of exit and shroud cavities with the leak-age flow are investigated. Time-accurate incompressible Navier-Stokes formulation with the overlapped grid scheme capability was evaluated by using MIT flapping foil experiment. The grid dependency, turbulence model effects, and the effect of order of differencing were investigated. Numerical results were compared against experimental data. The resulting procedure were applied to unsteady flapping foil calculations. Two upstream NACA 0025 foils perform high frequency synchronized motion and generate unsteady flow conditions to the downstream larger stationary foil. Comparison between unsteady experimental data and numerical results from two different moving boundary procedures will be presented.

Kiris, Cetin↗

Comparison of Numerical Schemes for a Realistic Computational Aeroacoustics Benchmark Problem

In this work, a nonlinear structured-multiblock CAA solver, the NASA GRC BASS code, will be tested on a realistic CAA benchmark problem. The purpose of this test is to ascertain what effect the high-accuracy solution methods used in CAA have on a realistic test problem, where both the mean flow and the unsteady waves are simultaneously computed on a fully curvilinear grid from a commercial grid generator. The proposed test will compare the solutions obtained using several finite-difference methods on identical grids to determine whether high-accuracy schemes have advantages for this benchmark problem.

Hixon, R.↗

Three-Dimensional Analysis and Modeling of a Wankel Engine

A new computer code, AGNI-3D, has been developed for the modeling of combustion, spray, and flow properties in a stratified-charge rotary engine (SCRE). The mathematical and numerical details of the new code are described by the first author in a separate NASA publication. The solution procedure is based on an Eulerian-Lagrangian approach where the unsteady, three-dimensional Navier-Stokes equations for a perfect gas-mixture with variable properties are solved in generalized, Eulerian coordinates on a moving grid by making use of an implicit finite-volume, Steger-Warming flux vector splitting scheme. The liquid-phase equations are solved in Lagrangian coordinates. The engine configuration studied was similar to existing rotary engine flow-visualization and hot-firing test rigs. The results of limited test cases indicate a good degree of qualitative agreement between the predicted and measured pressures. It is conjectured that the impulsive nature of the torque generated by the observed pressure nonuniformity may be one of the mechanisms responsible for the excessive wear of the timing gears observed during the early stages of the rotary combustion engine (RCE) development. It was identified that the turbulence intensities near top-dead-center were dominated by the compression process and only slightly influenced by the intake and exhaust processes. Slow mixing resulting from small turbulence intensities within the rotor pocket and also from a lack of formation of any significant recirculation regions within the rotor pocket were identified as the major factors leading to incomplete combustion. Detailed flowfield results during exhaust and intake, fuel injection, fuel vaporization, combustion, mixing and expansion processes are also presented. The numerical procedure is very efficient as it takes 7 to 10 CPU hours on a CRAY Y-MP for one entire engine cycle when the computations are performed over a 31 x16 x 20 grid.

Raju, M. S.↗

Numerical Simulation of Dynamic Stall Using Near-Body Adaptive Mesh Refinement

Time-dependent Reynolds-averaged Navier-Stokes simulations have been carried out for a NACA 0012 airfoil and a flexible UH-60A rotor undergoing dynamic stall. The OVERFLOW computational fluid dynamics (CFD) code is used to explore the use of near-body adaptive mesh refinement (NB-AMR) for the first time on a flexible helicopter rotor in forward flight. Emphasis is placed on understanding the temporal and spatial convergence of the solutions, their numerical efficiency and stability, and establishing a grid-converged solution. Some of the remaining differences between CFD and flight-test measurements are discussed. Time-dependent flow visualization is used to provide an improved understanding of the physical mechanisms involved with two-dimensional and three-dimensional dynamic stall, and the NB-AMR process.

Near-Body↗

Resonance Self-shielding Impact on Neutron Spectrum Determination for Missouri S and T Reactor

Neutron Activation Analysis is a useful technique for characterizing the neutron flux at nuclear facilities. It involves irradiating samples to produce unstable isotopes, which enables qualitative and quantitative analyses of unknown samples. In the process, samples are bombarded with neutrons; new elements then become radioactive isotopes due to absorption of neutrons. These emit different amounts of gamma energy during the decay process, and gamma rays are counted using HPGe. Each reaction product has a known half-life and cross-section, and by measuring the amount of gamma radiation energy emitted from a radioactive isotope, it is possible to identify the isotope and activity. However, the neutron spectrum is not precise or representative of the actual neutron environment because of the self-shielding effect. This occurs with large samples that have large cross-sections or thicknesses in the resonance range of the neutron flux spectrum. Several methods have been used to mitigate this problem, such as thinning or diluting alloys samples, calculating the self-shielding factor for each material, and generating a corrected library for the unfolding procedure code. The Missouri University of Science and Technology Reactor (MSTR) is a pool-type research reactor that operates up to 200 kW. It uses light-water for both moderation and natural convective heat removal. The reactor core consists of 15 fuel elements positioned in a 9 x 6 aluminum grid plate, with four elements used for reactor power control. The fuel element is 19.75% low enriched uranium in 18 curved fuel plates. The reactor has seen various configurations starting with configuration designation 101W and subsequent reconfiguration increases the numeric designation by 1. The current configuration, 120W, was implemented in 2010. The new configuration provides more in-core locations for activation experiments and uses more fuel elements than previous ones. One of the irradiation locations in the MSTR core is the source holder. It is desired to improve the documentation of neutronic characteristics of the MSTR which are cardinal to experiments supported by the facility. One of these characteristics is the neutron spectrum at various locations in the reactor. Successful efforts yielded the neutron spectrum at a known location through the unfolding of the activity of irradiated foils. The purpose of this current study is to investigate the impact of resonance self-shielding of some of the activated foils on the determined spectrum.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Demonstration of a Novel Technology to Manage Electricity Demand in Grid-Independent Military Microgrids

This research was conducted by the National Renewable Energy Laboratory (NREL) in collaboration with the S&C Electric Inc. through funding provided by the ESTCP. The project demonstrates use of cybersecure Automated Demand Response (ADR) technology to effectively manage microgrid loads during grid-independent, also known as "islanded," operation. When military microgrids become isolated from the main electrical grid, they are required to balance electricity supply and demand locally. Given that local generation may be constrained, the prevailing strategy involves shedding all but the most critical loads by tripping smart circuit breakers, which then necessitate manual resetting. This approach is generally implemented at the building level, which means that the buildings with mission-critical activities are exempt from load management and remain fully powered, whereas those deemed non-critical can experience a complete loss of service. In this research we developed a method that allows building automation systems to selectively control their assets in response to load shedding request from a microgrid controller, avoiding total loss of service in contrast to the conventional control approach. A commercial OpenADR client server by GridFabric is used for communication between the microgrid controller and the building management system (BMS). The microgrid controller monitors both generation capacity and various assets within the microgrid and issues a demand reduction request when necessary. This request is communicated to the OpenADR server via Modbus. Upon receiving the request, the OpenADR server forwards it to the BMS utilizing the OpenADR protocol. The BMS is pre-configured with various levels of load reduction strategies based on the controllable assets available, allowing for a nuanced approach to demand reduction. Both lab and field tests were performed that considered load shedding needed to achieve closed transition into island mode and to accommodate changing loads and power source availability while islanded. A commercial microgrid controller was used for these tests with normal programming within the expected constraints of the system capabilities. That is, the solution did not require any specialized modification to the code base of the controller. Given the latency of the round-trip communication path between the microgrid controller and the various devices involved with the load shed processes, there are certain scenarios for which the demonstrated solution are appropriate and some which are not. The methods described in this report can be used for load shedding/restoration during transitions between islanded and grid-tied modes of operation, as well as accommodating normal variations in load and the need to remove a power source from operation for maintenance. These methods should not be used for scenarios that require load shedding within a second or two such as sudden and unanticipated significant load increases or loss of power sources through equipment faults.

24 POWER TRANSMISSION AND DISTRIBUTION↗

SOC Synthetic Microstructure Bank

QUICK START: Start with property_library.html (can be found by typing the filename into the query box) and use the interactive table to filter, sort, and select a microstructure with the desired properties. Search for the alphabetic code to obtain the corresponding dataset. Full description: This is a bank of 1,970 unique 3-phase electrode microstructure files. When you account for reassigning phase IDs (e.g. declare that 1=Ni and 2=pore, instead of 1=pore and 2=Ni), it actually represents 5,910 unique electrode microstructures, each of which could be considered to be either an air or a fuel electrode (e.g. declare that the phase IDs correspond to pore, Ni, and YSZ; or that they correspond to pore, LSCF, and GDC; or whatever electron-conductor and ion-conductor combination is being studied). The voxel size is 50 nm and each electrode file contains a 4x4 grid of (12.5 micron)^3 sub-volumes. If placed together in a grid, they comprise a 50x50x12.5 micron electrode (note that the interfaces between sub-volumes will be sharp; this can be mitigated via simulating annealing/relaxation). The sub-volumes can also be used individually for a reasonably sized 12.5 micron cubic region-of-interest. A user can simply consider the voxel size to be a different value to rescale the volumes (and all of their morphological features, including particle size) as desired. These microstructures were generated using DREAM3D. The general procedure is outlined in https://doi.org/10.1016/j.jpowsour.2018.03.025 The file names are an alphabetic code having to do with the input parameters used in DREAM3D when they were generated. Most users would be best served by starting with the file property_library.html or property_library_subvols.html (which lists properties for each individual subvolume). These files contain a catalogue of the actual, measured properties of every microstructure in the database. Any combination of property values can be filtered and sorted until a desired electrode is found, at which point the user can find the file corresponding to that alphabetic code. The properties in the catalogue include connected TPB density, and for each phase: phase fraction, average particle size, polydispersity of particle size, tortuosity, and connected pair-wise interfacial area. They also include what fraction of each property is connected through to the interfaces of the volume. If the desired combination of properties is not found at first, remember that the phase IDs can be re-assigned arbitrarily, e.g. swapping 1s and 2s. In fact, the database was generated with this in mind so as not to generate redundant microstructures. If the database does not contain the desired property combinations, try to search for the other possible permutations of those properties with re-assigned phase IDs. Please cite https://doi.org/10.1149/10301.0909ecst for use. Please contact the maintainer, William K. Epting, for additional information or assistance.

3D microstructure↗

Highly Available Nuclear Power for Mission-Critical Applications

Some energy consumers require power on an anytime, all-year-round basis with a high level of certainty, including defense installations, isolated communities, and some industrial processes. For these customers, interruptions in electricity or heat can mean substantial financial loss or even loss of life. In the absence of grid-scale energy storage, a high level of power availability can be accomplished only through the robustness and redundancy of power generators. The NuScale small modular reactor design is well suited to provide highly available power because of several features related to both the nuclear steam supply system and the overall plant design. In analogy to Redundant Array of Independent Disks (RAID) systems used to provide highly reliable data storage, a NuScale plant can assure sustained power generation by virtue of its Redundant Array of Integral Reactors (RAIR).This paper describes the NuScale RAIR plant features and summarizes the results of a rigorous analysis of RAIR availability as a function of power or, conversely, the RAIR plant output power as a function of power availability. The analysis utilized the Matrix Laboratory code (MATLAB) and included probability distributions for the frequency and duration of module outages due to planned and unplanned events. Further, the study evaluated the impact of implementing turbine bypass rather than module shutdown and using one or more modules to supply house loads in the case of loss of off-site power (LOOP). Availability results are presented for a 12-module RAIR plant with and without turbine bypass enabled during a LOOP and for different possible connections to the off-site power distribution grid and dedicated service loads. Results indicate that a very high level of availability can be achieved at relatively high power output levels, regardless of turbine bypass and dedicated load connection, compared to the operating fleet.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part I: Lagrangian and indirect Euler AMR algorithms

Many applications of physics and engineering involve wide ranges of time and spatial scales. The numerical simulation of localized small scales such as shock waves and material interfaces requires a large number of computational cells in these regions. For these applications, Lagrangian and Arbitrary-Lagrangian-Eulerian (ALE) related methods are engaging since the moving mesh feature naturally brings mesh cells on shock discontinuities and material interfaces are carefully captured. In addition, Adaptive-Mesh-Refinement (AMR) strategies aim to optimize computational resources by concentrating finer mesh cells only in areas of interest while using coarser cells elsewhere. A key but challenging AMR requirement consists in efficiently distributing the computational effort to achieve high accuracy without the prohibitive computational costs associated with uniformly fine grids. Here, in this document, the coupling of the p4est AMR library with a cell-centered Lagrangian scheme is presented with the goal to perform reliable 3D Lagrangian-AMR and indirect Euler-AMR multi-material simulations. In particular, it is shown that starting from a 3D indirect ALE code, the memory management and load balancing requirements can be delegated to an external library (here the p4est library) to unlock ALE-AMR capabilities. First, we present a strategy to transcribe the octant-based connectivity of the 3D AMR framework with that of an unstructured mesh of polygonal cells used in Lagrangian hydrodynamics. Then, we show how refinement and coarsening operations must be adapted to the particular Lagrangian framework to ensure the conservation of volume during those steps. Finally, several numerical test cases are presented that demonstrate the capabilities of the Lagrangian-AMR and indirect Euler-AMR algorithms.

3D cell-centered Lagrangian numerical scheme↗

Decoding Resilience by Modeling Outage and Restoration Processes in Distribution Grids: A Pittsburgh Case Study

Climate-induced extreme weather events, such as floods and heatwaves, pose significant challenges to the resilience of urban power distribution grids. This paper examines the outage and restoration dynamics of Pittsburgh's power grid during the flood event in April 2024 and three heatwaves in June, July, and August 2024. We introduce a comprehensive modeling framework that integrates outage and restoration processes, enabling the quantification of resilience metrics, including total customer-hours of power outage, maximum residual values, and restoration durations across various ZIP codes. Our analysis highlights spatial disparities in outage impacts and restoration efficiencies, with ZIP Code 15222 experiencing the highest cumulative disruptions, while ZIP Code 15217 shows lower susceptibility. By linking statistical trends with weather events, this study underscores the critical need for targeted infrastructure upgrades and advanced restoration strategies. Furthermore, the proposed framework offers valuable insights for planning and managing resilient power systems in the face of increasing climate stress.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Particle-in-cell Simulations of Relativistic Magnetic Reconnection with Advanced Maxwell Solver Algorithms

Abstract Relativistic magnetic reconnection is a nonideal plasma process that is a source of nonthermal particle acceleration in many high-energy astrophysical systems. Particle-in-cell (PIC) methods are commonly used for simulating reconnection from first principles. While much progress has been made in understanding the physics of reconnection, especially in 2D, the adoption of advanced algorithms and numerical techniques for efficiently modeling such systems has been limited. With the GPU-accelerated PIC code WarpX, we explore the accuracy and potential performance benefits of two advanced Maxwell solver algorithms: a nonstandard finite-difference scheme (CKC) and an ultrahigh-order pseudo-spectral method (PSATD). We find that, for the relativistic reconnection problem, CKC and PSATD qualitatively and quantitatively match the standard Yee-grid finite-difference method. CKC and PSATD both admit a time step that is 40% longer than that of Yee, resulting in a ∼40% faster time to solution for CKC, but no performance benefit for PSATD when using a current deposition scheme that satisfies Gauss’s law. Relaxing this constraint maintains accuracy and yields a 30% speedup. Unlike Yee and CKC, PSATD is numerically stable at any time step, allowing for a larger time step than with the finite-difference methods. We found that increasing the time step 2.4–3 times over the standard Yee step still yields accurate results, but it only translates to modest performance improvements over CKC, due to the current deposition scheme used with PSATD. Further optimization of this scheme will likely improve the effective performance of PSATD.

79 ASTRONOMY AND ASTROPHYSICS↗

County Level Annual Population Projections for SSP3 and SSP5

This dataset consist of annual (2020-2100) county-level population projections for the United States (U.S.) for Shared Socioeconomic Pathways (SSPs) 3 and 5. The original decadal state-level data is used as an input to the gridded population data model to downscale the state-level projections for each SSP to a 1 km grid. The 1 km data is then aggregated to the county-level using the 2020 TIGER/Line county shapefiles from the U.S. Census Bureau. The two data files are shared as .csv files with the following structure: Rows = Data for individual counties which are identified using their Federal Information Processing Standard (FIPS) code. The FIPS code is stored in the first column. Columns = Each year from 2020-2100 is an individual column. For easier reference, the state name associated with each row is stored in the last column. Values in each cell are the downscaled estimated population for that county-year combination. Values are fractional due to the weighting scheme used in the downscaling. The paper and model source code cited in the "Related Works" section describes the initial source of the population projections.

FIPS↗

Boundary-fitted coordinate systems for arbitrary computational regions

The method of Smith and Wiegel was used to generate meshes for mixer lobes and subsonic inlets that are compatible with flow analysis codes requiring a boundary fitted coordinate system. Successful application of this mesh generator required development of procedures to distribute the mesh points along the boundaries, to regulate the dependence of the connecting function to the local boundary slope, to concentrate the mesh into regions of special interest, and to modify the mesh grid so that it possessed a smooth progression of cell metrics and cell volumes in all directions. The method of Smith and Wiegel when coupled with the extensions mentioned above has proven to be easy to use and control for the inlet and mixer lobe geometries investigated.

Kowalski, E. J.↗

Transonic flutter calculations using the Euler equations

In transonic flutter problems where shock motion plays an important part, it is believed that accurate predictions of the flutter boundaries will require the use of codes based on the Euler equations. Only Euler codes can obtain the correct shock location and shock strength, and the crucially important shock excursion amplitude and phase lag. The present study is based on the finite volume scheme developed by Jameson and Venkatakrishnan for the 2-D unsteady Euler equations. The equations are solved in integral form on a moving grid. The variable are pressure, density, Cartesian velocity components, and total energy.

Bendiksen, Oddvar O.↗

An efficient DSMC algorithm applied to a delta wing

A new algorithm for 3D direct simulation Monte Carlo (DSMC) is tested and numerical results are compared with wind tunnel data and results obtained earlier with a more traditional DSMC code. The test case is the flowfield around a delta wing at incidence at Knudsen number of 0.016 and Mach number of 20.2. The results are shown to compare favorably with both experimental and earlier numerical results. The new algorithm is described with special emphasis placed on its distinctive features: Cartesian/unstructured combination grid, special body surface definition, discretization in physical space.

Rault, D. F. G.↗

Multigrid methods for flow transition in three-dimensional boundary layers with surface roughness

The efficient multilevel adaptive method has been successfully applied to perform direct numerical simulations (DNS) of flow transition in 3-D channels and 3-D boundary layers with 2-D and 3-D isolated and distributed roughness in a curvilinear coordinate system. A fourth-order finite difference technique on stretched and staggered grids, a fully-implicit time marching scheme, a semi-coarsening multigrid method associated with line distributive relaxation scheme, and an improved outflow boundary-condition treatment, which needs only a very short buffer domain to damp all order-one wave reflections, are developed. These approaches make the multigrid DNS code very accurate and efficient. This allows us not only to be able to do spatial DNS for the 3-D channel and flat plate at low computational costs, but also to do spatial DNS for transition in the 3-D boundary layer with 3-D single and multiple roughness elements, which would have extremely high computational costs with conventional methods. Numerical results show good agreement with the linear stability theory, the secondary instability theory, and a number of laboratory experiments. The contribution of isolated and distributed roughness to transition is analyzed.

Liu, Chaoqun↗

Coupled Ablation, Heat Conduction, Pyrolysis, Shape Change and Spallation of the Galileo Probe

The Galileo probe enters the atmosphere of Jupiter in December 1995. This paper presents numerical methodology and detailed results of our final pre-impact calculations for the heat shield response. The calculations are performed using a highly modified version of a viscous shock layer code with massive radiation coupled with a surface thermochemical ablation and spallation model and with the transient in-depth thermal response of the charring and ablating heat shield. The flowfield is quasi-steady along the trajectory, but the heat shield thermal response is dynamic. Each surface node of the VSL grid is coupled with a one-dimensional thermal response calculation. The thermal solver includes heat conduction, pyrolysis, and grid movement owing to surface recession. Initial conditions for the heat shield temperature and density were obtained from the high altitude rarefied-flow calculations of Haas and Milos. Galileo probe surface temperature, shape, mass flux, and element flux are all determined as functions of time along the trajectory with spallation varied parametrically. The calculations also estimate the in-depth density and temperature profiles for the heat shield. All this information is required to determine the time-dependent vehicle mass and drag coefficient which are necessary inputs for the atmospheric reconstruction experiment on board the probe.

Milos, Frank S.↗