Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “VISCOSITY”

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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Lattice Boltzmann model for compressible fluids

A lattice Boltzmann model is derived which simulates compressible fluids. By choosing the parameters of the equilibrium distribution appropriately, the sound speed (which may be set arbitrarily low), bulk viscosity, and kinematic viscosity can be selected. This model simulates compressible flows and can include shocks. With a proper rescaling and zero-sound speed, this model simulates Burgers's equation. The viscosity determined by a Chapman-Enskog expansion compares well with that measured form simulations. The exact solutions of Burgers's equation on the unit circle are compared to solutions of lattice Boltzmann model finding reasonable agreement.

Alexander, F. J.↗
NASA NTRS1993

Conductivity As A Measure Of Degree Of Polymerization

In improved method of dielectric monitoring of process of polymerization of mixed ingredients of polymeric material, emphasis placed on measurement of conductivity rather than permittivity. Conductivity tends to change more during curing process than permittivity does, is less dependent on frequency, and more-sensitive and more-reliable indicator of progress of cure. Conductivity used to compute quantity called "ion viscosity," which relates to classical viscosity of fluid. During cure of typical polymer, ion viscosity varies over wide range, eventually rising to plateau toward end of cure.

Dean, David L.↗
NASA NTRS1992

What can we learn about impact mechanics from large craters on Venus?

More than 50 unequivocal peak-ring craters and multiringed impact basins have been identified on Venus from Earth-based Arecibo, Venera 15/16, and Magellan radar images. These ringed craters are relatively pristine, and so serve as an important new dataset that will further understanding of the structural and rheological properties of the venusian surface and of impact mechanics in general. They are also the most direct analogues for craters formed on the Earth in Phanerozoic time. Finite-element simulations of basin collapse and ring formation were undertaken in collaboration with V. J. Hillgren (University of Arizona). These calculations used an axisymmetric version of the viscoelastic finite element code TECTON, modeled structures on the scale of Klenova or Meitner, and demonstrated two major points. First, viscous flow and ring formation are possible on the timescale of crater collapse for the sizes of multiringed basins seen on Venus and heat flows appropriate to the plant. Second, an elastic lithosphere overlying a Newtonian viscous asthenosphere results mainly in uplift beneath the crater. Inward asthenospheric flow mainly occurs at deeper levels. Lithospheric response is dominantly vertical and flexural. Tensional stress maxima occur and ring formation by normal faulting is predicted in some cases, but these predicted rings occur too far out to explain observed ring spacings on Venus (or on the Moon). Overall, these estimates and models suggest that multiringed basin formation is indeed possible at the scales observed on Venus. Furthermore, due to the strong inverse dependence of solid-state viscosity on stress, the absence of Cordilleran-style ring faulting in craters smaller than Meitner or Klenova makes sense. The apparent increase in viscosity of shock-fluidized rock with crater diameter, greater interior temperatures accessed by larger, deeper craters, and decreased non-Newtonian viscosity associated with larger craters may conspire to make the transition with diameter from peak-ring crater to Orientale-type multiringed basin rather abrupt.

Mckinnon, William B.↗
NASA NTRS1992

The formation of ore mineral deposits on the Moon: A feasibility study

Most of the ore deposits on Earth are the direct result of formation by hydrothermal solutions. Analogous mineral concentrations do not occur on the Moon, however, because of the absence of water. Stratified ore deposits form in layered instrusives on Earth due to fractional crystallization of magma and crystal settling of high-density minerals, particularly chromium in the mineral chromite. We have evaluated the possibility of such mineral deposition on the Moon, based upon considerations of 'particle settling velocities' in lunar vs. terrestrial magmas. A first approximation of Stoke's Law would seem to indicate that the lower lunar gravity (1/6 terrestrial) would result in slower crystal settling on the Moon. However, the viscosity of the silicate melt is the most important factor affecting the settling velocity. The viscosities of typical lunar basaltic melts are 10-100 times less than their terrestrial analogs. These lower viscosities result from two factors: (1) lunar basaltic melts are typically higher in FeO and lower in Al2O3, Na2O, and K2O than terrestrial melts; and (2) lunar igneous melts and phase equilibria tend to be 100-150 C higher than terrestrial, largely because of the general paucity of water and other volatile phases on the Moon. Therefore, particle settling velocities on the Moon are 5-10 times greater than those on Earth. It is highly probable that stratiform ore deposits similar to those on Earth exist on the Moon. The most likely ore minerals involved are chromite, ilmenite, and native FeNi metal. In addition, the greater settling velocities of periodotite in lunar magmas indicate that the buoyancy effects of the melt are less than on Earth. Consequently, the possibility is considerably less than on Earth of deep-seated volcanism transporting upper mantle/lower crustal xenoliths to the surface of the Moon, such as occurs in kimberlites on Earth.

Taylor, Lawrence A.↗
NASA NTRS1992

Venusian sinuous rilles

After a preliminary assessment of venusian channels, it now seems to be clear that the channels have distinctive classes, which imply a wide range of formation parameters and formation mechanisms. They include outflow channels mainly formed by mechanical erosion from very high discharge flow, and canali-type channels requiring either constructional process or mechanical erosion by rather exotic low-viscosity lava such as carbonatite or sulfur. Here we focus on venusian sinuous rilles. Venusian sinuous rilles are generally simple, and originate from a collapsed source. They are shallow and narrow downstream. The venusian sinuous rilles are distinct from canali-type channels, which exhibit almost constant morphologies throughout their entire length, and from outflow channels, which are characterized by wide anastomosing reaches. The lunar sinuous rilles could have been formed initially as constructional channels. However, incision was caused by the long flow duration and high temperatures of eruption, along with relatively large discharge rates, possibly assisted by a low viscosity of the channel-forming lava. Channel narrowing and levee formation suggest relatively fast cooling. The venusian channels could have had a similar sequence of formation including rapid cooling. Assuming the substrate is typical tholeiitic lava, the flowing lavas' temperatures have to be higher than the melting temperature of the substrate. The flow should have a low viscosity to cause turbulence and keep a high Reynolds number to sustain efficient thermal erosion. Determining eruption conditions also provide insights to estimate lava composition. Assuming a channel is formed mostly by thermal erosion, the channel's length and longitudinal profile are functions of lava properties. The depth profiles of the channel are measured by radar foreshortening methods and stereo images. Eruption conditions of channel forming lava can be estimated by the methods developed by Hulme.

Komatsu, G.↗
NASA NTRS1992

Numerical Modeling of a Vortex Stabilized Arcjet

Arcjet thrusters are being actively considered for use in Earth orbit maneuvering applications. Experimental studies are currently the chief means of determining an optimal thruster configuration. Earlier numerical studies have failed to include all of the effects found in typical arcjets including complex geometries, viscosity, and swirling flow. Arcjet geometries are large area ratio converging nozzles with centerbodies in the subsonic portion of the nozzle. The nozzle walls serve as the anode while the centerbody functions as the cathode. Viscous effects are important because the Reynolds number, based on the throat radius, is typically less than 1,000. Experimental studies have shown that a swirl or circumferential velocity component stabilizes a constricted arc. This dissertation describes the equations governing flow through a constricted arcjet thruster. An assumption that the flowfield is in local thermodynamic equilibrium leads to a single fluid plasma temperature model. An order of magnitude analysis reveals the governing fluid mechanics equations are uncoupled from the electromagnetic field equations. A numerical method is developed to solve the governing fluid mechanics equations, the Thin Layer Navier-Stokes equations. A coordinate transformation is employed in deriving the governing equations to simplify the application of boundary conditions in complex geometries. An axisymmetric formulation is employed to include the swirl velocity component as well as the axial and radial velocity components. The numerical method is an implicit finite-volume technique and allows for large time steps to reach a converged steady-state solution. The inviscid fluxes are flux-split, and Gauss-Seidel line relaxation is used to accelerate convergence. Converging-diverging nozzles with exit-to-throat area ratios up to 100:1 and annular nozzles were examined. Quantities examined included Mach number and static wall pressure distributions, and oblique shock structures. As the level of swirl and viscosity in the flowfield increased the mass flow rate and thrust decreased. The technique was used to predict the flow through a typical arcjet thruster geometry. Results indicate swirl and viscosity play an important role in the complex geometry of an arcjet.

Pawlas, Gary E.↗
NASA NTRS1992

On the nonlinear three dimensional instability of Stokes layers and other shear layers to pairs of oblique waves

The nonlinear evolution of a pair of initially oblique waves in a high Reynolds Number Stokes layer is studied. Attention is focused on times when disturbances of amplitude epsilon have O(epsilon(exp 1/3)R) growth rates, where R is the Reynolds number. The development of a pair of oblique waves is then controlled by nonlinear critical-layer effects. Viscous effects are included by studying the distinguished scaling epsilon = O(R(exp -1)). This leads to a complicated modification of the kernel function in the integro-differential amplitude equation. When viscosity is not too large, solutions to the amplitude equation develop a finite-time singularity, indicating that an explosive growth can be introduced by nonlinear effects; we suggest that such explosive growth can lead to the bursts observed in experiments. Increasing the importance of viscosity generally delays the occurrence of the finite-time singularity, and sufficiently large viscosity may lead to the disturbance decaying exponentially. For the special case when the streamwise and spanwise wavenumbers are equal, the solution can evolve into a periodic oscillation. A link between the unsteady critical-layer approach to high-Reynolds-number flow instability, and the wave vortex approach is identified.

Wu, Xuesong↗
NASA NTRS1992

Turbulent transport of a passive-scalar field by using a renormalization-group method

A passive-scalar field is considered to evolve under the influence of a turbulent fluid governed by the Navier-Stokes equation. Turbulent-transport coefficients are calculated by small-scale elimination using a renormalization-group method. Turbulent processes couple both the viscosity and the diffusivity. In the absence of any correlation between the passive-scalar fluctuations and any component of the fluid velocity, the renormalized diffusivity is essentially the same as if the fluid velocity were frozen, although the renormalized equation does contain higher-order nonlinear terms involving viscosity. This arises due to the nonlinear interaction of the velocity with itself. In the presence of a finite correlation, the turbulent diffusivity becomes coupled with both the velocity field and the viscosity. There is then a dependence of the turbulent decay of the passive scalar on the turbulent Prandtl number.

Hossain, Murshed↗
NASA NTRS1992

Axisymmetric spherical shell models of mantle convection with variable properties and free and rigid lids

Axisymmetric spherical shell numerical simulations of mantle convection were carried out to investigate the influence of two end-member surface stress conditions: stress-free and rigid. These correspond approximately to a subducting or a rigid lithosphere and can be seen as end-member models of the surface of Venus. Our model assumed an effective Rayleigh number of 3 x 10 exp 6, similar to that for earth, and included uniform internal heating and depth-dependent thermal expansivity and thermal conductivity. The simulations utilized a Newtonian viscosity which was constant or varied with depth and/or temperature. We show how the temperature, speed, and vorticity fields change qualitatively and quantitatively with surface temperature, surface stress condition, internal heating and viscosity distribution. We find that a rigid lid and viscosity which increases with depth both promote steady large-scale circulation with smaller-scale circulation in the upper mantle.

Leitch, A. M.↗
NASA NTRS1993

Deep mantle viscous structure with prior estimate and satellite constraint

A radially stratified and incompressible earth model with secular variation of the second degree gravity field J-dot(2) is used here to test sensitivity of data to viscosity increases with depth and convective boundary layer structure. Prior estimates and observed nontidal J-dot(23)(-C-dot(20)) are consistent with a layered lower mantle viscosity. Details of this layering are examined by comparing predicted and observed J-dot(3), J-dot(4). Speculation that a high-velocity layer exists above D-double prime is considered. With a 650-km thick deep high-viscosity layer, the remaining lower mantle is in one of two ranges: 1.5 to 3.5 x 10 exp 20 or 3.5 to about 10 exp 22 Pa s.

Ivins, E. R.↗
NASA NTRS1993

The structure and stability of transonic accretion disks surrounding black holes

Stationary transonic alpha-viscosity models of accretion disks surrounding nonrotating black holes have been investigated. The viscosity is modified such that it vanishes in the supersonic region to ensure its effect does not violate the causality condition. In contrast to previous studies, the viscous stress is taken to be explicitly proportional to the angular velocity gradient and is not assumed to depend solely on the local pressure in the disk. The numerical results reveal that the structure of the innermost regions of the disk are more sensitive to the modified form of the viscosity than to the form of the viscous stress. The critical sonic point is located inside the innermost stable circular orbit of a test particle at 3 Schwarzschild radii. In these solutions, the transition from subsonic to supersonic flow results from pressure effects and not viscous effects. The linear stability of these disks has been examined in the local approximation. It is found that radiative energy transport and viscous stresses in the radial direction can have important effects. As a result, it is shown that the growth rate of the inertial-acoustic mode reaches a maximum at a critical wavelength.

Chen, Xingming↗
NASA NTRS1993

Central difference TVD schemes for time dependent and steady state problems

A Runge-Kutta formula in time is presently used to advance schemes in which central differences are used to solve the time-dependent Euler equations; a second difference is added near shocks as an artificial viscosity to reduce the given scheme to a first-order upwind one at shocks. A matrix-valued dissipation is introduced and compared with the scalar viscosity; a connection is shown between this artificial viscosity and flux limiters. The use of various flux limiters for this central difference scheme is compared.

Jorgenson, P.↗
NASA NTRS1950

Foaming of mixtures of pure hydrocarbons

Mixtures of pure liquid hydrocarbons are capable of foaming. Nine hydrocarbons were mixed in pairs, in all possible combinations, and four proportions of each combination. These mixtures were sealed in glass tubes, and the foaming was tested by shaking. Mixtures of aliphatic with other aliphatic hydrocarbons, or of alkyl benzenes with other alkyl benzenes, did not foam. Mixtures of aliphatic hydrocarbons with alkyl benzenes did foam. The proportions of the mixtures greatly affected the foaming, the maximum foaming of 12 of 20 pairs being at the composition 20 percent aliphatic hydrocarbon, 80 percent alkyl benzene. Six seconds was the maximum foam lifetime of any of these mixtures. Aeroshell 120 lubricating oil was fractionated into 52 fractions and a residue by extraction with acetone in a fractionating extractor. The index of refraction, foam lifetime, color, and viscosity of these fractions were measured. Low viscosity and high index fractions were extracted first. The viscosity of the fractions extracted rose and the index decreased as fractionation proceeded. Foam lifetimes and color were lowest in the middle fractions. Significance is attached to the observation that none of the foam lifetimes of the fractions or residue is as high as the foam lifetime of the original Aeroshell, indicating that the foaming is not due to a particular foaming constituent, but rather to the entire mixture.

LUBRICATING OILS - DISTILLATION↗
NASA NTRS1955

Effect of Surface-active Additives on Physical Properties of Slurries of Vapor-process Magnesium

The presence of 3 to 5 percent surface-active additive gave the lowest Brookfield apparent viscosity, plastic viscosity, and yield value that were obtained for slurry fuels containing approximately 50 percent vapor-process magnesium in JP-1 fuel. The slurries settled little and were easily remixed. A polyoxyethylene dodecyl alcohol was the most effective of 13 additives tested in reducing the Brookfield apparent viscosity and the yield value of the slurry. The seven most effective additives all had a hydroxyl group plus an ester or polyoxethylene group in the molecule. The densities of some of the slurries were measured.

Pinns, Murray L↗
NASA NTRS1993

Hypersonic flow computations around re-entry vehicles

The development of an algorithm for the solution of the compressible Euler equations at high Mach numbers on unstructured tetrahedral meshes is described. The basic algorithm is constructed in the form of a central difference scheme plus an explicit added artificial viscosity based upon fourth order differences of the solution. The stability of the solution in the vicinity of strong gradients is preserved by the incorporation of an additional artificial viscosity based upon a second order difference. Higher order accuracy is regained by using the ideas of flux corrected transport to limit the amount of added viscosity. The solution is advanced to steady state by means of an explicit multi-stage time-stepping method. The computational efficiency of the complete process is improved by incorporating an unstructured multigrid acceleration procedure. A number of flows of practical interest are analyzed to demonstrate the numerical performance of the proposed approach.

Peraire, J.↗
NASA NTRS1993

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 1: Analysis description

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. The governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models are described in detail.

Towne, Charles E.↗
NASA NTRS1993

Megaplumes on Venus

The geoid and topography heights of Atla Regio and Beta Regio, both peaks and slopes, appear explicable as steady-state plumes, if non-linear viscosity eta(Tau, epsilon) is taken into account. Strongly constrained by the data are an effective plume depth of about 700 km, with a temperature anomaly thereat of about 30 degrees, leading to more than 400 degrees at the plume head. Also well constrained is the combination Q(eta)/s(sup 4)(sub 0) = (volume flow rate)(viscosity)/(plume radius): about 11 Pa/m/sec. The topographic slopes dh/ds constrain the combination Q/A, where A is the thickness of the spreading layer, since the slope varies inversely with velocity. The geoid slopes dN/ds require enhancement of the deeper flow, as expected from non-linear viscosity. The Beta data are best fit by Q = 500 m(sup 3)/sec and A equals 140 km; the Atla, by Q equals 440 m(exp 3)/sec and A equals 260 km. The dynamic contribution to the topographic slope is minor.

Kaula, W. M.↗
NASA NTRS1993

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A computer code called Proteus 3D has been developed to solve the three dimensional, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort has been to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation have been emphasized. The governing equations are solved in generalized non-orthogonal body-fitted coordinates by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. It describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗