Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INCOMPRESSIBLE FLUID”

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 163 records · Page 9

In situ measurements of the plasma bulk velocity near the Io flux tube

The flow around the Io flux tube was studied by analyzing the eleven spectra taken by the Voyager 1 Plasma Science (PLS) experiment in its vicinity. The bulk plasma parameters were determined using a procedure that uses the full response function of the instrument and the data in all four PLS sensors. The mass density of the plasma in the vicinity of Io is found to be 22,500 + or - 2,500 amu/cu cm and its electron density is found to be 1500 + or - 200/cu cm. The Alfven speed was determined using three independent methods; the values obtained are consistent and taken together yield V sub A = 300 + or - 50 km/sec, corresponding to an Alfven Mach number of 0.19 + or - 0.02. For the flow pattern, good agreement was found with the model of Neubauer (1980), and it was concluded that the plasma flows around the flux tube with a pattern similar to the flow of an incompressible fluid around a long cylinder obstacle of radius 1.26 + or - 0.1 R sub Io.

Barnett, A.↗

Calculations of the stability of some axisymmetric flows proposed as a model of vortex breakdown

The term vortex breakdown refers to the abrupt and drastic changes of structure that can sometimes occur in swirling flows. It was conjectured that the bubble type of breakdown can be viewed as an axisymmetric wave traveling upstream in a primarily columnar vortex flow. In this scenario the wave's upstream progress is impeded only when it reaches a critical amplitude and it loses stability to some nonaxisymmetric disturbance. The stability of some axisymmetric wavy flows to three dimensional disturbances, viewing the amplitude of the wave as a bifurcation parameter is examined. The stability of a set of related columnar vortex flows, constructed by taking the two dimensional flow at a single axial location and extending it throughout the domain without variation, is investigated. The method used will be to expand the perturbation velocity in a series of divergence free vectors which ensures that the continuity equation for the incompressible fluid is satisfied exactly by the computed velocity field. Projections of the stability equation onto the space of inviscid vector fields eliminated the pressure term from the equation and reduces the differential eigen problem to a generalized matrix eigen problem. Results are presented both for the one dimensional, columnar vortex flows and also for the wavy bubble flow.

Mhuiris, N. M. G.↗

Classical free-streamline flow over a polygonal obstacle

In classical Kirchhoff flow, an ideal incompressible fluid flows past an obstacle and around a motionless wake bounded by free streamlines. Since 1869 it has been known that in principle, the two-dimensional Kirchhoff flow over a polygonal obstacle can be determined by constructing a conformal map onto a polygon in the log-hodograph plane. In practice, however, this idea has rarely been put to use except for very simple obstacles, because the conformal mapping problem has been too difficult. This paper presents a practical method for computing flows over arbitrary polygonal obstacles to high accuracy in a few seconds of computer time. We achieve this high speed and flexibility by working with a modified Schwarz-Christoffel integral that maps onto the flow region directly rather than onto the log-hodograph polygon. This integral and its associated parameter problem are treated numerically by methods developed earlier by Trefethen for standard Schwarz-Christoffel maps.

Elcrat, A. R.↗

In situ measurements of the plasma bulk velocity near the Io flux tube

Data obtained with the Voyager Plasma Science Experiment (PSE) during a flyby of the Jovian moon Io are analyzed to quantify the bulk plasma parameters near the Io flux tube. The PSE contains four modulated grid Faraday cups, three pentagonally shaped and the other circular. The path of the spacecraft around Jupiter and past Io ensured that the maximum flux was passing perpendicular to the spacecraft instruments when Voyager was closest to Io (20,000 km). A total of 11 high resolution (M mode) plasma spectra were obtained, revealing 6 species of ions: H(+), O(+), O(2+), S(+), S(2+) and S(3+). All species had the same bulk velocity and convected maxwellian distributions; the heavy ions species were all the same temperature. The mass density of the plasma near Io was 20,000-25,000 amu/cu cm and the electron density was 1300-1700/cu cm. Calculations of the Alfven speed yielded a value ranging from 250-350 km/sec and an Alfven Mach no. of 0.17-0.21. The data were equivalent to that expected from the flow of an incompressible fluid around a long cylinder with a radius about 1.25 that of Io's.

Barnett, A.↗

Viscous starting jets

The transient motion which is produced when a viscous incompressible fluid is forced from an initial state of rest is studied. The equations for unsteady particle paths, written in terms of similarity variables, are analyzed as a quasi-autonomous system with the Reynolds number treated as a parameter. By finding and classifying critical points in the system's phase portrait, the flow structure is examined. It is shown that: (1) bifurcations in the phase portrait occur at specific values of the Reynolds number of the flow in question, and (2) the exact solutions of the Stokes equations for the low-Reynolds-number limit contain two critical Reynolds numbers and three distinct states of motion which culminate in the onset of a vortex roll-up.

Cantwell, Brian J.↗

Tail modeling in a stretched magnetosphere. I - Methods and transformations

A new method is developed for representing the magnetospheric field B as a distorted dipole field. Because Delta-B = 0 must be maintained, such a distortion may be viewed as a transformation of the vector potential A. The simplest form is a one-dimensional 'stretch transformation' along the x axis, concisely represented by the 'stretch function' f(x), which is also a convenient tool for representing features of the substorm cycle. One-dimensional stretch transformations are extended to spherical, cylindrical, and parabolic coordinates and then to arbitrary coordinates. It is shown that distortion transformations can be viewed as mappings of field lines from one pattern to another; the final result only requires knowledge of the field and not of the potentials. General transformations in Cartesian and arbitrary coordinates are derived, and applications to field modeling, field line motion, MHD modeling, and incompressible fluid dynamics are considered.

Stern, David P.↗

Solitary waves in the resonant phenomenon between a surface gravity wave packet and an internal gravity wave

A two-layer inviscid incompressible fluid system of intermediate depth is considered. A multiple-scales perturbation technique is applied to the basic equations and boundary conditions for a two-layer fluid system to derive a system of weakly nonlinear partial integrodifferential equations governing the resonant interaction between a surface gravity wave packet and an internal gravity wave at an intermediate depth, providing a bridge between the existing shallow and deep fluid theories. The convolution integral term in these equations accounts for the dispersion in the lower-layer fluid. An iterative fast Fourier transform scheme is developed to find solitary wave solutions to this system of equations. The overtaking collision of two pairs of solitary waves, simulated using a spectral method, is found to be inelastic. It is found that the amplitude of the solitary waves changes slightly after the collision. The phase shifts these solitary waves undergo was calculated numerically.

Sepulveda, Nicasio↗

Viscous damped space structure for reduced jitter

A technique to provide modal vibration damping in high performance space structures was developed which uses less than one once of incompressible fluid. Up to 50 percent damping can be achieved which can reduce the settling times of the lowest structural mode by as much as 50 to 1. This concept allows the designers to reduce the weight of the structure while improving its dynamic performance. Damping by this technique is purely viscous and has been shown by test to be linear over 5 orders of input magnitude. Amplitudes as low as 0.2 microinch were demonstrated. Damping in the system is independent of stiffness and relatively insensitive to temperature.

Wilson, James F.↗

Spectral element methods: Algorithms and architectures

Spectral element methods are high-order weighted residual techniques for partial differential equations that combine the geometric flexibility of finite element methods with the rapid convergence of spectral techniques. Spectral element methods are described for the simulation of incompressible fluid flows, with special emphasis on implementation of spectral element techniques on medium-grained parallel processors. Two parallel architectures are considered: the first, a commercially available message-passing hypercube system; the second, a developmental reconfigurable architecture based on Geometry-Defining Processors. High parallel efficiency is obtained in hypercube spectral element computations, indicating that load balancing and communication issues can be successfully addressed by a high-order technique/medium-grained processor algorithm-architecture coupling.

Fischer, Paul↗

Star-forming instabilities of a decelerating plane-parallel slab of finite thickness

In order to study how supernova blast waves might catalyze star formation, the stability of a slab of decelerating gas of finite thickness is explored. Previous analyses are extended by applying shock-like boundary conditions to the leading edge of the slab and by studying the effects of arbitrary deceleration. Contrary to some earlier claims, it is found that blast waves can indeed accelerate the rate of star formation in the interstellar medium by compressing the interstellar gas and enabling it to cool. Also, it is demonstrated that in an incompressible fluid, the symmetric and antisymmetric modes in the case of zero acceleration transform continuously into Rayleigh-Taylor and gravity-wave modes as acceleration grows more important. The dynamical overstability of an isothermal shock wave and its implications for star formation are discussed.

Voit, G. Mark↗

Diode-Laser Doppler Velocimeter

Diode-laser Doppler velocimeter measures nonintrusively flow of incompressible fluid in narrow tube. New velocimeter rugged, compact, and competitive in cost. Includes three-section optical head mounted on tube containing flow. In slightly different version, beam splitter and mirror used to split laser beam into two beams.

Getzer, Gregory J.↗

Spectral element methods - Algorithms and architectures

Spectral element methods are high-order weighted residual techniques for partial differential equations that combine the geometric flexibility of finite element methods with the rapid convergence of spectral techniques. Spectral element methods are described for the simulation of incompressible fluid flows, with special emphasis on implementation of spectral element techniques on medium-grained parallel processors. Two parallel architectures are considered; the first, a commercially available message-passing hypercube system; the second, a developmental reconfigurable architecture based on Geometry-Defining Processors. High parallel efficiency is obtained in hypercube spectral element computations, indicating that load balancing and communication issues can be successfully addressed by a high-order technique/medium-grained processor algorithm-architecture coupling.

Fischer, Paul↗

Numerical investigation of the control of the secondary instability process in boundary layers

A numerical model has been developed for investigating boundary layer transition control for a flat plate boundary layer. Active control of a periodically forced boundary layer in an incompressible fluid is studied using surface heating techniques. The spatially evolving boundary layer is simulated. The Navier-Stokes and energy equations are integrated using a fully implicit finite difference/spectral method. Temperature perturbations are introduced locally along finite heater strips to directly attenuate the instability waves in the flow. A feedback control loop is employed in which a downstream sensor is used to monitor wall shear stress fluctuations. Active control of small amplitude two-dimensional and three-dimensional disturbances is numerically simulated. With proper phase control, in-phase reinforcement and out-of-phase attenuation are demonstrated. A receptivity study of the localized temperature perturbations is made. It is shown that narrow heater strips are more receptive in that they maximize the amplitude level of the disturbances in the flow. Active control of the early stages of the fundamental breakdown process is also numerically simulated. Control is achieved with either two-dimensional or three-dimensional control inputs.

Kral, L. D.↗

Conforming Chebyshev spectral collocation methods for the solution of laminar flow in a constricted channel

The numerical simulation of steady planar two-dimensional, laminar flow of an incompressible fluid through an abruptly contracting channel using spectral domain decomposition methods is described. The key features of the method are the decomposition of the flow region into a number of rectangular subregions and spectral approximations which are pointwise C(1) continuous across subregion interfaces. Spectral approximations to the solution are obtained for Reynolds numbers in the range 0 to 500. The size of the salient corner vortex decreases as the Reynolds number increases from 0 to around 45. As the Reynolds number is increased further the vortex grows slowly. A vortex is detected downstream of the contraction at a Reynolds number of around 175 that continues to grow as the Reynolds number is increased further.

Karageorghis, Andreas↗

Direct Finite-Difference Simulations Of Turbulent Flow

Report discusses use of upwind-biased finite-difference numerical-integration scheme to simulate evolution of small disturbances and fully developed turbulence in three-dimensional flow of viscous, incompressible fluid in channel. Involves use of computational grid sufficiently fine to resolve motion of fluid at all relevant length scales.

Rai, Man Mohan↗

The effect of wall compliance on the Goertler vortex instability

The stability of the flow of a viscous incompressible fluid over a curved compliant wall to longitudinal Goertler vortices is investigated. The compliant wall is modeled by a particularly simple equation relating the induced wall displacement to the pressure in the overlying fluid. Attention is restricted to the large Goertler number regime; this regime being appropriate to the most unstable Goertler mode. The effect of wall compliance on this most unstable mode is investigated.

Denier, J. P.↗

Multitasking the INS3D-LU code on the Cray Y-MP

This paper presents the results of multitasking the INS3D-LU code on eight processors. The code is a full Navier-Stokes solver for incompressible fluid in three dimensional generalized coordinates using a lower-upper symmetric-Gauss-Seidel implicit scheme. This code has been fully vectorized on oblique planes of sweep and parallelized using autotasking with some directives and minor modifications. The timing results for five grid sizes are presented and analyzed. The code has achieved a processing rate of over one Gflops.

Fatoohi, Rod↗

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.↗