SEARCH · Search NASA

Results for “hydrodynamics methods”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Removal of spurious modes encountered in solving stability problems by spectral methods

A technique based on the Galerkin approximation is developed to remove spurious roots arising when Chebyshev spectral methods are used to solve eigenvalue problems in hydrodynamic stability. The derivation of Galerkin-Chebyshev approximations is explained, and numerical results for the Orr-Sommerfeld equations of plane Poiseuille flow and a Blasius profile are presented in tables and compared with those obtained by the method of Zebib (1984). It is pointed out that the present method does not increase the size of the algebraic system to be solved.

Zebib, Abdelfattah↗

NASA NTRS1993

A mixed method Poisson solver for three-dimensional self-gravitating astrophysical fluid dynamical systems

A key ingredient in the simulation of self-gravitating astrophysical fluid dynamical systems is the gravitational potential and its gradient. This paper focuses on the development of a mixed method multigrid solver of the Poisson equation formulated so that both the potential and the Cartesian components of its gradient are self-consistently and accurately generated. The method achieves this goal by formulating the problem as a system of four equations for the gravitational potential and the three Cartesian components of the gradient and solves them using a distributed relaxation technique combined with conventional full multigrid V-cycles. The method is described, some tests are presented, and the accuracy of the method is assessed. We also describe how the method has been incorporated into our three-dimensional hydrodynamics code and give an example of an application to the collision of two stars. We end with some remarks about the future developments of the method and some of the applications in which it will be used in astrophysics.

Duncan, Comer↗

NASA NTRS1952

Theoretical damping in roll and rolling moment due to differential wing incidence for slender cruciform wings and wing-body combinations

A method of analysis based on slender-wing theory is developed to investigate the characteristics in roll of slender cruciform wings and wing-body combinations. The method makes use of the conformal mapping processes of classical hydrodynamics which transform the region outside a circle and the region outside an arbitrary arrangement of line segments intersecting at the origin. The method of analysis may be utilized to solve other slender cruciform wing-body problems involving arbitrarily assigned boundary conditions. (author)

Adams, Gaynor J↗

NASA NTRS1981

Effect of Surface Roughness on Hydrodynamic Bearings

A theoretical analysis on the performance of hydrodynamic oil bearings is made considering surface roughness effect. The hydrodynamic as well as asperity contact load is found. The contact pressure was calculated with the assumption that the surface height distribution was Gaussian. The average Reynolds equation of partially lubricated surface was used to calculate hydrodynamic load. An analytical expression for average gap was found and was introduced to modify the average Reynolds equation. The resulting boundary value problem was then solved numerically by finite difference methods using the method of successive over relaxation. The pressure distribution and hydrodynamic load capacity of plane slider and journal bearings were calculated for various design data. The effects of attitude and roughness of surface on the bearing performance were shown. The results are compared with similar available solution of rough surface bearings. It is shown that: (1) the contribution of contact load is not significant; and (2) the hydrodynamic and contact load increase with surface roughness.

Majumdar, B. C.↗

NASA NTRS1975

Intrinsic hydrodynamics with applications to space-time fluids

Cartan's methods of exterior differential forms are utilized to derive conformal and absolute conservation equations for various field quantities in relativistic hydromechanics. In particular, Fridman's results for the conservation of vector lines are extended to space time situations. The invariance concepts of Helmholtz and Bernoulli are formulated for arbitrary frames of reference, and the separate conservation concepts of total mass and total number in space time are distinguished and associated with the concept of turbulence.

Kiehn, R. M.↗

NASA NTRS1981

Analysis of starvation effects on hydrodynamic lubrication in nonconforming contacts

The effects of lubricant starvation on minimum film thickness, under conditions of a hydrodynamic point contact, are determined by numerical methods where: (1) starvation is effected by varying the fluid inlet level; (2) the Reynolds boundary conditions are applied at the cavitation boundary; and (3) zero pressure is stipulated at the meniscus or inlet boundary. Seventy-four cases were used to numerically determine a minimum-film-thickness equation, as a function of the ratio of dimensionless load to dimensionless speed for varying degrees of starvation. A film reduction factor was in turn determined as a function of the fluid inlet level, and a starved, fully-flooded boundary was defined along with an expression determining the onset of starvation. It is found that as the degree of starvation increases, the minimum film thickness decreases gradually until the fluid inlet becomes critical. The changes in the inlet pressure buildup due to changing the available lubricant supply are presented in the form of three-dimensional isometric plots.

Brewe, D. E.↗

NASA NTRS1993

Simplified Computation Of Blade-Surface Coordinates

Coordinates of points on complicated three-dimensional surface of turbo-machine blade computed by method partly simplifying three-dimensional problem to two-dimensional problem. Method implemented in computer program enabling rapid hydrodynamical evaluation and optimization of blade surface prior to definition of final design of, and analysis of stresses in blade. By incorporating program into procedure for designing turbomachine, one reduces time spent on design iterations.

Chen, Wei-Chung↗

NASA NTRS1994

X-ray clusters from a high-resolution hydrodynamic PPM simulation of the cold dark matter universe

A new three-dimensional hydrodynamic code based on the piecewise parabolic method (PPM) is utilized to compute the distribution of hot gas in the standard Cosmic Background Explorer (COBE)-normalized cold dark matter (CDM) universe. Utilizing periodic boundary conditions, a box with size 85 h(exp-1) Mpc, having cell size 0.31 h(exp-1) Mpc, is followed in a simulation with 270(exp 3)=10(exp 7.3) cells. Adopting standard parameters determined from COBE and light-element nucleosynthesis, Sigma(sub 8)=1.05, Omega(sub b)=0.06, we find the X-ray-emitting clusters, compute the luminosity function at several wavelengths, the temperature distribution, and estimated sizes, as well as the evolution of these quantities with redshift. The results, which are compared with those obtained in the preceding paper (Kang et al. 1994a), may be used in conjuction with ROSAT and other observational data sets. Overall, the results of the two computations are qualitatively very similar with regard to the trends of cluster properties, i.e., how the number density, radius, and temeprature depend on luminosity and redshift. The total luminosity from clusters is approximately a factor of 2 higher using the PPM code (as compared to the 'total variation diminishing' (TVD) code used in the previous paper) with the number of bright clusters higher by a similar factor. The primary conclusions of the prior paper, with regard to the power spectrum of the primeval density perturbations, are strengthened: the standard CDM model, normalized to the COBE microwave detection, predicts too many bright X-ray emitting clusters, by a factor probably in excess of 5. The comparison between observations and theoretical predictions for the evolution of cluster properties, luminosity functions, and size and temperature distributions should provide an important discriminator among competing scenarios for the development of structure in the universe.

Bryan, Greg L.↗

NASA NTRS1994

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 2: Homogeneous medium at Mach 3 with gamma = 5/3

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing spheres of varying sizes (from 10 down to 0.01 accretion radii) move at Mach 3 relative to a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3). To accommodate the long-range gravitational forces, the extent of the computational volume is 32(exp 3) accretion radii. We examine the influence of numerical procedure on physical behavior. The hydrodynamics is modeled by the 'piecewise parabolic method.' No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (5-10) grids around the sphere, each finer grid being a factor of 2 smaller in zone dimension that the next coarser grid. The largest dynamic range (ratio of size of the largest grid to size of the finest zone) is 16,384. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid, while at the same time evolving the gas on the coarser grids. Initially (at time t = 0-10), a shock front is set up, a Mach cone develops, and the accretion column is observable. Eventually the flow becomes unstable, destroying axisymmetry. This happens approximately when the mass accretion rate reaches the values (+/- 10%) predicted by the Bondi-Hoyle accretion formula (factor of 2 included). However, our three-dimensional models do not show the highly dynamic flip-flop flow so prominent in two-dimensional calculations performed by other authors. The flow, and thus the accretion rate of all quantities, shows quasi-periodic (P approximately equals 5) cycles between quiescent and active states. The interpolation formula proposed in an accompanying paper is found to follow the collected numerical data to within approximately 30%. The specific angular momentum accreted is of the same order of magnitude as the values previously found for two-dimensional flows.

Ruffert, Maximilian↗

NASA NTRS1996

Numerical Investigation of Two-Phase Flows With Charged Droplets in Electrostatic Field

A numerical method to solve two-phase turbulent flows with charged droplets in an electrostatic field is presented. The ensemble-averaged Navier-Stokes equations and the electrostatic potential equation are solved using a finite volume method. The transitional turbulence field is described using multiple-time-scale turbulence equations. The equations of motion of droplets are solved using a Lagrangian particle tracking scheme, and the inter-phase momentum exchange is described by the Particle-In-Cell scheme. The electrostatic force caused by an applied electrical potential is calculated using the electrostatic field obtained by solving a Laplacian equation and the force exerted by charged droplets is calculated using the Coulombic force equation. The method is applied to solve electro-hydrodynamic sprays. The calculated droplet velocity distributions for droplet dispersions occurring in a stagnant surrounding are in good agreement with the measured data. For droplet dispersions occurring in a two-phase flow, the droplet trajectories are influenced by aerodynamic forces, the Coulombic force, and the applied electrostatic potential field.

Kim, Sang-Wook↗

NASA NTRS2000

Remobilizing the Interfaces of Thermocapillary Driven Bubbles Retarded by the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is applied to the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow which propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillarity to direct the movement of bubbles in space is the fact that surfactant impurities which are unavoidably present in the continuous phase can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature tension gradient, and decreases to near zero the thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have been retarded by the adsorption onto the bubble surface of a surfactant impurity, Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition the rapid kinetic exchange and high bulk concentration maintain a saturated surface with a uniform surface concentrations. This prevents retarding surface tension gradients and keeps the velocity high. In our first report last year, we detailed experimental results which verified the theory of remobilization in ground based experiments in which the steady velocity of rising bubbles was measured in a continuous phase consisting of a glycerol/water mixture containing a polyethylene glycol surfactant C12E6 (CH3(CH2)11(OCH2CH2)6OH). In our report this year, we detail our efforts to describe theoretically the remobilization observed. We construct a model in which a bubble rises steadily by buoyancy in a continuous (Newtonian) viscous fluid containing surfactant with a uniform far field bulk concentration. We account for the effects of inertia as well as viscosity in the flow in the continuous phase caused by the bubble motion (order one Reynolds number), and we assume that the bubble shape remains spherical (viscous and inertial forces are smaller than capillary forces, i e. small Weber and capillary numbers). The surfactant distribution is calculated by solving the mass transfer equations including convection and diffusion in the bulk, and finite kinetic exchange the bulk and the surface. Convective effects dominate diffusive mass transfer in the bulk of the liquid (high Peclet numbers) except in a thin boundary layer near the surface. A finite volume method is used to numerically solve the hydrodynamic and mass transfer equations on a staggered grid which accounts specifically for the thin boundary layer. We present the results of the nondimensional drag as a function of the bulk concentration of surfactant for different rates of kinetic exchange, from which we develop criteria for the concentration necessary to develop a prescribed degree of remobilization. The criteria compare favorably with the experimental results.

Palaparthi, Ravi↗

NASA NTRS1970

Lagrangian description of warm plasmas

Efforts are described to extend the averaged Lagrangian method of describing small signal wave propagation and nonlinear wave interaction, developed by earlier workers for cold plasmas, to the more general conditions of warm collisionless plasmas, and to demonstrate particularly the effectiveness of the method in analyzing wave-wave interactions. The theory is developed for both the microscopic description and the hydrodynamic approximation to plasma behavior. First, a microscopic Lagrangian is formulated rigorously, and expanded in terms of perturbations about equilibrium. Two methods are then described for deriving a hydrodynamic Lagrangian. In the first of these, the Lagrangian is obtained by velocity integration of the exact microscopic Lagrangian. In the second, the expanded hydrodynamic Lagrangian is obtained directly from the expanded microscopic Lagrangian. As applications of the microscopic Lagrangian, the small-signal dispersion relations and the coupled mode equations are derived for all possible waves in a warm infinite, weakly inhomogeneous magnetoplasma, and their interactions are examined.

Kim, H.↗

NASA NTRS1992

A new scheme for multidimensional line transfer. II - ETLA method in one dimension with application to iron K-alpha lines

The 1D version of the radiation hydrodynamics computer code ALTAIR is discribed. It is an implementation of the equivalent two-level atom method, and also incorporates a variety of other iterative techniques to achieve a fast, accurate solution of the problem of coupled atomic kinetics and radiative transfer. All the techniques are described in sufficient detail to enable their use by others. The application of the program is illustrated by a set of calculations of the iron X-ray emission spectrum from intercloud hot gas that may exist in the broad-line region of active galactic nuclei. These calculations demonstrate that such spectra may be a rich load to mine for information about the hot gas in AGNs, and also that the concept of 'the iron K-alpha line' inadequately depicts the spectrum.

Castor, John I.↗

NASA NTRS1984

Real Surface Effects in Elastohydrodynamic Lubrication

Following a general discussion of real effects in EHL, consideration is given in detail to the role played by stochastic roughness superimposed on the nominal solid boundaries. Particular attention is given to the full-film EHL regime where incipient asperity contact bears a negligible fraction of the load. Flow, from which an averaged Reynolds equation can be formed, is nonetheless modified by the amplitude and texture of the roughness patterns. In describing these effects by means of certain well-defined flow factors, it is found that only two extra parameters are needed - the rms surface height and the ratio of the correlation lengths in the two principal roughness directions. In a lowest order perturbation expansion of these flow factors in powers of the ratio of rms roughness to nominal film thickness, no other properties of the roughness appear. In most cases of practical interest, the factors describing Poiseuille flow are separable into the sum of two single-surface flow factors which means that a combination of a single equivalent rough surface versus an ideal smooth surface can always be found. For flow entrained by slip velocity the factors separate instead into a difference of single-surface factors and it becomes significant which of the two surfaces carries the equivalent roughness. Results are discussed of some applications of the averaged Reynolds equation based on the flow factor method to the EHL line contact problem. Finally, the partial EHL regime is considered where comparable load fractions are carried by the hydrodynamic film and by incipient mechanical contact. An extension of the method into this regime by combining it with asperity contact models appears most encouraging.

John H Tripp↗

NASA NTRS1975

Evolution of rotating interstellar clouds. I - Numerical techniques

A method is described for the numerical calculation of the hydrodynamic evolution of a self-gravitating configuration in two space dimensions with assumed axial symmetry. The calculation is formulated in cylindrical coordinates with respect to a moving Eulerian grid and is solved using explicit hydrodynamics combined with implicit radiative transfer. The physics included is appropriate for calculation of the collapse of a rotating protostellar cloud. The gravitational field is obtained by means of an alternating-direction iterative technique. Numerical tests to demonstrate the correctness of the method are presented for special cases.

Black, D. C.↗

NASA NTRS1985

Flight vehicle propellant surface motion by numerical method

The marker and cell technique is applied to three idealized hydrodynamic situations of interest to the Space Shuttle External Tank (ET). These applications illustrate the versatility of this technique as well as its usefulness in solving similar problems in the aerospace industry. The codes used are the simplified marker and cell (SMAC) code and the variable grid simplified marker and cell code named ERIE. The examples focus on the dropout height in the LH2 tank, the slosh behavior of propellant in the bottom dome of the ET during thrust termination, both of which are required to accurately determine the amount of usable propellants, and the effect on liquid surface due to LO2 tank wall motions occurring during Space Shuttle lift-off. Comparisons with experimental results and linear theory are also presented.

Mehta, G.↗

NASA NTRS1980

Protostellar formation in rotating interstellar clouds. I - Numerical methods and tests

Attention is given to numerical methods and tests of a series of gravitational hydrodynamics computer codes constructed in order to numerically follow the dynamic collapse of spherically symmetric (1-D), axisymmetric (2-D), and non-axisymmetric (3-D) isothermal interstellar clouds. A spherical harmonic expansion is used to solve the Poisson equation for the gravitational potential. It is shown that the use of explicit donor-cell hydrodynamics on a moving spherical coordinate grid ensures mass and momentum conservation and allows the grid to follow the collapse of the fluid. Finally, the performance of the codes is examined.

Boss, A. P.↗

NASA NTRS1983

Filtering flight data prior to aerodynamic system identification

An algorithm for processing flight-test data to provide state estimates and instrument calibrations for aerodynamic or hydrodynamic system identification by the equation error estimation method is developed and demonstrated on synthesized data. The extended-Kalman-filter algorithm employs a locally level, north-pointing frame of reference, accounts for rotating ellipsoidal earth effects, and estimates sensor bias, scale factors, wind components, and process noise levels by maximum-likelihood parameters. The method is found to be most effective with navigation quality inertial input data. The algorithm is applied to data from a six-degree-of-freedom F-4 aircraft simulation and shown to produce state estimates in good agreement with the simulation values.

Trankle, T. L.↗