Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “oscillatory flows”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

Atmospheric Pressure Effects on Cryogenic Storage Tank Boil-Off

The Cryogenics Test Laboratory (CTL) at the Kennedy Space Center (KSC) routinely utilizes cryostat test hardware to evaluate comparative and absolute thermal conductivities of a wide array of insulation systems. The test method is based on measurement of the flow rate of gas evolved due to evaporative boil-off of a cryogenic liquid. The gas flow rate typically stabilizes after a period of a couple of hours to a couple of days, depending upon the test setup. The stable flow rate value is then used to calculate the thermal conductivity for the insulation system being tested. The latest set of identical cryostats, 1,000-L spherical tanks, exhibited different behavior. On a macro level, the flow rate did stabilize after a couple of days; however the stable flow rate was oscillatory with peak to peak amplitude of up to 25 percent of the nominal value. The period of the oscillation was consistently 12 hours. The source of the oscillation has been traced to variations in atmospheric pressure due to atmospheric tides similar to oceanic tides. This paper will present analysis of this phenomenon, including a calculation that explains why other cryostats are not affected by it.

Sass, J. P.↗

The nonlinear rheology of complex yield stress foods

Many foods have a yield stress that allows them to retain a desired shape at rest, but transition into a viscous fluid when being served or consumed. The determination of the yield stress of the food dictates how the foods are formed and packaged, how they are served, and how they are perceived when being eaten. Oscillatory shearing provides an ideal test protocol to map the rheology across a range of time and flow strength scales. We couple oscillatory shearing and an iterative recovery procedure to show that the yielding process is a continuous transition for two common yield stress foods. We show that unrecoverable processes from oscillatory tests are equivalent to the steady shear flow behavior. We show that this yielding behavior can be well approximated by a recently published model that treats yield stress materials as continuous viscoelastic fluids with a rate-dependent relaxation time and has parameters that can be obtained from the linear viscoelastic oscillatory frequency sweep and the steady shear flow curve.

Griebler, James J.↗

Convection in a Very Compressible Fluid: Comparison of Simulations With Experiments

The time profile (Delta)T(t) of the temperature difference, measured across a very compressible fluid layer of supercritical He-3 after the start of a heat flow, shows a damped oscillatory behavior before steady state convection is reached. The results for (Delta)T(t) obtained from numerical simulations and from laboratory experiments are compared over a temperature range where the compressibility varies by a factor of approx. = 40. First the steady-state convective heat current j(sup conv) as a function of the Rayleigh number R(alpha) is presented, and the agreement is found to be good. Second, the shape of the time profile and two characteristic times in the transient part of (Delta)T(t) from simulations and experiments are compared, namely: 1) t(sub osc), the oscillatory period and 2) t(sub p), the time of the first peak after starting the heat flow. These times, scaled by the diffusive time tau(sub D) versus R(alpha), are presented. The agreement is good for t(sup osc)/tau(sub D), where the results collapse on a single curve showing a powerlaw behavior. The simulation hence confirms the universal scaling behavior found experimentally. However for t(sub p)/tau(sub D), where the experimental data also collapse on a single curve, the simulation results show systematic departures from such a behavior. A possible reason for some of the disagreements, both in the time profile and in t(sub p) is discussed. In the Appendix a third characteristic time, t(sub m), between the first peak and the first oscillation minimum is plotted and a comparison between the results of experiments and simulations is made.

Meyer, H.↗

Numerical Studies of a Fluidic Diverter for Flow Control

The internal flow structure in a specific fluidic diverter is studied over a range from low subsonic to sonic inlet conditions by a time-dependent numerical analysis. The understanding will aid in the development of fluidic diverters with minimum pressure losses and advanced designs of flow control actuators. The velocity, temperature and pressure fields are calculated for subsonic conditions and the self-induced oscillatory behavior of the flow is successfully predicted. The results of our numerical studies have excellent agreement with our experimental measurements of oscillation frequencies. The acoustic speed in the gaseous medium is determined to be a key factor for up to sonic conditions in governing the mechanism of initiating the oscillations as well as determining its frequency. The feasibility of employing plasma actuation with a minimal perturbation level is demonstrated in steady-state calculations to also produce oscillation frequencies of our own choosing instead of being dependent on the fixed-geometry fluidic device.

Gokoglu, Suleyman A.↗

Control of Cavity Resonance Using Oscillatory Blowing

The near-zero net mass oscillatory blowing control of a subsonic cavity flow has been experimentally investigated. An actuator was designed and fabricated to provide both steady and oscillatory blowing over a range of blowing amplitudes and forcing frequencies. The blowing was applied just upstream of the cavity front Wall through interchangeable plate configurations These configurations enabled the effects of hole size, hole shape, and blowing angle to be examined. A significant finding is that in terms of the blowing amplitude, the near zero net mass oscillatory blowing is much more effective than steady blowing; momentum coefficients Lip two orders of magnitude smaller than those required for steady blowing are sufficient to accomplish the same control of cavity resonance. The detailed measurements obtained in the experiment include fluctuating pressure data within the cavity wall, and hot-wire measurements of the cavity shear layer. Spectral and wavelet analysis techniques are applied to understand the dynamics and mechanisms of the cavity flow with control. The oscillatory blowing, is effective in enhancing the mixing in the cavity shear layer and thus modifying the feedback loop associated with the cavity resonance. The nonlinear interactions in the cavity flow are no longer driven by the resonant cavity modes but by the forcing associated with the oscillatory blowing. The oscillatory blowing does not suppress the mode switching behavior of the cavity flow, but the amplitude modulation is reduced.

Scarfe, Alison Lamp↗

Convective effects in float-zone and Czochralski melts

The hydrodynamics of crystal-growth melts is a relatively new research area. Numerical modeling of these processes is necessary. The work discussed herein is in two parts: numerical simulations of the flow in a Czochralski melt, and also of that in a float zone. In addition, for the float-zone case, energy stability theory will be used to determine stability bounds for the onset of oscillatory thermo-capillary flow. Convective effects in crystal-growth melts arise from a variety of mechanisms. Temperature gradients both in the direction of gravity and normal to it give rise to convection due to buoyancy effects. Rotation of the crucible and/or crystal causes a forced convection which may augment or oppose the buoyancy-driven flow. Finally, thermo-capillary forces (due to the variation of surface tension with temperature) drive surface motions which in turn generate convection in the bulk fluid. All of these mechanisms are present in either Czochralski or float-zone growth. The objective of the Czochralski modeling is to develop an accurate numerical simulation of the flow in a Czochralski silicon melt and to investigate the effects of various parameters on the flow properties. Like some earlier investigations, the intent is to simulate the effects of buoyancy, forced and thermo-capillary convection, including unsteady effects. Unlike earlier work, the aim is to include the effects of a variable free surface and freezing interface and, possibly incorporate nonaxisymmetric effects.

Neitzel, G. P.↗

Hydrodynamic irreversibility of non-Brownian suspensions in highly confined duct flow

The irreversible behaviour of a highly confined non-Brownian suspension of spherical particles at low Reynolds number in a Newtonian fluid is studied experimentally and numerically. In the experiment, the suspension is confined in a thin rectangular channel that prevents complete particle overlap in the narrow dimension and is subjected to an oscillatory pressure-driven flow. In the small cross-sectional dimension, particles rapidly separate to the walls, whereas in the large dimension, features reminiscent of shear-induced migration in bulk suspensions are recovered. Furthermore, as a consequence of the channel geometry and the development and application of a single-camera particle tracking method, three-dimensional particle trajectories are obtained that allow us to directly associate relative particle proximity with the observed migration. Companion simulations of a steadily flowing suspension highly confined between parallel plates are conducted using the force coupling method, which also show rapid migration to the walls as well as other salient features observed in the experiment. While we consider relatively low volume fractions compared to most prior work in the area, we nevertheless observe significant and rapid migration, which we attribute to the high degree of confinement.

42 ENGINEERING↗

Concentration Measurements in Self-Excited, Momentum-Dominated Helium Jets

Flow structure of momentum-dominated pure helium jets discharged vertically into ambient air was investigated using high-speed rainbow schlieren deflectometry (RSD) technique. Effects of the operating parameters, i.e., Reynolds number (Re) and Richardson number (Ri), on the oscillatory behavior of the flow were examined over a range of experimental conditions. To seek the individual effect of these parameters, one of them was fixed and the other was varied with certain constraints. Measurements revealed highly periodic oscillations in the laminar region as well as high regularity in transition and turbulent regions. Maximum spectral power profiles at different axial locations indicated the oscillation amplitude increasing until the breakdown of the jet in the turbulent regime. The transition from the laminar to turbulent flow was also investigated. Fast Fourier transform analysis performed in the transition regime showed that the flow oscillates at a unique frequency, which was the same in the upstream laminar flow region. Measured deflection angle data were used in Abel inversion algorithm to construct the helium concentration fields. Instantaneous helium concentration contours revealed changes in the flow structure and evolution of vortical structures during an oscillation cycle. Temporal evolution plots of helium concentration at different axial location showed repeatable oscillations at all axial and radial locations up to the turbulent regime. A cross-correlation technique, applied to find the spatial displacements of the vortical structures, provided correlation coefficient peaks between consecutive schlieren images. Results show that the vortical structure convected and accelerated only in the axial direction.

Yildirim, Bekir Sedat↗

Developmental testing resulting in a simplified liquid oxygen check valve for the Space Shuttle Main Propulsion System

The coil spring in a Space Shuttle liquid oxygen check valve failed due to cyclic fatigue in September, 1991. The dual-flapper, swing check valve is used to prevent reverse flow to the Space Shuttle Main Engines. Upon inspection of the failed component, the spring tangs were found to be missing and heavy wear was observed on the inner diameter of the spring coils. The fracture surfaces revealed that the metal had been steadily worn away until a simple overload caused the final fracture. A series of flow tests using water and a water/gas mixture was conducted to identify the flow phenomenon which produced the cyclic wear. A Plexiglas outlet housing was utilized to view the flapper behavior under different flow conditions and to aid in high speed photography. The tests revealed that flow instabilities induced two oscillatory flapper responses: a rocking mode and a chattering mode. Initially, attempts were made to reduce the spring-flapper oscillations. However, the final solution to the problem was a springless configuration which satisfied the valve's design requirements and eliminated the oscillations. The springless design relied on the inherent ability of the reverse flow momentum to close the flappers.

Aber, Gregory S.↗

Effect of nozzle translation on a bounded vortex flow

A bounded vortex flow is generated by a nozzle that combines azimuthally tilted downward jets in a circular array and a central suction port. When the nozzle is directed toward a bottom surface (called the impingement surface), the flow develops a strong intake vortex spanning from the impingement surface to the suction port, which acts to generate strong shear stress on the impingement surface. This strong shear stress can aid in removing particles attached to the impingement surface or in enhancement of heat or mass transfer, where in typical applications, the nozzle would be translated along the impingement surface. This paper reports on a computational study of the effect of nozzle translation on the bounded vortex flow field. The nozzle translation induces a bending of the bounded vortex in the downstream flow direction, which becomes larger with increase in the nozzle translation speed. We identify an interesting oscillatory phenomenon by which the bounded vortex can resist the bending induced by the translational flow and rotate about a mean position downstream of the suction port. At sufficiently high translation speeds, the oscillatory process can no longer resist the downwash effect from nozzle translation and the vortex eventually detaches from the impingement surface, limiting the effectiveness of this flow device.

42 ENGINEERING↗

Experimental Investigation of Fuel Staging Effect on Modal Dynamics of Thermoacoustic Azimuthal Instabilities in a Multi-Nozzle Can Combustor

Abstract This paper analyzes the dynamics of unstable azimuthal thermoacoustic modes in a lean premixed combustor. Azimuthal modes can be decomposed into two counter rotating waves where they can either compete and potentially suppress one of them (spinning) or coexist (standing), depending on the operating conditions. This paper describes experimental results of the dynamical behaviors of these two waves. The experimental data were taken at different mass flow rates as well as different azimuthal fuel staging in a multi-nozzle can combustor. It is shown that at a low flow rate with uniform fuel distribution, the two waves have similar amplitudes, giving rise to a standing wave. However, the two amplitudes are slowly oscillating out of phase to each other, and the phase difference between the two waves also shows oscillatory behavior. For an intermediate flow rate, the dynamics show intermittency between standing and spinning waves, indicating that the system is bistable. In addition, the phase difference dramatically shifts when the mode switches between standing and spinning waves. For a high flow rate, the system stabilizes at a spinning wave most of the time. These experimental observations demonstrate that not only the amplitudes of two waves but also the phase difference plays an important role in the dynamics of azimuthal mode. For non-uniform azimuthal fuel staging, the modal dynamics exhibit only an oscillatory standing wave behavior regardless of the mass flow rate. Compared to the uniform fuel staging, however, the pressure magnitude is considerably reduced, which provides a potential strategy to mitigate and/or suppress the instabilities.

Kim, J.↗

Spin-Up Instability of a Levitated Molten Drop in MHD-Flow Transition to Turbulence

When an alternating magnetic field interacts with induced eddy currents in a conducting body, there will be a repulsive force between the body and the driving coil system generating the field. This repulsive force is the basis of electromagnetic levitation, which allows containerless processing of different materials. The eddy currents in the conducting body also generate Joule heating. Axial rotation of electromagnetically levitated objects is a common observation in levitation systems and often an undesirable side effect of such experiments on 1-g and -g. There have been recent efforts to use magnetic damping and suppress this tendency of body rotation. The first report of rotation in EML drops was attributed to a slight asymmetry of the shape and location of the levitation coils could change the axis and speed of rotation. Other theories of sample rotation include a frequency difference in the traveling electromagnetic waves and a phase difference in two different applied fields of the same frequency. All of these different mechanisms share the following characteristics: the torque is small, constant for constant field strength, and very weakly dependent on the sample's temperature and phase (solid or liquid). During experiments on the MSL-1 (First Microgravity Science Laboratory) mission of the Space Shuttle (STS-83 and STS-94, April and July 1997), a droplet of palladium-silicon alloy was electromagnetically levitated for viscosity measurements. For the non-deforming droplet, the resultant MHD flow inside the drop is inferred from motion of impurities on the surface. These observations indicate formation of a pair of co-rotating toroidal flow structures inside the spheroidal levitated drop that undergo secondary flow instabilities. As rise in the fluid temperature rises, the viscosity falls and the internal flow accelerates and becomes oscillatory; and beyond a point in the experiments, the surface impurities exhibit non-coherent chaotic motion signifying emergence of turbulence inside the drop. In this work, a background of these set of observations will be given followed by a presentation of our results on the digital particle tracking analysis that has been performed on a number of available videos. The analysis indicates that the levitated drop attains a constant rotational speed during the melting phase and formation of the co-rotating axi-symmetric laminar toroidal structures. However, the rate of axial rotation increases dramatically during the deformation of the toroidal structures anti their breakup into chaotic entities. This new data suggests an interaction between the flow inside the levitated molten drop and the driving coils in the experiments. Possible mechanisms for this interaction will be reviewed. The data will also be used to make an assessment of existing theories on droplet rotation.

Abedian, B.↗

An oscillatory kernel function method for lifting surfaces in mixed transonic flow

A study was conducted on the use of combined subsonic and supersonic linear theory to obtain economical and yet realistic solutions to unsteady transonic flow problems. With some modification, existing linear theory methods were combined into a single computer program. The method was applied to problems for which measured steady Mach number distributions and unsteady pressure distributions were available. By comparing theory and experiment, the transonic method showed a significant improvement over uniform flow methods. The results also indicated that more exact local Mach number effects and normal shock boundary conditions on the perturbation potential were needed. The validity of these improvements was demonstrated by application to steady flow.

Cunningham, A. M., Jr.↗

Multiple crossings of a very thin plasma sheet in the earth's magnetotail

High resolution magnetic field, thermal plasma, and energetic particle data measured on the IMP 8 spacecraft are used to demonstrate that a thin plasma sheet with rapidly flowing plasma undergoes violent oscillatory motion that causes the entire thickness of the plasma sheet to pass over the spacecraft in times which are sometimes as short as 10 s. Ions with 50 and 290 keV energy are investigated since their gyroradii are comparable to the plasma sheet thickness. As the spacecraft moves from the plasma sheet into the lobe, the energetic particle intensities decrease, but the 290 keV ions are relatively more numerous than 50 keV ions on field lines located deeper in the lobe.

Fairfield, D. H.↗

Unsteady convection in tin in a Bridgman configuration

When a quiescent fluid is heated sufficiently from below, steady convection will begin. Further heating will cause oscillatory and then turbulent flow. Theoretical results predict that the frequency of oscillation will depend on the square root of the Rayleigh number in the fluid. In the current work, liquid tin was heated from below for three aspect ratios, h/R = 3.4, 5.3, and 7.0. The experimental results are curve-fit for the square-root relation and also for a linear relation. The fit of the expression is evaluated using a correlation coefficient. An estimate for the first critical Rayleigh number (onset of steady convection) is obtained for both expressions. These values are compared to previous experimental results.

Knuteson, David J.↗

Practical aspects of spatially high accurate methods

The computational qualities of high order spatially accurate methods for the finite volume solution of the Euler equations are presented. Two dimensional essentially non-oscillatory (ENO), k-exact, and 'dimension by dimension' ENO reconstruction operators are discussed and compared in terms of reconstruction and solution accuracy, computational cost and oscillatory behavior in supersonic flows with shocks. Inherent steady state convergence difficulties are demonstrated for adaptive stencil algorithms. An exact solution to the heat equation is used to determine reconstruction error, and the computational intensity is reflected in operation counts. Standard MUSCL differencing is included for comparison. Numerical experiments presented include the Ringleb flow for numerical accuracy and a shock reflection problem. A vortex-shock interaction demonstrates the ability of the ENO scheme to excel in simulating unsteady high-frequency flow physics.

Godfrey, Andrew G.↗

Combustion wave instabilities in turbulent flows

Navier-Stokes solutions are used to calculate oscillatory components of pressure, velocity, and density, which in turn provide necessary data to compute energy growth factors to determine combustion instability. It is shown that wave instabilities are associated with changes in entropy and the space and time averages of oscillatory components of pressure, velocity and density, together with the mean flow field in the energy equation. Compressible laminar and turbulent flows and reacting flows with hydrogen/oxygen combustion are considered in this paper. The space shuttle main engine combustion/thrust chamber is used for illustration of the theory. The analysis shows that the increase of mean pressure and disturbances consistently results in the increase in instability. It is shown that adequate combustion instability analysis requires at least third order nonlinearity in energy growth or decay.

Yoon, W. S.↗

Mean Flow Augmented Acoustics in Rocket Systems

Oscillatory motion in solid rocket motors and liquid engines has long been a subject of concern. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process and gas dynamics. The customary approach to modeling acoustic waves inside a rocket chamber is to apply the classical inhomogeneous wave equation to the combustion gas. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while the acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The converging section of a rocket nozzle, where gradients in pressure, density, and velocity become large, is a notable region where this approach is not applicable. The expulsion of unsteady energy through the nozzle of a rocket is identified as the predominate source of acoustic damping for most rocket systems. An accurate model of the acoustic behavior within this region where acoustic modes are influenced by the presence of a steady mean flow is required for reliable stability predictions. Recently, an approach to address nozzle damping with mean flow effects was implemented by French [1]. This new approach extends the work originated by Sigman and Zinn [2] by solving the acoustic velocity potential equation (AVPE) formulated by perturbing the Euler equations [3]. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, (del squared)(psi) − (lambda/c)(exp 2)(psi) − M(dot)[M(dot)(del)(del(psi))] − 2(lambda(M/c) + (M(dot)del(M))(dot)del(psi)−2(lambda)(psi)[M(dot)del(1/c)]=0 (1) with M as the Mach vector, c as the speed of sound, and lambda as the complex eigenvalue. French apply the finite volume method to solve the steady flow field within the combustion chamber and nozzle with inviscid walls. The complex eigenvalues and eigenvector are determined with the use of the ARPACK eigensolver. The present study employs the COMSOL Multphysics framework to solve the coupled eigenvalue problem using the finite element approach. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluated the gas flow inside of a solid rocket motor using St. Robert's law modeling solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are used by the coefficient form of the mathematics module to determine the eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. Pertinent results from these analyses are the complex valued eigenvalue and eigenvectors. Comparisons are made to the French results to evaluate the modeling approach. A comparison of the French results with that of the present analysis is displayed in figures 1 and 2, respectively. The graphic shows the first tangential eigenvector's real (a) and imaginary (b) values.

Fischbach, Sean R.↗