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 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Three-dimensional modeling of diesel engine intake flow, combustion and emissions-2

A three-dimensional computer code, KIVA, is being modified to include state-of-the-art submodels for diesel engine flow and combustion. Improved and/or new submodels which have already been implemented and previously reported are: wall heat transfer with unsteadiness and compressibility, laminar-turbulent characteristic time combustion with unburned HC and Zeldo'vich NO(x), and spray/wall impingement with rebounding and sliding drops. Progress on the implementation of improved spray drop drag and drop breakup models, the formulation and testing of a multistep kinetics ignition model, and preliminary soot modeling results are described. In addition, the use of a block structured version of KIVA to model the intake flow process is described. A grid generation scheme was developed for modeling realistic (complex) engine geometries, and computations were made of intake flow in the ports and combustion chamber of a two-intake-value engine. The research also involves the use of the code to assess the effects of subprocesses on diesel engine performance. The accuracy of the predictions is being tested by comparisons with engine experiments. To date, comparisons were made with measured engine cylinder pressure, temperature and heat flux data, and the model results are in good agreement with the experiments. Work is in progress that will allow validation of in-cylinder flow and soot formation predictions. An engine test facility is described that is being used to provide the needed validation data. Test results were obtained showing the effect of injection rate and split injections on engine performance and emissions.

Reitz, R. D.↗

NASA NTRSDate not supplied

Aeroheating Predictions for a Hypersonic, Turbulent Near-Wake

The accuracy of heating predictions using various turbulence models is examined for an axisymmetric near-wake at Mach 6. The CFD predictions are compared with experimental data collected under AGARD Working Group 18 on the wake of a 70-degree sphere-cone. The impact of grid resolution and discretization error is estimated, which allows a comparison of stacked-block and conventional structured meshes. The accuracy of steady Reynolds-averaged Navier-Stokes (RANS) models is contrasted with that of a hybrid RANS/Large-Eddy Simulation model. The predictions are made with three different CFD codes (LAURA, FUN3D, and HyperSolve), to demonstrate the code-to-code variation in the results. Steady SST models capture the qualitative nature of the heating in the wake, including the increasing peak heating with increasing Reynolds number. Spalart-Allmaras models, including SA-Catris, under-predicted the peak heating in the wake. Hybrid RANS/LES models improve upon the SA results but have their own modeling difficulties near the shear layer impingement. These results are generally consistent across solvers and grid topologies.

RANS↗

NASA NTRSDate not supplied

Aeroheating Predictions for a Hypersonic, Turbulent Near-Wake

RANS↗

NASA NTRS1985

Unsteady transonic flow calculations for interfering lifting surface configurations

Unsteady transonic flow calculations are presented for aerodynamically interfering lifting surface configurations. Calculations are performed by extending the XTRAN3S (Version 1.5) unsteady transonic small-disturbance code to allow the treatment of an additional lifting surface. The research was conducted as a first-step toward developing the capability to treat a complete flight vehicle. Grid generation procedures for swept tapered interfering lifting surface applications of XTRAN3S are described. Transonic calculations are presented for wing-tail and canard-wing configurations for several values of mean angle of attack. The effects of aerodynamic interference on transonic steady pressure distributions and steady and oscillatory spanwise lift distributions are demonstrated. Results due to wing, tail, or canard pitching motions are presented and discussed in detail.

Batina, J. T.↗

NASA NTRS1985

Unsteady transonic flow calculations for interfering lifting surface configurations

Batina, J. T.↗

NASA NTRS1990

Prediction of the aerodynamic and thermal environment in turbines

A numerical study of the steady, viscous flow prediction capabilities of the three-dimensional turbine stage code ROTOR3 is presented. Computations were performed with RAI3DC, a cascade version of ROTOR3 capable of being run in a planar or annular mode. Computed results are compared with experimental data obtained for Hodson's cascade, Kopper's cascade, and United Technologies Research Center's Large Scale Rotating Rig (LSRR) first-stage stator. The code's predictive capability is assessed in terms of the accuracy of predicted airfoil loadings, performance (including secondary flows in the LSRR case), boundary layers, and heat transfer. A grid refinement study was conducted in the LSRR case in an effort to more accurately model the boundary layers on the airfoil and endwall surfaces. The effects of the inlet total pressure profile in secondary flow prediction were also assessed.

Griffin, Lisa W.↗

NASA NTRS1992

Efficient simulation of incompressible viscous flow over single and multi-element airfoils

Incompressible viscous turbulent flows over single- and multiple-element airfoils are numerically simulated in an efficient manner by solving the incompressible Navier-Stokes equations. The solution algorithm uses the method of pseudocompressibility with an upwind-differencing scheme for the convective fluxes and an implicit line-relaxation scheme to study high-lift take-off and landing configurations and to compute lift and drag at various angles of attack up to stall. Two different turbulence models are tested in computing the flow over an NACA 4412 airfoil. The approach used for multiple-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.↗

NASA NTRS2012

LSPRAY-IV: A Lagrangian Spray Module

LSPRAY-IV is a Lagrangian spray solver developed for application with parallel computing and unstructured grids. It is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo Probability Density Function (PDF) solvers. The solver accommodates the use of an unstructured mesh with mixed elements of either triangular, quadrilateral, and/or tetrahedral type for the gas flow grid representation. It is mainly designed to predict the flow, thermal and transport properties of a rapidly vaporizing spray. Some important research areas covered as a part of the code development are: (1) the extension of combined CFD/scalar-Monte- Carlo-PDF method to spray modeling, (2) the multi-component liquid spray modeling, and (3) the assessment of various atomization models used in spray calculations. The current version contains the extension to the modeling of superheated sprays. The manual provides the user with an understanding of various models involved in the spray formulation, its code structure and solution algorithm, and various other issues related to parallelization and its coupling with other solvers.

Raju, M. S.↗

NASA NTRS2014

A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de

We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.

electric propulsion↗

NASA NTRS1992

Development of a multigrid transonic potential flow code for cascades

Finite-volume methods for discretizing transonic potential flow equations have proven to be very flexible and accurate for both two and three dimensional problems. Since they only use local properties of the mapping, they allow decoupling of the grid generation from the rest of the problem. A very effective method for solving the discretized equations and converging to a solution is the multigrid-ADI technique. It has been successfully applied to airfoil problems where O type, C type and slit mappings have been used. Convergence rates for these cases are more than an order of magnitude faster than with relaxation techniques. In this report, we describe a method to extend the above methods, with the C type mappings, to airfoil cascade problems.

Steinhoff, John↗

NASA NTRS1981

Transonic flow over wing/fuselage configurations

Wing Body Code (WIBCO) program simulates flow-field configurations for reduction of design cost and improvement of aircraft performance. Inputs to WIBCO consist of ambient flow conditions and geometric configuration data; grid control and relaxation parameters are internally set. Outputs include input data echo, grid system verification, relaxation-solution convergence history, and computed velocities, pressures, forces, moments, reference lengths, and areas. Program is written in FORTRAN IV for batch execution.

Boppe, C. W.↗

NASA NTRS1989

Analysis of the SSME HPOTP bearing inlet cavity

Analysis of the flow in the Space Shuttle Main Engine (SSME) high pressure oxygen turbopump (HPOTP) bearing no. 1 inlet cavity was completed in support of return-to-flight. With the incorporation of several design changes in the Phase 2 turbopump, rotordynamic stability of the pumps was enhanced, but the durability and life of the LOX-cooled bearings has decreased. During the post-Challenger SSME recertification, the causes of limited bearing durability were investigated. One topic addressed was the flow environment upstream of the pump-end bearing and the effect of seal exit swirl and a cavity anti-vortex rib on the bearing environment and life. The objective is to define the hydrodynamic environment upstream of the pump-end bearing and determine the effect of seal exit swirl and the anti-vortex rib on bearing inlet swirl. The problem was posed as an axisymmetric cavity flow with the computational domain extending from the seal exit to the bearing inlet. This domain was discretized with 22800 grid points. Boundary conditions were obtained from a 1-D model of the SSME coolant path. The inlet Mach number was 0.19 and the problem was solved with the CMINT code utilizing the Briley-McDonald/Beam-Warming algorithm with preconditioning to speed convergence at low Mach numbers. Three parametric cases with inlet swirl of 50 percent shaft speed (labyrinth seal), 20 percent shaft speed (damping seal), and no inlet swirl were considered. Computational results indicate large vortical flow structures in the cavity, with the labyrinth, damping, and no-swirl cases yielding bearing inlet swirl rates of 14, 10, and 9 percent of shaft speed, respectively. When these results were used as input to the SHABRETH bearing model, limited durability could not be explained by these small differences in swirl. Also, based on these results, a proposed design change for the cavity anti-vortex rib was not implemented by the SSME chief engineer.

Mcconnaughey, P. K.↗

NASA NTRS1987

Design of Optimum Ducts Using an Efficient 3-D Viscous Computational Flow Analysis

Design of fluid dynamically efficient ducts is addressed through the combination of an optimization analysis with a three-dimensional viscous fluid dynamic analysis code. For efficiency, a parabolic fluid dynamic analysis was used. Since each function evaluation in an optimization analysis is a full three-dimensional viscous flow analysis requiring 200,000 grid points, it is important to use both an efficient fluid dynamic analysis and an efficient optimization technique. Three optimization techniques are evaluated on a series of test functions. The Quasi-Newton (BFGS, eta = .9) technique was selected as the preferred technique. A series of basic duct design problems are performed. On a two-parameter optimization problem, the BFGS technique is demonstrated to require half as many function evaluations as a steepest descent technique.

Madabhushi, Ravi K.↗

NASA NTRS1995

Technical report series on global modeling and data assimilation. Volume 5: Documentation of the AIRES/GEOS dynamical core, version 2

A detailed description of the numerical formulation of Version 2 of the ARIES/GEOS 'dynamical core' is presented. This code is a nearly 'plug-compatible' dynamics for use in atmospheric general circulation models (GCMs). It is a finite difference model on a staggered latitude-longitude C-grid. It uses second-order differences for all terms except the advection of vorticity by the rotation part of the flow, which is done at fourth-order accuracy. This dynamical core is currently being used in the climate (ARIES) and data assimilation (GEOS) GCMs at Goddard.

Suarez, Max J.↗

NASA NTRS2001

Method Improvements in Thermal Analysis of Mach 10 Leading Edges

Several improvements have recently been made in the thermal analysis methods for leading edges of a hypersonic vehicle. The leading edges of this vehicle undergo exceptionally high heat loads that incorporate extreme spatial gradients as well as severe transients. Due to the varying flight conditions, complex geometry, and need for thermal loads at many points along the trajectory, full computational fluid dynamics (CFD) analysis of the aeroheating; loads is not feasible. Thus, engineering methods must be used to determine the aeroheating on the vehicle surfaces, and that must be utilized in the thermal analysis. Over the last year, the thermal analysis of a hypersonic vehicle has been enhanced in several ways. Two different engineering codes are used to predict aeroheating loads: one over the curve near the stagnation point, and the other on flat surfaces downstream of the leading edge. These two are matched together at the intersection point using a method that allows closer approximation of CFD results. User-developed FORTRAN, which is part of the thermal solver PATRAN Thermal, is used to accomplish this. The customizable FORTRAN code also allows use of many different time- and space-dependent factors, interpolation of the heat load in time and space, and inclusion of both highly swept and unswept grid structures. This FORTRAN is available to other PATRAN users who may want to accomplish a similar objective in analysis. Flux, rather than convective coefficient, is used to define heat loads, which allows more accurate analysis as well as better application of margins. Improvements have also been made in more efficient utilization of imported CAD geometry, by creating faces on solids to facilitate load application.

Amundsen, Ruth M.↗

NASA NTRS2002

Design and Analysis of Turbomachinery for Space Applications

This presentation provides an overview of CORSAIR, a three dimensional computational fluid dynamics software code for the analysis of turbomachinery components available from NASA, and discusses its potential use in the design of these parts. Topics covered include: time-dependent equations of motion, grid topology, turbulence models, boundary conditions, parallel simulations and miscellaneous capabilities.

Dorney, D.↗

NASA NTRS2018

Comparison of Staggered Grid Finite Difference Schemes for Ultrasound Simulation in Curving Composites

The optimization of ultrasonic nondestructive evaluation (NDE) simulation tools for composites has the potential to reduce both individual part inspection time and overall certification time for composite parts and structures. Inspection guidance based on simulation provides increased confidence in the veracity of inspection results in addition to time reductions. This paper outlines ongoing work targeted to advance this objective through the use of finite difference (FD) simulation techniques formulated for composite structures with realistic geometries. Two staggered grid explicit FD schemes which show promise for this purpose are assessed: the Lebedev FD scheme and the rotated staggered grid (RSG) FD scheme. Algorithmic points which provide challenges for complex geometries are addressed, in particular handling of traction free surfaces and bi-material interfaces present at lamina boundaries. Code execution time estimates are performed as well to guide feasible domain sizes relative to algorithm choice and available hardware. Three test cases are simulated: a delaminated plate, a cylinder, and a triclinic lamina. These tests demonstrate that the Lebedev FD scheme needs additional work to handle inter-laminar interfaces and traction free boundaries in the presence of stair-stepping approximations. In contrast, the simple structure of the RSG unit cell makes it more straightforward to construct a 3D simulation technique for curved composite laminates.

Frankforter, Erik↗

NASA NTRSDate not supplied

NASA LaRC Contribution to the High Angle Working Group of the Third Aeroelastic Prediction Workshop: BSCW Shock Buffet

FUN3D Core Capabilities - Established as a research code in late 1980s; now supports numerous internal and external efforts across the speed range - Solves 2D/3D steady and unsteady Euler and RANS equations on node-based mixed element grids for compressible and incompressible flows - General dynamic mesh capability: any combination of rigid / overset / morphing grids, including 6-DOF effects - Aeroelastic modeling using mode shapes, full FEM, etc. - Constrained / multipoint adjoint-based design and mesh adaptation - Distributed development team using agile/extreme software practices including 24/7 regression, performance testing - Capabilities fully integrated, online documentation, training videos, tutorials

Pawel Chwalowski↗