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,099 records · Page 61

Navier-Stokes computations and experimental comparisons for multielement airfoil configurations

A two-dimensional unstructured Navier-Stokes code is utilized for computing the flow around multielement configurations. Comparisons are shown for a landing configuration with an advanced-technology flap as well as for a takeoff configuration. For each, a grid convergence study is conducted to assess inaccuracies caused by inadequate grid resolution. Although adequate resolution is obtained for determining the pressure distributions, further refinement is needed to sufficiently resolve the velocity profiles at high angles of attack. For the advanced flap configuration, comparisons of pressure distributions and lift are made with experimental data. Here, two flap riggings and two Reynolds numbers are considered. In general, the trends caused by variations in these quantities are well predicted by the computations, although the angle of attack for maximum lift is overpredicted. For the takeoff configuration, comparisons between pressures, velocity profiles, and Reynolds stresses are made at a single angle of attack with generally good agreement.

Anderson, W. K.↗

Development of an Atomization Methodology for Spray Combustion

In liquid rocket propulsion, the knowledge and the understanding of liquid-gas interfacial phenomena are very important. This is important for predicting the onset of cavitation occurring in swirl injection elements used in STME, as well as atomization processes in shear-induced injectors (co-axial) and impinging injector elements. From the fact that all the physical processes including droplet size distribution, droplet dispersion, mixing and combustion are controlled by atomization processes, it is expected that the successful incorporation of the volume of fraction (VOF) will greatly enhance the analytical capability of predicting spray combustion processes in liquid-fueled engines. In this paper, a methodology is developed to define and track interfaces between two fluids in non-orthogonal, body-fitted grids using a single fractional volume of fluid (VOF) variable to describe the distribution of the liquid phase in a gas-liquid flow field. This method was implemented in a mature CFD code MAST (Multiphase All-Speed Transient) utilizing the general PISO-C algorithm. For the preliminary study on the analysis of spray combustion and tracking of the interface between two phases, we will report on the progress of the simulation of the instability on the liquid column; the surface wave instability and the droplet breakup from the liquid surface.

Seung, S. P.↗

SCISEAL: A CFD Code for Analysis of Fluid Dynamic Forces in Seals

A 3D CFD code, SCISEAL, has been developed and validated. Its capabilities include cylindrical seals, and it is employed on labyrinth seals, rim seals, and disc cavities. State-of-the-art numerical methods include colocated grids, high-order differencing, and turbulence models which account for wall roughness. SCISEAL computes efficient solutions for complicated flow geometries and seal-specific capabilities (rotor loads, torques, etc.).

Althavale, Mahesh M.↗

Time-Dependent Simulations of Incompressible Flow in a Turbopump Using Overset Grid Approach

This viewgraph presentation provides information on mathematical modelling of the SSME (space shuttle main engine). The unsteady SSME-rig1 start-up procedure from the pump at rest has been initiated by using 34.3 million grid points. The computational model for the SSME-rig1 has been completed. Moving boundary capability is obtained by using DCF module in OVERFLOW-D. MPI (Message Passing Interface)/OpenMP hybrid parallel code has been benchmarked.

Kiris, Cetin↗

New Software for Predicting Charging of Spacecraft

Charging System Analyzer Program (Nascap-2K) is a comprehensive update, revision, and extension of several NASA and Air Force codes for predicting electrical charging of spacecraft. Nascap-2K integrates the capabilities and models included in four independent programs: NASCAP/LEO for low-Earth orbits, NASCAP/GEO for geosynchronous orbits, POLAR for auroral charging in polar orbits, and DynaPAC (Dynamic Plasma Analysis Code) for time-dependent plasma interactions. While each of the earlier codes works well for the range of problems for which it was designed, by today s standards these codes are difficult to learn, cumbersome to use, and overly restrictive in their geometric modeling capabilities. Nascap-2K incorporates these models into a single software package that includes spacecraft surface modeling, spatial gridding, environmental specifications, calculating scripting, and post-processing analysis and visualization. The provided material properties database includes values from earlier programs as well as values from recent measurements. Development of Nascap-2K continues with future capabilities to include interactions with dense plasma such as those produced by electric propulsion.

Mandell, Myron↗

Numerical Analysis of a Rotating Detonation Engine in the Relative Reference Frame

A two-dimensional, computational fluid dynamic (CFD) simulation of a semi-idealized rotating detonation engine (RDE) is described. The simulation operates in the detonation frame of reference and utilizes a relatively coarse grid such that only the essential primary flow field structure is captured. This construction yields rapidly converging, steady solutions. Results from the simulation are compared to those from a more complex and refined code, and found to be in reasonable agreement. The performance impacts of several RDE design parameters are then examined. Finally, for a particular RDE configuration, it is found that direct performance comparison can be made with a straight-tube pulse detonation engine (PDE). Results show that they are essentially equivalent.

pulse detonation engines↗

Numerical Analysis of a Rotating Detonation Engine in the Relative Reference Frame

A two-dimensional, computational fluid dynamic (CFD) simulation of a semi-idealized rotating detonation engine (RDE) is described. The simulation operates in the detonation frame of reference and utilizes a relatively coarse grid such that only the essential primary flow field structure is captured. This construction yields rapidly converging, steady solutions. Results from the simulation are compared to those from a more complex and refined code, and found to be in reasonable agreement. The performance impacts of several RDE design parameters are then examined. Finally, for a particular RDE configuration, it is found that direct performance comparison can be made with a straight-tube pulse detonation engine (PDE). Results show that they are essentially equivalent.

propulsion↗

The 4th order GISS model of the global atmosphere

The new GISS 4th order model of the global atmosphere is described. It is based on 4th order quadratically conservative differences with the periodic application of a 16th order filter on the sea level pressure and potential temperature equations, a combination which is approximately enstrophy conserving. Several short range forecasts indicate a significant improvement over 2nd order forecasts with the same resolution (approximately 400 km). However the 4th order forecasts are somewhat inferior to 2nd order forecasts with double resolution. This is probably due to the presence of short waves in the range between 1000 km and 2000 km, which are computed more accurately by the 2nd order high resolution model. An operation count of the schemes indicates that with similar code optimization, the 4th order model will require approximately the same amount of computer time as the 2nd order model with the same resolution. It is estimated that the 4th order model with a grid size of 200 km provides enough accuracy to make horizontal truncation errors negligible over a period of a week for all synoptic scales (waves longer than 1000 km).

Kalnay-Rivas, E.↗

Computation of transonic separated wing flows using an Euler/Navier-Stokes zonal approach

A computer program called Transonic Navier Stokes (TNS) has been developed which solves the Euler/Navier-Stokes equations around wings using a zonal grid approach. In the present zonal scheme, the physical domain of interest is divided into several subdomains called zones and the governing equations are solved interactively. The advantages of the Zonal Grid approach are as follows: (1) the grid for any subdomain can be generated easily; (2) grids can be, in a sense, adapted to the solution; (3) different equation sets can be used in different zones; and, (4) this approach allows for a convenient data base organization scheme. Using this code, separated flows on a NACA 0012 section wing and on the NASA Ames WING C have been computed. First, the effects of turbulence and artificial dissipation models incorporated into the code are assessed by comparing the TNS results with other CFD codes and experiments. Then a series of flow cases is described where data are available. The computed results, including cases with shock-induced separation, are in good agreement with experimental data. Finally, some futuristic cases are presented to demonstrate the abilities of the code for massively separated cases which do not have experimental data.

Kaynak, Uenver↗

A numerical investigation of turbulent baseflow

The turbulent base flow of a cone with a sting in hypersonic flow is investigated using a time-dependent Navier-Stokes code called UWIN. A PNS solver provides the upstream flow conditions which drive the development of the base flow. The computational method is described along with its application to the base flow problem. The grid generation and boundary conditions are described. The results and the summary conclude this note.

Chaussee, Denny S.↗

ISOLEV - A level surface cutting plane program for fluid flow data

A computer program, ISOLEV, allows interactive visualization of CFD scalar and vector functions by means of user-specified sections of the data. The application is based on table lookups which govern both isosurface generation on hexahedral grid cells and recursive subdivision of the cells. The program supports Gouraud-shaded color maps of the data, surface-on-surface maps, and deformation surfaces of vector fields. The execution of the code for animated sweeps is improved by resorting to the cells in the database and maintaining an active set of cells to be rendered.

Kerlick, G. David↗

Unsteady-flow-field predictions for oscillating cascades

The unsteady flow field around an oscillating cascade of flat plates with zero stagger was studied by using a time marching Euler code. This case had an exact solution based on linear theory and served as a model problem for studying pressure wave propagation in the numerical solution. The importance of using proper unsteady boundary conditions, grid resolution, and time step size was shown for a moderate reduced frequency. Results show that an approximate nonreflecting 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 to the airfoils than when reflective boundaries are used. Stretching the boundary to dampen the unsteady waves is another way to minimize reflections. Grid clustering near the plates captures the unsteady flow field better than when uniform grids are used as long as the 'Courant Friedrichs Levy' (CFL) number is less than 1 for a sufficient portion of the grid. Finally, a solution based on an optimization of grid, CFL number, and boundary conditions shows good agreement with linear theory.

Huff, Dennis L.↗

Computations of spray, fuel-air mixing, and combustion in a lean-premixed-prevaporized combustor

A code was developed for computing the multidimensional flow, spray, combustion, and pollutant formation inside gas turbine combustors. The code developed is based on a Lagrangian-Eulerian formulation and utilizes an implicit finite-volume method. The focus of this paper is on the spray part of the code (both formulation and algorithm), and a number of issues related to the computation of sprays and fuel-air mixing in a lean-premixed-prevaporized combustor. The issues addressed include: (1) how grid spacings affect the diffusion of evaporated fuel, and (2) how spurious modes can arise through modelling of the spray in the Lagrangian computations. An upwind interpolation scheme is proposed to account for some effects of grid spacing on the artificial diffusion of the evaporated fuel. Also, some guidelines are presented to minimize errors associated with the spurious modes.

Dasgupta, A.↗

Hypersonic Shock Interactions About a 25 deg/65 deg Sharp Double Cone

This paper presents the results of a numerical study of shock interactions resulting from Mach 10 air flow about a sharp double cone. Computations are made with the direct simulation Monte Carlo (DSMC) method by using two different codes: the G2 code of Bird and the DAC (DSMC Analysis Code) code of LeBeau. The flow conditions are the pretest nominal free-stream conditions specified for the ONERA R5Ch low-density wind tunnel. The focus is on the sensitivity of the interactions to grid resolution while providing information concerning the flow structure and surface results for the extent of separation, heating, pressure, and skin friction.

Moss, James N.↗

Development Status of the Phantom Code for Turbomachinery

Contents include the following: Motivation. Code details. Test cases: aardvark (2D), phantom (3D). Algorithm basics. New feature I, II, and III. UTRC LSRR 1 - 1/2 stage turbine. Supersonic turbine. UTRC LSRR single-stage turbine. LOX elbow duct. Computational grid. High-pressure fuel pump inducer. Development plan. Summary. Appendix - Background.

Dorney, D. J.↗

Influence of Turbulence on the Restraint of Liquid Jets by Surface Tension in Microgravity Investigated

Microgravity poses many challenges to the designer of spacecraft tanks. Chief among these are the lack of phase separation and the need to supply vapor-free liquid or liquidfree vapor to the spacecraft processes that require fluid. One of the principal problems of phase separation is the creation of liquid jets. A jet can be created by liquid filling, settling of the fluid to one end of the tank, or even closing a valve to stop the liquid flow. Anyone who has seen a fountain knows that jets occur in normal gravity also. However, in normal gravity, the gravity controls and restricts the jet flow. In microgravity, with gravity largely absent, surface tension forces must contain jets. To model this phenomenon, a numerical method that tracks the fluid motion and the surface tension forces is required. Jacqmin has developed a phase model that converts the discrete surface tension force into a barrier function that peaks at the free surface and decays rapidly away. Previous attempts at this formulation were criticized for smearing the interface. This can be overcome by sharpening the phase function, double gridding the fluid function, and using a higher order solution for the fluid function. The solution of this equation can be rewritten as two coupled Poisson equations that also include the velocity. After the code was implemented in axisymmetric form and verified by several test cases at the NASA Glenn Research Center, the drop tower runs of Aydelott were modeled. Work last year with a laminar model was found to overpredict Aydelott's results, except at the lowest Reynolds number conditions of 400. This year, a simple turbulence model was implemented by adding a turbulent viscosity based on the mixing-length hypothesis and empirical measurements of previous works. Predictions made after this change was implemented have been much closer to experimentally observed flow patterns and geyser heights. Two model runs is shown. The first, without any turbulence correction, breaks through the free surface and strikes the far end of the tank. In the second, the turbulence spreads the jet momentum over more of the free surface, enabling the surface tension forces to turn the jet back into the bulk liquid. The model geyser height with the second model is 1.1 cm. This is quite close to the 1.5-cm geyser height measured by Aydelott.

Chato, David J.↗

NASA Langley FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

This work summarizes the results of FUN3D analyses conducted for the 1st AIAA Stability and Control Workshop on behalf of participants from the NASA Langley Research Center. The workshop was created to establish best practices for the prediction of stability and control derivatives using computational fluid dynamics and assess the limitations of these methods when those best practices are applied. The inaugural workshop considered the ONERA version of the NASA/Boeing Common Research Model, which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel data at small sideslip angles remain unpublished and served as ‘blind’ data for computational comparisons. The present research generated workshop test case data using the NASA FUN3D code, which is a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Steady- state numerical simulations were conducted for workshop test cases investigating the following: grid convergence, Mach number effect on static stability, wind tunnel sting increments, static stability-derivative calculations, and a sideslip angle sweep. Results were generated for two series of unstructured, mixed-element grids, one set provided by the workshop and another set created using the HeldenMesh grid generation software. The results provided include total- and component-level breakdowns of the force and moment coefficients, in addition to sectional pressure distributions for the wing and tail components for comparisons to wind tunnel data.

CFD↗

Prediction of Business Jet Airloads Using The Overflow Navier-Stokes Code

The objective of this work is to evaluate the application of Navier-Stokes computational fluid dynamics technology, for the purpose of predicting off-design condition airloads on a business jet configuration in the transonic regime. The NASA Navier-Stokes flow solver OVERFLOW with Chimera overset grid capability, availability of several numerical schemes and convergence acceleration techniques was selected for this work. A set of scripts which have been compiled to reduce the time required for the grid generation process are described. Several turbulence models are evaluated in the presence of separated flow regions on the wing. Computed results are compared to available wind tunnel data for two Mach numbers and a range of angles-of-attack. Comparisons of wing surface pressure from numerical simulation and wind tunnel measurements show good agreement up to fairly high angles-of-attack.

Bounajem, Elias↗