''Pressure-drop'' oscillations in forced convection flow with boiling.
Pressure drop oscillations in forced convection flow of Freon-11 near onset of boiling, noting flow rate for instability and limit cycles
SEARCH · Search NASA
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.
Pressure drop oscillations in forced convection flow of Freon-11 near onset of boiling, noting flow rate for instability and limit cycles
A test rig designed to generate heat exchanger pressure drop information under oscillating flow conditions is described. This oscillating flow rig is based on a variable stroke and variable frequency linear drive motor. A frequency capability of 120 hertz and a mean test pressure up to 15 mPA (2200 psi) allows for testing at flow conditions found in modern high specific power Stirling engines. An important design feature of this rig is that it utilizes a single close coupled dynamic pressure transducer to measure the pressure drop across the test sample. This eliminates instrumentation difficulties associated with the pressure sensing lines common to differential pressure transducers. Another feature of the rig is that it utilizes a single displacement piston. This allows for testing of different sample lengths and configurations without hardware modifications. All data acquisition and reduction for the rig is performed with a dedicated personal computer. Thus the overall system design efficiently integrates the testing and data reduction procedures. The design methodology and details of the test rig is described.
A test rig designed to generate heat exchanger pressure drop information under oscillating flow conditions is described. This oscillating flow rig is based on a variable stroke and variable frequency linear drive motor. A frequency capability of 120 hertz and a mean test pressure up to 15 mPA (2200 psi) allows for testing at flow conditions found in modern high specific power Stirling engines. An important design feature of this rig is that it utilizes a single close coupled dynamic pressure transducer to measure the pressure drop across the test sample. This eliminates instrumentation difficulties associated with the pressure sensing lines common to differential pressure transducers. Another feature of the rig is that it utilizes a single displacement piston. This allows for testing of different sample lengths and configurations without hardware modifications. All data acquisition and reduction for the rig is performed with a dedicated personal computer. Thus the overall system design efficiently integrates the testing and data reduction procedures. The design methodology and details of the test rig are described.
The development of systems such as acoustic levitation chambers will allow for the positioning and manipulation of material samples (drops) in a microgravity environment. This provides the capability for fundamental studies in droplet dynamics as well as containerless processing work. Such systems use acoustic radiation pressure forces to position or to further manipulate (e.g., oscillate) the sample. The primary objective was to determine the effect of a viscous acoustic field/tangential radiation pressure forcing on drop oscillations. To this end, the viscous acoustic field is determined. Modified (forced) hydrodynamic field equations which result from a consistent perturbation expansion scheme are solved. This is done in the separate cases of an unmodulated and a modulated acoustic field. The effect of the tangential radiation stress on the hydrodynamic field (drop oscillations) is found to manifest as a correction to the velocity field in a sublayer region near the drop/host interface. Moreover, the forcing due to the radiation pressure vector at the interface is modified by inclusion of tangential stresses.
Density wave, pressure drop and thermal flow oscillation modes in forced-convection boiling Freon 11
A series of experiments have been performed to investigate the oscillatory flow pressure and pressure drop through complex geometries. These experiments were conducted at the CSU-SLRE facility which is a horizontally opposed, two-piston, single-acting engine with a split crankshaft driving mechanism. Flow through a rectangular duct, with no insert (obstruction), was studied first. Then four different inserts were examined: Abrupt, Manifold, Diverging Short and Diverging Long. The inserts were mounted in the center of the rectangular duct to represent different type of geometries that could be encountered in Stirling machines. The pressure and pressure drop of the oscillating flow was studied for: 1) different inserts, 2) different phase angle between the two pistons of the engine (zero, 90 lead, 180, and 90 lag), and 3) for different piston frequencies (5, 10, 15, and 20 Hz). It was found that the pressure drop of the oscillatory flow increases with increasing Reynolds number. The pressure drop was shown to be mainly due to the gas inertia for the case of oscillatory flow through a rectangular duct with no insert. On the other hand, for the cases with different inserts into the rectangular duct, the pressure drop has three sources: inertia, friction, and local losses. The friction pressure drop is only a small fraction of the total pressure drop. It was also shown that the dimensionless pressure drop decreases with increasing kinetic Reynolds number.
The creeping flow of a dilute (0.025 wt%) monodisperse polystyrene/polystyrene Boger fluid through a 4:1:4 axisymmetric contraction/expansion is experimentally observed for a wide range of Deborah numbers. Pressure drop measurements across the orifice plate show a large extra pressure drop that increases monotonically with Deborah number above the value observed for a similar Newtonian fluid at the same flow rate. This enhancement in the dimensionless pressure drop is not associated with the onset of a flow instability, yet it is not predicted by existing steady-state or transient numerical computations with simple dumbbell models. It is conjectured that this extra pressure drop is the result of an additional dissipative contribution to the polymeric stress arising from a stress-conformation hysteresis in the strong non-homogeneous extensional flow near the contraction plane. Such a hysteresis has been independently measured and computed in recent studies of homogeneous transient uniaxial stretching of PS/PS Boger fluids. Flow visualization and velocity field measurements using digital particle image velocimetry (DPIV) show large upstream growth of the corner vortex with increasing Deborah number. At large Deborah numbers, the onset of an elastic instability is observed, first locally as small amplitude fluctuations in the pressure measurements, and then globally as an azimuthal precessing of the upstream corner vortex accompanied by periodic oscillations in the pressure drop across the orifice.
A high efficiency Stirling Radioisotope Generator (SRG) is being developed for possible use in long-duration space science missions. NASA s advanced technology goals for next generation Stirling convertors include increasing the Carnot efficiency and percent of Carnot efficiency. To help achieve these goals, a multi-dimensional Computational Fluid Dynamics (CFD) code is being developed to numerically model unsteady fluid flow and heat transfer phenomena of the oscillating working gas inside Stirling convertors. In the absence of transient pressure drop data for the zero mean oscillating multi-dimensional flows present in the Technology Demonstration Convertors on test at NASA Glenn Research Center, unidirectional flow pressure drop test data is used to compare against 2D and 3D computational solutions. This study focuses on tracking pressure drop and mass flow rate data for unidirectional flow though a Stirling heater head using a commercial CFD code (CFD-ACE). The commercial CFD code uses a porous-media model which is dependent on permeability and the inertial coefficient present in the linear and nonlinear terms of the Darcy-Forchheimer equation. Permeability and inertial coefficient were calculated from unidirectional flow test data. CFD simulations of the unidirectional flow test were validated using the porous-media model input parameters which increased simulation accuracy by 14 percent on average.
The dynamic behavior of the vapor flow in heat pipes is investigated at startup and during operational transients. The vapor is modeled as two-dimensional, compressible viscous flow in an enclosure with inflow and outflow boundary conditions. For steady-state and operating transients, the SIMPLER method is used. In this method a control volume approach is employed on a staggered grid which makes the scheme very stable. It is shown that for relatively low input heat fluxes the compressibility of the vapor flow is low and the SIMPLER scheme is suitable for the study of transient vapor dynamics. When the input heat flux is high or the process under a startup operation starts at very low pressures and temperatures, the vapor is highly compressible and a shock wave is created in the evaporator. It is shown that for a wide range of input heat fluxes, the standard methods, including the SIMPLER scheme, are not suitable. A nonlinear filtering technique, along with the centered difference scheme, are then used for shock capturing as well as for the solution of the cell Reynolds-number problem. For high heat flux, the startup transient phase involves multiple shock reflections in the evaporator region. Each shock reflection causes a significant increase in the local pressure and a large pressure drop along the heat pipe. Furthermore, shock reflections cause flow reversal in the evaporation region and flow circulations in the adiabatic region. The maximum and maximum-averaged pressure drops in different sections of the heat pipe oscillate periodically with time because of multiple shock reflections. The pressure drop converges to a constant value at steady state. However, it is significantly higher than its steady-state value at the initiation of the startup transient. The time for the vapor core to reach steady-state condition depends on the input heat flux, the heat pipe geometry, the working fluid, and the condenser conditions. However, the vapor transient time, for an Na-filled heat pipe is on the order of seconds. Depending on the time constant for the overall system, the vapor transient time may be very short. Therefore, the vapor core may be assumed to be quasi-steady in the transient analysis of a heat pipe operation.
Four configurations of piloted flameholders were tested. The range of flame stabilization, flame propagation, pressure oscillation during ignition, and pressure drop of the configurations were determined. Some tests showed a very strong effect of inlet flow velocity profile and flameholder geometry on flame stabilization. These tests led to the following conclusions. (1) The use of a piloted flameholder in the turbofan augmentor may minimize the peak pressure rise during ignition. At the present experimental conditions, delta P/P asterisk over 2 is less than 10 percent; therefore, the use of a piloted flameholder is a good method to realize soft ignition. (2) The geometry of the piloted flameholder and the amount of fuel injected into the flameholder have a strong effect on the pressure oscillation during ignition of the fuel-air mixture in the secondary zone. (3) Compared with the V-gutter flameholder with holes in its wall, the V-gutter flameholder without holes not only has advantages such as simple structure and good rigidity but offers a wide combustion stability limit and a high capability of igniting the fuel-air mixture of the secondary zone.
The development of acoustic levitation systems has provided a technology with which to undertake droplet studies as well as do containerless processing experiments in a microgravity environment. Acoustic levitation chambers utilize radiation pressure forces to position/manipulate the drop. Oscillations can be induced via frequency modulation of the acoustic wave, with the modulated acoustic radiation vector acting as the driving force. To account for tangential as well as radial forcing, it is necessary that the viscous effects be included in the acoustic field. The method of composite expansions is employed in the determination of the acoustic field with viscous effects.
Similarity parameters for characterizing the effect of flow oscillation on wall shear stress, viscous dissipation, pressure drop and heat transfer rates are proposed. They are based on physical agruments and are derived by normalizing the governing equations. The literature on oscillating duct flows, regenerator and porous media flows is surveyed. The operating characteristics of the heat exchanger of eleven Stirling engines are discribed in terms of the similarity parameters. Previous experimental and analytical results are discussed in terms of these parameters and used to estimate the nature of the oscillating flow under engine operating conditions. The operating points for many of the modern Stirling engines are in or near the laminar to turbulent transition region. In several engines, working fluid does not pass entirely through heat exchangers during a cycle. Questions that need to be addressed by further research are identified.
Although stable operating regimes for cryogenic coaxial injectors have been empirically determined, there is no knowledge of the spray characteristics corresponding to stable operation, or the physical processes which produce the atomization patterns that result in stable or unstable operation. The current engineering method for determining the stable operating regime of a cryogenic coaxial injector is the 'hydrogen temperature ramping' method, however there is no definitive knowledge of whether the hydrogen temperature influences the chamber stability by decreasing the injected gas velocity, by affecting a recirculation region at the base of the LOX Post, or by changing the pressure drop across the injector, allowing chamber pressure oscillations to couple to the fuel feed system. Results for the injector response from a linearized lumped-element model are presented as a function of temperature and frequency. LDV measurements in the recess region at the base of the LOX post show reverse flow indicative of a recirculation region. Finally, Phase Doppler Particle Analyzer (PDPA) measurements of droplet size and velocity distributions are discussed for a coaxial injector element similar to the SSME preburner element operating with water and air at atmospheric pressure and liquid and gaseous nitrogen at 20 bars.
The design and the performance of an electrostatic-acoustic hybrid system and its application to a charge drop rotation experiment are presented. This system can levitate a charged drop electrostatically and induce drop rotation or oscillation by imposing an acoustic torque or an oscillating acoustic pressure. Using this system, the equilibrium shapes and stability of a rotating charged drop were experimentally investigated. A 3 mm size water drop was rotated as a rigid body and its gyrostatic equilibrium shapes were observed. Families of axisymmetric shapes, two-lobed shapes, and eventual fissioning have been observed. With the assumption of 'effective surface tension' in which the surface charge simply modified the surface tension of neutral liquid, the results agree exceptionally well with the Brown and Scriven's (1980) prediction for uncharged drops.
We have developed a calculational model that treats all the components of an orifice pulse tube cooler. We base our analysis on 1-dimensional thermodynamic equations for the regenerator and we assume that all mass flows, pressure oscillations and temperature oscillations are small and sinusoidal. Non-linear pressure drop effects are included in the regenerator to account for finite pressure amplitude effects. The resulting mass flows and pressures are matched at the boundaries with the other components of the cooler: compressor, aftercooler, cold heat exchanger, pulse tube, hot heat exchanger, orifice and reservoir. The results of the calculation are oscillating pressures, mass flows and enthalpy flows in the main components of the cooler. By comparing with the calculations of other available models, we show that our model is very similar to REGEN 3 from NIST and DeltaE from Los Alamos National Lab. Our model is much easier to use than other available models because of its simple graphical interface and the fact that no guesses are required for the operating pressures or mass flows. In addition, the model only requires a few minutes of running time allowing many parameters to be optimized in a reasonable time. A version of the model is available for use over the World Wide Web at http://irtek.arc.nasa.gov. Future enhancements include adding a bypass orifice and including second order terms in steady mass streaming and steady heat transfer. A two-dimensional anelastic approximation of the fluid equations will be used as the basis for the latter analysis. Preliminary results are given in dimensionless numbers appropriate for oscillating compressible flows. The model shows how transverse heat transfer reduces enthalpy flow, particularly for small pulse tubes. The model also clearly shows mass recirculation in the open tube on the order of the tube length. They result from the higher order Reynolds stresses. An interesting result of the linearized approach is that the steady mass streaming does not affect the enthalpy flow at second order. The major effect of recirculating mass streaming is to increase transverse temperature gradients, which leads to higher entropy production and reduced efficiency.
Large amplitude oscillations of drops and bubbles immersed in an immiscible liquid host have been investigated using ultrasonic radiation pressure techniques.
The NASA Marshall Space Flight Center's electrostatic levitation (ESL) laboratory has been recently upgraded with an oxygen partial pressure controller. This system allows the oxygen partial pressure within the vacuum chamber to be measured and controlled, theoretically in the range from 10‐36 to 100 bar. The oxygen control system installed in the ESL laboratory's main chamber consists of an oxygen sensor, oxygen pump, and a control unit. The sensor is a potentiometric device that determines the difference in oxygen activity in two gas compartments (inside the chamber and the air outside of the chamber) separated by an electrolyte, which is yttria‐stabilized zirconia. The pump utilizes coulometric titration to either add or remove oxygen. The system is controlled by a desktop control unit, which can also be accessed via a computer. The controller performs temperature control for the sensor and pump, PID-based current loop, and a control algorithm. Oxygen partial pressure has been shown to play a significant role in the surface tension of liquid metals. Oxide films or dissolved oxygen may lead to significant changes in surface tension. The effects of oxygen partial pressure on the surface tension of undercooled liquid nickel will be analyzed, and the results will be presented. The surface tension will be measured at several different oxygen partial pressures while the sample is undercooled. Surface tension will be measured using the oscillating drop method. While undercooled, each sample will be oscillated several times consecutively to investigate how the surface tension behaves with time while at a particular oxygen partial pressure.
A transient model of a hybrid motor was formulated to study the cause and elimination of nonacoustic combustion instability. The transient model was used to simulate four key tests out of a series of seventeen hybrid motor tests conducted by Thiokol, Rocketdyne, and Martin Marietta at NASA Marshall Space Flight Center (MSFC). These tests were performed under the Hybrid Propulsion Technology for Launch Vehicle Boosters (HPTLVB) program. The first test resulted in stable combustion. The second test resulted in large-amplitude, 6.5-Hz chamber pressure oscillations that gradually damped away by the end of the test. The third test resulted in large-amplitude, 7.5-Hz chamber pressure oscillations that were sustained throughout the test. The seventh test resulted in elimination of combustion instability with the installation of an orifice immediately upstream of the injector. Formulation and implementation of the model are the scope of this presentation. The current model is an independent continuation of modeling presented previously by joint Thiokol-Rocketdyne collaborators Boardman, Hawkins, Wassom. and Claflin. The previous model simulated an unstable independent research and development (IR&D) hybrid motor test performed by Thiokol. There was very good agreement between the model and test data. Like the previous model, the current model was developed using Matrix-x simulation software. However, tests performed at MSFC under the HPTLVB program were actually simulated. ln the current model, the hybrid motor, consisting of the liquid oxygen (lox) injector, the multiport solid fuel grain, and nozzle, was simulated. The lox feedsystem, consisting of the tank, venturi. valve, and feed lines, was also simulated in the model. All components of the hybrid motor and lox feedsystem are treated by a lumped-parameter approach. Agreement between the results of the transient model and actual test data was very good. This agreement between simulated and actual test data indicated that the combustion instability in the hybrid motor was due to two causes: 1. a lox feed system of insufficient stiffness, and 2. a lox injector with an impedance (it pressure drop that was too low to provide damping against the feed system oscillations. Also, it was discovered that testing with a new grain of solid fuel sustained the combustion instability. However, testing with a used grain of solid fuel caused the combustion instability to gradually decay.