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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

NASA NTRS1986

A comparison of computational methods for three-dimensional, turbulent turbomachinery flow fields

A space-marching method and a time-marching method have been used to compute the three-dimensional turbulent flow in an end wall cascade of airfoils. Using an identical grid and turbulence model, the two codes were used to predict a variety of flow quantities. Predictions by the two methods are compared to each other and to experimental data. In general both methods predict measured quantities well, with a small edge in prediction accuracy going to the space-marching method. Secondary flow comparisons show the time-marching solution more accurately predicting the underturning of the flow in the outer portion of the end wall boundary layer while the space-marching method more accurately predicted the overturning of the flow very near the end wall. The prediction comparisons are discussed along with computational details and other attributes of the two methods.

Kirtley, K. R.↗

NASA NTRS1990

Implicit solution of three-dimensional internal turbulent flows

The scalar form of the approximate factorization method was used to develop a new code for the solution of three-dimensional internal laminar and turbulent compressible flows. The Navier-Stokes equations in their Reynolds-averaged form are iterated in time until a steady solution is reached. Evidence is given to the implicit and explicit artificial damping schemes that proved to be particularly efficient in speeding up convergence and enhancing the algorithm robustness. A conservative treatment of these terms at domain boundaries is proposed in order to avoid undesired mass and/or momentum artificial fluxes. Turbulence effects are accounted for by the zero-equation Baldwin-Lomax turbulence model and the q-omega two-equation model. For the first, an investigation on the model behavior in case of multiple boundaries is performed. The flow in a developing S-duct is then solved in the laminar regime at Reynolds number (Re) 790 and in the turbulent regime at Re=40,000 using the Baldwin-Lomax model . The Stanitz elbow is then solved using an inviscid version of the same code at M(sub inlet)=0.4. Grid dependence and convergence rate are investigated showing that for this solver the implicit damping scheme may play a critical role for convergence characteristics. The same flow at Re=2.5x10(exp 6) is solved with the Baldwin-Lomax and the q-omega models. Both approaches showed satisfactory agreement with experiments, although the q-omega model is slightly more accurate.

Michelassi, V.↗

NASA NTRS1990

Prediction of the Thrust Performance and the Flowfield of Liquid Rocket Engines

In an effort to improve the current solutions in the design and analysis of liquid propulsive engines, a computational fluid dynamics (CFD) model capable of calculating the reacting flows from the combustion chamber, through the nozzle to the external plume, was developed. The Space Shuttle Main Engine (SSME) fired at sea level, was investigated as a sample case. The CFD model, FDNS, is a pressure based, non-staggered grid, viscous/inviscid, ideal gas/real gas, reactive code. An adaptive upwinding differencing scheme is employed for the spatial discretization. The upwind scheme is based on fourth order central differencing with fourth order damping for smooth regions, and second order central differencing with second order damping for shock capturing. It is equipped with a CHMQGM equilibrium chemistry algorithm and a PARASOL finite rate chemistry algorithm using the point implicit method. The computed flow results and performance compared well with those of other standard codes and engine hot fire test data. In addition, the transient nozzle flowfield calculation was also performed to demonstrate the ability of FDNS in capturing the flow separation during the startup process.

Wang, T.-S.↗

NASA NTRS1992

Sonic boom prediction for the Langley Mach 2 low-boom configuration

Sonic boom pressure signatures and aerodynamic force data for the Langley Mach 2 low sonic boom configuration were computed using the TranAir full-potential code. A solution-adaptive Cartesian grid scheme is utilized to compute off-body flow field data. Computations were performed with and without nacelles at several angles of attack. Force and moment data were computed to measure nacelle effects on the aerodynamic characteristics and sonic boom footprints of the model. Pressure signatures were computed both on and off ground-track. Near-field pressure signature computations on ground-track were in good agreement with experimental data. Computed off ground-track signatures showed that maximum pressure peaks were located off ground-track and were significantly higher than the signatures on ground-track. Bow shocks from the nacelle inlets increased lift and drag, and also increased the magnitude of the maximum pressure both on and off ground-track.

Madson, Michael D.↗

NASA NTRS1993

Averaging techniques for steady and unsteady calculations of a transonic fan stage

It is often desirable to characterize a turbomachinery flow field with a few lumped parameters such as total pressure ratio or stage efficiency. Various averaging schemes may be used to compute these parameters. The momentum, energy, and area averaging schemes are described and compared. The schemes were compared for two computed solutions of the midspan section of a transonic fan stage: a steady averaging-plane solution in which average rotor outflow conditions were used as stator inflow conditions, and an unsteady rotor-stator interaction solution. The solutions were computed on identical grids using similar Navier-Stokes codes and an algebraic turbulence model. The unsteady solution is described, some unsteady flow phenomena are discussed, and the steady pressure distributions are compared. Despite large unsteady pressure fluctuations on the stator surface, the steady pressure distribution matched the average unsteady distribution almost exactly. Stator wake profiles, stator loss coefficient, and stage efficiency were computed for the two solutions with the three averaging schemes and are compared. In general, the energy averaging scheme gave good agreement between the averaging-plane solution and the time-averaged unsteady solution, even though certain phenomena due to unsteady wake migration were neglected.

Wyss, M. L.↗

NASA NTRS1995

Automatic structured grid generation using Gridgen (some restrictions apply)

The authors have noticed in the recent grid generation literature an emphasis on the automation of structured grid generation. The motivation behind such work is clear; grid generation is easily the most despised task in the grid-analyze-visualize triad of computational analysis (CA). However, because grid generation is closely coupled to both the design and analysis software and because quantitative measures of grid quality are lacking, 'push button' grid generation usually results in a compromise between speed, control, and quality. Overt emphasis on automation obscures the substantive issues of providing users with flexible tools for generating and modifying high quality grids in a design environment. In support of this paper's tongue-in-cheek title, many features of the Gridgen software are described. Gridgen is by no stretch of the imagination an automatic grid generator. Despite this fact, the code does utilize many automation techniques that permit interesting regenerative features.

Chawner, John R.↗

NASA NTRS1998

Euler Technology Assessment for Preliminary Aircraft Design-Unstructured/Structured Grid NASTD Application for Aerodynamic Analysis of an Advanced Fighter/Tailless Configuration

This study supports the NASA Langley sponsored project aimed at determining the viability of using Euler technology for preliminary design use. The primary objective of this study was to assess the accuracy and efficiency of the Boeing, St. Louis unstructured grid flow field analysis system, consisting of the MACGS grid generation and NASTD flow solver codes. Euler solutions about the Aero Configuration/Weapons Fighter Technology (ACWFT) 1204 aircraft configuration were generated. Several variations of the geometry were investigated including a standard wing, cambered wing, deflected elevon, and deflected body flap. A wide range of flow conditions, most of which were in the non-linear regimes of the flight envelope, including variations in speed (subsonic, transonic, supersonic), angles of attack, and sideslip were investigated. Several flowfield non-linearities were present in these solutions including shock waves, vortical flows and the resulting interactions. The accuracy of this method was evaluated by comparing solutions with test data and Navier-Stokes solutions. The ability to accurately predict lateral-directional characteristics and control effectiveness was investigated by computing solutions with sideslip, and with deflected control surfaces. Problem set up times and computational resource requirements were documented and used to evaluate the efficiency of this approach for use in the fast paced preliminary design environment.

Michal, Todd R.↗

NASA NTRS1999

Overflow: Facts on Friction

This paper presents results of three minor studies into the behavior of the OVERFLOW with respect to the prediction of skin friction drag on wing bodies at cruise Mach number and wind tunnel Reynolds number. The studies include a preliminary assessment of the behavior of the two new 2-equation turbulence models introduced with the latest version of OVERFLOW (v. 1.8f), an investigation into potential improvements in the matrix dissipation scheme currently implemented in OVERFLOW, and an analysis of the observed sensitivity of the code's skin friction predictions to grid stretching at solid surface boundaries.

Lawrence, Scott↗

NASA NTRS1999

Installed Transonic 2D Nozzle Nacelle Boattail Drag Study

The Transonic Nozzle Boattail Drag Study was initiated in 1995 to develop an understanding of how external nozzle transonic aerodynamics effect airplane performance and how strongly those effects are dependent on nozzle configuration (2D vs. axisymmetric). MDC analyzed the axisymmetric nozzle. Boeing subcontracted Northrop-Grumman to analyze the 2D nozzle. AU participants analyzed the AGARD nozzle as a check-out and validation case. Once the codes were checked out and the gridding resolution necessary for modeling the separated flow in this region determined, the analysis moved to the installed wing/body/nacelle/diverter cases. The boat tail drag validation case was the AGARD B.4 rectangular nozzle. This test case offered both test data and previous CFD analyses for comparison. Results were obtained for test cases B.4.1 (M=0.6) and B.4.2 (M=0.938) and compared very well with the experimental data. Once the validation was complete a CFD grid was constructed for the full Ref. H configuration (wing/body/nacelle/diverter) using a combination of patched and overlapped (Chimera) grids. This was done to ensure that the grid topologies and density would be adequate for the full model. The use of overlapped grids allowed the same grids from the full configuration model to be used for the wing/body alone cases, thus eliminating the risk of grid differences affecting the determination of the installation effects. Once the full configuration model was run and deemed to be suitable the nacelle/diverter grids were removed and the wing/body analysis performed. Reference H wing/body results were completed for M=0.9 (a=0.0, 2.0, 4.0, 6.0 and 8.0), M=1.1 (a=4.0 and 6.0) and M=2.4 (a=0.0, 2.0, 4.4, 6.0 and 8.0). Comparisons of the M=0.9 and M=2.4 cases were made with available wind tunnel data and overall comparisons were good. The axi-inlet/2D nozzle nacelle was analyzed isolated. The isolated nacelle data coupled with the wing/body result enabled the interference effects of the installed nacelles to be determined. Isolated nacelle mm were made at M=0.9 and M=1.1 for both the supersonic and transonic nozzle settings. AU of the isolated nacelle cases were run at alpha=0. Full configuration runs were to be made at Mach numbers of 0.9, 1.1, and 2.4 (the same as the wing/body and isolated nacelles). Both the isolated nacelles and installed nacelles were run with inlet conditions designed to give zero spillage. This was to be done in order to isolate the boattail effects as much as possible. Full configuration runs with the supersonic nozzles were completed for M=0.9 and 1.1 at a=4.0 and 6.0 (4 runs total) and with the transonic nozzles at M=0.9 and 1.1 at a=2.0, 4.0 and 6.0 (6 runs total). Drag breakdowns were completed for the M=0.9 and M= 1.1 showing favorable interference drag for both cases.

Malone, Michael B.↗

NASA NTRS2010

Micro-Ramps for External Compression Low-Boom Inlets

The application of vortex generators for flow control in an external compression, axisymmetric, low-boom concept inlet was investigated using RANS simulations with three-dimensional (3-D), structured, chimera (overset) grids and the WIND-US code. The low-boom inlet design is based on previous scale model 1- by 1-ft wind tunnel tests and features a zero-angle cowl and relaxed isentropic compression centerbody spike, resulting in defocused oblique shocks and a weak terminating normal shock. Validation of the methodology was first performed for micro-ramps in supersonic flow on a flat plate with and without oblique shocks. For the inlet configuration, simulations with several types of vortex generators were conducted for positions both upstream and downstream of the terminating normal shock. The performance parameters included incompressible axisymmetric shape factor, separation area, inlet pressure recovery, and massflow ratio. The design of experiments (DOE) methodology was used to select device size and location, analyze the resulting data, and determine the optimal choice of device geometry. The optimum upstream configuration was found to substantially reduce the post-shock separation area but did not significantly impact recovery at the aerodynamic interface plane (AIP). Downstream device placement allowed for fuller boundary layer velocity profiles and reduced distortion. This resulted in an improved pressure recovery and massflow ratio at the AIP compared to the baseline solid-wall configuration.

Rybalko, Michael↗

NASA NTRS2013

A Simple GPU-Accelerated Two-Dimensional MUSCL-Hancock Solver for Ideal Magnetohydrodynamics

We describe our experience using NVIDIA's CUDA (Compute Unified Device Architecture) C programming environment to implement a two-dimensional second-order MUSCL-Hancock ideal magnetohydrodynamics (MHD) solver on a GTX 480 Graphics Processing Unit (GPU). Taking a simple approach in which the MHD variables are stored exclusively in the global memory of the GTX 480 and accessed in a cache-friendly manner (without further optimizing memory access by, for example, staging data in the GPU's faster shared memory), we achieved a maximum speed-up of approx. = 126 for a sq 1024 grid relative to the sequential C code running on a single Intel Nehalem (2.8 GHz) core. This speedup is consistent with simple estimates based on the known floating point performance, memory throughput and parallel processing capacity of the GTX 480.

graphics processing units↗

NASA NTRS1986

High speed viscous flow calculations about complex configurations

Applications of the NASA Ames Parabolized Navier-Stokes (PNS) code to a variety of complex generic configurations is presented. The algorithm, boundary conditions, initial conditions, and grid generators are discussed as applied to these configurations. The PNS code was used as the mainline procedure to numerically simulate the viscous supersonic flow over these generic configurations. The turbulence model that was used in this study is the Baldwin-Lomax model. The boundary conditions are the usual viscous no slip at the wall, and a characteristic procedure is used to fit the bow shock wave which is the outermost boundary. An elliptic grid generator is employed to discretize the flow domain. In addition, an equilibrium air capability has been incorporated into the code. It uses the curve fits of Tannehill, et al. The flow regimes vary from a Mach number of 2 up to 25. Both laminar and turbulent flow are considered. Varying angles of attack have also been computed. Configurations vary from simple cone-type bodies to lifting winged bodies, such as the space shuttle or the generic supersonic cruise fighter.

Chaussee, D. S.↗

DOE OSTI2018

Investigating the Unusual Spectroscopic Time Evolution in SN 2012fr

The type Ia supernova (SN) 2012fr displayed a curious combination of its Si ii λλ5972, 6355 features. This includes the ratio of their pseudo-equivalent widths, placing it at the border of the shallow silicon (SS) and core normal (CN) spectral subtype in the Branch diagram, while the Si ii λ6355 expansion velocities place it as a high-velocity (HV) object in the Wang et al. spectral type that most intriguingly evolves slowly, placing it in the low-velocity gradient (LVG) typing of Benetti et al. Only 5% of SNe Ia are HV and located in the SS+CN portion of the Branch diagram, and fewer than 10% of SNe Ia are both HV and LVG. These features point toward SN 2012fr being quite unusual, similar in many ways to the peculiar SN 2000cx. We modeled the spectral evolution of SN 2012fr to see if we could gain some insight into its evolutionary behavior. We use the parameterized radiative transfer code SYNOW to probe the abundance stratification of SN 2012fr at pre-maximum, maximum, and post-maximum light epochs. We also use a grid of W7 models in the radiative transfer code PHOENIX to probe the effect of different density structures on the formation of the Si ii λ6355 absorption feature at post-maximum epochs. We find that the unusual features observed in SN 2012fr are likely due to a shell-like density enhancement in the outer ejecta. We report on possible reasons for atypical Ca ii absorption features, and suggest that they are related to the Si ii features.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1987

Phase space simulation of collisionless stellar systems on the massively parallel processor

A numerical technique for solving the collisionless Boltzmann equation describing the time evolution of a self gravitating fluid in phase space was implemented on the Massively Parallel Processor (MPP). The code performs calculations for a two dimensional phase space grid (with one space and one velocity dimension). Some results from calculations are presented. The execution speed of the code is comparable to the speed of a single processor of a Cray-XMP. Advantages and disadvantages of the MPP architecture for this type of problem are discussed. The nearest neighbor connectivity of the MPP array does not pose a significant obstacle. Future MPP-like machines should have much more local memory and easier access to staging memory and disks in order to be effective for this type of problem.

White, Richard L.↗

NASA NTRS1994

A 3D-CFD code for accurate prediction of fluid flows and fluid forces in seals

Current and future turbomachinery requires advanced seal configurations to control leakage, inhibit mixing of incompatible fluids and to control the rotodynamic response. In recognition of a deficiency in the existing predictive methodology for seals, a seven year effort was established in 1990 by NASA's Office of Aeronautics Exploration and Technology, under the Earth-to-Orbit Propulsion program, to develop validated Computational Fluid Dynamics (CFD) concepts, codes and analyses for seals. The effort will provide NASA and the U.S. Aerospace Industry with advanced CFD scientific codes and industrial codes for analyzing and designing turbomachinery seals. An advanced 3D CFD cylindrical seal code has been developed, incorporating state-of-the-art computational methodology for flow analysis in straight, tapered and stepped seals. Relevant computational features of the code include: stationary/rotating coordinates, cylindrical and general Body Fitted Coordinates (BFC) systems, high order differencing schemes, colocated variable arrangement, advanced turbulence models, incompressible/compressible flows, and moving grids. This paper presents the current status of code development, code demonstration for predicting rotordynamic coefficients, numerical parametric study of entrance loss coefficients for generic annular seals, and plans for code extensions to labyrinth, damping, and other seal configurations.

Athavale, M. M.↗

NASA NTRSDate not supplied

Rapid Hypersonic Simulations using US3D and Pointwise

For hypersonic simulations, unstructured flow solvers typically have problems predicting surface heat fluxes when strong shocks are present. To address these issues, this paper outlines a workflow that applies best practices developed for structured grids to unstructured meshes. In addition, unstructured grid generation can significantly reduce the time required to create quality grids for complex geometries. Several examples are computed using DPLR, a structured grid flow solver, and US3D flow, an unstructured mesh solver. Results from the two codes are compared, and they show excellent agreement. Overall, the unstructured grid workflow offers a viable and attractive alternative for hypersonic simulations.

C. Tang↗

NASA NTRSDate not supplied

Rapid Hypersonic Simulations using US3D and Pointwise

Chun Tang↗

NASA NTRS1998

Numerical Study of Interaction of a Vortical Density Inhomogeneity with Shock and Expansion Waves

We studied the interaction of a vortical density inhomogeneity (VDI) with shock and expansion waves. We call the VDI the region of concentrated vorticity (vortex) with a density different from that of ambiance. Non-parallel directions of the density gradient normal to the VDI surface and the pressure gradient across a shock wave results in an additional vorticity. The roll-up of the initial round VDI towards a non-symmetrical shape is studied numerically. Numerical modeling of this interaction is performed by a 2-D Euler code. The use of an adaptive unstructured numerical grid makes it possible to obtain high accuracy and capture regions of induced vorticity with a moderate overall number of mesh points. For the validation of the code, the computational results are compared with available experimental results and good agreement is obtained. The interaction of the VDI with a propagating shock wave is studied for a range of initial and induced circulations and obtained flow patterns are presented. The splitting of the VDI develops into the formation of a non-symmetrical vortex pair and not in a set of vortices. A method for the analytical computation of an overall induced circulation Gamma(sub 1) as a result of the interaction of a moving VDI with a number of waves is proposed. Simplified, approximated, expressions for Gamma(sub 1) are derived and their accuracy is discussed. The splitting of the VDI passing through the Prandtl-Meyer expansion wave is studied numerically. The obtained VDI patterns are compared to those for the interaction of the VDI with a propagating shock wave for the same values of initial and induced circulations. These patterns have similar shapes for corresponding time moments.

Povitsky, A.↗