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,513 records · Page 84

Towards 100,000 CPU Cycle-Scavenging by Genetic Algorithms

We examine a web-centric design using standard tools such as web servers, web browsers, PHP, and mySQL. We also consider the applicability of Information Power Grid tools such as the Globus (no relation to the author) Toolkit. We intend to implement this architecture with JavaGenes running on at least two cycle-scavengers: Condor and United Devices. JavaGenes, a genetic algorithm code written in Java, will be used to evolve multi-species reactive molecular force field parameters.

Globus, Al↗

A Two-Stage Procedure Toward the Efficient Implementation of PANS and Other Hybrid Turbulence Models

The main objective of this article is to introduce and to show the implementation of a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for Partial Averaged Navier-Stokes (PANS) and other hybrid models. It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The first step is to solve the unsteady or steady Reynolds Averaged Navier-Stokes (URANS/RANS) equations. From this preprocessing step, the turbulence length-scale field is obtained. This is then used to compute the characteristic length-scale ratio between the turbulence scale and the grid spacing. Based on this ratio, we can assess the finest scale resolution that a given grid for a given flow can support. Along with other additional criteria, we are able to analytically identify the appropriate hybrid solver resolution for different regions of the flow. This procedure removes the grid dependency issue that affects the results produced by different hybrid procedures in solving unsteady flows. The formulation, implementation methodology, and validation example are presented. We implemented this capability in a production Computational Fluid Dynamics (CFD) code, PAB3D, for the simulation of unsteady flows.

Abdol-Hamid, Khaled S.↗

USM3D Simulations for Second Sonic Boom Workshop

The NASA Tetrahedral Unstructured Software System with the USM3D flow solver was used to compute test cases for the Second AIAA Sonic Boom Prediction Workshop. The intent of this report is to document the USM3D results for SBPW2 test cases. The test cases included an axisymmetric equivalent area body, a JAXA wing body, a NASA low boom supersonic configuration modeled with flow through nacelles and engine boundary conditions. All simulations were conducted for a free stream Mach number of 1.6, zero degrees angle of attack, and a Reynolds number of 5.7 million per meter. Simulations were conducted on tetrahedral grids provided by the workshop committee, as well as a family of grids generated by an in-house approach for sonic boom analyses known as BoomGrid using current best practices. The near-field pressure signatures were extracted and propagated to the ground with the atmospheric propagation code, sBOOM. The USM3D near-field pressure signatures, corresponding sBOOM ground signatures, and loudness levels on the ground are compared with mean values from other workshop participants.

Elmiligui, Alaa↗

USM3D Simulations for Second Sonic Boom Workshop

The NASA Tetrahedral Unstructured Software System with the USM3D flow solver was used to compute the test cases for the Second AIAA Sonic Boom Prediction Workshop. The intent of this article is to document the USM3D results for SBPW2 test cases. The test cases included an axisymmetric equivalent area body, a JAXA wing body, a NASA low boom supersonic configuration modeled with flow-through nacelles and engine boundary conditions. All simulations were conducted for a free stream Mach number of 1.6, zero degrees angle of attack, and a Reynolds number of 5.7 million per meter. Simulations were conducted on the tetrahedral grids provided by the workshop committee, as well as a family of grids generated by an in-house approach for sonic boom analyses known as BoomGrid using current best practices. The near-field pressure signatures were extracted and propagated to the ground with the atmospheric propagation code, sBOOM. The USM3D near-field pressure signatures, corresponding sBOOM ground signatures, and loudness levels on the ground are compared with mean values from other workshop participants.

Elmiligui, Alaa↗

Review of WEC-Sim development and applications

WEC-Sim (Wave Energy Converter Simulator) is an open-source code for simulating wave energy converters, which has been actively developed and applied to simulate a wide variety of device archetypes, and has become a popular tool since its release. This paper reviewed the development efforts and the usage of WEC-Sim. The publications considered in this study have been broken down into six topic areas, namely feature development, experimental validation, device modeling, control modeling, PTO and grid modeling, and novel applications, which even includes some non-wave energy applications. This review paper has also attempted to recognize the contributions of the broader WEC-Sim development effort, meaning not only the internal WEC-Sim development team but also the external efforts from the academia researchers and technology developers around the world. The growing trend of external applications of WEC-Sim has demonstrated the broader acceptance of the open-source code, and how WEC-Sim has been used in a certain topic area also highlights the potential future development needs.

16 TIDAL AND WAVE POWER↗

Development and application of computational aerothermodynamics flowfield computer codes

Presented is a collection of papers on research activities carried out during the funding period of October 1991 to March 1992. Topics covered include: blunt body flows in thermochemical equilibrium; thermochemical relaxation in high enthalpy nozzle flow; single expansion ramp nozzle simulations; lunar return aerobraking; line boundary problem for three dimensional grids; and unsteady shock induced combustion.

Venkatapathy, Ethiraj↗

Model atmospheres for M (sub)dwarf stars. 1: The base model grid

We have calculated a grid of more than 700 model atmospheres valid for a wide range of parameters encompassing the coolest known M dwarfs, M subdwarfs, and brown dwarf candidates: 1500 less than or equal to T(sub eff) less than or equal to 4000 K, 3.5 less than or equal to log g less than or equal to 5.5, and -4.0 less than or equal to (M/H) less than or equal to +0.5. Our equation of state includes 105 molecules and up to 27 ionization stages of 39 elements. In the calculations of the base grid of model atmospheres presented here, we include over 300 molecular bands of four molecules (TiO, VO, CaH, FeH) in the JOLA approximation, the water opacity of Ludwig (1971), collision-induced opacities, b-f and f-f atomic processes, as well as about 2 million spectral lines selected from a list with more than 42 million atomic and 24 million molecular (H2, CH, NH, OH, MgH, SiH, C2, CN, CO, SiO) lines. High-resolution synthetic spectra are obtained using an opacity sampling method. The model atmospheres and spectra are calculated with the generalized stellar atmosphere code PHOENIX, assuming LTE, plane-parallel geometry, energy (radiative plus convective) conservation, and hydrostatic equilibrium. The model spectra give close agreement with observations of M dwarfs across a wide spectral range from the blue to the near-IR, with one notable exception: the fit to the water bands. We discuss several practical applications of our model grid, e.g., broadband colors derived from the synthetic spectra. In light of current efforts to identify genuine brown dwarfs, we also show how low-resolution spectra of cool dwarfs vary with surface gravity, and how the high-regulation line profile of the Li I resonance doublet depends on the Li abundance.

Allard, France↗

Potential Research and Development Opportunities for Light Water Reactor Spent Nuclear Fuel at INL

This report documents the commercial spent nuclear fuel currently in storage at Idaho National Laboratory with a specific focus on the potential research and development (R&D) opportunities for this material. Commercial spent nuclear fuel is stored in dry casks on the INL site at the Idaho Nuclear Technology and Engineering Center (INTEC, formerly the Idaho Chemical Processing Plant or ICPP). This spent fuel was part of the DOE-industry cooperative dry cask demonstration project that began in 1983 and from historical fuel testing and treatment activities that had been performed at the Test Area North (TAN)-607 facility. Most of the fuel is “low burnup,” but given that it was removed from the respective reactors between 1970 and 1985, it was representative for its time. The majority of the fuel is of the 15×15 Westinghouse PWR variety, though additional fuel types (e.g., TN-BRP) and configurations are also stored. Potential R&D opportunities identified are (1) the REA-2023 cask, (2) the TN-BRP cask, and (3) the Castor V/21 cask. Fuel from the TN-24P may be an alternative to the intact assemblies in the Castor V/21. Specifically, the REA-2023 cask contains 9 rods that have been exposed to air since 2005. These rods have a detailed pre-loading characterization and thus any changes in corrosion or geometry can most likely be attributed to the storage environment. In addition, similar rods to the nine in REA-2023 reside in the Castor V/21 cask whose environment has been maintained with inert nitrogen. The Castor V/21 cask also contains intact PWR assemblies that could be used to perform limited confirmation of the pyroprocessing approach. An alternative to the use of the Castor assemblies for pyroprocessing is the consolidated fuel that is stored in the TN-24P, which would not require removal from the assembly grid since it has already been packaged in accessible two-piece shell strongback containers. The TN-BRP cask contains a number of Gd-doped, low burnup assemblies that could be of use to validate code packages used to justify burnup-up credit for boiling water reactors (BWRs). The INL dry casks are stored on the CPP-2707 pad at INTEC, adjacent to the CPP-603 Irradiated Fuel Storage Facility, which incorporates a dry hot cell that has the potential to receive a full scale commercial cask of the types identified. Several evaluations have been made regarding this task, including full scale mockup testing of the physical clearances required to remove and handle fuel assemblies and components. A proposed schedule of 14 months and cost of $5.3 million was the outcome of analysis performed by Fluor Idaho and is included as Appendix A & B. In addition to the spent fuel at INL, there is also a full-scale mockup assembly that was fabricated for the Protypical Consolidation Demonstration Project to validate fuel rod consolidation approaches. The assembly was fabricated by Westinghouse as a 15 x 15 PWR, using Zircaloy-4 tubing and steel end fittings in accordance with product specifications. The only departure from specification is the use of copper metal in lieu of uranium oxide pellets to achieve a weight that is within 10% of the standard.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Hanford Site Composite Analysis: LLBG-200-W B Vadose Zone Model

The objectives of the vadose modeling for the updated Hanford Site Composite Analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates to the CA saturated zone model (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort: carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP 1 ) simulator, officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP), is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the saturated zone model. The rates will be summed over the 100 by 100 m saturated zone model grid cells that fall within the vadose zone model source domain.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Magnetic Reconnection

Magnetic reconnection is a fundamental plasma physics process ubiquitous in astrophysics, and important in both magnetic confinement fusion and space weather. The MARZ fundamental science program was recently established on Z to enable the first laboratory astrophysics platform able to access and study the strongly radiatively cooled magnetic reconnection regime. Simulations of this system have successfully used a resistive-MHD approach, but in some regions of parameter space Hall physics has the potential to be important. We describe implementation of a Hall method on a staggered grid resistive-MHD method (compatible with the approach used to model MARZ experiments. We then present a different Hall method based on cell-centered field quantities. Both approaches have been implemented in the Sandia KRAKEN code, to enable us to contrast different numerical Hall-MHD methods within the same HED code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The ICON-A model for direct QBO simulations on GPUs (version icon-cscs:baf28a514)

Abstract. Classical numerical models for the global atmosphere, as used for numerical weather forecasting or climate research, have been developed for conventional central processing unit (CPU) architectures. This hinders the employment of such models on current top-performing supercomputers, which achieve their computing power with hybrid architectures, mostly using graphics processing units (GPUs). Thus also scientific applications of such models are restricted to the lesser computer power of CPUs. Here we present the development of a GPU-enabled version of the ICON atmosphere model (ICON-A), motivated by a research project on the quasi-biennial oscillation (QBO), a global-scale wind oscillation in the equatorial stratosphere that depends on a broad spectrum of atmospheric waves, which originates from tropical deep convection. Resolving the relevant scales, from a few kilometers to the size of the globe, is a formidable computational problem, which can only be realized now on top-performing supercomputers. This motivated porting ICON-A, in the specific configuration needed for the research project, in a first step to the GPU architecture of the Piz Daint computer at the Swiss National Supercomputing Centre and in a second step to the JUWELS Booster computer at the Forschungszentrum Jülich. On Piz Daint, the ported code achieves a single-node GPU vs. CPU speedup factor of 6.4 and allows for global experiments at a horizontal resolution of 5 km on 1024 computing nodes with 1 GPU per node with a turnover of 48 simulated days per day. On JUWELS Booster, the more modern hardware in combination with an upgraded code base allows for simulations at the same resolution on 128 computing nodes with 4 GPUs per node and a turnover of 133 simulated days per day. Additionally, the code still remains functional on CPUs, as is demonstrated by additional experiments on the Levante compute system at the German Climate Computing Center. While the application shows good weak scaling over the tested 16-fold increase in grid size and node count, making also higher resolved global simulations possible, the strong scaling on GPUs is relatively poor, which limits the options to increase turnover with more nodes. Initial experiments demonstrate that the ICON-A model can simulate downward-propagating QBO jets, which are driven by wave–mean flow interaction.

54 ENVIRONMENTAL SCIENCES↗

Simplified clustering of nonorthogonal grids generated by elliptic partial differential equations

A simple clustering transformation is combined with the Thompson, Thames, and Mastin (TTM) method of generating computational grids to produce controlled mesh spacings. For various practical grids, the resulting hybrid scheme is easier to apply than the inhomogeneous clustering terms included in the TTM method for this purpose. The technique is illustrated in application to airfoil problems, and listings of a FORTRAN computer code for this usage are included.

Sorenson, R. L.↗

Self-adjusting grid methods for one-dimensional hyperbolic conservation laws

The automatic adjustment of a grid which follows the dynamics of the numerical solution of hyperbolic conservation laws is given. The grid motion is determined by averaging the local characteristic velocities of the equations with respect to the amplitudes of the signals. The resulting algorithm is a simple extension of many currently popular Godunov-type methods. Computer codes using one of these methods can be easily modified to add the moving mesh as an option. Numerical examples are given that illustrate the improved accuracy of Godunov's and Roe's methods on a self-adjusting mesh. Previously announced in STAR as N83-15008

Harten, A.↗

A locally implicit method for fluid flow problems

The fluid flow inside the space shuttle main engine (SSME) traverses through a complex geometrical configuration. The flow is compressible, viscous, and turbulent with pockets of separated regions. Several computer codes are being developed to solve three dimensional Navier-Stokes equations with different turbulence models for analyzing the SSME internal flow. The locally implicit scheme is a computationally efficient scheme which converges rapidly in multi-grid modes for elliptic problems. It has the promise of providing a rapidly converging algorithm for steady-state viscous flow problems.

Reddy, K. C.↗

Three-dimensional multigrid Navier-Stokes computations for turbomachinery applications

The fully three-dimensional, time-dependent compressible Navier-Stokes equations in cylindrical coordinates are presently used, in conjunction with the multistage Runge-Kutta numerical integration scheme for solution of the governing flow equations, to simulate complex flowfields within turbomechanical components whose pertinent effects encompass those of viscosity, compressibility, blade rotation, and tip clearance. Computed results are presented for selected cascades, emphasizing the code's capabilities in the accurate prediction of such features as airfoil loadings, exit flow angles, shocks, and secondary flows. Computations for several test cases have been performed on a Cray-YMP, using nearly 90,000 grid points.

Subramanian, S. V.↗

Experimental and computational studies of hovering rotor flows

HELIX-I, an essentially standard full-potential CFD helicopter rotor code, is unique in its use of the Clebisch kinematical flow description to specify a freely convecting wake and its capacity for predicting hover performance. A study is presently performed to assess the method's sensitivity to grids and solution-starting techniques. The effects of these parameters on thrust, power, load distribution, and wake geometry are ascertained and compared with an extensive rotor data base. The use of a fairly accurate starting solution yields no obvious advantage over the use of a novel starting method which employs a succession of diminishing artificial flows.

Nsi Mba, M.↗

Computational study of generic hypersonic vehicle flow fields

The geometric data of the generic hypersonic vehicle configuration included body definitions and preliminary grids for the forebody (nose cone excluded), midsection (propulsion system excluded), and afterbody sections. This data was to be augmented by the nose section geometry (blunt conical section mated with the noncircular cross section of the forebody initial plane) along with a grid and a detailed supersonic combustion ramjet (scramjet) geometry (inlet and combustor) which should be merged with the nozzle portion of the afterbody geometry. The solutions were to be obtained by using a Navier-Stokes (NS) code such as TUFF for the nose portion, a parabolized Navier-Stokes (PNS) solver such as the UPS and STUFF codes for the forebody, a NS solver with finite rate hydrogen-air chemistry capability such as TUFF and SPARK for the scramjet and a suitable solver (NS or PNS) for the afterbody and external nozzle flows. The numerical simulation of the hypersonic propulsion system for the generic hypersonic vehicle is the major focus of this entire work. Supersonic combustion ramjet is such a propulsion system, hence the main thrust of the present task has been to establish a solution procedure for the scramjet flow. The scramjet flow is compressible, turbulent, and reacting. The fuel used is hydrogen and the combustion process proceeds at a finite rate. As a result, the solution procedure must be capable of addressing such flows.

Narayan, Johnny R.↗

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 1: Code validation and stationary accretors

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing stationary spheres of varying sizes (from 10.0 down to 0.02 Bondi radii) accrete matter from a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3 or 1.2). To accommodate the long-range gravitational forces, the extent of the computational volume is typically a factor of 100 larger than the radius of the accretor. We compare the numerical mass accretion rates with the theoretical predictions of Bondi, to assess the validity of the code. The hydrodynamics is modeled by the piecewise parabolic method. No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (6-8) grids around the stationary sphere, each finer grid being a factor of 2 smaller spatially than the next coarser grid. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid while at the same time evolving the gas on the coarser grids. The accretion rates derived numerically are in in very good agreement (to about 10% over several orders of magnitude) with the values given by Bondi for a stationary accretor within a hydrodynamic medium. However, the equations have to be changed in order to include the finite size of the accretor (in some cases very large compared to the sonic point or even to the Bondi radius).

Ruffert, Maximilian↗