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.

SourceResearch connectionFrom yearThrough year

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

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS1981

Flow aerodynamics modeling of an MHD swirl combustor - Calculations and experimental verification

The paper describes a computer code for calculating the flow dynamics of a constant-density flow in the second-stage trumpet shaped nozzle section of a two-stage MHD swirl combustor for application to a disk generator. The primitive pressure-velocity variable, finite-difference computer code has been developed for the computation of inert nonreacting turbulent swirling flows in an axisymmetric MHD model swirl combustor. The method and program involve a staggered grid system for axial and radial velocities, and a line relaxation technique for the efficient solution of the equations. The code produces as output the flow field map of the nondimensional stream function, axial and swirl velocity. It was found that the best location for seed injection to obtain a uniform distribution at the combustor exit is in the central location for seed injected at the entrance to the second stage combustor.

Gupta, A. K.↗

NASA NTRS2010

Statistical Analysis of CFD Solutions from the Fourth AIAA Drag Prediction Workshop

A graphical framework is used for statistical analysis of the results from an extensive N-version test of a collection of Reynolds-averaged Navier-Stokes computational fluid dynamics codes. The solutions were obtained by code developers and users from the U.S., Europe, Asia, and Russia using a variety of grid systems and turbulence models for the June 2009 4th Drag Prediction Workshop sponsored by the AIAA Applied Aerodynamics Technical Committee. The aerodynamic configuration for this workshop was a new subsonic transport model, the Common Research Model, designed using a modern approach for the wing and included a horizontal tail. The fourth workshop focused on the prediction of both absolute and incremental drag levels for wing-body and wing-body-horizontal tail configurations. This work continues the statistical analysis begun in the earlier workshops and compares the results from the grid convergence study of the most recent workshop with earlier workshops using the statistical framework.

Statistical analysis↗

NASA NTRS1990

Relative efficiency and accuracy of two Navier-Stokes codes for simulating attached transonic flow over wings

In the present study, two codes which solve the three-dimensional Thin-Layer Navier-Stokes (TLNS) equations are used to compute the steady-state flow for two test cases representing typical finite wings at transonic conditions. Several grids of C-O topology and varying point densities are used. After a description of each code and test case, standards for determining code efficiency and accuracy are defined and applied to determine the relative performance of the two codes in predicting turbulent transonic wing flows. Comparisons of computed surface pressure distributions with experimental data are made.

Bonhaus, Daryl L.↗

NASA NTRSDate not supplied

TPSAS-NF1676L-10767-DND

Version 2.1ae of the NASA CFD code OVERFLOW has been used to simulate the flow past tandem cylinders. The multiblock grid system consisted of 15 overset meshes with a total grid count of 15.7 million points. The grid extended 3 cylinder diameters in the spanwise direction with periodic boundary conditions at the ends. Riemann conditions were used at the outer boundaries, and a no-slip condition on the cylinder surfaces. The 3rd order HLLC upwind spatial discretization scheme was used with a dual-time stepping algorithm. A modified version of the Delayed Detached Eddy Simulation hybrid RANS/LES turbulence model proposed by Spalart was used in a fully turbulent mode. The code was run on 53 Core 2 Duo nodes using MPI across nodes and OPENMP within each node. Simulation results will be compared with the available experimental data.

David P. Lockard↗

NASA NTRS2008

Efficient Unstructured Grid Adaptation Methods for Sonic Boom Prediction

This paper examines the use of two grid adaptation methods to improve the accuracy of the near-to-mid field pressure signature prediction of supersonic aircraft computed using the USM3D unstructured grid flow solver. The first method (ADV) is an interactive adaptation process that uses grid movement rather than enrichment to more accurately resolve the expansion and compression waves. The second method (SSGRID) uses an a priori adaptation approach to stretch and shear the original unstructured grid to align the grid with the pressure waves and reduce the cell count required to achieve an accurate signature prediction at a given distance from the vehicle. Both methods initially create negative volume cells that are repaired in a module in the ADV code. While both approaches provide significant improvements in the near field signature (< 3 body lengths) relative to a baseline grid without increasing the number of grid points, only the SSGRID approach allows the details of the signature to be accurately computed at mid-field distances (3-10 body lengths) for direct use with mid-field-to-ground boom propagation codes.

Campbell, Richard L.↗

NASA NTRS2004

Agglomeration Multigrid for an Unstructured-Grid Flow Solver

An agglomeration multigrid scheme has been implemented into the sequential version of the NASA code USM3Dns, tetrahedral cell-centered finite volume Euler/Navier-Stokes flow solver. Efficiency and robustness of the multigrid-enhanced flow solver have been assessed for three configurations assuming an inviscid flow and one configuration assuming a viscous fully turbulent flow. The inviscid studies include a transonic flow over the ONERA M6 wing and a generic business jet with flow-through nacelles and a low subsonic flow over a high-lift trapezoidal wing. The viscous case includes a fully turbulent flow over the RAE 2822 rectangular wing. The multigrid solutions converged with 12%-33% of the Central Processing Unit (CPU) time required by the solutions obtained without multigrid. For all of the inviscid cases, multigrid in conjunction with an explicit time-stepping scheme performed the best with regard to the run time memory and CPU time requirements. However, for the viscous case multigrid had to be used with an implicit backward Euler time-stepping scheme that increased the run time memory requirement by 22% as compared to the run made without multigrid.

Frink, Neal↗

NASA NTRSDate not supplied

Coupling a Particle Trajectory Capability to the US3D CFD Flow Solver

The purpose of the current research effort is to incorporate particle trajectory equations into the US3D CFD flow solver. The integrated code will allow higher-fidelity modeling of particle-fluid interactions and allow for the simulation of physical effects such as surface heating augmentation due to particle impacts. The key elements underpinning the current work are a robust tracking algorithm that can identify the encompassing mesh elements for individual particles and a point to point MPI frame-work. These developments combined with accurate property models and high-order time-stepping deliver realistic simulations on larger computational grids for a wide range of particle mass loadings. The coupled US3D-particle code will be applied to the Schiarparelli capsule and Earth-based experimental test cases.

A Sahai↗

NASA NTRS1984

Numerical simulation of the viscous flow fields over three-dimensional complicated geometries

A 'thin-layer' Navier-Stokes code capable of predicting steady state viscous flows is applied to complicated three-dimensional flow fields. The code is written in a generalized coordinate system, and a recently developed grid generation procedure is used for the flow-field discretization. Application is made to the vortical flow over a delta wing at high angle of attack, and the computed results are compared with experimental results. The results indicate that the present method can capture the physical phenomenon well. It is observed that a leading-edge separation vortex is formed over the wing, as is a secondary separation vortex near the leading edge. The flow field behind the trailing edge is also well described. Application is also made to the transonic flow over the Shuttle configuration. The result appears to be reasonable even though no experimental data are available for comparison. These results indicate that the present approach is capable of computing complicated three-dimensional flow fields.

Fujii, K.↗

NASA NTRS1984

Digital image processing of Seabeam bathymetric data for structural studies of seamounts near the East Pacific Rise

The problem of displaying information on the seafloor morphology is attacked by utilizing digital image processing techniques to generate images for Seabeam data covering three young seamounts on the eastern flank of the East Pacific Rise. Errors in locations between crossing tracks are corrected by interactively identifying features and translating tracks relative to a control track. Spatial interpolation techniques using moving averages are used to interpolate between gridded depth values to produce images in shaded relief and color-coded forms. The digitally processed images clarify the structural control on seamount growth and clearly show the lateral extent of volcanic materials, including the distribution and fault control of subsidiary volcanic constructional features. The image presentations also clearly show artifacts related to both residual navigational errors and to depth or location differences that depend on ship heading relative to slope orientation in regions with steep slopes.

Edwards, M. H.↗

NASA NTRS1987

Inviscid transonic wing design using inverse methods in curvilinear coordinates

An inverse wing design method has been developed around an existing transonic wing analysis code. The original analysis code, TAWFIVE, has as its core the numerical potential flow solver, FLO30, developed by Jameson and Caughey. Features of the analysis code include a finite-volume formulation; wing and fuselage fitted, curvilinear grid mesh; and a viscous boundary layer correction that also accounts for viscous wake thickness and curvature. The development of the inverse methods as an extension of previous methods existing for design in Cartesian coordinates is presented. Results are shown for inviscid wing design cases in super-critical flow regimes. The test cases selected also demonstrate the versatility of the design method in designing an entire wing or discontinuous sections of a wing.

Gally, Thomas A.↗

NASA NTRS1988

Analysis of the convergence history of flow through nozzles with shocks

Acceleration techniques such as Wynn's (1986) epsilon algorithm and analysis techiques such as eigensystem analysis are used here to study numerically the convergence properties of an iterative scheme applied to the quasi-one-dimensional Euler and Navier-Stokes equations for flow through nozzles with shocks. The convergence and stability properties are studied by analyzing the dependence of convergence of the code on the discretization technique, boundary conditions, time-step, number of grid points, and the physics of the problem.

Cheer, A. Y.↗

NASA NTRS1989

Development of direct-inverse 3-D methods for applied transonic aerodynamic wing design and analysis

An inverse wing design method was developed around an existing transonic wing analysis code. The original analysis code, TAWFIVE, has as its core the numerical potential flow solver, FLO30, developed by Jameson and Caughey. Features of the analysis code include a finite-volume formulation; wing and fuselage fitted, curvilinear grid mesh; and a viscous boundary layer correction that also accounts for viscous wake thickness and curvature. The development of the inverse methods as an extension of previous methods existing for design in Cartesian coordinates is presented. Results are shown for inviscid wing design cases in super-critical flow regimes. The test cases selected also demonstrate the versatility of the design method in designing an entire wing or discontinuous sections of a wing.

Carlson, Leland A.↗

NASA NTRS1990

Simulation of disk-disk encounters with co-moving polar grids

The two-grid simulation method combining advantages of both polar and Cartesian mesh-codes is described. In addition to the stellar component reacting solely to gravitational forces, the gas component is included with dissipatively colliding particles. This allows fairly realistic simulation of planar encounters where both systems contain star plus gas disks.

Salo, Heikki↗

NASA NTRS1992

Application of engineering codes to the assessment of spacecraft charging induced hazards

One of the major hurdles to developing a better understanding of the environmentally induced interactions in spacecraft systems (spacecraft charging) is the lack of simple tools to analyze the behavior of the spacecraft in a sub-storm environment. The available tool, NASA Charging Analyzer Program (NASCAP), is expensive to use and treats only half of the problem - surface effects only. In response to this need, we have developed, under our IRAD Program, a set of engineering tools that operate on personal computers and addresses both the surface and bulk charging aspects of the interactions. These engineering codes utilize circuit analysis techniques and model simplifications to allow rapid computations. The codes were compared to both NASCAP predictions and to bulk charging data with excellent agreement for both phenomena. The models used in our engineering codes are shown, and in the engineering codes, the models are defined in terms of plates or triangles in an unconfined space. The predictions are based upon current balances to each plate. The code has routines to allow interactions between the plates to simulate 'barrier' effects. The agreement for the predicted differential surface voltages for the model is given. In conducting the comparison of engineering codes to NASCAP predictions, it was found that some inconsistencies existed in the NASCAP code that are not obvious to the user. The use of the grid imposes constraints that voltage gradients exist between the grid points. This allows the existence of voltage gradients in conductors and influences surface resistivity predictions. These points are discussed along with the predictions for treating these inconsistencies properly. The model predictions for a variable intensity substorm is also presented.

Stevens, N. J.↗

NASA NTRS1992

Full Navier-Stokes analysis of a two-dimensional mixer/ejector nozzle for noise suppression

A three-dimensional full Navier-Stokes (FNS) analysis was performed on a mixer/ejector nozzle designed to reduce the jet noise created at takeoff by a future supersonic transport. The PARC3D computational fluid dynamics (CFD) code was used to study the flow field of the nozzle. The grid that was used in the analysis consisted of approximately 900,000 node points contained in eight grid blocks. Two nozzle configurations were studied: a constant area mixing section and a diverging mixing section. Data are presented for predictions of pressure, velocity, and total temperature distributions and for evaluations of internal performance and mixing effectiveness. The analysis provided good insight into the behavior of the flow.

Debonis, James R.↗

NASA NTRS1992

Full Navier-Stokes analysis of a two-dimensional mixer/ejector nozzle for noise suppression

Debonis, James R.↗

NASA NTRS1992

Gigaflop (billion floating point operations per second) performance for computational electromagnetics

Accurate and rapid evaluation of radar signature for alternative aircraft/store configurations would be of substantial benefit in the evolution of integrated designs that meet radar cross-section (RCS) requirements across the threat spectrum. Finite-volume time domain methods offer the possibility of modeling the whole aircraft, including penetrable regions and stores, at longer wavelengths on today's gigaflop supercomputers and at typical airborne radar wavelengths on the teraflop computers of tomorrow. A structured-grid finite-volume time domain computational fluid dynamics (CFD)-based RCS code has been developed at the Rockwell Science Center, and this code incorporates modeling techniques for general radar absorbing materials and structures. Using this work as a base, the goal of the CFD-based CEM effort is to define, implement and evaluate various code development issues suitable for rapid prototype signature prediction.

Shankar, V.↗

NASA NTRS1993

Carbon and nitrogen in Type 2 supernova diamonds

Abundant diamonds found in meteorites seem either to have condensed within supernova interiors during their expansions and coolings or to have been present around those explosions. Either alternative allows implantation of Xe-HL prior to interstellar mixing. A puzzling feature is the near normalcy of the carbon isotopes, considering that the only C-rich matter, the He-burning shell, is pure C-12 in that region. That last fact has caused many to associate supernova carbon with C-12 carbon, so that its SUNOCONS have been anticipated as very C-12-rich. We show that this expectation is misleading because the C-13-rich regions of Type 2's have been largely overlooked in this thinking. We here follow the idea that the diamonds nucleated in the C-12-rich He shell, the only C-rich site for nucleation, but then attached C-13-rich carbon during turbulent encounters with overlying C-13-rich matter. That is, the initial diamonds continued to grow during the same collisional encounters that cause the Xe-HL implantation. Instead of interacting with the small carbon mass having 13/12 = 0.2 in the upper He zone, however, we have calculated the remnants of the initial H-burning core, which left behind C-13-rich matter as it receded during core hydrogen burning. Howard et al. described why the velocity mixing would be essential to understanding the implantation of both the Xe-H and Xe-L components. Velocity mixing is now known to occur from the X-ray and gamma-ray light curves of supernova 1987A. Using the stellar evolution code developed at Goettingen, we calculated at Clemson the evolution of a grid of massive stars up to the beginning of core He burning. We paid attention to all H-burning reactions throughout the star, to the treatment of both convection and semiconvection, and to the recession of the outer boundary of the convective H-burning core as the star expands toward a larger redder state. This program was to generate a careful map of the CNO isotope distribution as He burning begins. Our result for the 30 solar mass star is shown.

Clayton, Donald D.↗