Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid codes”

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems. Task 5: Unsteady counterrotation ducted propfan analysis. Computer program user's manual

The primary objective of this study was the development of a time-marching three-dimensional Euler/Navier-Stokes aerodynamic analysis to predict steady and unsteady compressible transonic flows about ducted and unducted propfan propulsion systems employing multiple blade rows. The computer codes resulting from this study are referred to as ADPAC-AOACR (Advanced Ducted Propfan Analysis Codes-Angle of Attack Coupled Row). This report is intended to serve as a computer program user's manual for the ADPAC-AOACR codes developed under Task 5 of NASA Contract NAS3-25270, Unsteady Counterrotating Ducted Propfan Analysis. The ADPAC-AOACR program is based on a flexible multiple blocked grid discretization scheme permitting coupled 2-D/3-D mesh block solutions with application to a wide variety of geometries. For convenience, several standard mesh block structures are described for turbomachinery applications. Aerodynamic calculations are based on a four-stage Runge-Kutta time-marching finite volume solution technique with added numerical dissipation. Steady flow predictions are accelerated by a multigrid procedure. Numerical calculations are compared with experimental data for several test cases to demonstrate the utility of this approach for predicting the aerodynamics of modern turbomachinery configurations employing multiple blade rows.

Hall, Edward J.↗
NASA NTRS1993

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems. Task 5: Unsteady counterrotation ducted propfan analysis

The primary objective of this study was the development of a time-marching three-dimensional Euler/Navier-Stokes aerodynamic analysis to predict steady and unsteady compressible transonic flows about ducted and unducted propfan propulsion systems employing multiple blade rows. The computer codes resulting from this study are referred to as ADPAC-AOAR\CR (Advanced Ducted Propfan Analysis Codes-Angle of Attack Coupled Row). This document is the final report describing the theoretical basis and analytical results from the ADPAC-AOACR codes developed under task 5 of NASA Contract NAS3-25270, Unsteady Counterrotating Ducted Propfan Analysis. The ADPAC-AOACR Program is based on a flexible multiple blocked grid discretization scheme permitting coupled 2-D/3-D mesh block solutions with application to a wide variety of geometries. For convenience, several standard mesh block structures are described for turbomachinery applications. Aerodynamic calculations are based on a four-stage Runge-Kutta time-marching finite volume solution technique with added numerical dissipation. Steady flow predictions are accelerated by a multigrid procedure. Numerical calculations are compared with experimental data for several test cases to demonstrate the utility of this approach for predicting the aerodynamics of modern turbomachinery configurations employing multiple blade rows.

Hall, Edward J.↗
NASA NTRS2000

Full 3D Analysis of the GE90 Turbofan Primary Flowpath

The multistage simulations of the GE90 turbofan primary flowpath components have been performed. The multistage CFD code, APNASA, has been used to analyze the fan, fan OGV and booster, the 10-stage high-pressure compressor and the entire turbine system of the GE90 turbofan engine. The code has two levels of parallel, and for the 18 blade row full turbine simulation has 87.3 percent parallel efficiency with 121 processors on an SGI ORIGIN. Grid generation is accomplished with the multistage Average Passage Grid Generator, APG. Results for each component are shown which compare favorably with test data.

Turner, Mark G.↗
NASA NTRS1999

Helicon Plasma Source Optimization Studies for VASIMR

A helicon plasma source at Oak Ridge National Laboratory is being used to investigate operating scenarios relevant to the VASIMR (VAriable Specific Impulse Magnetoplasma Rocket). These include operation at high magnetic field (> = 0.4 T), high frequency (<= 30 MHz), high power (< = 3 kW), and with light ions (He+, H+). To date, He plasmas have been produced with n(sub e0) = 1.7 x 10(exp 19)/cu m (measured with an axially movable 4mm microwave interferometer), with Pin = I kW at f = 13.56 MHz and absolute value of B(sub 0) = 0.16 T. In the near future, diagnostics including a mass flow meter and a gridded energy analyzer array will be added to investigate fueling efficiency and the source power balance. The latest results, together with modeling results using the EMIR rf code, will be presented.

Goulding, R. H.↗
NASA NTRS2006

Multi-Megawatt Gas Turbine Power Systems for Lunar Colonies

A concept for development of second generation 10 MWe prototype lunar power plant utilizing a gas cooled fission reactor supplying heated helium working fluid to two parallel 5 MWe closed cycle gas turbines is presented. Such a power system is expected to supply the energy needs for an initial lunar colony with a crew of up to 50 persons engaged in mining and manufacturing activities. System performance and mass details were generated by an author developed code (BRMAPS). The proposed pilot power plant can be a model for future plants of the same capacity that could be tied to an evolutionary lunar power grid.

Juhasz, Albert J.↗
NASA NTRS2002

Isotopic Effects in Nuclear Fragmentation and GCR Transport Problems

Improving the accuracy of the galactic cosmic ray (GCR) environment and transport models is an important goal in preparing for studies of the projected risks and the efficiency of potential mitigations methods for space exploration. In this paper we consider the effects of the isotopic composition of the primary cosmic rays and the isotopic dependence of nuclear fragmentation cross sections on GCR transport models. Measurements are used to describe the isotopic composition of the GCR including their modulation throughout the solar cycle. The quantum multiple-scattering approach to nuclear fragmentation (QMSFRG) is used as the data base generator in order to accurately describe the odd-even effect in fragment production. Using the Badhwar and O'Neill GCR model, the QMSFRG model and the HZETRN transport code, the effects of the isotopic dependence of the primary GCR composition and on fragment production for transport problems is described for a complete GCR isotopic-grid. The principle finding of this study is that large errors ( 100%) will occur in the mass-flux spectra when comparing the complete isotopic-grid (141 ions) to a reduced isotopic-grid (59 ions), however less significant errors 30%) occur in the elemental-flux spectra. Because the full isotopic-grid is readily handled on small computer work-stations, it is recommended that they be used for future GCR studies.

Cucinotta, Francis A.↗
NASA NTRS2015

Aeroacoustic Simulations of a Nose Landing Gear Using FUN3D on Pointwise Unstructured Grids

Numerical simulations have been performed for a partially-dressed, cavity-closed (PDCC) nose landing gear configuration that was tested in the University of Florida's open-jet acoustic facility known as the UFAFF. The unstructured-grid flow solver FUN3D is used to compute the unsteady flow field for this configuration. Mixed-element grids generated using the Pointwise(TradeMark) grid generation software are used for these simulations. Particular care is taken to ensure quality cells and proper resolution in critical areas of interest in an effort to minimize errors introduced by numerical artifacts. A hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) turbulence model is used for these simulations. Solutions are also presented for a wall function model coupled to the standard turbulence model. Time-averaged and instantaneous solutions obtained on these Pointwise grids are compared with the measured data and previous numerical solutions. The resulting CFD solutions are used as input to a Ffowcs Williams-Hawkings noise propagation code to compute the farfield noise levels in the flyover and sideline directions. The computed noise levels compare well with previous CFD solutions and experimental data.

Vatsa, Veer N.↗
NASA NTRS2018

Low Reynolds Number Airfoil Evaluation for the Mars Helicopter Rotor

The present research is aimed at providing a performance model for the Mars Helicopter (MH), to understand the complexity of the flow, and future regions of improvement. The Martian atmosphere's low density and the MH's relatively small rotor result in very low chord-based Reynolds number flows, Rec = O(10(exp 3)-10(exp 4)). The low density and subcritical Reynolds number reduce the lifting force and lifting efficiency, respectively. The high drag coefficients in subcritical flow, especially for thicker sections, are attributed to laminar separation from the rear of the airfoil. The goal is to generate a performance model for the MH rotor for a free wake analysis, since the computational budget for a complete Navier-Stokes solution for a rotating body-fitted rotor is substantial. In this study, a RANS-based approach is used to generate the airfoil deck using OVERFLOW with stitched experimental data for very high angles of attack. The model is presented through airfoil data tables (C81 files) that are used by comprehensive rotor analysis codes such as CAMRADII, or the mid-fidelity CFD solver RotCFD. These codes have proven to provide accurate performance predictions for all rotor operations at only a fraction of the computational expense of three- dimensional body-fitted viscous grids.

Koning, Witold J. K.↗
NASA NTRS1994

Automated Instrumentation, Monitoring and Visualization of PVM Programs Using AIMS

We present views and analysis of the execution of several PVM codes for Computational Fluid Dynamics on a network of Sparcstations, including (a) NAS Parallel benchmarks CG and MG (White, Alund and Sunderam 1993); (b) a multi-partitioning algorithm for NAS Parallel Benchmark SP (Wijngaart 1993); and (c) an overset grid flowsolver (Smith 1993). These views and analysis were obtained using our Automated Instrumentation and Monitoring System (AIMS) version 3.0, a toolkit for debugging the performance of PVM programs. We will describe the architecture, operation and application of AIMS. The AIMS toolkit contains (a) Xinstrument, which can automatically instrument various computational and communication constructs in message-passing parallel programs; (b) Monitor, a library of run-time trace-collection routines; (c) VK (Visual Kernel), an execution-animation tool with source-code clickback; and (d) Tally, a tool for statistical analysis of execution profiles. Currently, Xinstrument can handle C and Fortran77 programs using PVM 3.2.x; Monitor has been implemented and tested on Sun 4 systems running SunOS 4.1.2; and VK uses X11R5 and Motif 1.2. Data and views obtained using AIMS clearly illustrate several characteristic features of executing parallel programs on networked workstations: (a) the impact of long message latencies; (b) the impact of multiprogramming overheads and associated load imbalance; (c) cache and virtual-memory effects; and (4significant skews between workstation clocks. Interestingly, AIMS can compensate for constant skew (zero drift) by calibrating the skew between a parent and its spawned children. In addition, AIMS' skew-compensation algorithm can adjust timestamps in a way that eliminates physically impossible communications (e.g., messages going backwards in time). Our current efforts are directed toward creating new views to explain the observed performance of PVM programs. Some of the features planned for the near future include: (a) ConfigView, showing the physical topology of the virtual machine, inferred using specially formatted IP (Internet Protocol) packets; and (b) LoadView, synchronous animation of PVM-program execution and resource-utilization patterns.

Mehra, Pankaj↗
DOE OSTI2022

Experimental measurements of fluid flow in an 84-pin hexagonal rod bundle with spacer grid for a gas-cooled fast modular reactor

A 50 MW e helium-cooled fast modular reactor (FMR) is under development, and as part of the Department of Energy Integrated Research Project (IRP), Texas A&M University is conducting the thermal-hydraulic characterization of its baseline core configuration. Here, we experimentally investigated the axial and cross flow fields characteristics in the hexagonal fuel rod bundle composed of 84 rods, a central rod, and spacer grids at a Reynolds number of 12,000. A fully transparent experimental facility resembling of one unit of the fuel assembly was constructed. Time-resolved particle image velocimetry (TR-PIV) measurements were performed to characterize hydraulic behavior downstream the spacer grid. Velocity measurements were conducted in the axial and radial direction of the rod bundle. From the PIV velocity vector fields, the full-field flow statistics were computed for the mean velocity, vorticity, and Reynolds stresses. An analysis of the energy decay downstream of the spacer grid was conducted with calculations of secondary-flow intensities, turbulent kinetic energy, power spectrum, and spatial-temporal velocity cross correlations. The vorticity field was obtained and pairs of counter-rotating vortexes were identified using a 3D reconstruction of several measurement planes. The data from this experimental campaign will be used to validate numerical models based on Computational Fluid Dynamics and other codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Continuum gyrokinetic simulations of edge plasmas in single-null geometries

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. Furthermore, in order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

H-AMR: A New GPU-accelerated GRMHD Code for Exascale Computing with 3D Adaptive Mesh Refinement and Local Adaptive Time Stepping

General relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black hole accretion. Here, we present a GPU-accelerated GRMHD code H-AMR with multifaceted optimizations that, collectively, accelerate computation by 2–5 orders of magnitude for a wide range of applications. First, it introduces a spherical grid with 3D adaptive mesh refinement that operates in each of the three dimensions independently. This allows us to circumvent the Courant condition near the polar singularity, which otherwise cripples high-resolution computational performance. Second, we demonstrate that local adaptive time stepping on a logarithmic spherical-polar grid accelerates computation by a factor of ≲10 compared to traditional hierarchical time-stepping approaches. Jointly, these unique features lead to an effective speed of ~10 9 zone cycles per second per node on 5400 NVIDIA V100 GPUs (i.e., 900 nodes of the OLCF Summit supercomputer). We illustrate H-AMR's computational performance by presenting the first GRMHD simulation of a tilted thin accretion disk threaded by a toroidal magnetic field around a rapidly spinning black hole. With an effective resolution of 13,440 × 4608 × 8092 cells and a total of ≲22 billion cells and ~0.65 × 10 8 time steps, it is among the largest astrophysical simulations ever performed. We find that frame dragging by the black hole tears up the disk into two independently precessing subdisks. The innermost subdisk rotation axis intermittently aligns with the black hole spin, demonstrating for the first time that such long-sought alignment is possible in the absence of large-scale poloidal magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Ensemble-Based Characterization of Historical Storm Types in ERA5 over CONUS

We present Storm-Type Labeled Precipitation, a gridded dataset that classifies ERA5 precipitation over CONUS (20–50° N, 125–66° W) by storm type at 6-hourly resolution. Each grid cell/time step is assigned to one of five classes—mesoscale convective system (MCS), extratropical cyclone (ETC), hurricane (HUR), atmospheric river (AR), or other convective, with an “unidentified” code reserved for cases with no detected type. Two products are provided: (i) a TempestExtremes (TE)–only version, and (ii) a hybrid TRACK–TE version that uses TRACK for ETCs and HURs and TE for MCSs and ARs. Outputs are integer masks aligned to ERA5 precipitation, supporting storm-type attribution of rainfall, event compositing, trend analysis, and model evaluation. Future work will extend detection to high-resolution downscaled Earth system model projections to support evaluations of projected changes in hydrologic hazards that may threaten critical water and energy infrastructures.

Rastogi, Deeksha [ORNL] (ORCID:0000000204624027)↗
NASA NTRS1984

Implicit treatment of the unsteady full potential equation in conservation form

An implicit, conservative treatment for the unsteady full potential equation in two-dimensions is presented. The method employs a local time linearization for density, and introduces flux biasing concepts based on sonic conditions for the generation of artificial viscosity to capture shocks without any overshoots. The boundary condition is treated implicitly using a splitting procedure consistent with the approximate factorization scheme. This allows for extremely large Courant numbers, even for nonorthogonal grid at the body. The method has application not only to unsteady problems, but also to generate the starting blunt body solution for a supersonic full potential marching code. Results are presented for flows over cylinders, spheres and airfoils. Comparisons are made with available Euler and full potential results, and are in excellent agreement.

Shankar, V.↗
NASA NTRS1992

Pressure wave propagation studies for oscillating cascades

The unsteady flow field around an oscillating cascade of flat plates is studied using a time marching Euler code. Exact solutions based on linear theory serve as model problems to study pressure wave propagation in the numerical solution. The importance of using proper unsteady boundary conditions, grid resolution, and time step is demonstrated. Results show that an approximate non-reflecting boundary condition based on linear theory does a good job of minimizing reflections from the inflow and outflow boundaries and allows the placement of the boundaries to be closer than cases using reflective boundary conditions. Stretching the boundary to dampen the unsteady waves is another way to minimize reflections. Grid clustering near the plates does a better job of capturing the unsteady flow field than cases using uniform grids as long as the CFL number is less than one for a sufficient portion of the grid. Results for various stagger angles and oscillation frequencies show good agreement with linear theory as long as the grid is properly resolved.

Huff, Dennis L.↗
NASA NTRS1992

Pressure wave propagation studies for oscillating cascades

The unsteady flowfield around an oscillating cascade of flat plates is studied using a time marching Euler code. Exact solutions based on linear theory serve as model problems to study pressure wave propagation in the numerical solution. The importance of using proper unsteady boundary conditions, grid resolution, and time step is demonstrated. Results show that an approximate non-reflecting boundary condition based on linear theory does a good job of minimizing reflections from the inflow and outflow boundaries and allows the placement of the boundaries to be closer than cases using reflective boundary conditions. Stretching the boundary to dampen the unsteady waves is another way to minimize reflections. Grid clustering near the plates does a better job of capturing the unsteady flowfield than cases using uniform grids as long as the CFL number is less than one for a sufficient portion of the grid. Results for various stagger angles and oscillation frequencies show good agreement with linear theory as long as the grid is properly resolved.

Huff, Dennis L.↗
NASA NTRS1993

Navier-Stokes analysis of airfoils with leading edge ice accretions

A numerical analysis of the flowfield characteristics and the performance degradation of an airfoil with leading edge ice accretions was performed. The important fluid dynamic processes were identified and calculated. Among these were the leading edge separation bubble at low angles of attack, complete separation on the low pressure surface resulting in premature shell, drag rise due to the ice shape, and the effects of angle of attack on the separated flow field. Comparisons to experimental results were conducted to confirm these calculations. A computer code which solves the Navier-Stokes equations in two dimensions, ARC2D, was used to perform the calculations. A Modified Mixing Length turbulence model was developed to produce grids for several ice shape and airfoil combinations. Results indicate that the ability to predict overall performance characteristics, such as lift and drag, at low angles of attack is excellent. Transition location is important for accurately determining separation bubble shape. Details of the flowfield in and downstream of the separated regions requires some modifications. Calculations for the stalled airfoil indicate periodic shedding of vorticity that was generated aft of the ice accretion. Time averaged pressure values produce results which compare favorably with experimental information. A turbulence model which accounts for the history effects in the flow may be justified.

Potapczuk, Mark G.↗
NASA NTRS1993

Efficient simulation of incompressible viscous flow over multi-element airfoils

The incompressible, viscous, turbulent flow over single and multi-element airfoils is numerically simulated in an efficient manner by solving the incompressible Navier-Stokes equations. The solution algorithm employs the method of pseudo compressibility and utilizes an upwind differencing scheme for the convective fluxes, and an implicit line-relaxation scheme. The motivation for this work includes interest in studying high-lift take-off and landing configurations of various aircraft. In particular, accurate computation of lift and drag at various angles of attack up to stall is desired. Two different turbulence models are tested in computing the flow over an NACA 4412 airfoil; an accurate prediction of stall is obtained. The approach used for multi-element airfoils involves the use of multiple zones of structured grids fitted to each element. Two different approaches are compared; a patched system of grids, and an overlaid Chimera system of grids. Computational results are presented for two-element, three-element, and four-element airfoil configurations. Excellent agreement with experimental surface pressure coefficients is seen. The code converges in less than 200 iterations, requiring on the order of one minute of CPU time on a CRAY YMP per element in the airfoil configuration.

Rogers, Stuart E.↗