Search NASA⌕ Search

SEARCH · Search NASA

Results for “Navier Stokes solver”

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 253 records · Page 14

Computation of the flow field in a centrifugal impeller with splitter blades

To support the design effort of the Space Transportation Main Engine (STME) Fuel Pump Stage, viscous flow calculations were performed in a centrifugal impeller with splitter blades. These calculations were carried out with a Navier-Stokes solver (MINT), which employs a linearized block-implicit Alternating Direction Implicit (ADI) procedure to iteratively solve a finite difference form of the system of conservation equations of mass, momentum, and energy in body-fitted coordinates. A computational grid was generated algebraically for the 'channel' between two main blades of the impeller and extended both upstream of the impeller inlet and downstream of the impeller exit so that the appropriate boundary conditions could be applied. The results of the calculations show that although the overall level of flow distortion near the impeller exit is not very large, there is a noticeable difference between the flow patterns in the two 'passages' (one passage between the pressure side of the splitter blade and the suction side of the next full blade).

Dejong, Frederik J.↗

Impeller tandem blade study with grid embedding for local grid refinement

Flow non-uniformity at the discharge of high power density impellers can result in significant unsteady interactions between impeller blades and downstream diffuser vanes. These interactions result in degradation of both performance and pump reliability. The MSFC Pump Technology Team has recognized the importance of resolving this problem and has thus initiated the development and testing of a high head coefficient impeller. One of the primary goals of this program is to improve impeller performance and discharge flow uniformity. The objective of the present work is complimentary. Flow uniformity and performance gains were sought through the application of a tandem blade arrangement. The approach adopted was to numerically establish flow characteristics at the impeller discharge for the baseline MSFC impeller and then parametrically evaluate tandem blade configurations. A tandem design was sought that improves both impeller performance and discharge uniformity. The Navier-Stokes solver AEROVISC was used to conduct the study. Grid embedding is used to resolve local gradients while attempting to minimize model size. Initial results indicate that significant gains in flow uniformity can be achieved through the tandem blade concept and that blade clocking rather than slot location is the primary driver for flow uniformity.

Bache, George↗

Numerical prediction of wakes in cascades and compressor rotors including the effects of mixing. II - Rotor passage flow and wakes including the effects of spanwise mixing

The results of a numerical investigation to predict the flow field including wakes and mixing in axial-flow compressor rotors are presented. The wake behavior in a moderately loaded compressor rotor is studied numerically using a 3D incompressible Navier-Stokes solver with a high Reynolds number form of a turbulence model. The equations are solved using a time dependent implicit technique. The agreement between the measured data and the predictions is good; including the blade boundary-layer profiles, wake mean-velocity profiles, and decay. The ability of the pseudocompressibility scheme to predict the entire flow field including the near and far wake profiles and its decay characteristics, effect of loading, and the viscous losses of a 3D rotor flow field are demonstrated. The mixing in the downstream regions away from the hub and annulus walls is dominated by wake diffusion. In regions away from the walls the radial mixing is predominantly caused by the transport of mass, momentum, and energy by the radial component of velocity in the wake.

Lakshminarayana, B.↗

Turbulent flow calculations using unstructured and adaptive meshes

A method of efficiently computing turbulent compressible flow over complex two dimensional configurations is presented. The method makes use of fully unstructured meshes throughout the entire flow field, thus enabling the treatment of arbitrarily complex geometries and the use of adaptive meshing techniques throughout both viscous and inviscid regions of the flow field. Mesh generation is based on a locally mapped Delaunay technique in order to generate unstructured meshes with highly-stretched elements in the viscous regions. The flow equations are discretized using a finite element Navier-Stokes solver, and rapid convergence to steady-state is achieved using an unstructured multigrid algorithm. Turbulence modeling is performed using an inexpensive algebraic model, implemented for use on unstructured and adaptaive meshes. Compressible turbulent flow solutions about multiple-element airfoil geometries are computed and compared with expeirmental data.

Mavriplis, Dimitri J.↗

An assumed pdf approach for the calculation of supersonic mixing layers

In an effort to predict the effect that turbulent mixing has on the extent of combustion, a one-equation turbulence model is added to an existing Navier-Stokes solver with finite-rate chemistry. To average the chemical-source terms appearing in the species-continuity equations, an assumed pdf approach is also used. This code was used to analyze the mixing and combustion caused by the mixing layer formed by supersonic coaxial H2-air streams. The chemistry model employed allows for the formation of H2O2 and HO2. Comparisons are made with recent measurements using laser Raman diagnostics. Comparisons include temperature and its rms, and concentrations of H2, O2, N2, H2O, and OH. In general, good agreement with experiment was noted.

Baurle, R. A.↗

Turbulent reacting flow computations including turbulence-chemistry interactions

A two-equation (k-epsilon) turbulence model has been extended to be applicable for compressible reacting flows. A compressibility correction model based on modeling the dilatational terms in the Reynolds stress equations has been used. A turbulence-chemistry interaction model is outlined. In this model, the effects of temperature and species mass concentrations fluctuations on the species mass production rates are decoupled. The effect of temperature fluctuations is modeled via a moment model, and the effect of concentration fluctuations is included using an assumed beta-pdf model. Preliminary results obtained using this model are presented. A two-dimensional reacting mixing layer has been used as a test case. Computations are carried out using the Navier-Stokes solver SPARK using a finite rate chemistry model for hydrogen-air combustion.

Narayan, J. R.↗

Single expansion ramp nozzle simulations

The single-expansion-ramp-nozzle (SERN) experiment underway at NASA Ames Research Center simulates the National Aerospace Plane propulsive jet-plume flow. Recently, limited experimental data has become available from an experiment with a generic nozzle/afterbody model in a hypersonic wind tunnel. The present paper presents full three-dimensional solutions obtained with the implicit Navier-Stokes solver, FL3D, for the baseline model and a version of the model with side extensions. Analysis of the computed flow clearly shows the complex 3-D nature of the flow, critical flow features, and the effect of side extensions on the plume flow development. Flow schematics appropriate for the conditions tested are presented for the baseline model and the model with side extensions. The computed results show excellent agreement with experimental shadowgraph and with surface pressure measurements. The computed and experimental surface oil-flows show the same features but may be improved by appropriate turbulence modeling.

Ruffin, Stephen M.↗

Numerical investigation of performance degradation of wings and rotors due to icing

The aerodynamic load characteristics and the performance degradation of moderate aspect ratio wings and rotors with simulated glaze leading-edge ice have been studied using a three-dimensional, compressible Navier-Stokes solver. The effect of a splitter plate at the wing root on both clean and iced wing configurations has been studied and the results are compared with the experiment. A significant difference has been observed with and without splitter plates in the magnitude of flow separation and aerodynamic loading at the inboard stations for the iced wing at 8-deg angle of attack. Inviscid calculations were performed and compared with viscous calculations to investigate whether the performance of iced swept wings can be inexpensively predicted using Euler methods. It is shown that inviscid calculations predict higher aerodynamic loading than viscous calculations, and cannot model separation effects. A typical nonlifting helicopter rotor in forward flight condition is also studied, and the penalty due to the leading-edge ice formation on the required torque is numerically demonstrated.

Kwon, Oh J.↗

Numerical simulation of vortex unsteadiness on a slender body at high incidence

Numerical results obtained with the present3D thin-layer Navier-Stokes solver are shown to qualitatively follow experimental results. They also indicate the existence of at least two main frequencies on an ogive cylinder at 40 and 60 deg angles-of-attack: (1) a low frequency associated with fluctuations of the primary vortices, and leads to vortex shedding, and (2) a high frequency associated with shear-layer fluctuations that are concentrated in several small areas above the leeward side of the body's cylindrical portion.

Degani, David↗

An Eulerian/Lagrangian method for computing blade/vortex impingement

A combined Eulerian/Lagrangian approach to calculating helicopter rotor flows with concentrated vortices is described. The method computes a general evolving vorticity distribution without any significant numerical diffusion. Concentrated vortices can be accurately propagated over long distances on relatively coarse grids with cores only several grid cells wide. The method is demonstrated for a blade/vortex impingement case in 2D and 3D where a vortex is cut by a rotor blade, and the results are compared to previous 2D calculations involving a fifth-order Navier-Stokes solver on a finer grid.

Steinhoff, John↗

Characteristics of the Shuttle Orbiter leeside flow during a reentry condition

A study of the leeside flow characteristics of the Shuttle Orbiter is presented for a reentry flight condition. The flow is computed using a point-implicit, finite-volume scheme known as the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA). LAURA is a second-order accurate, laminar Navier-Stokes solver, incorporating finite-rate chemistry with a radiative equilibrium wall temperature distribution and finite-rate wall catalysis. The resulting computational solution is analyzed in terms of salient flow features and the surface quantities are compared with flight data.

Kleb, William L.↗

Navier-Stokes turbine heat transfer predictions using two-equation turbulence closures

Navier-Stokes calculations were carried out in order to predict the heat-transfer rates on turbine blades. The calculations were performed using TRAF2D which is a k-epsilon, explicit, finite volume mass-averaged Navier-Stokes solver. Turbulence was modeled using Coakley's q-omega and Chien's k-epsilon two-equation models and the Baldwin-Lomax algebraic model. The model equations along with the flow equations were solved explicitly on a nonperiodic C grid. Implicit residual smoothing (IRS) or a combination of multigrid technique and IRS was applied to enhance convergence rates. Calculations were performed to predict the Stanton number distributions on the first stage vane and blade row as well as the second stage vane row of the SSME high-pressure fuel turbine. The comparison serves to highlight the weaknesses of the turbulence models for use in turbomachinery heat-transfer calculations.

Ameri, Ali A.↗

Numerical modeling of fluid and electromagnetic phenomena in an arcjet

An explicit numerical technique is used to solve the axisymmetric reduced electromagnetic field equation. The effect of an electrical arc on a viscous, axisymmetric flow is approximated using an implicit thin layer Navier-Stokes solver with additional electromagnetic source terms in conjunction with the explicit finite difference code.

Flowe, Anita C.↗

Three dimensional flows in multistage compressors

Current multi-stage turbomachinery design/analysis methods are based on a time-averaged, axisymmetric representation of the flow field. The actual flow field is asymmetric and unsteady due to blade row interactions. The Reynolds averaged Navier-Stokes solvers are limited to single-stage machines for existing computers. Therefore, advanced multi-stage compressors will operate far off-design for portions of the flight regime. The objectives are to provide an experimentally validated average-passage calculation of multistage compressor blade row interactions and an experimentally validated time-accurate calculation of multi-stage compressor blade row interactions. Various aspects of this investigation are presented in viewgraph form.

Hathaway, Michael D.↗

Navier-Stokes turbine heat transfer predictions using two-equation turbulence

Navier-Stokes calculations were carried out in order to predict the heat transfer rates on turbine blades. The calculations were performed using TRAF2D which is a two-dimensional, explicit, finite volume mass-averaged Navier-Stokes solver. Turbulence was modeled using q-omega and k-epsilon two-equation models and the Baldwin-Lomax algebraic model. The model equations along with the flow equations were solved explicitly on a non-periodic C grid. Implicit residual smoothing (IRS) or a combination of multigrid technique and IRS was applied to enhance convergence rates. Calculations were performed to predict the Stanton number distributions on the first stage vane and blade row as well as the second stage vane row of the Rocketdyne Space Shuttle Main Engine (SSME) high pressure fuel turbine. The comparison with the experimental results, although generally favorable, serves to highlight the weaknesses of the turbulence models and the possible areas of improving these models for use in turbomachinery heat transfer calculations.

Ameri, Ali A.↗

Navier-Stokes computation of compressible turbulent flows with a second order closure, part 1

A second order closure turbulence model for compressible flows is developed and implemented in a 2D Reynolds-averaged Navier-Stokes solver. From the beginning where a kappa-epsilon turbulence model was implemented in the bidiagonal implicit method of MACCORMACK (referred to as the MAC3 code) to the final stage of implementing a full second order closure in the efficient line Gauss-Seidel algorithm, numerous work was done, individually and collectively. Besides the collaboration itself, the final product of this work is a second order closure derived from the Launder, Reece, and Rodi model to account for near wall effects, which has been called FRAME model, which stands for FRench-AMerican-Effort. During the reporting period, two different problems were worked out. The first was to provide Ames researchers with a reliable compressible boundary layer code including a wide collection of turbulence models for quick testing of new terms, both in two equations and in second order closure (LRR and FRAME). The second topic was to complete the implementation of the FRAME model in the MAC5 code. The work related to these two different contributions is reported. dilatation in presence of stron shocks. This work, which has been conducted during a work at the Center for Turbulence Research with Zeman aimed also to cros-check earlier assumptions by Rubesin and Vandromme.

Haminh, Hieu↗

Grid adaption using Chimera composite overlapping meshes

The objective of this paper is to perform grid adaptation using composite over-lapping meshes in regions of large gradient to capture the salient features accurately during computation. The Chimera grid scheme, a multiple overset mesh technique, is used in combination with a Navier-Stokes solver. The numerical solution is first converged to a steady state based on an initial coarse mesh. Solution-adaptive enhancement is then performed by using a secondary fine grid system which oversets on top of the base grid in the high-gradient region, but without requiring the mesh boundaries to join in any special way. Communications through boundary interfaces between those separated grids are carried out using tri-linear interpolation. Applications to the Euler equations for shock reflections and to a shock wave/boundary layer interaction problem are tested. With the present method, the salient features are well resolved.

Kao, Kai-Hsiung↗

Calculation of fully three-dimensional separated flow with an unsteady viscous-inviscid interaction method

Three previous papers have shown that the viscous inviscid numerical methods were capable of calculating separated turbulent flows. The 'Semi-inverse method' and the models of the first author allow the computation of massive separations, stalled flows, and shock wave-boundary layer interactions, in two-dimensional or quasi-three-dimensional conditions, from low speeds to supersonic speeds. The 'Semi-implicit' method for time-consistent coupling allows for the computation of time-accurate transonic separated flow and buffer computations in two-dimensions. The present paper shows that the viscous-inviscid interaction approach is also able to compute the fully three-dimensional flow separation phenomena. The method is based on a thin-layer approximation of the theory of 'Defect-Formulation' that provides the viscous-inviscid splitting of the Navier-Stokes equations. A parametric analytical modelling of the 3D-turbulent velocity profiles is involved. Numerically, the 3D-velocity profiles are discretized in the normal z-direction and driven by parametric integral 3D-equations in direct or inverse modes in the x-direction. The viscous-inviscid coupling is fully 3D and solved the time-consistency problem with an extension of the 'Semi-implicit' method previously suggested in two-dimensions. A 3D inviscid subroutine with TSP approximation is used. Results are obtained for transonic steady flows over wings with shock-induced transonic separation. The method provides results for highly three-dimensional flow separations, such as induced by a three-dimensional through at the wall of a flat plate. The 3D viscous-inviscid coupling and the 3D model of the velocity field provide three-dimensional instantaneous skin-friction lines whose patterns exhibit the same complex topology as Navier-Stokes solvers, with foci, nodes, and saddle-points.

Leballeur, J. C.↗