CFD analysis of high frequency miniature pulse tube refrigerators for space applications with thermal non-equilibrium model

High frequency, miniature, pulse tube cryocoolers are extensively used in space applications because of their simplicity. Parametric studies of inertance type pulse tube cooler are performed with different length-to-diameter ratios of the pulse tube with the help of the FLUENT (R) package. The local thermal non-equilibrium of the gas and the matrix is taken into account for the modeling of porous zones, in addition to the wall thickness of the components. Dynamic characteristics and the actual mechanism of energy transfer in pulse are examined with the help of the pulse tube wall time constant. The heat interaction between pulse tube wall and the oscillating gas, leading to surface heat pumping, is quantified. The axial heat conduction is found to reduce the performance of the pulse tube refrigerator. The thermal non-equilibrium predicts a higher cold heat exchanger temperature compared to thermal equilibrium. The pressure drop through the porous medium has a strong non-linear effect due to the dominating influence of Forchheimer term over that of the linear Darcy term at high operating frequencies. The phase angle relationships among the pressure, temperature and the mass flow rate in the porous zones are also important in determining the performance of pulse tuberefrigerator.

Effect of Bottom Bubbling Conditions on Surface Reaction Rate in Oxygen-Water System

In secondary steelmaking, the enhancement of the reaction rate in the low carbon period during the decarburization of steel is considered the most effective method to produce ultralow carbon steel. In a previous study, it was revealed that the surface reaction is dominant during the final stage of the actual refining process. In order to improve the surface reaction rate, it is necessary to enlarge the reaction region, which is usually achieved by increasing the plume eye area. In this study, water model experiments were carried out to estimate the influence of bottom stirring conditions on the gas-liquid reaction rate; for this purpose, the deoxidation rate during the bottom bubbling process was measured. Five types of nozzle configurations were used to study the effect of the plume eye area on the reaction rate at various gas flow rates. The results reveal that the surface reaction rate is influenced by the gas flow rate and the plume eye area. An empirical correlation was developed for the reaction rate and the plume eye area. This correlation was applied to estimate the gas-liquid reaction rate mat the bath surface.

Studies in bubble formation — III

Bubble formation from single horizontal orifices submerged in Newtonian liquids has been investigated for such chamber volumes that both the pressure inside the chamber and flow rate into the bubble are time dependent. The data collected show that under these conditions the bubble volume decreases exponentially with increase in orifice submergence. The equations for the generalized two stage model of bubble formation, taking the variation of gas flow rate with time into account, have been derived. These equations reduce to the cases of constant gas flow rate and constant pressure when adequate constraints are imposed. The results obtained under intermediate conditions have been quantitatively explained on the basis of these equations.

Observations on transition in plane bubble plumes

Flow visualization studies of plane laminar bubble plumes have been conducted to yield quantitative data on transition height, wavelength and wave velocity of the most unstable disturbance leading to transition. These are believed to be the first results of this kind. Most earlier studies are restricted to turbulent bubble plumes. In the present study, the bubble plumes were generated by electrolysis of water and hence very fine control over bubble size distribution and gas flow rate was possible to enable studies with laminar bubble plumes. Present observations show that (a) the dominant mode of instability in plane bubble plumes is the sinuous mode, (b) transition height and wavelength are related linearly with the proportionality constant being about 4, (c) wave velocity is about 40 % of the mean plume velocity, and (d) normalized transition height data correlate very well with a source Grashof number. Some agreement and some differences in transition characteristics of bubble plumes have been observed compared to those for similar single-phase flows.

Studies on helium flow in experimental loops with different impedances operated with fountain effect pumps

An experimental flow loop with He II flow driven by fountain effect pumps (FEPs) is studied with respect to operation at different flow impedances and with thermal loads applied at different positions. The measured values of temperature, flow rate and pressure drop are compared with calculations resulting from a simplified model which assumes ideal performance of the porous plug and of the heat exchangers and which does not take into account Gorter-Mellink (GM) conduction. The main features of the loop are shown to be well described by this model. Refined calculations with a more complex model, including GM conduction of the He II, are only required for predicting the temperature distribution in some discrete regions of the loop.

Low cycle fatigue behaviour of a nitrogen modified 316L stainless steel at 823 K

Strain controlled low cycle fatigue tests on solution annealed nitrogen modified 316L stainless steel have been conducted in air at 823 K to ascertain the influence of strain rate and strain amplitude. Effect of strain rate was examined from 3x10(-5) s(-1) to 3 x 10(-2) at a fixed strain amplitude of +/- 0.6%. The influence of strain amplitude was evaluated between +/- 0.25 % and +/- 1.0% at a constant strain rate of 3x10(-3) s(-1). The cyclic stress response at all testing conditions is characterized by an initial hardening followed by saturation. Serrated flow, a characteristic feature of dynamic strain ageing (DSA) was seen at strain rates lower than 3x10(-3) s(-1). Fatigue life was found to decrease with decrease in strain rate. The reduction in fatigue resistance is attributed mainly to the detrimental effects associated with DSA.

Work fluctuations in an elastic dumbbell model of polymers in planar elongational flow

We use a path-integral approach to calculate the distribution P(w, t) of the fluctuations in the work W at time t of a polymer molecule (modeled as an elastic dumbbell in a viscous solvent) that is acted on by an elongational flow field having a flow rate (gamma) over dot. We find that P(w, t) is non-Gaussian and that, at long times, the ratio P(w, t)/ P (-w, t) is equal to expw/(k(B)T)], independent of (gamma) over dot. On the basis of this finding, we suggest that polymers in elongational flows satisfy a fluctuation theorem.

Flow of liquid/liquid dispersions in a Hele-Shaw cell

Flow of liquid/liquid dispersions have been investigated in a Hele-Shaw cell which contained a thin disk held between two parallel plates. This device offers a well defined flow field and also permits visual observation of the dispersed drop movement. The dispersed drops coalesce with the disk for the systems where the dispersed phase wets the disk surface. The dispersed phase accumulate at the downstream end of the disk and they detach from there as blobs. Through an accurate measurement of accumulated dispersed phase volume, the coalescence rate was determined. The coalescence efficiency in the Hele Shaw cell is determined by dividing the coalescence hate by the undisturbed flow rate of the dispersed phase through an area equal to the projected area of the disk on a plane normal to the flow direction. The coalescence efficiency first increases and then decreases with the flow rate of dispersion. The coalescence rate/disk dimensions increases with the decrease in the disk dimensions. The rate of coalescence increases with the increase in the dispersed drop diameter and it decreases with the increase in the continuous phase viscosity. The presence of surfactants reduces the coalescence rate. All these results are quantitatively explained through a model, which takes into account several important features like various mechanism of drainage, the roles of dispersion and continuous phase viscosities, and the drop deformation.

Profile-Based Load Scheduling in Wireless Energy Harvesting Sensors for Data Rate Maximization

In this paper, power management algorithms for energy harvesting sensors (EHS) that operate purely based on energy harvested from the environment are proposed. To maintain energy neutrality, EHS nodes schedule their utilization of the harvested power so as to save/draw energy into/from an inefficient battery during peak/low energy harvesting periods, respectively. Under this constraint, one of the key system design goals is to transmit as much data as possible given the energy harvesting profile. For implementational simplicity, it is assumed that the EHS transmits at a constant data rate with power control, when the channel is sufficiently good. By converting the data rate maximization problem into a convex optimization problem, the optimal load scheduling (power management) algorithm that maximizes the average data rate subject to energy neutrality is derived. Also, the energy storage requirements on the battery for implementing the proposed algorithm are calculated. Further, robust schemes that account for the insufficiency of battery storage capacity, or errors in the prediction of the harvested power are proposed. The superior performance of the proposed algorithms over conventional scheduling schemes are demonstrated through computations using numerical data from solar energy harvesting databases.

Particle dynamics in a turbulent particle–gas suspension at high Stokes number. Part 1. Velocity and acceleration distributions

The effect of fluid velocity fluctuations on the dynamics of the particles in a turbulent gas–solid suspension is analysed in the low-Reynolds-number and high Stokes number limits, where the particle relaxation time is long compared with the correlation time for the fluid velocity fluctuations, and the drag force on the particles due to the fluid can be expressed by the modified Stokes law. The direct numerical simulation procedure is used for solving the Navier–Stokes equations for the fluid, the particles are modelled as hard spheres which undergo elastic collisions and a one-way coupling algorithm is used where the force exerted by the fluid on the particles is incorporated, but not the reverse force exerted by the particles on the fluid. The particle mean and root-mean-square (RMS) fluctuating velocities, as well as the probability distribution function for the particle velocity fluctuations and the distribution of acceleration of the particles in the central region of the Couette (where the velocity profile is linear and the RMS velocities are nearly constant), are examined. It is found that the distribution of particle velocities is very different from a Gaussian, especially in the spanwise and wall-normal directions. However, the distribution of the acceleration fluctuation on the particles is found to be close to a Gaussian, though the distribution is highly anisotropic and there is a correlation between the fluctuations in the flow and gradient directions. The non-Gaussian nature of the particle velocity fluctuations is found to be due to inter-particle collisions induced by the large particle velocity fluctuations in the flow direction. It is also found that the acceleration distribution on the particles is in very good agreement with the distribution that is calculated from the velocity fluctuations in the fluid, using the Stokes drag law, indicating that there is very little correlation between the fluid velocity fluctuations and the particle velocity fluctuations in the presence of one-way coupling. All of these results indicate that the effect of the turbulent fluid velocity fluctuations can be accurately represented by an anisotropic Gaussian white noise.

Microwave ECR Plasma Assisted MOCVD of Y(2)O(3) Thin Films Using Y(tod)(3) Precursor and Their Characterization

Yttrium oxide (Y(2)O(3)) thin films were deposited by microwave electron cyclotron resonance (ECR) plasma assisted metal organic chemical vapour deposition (MOCVD) process using indigenously developed metal organic precursors Yttrium 2,7,7-trimethyl-3,5-octanedionates, commonly known as Y(tod)(3) which were synthesized by an ultrasound method. A series of thin films were deposited by varying the oxygen flow rate from 1-9 sccm, keeping all other parameters constant. The deposited coatings were characterized by X-ray photoelectron spectroscopy, glancing angle X-ray diffraction and infrared spectroscopy. Thickness and roughness for the films were measured by stylus profilometry. Optical properties of the coatings were studied by the spectroscopic ellipsometry. Hardness and elastic modulus of the films were measured by nanoindentation technique. Being that microwave ECR CVD process is operating-pressure-sensitive, optimum oxygen activity is very essential for a fixed flow rate of precursor, in order to get a single phase cubic yttrium oxide in the films. To the best of our knowledge, this is the first effort that describes the use of Y(tod)(3) precursor for deposition of Y(2)O(3) films using plasma assisted CVD process.

Single Particle and Packed Bed Combustion in Modern Gasifier Stoves—Density Effects

This article is concerned with a study of an unusual effect due to density of biomass pellets in modern stoves based on close-coupled gasification-combustion process. The two processes, namely, flaming with volatiles and glowing of the char show different effects. The mass flux of the fuel bears a constant ratio with the air flow rate of gasification during the flaming process and is independent of particle density; char glowing process shows a distinct effect of density. The bed temperatures also have similar features: during flaming, they are identical, but distinct in the char burn (gasification) regime. For the cases, wood char and pellet char, the densities are 350, 990 kg/m(3), and the burn rates are 2.5 and 3.5 g/min with the bed temperatures being 1380 and 1502 K, respectively. A number of experiments on practical stoves showed wood char combustion rates of 2.5 +/- 0.5 g/min and pellet char burn rates of 3.5 +/- 0.5 g/min. In pursuit of the resolution of the differences, experimental data on single particle combustion for forced convection and ambient temperatures effects have been obtained. Single particle char combustion rate with air show a near-d(2) law and surface and core temperatures are identical for both wood and pellet char. A model based on diffusion controlled heat release-radiation-convection balance is set up. Explanation of the observed results needs to include the ash build-up over the char. This model is then used to explain observed behavior in the packed bed; the different packing densities of the biomass chars leading to different heat release rates per unit bed volume are deduced as the cause of the differences in burn rate and bed temperatures.

Filling Characteristics for an Activated Carbon Based Adsorbed Natural Gas Storage System

The storage capacity of an activated carbon bed is studied using a 2D transport model with constant inlet flow conditions. The predicted filling times and variation in bed pressure and temperature are in good agreement with experimental observations obtained using a 1.82 L prototype ANG storage cylinder. Storage efficiencies based on the maximum achievable V/V (volume of gas/volume of container) and filling times are used to quantify the performance of the charging process. For the high permeability beds used in the experiments, storage efficiencies are controlled by the rate of heat removal. Filling times, defined as the time at which the bed pressure reaches 3.5 MPa, range from 120 to 3.4 min for inlet flow rates of 1.0 L min(-1) and 30.0 L min(-1), respectively. The corresponding storage efficiencies, eta(s), vary from 90% to 76%, respectively. Simulations with L/D ratios ranging from 0.35 to 7.8 indicate that the storage efficiencies can be improved with an increase in the LID ratios and/or with water cooled convection. Thus for an inlet flow rate of 30.0 L min(-1), an eta(s) value of 90% can be obtained with water cooling for an L/D ratio of 7.8 and a filling time of a few minutes. In the absence of water cooling the eta(s) value reduces to 83% at the same L/D ratio. Our study suggests that with an appropriate choice of cylinder dimensions, solutions based on convective cooling during adsorptive storage are possible with some compromise in the storage capacity.

Softening and dynamic recrystallization in magnesium single crystals during c-axis compression

Commercial purity (99.8%) magnesium single crystals were subjected to plane strain compression (PSC) along the c-axis at 200 and 370 degrees C and a constant strain rate of 10(-3) s(-1). Extension was confined to the < 1 1 (2) over bar 0 > direction and the specimens were strained up to a logarithmic true strain of -1. The initial rapid increase in flow stress was followed by significant work softening at different stresses and comparable strains of about -0.05 related to macroscopic twinning events. The microstructure of the specimen after PSC at 200 degrees C was characterized by a high density of {1 0 (1) over bar 1} and {1 0 (1) over bar 3} compression twins, some of which were recrystallized. After PSC at 370 degrees C, completely recrystallized twin bands were the major feature of the observed microstructure. All new grains in these bands retained the same c-axis orientation of their compression twin hosts. The basal plane in these grains was randomly rotated around the c-axis, forming a fiber texture component. The obtained results are discussed with respect to the mechanism of recrystallization, the specific character of the boundaries between new grains and the initial matrix, and the importance of the dynamically recrystallized bands for strain accommodation in these deformed magnesium single crystals. (C) 2011 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Flow from T-spout ladles

The ladle constitutes a crucial element in the pouring system for developing process consistency and quality in the manufacture of castings. Flow of molten metal from T-spout ladles, with the spout projecting at an angle from the shell of the ladle, was investigated by simulation, modeling and analysis. This was followed by experimental validation for water flow in translucent models, and verification in a cast iron foundry. Key parameters in the design of a ladle (for its geometry) have been identified as the ladle size, spout diameter and its angle. Velocity and flow rate from a T-spout ladle are governed by and increase with increase in angle of tilt apart from the aforementioned design parameters. Cross section and profile of the issuing jet displaying a stable twisting pattern, during its free-fall of model fluid water and molten metal iron, are notable aspects in the design of the ladle, particularly the spout exit.