Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stokes flow”

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 325 records · Page 18

Optimal Growth in Hypersonic Boundary Layers

The linear form of the parabolized linear stability equations is used in a variational approach to extend the previous body of results for the optimal, nonmodal disturbance growth in boundary-layer flows. This paper investigates the optimal growth characteristics in the hypersonic Mach number regime without any high-enthalpy effects. The influence of wall cooling is studied, with particular emphasis on the role of the initial disturbance location and the value of the spanwise wave number that leads to the maximum energy growth up to a specified location. Unlike previous predictions that used a basic state obtained from a self-similar solution to the boundary-layer equations, mean flow solutions based on the full Navier-Stokes equations are used in select cases to help account for the viscous- inviscid interaction near the leading edge of the plate and for the weak shock wave emanating from that region. Using the full Navier-Stokes mean flow is shown to result in further reduction with Mach number in the magnitude of optimal growth relative to the predictions based on the self-similar approximation to the base flow.

Navier Stokes equations↗

Analysis of a high speed civil transport configuration at subsonic flow conditions using a Navier-Stokes solver

Computations of three dimensional vortical flows over a generic High Speed Civil Transport (HSCT) configuration with an aspect ratio of 3.04 are performed using a thin-layer Navier-Stokes solver. The HSCT cruise configuration is modeled without leading or trailing edge flap deflections and without engine nacelles. The flow conditions, which correspond to tests done in the NASA Langley 8-Foot Transonic Pressure Tunnel (TPT), are a subsonic Mach number of 0.3 and Reynolds number of 4.4 million for a range-of-attack (-.23 deg to 17.78 deg). The effects of the farfield boundary location with respect to the body are investigated. The boundary layer is assumed turbulent and simulated using an algebraic turbulence model. The key features of the vortices and their interactions are captured. Grid distribution in the vortex regions is critical for predicting the correct induced lift. Computed forces and surface pressures compare reasonably well with the experimental TPT data.

Lessard, Victor R.↗

Navier-Stokes computations of lee-side flows over delta wings

Solutions to the Navier-Stokes equations for the flow over delta wings are computed with emphasis on the separated vortical flows developing on the lee side at high angles of attack. A recently developed implicit algorithm is used which employs upwind differencing for the pressure and convection terms and central differencing for the shear stress and heat transfer terms. Solutions to both the three-dimensional equations and the approximate conical flow equations are compared parametrically with an extensive experimental data base at supersonic speeds. The computations indicate that the conical flow approximation provides results in close agreement with the three-dimensional equations, even to angles of attack as high as 20 degrees. Good agreement with experimentally measured pressures and vapor screen photographs is obtained for the conditions investigated. The method predicts the classical pattern of vortical flow over a delta wing and transition to other flow patterns as the leading edge sweep angle and leading edge normal Mach number are varied.

Thomas, J. L.↗

A three-dimensional Navier-Stokes stage analysis of the flow through a compact radial turbine

A steady, three-dimensional Navier-Stokes average passage computer code is used to analyze the flow through a compact radial turbine stage. The code is based upon the average passage set of equations for turbomachinery, whereby the flow fields for all passages in a given blade row are assumed to be identical while retaining their three-dimensionality. A stage solution is achieved by alternating between stator and rotor calculations, while coupling the two solutions by means of a set of axisymmetric body forces which model the absent blade row. Results from the stage calculation are compared with experimental data and with results from an isolated rotor solution having axisymmetric inlet flow quantities upstream of the vacated stator space. Although the mass-averaged loss through the rotor is comparable for both solutions, the details of the loss distribution differ due to stator effects. The stage calculation predicts smaller spanwise variations in efficiency, in closer agreement with the data. The results of the study indicate that stage analyses hold promise for improved prediction of loss mechanisms in multi-blade row turbomachinery, which could lead to improved designs through the reduction of these losses.

Heidmann, James D.↗

Solutions to Three-Dimensional Thin-Layer Navier-Stokes Equations in Rotating Coordinates for Flow Through Turbomachinery

The viscous, Navier-Stokes solver for turbomachinery applications, MSUTC has been modified to include the rotating frame formulation. The three-dimensional thin-layer Navier-Stokes equations have been cast in a rotating Cartesian frame enabling the freezing of grid motion. This also allows the flow-field associated with an isolated rotor to be viewed as a steady-state problem. Consequently, local time stepping can be used to accelerate convergence. The formulation is validated by running NASA's Rotor 67 as the test case. results are compared between the rotating frame code and the absolute frame code. The use of the rotating frame approach greatly enhances the performance of the code with respect to savings in computing time, without degradation of the solution.

Ghosh, Amrit Raj↗

A streamwise upwind algorithm for the Euler and Navier-Stokes equations applied to transonic flows

A new algorithm was developed for the Euler and Navier-Stokes equations that uses upwind differencing based on the streawise direction. This algorithm is time accurate and can be used in codes for calculating unsteady transonic flows over wings. Such codes can be used for the flutter analysis of wings. In this algorithm, the coordinate system is locally rotated to align with the streamwise direction. For differencing the convective terms in the streamwise direction, a new form of flux splitting is employed, in which the biasing depends on the local Mach number. In the plane perpendicular to the stream direction, the new flux splitting uses the condition of no flow in that local plane. By using a locally rotated coordinate system, the convective flux vector biasing depends on the total Mach number. Hence, the switching of the flux vector biasing occurs across shock waves and the proper domain of dependence is used in supersonic regions. For comparison, many other upwind methods switch differencing based on Mach number of shock waves in multidimensional flows. The formulas for the convective flux vector differencing do not contain any user specified parameters. So, the amount of numerical dissipation is automatically determined.

Goorjian, Peter M.↗

A streamwise upwind algorithm for the Euler and Navier-Stokes equations applied to transonic flows

A new algorithm was developed for the Euler and Navier-Stokes equations that uses upwind differencing based on the streamwise direction. This algorithm is time accurate and can be used in codes for calculating unsteady transonic flows over wings. Such codes can be used for the flutter analysis of wings. In this algorithm, the coordinate system is locally rotated to align with the streamwise direction. For differencing the convective terms in the streamwise direction, a new form of flux splitting is employed, in which the biasing depends on the local Mach number. In the plane perpendicular to the stream direction, the new flux splitting uses the condition of no flow in that local plane. By using a locally rotated coordinate system, the convective flux vector biasing depends on the total Mach number. Hence, the switching of the flux vector biasing occurs across shock waves and the proper domain of dependence is used in supersonic regions. For comparison, many other upwind methods switch differencing based on Mach number of shock wave in multidimensional flows. The formulas for the convective flux vector differencing do not contain any user specified parameters. So, the amount of numerical dissipation is automatically determined.

Goorjian, P. M.↗

A new algorithm for the Navier-Stokes equations applied to transonic flows over wings

A new algorithm has been developed for the Navier-Stokes equations. For the convective terms, this algorithm employs flux vector splitting together with a locally rotated coordinate system that is aligned with the streamwise direction. The algorithm has been implemented into a Navier-Stokes code, the NASA Ames Research Center's TNS code, and several cases of steady flow have been calculated. The results show improvements over the original method in the code. First, in a case of separated flow, there is better agreement with the experimental results and also there is an absence of numerical oscillations that occur with the original method. Second, in a case with a strong supersonic to subsonic shock wave, comparisons show an improvement in the resolution of the shock wave and the reexpansion singularity.

Goorjian, Peter M.↗

An efficient and robust algorithm for two dimensional time dependent incompressible Navier-Stokes equations: High Reynolds number flows

An algorithm is presented for unsteady two-dimensional incompressible Navier-Stokes calculations. This algorithm is based on the fourth order partial differential equation for incompressible fluid flow which uses the streamfunction as the only dependent variable. The algorithm is second order accurate in both time and space. It uses a multigrid solver at each time step. It is extremely efficient with respect to the use of both CPU time and physical memory. It is extremely robust with respect to Reynolds number.

Goodrich, John W.↗

An efficient and robust algorithm for two dimensional time dependent incompressible Navier-Stokes equations - High Reynolds number flows

An algorithm is presented for unsteady two-dimensional incompressible Navier-Stokes calculations. This algorithm is based on the fourth order partial differential equation for incompressible fluid flow which uses the streamfunction as the only dependent variable. The algorithm is second order accurate in both time and space. It uses a multigrid solver at each time step. It is extremely efficient with respect to the use of both CPU time and physical memory. It is extremely robust with respect to Reynolds number.

Goodrich, John W.↗

A compressible Navier-Stokes code for turbulent flow modeling

An implicit, finite volume code for solving two dimensional, compressible turbulent flows is described. Second order upwind differencing of the inviscid terms of the equations is used to enhance stability and accuracy. A diagonal form of the implicit algorithm is used to improve efficiency. Several zero and two equation turbulence models are incorporated to study their impact on overall flow modeling accuracy. Applications to external and internal flows are discussed.

Coakley, T. J.↗

Navier-Stokes computations of vortical flows over low aspect ratio wings

An upwind-biased finite-volume algorithm is applied to the low-speed flow over a low aspect ratio delta wing from zero to forty degrees angle of attack. The differencing is second-order accurate spatially, and a multigrid algorithm is used to promote convergence to the steady state. The results compare well with the detailed experiments of Hummel (1983) and others for a Re(L) of 0.95 x 10 to the 6th. The predicted maximum lift coefficient of 1.10 at thirty-five degrees angle of attack agrees closely with the measured maximum lift of 1.06 at thirty-three degrees. At forty degrees angle of attack, a bubble type of vortex breakdown is evident in the computations, extending from 0.6 of the root chord to just downstream of the trailing edge.

Thomas, J. L.↗

Incompressible Navier-Stokes computations of vortical flows over double-delta wings

An implicit flux-difference splitting scheme is used to compute incompressible laminar vortical flows about a thin round-edged double-delta wing with 80 deg and 60 deg leading-edged sweep for the strake and the main wing respectively. The numerical scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. It is second-order accurate spatially by applying TVD-like upwind discretization to the inviscid fluxes and central differencing to the viscous shear fluxes. The interaction between the two primary vortices emanating from the apex and kink leading-edges is successfully simulated. The computed trajectories of vortical cores compare well with experimental data.

Hsu, Chung-Hao↗

Navier-Stokes calculations of transonic flows past open and transitional cavities

Computational simulations of deep and transitional cavities at transonic regimes were performed. The three-dimensional unsteady separation on the cavity floor was analyzed, and computational flow visualization inside the cavity was done. A parametric study was performed to investigate the effects of varying Mach number and the ratio delta/D, where delta is the instantaneous boundary layer thickness at the upstream lip and D is the depth. Both time-averaged and instantaneous solutions were obtained. Comparison of the solutions with experimental data were generally favorable; some discrepancies are explained.

Srinivasan, S.↗