Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stokes”

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,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

The dynamics of oceanic fronts. I - The Gulf Stream

The establishment and maintenance of the mean hydrographic properties of large-scale density fronts in the upper ocean is considered. The dynamics is studied by posing an initial value problem starting with a near-surface discharge of buoyant water with a prescribed density deficit into an ambient stationary fluid of uniform density; full time dependent diffusion and Navier-Stokes equations are then used with constant eddy diffusion and viscosity coefficients, together with a constant Coriolis parameter. Scaling analysis reveals three independent scales of the problem including the radius of deformation of the inertial length, buoyancy length, and diffusive length scales. The governing equations are then suitably scaled and the resulting normalized equations are shown to depend on the Ekman number alone for problems of oceanic interest. It is concluded that the mean Gulf Stream dynamics can be interpreted in terms of a solution of the Navier-Stokes and diffusion equations, with the cross-stream circulation responsible for the maintenance of the front; this mechanism is suggested for the maintenance of the Gulf Stream dynamics.

Kao, T. W.↗
NASA NTRS1980

Computation of supersonic turbulent flows over an inclined ogive-cylinder-flare

A supersonic turbulent flow over an ogive-cylinder-flare has been solved numerically. Initially, the parabolized Navier-Stokes equations are solved for the ogive cylinder back to a location upstream of the shock-wave and boundary-layer interaction. Then, the time-dependent Navier-Stokes equations with a thin-layer approximation are solved for the remaining cylinder-flare portion. Results for a Mach number of 2.9 and a unit Reynolds number of 11.42 x 10 to the 6th/m are obtained for angles of attack alpha = 0, 4, and 8 deg. Good agreement has been found between computed and experimental results of the surface pressure on the ogive-cylinder portion, and for the interaction region at alpha = 0 and 4 deg. The role of circumferential communication in a three-dimensional shock-wave and boundary-layer interaction flow field is discussed.

Hung, C. M.↗
NASA NTRS1981

Ellipsoidal corrections for geoid undulation computations

The computation of accurate geoid undulations is usually done combining potential coefficient information and terrestrial gravity data in a cap surrounding the computation point. In doing this a spherical approximation is made that can cause the errors that were investigated. The equations dealing with ellipsoidal corrections developed by Lelgemann and by Moritz were used to develop a computational procedure considering the ellipsoid as a reference surface. Terms in the resulting expression for the geoid undulation are identified as ellipsoidal correction terms. These equations were developed for the case where the Stokes function is used, and for the case where the modified Stokes function is used. For a cap of 20 deg the correction can reach -33 cm.

Rapp, R. H.↗
NASA NTRS1981

A Quasi-Parabolic technique for computation of three-dimensional viscous flows

A computational technique is presented for obtaining flowfield solutions to a parabolic form of the Navier-Stokes equations. The point of departure is the General Interpolant Method (GIM) which provides a discretization for partial differential equations on arbitrary three-dimensional geometries. The new scheme, termed Quasi-Parabolic, treats the parabolized equations but with 'time-like' terms appended. Addition of these extra terms, which are relaxed by iteration, avoid many of the singularities inherent in classical parabolic Navier-Stokes methods. Streamwise derivatives are approximated by three-point backward differences and the cross plane operators use an alternating forward-backward sweep. A two-step sequence is used to implement the difference scheme in the spatial dimensions and a time-like relaxation converges the quasi-parabolic procedure at each plane. Solutions are presented for flows in two and three dimensions. Inviscid flows are solved for internal and external applications and viscous flows in boundary layers and free shear layers are also computed with the GIM/Quasi-Parabolic scheme.

Spradley, L. W.↗
NASA NTRS1981

High-resolution interferometric spectrophotopolarimetry

Spectrophotopolarimetric capability can be added to a laboratory interferometer-spectrometer by use of a specially designed module described herein. With the instrument so augmented, high-resolution spectra can be obtained of the Stokes parameters of the reference beam and the beams diffusely reflected or transmitted by a sample medium of interest. For any such beam, the exponential Fourier transforms of the two interferograms obtained with a polarizer-analyzer oriented along the 0 deg and the 90 deg directions provide the spectra of I and Q, separately. Within experimental (and numerical) noise, this I spectrum should be the same as the one obtained with the polarizer removed. The remaining Stokes parameters U and V are obtained with a third interferogram recorded with the polarizer along the 45 deg direction. The complete theory of this instrument is described including the detailed analysis of the polarization-interferograms it provides.

Fymat, A. L.↗
NASA NTRS1981

Algol - A further study of the optical polarization

Results of a study of the phase-locked optical polarization of the Algol triple star system are presented. Polarimetry was obtained in the B filter during 75 nights from 1979 to 1980, and results were folded on the AB orbital period for the equatorial Stokes parameters Q and U. As expected for a nearly edge-on system, the polarization variations are predominantly in one direction of the Stokes-vector plane, and are mainly second-harmonic in character. The variable polarization indicates an eigendirection for the AB system at an angle of 137 + or - 4 deg, and a projected major axis of the third star at 133 + or - 3 deg with respect to the AB system. Considerations of the lowest-order theories of polarizing mechanisms, including limb polarization and reflection, indicate that the AB and AB-C orbital planes are almost normal to each other, rather than coplanar. An angular splitting in the Q-U curve between the two quadrature phases 0.25 and 0.75 is also found which can be explained by centrifugal flattening of the secondary's Roche lobe, or reflection from an optically thick, flattened gas cloud near the L1 point.

Kemp, J. C.↗
NASA NTRS1982

A viscous-inviscid interaction algorithm for three-dimensional turbulent subsonic aerodynamic juncture region flow

An order-of magnitude analysis of the subsonic, three-dimensional, steady time-averaged Navier-Stokes equations for semibounded aerodynamic juncture geometries yields the parabolic Navier-Stokes simplification. The numerical solution of the resultant pressure Poisson equation is cast into complementary and particular parts, yielding an iterative interaction algorithm with an exterior three-dimensional potential flow solution. A parabolic transverse momentum equation set is constructed, wherein robust enforcement first-order continuity effects is accomplished using a penalty differential constraint concept within a finite element solution algorithm. A Reynolds stress constitutive equation, with low turbulence Reynolds number wall functions, is employed for closure, using parabolic forms of the two-equation turbulent kinetic energy-dissipation equation system. The algorithm is employed to predict the three-dimensional turbulent flowfield in the juncture region formed by two intersecting parabolic arcs at a free-stream Mach-number of 0.08, and Re/C = 600,000/m.

Baker, A. J.↗
NASA NTRS1982

Numerical solution of Space Shuttle Orbiter flow field

The supersonic, viscous laminar flow around the Space Shuttle Orbiter forebody has been computed with a parabolized Navier-Stokes code using a generalized coordinate transformation. The initial solution for the nose part of the Orbiter geometry was obtained with a three-dimensional time-dependent Navier-Stokes solver. It was necessary to employ a wind axis oriented coordinate system to obtain the initial solution with the time-dependent code. The generalized PNS technique used in this study allows the solution to be marched from the given initial data surface to any desired surface downstream. A grid point clustering scheme was employed to accurately describe the body shape by clustering points at the wing tip and at the wing body juncture. The computed heat transfer coefficients, pressure coefficients, and shock shapes are compared with the available experimental data for 0 degrees and 30 degrees angle of attack.

Venkatapathy, E.↗
NASA NTRS1982

The CMC:3DPNS computer program for prediction of three-dimensional, subsonic, turbulent aerodynamic juncture region flow. Volume 1: Theoretical

An order-of-magnitude analysis of the subsonic three dimensional steady time averaged Navier-Stokes equations, for semibounded aerodynamic juncture geometries, yields the parabolic Navier-Stokes simplification. The numerical solution of the resultant pressure Poisson equation is cast into complementary and particular parts, yielding an iterative interaction algorithm with an exterior three dimensional potential flow solution. A parabolic transverse momentum equation set is constructed, wherein robust enforcement of first order continuity effects is accomplished using a penalty differential constraint concept within a finite element solution algorithm. A Reynolds stress constitutive equation, with low turbulence Reynolds number wall functions, is employed for closure, using parabolic forms of the two-equation turbulent kinetic energy-dissipation equation system. Numerical results document accuracy, convergence, and utility of the developed finite element algorithm, and the CMC:3DPNS computer code applied to an idealized wing-body juncture region. Additional results document accuracy aspects of the algorithm turbulence closure model.

Baker, A. J.↗
NASA NTRS1983

Numerical studies of laminar and turbulent drag reduction, part 2

The flow over wave shaped surfaces is studied using a Navier Stokes solver. Detailed comparisons with theoretical results are presented, including the stability of a laminar flow over wavy surfaces. Drag characteristics of nonplanar surfaces are predicted using the Navier-Stokes solver. The secondary instabilities of wall bounded and free shear flows are also discussed.

Balasubramanian, R.↗
NASA NTRS1982

The indirect effects on the computation of geoid undulations

The indirect effects on the geoid computation due to the second method of Helmert's condensation were studied. when Helmert's anomalies are used in Stokes's equation, there are three types of corrections to the free air geoid. The first correction, the indirect effect on geoid undulation due to the potential change in Helmert's reduction, had a maximum value of 0.51 meters in the test area covering the United States. The second correction, the attraction change effect on geoid undulation, had a maximum value of 9.50 meters when the 10 deg cap was used in Stokes' equation. The last correction, the secondary indirect effect on geoid undulatin, was found negligible in the test area. The corrections were applied to uncorrected free air geoid undulations at 65 Doppler stations in the test area and compared with the Doppler undulations. Based on the assumption that the Doppler coordinate system has a z shift of 4 meters with respect to the geocenter, these comparisons showed that the corrections presented in this study yielded improved values of gravimetric undulations.

Wichiencharoen, C.↗
NASA NTRS1983

Turbulent drag characteristic of small amplitude rigid surface waves

The paper discusses an experimental and theoretical investigation of low speed turbulent flow over transverse rigid wavy surfaces having a wavelength on the order of the boundary layer thickness. In addition the results of theoretical investigations of laminar flow over short wavelength waves and turbulent flow over long wavelength (wavelength much greater than the boundary layer thickness) wavy walls are presented. Surface pressure and net drag measurements obtained for the short wavelength symmetric and asymmetric waves in a turbulent flow agreed with the results obtained with a Navier-Stokes spectral code and indicated no significant drag reduction. Several asymmetric wave configurations did significantly reduce the pressure drag compared to sine wave geometries, but these waves still did not produce a net drag reduction. The Navier-Stokes calculations for laminar flow indicated a drag reduction of up to 17 percent for certain short wavelength waves. Boundary layer calculations for turbulent flow over long wavelength surface waves indicated that net drag reductions on the order of 10 percent may be possible for large-radius waisted bodies.

Lin, J. C.↗
NASA NTRS1983

Calculation of viscous supersonic flows over finned bodies

The parabolized Navier-Stokes (PNS) equations are used to calculate the viscous, supersonic flow fields about a six-finned projectile and a generic four-finned missile at angles of attack. Since current computer speeds and storage preclude a fully three-dimensional calculation using the unsteady, Reynolds-averaged, Navier-Stokes equations, the applicability of the PNS equations to the above flow fields is of considerable interest. Two important aspects of the calculation are grid generation and the type of smoothing used to prevent nonphysical solutions. This paper includes a description of the grid-generation process. Results in the form of density contours and velocity vector plots are presented for the two configurations. The applicability of the PNS equations to the complicated flow fields considered is successfully demonstrated.

Rai, M. M.↗
NASA NTRS1983

Finite-difference simulation of transonic separated flow using a full potential boundary layer interaction approach

A new, fast, direct-inverse, finite-difference boundary-layer code has been developed and coupled with a full-potential transonic airfoil analysis code via new inviscid-viscous interaction algorithms. The resulting code has been used to calculate transonic separated flows. The results are in good agreement with Navier-Stokes calculations and experimental data. Solutions are obtained in considerably less computer time than Navier-Stokes solutions of equal resolution. Because efficient inviscid and viscous algorithms are used, it is expected this code will also compare favorably with other codes of its type as they become available.

Van Dalsem, W. R.↗
NASA NTRS1983

Numerical study of the effect of an embedded surface-heat source on the separation bubble of supersonic flow

A numerical study of the conjugated problem of a separated supersonic flow field and a conductive solid wall with an embedded heat source is presented. Implicit finite-difference schemes were used to solve the two-dimensional time-dependent compressible Navier-Stokes equations and the time-dependent heat-conduction equation for the solid in both general coordinate systems. A detailed comparison between the thin-layer and Navier-Stokes models was made for steady and unsteady supersonic flow and showed insignificant differences. Steady-state and transient cases were computed and the results show that a temperature pulse at the solid-fluid interface can be used to detect the flow direction near the wall in the vicinity of separation without significant distortion of the flow field.

Degani, D.↗
NASA NTRS1983

Iterative PNS method for attached flows with upstream influence

A stable global-iteration procedure is developed by utilizing successive sweeps, from inflow to outflow boundaries, of a parabolized Navier-Stokes code for attached, steady supersonic flow. It is shown that the procedure converges in about ten or fewer iterations, and allows for the upstream influence within the subsonic region of a supersonic boundary layer. An implicit forward difference is employed in the subsonic region to evaluate the pressure-gradient term in the streamwise momentum equation. The pressure-gradient term is normally approximated in standard single-sweep parabolized methods, suppressing the upstream influence. The iterative results obtained demonstrate the validity of the single-sweep method for weak interactions. For hypersonic viscous interaction on a flat plate at a Mach number of 5.8, it is found that the single-sweep method has a small error that vanishes as the Reynolds number of the flow is increased. However, the skin friction has an error of about 10 percent for low Reynolds numbers and for hot wall conditions. The application of this method to laminar two-dimensional flow over weak expansion and compression corners shows that the pressure and skin friction results in the vicinity of the expansion corner compare favorably with a time-dependent Navier-Stokes numerical solution.

Rakich, J. V.↗
NASA NTRS1982

Turbulent flow simulation - A large eddy simulator's viewpoint

The state of the art of the data base on turbulent structures and numerical turbulence modelling is assessed. The failure of statistical methods to produce accurate results revealed that turbulence is more than simply random noise. Both randomness and coherence may be characteristics of turbulence, which may also feature noise that is a remnant of previous coherence. Five methods are presently used for numerically predicting turbulent flow behavior: correlations, integral methods, the k-epsilon model, and phenomenological and exact solutions to the Navier-Stokes equations. The higher level models include large eddy simulation and full simulation with numerical solutions to the Navier-Stokes equations; the two methods are applicable to experimental, rather than design work.

Ferziger, J. H.↗
NASA NTRS1984

Physical aspects of computing the flow of a viscous fluid

One of the main themes in fluid dynamics at present and in the future is going to be computational fluid dynamics with the primary focus on the determination of drag, flow separation, vortex flows, and unsteady flows. A computation of the flow of a viscous fluid requires an understanding and consideration of the physical aspects of the flow. This is done by identifying the flow regimes and the scales of fluid motion, and the sources of vorticity. Discussions of flow regimes deal with conditions of incompressibility, transitional and turbulent flows, Navier-Stokes and non-Navier-Stokes regimes, shock waves, and strain fields. Discussions of the scales of fluid motion consider transitional and turbulent flows, thin- and slender-shear layers, triple- and four-deck regions, viscous-inviscid interactions, shock waves, strain rates, and temporal scales. In addition, the significance and generation of vorticity are discussed. These physical aspects mainly guide computations of the flow of a viscous fluid.

Mehta, U. B.↗