Numerical analysis on thermoacoustic prime movers for development of pulse tube cryocoolers

Thermo Acoustic Prime Movers (TAPMs) are being considered as the ideal choice for driving the Pulse Tube Cryocoolers replacing the conventional compressors. The advantages are the absence of moving components and they can be driven by low grade energy as such as fuel, gas, solar energy, waste heat etc. While the development of such TAPMs is in progress in our laboratory, their design and fabrication should be guided by numerical modeling and this may be done by several methods such as solving the energy equation 1], enthalpy flow model 2], CFD 3], etc. We have used CFD technique, since it provides a better insight into the velocity and temperature profiles. The analysis is carried out by varying parameters such as (a) temperature difference across the stack, (b) stack and resonator lengths and (c) different working fluids such as air, nitrogen, argon etc. The theoretical results are compared with the experimental data wherever possible and they are in reasonably good agreement with each other. The analysis indicate that (i) larger temperature difference across the stack leads to increased acoustic amplitude, (ii) longer resonator leads to decrease in frequency with lesser amplitude and (iii) there exists an optimal stack length for the best performance of TAPM. These results are presented here.

Development of moving magnet type linear motor for dual piston compressor for pulse tube cryocooler

This paper describes the design, fabrication and testing of a moving magnet type linear motor of dual piston configuration for a pulse tube cryocooler for ground applications. Eight radially magnetized segmented magnets were used to form one set of a magnet ring. Four magnet rings of such type were constructed, in which one pair of rings has north-pole on its outer diameter and south-pole on inner diameter, while the other pair is it's complementary. The magnets were mounted with opposite poles together on the magnet holder with an axial moving shaft having a piston mounted on both ends of the shaft. The shaft movement was restricted to the axial direction by using C-clamp type flexures, mounted on both sides of the shaft. The force requirement for driving the compressor was calculated based on which the electrical circuit of motor is designed by proper selection of wire gauge and Ampere-turns. The flexure spring force estimation was done through simulation using ANSYS 11.0 and was verified experimentally; while the magnet spring force was determined experimentally. The motor with mounted piston was tested using a variable voltage and variable frequency power supply capable of driving 140 watts of load.

Performances of single and two-stage pulse tube cryocoolers under different vacuum levels with and without thermal radiation shields

Single and two-stage Pulse Tube Cryocoolers (PTC) have been designed, fabricated and experimentally studied. The single stage PTC reaches a no-load temperature of similar to 29 K at its cold end, the two-stage PTC reaches similar to 2.9 K in its second stage cold end and similar to 60 K in its first stage cold end. The two-stage Pulse Tube Cryocooler provides a cooling power of similar to 250 mW at 4.2 K. The single stage system uses stainless steel meshes along with Pb granules as its regenerator materials, while the two-stage PTC uses combinations of Pb along with Er3Ni/HoCu2 as the second stage regenerator materials. Normally, the above systems are insulated by thermal radiation shields and mounted inside a vacuum chamber which is maintained at high vacuum. To evaluate the performance of these systems in the possible conditions of loss of vacuum with and without radiation shields, experimental studies have been performed. The heat-in-leak under such severe conditions has been estimated from the heat load characteristics of the respective stages. The experimental results are analyzed to obtain surface emissivities and effective thermal conductivities as a function of interspace pressure.

High-temperature kinetics of the reaction between CN and hydrocarbons using a novel high-enthalpy flow tube

More than 70 molecules of varied nature have been identified in the envelopes of carbon-rich stars through their spectral fingerprints in the microwave or far infrared regions. Many of them are carbon chain molecules and radicals, and a significant number are unique to the circumstellar medium. The determination of relevant laboratory kinetics data is critical to keep up with the development of the high spectral and spatial resolution observations and of the refinement of chemical models. Neutralneutral reactions of the CN radical with unsaturated hydrocarbons could be a dominant route in the formation of cyanopolyynes, even at low temperatures and deserve a detailed laboratory investigation. The approach we have developed aims to bridge the temperature gap between resistively heated flow tubes and shock tubes. The present kinetic measurements are obtained using a new reactor combining a high-enthalpy source with a flow tube and a pulsed laser photolysislaser-induced fluorescence system to probe the undergoing chemical reactions. The high-enthalpy flow tube has been used to measure the rate constant of the reaction of the CN radical with propane (C3H8), propene (C3H6), allene (C3H4), 1,3-butadiene (1,3-C4H6), and 1-butyne (C4H6) over a temperature range extending from 300 to 1200 K. All studied reactions of CN with unsaturated hydrocarbons are rapid, with rate coefficients greater than 10-10 cm3 center dot molecule-1 center dot s-1 and exhibit slight negative temperature dependence above room temperature. (c) 2012 Wiley Periodicals, Inc. Int J Chem Kinet 44: 753766, 2012

Theoretical-experimental study of shock wave-assisted metal forming process using a diaphragmless shock tube

The use of high-velocity sheet-forming techniques where the strain rates are in excess of 10(2)/s can help us solve many problems that are difficult to overcome with traditional metal-forming techniques. In this investigation, thin metallic plates/foils were subjected to shock wave loading in the newly developed diaphragmless shock tube. The conventional shock tube used in the aerodynamic applications uses a metal diaphragm for generating shock waves. This method of operation has its own disadvantages including the problems associated with repeatable and reliable generation of shock waves. Moreover, in industrial scenario, changing metal diaphragms after every shot is not desirable. Hence, a diaphragmless shock tube is calibrated and used in this study. Shock Mach numbers up to 3 can be generated with a high degree of repeatability (+/- 4 per cent) for the pressure jumps across the primary shock wave. The shock Mach number scatter is within +/- 1.5 per cent. Copper, brass, and aluminium plates of diameter 60 mm and thickness varying from 0.1 to 1 mm are used. The plate peak over-pressures ranging from 1 to 10 bar are used. The midpoint deflection, circumferential, radial, and thickness strains are measured and using these, the Von Mises strain is also calculated. The experimental results are compared with the numerical values obtained using finite element analysis. The experimental results match well with the numerical values. The plastic hinge effect was also observed in the finite element simulations. Analysis of the failed specimens shows that aluminium plates had mode I failure, whereas copper plates had mode II failure.

Minimum power consumption for drag reduction on a circular cylinder by tangential surface motion

We investigate the effect of a prescribed tangential velocity on the drag force on a circular cylinder in a spanwise uniform cross flow. Using a combination of theoretical and numerical techniques we make an attempt at determining the optimal tangential velocity profiles which will reduce the drag force acting on the cylindrical body while minimizing the net power consumption characterized through a non-dimensional power loss coefficient (C-PL). A striking conclusion of our analysis is that the tangential velocity associated with the potential flow, which completely suppresses the drag force, is not optimal for both small and large, but finite Reynolds number. When inertial effects are negligible (R e << 1), theoretical analysis based on two-dimensional Oseen equations gives us the optimal tangential velocity profile which leads to energetically efficient drag reduction. Furthermore, in the limit of zero Reynolds number (Re -> 0), minimum power loss is achieved for a tangential velocity profile corresponding to a shear-free perfect slip boundary. At finite Re, results from numerical simulations indicate that perfect slip is not optimum and a further reduction in drag can be achieved for reduced power consumption. A gradual increase in the strength of a tangential velocity which involves only the first reflectionally symmetric mode leads to a monotonic reduction in drag and eventual thrust production. Simulations reveal the existence of an optimal strength for which the power consumption attains a minima. At a Reynolds number of 100, minimum value of the power loss coefficient (C-PL = 0.37) is obtained when the maximum in tangential surface velocity is about one and a half times the free stream uniform velocity corresponding to a percentage drag reduction of approximately 77 %; C-PL = 0.42 and 0.50 for perfect slip and potential flow cases, respectively. Our results suggest that potential flow tangential velocity enables energetically efficient propulsion at all Reynolds numbers but optimal drag reduction only for Re -> infinity. The two-dimensional strategy of reducing drag while minimizing net power consumption is shown to be effective in three dimensions via numerical simulation of flow past an infinite circular cylinder at a Reynolds number of 300. Finally a strategy of reducing drag, suitable for practical implementation and amenable to experimental testing, through piecewise constant tangential velocities distributed along the cylinder periphery is proposed and analysed.

A constant factor approximation algorithm for boxicity of circular arc graphs

Boxicity of a graph G(V, E) is the minimum integer k such that G can be represented as the intersection graph of k-dimensional axis parallel boxes in Rk. Equivalently, it is the minimum number of interval graphs on the vertex set V such that the intersection of their edge sets is E. It is known that boxicity cannot be approximated even for graph classes like bipartite, co-bipartite and split graphs below O(n0.5-ε)-factor, for any ε > 0 in polynomial time unless NP = ZPP. Till date, there is no well known graph class of unbounded boxicity for which even an nε-factor approximation algorithm for computing boxicity is known, for any ε < 1. In this paper, we study the boxicity problem on Circular Arc graphs - intersection graphs of arcs of a circle. We give a (2+ 1/k)-factor polynomial time approximation algorithm for computing the boxicity of any circular arc graph along with a corresponding box representation, where k ≥ 1 is its boxicity. For Normal Circular Arc(NCA) graphs, with an NCA model given, this can be improved to an additive 2-factor approximation algorithm. The time complexity of the algorithms to approximately compute the boxicity is O(mn+n2) in both these cases and in O(mn+kn2) which is at most O(n3) time we also get their corresponding box representations, where n is the number of vertices of the graph and m is its number of edges. The additive 2-factor algorithm directly works for any Proper Circular Arc graph, since computing an NCA model for it can be done in polynomial time.

Manually operated piston-driven shock tube

A simple hand-operated shock tube capable of producing Mach 2 shock waves is described. Performance of this miniature shock tube using compressed high pressure air created by a manually operated piston in the driver section of the shock tube as driver gas with air at 1 atm pressure as the test gas in the driven tube is presented. The performance of the shock tube is found to match well with the theoretically estimated values using normal shock relations. Applications of this shock tube named Reddy tube, include study of blast-induced traumatic brain injuries and high temperature chemical kinetics.

Tuning of the extended concentric tube resonators

It has been shown recently that the acoustic performance of the extended tube expansion chambers can be improved substantially by making the extended inlet and outlet equal to half and quarter chamber lengths, duly incorporating the end corrections due to the evanescent higher order modes that would be generated at the discontinuities. Such chambers however suffer from the disadvantages of high back pressure and generation of aerodynamic noise at the area discontinuities. These two disadvantages can be overcome by means of a perforated bridge between the extended inlet and extended outlet. This paper deals with design or tuning of these extended concentric tube resonators.

Asymptotic wavenumber expansions of a strongly orthotropic fluid-filled circular cylindrical shell

In this paper, a suitable nondimensional `orthotropy parameter' is defined and asymptotic expansions are found for the wavenumbers in in vacuo and fluid-filled orthotropic circular cylindrical shells modeled by the Donnell-Mushtari theory. Here, the elastic moduli in the two directions are greatly different; the particular case of E-x >> E-theta is studied in detail, i.e., the elastic modulus in the longitudinal direction is much larger than the elastic modulus in the circumferential direction. These results are compared with the corresponding results for a `slightly orthotropic' shell (E-x approximate to E-theta) and an isotropic shell. The novelty of this presentation lies in obtaining closed-form expansions for the in vacuo and coupled wavenumbers in an orthotropic shell using perturbation methods aiding in a better physical understanding of the problem.

Computational fluid dynamics simulation of open ended thermoacoustic prime mover with straight and tapered resonator

A thermoacoustic refrigerator driven by a thermoacoustic primemover is an effective way to produce durable and long lasting refrigeration due to high reliability, no exotic materials, and no moving parts. Resonator geometry is also one of the important factors that influence the performance of a thermoacoustic prime mover, namely, frequency. Computational fluid dynamics simulation of performance comparison of thermoacoustic prime mover with a straight and tapered resonator is chosen for the present study under an identical stack condition with the air as a working fluid. The frequency and pressure amplitude of oscillations obtained from simulation results were found to be more in the tapered resonator than the straight resonator. Apart from computational fluid dynamics simulation, the simulation studies have also been conducted using design environment for low-amplitude thermoacoustic energy conversion, which predicts the performance of thermoacoustic primemover comparatively well. Simulation results from computational fluid dynamics and design environment for low-amplitude thermoacoustic energy conversion were compared and found to be matching well, representing the good validity of computational fluid dynamics modeling.

Recent advances in muffler acoustics

Exhaust noise in engines has always been a major source of automotive noise. Challenges for muffler design have been constraints on size, back pressure, and, of course, the cost. Designing for sufficient insertion loss at the engine firing frequency and the first few harmonics has been the biggest challenge. Most advances in the design of efficient mufflers have resulted from linear plane wave theory, making use of the transfer matrix method. This review paper deals with evaluating approximate source characteristics required for prediction of the unmuffled intake and exhaust noise, making use of the electroacoustical analogies. In the last few years, significant advances have been made in the analysis of variable area perforated ducts, transverse plane wave analysis of short elliptical as well as circular chambers, double-tuned expansion chambers and concentric tube resonators, catalytic converters, diesel particulate filters, air cleaners, etc. The development of long strand fibrous materials that can be used in hot exhaust systems without binders has led to the use of combination mufflers in exhaust systems. Breakthroughs have been achieved in the prediction and control of breakout noise from the elliptical and circular muffler shell as well as the end plates of typical mufflers. Diesel particulate filters and inlet air cleaners have also been modeled acoustically. Some of these recent advances are the subject of this review paper.

Interfacial growth as a model of tube-width heterogeneities in concentrated solutions of stiff polymers

Recent experimental measurements of the distribution P(w) of transverse chain fluctuations w in concentrated solutions of F-actin filaments B. Wang, J Guan, S. M. Anthony, S. C. Bae, K. S. Schweizer, and S. Granick, Phys. Rev. Lett. 104, 118301 (2010); J. Glaser, D. Chakraborty, K. Kroy, I. Lauter, M. Degawa, N. Kirchgessner, B. Hoffmann, R. Merkel, and M. Giesen, Phys. Rev. Lett. 105, 037801 (2010)] are shown to be well-fit to an expression derived from a model of the conformations of a single harmonically confined weakly bendable rod. The calculation of P(w) is carried out essentially exactly within a path integral approach that was originally applied to the study of one-dimensional randomly growing interfaces. Our results are generally as successful in reproducing experimental trends as earlier approximate results obtained from more elaborate many-chain treatments of the confining tube potential.

Asymptotic expansions for the structural wavenumbers of isotropic and orthotropic fluid-filled circular cylindrical shells in the intermediate frequency range

We consider wavenumbers in in vacuo and fluid-filled isotropic and orthotropic shells. Using the Donnell-Mushtari (DM) theory we find compact and elegant asymptotic expansions for the wavenumbers in the intermediate frequency range, i.e., around the ring frequency. This frequency range corresponds to the frequencies where there is a rapid change in the values of bending wavenumbers and is found to exist in isotropic and orthotropic shells (in vacua and fluid-filled) for low circumferential orders n only. The same is first identified using the n=0 mode of an orthotropic shell. Following this, using the expression for the intermediate frequency, asymptotic expansions are found for other cases. Here, in order to get compact expansions we consider slight orthotropy (epsilon << 1) and light fluid loading (mu << 1). Thus, the orthotropy parameter epsilon and the fluid loading parameter mu are used as asymptotic parameters along with the non-dimensional thickness parameter beta. The methodology can be extended to any order of epsilon, only the expansions become unwieldy. The expansions are matched with the numerical solutions of the corresponding dispersion relation. The match is found to be good.

Stability of a long unsupported circular tunnel in soils with seismic forces

By using the lower-bound finite element limit analysis, the stability of a long unsupported circular tunnel has been examined with an inclusion of seismic body forces. The numerical results have been presented in terms of a non-dimensional stability number (gamma H/c) which is plotted as a function of horizontal seismic earth pressure coefficient (k (h)) for different combinations of H/D and I center dot; where (1) H is the depth of the crest of the tunnel from ground surface, (2) D is the diameter of the tunnel, (3) k (h) is the earthquake acceleration coefficient and (4) gamma, c and I center dot define unit weight, cohesion and internal friction angle of soil mass, respectively. The stability numbers have been found to decrease continuously with an increase in k (h). With an inclusion of k (h), the plastic zone around the periphery of the tunnel becomes asymmetric. As compared to the results reported in the literature, the present analysis provides a little lower estimate of the stability numbers. The numerical results obtained would be useful for examining the stability of unsupported tunnel under seismic forces.