SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1981

Modifying Stokes' function to reduce the error of geoid undulation computations

Errors in computations of the geoid undulation, which determines the position of the geoid with respect to a reference ellipsoid, introduced by various modifications of the Stokes' function in the Stokes integral formula for the undulation are examined in light of currently available harmonic coefficients. The variants considered include the conventional Molodenskii truncation theory, in which integration is limited to a spherical cap, and the variations suggested by Molodenskii (1958) and Meissl (1971). In computations of the undulation from gravity anomalies in a spherical cap and terrestrial and satellite harmonic potential coefficients (GEM 9 and GEM 10B), it is shown that either of the modifications to the conventional Molodenskii theory results in substantial reductions in truncation errors, with the simpler modification of Meissl even providing superior results in certain circumstances.

Jekeli, C.↗

NASA NTRS1981

The Stokes problem for the ellipsoid using ellipsoidal kernels

A brief review of Stokes' problem for the ellipsoid as a reference surface is given. Another solution of the problem using an ellipsoidal kernel, which represents an iterative form of Stokes' integral, is suggested with a relative error of the order of the flattening. On studying of Rapp's method in detail the procedures of improving its convergence are discussed.

Zhu, Z.↗

NASA NTRS1981

Turbulent solution of the Navier-Stokes equations for uniform shear flow

To study the nonlinear physics of uniform turbulent shear flow, the unaveraged Navier-Stokes equations are solved numerically. This extends our previous work in which mean gradients were absent. For initial conditions, modified three-dimensional-cosine velocity fluctuations are used. The boundary conditions are modified periodic conditions on a stationary three-dimensional numerical grid. A uniform mean shear is superimposed on the initial and boundary conditions. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. As in the case of no shear the initially nonrandom flow develops into an apparently random turbulence at higher Reynolds number. Thus, randomness or turbulence can apparently arise as a consequence of the structure of the Navier-Stokes equations. Except for an initial period of adjustment, all fluctuating components grow with time. The initial equality of the three intensity components is destroyed by the shear, the transverse components becoming smaller than the longitudinal one, in agreement with experiment. Also, the shear creates a small-scale structure in the turbulence. The nonlinear solutions are compared with linearized ones.

Deissler, R. G.↗

NASA NTRS1983

Parabolized Navier-Stokes solutions of separation and trailing-edge flows

A robust, iterative solution procedure is presented for the parabolized Navier-Stokes or higher order boundary layer equations as applied to subsonic viscous-inviscid interaction flows. The robustness of the present procedure is due, in part, to an improved algorithmic formulation. The present formulation is based on a reinterpretation of stability requirements for this class of algorithms and requires only second order accurate backward or central differences for all streamwise derivatives. Upstream influence is provided for through the algorithmic formulation and iterative sweeps in x. The primary contribution to robustness, however, is the boundary condition treatment, which imposes global constraints to control the convergence path. Discussed are successful calculations of subsonic, strong viscous-inviscid interactions, including separation. These results are consistent with Navier-Stokes solutions and triple deck theory.

Brown, J. L.↗

NASA NTRS1981

Anti-Stokes scattering as an XUV radiation source

Harris (1977) has shown that laser induced scattering from atoms stored in a metastable level may be used to produce an XUV radiation source with several unique properties, including narrow linewidth, tunability, linear polarization, picosecond pulsewidth, and relatively high peak spectral brightness. The present investigation provides a summary concerning the progress made on the use of this radiation source as a unique instrument for XUV spectroscopy, and as a possible flashlamp for a 200 A laser. Zych et al. (1978) demonstrated that the maximum spectral brightness of the source is that of a blackbody with a temperature of the metastable storage level. The peak count rate and peak spectral brightness of the anti-Stokes emission at 569 A are 17 times and 139 times greater than those of the 584 A resonance line. The anti-Stokes process was employed as an external source of radiation for high resolution spectroscopy of the 3p6 shell of potassium. Attention is also given to emission spectroscopy of lithium, metastable levels which are appropriate for storage, and a hollow cathode discharge.

Harris, S. E.↗

NASA NTRS1983

Calculation of a simulated 3-D high speed inlet using the Navier-Stokes equations

A hybrid numerical algorithm, developed to solve the full three-dimensional Navier-Stokes equations, is applied to the computation of the flowfield in a simulated three-dimensional high speed aircraft inlet at a Mach number of 2.5 and Reynolds number of 1.4 x 10 to the 7th based on inlet length. The numerical algorithm incorporates a coordinate transformation in order to handle general flow geometries, and utilizes the algebraic turbulent eddy viscosity model of Baldwin and Lomax. The hybrid algorithm has been vectorized on the CDC CYBER 203 computer using the SL/1 vector programming language developed at NASA Langley. The computed results are compared with experimental measurements of the ramp and cowl static pressures, and boundary layer pitot profiles. The results are also compared with a previous two-dimensional Navier-Stokes computation of the same configuration. The agreement with the experimental data is generally good; however, additional improvements in turbulence modeling are needed.

Knight, D. D.↗

NASA NTRS1983

A fully vectorized numerical solution of the incompressible Navier-Stokes equations

A vectorizable algorithm is presented for the implicit finite difference solution of the incompressible Navier-Stokes equations in general curvilinear coordinates. The unsteady Reynolds averaged Navier-Stokes equations solved are in two dimension and non-conservative primitive variable form. A two-layer algebraic eddy viscosity turbulence model is used to incorporate the effects of turbulence. Two momentum equations and a Poisson pressure equation, which is obtained by taking the divergence of the momentum equations and satisfying the continuity equation, are solved simultaneously at each time step. An elliptic grid generation approach is used to generate a boundary conforming coordinate system about an airfoil. The governing equations are expressed in terms of the curvilinear coordinates and are solved on a uniform rectangular computational domain. A checkerboard SOR, which can effectively utilize the computer architectural concept of vector processing, is used for iterative solution of the governing equations.

Patel, N.↗

NASA NTRS1984

Application of the implicit MacCormack scheme to the parabolized Navier-Stokes equations

MacCormack's implicit finite-difference scheme was used to solve the two-dimensional parabolized Navier-Stokes (PNS) equations. This method for solving the PNS equations does not require the inversion of block tridiagonal systems of algebraic equations and permits the original explicit MacCormack scheme to be employed in those regions where implicit treatment is not needed. The advantages and disadvantages of the present adaptation are discussed in relation to those of the conventional Beam-Warming scheme for a flat plate boundary layer test case. Comparisons are made for accuracy, stability, computer time, computer storage, and ease of implementation. The present method was also applied to a second test case of hypersonic laminar flow over a 15% compression corner. The computed results compare favorably with experiment and a numerical solution of the complete Navier-Stokes equations.

Lawrence, J. L.↗

NASA NTRS1983

Accuracy of approximations to the Navier-Stokes equations

The systems of truncated differential equations that have been proposed to reduce the complexity and large computational costs of solutions to the full Navier-Stokes equations are considered. These systems are computationally efficient and capture all the physically relevant behavior. The systems follow a certain hierarchy: (1) the classical boundary-layer equations with specified edge properties (usually the streamwise pressure distribution); (2) the coupled boundary-layer/inviscid equations; (3) the so-called thin-layer equations that discard streamwise diffusion; and (4) the Navier-Stokes equations. Consideration is given to each of these approximations applied to an incompressible, laminar-separating flow at low and moderate Reynolds numbers. It is pointed out that for any flow or region of flow for which viscous-inviscid interaction effects are small, classical boundary-layer equations will provide a satisfactory description of the viscous flow at a fraction of the computational cost of any higher approximations.

Murphy, J. D.↗

NASA NTRS1984

Navier-Stokes calculations for unsteady three-dimensional vortical flows in unbounded domains

Finite-difference Navier-Stokes calculations for unsteady, three-dimensional, incompressible, viscous flows induced by initial vorticity distributions are presented and discussed in this paper. The initial vorticity distributions are assumed to be embedded in a flow field of infinite extent that is quiescent at infinity. These vorticity distributions are typical of vortex rings and other closed vortical tubes or structures. Such structures are important elements in fluid flows such as jets, atmospheric convection and the far-field wakes of aircraft; studies of their interaction may aid in an understanding of complex fluid flows. The calculations employ a method recently proposed by Ting to approximate the infinite-domain boundary value problem with a finite boundary computational domain, and this method is shown to yield accurate three-dimensional results for reasonable expenditures of computer time. Because of the efficiency of the boundary condition technique and the resulting Navier-Stokes code, a 16-bit minicomputer with virtual memory was capable of performing the calculations for the unsteady motion of two obliquely colliding vortex rings. The results of these calculations are presented in the paper.

Chamberlain, J. P.↗

NASA NTRS1984

Dynamic response of shock waves in transonic diffuser and supersonic inlet - An analysis with the Navier-Stokes equations and adaptive grid

An existing method which solves the multi-dimensional ensemble-averaged compressible time-dependent Navier-Stokes equations in conjunction with mixing length turbulence model and shock capturing technique has been extended to include the shock-tracking adaptive grid systems. The numerical scheme for solving the governing equations is based on a linearized block implicit approach. The effects of grid-motion and grid-distribution on the calculated flow solutions have been studied in relative detail and this is carried out in the context of physically steady, shocked flows computed with non-stationary grids. Subsequently, the unsteady dynamics of the flows occurring in a supercritically operated transonic diffuser and a mixed compression supersonic inlet have been investigated with the adaptive grid systems by solving the Navier-Stokes equations.

Liu, N.-S.↗

NASA NTRS1984

Three-dimensional Navier-Stokes calculations of multiple interacting vortex rings

Results from a finite-difference Navier-Stokes code for three-dimensional, unsteady, vortical flows in unbounded domains are presented and analyzed in this paper. The vortical flows presented are representative of vortex rings and other closed vortical tubes or structures in fluid mechanics. Such structures are important elements in fluid flows such as jets, atmospheric turbulence, and the far-field wakes of aircraft, and studies of their interaction may aid in an understanding of complex fluid flows. The paper demonstrates that computational methods can be used as a viable alternative or supplement to experimental techniques for studying the physics of vortex flows. The separate visualization of vortex stretching, convection, and diffusion is presented in this paper for a single elliptical vortex ring.The calculations employ a truncated series expansion technique to simulate the unbounded nature of the fluid flow with a finite computational domain, which is a more accurate technique than the conventional freestream boundary specification. The numerical divergence of the three-dimensional vorticity field is considered as a useful estimate of truncation error, and the use of a kinetic energy decay law as a calculation check is demonstrated. Results from the Navier-Stokes code are presented for the unsteady motion of two and four vortex rings along parallel axes, and the results agree qualitatively with experimental flow visualization.

Chamberlain, J. P.↗

NASA NTRS1982

Comparison of solutions of the Navier-Stokes and interacting boundary layer equations for separated turbulent flow

Comparisons are presented between numerical solutions of the Navier-Stokes and interacting boundary-layer equations for the subsonic separated turbulent flow over an axisymmetric boattail configuration for which experimental data was available. Comparisons were made using a standard algebraic turbulence model and its modification using a streamwise relaxation model. Grid size studies were performed to provide an approximate assessment of the sensitivity of each method to mesh size. In general, good agreement was observed between the solutions of the Navier-Stokes and interacting boundary-layer equations; however, the interaction approach responded more strongly to the modified turbulence model thereby giving better agreement with experimental data.

Vatsa, V. N.↗

NASA NTRS1984

A note on the variational method of Stokes and DeMarcus for radiative transfer in planetary atmospheres

A generalization of the variational principle of Stokes and DeMarcus (1971) by the construction of two bifunctionals is proposed. The equation of radiative transfer is obtained by imposing the stationary condition on such functionals, and the reflectivity and transmissivity at arbitrary angle are given by their extremum values. Further simplification is possible when trail functions are used that are linear in the variational parameters, in which case the source function derived from the method of Stokes and DeMarcus can give the reflectivity at arbitrary emergent angle to second order. The usefulness and limitations of these variational priinciples are explored with extensive numerical results.

Chow, K. S. K.↗

NASA NTRS1985

Hydrodynamic stability problem formulated by numerical solutions of the Navier-Stokes equations

The problem of hydrodynamic stability and the transition from laminar to turbulent flows are reformulated by seeking numerical solutions of the full, unapproximated Navier-Stokes equations. This method differs significantly from the well known Orr-Sommerfeld equation approach. The oncoming laminar flow is disturbed by forced, time dependent perturbations. The magnitudes of these perturbations are arbitrary. Then, the ensuing spatial and temporal development of the imposed perturbations on the basic flow is calculated by direct numerical solutions of the time dependent Navier-Stokes equations. One of the main advantages of this method is its ability to simulate nonlinear processes. As a specific application of this technique to SSME (Space Shuttle Main Engine) flow configurations, computer programs have been written for the two dimensional flow over a backward facing step. This numerical code will be tested for operational use as part of a continued research collaboration effort with the NASA/MSFC counterparts.

Hyun, J. M.↗

NASA NTRS1985

The San Fernando Observatory video Stokes polarimeter

A study was conducted to determine the suitability of the San Fernando Observatory's 61 cm (24 inch) aperture vacuum solar telescope and 3 m (118 inch) focal length vacuum spectroheliograph for Stokes Polarimetry measurements. The polarization characteristics of these two instruments was measured by determining their Mueller matrices as a function of telescope orientation, field angle, wavelength, grating type, and position of the measuring beam in the telescope entrance window. In general, the polarizing and depolarizing properties are small so that inversion of the system Mueller matrix will permit the accurate measurement of Stokes profiles for vector magnetic field determination. A proposed polarimeter design based on the use of a TV camera system to simultaneously scan six different polarization components of a given line profile is described. This design, which uses no rotating optics or electronic modulators and makes efficient use of the available irradiance, promises to yield high quality vector magnetograms.

Richter, P. H.↗

NASA NTRS1985

Least squares inversion of Stokes profiles in the presence of velocity gradients

The Auer, Heasley and House Stokes inversion procedure in use at High Altitude Observatory is based on the analytic solution of the equation of transfer for polarized light where the representation of the thermodynamic and magnetic structure of the atmosphere is assumed to have a high degree of invariance, namely, a Milne-Eddington (ME) structure with a constant magnetic field. In the presence of invariance breaking gradients the resultant Stokes profiles are represented only approximately, if at all, by analytic forms. The accuracy of the inversion parameters and their significance as measures of actual structure are explored for the ME and the Landman-Finn sunspot models under the effects of velocity gradients. The resulting field parameters are good to a few percent and prove to be insensitive to the errors committed by the use of a ME-representation, but the resulting ME parameters yield a less precise measure of thermal structure.

Skumanich, A.↗

NASA NTRS1985

Time-dependent Navier-Stokes computations of separated flows over airfoils

The application of both a central-difference and an upwind implicit approximate factorization Navier-Stokes algorithm to highly separated flow is described. Using the thin-layer approximation, both algorithms are employed to solve the low Reynolds number laminar flow around a circular cylinder with periodic shedding. Results agree with experimental Strouhal numbers to within about 10 percent. A grid refinement study is conducted using the upwind code for both the thin-layer and complete Navier-Stokes approximations. The circular cylinder solution is shown to be highly dependent upon both grid density outward from the body and grid extent. Finally, the thin-layer central-difference code is used to predict the time-accurate unsteady flow about two stalled airfoils. The NACA 0012 airfoil turbulent solutions are periodic with Strouhal numbers of about 0.1; a strong vortex is periodically shed from the leading edge. The LS(1)-0417 airfoil turbulent solution shows a high frequency oscillatory trailing edge separation. The 20 degree NACA 0012 airfoil solution is massively separated, but shows no evidence of periodicity.

Rumsey, C. L.↗