Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid instabilities”

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 631 records · Page 35

Surface Tension and Fingering of Miscible Interfaces

Experiments on miscible, buoyantly unstable reaction-diffusion fronts and non-reacting displacement fronts in Hele-Shaw cells show a fingering-type instability whose wavelengths (lambda*) are consistent with an interfacial tension (sigma) at the front caused by the change in chemical composition, even though the solutions are miscible in all proportions. In conjunction with the Saffman-Taylor model, the relation sigma = K/tau, where tau is the interface thickness and K approximately equal 4 +/- 2 x 10(exp -6) dyne, enables prediction of our measured values of lambda* as well as results from prior experiments on miscible interfaces. These results indicate that even for miscible fluids, surface tension is generally a more significant factor than diffusion in interfacial stability and flow characteristics.

Abib, Mohammed↗

Electron acceleration by Alfven waves in the magnetosphere

The self-consistent electron kinetics of Alfven waves on the electron inertial scale is studied using a 2D hybrid-kinetic description. The ions follow a fluid description for Alfven waves at frequencies below the ion cyclotron frequency. The parallel electron dynamics is treated kinetically using particle-in-cell techniques. In this model, the electron plasma mode is eliminated, and only the physics of the Alfven waves is retained. At sufficiently large amplitudes, it is found that oblique Alfven waves break due to finite electron inertia in a cold plasma. The consequence of wave breaking is the formation of an electron beam which can be unstable to the beam-plasma instability. The electrons supporting the parallel current thermalize into a non-Maxwellian distribution with an energetic tail up to several keV, assuming a reasonable magnetospheric Alfven speed. In hot plasma simulations, electron trapping is the principal mechanism of electron acceleration. It is proposed that wave breaking or electron trapping of oblique Alfven waves at 1 R(E) can result in electron acceleration and may explain some observed auroral phenomena.

Hui, C.-H.↗

Inlet-Compressor Analysis using Coupled CFD Codes

Propulsion performance and operability are key factors in the development of a successful aircraft. For high-speed supersonic aircraft, mixed-compression inlets offer high performance but are susceptible to an instability referred to as unstart. An unstart occurs when a disturbance originating in the atmosphere or the engine causes the shock system to be expelled from the inlet. This event can have adverse effects on control of the aircraft, which is unacceptable for a passenger plane such as the high speed civil transport (HSCT). The ability to predict the transient response of such inlets to flow perturbations is, therefore, important to the proper design of the inlet and the control measures used to prevent unstart. Computational fluid dynamics (CFD) is having an increasing role in the analysis of individual propulsion components. Isolated inlet studies are relatively easy to perform, but a major uncertainty is the boundary condition used at the inlet exit to represent the engine - the so-called compressor face boundary condition. A one-dimensional (1-D) Euler inlet simulation showed that the predicted inlet unstart tolerance to free-stream pressure perturbations can vary by as much as a factor of about six, depending on the boundary condition used. Obviously a thorough understanding of dynamic interactions between inlets and compressors/fans is required to provide the proper boundary condition. To aid in this understanding and to help evaluate possible boundary conditions, an inlet-engine experiment was conducted at the University of Cincinnati. The interaction of acoustic pulses, generated in the inlet, with the engine were investigated. Because of the availability of experimental data for validation, it was decided to simulate the experiment using CFD. The philosophy here is that the inlet-engine system is best simulated by coupling (existing) specialized CFD component-codes. The objectives of this work were to aid in a better understanding of inlet-compressor interaction physics and the formulation of a more realistic compressor-face boundary condition for time-accurate CFD simulations of inlets. Previous simulations have used 1-D Euler engine simulations in conjunction with 1-D Euler and axisymmetric Euler inlet simulations. This effort is a first step toward CFD simulation of an entire engine by coupling multidimensional component codes.

Cole, Gary↗

Inlet-Compressor Analysis Using Coupled CFD Codes

Propulsion performance and operability are key factors in the development of a successful aircraft. For high-speed supersonic aircraft, mixed-compression inlets offer high performance but are susceptible to an instability referred to as unstart. An unstart occurs when a disturbance originating in the atmosphere or the engine causes the shock system to be expelled from the inlet. This event can have adverse effects on control of the aircraft, which is unacceptable for a passenger plane such as the high speed civil transport (HSCT). The ability to predict the transient response of such inlets to flow perturbations is, therefore, important to the proper design of the inlet and the control measures used to prevent unstart. Computational fluid dynamics (CFD) is having an increasing role in the analysis of individual propulsion components. isolated inlet studies are relatively easy to perform, but a major uncertainty is the boundary condition used at the inlet exit to represent the engine - the so-called compressor face boundary condition. A one-dimensional (I-D) Euler inlet simulation (ref. 1) showed that the predicted inlet unstart tolerance to free-stream pressure perturbations can vary by as much as a factor of about six, depending on the boundary condition used. Obviously, a thorough understanding of dynamic interactions between inlets and compressors/fans is required to provide the proper boundary condition.

Cole, Gary↗

Concept Demonstration of Dopant Selective Reactive Etching (DSRIE) in Silicon Carbide

Accurate quantification of combustor pressure dynamics for the primary purpose of experimental validation of computational fluid dynamics (CFD) codes requires the use of robust, reliable and sensitive pressure sensors that can resolve sub--pound-per-square-inch pressure levels in high temperature environments (i.e., combustor). The state of the art microfabricated piezoresistive silicon carbide (SiC) pressure sensors that we have developed are capable of operating reliably at 600 degrees Centigrade. This technology was used in support of the ARMD ISRP-ERA (NASA's Aeronautics Research Mission Directorate, Integrated System Research Project - Environmentally Responsible Aviation) program to quantify combustor thermoacoustic instabilities. The results showed that while the SiC pressure sensors survived the high temperature and measured instabilities, the diaphragm (force collector) was not thin enough to be sensitive in resolving sub-pound-per-square-inch pressures; 30 meters is the thinnest diaphragm achievable with conventional reactive ion etching (RIE) processes. Therefore, this precludes its use for sub-pound-per-square-inch pressure measurement with high fidelity. In order to effectively resolve sub-pound-per-square-inch pressures, a thinner more sensitive diaphragm (10 meters) is needed. To achieve this would require a new and innovative fabrication process technique.

static↗

Unsteady aerodynamic loading of delta wings for low and high angles of attack

Experimental and theoretical investigations dealing with unsteady flow phenomena are surveyed, with the emphasis on the pattern of vortices which originate from flow separation at sharp leading edges. It is concluded that these vortices exhibit quasi-steady behavior when the alpha-vibrations are such that bursting instability does not occur above the wing surface. A selection of test results from Jarrah (1988) is presented and discussed. For sharp-edged delta models at low speeds, the aerodynamic loads which are plotted quantify the role of parameters AR and K for three ranges of alpha-variation. An extremely approximate and empirical 'theory' is offered, with data on crossflow drag and burst location, to reproduce the behavior of these airloads up to 90 deg. Recent attempts to apply the more sophisticated tools of computational fluid dynamics to the combination of unsteadiness and very high alpha are shown to be deficient.

Ashley, H.↗

Towards an Entropy Stable Spectral Element Framework for Computational Fluid Dynamics

Entropy stable (SS) discontinuous spectral collocation formulations of any order are developed for the compressible Navier-Stokes equations on hexahedral elements. Recent progress on two complementary efforts is presented. The first effort is a generalization of previous SS spectral collocation work to extend the applicable set of points from tensor product, Legendre-Gauss-Lobatto (LGL) to tensor product Legendre-Gauss (LG) points. The LG and LGL point formulations are compared on a series of test problems. Although being more costly to implement, it is shown that the LG operators are significantly more accurate on comparable grids. Both the LGL and LG operators are of comparable efficiency and robustness, as is demonstrated using test problems for which conventional FEM techniques suffer instability. The second effort generalizes previous SS work to include the possibility of p-refinement at non-conforming interfaces. A generalization of existing entropy stability machinery is developed to accommodate the nuances of fully multi-dimensional summation-by-parts (SBP) operators. The entropy stability of the compressible Euler equations on non-conforming interfaces is demonstrated using the newly developed LG operators and multi-dimensional interface interpolation operators.

Carpenter, Mark H.↗

Benard and Marangoni convection in multiple liquid layers

Convective fluid dynamics of immiscible double and triple liquid layers are considered. First results on multilayer convective flow, in preparation for spaceflight experiment aboard IML-2 (International Microgravity Laboratory), are discussed. Convective flow in liquid layers with one or two horizontal interfaces with heat flow applied parallel to them is one of the systems investigated. The second system comprises two horizontally layered immiscible liquids heated from below and cooled from above, that is, heat flow orthogonal to the interface. In this system convection results due to the classical Benard instability.

Koster, Jean N.↗

Towards an "All Speed" Unstructured Upwind Scheme

In the authors previous studies [1], a time-accurate, upwind finite volume method (ETAU scheme) for computing compressible flows on unstructured grids was proposed. The scheme is second order accurate in space and time and yields high resolution in the presence of discontinuities. The scheme features a multidimensional limiter and multidimensional numerical dissipation. These help to stabilize the numerical process and to overcome the annoying pathological behaviors of upwind schemes. In the present paper, it will be further shown that such multidimensional treatments also lead to a nearly all-speed or Mach number insensitive upwind scheme. For flows at very high Mach number, e.g., 10, local numerical instabilities or the pathological behaviors are suppressed, while for flows at very low Mach number, e.g., 0.02, computation can be directly carried out without invoking preconditioning. For flows in different Mach number regimes, i.e., low, medium, and high Mach numbers, one only needs to adjust one or two parameters in the scheme. Several examples with low and high Mach numbers are demonstrated in this paper. Thus, the ETAU scheme is applicable to a broad spectrum of flow regimes ranging from high supersonic to low subsonic, appropriate for both CFD (computational fluid dynamics) and CAA (computational aeroacoustics).

Loh, Ching Y.↗

Hydrodynamics of coalescing binary neutron stars: Ellipsoidal treatment

We employ an approximate treatment of dissipative hydrodynamics in three dimensions to study the coalescence of binary neutron stars driven by the emission of gravitational waves. The stars are modeled as compressible ellipsoids obeying a polytropic equation of state; all internal fluid velocities are assumed to be linear functions of the coordinates. The hydrodynamics equations then reduce to a set of coupled ordinary differential equations for the evolution of the principal axes of the ellipsoids, the internal velocity parameters, and the binary orbital parameters. Gravitational radiation reaction and viscous dissipation are both incorporated. We set up exact initial binary equilibrium configurations and follow the transition from the quasi-static, secular decay of the orbit at large separation to the rapid dynamical evolution of the configurations just prior to contact. A hydrodynamical instability resulting from tidal interactions significantly accelerates the coalescence at small separation, leading to appreciable radial infall velocity and tidal lag angles near contact. This behavior is reflected in the gravitational waveforms and may be observable by gravitational wave detectors under construction. In cases where the neutron stars have spins which are not aligned with the orbital angular momentum, the spin-induced quadrupole moment can lead to precession of the orbital plane and therefore modulation of the gravitational wave amplitude even at large orbital radius. However, the amplitude of the modulation is small for typical neutron star binaries with spins much smaller than the orbital angular momentum.

Lai, Dong↗

Turbulent Transition in Electromagnetically Levitated Drops

Electromagnetic levitation (EML) is an important tool in materials research. Because a sample can be processed without contact with a container, experiments may be performed on high temperature, highly reactive, and undercooled liquid metals. Many of these experiments are affected by fluid flow in the sample, driven by the electromagnetic positioning force. Despite the importance of convection in these experiments, the transition to turbulence is not well understood in this system. However, we have observed a transition from laminar to turbulent flow in EML droplets in the course of microgravity experiments in TEMPUS on the Space Shuttle (STS-94). The transition occurs repeatably and over a narrow range of conditions. These experimental observations are compared with two competing theories about the transition to turbulence. Also, the results of a particle tracking study of the instabilities leading up to the transition to turbulence are presented.

Hyers, Robert W.↗

Findings from the Supersonic Qualification Program of the Mars Science Laboratory Parachute System

In 2012, the Mars Science Laboratory Mission (MSL) will deploy NASA's largest extra-terrestrial parachute, a technology integral to the safe landing of its advanced robotic explorer on the surface. The supersonic parachute system is a mortar deployed 21.5 m disk-gap-band (DGB) parachute, identical in geometric scaling to the Viking era DGB parachutes of the 1970's. The MSL parachute deployment conditions are Mach 2.3 at a dynamic pressure of 750 Pa. The Viking Balloon Launched Decelerator Test (BLDT) successfully demonstrated a maximum of 700 Pa at Mach 2.2 for a 16.1 m DGB parachute in its AV4 flight. All previous Mars deployments have derived their supersonic qualification from the Viking BLDT test series, preventing the need for full scale high altitude supersonic testing. The qualification programs for Mars Pathfinder, Mars Exploration Rover, and Phoenix Scout Missions were all limited to subsonic structural qualification, with supersonic performance and survivability bounded by the BLDT qualification. The MSL parachute, at the edge of the supersonic heritage deployment space and 33% larger than the Viking parachute, accepts a certain degree of risk without addressing the supersonic environment in which it will deploy. In addition, MSL will spend up to 10 seconds above Mach 1.5, an aerodynamic regime that is associated with a known parachute instability characterized by significant canopy projected area fluctuation and dynamic drag variation. This aerodynamic instability, referred to as "area oscillations" by the parachute community has drag performance, inflation stability, and structural implications, introducing risk to mission success if not quantified for the MSL parachute system. To minimize this risk and as an alternative to a prohibitively expensive high altitude test program, a multi-phase qualification program using computation simulation validated by subscale test was developed and implemented for MSL. The first phase consisted of 2% of fullscale supersonic wind tunnel testing of a rigid DGB parachute with entry-vehicle to validate two high fidelity computational fluid dynamics (CFD) tools. The computer codes utilized Large Eddy Simulation and Detached Eddy Simulation numerical approaches to accurately capture the turbulent wake of the entry vehicle and its coupling to the parachute bow-shock. The second phase was the development of fluid structure interaction (FSI) computational tools to predict parachute response to the supersonic flow field. The FSI development included the integration of the CFD from the first phase with a finite element structural model of the parachute membrane and cable elements. In this phase, a 4% of full-scale supersonic flexible parachute test program was conducted to provide validation data to the FSI code and an empirical dataset of the MSL parachute in a flight-like environment. The final phase is FSI simulations of the full-scale MSL parachute in a Mars type deployment. Findings from this program will be presented in terms of code development and validation, empirical findings from the supersonic testing, and drag performance during supersonic operation.

EDL↗

Some heat transfer and hydrodynamic problems associated with superconducting cables (SPTL)

To study some effects of thermogravitation on (CIIK-SPTL) systems, a heated tube experiment was set up at Krzhizhanovsky Power Engineering Institute Moscow, U.S.S.R. Heat transfer data were taken with fluid helium flowing through a 2.85 m, 19 mm diameter uniformly heated horizontal tube. Temperatures were measured on the top and bottom of the tube at six axial locations with three other circumferential measurements made at (X/L) =57. Typical temperature profiles show significant variations both axially and circumferentially. The data are grouped using reduced Nusselt number (NuR) and the bulk expansion parameter for each axial location. The average data for 0.26 less than or equal to X/L less than or equal to 0.76 follow a power law relation with the average expansion parameter. System instabilities are noted and discussed. Future work including heat transfer in coaxial cylinders is discussed.

Hendricks, R. C.↗

Progress and supercomputing in computational fluid dynamics; Proceedings of U.S.-Israel Workshop, Jerusalem, Israel, December 1984

Current developments and future trends in the application of supercomputers to computational fluid dynamics are discussed in reviews and reports. Topics examined include algorithm development for personal-size supercomputers, a multiblock three-dimensional Euler code for out-of-core and multiprocessor calculations, simulation of compressible inviscid and viscous flow, high-resolution solutions of the Euler equations for vortex flows, algorithms for the Navier-Stokes equations, and viscous-flow simulation by FEM and related techniques. Consideration is given to marching iterative methods for the parabolized and thin-layer Navier-Stokes equations, multigrid solutions to quasi-elliptic schemes, secondary instability of free shear flows, simulation of turbulent flow, and problems connected with weather prediction.

Murman, E. M.↗

Vortex rings

The vortex-ring problem in fluid mechanics is examined generally in terms of formation, the steady state, the duration of the rings, and vortex interactions. The formation is studied by examining the generation of laminar and turbulent vortex rings and their resulting structures with attention given to the three stages of laminar ring development. Inviscid dynamics is addressed to show how core dynamics affects overall ring motion, and laminar vortex structures are described in two dimensions. Viscous and inviscid structures are related in terms of 'leapfrogging', head-on collisions, and collisions with a no-slip wall. Linear instability theory is shown to successfully describe observational data, although late stages in the breakdown are not completely understood. This study of vortex rings has important implications for key aerodynamic issues including sound generation, transport and mixing, and vortex interactions.

Shariff, Karim↗

Laboratory and theoretical studies of baroclinic processes

An understanding is being developed for processes which may be important in the atmosphere, and the definition and analysis of baroclinic experiments utilizing the geophysical fluid flow cells (GFFC) apparatus in microgravity space flights. Included are studies using numerical codes, theoretical models, and terrestrial laboratory experiments. The numerical modeling is performed in three stages: calculation of steady axisymmetric flow, calculation of fastest-growing linear eigenmodes, and nonlinear effects (first, wave-mean flow interactions, then wave-wave interactions). The code can accommodate cylindrical, spherical, or channel geometry. It uses finite differences in the vertical and meridional directions, and is spectral in the azimuthal. The theoretical work was mostly in the area of effects of topography upon the baroclinic instability problem. The laboratory experiments are performed in a cylindrical annulus which has a temperture gradient imposed upon the lower surface and an approximately isothermal outer wall, with the upper and inner surfaces being nominally thermally insulating.

Miller, Timothy↗

Ground based preparation for microgravity growth of alloy semiconductors

Ground-based research conducted in order to prepare a microgravity space flight experiment is presented. The thermophysical properties of a PbSnTe alloy used for semiconductors are investigated, and furnace calibration and fluid-flow measurements are performed. The alloy has a zero energy crossing at approximately 40 percent SnTe in its band-gap vs composition diagram, which facilitates the design of long-wavelength IR detectors and lasers. The uniformity of devices made from this material depends on the ratio of PbTe and SnTe and requires the composition of the crystal growth to be closely controlled. The main obstacle to such control is the fact that liquid of this material is always solutally or thermally unstable, and, in a high-temperature gradient, the double convective instability cannot be made stable by balancing thermal and solutal expansion. In order to extend the science of crystal growth, the limits of suppression of convection have to be tested in low earth orbit.

Fripp, Archibald L.↗

Experimental determination of gap flow-conditioned forces at turbine stages and their effect on the running stability of simple rotors

Instabilities in turbine operation can be caused by forces which are produced in connection with motions involving the oil film in the bearings. An experimental investigation regarding the characteristics of such forces in the case of three typical steam turbine stages is conducted, taking into account the effect of various parameters. Supplementary kinetic tests are carried out to obtain an estimate of the flow forces which are proportional to the velocity. The measurements are based on the theoretical study of the damping characteristics of a vibrational model. A computational analysis of the effect of the measured fluid forces on the stability characteristics of simple rotor model is also conducted.

Wohlrab, R.↗