Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rotating 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 739 records · Page 41

Fragmentation in stellar collapse

A three dimensional Cartesian code was developed to model molecular cloud collapse and protostellar formation. The basic Eulerian scheme employs spatially second-order accurate, finite difference methods to advance the fluid variables in time. The hydrodynamic equations governing the collapse include the effects of self-gravity, rotation, and radiative transfer. Supplementary equations include Poisson's equation which relates the gravitational potential to the density distribution and a mean intensity equation derived with the Eddington approximation for radiative transfer. The gas pressures and internal energies are determined from equations of state derived for a molecular cloud composed of hydrogen, helium, and heavier elements.

Myhill, E. M.↗

A hierarchy for modeling high speed propulsion systems

General research efforts on reduced order propulsion models for control systems design are overviewed. Methods for modeling high speed propulsion systems are discussed including internal flow propulsion systems that do not contain rotating machinery, such as inlets, ramjets, and scramjets. The discussion is separated into four areas: (1) computational fluid dynamics models for the entire nonlinear system or high order nonlinear models; (2) high order linearized models derived from fundamental physics; (3) low order linear models obtained from the other high order models; and (4) low order nonlinear models (order here refers to the number of dynamic states). Included in the discussion are any special considerations based on the relevant control system designs. The methods discussed are for the quasi-one-dimensional Euler equations of gasdynamic flow. The essential nonlinear features represented are large amplitude nonlinear waves, including moving normal shocks, hammershocks, simple subsonic combustion via heat addition, temperature dependent gases, detonations, and thermal choking. The report also contains a comprehensive list of papers and theses generated by this grant.

Hartley, Tom T.↗

A hierarchy for modeling high speed propulsion systems

General research efforts on reduced order propulsion models for control systems design are overviewed. Methods for modeling high speed propulsion systems are discussed including internal flow propulsion systems that do not contain rotating machinery such as inlets, ramjets, and scramjets. The discussion is separated into four sections: (1) computational fluid dynamics model for the entire nonlinear system or high order nonlinear models; (2) high order linearized model derived from fundamental physics; (3) low order linear models obtained from other high order models; and (4) low order nonlinear models. Included are special considerations on any relevant control system designs. The methods discussed are for the quasi-one dimensional Euler equations of gasdynamic flow. The essential nonlinear features represented are large amplitude nonlinear waves, moving normal shocks, hammershocks, subsonic combustion via heat addition, temperature dependent gases, detonation, and thermal choking.

Hartley, Tom T.↗

Analysis of thermo-chemical nonequilibrium models for carbon dioxide flows

The aerothermodynamics of thermochemical nonequilibrium carbon dioxide flows is studied. The chemical kinetics models of McKenzie and Park are implemented in separate three-dimensional computational fluid dynamics codes. The codes incorporate a five-species gas model characterized by a translational-rotational and a vibrational temperature. Solutions are obtained for flow over finite length elliptical and circular cylinders. The computed flowfields are then employed to calculate Mach-Zehnder interferograms for comparison with experimental data. The accuracy of the chemical kinetics models is determined through this comparison. Also, the methodology of the three-dimensional thermochemical nonequilibrium code is verified by the reproduction of the experiments.

Rock, Stacey G.↗

Kinetic structure of intermediate shocks: Implications for the magnetopause

A general study of the structure and stability of intermediate shocks (IS) in an isotropic plasma is presented using a hybrid as well as a resistive Hall MHD code. Special emphasis is put on the question of whether the rotational layers observed at the magnetopause can be intermediate shocks. The shocks are formed dynamically by the interaction between a flowing plasma and a stationary piston. Coplanar ISs (both strong and weak) are found to be stable in a collisionless plasma. The existence of slow shocks in a high beta plasma is also established for the first time. Noncoplanar ISs are found to be time-dependent, evolving toward a rotational discontinuity (RD) after some characteristic time tau which can be quite long (1000 Omega(exp -1), where Omega is the ion gyrofrequency). The value tau is larger the closer the rotation angle is to 180 deg. Rotations larger than 180 deg are found to be unstable, decaying into a state of minimum shear (i.e., rotation angle less than 180 deg). There are various length scales associated with an IS in the kinetic regime. The shortest scale is found to be the length scale over which rotation of the transverse component of the magnetic field takes place. This scale can have a half width as small as one ion inertial length (c/omega(sub p)) for electron sense rotations and 3c/omega(sub p) for ion sense rotations, for an upstream ion beta of unity. Both of these scales are consistent with the observed thickness at the magnetopause and identical to the corresponding RD scales. A detailed study of the mode conversion of the Alfven ion cyclotron waves (A/IC) waves across both slow and intermediate shocks and the resulting downstream wave spectrum are presented. The possibility that the large number of relfected ions observed at the magnetopause may be due to the presence of strong ISs is considered. The identification of strong ISs and their distinction from RDs should be possible in observations due to significant differences that exist between jump conditions and overall structure of the two discontinuities. The jumps in the plasma parameters across a weak IS are typically small. This together with the fact that the weak ISs and RDs have very similar thickness and other overall properties makes the distinction between weak ISs and RDs in the observations largely inconsequential. However, at large noncoplanarity angles the weak IS approaches the RD limit in a relatively short time (approximately less than 100 Omega (exp -1)). Thus, magnetopause rotations with large noncoplanarity angles are most likely either RDs or strong ISs. Finally, direct comparisons between fluid (resistive Hall MHD) and kinetic simulations show that fluid theory is not applicable to study of ISs in a collisionless plasma.

Karimabadi, H.↗

Experimental study of vortex breakdown in a cylindrical, swirling flow

The stability of a steady, vortical flow in a cylindrical container with one rotating endwall has been experimentally examined to gain insight into the process of vortex breakdowwn. The dynamics of the flow are governed by the Reynolds number (Re) and the aspect ratio of the cylinder. Re is given by Omega R(sup 2)/nu, where Omega is the speed of rotation of the endwall, R is the cylinder radius, and nu is the kinematic viscosity of the fluid filling the cylinder. The aspect ratio is H/R, where H is the height of the cylinder. Numerical simulation studies disagree whether or not the steady breakdown is stable beyond a critical Reynolds number, Re(sub c). Previous experimental researches have considered the steady and unsteady flows near Re(sub c), but have not explored the stability of the steady breakdown structures beyond this value. In this investigation, laser induced fluorescence was utilized to observe both steady and unsteady vortex breakdown at a fixed H/R of 2.5 with Re varying around Re(sub c). When the Re of a steady flow was slowly increased beyond Re(sub c), the breakdown structure remained steady even though unsteadiness was possible. In addition, a number of hysteresis events involving the oscillation periods of the unsteady flow were noted. The results show that both steady and unsteady vortex breakdown occur for a limited range of Re above Re(sub c). Also, with increasing Re, complex flow transformations take place that alter the period at which the unsteady flow oscillates.

Stevens, J. L.↗

A Power Conversion Concept for the Jupiter Icy Moons Orbiter

The Jupiter Icy Moons Orbiter (JIMO) mission is currently under study by the Office of Space Science under the Project Prometheus Program. JIMO is examining the use of Nuclear Electric Propulsion (NEP) to carry scientific payloads to three Jovian moons. A potential power system concept includes dual 100 kWe Brayton converters, a deployable pumped loop heat rejection subsystem, and a 400 Vac Power Management and Distribution (PMAD) bus. Many trades were performed in aniving at this candidate power system concept. System-level studies examined design and off-design operating modes, determined startup requirements, evaluated subsystem redundancy options, and quantified the mass and radiator area of reactor power systems from 20 to 200 kWe. In the Brayton converter subsystem, studies were performed to investigate converter packaging options, and assess the induced torque effects on spacecraft dynamics due to rotating machinery. In the heat rejection subsystem, design trades were conducted on heat transport approaches, material and fluid options, and deployed radiator geometries. In the PMAD subsystem, the overall electrical architecture was defined and trade studies examined distribution approaches, voltage levels, and cabling options.

Mason, Lee S.↗

Obliquity histories of Earth and Mars: Influence of inertial and dissipative core-mantle coupling

For both the Earth and Mars, secular variations in the angular separation of the spin axis from the orbit normal are suspected of driving major climatic changes. There is considerable interest in determining the amplitude and timing of these obliquity variations. If the orientation of the orbital plane were inertially fixed, and the planet were to act as a rigid body in it response to precessional torques, the spin axis would simply precess around the orbit at a fixed obliquity and at a uniform angular rate. The precession rate parameter depends on the principal moments of inertia and rotation rate of the perturbed body, and on the gravitational masses and semiminor axes of the perturbing bodies. For Mars, the precession rate is not well known, but probably lies in the interval 8 to 10 arcsec/year. Gravitational interactions between the planets lead to secular motions of the orbit planes. In the rigid body case, the spin axis still attempts to precess about the instantaneous orbit normal, but now the obliquity varies. The hydrostatic figure of a planet represents a compromise between gravitation, which attempts to attain spherical symmetry, and rotation, which prefers cylindrical symmetry. Due to their higher mean densities the cores of the Earth and Mars will be more nearly spherical than the outer layers of these planets. On short time scales it is appropriate to consider the core to be an inviscid fluid constrained to move with the ellipsoidal region bounded by the rigid mantle. The inertial coupling provided by this mechanism is effective whenever the ellipticicy of the container exceeds the ratio of precessional to rotational rates. If the mantle were actually rigid, this would be an extremely effective type of coupling. However, on sufficiently long time scales, the mantle will deform viscously and can accommodate the motions of the core fluid. A fundamentally different type of coupling is provided by electromagnetic or viscous torques. This type of coupling is likely to be most important on longer time scales. In each case, the mantle exerts an equal and opposite torque on the core.

Bills, Bruce G.↗

Venus' Free Obliquity

The predicted orientation of Venus' rotation axis relative to its orbit can be uniquely determined given knowledge of its J2 gravity coefficient and polar moment of inertia C if its free obliquity is fully damped. This assumption seems warranted given the dominant damping mechanism: turbulent fluid friction at a core mantle boundary. The determination of the free obliquity is examined.

Venus↗

Shape-anisotropy inverses the behavior of emergent vortices in active chiral fluids

Active colloidal fluids exhibit spontaneous emergence of correlated states, characterized by complex collective dynamics and self-organization. In geometrically confined systems, activity modulations trigger robust polar state reversals of a macroscopic vortex formed by colloidal rollers. Here, we reveal that the shape anisotropy of dielectric rollers transforms the system into a chiral active fluid comprised of rollers of arbitrary handedness. The chiral rollers self-assemble into multiple freestanding vortices with a spontaneously selected sense of rotation. We demonstrate that upon reactivation of the system after a complete cessation of activity beyond all relevant timescales, the vortices simultaneously restore their previous chiral states in striking contrast to the chiral state reversals demonstrated by spherical rollers. The analysis reveals that shape-anisotropy modifies the collective state memory of the ensemble. The findings provide insights into the emergence of complex collective behavior in chiral colloidal fluids governed by an interplay between shape anisotropy, chiral motion, and activity modulations.

Colloids↗

Design and fabrication of the Brayton rotating unit

The Brayton rotating unit (BRU), operating on a gas bearing system, has been designed, fabricated, and demonstrated for use in a closed Brayton cycle space power conversion system. The BRU uses a binary mixture of xenon and helium (molecular weight, 83.8) as the cycle working fluid and bearing lubricating medium and was designed to produce from 2.25 to 10.5 kw sub e of 1200 Hz three-phase electrical power. The single-shaft rotating assembly operates at a design speed of 36,000 rpm and comprises a radial single-stage compressor, a four-pole Rice alternator rotor, and a radial inflow turbine. Four units, a dynamic simulator and three component research packages, were supplied to the NASA Lewis Research Center for performance testing and further development.

Davis, J. E.↗

Rotating bubble membrane radiator for space applications

An advanced radiator concept for heat rejection in space is described which uses a two-phase working fluid to radiate waste heat. The development of advanced materials and the large surface area per mass makes the Bubble Membrane Radiator an attractive alternative to both conventional heat pipes and liquid droplet radiators for mid to high temperature applications. A system description, a discussion of design requirements, and a mass comparison with heat pipes and liquid droplet radiators are provided.

Webb, Brent J.↗

The heat Capacity of (sup 4)He Under Rotation Near T(wavelength)

The Lambda-transition in liquid helium has become a model for phase transitions. The transition with no counterflow between normal fluid and superfluid is rather well studied and understood in terms of renormalization group theory. Broken gauge symmetry gives rise to a new thermodynamic degree of freedom below the superfluid transition.

transition↗

Numerical Investigation of the Spark Discharge Process in a Crossflow

The present study numerically investigates the spark discharge process under crossflow conditions using a thermal equilibrium plasma solver that fully couples the electromagnetic physics and fluid dynamics in a computational framework. Numerical results are validated by the comparison with experimental data. Here, the spark discharge experiment is performed in a constant volume vessel using an inductive coil ignition system for automotive applications, and the evolution of the spark channel is measured using high-speed imaging. The crossflow in the gap between the spark-plug electrodes is generated by a rotating fan with two different fan speeds, and the flow velocity across the gap is characterized by particle image velocimetry (PIV) measurement. A computational fluid dynamics (CFD) solver is employed to simulate the crossflow and provide the flow field variables (velocity, pressure, temperature) to the plasma solver. The crossflow velocity predicted in the flow simulation agrees well with the PIV data in that the non-uniform velocity profiles at monitoring points are reproduced by the CFD code. With the crossflow initialization in the plasma solver, the simulated spark discharge process from the breakdown to spark discharge matches the experimental data, including the voltage and circuit waveforms and the high-speed images of the spark channel evolution. The stretch of spark channel captured by plasma simulations agrees with the measured data. The plasma simulation reveals that the mean temperature of the spark channel is maintained at 5000 K during the discharge phase, and the temperature varies along the spark channel so that the highest value is obtained at the spark root on the center electrode. Overall, the results presented in this paper are meant to provide valuable information about the properties of the plasma generated by the spark discharge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Explicit simulation of the Brownian rotation of arbitrary shaped aerosol particles using quaternions

The shape of an aerosol particle strongly influences its mass and momentum transfer cross-sections, charging properties, and other physical properties. Here, we present an explicit time-stepping procedure to simulate the rotational Brownian motion of arbitrary shaped aerosol particles by solving Euler’s equation of rotation. A Langevin formulation of the rotation equations is used, wherein Brownian motion due to thermal collisions between a particle and background gas molecules is represented using a stochastic fluctuating torque and fluid resistance is included as a drag torque. To avoid singularities associated with describing the orientation of a shape with Euler angles, we employ a quaternion formulation that leads to first-order stochastic differential equations to describe the evolution of the angular position and angular velocity of a rigid body. We perform all the rotational dynamics calculations in the body-fixed frame of reference attached to the rotating shape whose basis vectors are the normalized eigenvectors of the inertia tensor of the particle. Numerical solutions to rotation under torque-free conditions, damped rotation without Brownian motion, and stochastic rotation for arbitrary shapes are presented and discussed. The presented method enables time-resolved simulation of Brownian rotation for direct comparison with experimentally measured trajectories or statistical measures. The second order accuracy of the used time-stepping procedure places a severe restriction on the timestep that can be used for obtaining accurate results. Animations of presented simulations are included for visualizing rotational motion at various gas pressures. To aid implementation, MATLAB ® codes are also provided. Extension to include translation Brownian motion is straightforward.

Roy, Mrittika↗

A View of the Lunar Interior Through Lunar Laser Range Analysis

Laser ranges between observatories on the Earth and retroreflectors on the Moon started in 1969 and continue to the present. Recent range accuracies are 2 cm while earliest ranges are an C, order of magnitude less certain. Four retroreflectors are ranged: three located at the Apollo 11, 14, and 15 sites and one on the Lunakhod 2 rover. Accurate analysis of the range data determines a number of lunar science parameters. The lunar interior variables include a fluid core parameter. The Lunar Laser Ranging effort is reviewed elsewhere. Many parameters are detected through their influence on rotation. Also detected are solid-body tides and accurate selenocentric reflector locations. Determined through the rotation are moment-of-inertia differences, gravitational harmonics, potential Love number, and dissipation effects due to tides and molten core. The rotation of the Moon is not at its minimum energy state; some recently active process has caused free librations. The moment differences contributed to the recent improvement of the Moon's moment of inertia from the Lunar Prospector gravity field. The Love numbers provide bulk elastic properties. Future possibilities for measurement include oblateness of the core-mantle boundary and core moment. A study of dissipation signatures in the rotation determines tidal Q vs. frequency and concludes that the Moon has a molten core. At 1 month the tidal Q is 37 and at 1 yr it is 60. The core radius is < or = 352 km for Fe and < or = 374 km for the Fe-FeS eutectic. The core detection exceeds 3x its uncertainty. The spin of the core is not aligned with the spin of the mantle and torque arises from the velocity difference at the boundary. Yoder's turbulent boundary layer theory is used to compute the radii. The present heat generation from tides and core interaction is minor compared to radiogenic heating. The heating for ancient times is more interesting. Peale and Cassen investigated lunar tidal heating while the lunar orbit expanded due to tides on Earth. Their calculations predate the measurement of Q and should be multiplied by 3.45 to match the lunar-laser-determined Love number and monthly Q. Tidal-heating computations depend on how fast the lunar orbit evolved and whether the tidal dissipation is localized. Neither is known, but plausible assumptions lead to early central region temperature increases of several hundred degrees. Most of the energy is deposited early in the Moon's history. The turbulent boundary layer theory allows a prediction of energy dissipated at the core-mantle boundary during orbit evolution. Under the assumption that the properties of the early core are the same as at present, the energy dissipated by core-mantle interaction is about the same as for tidal dissipation, but it is deposited in a smaller volume. This source of energy is capable of promoting convection in an early fluid core and driving a dynamo. This is a transient phase with a duration depending on the rate of orbit evolution. Plausible assumptions lead to a duration of a few hundred million years. Thus, the remnant magnetization of many lunar rocks is compatible with a brief global magnetic field powered by dynamical energy dissipation. Analysis of the lunar laser ranges is providing information on lunar geophysics. Future data will improve accuracies of present solution parameters, and several more interior effects should be detectable.

Williams, J. G.↗

Symmetric baroclinic instability of a Hadley cell

A symmetric baroclinic instability is examined in terms of a Boussinesq fluid contained between two horizontal plates to determine the effects of the Ekman and thermal layers. Governing equations are written for a rotating reference frame, taking into account the Rossby, Ekman, and Prandtl numbers. Equations are defined for the perturbation functions, treated as an eigenvalue problem, and a numerical integration of the full eighth order differential system is performed by a shooting technique. An instability is found to occur in the Hadley cell containing both Ekman and thermal boundary layers when the Richardson number is close to unity. If the Prandtl number is fixed the critical Richardson number decreases with an increasing Ekman number until the Ekman number reaches a certain value, at which time the fluid is stable.

Antar, B. N.↗

Steady-state temperature distribution within a Brayton rotating unit operating in a power conversion system using helium-xenon gas

The Brayton rotating unit (BRU), consisting of a turbine, an alternator, and a compressor, was tested as part of a Brayton cycle power conversion system over a side range of steady state operating conditions. The working fluid in the system was a mixture of helium-xenon gases. Turbine inlet temperature was varied from 1200 to 1600 F, compressor inlet temperature from 60 to 120 F, compressor discharge pressure from 20 to 45 psia, rotative speed from 32 400 to 39 600 rpm, and alternator liquid-coolant flow rate from 0.01 to 0.27 pound per second. Test results indicated that the BRU internal temperatures were highly sensitive to alternator coolant flow below the design value of 0.12 pound per second but much less so at higher values. The armature winding temperature was not influenced significantly by turbine inlet temperature, but was sensitive, up to 20 F per kVA alternator output, to varying alternator output. When only the rotational speed was changed (+ or - 10% of rated value), the BRU internal temperatures varied directly with the speed.

Johnsen, R. L.↗