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 37 records · Page 2

Numerical simulation of the tip vortex off a low-aspect-ratio wing at transonic speed

The viscous transonic flow around a low-aspect-ratio wing has been computed using an implicit, three-dimensional, 'thin-layer' Navier-Stokes solver. The grid around the geometry of interest is obtained numerically as a solution to a Dirichlet problem for the cube. The geometry chosen for this study is a low-aspect-ratio wing with large sweep, twist, taper, and camber. The topology chosen to wrap the mesh around the wing with good tip resolution is a C-O type mesh. Using this grid, the flow around the wing was computed for a free-stream Mach number of 0.82 at an angle of attack of 5 deg. At this Mach number, an oblique shock forms on the upper surface of the wing, and a tip vortex and three-dimensional flow separation off the wing surface are observed. Particle path lines indicate that the three-dimensional flow separation on the wing surface is part of the roots of the tip-vortex formation. The lifting of the tip vortex before the wing trailing edge is clearly observed by following the trajectory of particles released around the wing tip.

Mansour, N. N.↗

Comparison of secondary flows predicted by a viscous code and an inviscid code with experimental data for a turning duct

A comparison of the secondary flows computed by the viscous Kreskovsky-Briley-McDonald code and the inviscid Denton code with benchmark experimental data for turning duct is presented. The viscous code is a fully parabolized space-marching Navier-Stokes solver while the inviscid code is a time-marching Euler solver. The experimental data were collected by Taylor, Whitelaw, and Yianneskis with a laser Doppler velocimeter system in a 90 deg turning duct of square cross-section. The agreement between the viscous and inviscid computations was generally very good for the streamwise primary velocity and the radial secondary velocity, except at the walls, where slip conditions were specified for the inviscid code. The agreement between both the computations and the experimental data was not as close, especially at the 60.0 deg and 77.5 deg angular positions within the duct. This disagreement was attributed to incomplete modeling of the vortex development near the suction surface.

Schwab, J. R.↗

A parabolized Navier-Stokes algorithm for separated supersonic internal flows

A stable global-iteration procedure is applied to a fully implicit noniterative parabolized Navier-Stokes algorithm to improve the solutions of the single-sweep method as well as to resolve small separation bubbles due to shock-wave/boundary-layer interactions. A modified form of the FLARE approximation is employed in the separated regions and some requirements for stable iterations are discussed. The results of the global-iteration scheme for two test cases of supersonic internal flows are compared against the results of the single-sweep method as well as the results of the full Navier-Stokes solver.

Chitsomboon, T.↗

Numerical study of Reynolds stress in compressible flows

A second order closure has been implemented in an implicit Navier-Stokes solver to study the behavior of the Reynolds stresses under the influence of severe pressure gradients. In the boundary layer zone, the strongly sheared character of the mean flow dominates the turbulence generation mechanisms. However, the pressure gradients also play a very important role for these processes, but at different locations within the boundary layer.

Vandromme, D.↗

The behaviour of turbulence anisotropy through shock waves and expansions

A second order closure has been implemented in an implicit Navier-Stokes solver to study the behavior of the Reynolds stresses under the influence of severe pressure gradients. In the boundary layer zone, the strongly sheared character of the mean flow dominates the turbulence generation mechanisms. However, the pressure gradients play also a very important role for these processes, but at different locations within the boundary layer. This aspect may be emphasized by the analysis of turbulence anisotropy through shock waves and expansions.

Minh, H. H.↗

A Navier-Stokes boundary element solver

Using global interpolation functions (GIF's) boundary element solutions are obtained for two-dimensional laminar flows. Two schemes are proposed for handling the convective terms. The first treats convection as a forcing function, and converts the flow equations to pseudo-Poisson equations. In the second scheme, some convective effect is incorporated into the fundamental solution used in constructing the pertinent integral equations. The lid-driven cavity flow is selected as the benchmark problem.

Reddy, D. R.↗

Navier-Stokes computations for exotic airfoils

An efficient hyperbolic grid generator with improvements for handling sharp corners and concave surfaces is combined with an efficient and accurate Navier-Stokes flow solver. This combination is applied to some rather complex two-dimensional airfoil configurations. Steady separated flow about an iced leading edge of an airfoil is presented. Unsteady viscous separated flows past an airfoil at two angles of attack with a spoiler deployed at 60 degrees are presented and compared with experiment. The spoiler computations are performed with two different topological maps of the physical domain to the computational domain. Innovative graphical techniques for both the static and unsteady display of the flow fields are presented and discussed.

Barth, T. J.↗

Navier-Stokes calculations with a coupled strongly implicit method. Part 2: Spline solutions

A coupled strongly implicit method is combined with a deferred-corrector spline solver for the vorticity-stream function form of the Navier-Stokes equation. Solutions for cavity, channel and cylinder flows are obtained with the fourth-order spline 4 procedure. The strongly coupled spline corrector method converges as rapidly as the finite difference calculations and also allows for arbitrary large time increments for the Reynolds numbers considered. In some cases fourth-order smoothing or filtering is required in order to suppress high frequency oscillations.

Rubin, S. G.↗

Numerical approach for the aerodynamic analysis if airfoils with laminar separation

A numerical method for simultaneously and efficiently coupling an external subsonic potential flow and an interior viscous flow such that the two flows match at an interfacing boundary is discussed. Both a panel method and a simple point compressible vortex model are used for the outer potential field. The interior flow solvers which were used are the Navier-Stokes and Euler codes of T. J. Coakley and the Euler code of A. Verhoff. In order to test compatibility, the panel method is coupled to the less expensive Euler codes since the coupling procedure is identical with the Navier-Stokes code. The results show significant efficiency improvements can be obtained over the uncoupled approach. Results also indicate the outer potential flow is best represented by the simple point compressible vortex model. The panel method couples smoothly to Coakley's implicit code but is numerically incompatible as coupled with the explicit Euler code. An improved Navier-Stokes code is under initial development which extends the Euler code to include the necessary viscous terms. Results are shown for all infinite length channel with one wavy periodic wall with and without laminar separation.

Halt, D. W.↗

Navier-Stokes calculations with a coupled strongly implicit method. II Spline deferred-corrector solutions

The coupled strongly implicit (CSIP) method described previously is combined with a deferred-corrector spline solver for the vorticity-stream function form of the Navier-Stokes equations. Solutions for cavity, channel and cylinder flows are obtained with the fourth-order spline 4 procedure. The strongly coupled spline corrector method converges as rapidly as the finite difference calculations and also allows for arbitrary large time increments for the Reynolds numbers considered (equal to or less than 1000). In some cases fourth-order smoothing or filtering is required in order to suppress high frequency oscillations.

Rubin, S. G.↗

Calculation of Unsteady Transonic Flow over an Airfoil

An implicit finite-difference solver for either the Euler equations or the "thin-layer" Navier-Stokes equations was used to calculate a transonic flow over the NACA 64A010 airfoil pitching about its one-quarter chord. An unsteady automatic grid-generation procedure that will improve significantly the computational efficiency of various unsteady flow problems is described. The calculated results for both inviscid and viscous flows at Much number 0.8 over the airfoil oscillating with reduced frequency referenced to one-half chord, 0.2, are compared with experimental data measured in the Ames 11 x 11-ft Transonic Wind Tunnel. Nonlinear, unsteady effects of the flow on the surface pressure variations, shock-wave excursions, and overall airloads are examined. Good agreements between the results of computations and experiments were obtained. In the shock-wave region, however, the results of the viscous-flow computations showed closer agreement with the experimental data.

Chyu, W. J.↗

Calculation of Unsteady Transonic Flow Over an Airfoil

An implicit finite-difference solver for either the Euler equations or the "thin-layer" Navier-Stokes equations was used to calculate a transonic flow over the NACA 64A010 airfoil pitching about its one-quarter chord. An unsteady automatic grid-generation procedure that will improve significantly the computational efficiency of various unsteady flow problems is described. The calculated results for both inviscid and viscous flows at Mach number 0.8 over the airfoil oscillating with reduced frequency referenced to one-half chord, 0.2, are compared with experimental data measured in the Ames 11 x 11 ft Transonic Wind Tunnel. Nonlinear, unsteady effects of the flow on the surface pressure variations, shock-wave excursions, and overall airloads are examined. Good agreements between the results of computations and experiments were obtained. In the shock-wave region, however, the results of the viscous-flow computations showed closer agreement with the experimental data.

Chyu, W. J.↗

Approximate factorization with an elliptic pressure solver for incompressible flow

Two-dimensional curvilinear coordinates are used to solve the incompressible Navier-Stokes equations, in conjunction with approximate factorization for the solution of the momentum equation and the successive overrelaxation by lines method for the solution of a Poisson equation for the pressure. The combined algorithm, although not fully explicit, is marginally stable at Reynolds numbers lower than 10,000 and time increments of 0.01. Pressure distributions calculated for attack angles of zero and 6 deg are of the same shape as the experimental curves, but are shifted to one side.

Bernard, R. S.↗

A rapid solver for hyperbolic systems of equations

A numerical method is proposed for solving the time-dependent compressible Navier-Stokes equations in two dimensions. The equations are time-split into a hyperbolic part and a parabolic part. This paper describes an explicit numerical method for solving the hyperbolic (inviscid) part. The hyperbolic operator is explicit, conservative, uses characteristic relations to predict convection and pressure fields, and is stable under a specified condition.

Maccormack, R. W.↗

The NASA Turbulent Heat Flux Experiments: Summary and Lessons Learned

The Turbulent Heat Flux (THX) experiments were conducted at NASA Glenn Research Center (GRC) in order to collect measurements of velocities and temperatures for computational fluid dynamics (CFD) validation of heated flows, with a focus on propulsion system components. The experiments spanned 5 phases; four of which were conducted in the GRC AeroAcoustic Propulsion Laboratory (AAPL) using the Small Hot Jet Flow Rig (SHJAR). In addition to making velocity measurements with Particle Image Velocimetry (PIV), the THX experiments introduced a new Raman-scattering based capability to measure temperatures. Computational studies were also conducted for each of the experimental configurations, in order to provide a baseline of expected CFD results and conduct an assessment of the capability of various CFD approaches for calculating flows where the turbulent transport of heat was important. Two of the collected sets of data were used for American Institute of Aeronautics and Astronautics (AIAA) Propulsion Aerodynamic Workshops (PAWs). The data set from the 5th phase, collected for heated supersonic jets, was also used to construct new validation cases for the NASA Turbulence Model Resource (TMR). This paper provides an overview of the experiments and associated computations for each of the 5 test phases. Key experimental findings are presented. Lessons learned are provided concerning the effect of computational modeling choice on accuracy of predicting turbulent flows where thermal transport is important. Emphasis is placed on comparing Reynolds- averaged Navier-Stokes approaches with large-eddy simulation approaches. The benefits of utilizing a conjugate heat transfer method in conjunction with CFD solver for film cooling is demonstrated.

RANS↗

The NASA Turbulent Heat Flux (THX) Experiments: Summary and Lessons Learned

The Turbulent Heat Flux (THX) experiments were conducted at NASA Glenn Research Center (GRC) in order to collect measurements of velocities and temperatures for computational fluid dynamics (CFD) validation of heated flows, with a focus on propulsion system components. The experiments spanned 5 phases; four of which were conducted in the GRC AeroAcoustic Propulsion Laboratory (AAPL) using the Small Hot Jet Flow Rig (SHJAR). In addition to making velocity measurements with Particle Image Velocimetry (PIV), the THX experiments introduced a new Raman-scattering based capability to measure temperatures. Computational studies were also conducted for each of the experimental configurations, in order to provide a baseline of expected CFD results and conduct an assessment of the capability of various CFD approaches for calculating flows where the turbulent transport of heat was important. Two of the collected sets of data were used for American Institute of Aeronautics and Astronautics (AIAA) Propulsion Aerodynamic Workshops (PAWs). The data set from the 5th phase, collected for heated supersonic jets, was also used to construct new validation cases for the NASA Turbulence Model Resource (TMR). This paper provides an overview of the experiments and associated computations for each of the 5 test phases. Key experimental findings are presented. Lessons learned are provided concerning the effect of computational modeling choice on accuracy of predicting turbulent flows where thermal transport is important. Emphasis is placed on comparing Reynolds-averaged Navier-Stokes approaches with large-eddy simulation approaches. The benefits of utilizing a conjugate heat transfer method in conjunction with CFD solver for film cooling is demonstrated.

RANS↗

Semidirect calculation of steady two- and three-dimensional flows

This paper describes a semidirect method for rapidly solving steady-state viscous flows described by the complete Navier-Stokes equations. The current results are for two-dimensional incompressible flows in general channels at arbitrary Reynolds numbers, but work in progress on compressible and three-dimensional flows is also described. The basic concept of semidirect methods is to use the recently developed fast (direct, or noniterative) linear solvers to solve linearized equations, which are then iterated to solve the nonlinearity. The method used here is an extension of the Split NOS method (Roache, 1975)

Roache, P. J.↗

Recent computation of viscous effects in transonic flow

The paper describes some improvements to a method of solving the time-dependent Reynolds-averaged Navier-Stokes equations for two-dimensional compressible flows for describing the coupling between viscous and inviscid regions in transonic flows. Mesh patching is implemented at the sublayer-boundary layer interfaces and the boundary layer-inviscid flow interface. Mesh stretching is geometric and is small in regions of large gradients and large in regions of small gradients, thus preserving second-order accuracy. Mesh lines are arranged so as to accommodate MacCormack's rapid solver for hyperbolic systems and a wider class of turbulent transport models. Results of computations of transonic flows over lifting airfoils are compared with experimental data and other solutions.

Deiwert, G. S.↗