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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Computational Study of Low-Speed Boundary-Layer Laminar-Turbulent Transition over a 30º Swept Plate in a Channel

Reliable boundary-layer laminar-turbulent transition predictions are critical for the design of aircraft that use natural or controlled laminar flow to achieve low drag and better fuel efficiency. During cruise flight through low turbulence, the three-dimensional boundary layer over swept surfaces can transition prematurely due to stationary crossflow instabilities induced by surface roughness. Parabolized stability equations analyses have been successfully used to explain the different phases of the transition process under well controlled perturbations such as steps or discrete roughness elements. However, methodologies are still lacking to model the boundary layer receptivity to distributed random roughness, and the following non-linear instability stages leading to transition. For engineering applications, several semi-empirical methods based on linear amplification factors or local helicity have found growing adoption, but their reliability is limited by their level of simplification and the sparsity of data they are based on. To address this gap, a new study of the effect of distributed random surface roughness on crossflow transition was started, using a 30◦ swept plate as a test vehicle. The plate, which is mounted in the 2 Foot × 3 Foot Low-Speed Boundary-Layer Channel at the NASA Langley Research Center, has been tested with micro-scale roughness of varied geometry near the leading edge. In order to analyze the flow in the test section and prepare Direct Numerical Simulations with resolved roughness, Reynolds-Averaged Navier-Stokes (RANS) simulations of the flow over the plate installed in the wind tunnel test section were performed using the Langtry-Menter transition model with stationary crossflow extension. This paper presents the results of these simulations, with comparison to the experiment.

ARMD↗
NASA NTRS2016

Comparison of NACA 0012 Laminar Flow Solutions: Structured and Unstructured Grid Methods

In this paper we consider the solution of the compressible Navier-Stokes equations for a class of laminar airfoil flows. The principal objective of this paper is to demonstrate that members of this class of laminar flows have steady-state solutions. These laminar airfoil flow cases are often used to evaluate accuracy, stability and convergence of numerical solution algorithms for the Navier-Stokes equations. In recent years, such flows have also been used as test cases for high-order numerical schemes. While generally consistent steady-state solutions have been obtained for these flows using higher order schemes, a number of results have been published with various solutions, including unsteady ones. We demonstrate with two different numerical methods and a range of meshes with a maximum density that exceeds 8 × 106 grid points that steady-state solutions are obtained. Furthermore, numerical evidence is presented that even when solving the equations with an unsteady algorithm, one obtains steady-state solutions.

Swanson, R. C.↗
NASA NTRS1988

A block space-marching method for the Navier-Stokes equations

A solution procedure is presented for predominantly supersonic viscous flows. The procedure approximately solves the Navier-Stokes equations by marching blocks of grid points in the streamwise direction and solving the fully elliptic equations within each block. In this manner elliptic effects of limited streamwise extent may be accurately calculated. Results are presented for calculations of a Mach 2 laminar flat-plate shock/boundary-layer interaction, and for a Mach 10 hypervelocity interceptor.

Imlay, Scott T.↗
NASA NTRS1977

Computation of viscous transonic flow about a lifting airfoil

The Reynolds averaged Navier-Stokes equations are solved numerically for the viscous transonic flow about a stationary NACA 64A010 airfoil in free air. This paper presents descriptions of the numerical method, turbulence models employed, and boundary conditions appropriate to simulation of free-air flight. Computed results are presented for the airfoil at a free-stream Mach number of 0.8, angles of attack of 0 and 2 deg, and a Reynolds number based on a chord of 4 x 10 to the 6th. For the lifting case, unsteady periodic motion was calculated along the aft portion of the airfoil and in its wake. Recent experimental results obtained by Johnson indicate periodicity aft of the shock closely approximates the computed frequency, but the amplitude of the disturbances was significantly less than the calculated amplitude.

Walitt, L.↗
NASA NTRS1988

Investigation of two-dimensional chemically reacting and radiating supersonic channel flows

The two-dimensional time-dependent Navier-Stokes equations are solved in order to study supersonic flows with finite rate chemistry and radiation for hydrogen-air systems. The problem of the flow in a channel with a ten-degree compression-expansion ramp is solved using the finite volume technique of Jameson et al. (1981) and the unsplit finite difference scheme of MacCormack (1969). The problem of chemically reacting and radiating flows is considered for the flow of premixed hydrogen-air through a channel with parallel boundaries and a channel with a compression corner. Results suggest that radiative interaction can have a significant effect on the entire flowfield.

Mani, M.↗
NASA NTRS1992

Hypersonic code efficiency and validation studies

Renewed interest in hypersonic and supersonic flows spurred the development of the Compressible Navier-Stokes (CNS) code. Originally developed for external flows, CNS was modified to enable it to also be applied to internal high speed flows. In the initial phase of this study CNS was applied to both internal flow applications and fellow researchers were taught to run CNS. The second phase of this research was the development of surface grids over various aircraft configurations for the High Speed Research Program (HSRP). The complex nature of these configurations required the development of improved surface grid generation techniques. A significant portion of the grid generation effort was devoted to testing and recommending modifications to early versions of the S3D surface grid generation code.

Bennett, Bradford C.↗
NASA NTRS1977

Reynolds number effects on shock-wave turbulent boundary-layer interactions - A comparison of numerical and experimental results

An experiment is described that tests and guides computations of a shock-wave turbulent boundary-layer interaction flow over a 20-deg compression corner at Mach 2.85. Numerical solutions of the time-averaged Navier-Stokes equations for the entire flow field, employing various turbulence models, are compared with the data. Each model is critically evaluated by comparisons with the details of the experimental data. Experimental results for the extent of upstream pressure influence and separation location are compared with numerical predictions for a wide range of Reynolds numbers and shock-wave strengths.

Horstman, C. C.↗
NASA NTRS1993

A study of the accuracy of neutrally buoyant bubbles used as flow tracers in air

Research has been performed to determine the accuracy of neutrally buoyant and near neutrally buoyant bubbles used as flow tracers in air. Theoretical, computational, and experimental results are presented to evaluate the dynamics of bubble trajectories and factors affecting their ability to trace flow-field streamlines. The equation of motion for a single bubble was obtained and evaluated using a computational scheme to determine the factors which affect a bubble's trajectory. A two-dimensional experiment was also conducted to experimentally determine bubble trajectories in the stagnation region of NACA 0012 airfoil at 0 deg angle of attack using a commercially available helium bubble generation system. Physical properties of the experimental bubble trajectories were estimated using the computational scheme. These properties included the density ratio and diameter of the individual bubbles. the helium bubble system was then used to visualize and document the flow field about a 30 deg swept semispan wing with simulated glaze ice. Results were compared to Navier-Stokes calculations and surface oil flow visualization. The theoretical and computational analysis have shown that neutrally buoyant bubbles will trace even the most complex flow patterns. Experimental analysis revealed that the use of bubbles to trace flow patterns should be limited to qualitative measurements unless care is taken to ensure neutral buoyancy. This is due to the difficulty in the production of neutrally buoyant bubbles.

Kerho, Michael F.↗
NASA NTRS1986

Separated flows near the nose of a body of revolution

The solution of the Navier-Stokes equations for the problem of cross-flow separataion about a deforming cylinder was achieved by iteration. It was shown that the separation starts at the rear stagnation point and the point of primary separation moves upstram along the cylinder surface. A general method of linear stability analysis for nonparallel external flows was constructed, which consists of representing the eigenfunctions with complete orthogonal sets and forms characteristic equations with the Galerkin method. The method was applied to the Kovasznay flow which is an exact solution of the Navier-Stokes equation. The results show that when the critical parameter is exceeded, there are only a few isolated unstable eigen-frequencies. Another exact solution is shown to be absolutely and monotonically stable with respect to infinitesimal disturbances of all frequencies. The flow is also globally, asymptotically, and monotonically stable in the mean with respect o three-dimensional disturbances. This result forms the sound foundation of rigorous stability analysis for nonparallel flows, and provides an invaluable test ground for future studies of nonparallel flows in which the basic states do not posses exact solutions. The application of this method to the study of the formation of spiral vorticies near the nose of a rotating body of revolution is underway. The same method will be applied to the stability analysis of reversed flow over a plate with suction.

Lin, S. P.↗
NASA NTRS1976

Numerical simulation of supersonic and hypersonic turbulent compression corner flows using relaxation models

Relaxation turbulence eddy viscosity models are incorporated to solve the complete Navier-Stokes equations for supersonic and hypersonic flows over a two-dimensional compression corner. The system of equations is solved by a time-split, second-order accurate numerical scheme. Details of relaxation process are studied, and several relaxation models are tested and compared. Good improvement in the prediction of upstream pressure propagation is obtained for a Mach number of 2.96, and Reynolds number of 10,000,000, with wedge angle of 25 deg. However, the application of the relaxation models to hypersonic flow at Mach number 8.66 and Reynolds number of 22,000,000, with highly cooled wall, shows unfavorable effects on heat transfer and skin friction in the reattachment and recompression regions.

Hung, C. M.↗
NASA NTRS1987

Numerical solutions of Navier-Stokes equations for a Butler wing

The flow field is simulated on the surface of a Butler wing in a uniform stream. Results are presented for Mach number 3.5 and Reynolds number of 2,000,000. The simulation is done by integrating the viscous Navier-Stokes equations. These equations govern the unsteady, viscous, compressible and heat conducting flow of an ideal gas. The equations are written in curvilinear coordinates so that the wing surface is represented accurately. The O-type and H-type grids have been used for this study, and results are compared. The governing equations are solved by the MacCormack time-split method, and the results are compared with other theoretical and experimental results. The codes are written in FORTRAN, vectorized and currently run on the CDC Vector Processing System (VPS-32) computer.

Abolhassani, Jamshid S.↗
NASA NTRS1973

Finite element analysis and computer graphics visualization of flow around pitching and plunging airfoils

A general computational method for analyzing unsteady flow around pitching and plunging airfoils was developed. The finite element method was applied in developing an efficient numerical procedure for the solution of equations describing the flow around airfoils. The numerical results were employed in conjunction with computer graphics techniques to produce visualization of the flow. The investigation involved mathematical model studies of flow in two phases: (1) analysis of a potential flow formulation and (2) analysis of an incompressible, unsteady, viscous flow from Navier-Stokes equations.

Bratanow, T.↗
NASA NTRS1975

Numerical solution of the Navier-Stokes equations for arbitrary two-dimensional airfoils

A method for numerical solution of the Navier-Stokes equations for the flow about arbitrary airfoils or other bodies is presented. This method utilizes a numerically generated curvilinear coordinate system having a coordinate line coincident with the body contour. Streamlines, velocity profiles, and pressure and force coefficients for several airfoils and an arbitrary rock are given. Potential flow solutions are also presented. The procedure capable of treating multiple-element airfoils, and potential flow results are presented.

Thames, F. C.↗
NASA NTRS1975

A general algorithm using finite element method for aerodynamic configurations at low speeds

A finite element algorithm for numerical simulation of two-dimensional, incompressible, viscous flows was developed. The Navier-Stokes equations are suitably modelled to facilitate direct solution for the essential flow parameters. A leap-frog time differencing and Galerkin minimization of these model equations yields the finite element algorithm. The finite elements are triangular with bicubic shape functions approximating the solution space. The finite element matrices are unsymmetrically banded to facilitate savings in storage. An unsymmetric L-U decomposition is performed on the finite element matrices to obtain the solution for the boundary value problem.

Balasubramanian, R.↗
NASA NTRS1976

The relative merits of several numerical techniques for solving the compressible Navier-Stokes equations

Four explicit finite difference techniques designed to solve the time-dependent, compressible Navier Stokes equations are compared. These techniques are: (1) MacCormack, (2) modified Du Fort-Frankel, (3) modified hopscotch, and (4) Brailovskaya. The comparison was made numerically by solving the quasi-one dimensional Navier Stokes equations for the flow in a converging-diverging nozzle. Solutions with and without standing normal shock waves were computed for unit Reynolds numbers (based on total conditions) ranging from 45374 to 2269. The results indicate that all four techniques are comparable in accuracy; however, the modified hopscotch scheme is two to three times faster than the Brailovskaya and MacCormack schemes and three to six times faster than the modified Du Fort-Frankel scheme.

Holst, T. L.↗
NASA NTRS1978

Coefficient matrices for implicit finite difference solution of the inviscid fluid conservation law equations

Although the Navier-Stokes equations describe most flows of interest in aerodynamics, the inviscid conservation law equations may be used for small regions with viscous forces. Thus, Euler equations and several time-accurate finite difference procedures, explicit and implicit, are discussed. Although implicit techniques require more computational work, they permit larger time steps to be taken without instability. It is noted that the Jacobian matrices for Euler equations in conservation-law form have certain eigenvalue-eigenvector properties which may be used to construct conservative-form coefficient matrices. This reduces the computation time of several implicit and semiimplicit schemes. Extensions of the basic approach to other areas are suggested.

Steger, J. L.↗
NASA NTRS1979

Program of research in severe storms

Two modeling areas, the development of a mesoscale chemistry-meteorology interaction model, and the development of a combined urban chemical kinetics-transport model are examined. The problems associated with developing a three dimensional combined meteorological-chemical kinetics computer program package are defined. A similar three dimensional hydrostatic real time model which solves the fundamental Navier-Stokes equations for nonviscous flow is described. An urban air quality simulation model, developed to predict the temporal and spatial distribution of reactive and nonreactive gases in and around an urban area and to support a remote sensor evaluation program is reported.

Source record↗