当前位置:首页 >> 能源/化工 >>

Design


Design Guidelines For a Thermoacoustic Refrigerator
Ram C. Dhuley, M.D. Atrey Department of Mechanical Engineering Indian Institute of Technology Bombay, Powai Mumbai-400 076

Development of refrigerators based on Thermoacoustic technology is a novel solution to the present day need of cooling without causing environmental hazards. With added advantages of absence of moving parts and circulating refrigerants, these devices can attain very low temperatures maintaining a compact size. The present theoretical work is based on theory of linear thermoacoustics[1]. Under the short stack and invicid assumptions, an algorithm for design of a standing wave thermoacoustic refrigerator, with main focus on the stack, is described. A stack is designed for a given cooling requirement of the refrigerator and certain chosen operation parameters.

Key Words: design, thermoacoustic, refrigeration

Nomenclature p pressure (N m-2) T f ρ a λ k ω γ δk β K Cp Γ s temperature (K) frequency (Hz) density (kg m-3) sound speed (m s-1) wavelength (m) wave number (m-1) angular frequency (rad s-1) ratio of specific heats thermal penetration depth (m) thermal expansion coefficient (k-1) thermal conductivity (W m-1 k-1) isobaric specific heat (J kg-1 k-1) normalized temperature gradient specific entropy (J kg-1 k-1)

? Q c ? W
COP L Ls A l y0 П xs Subscripts m a 1 p

cooling power (W) acoustic power (W) coefficient of performance length of resonator (m) length of stack (m) cross section area of stack (m2) half plate thickness (m) half plate spacing (m) plate perimeter (m) stack centre position (m)

mean amplitude local amplitude per plate

INTRODUCTION ‘Thermoacoustics’ is the interaction between heat and sound. Of the several applications[1,2] of thermoacoustics in energy conversion, the thermoacoustic refrigerator is a well known device. A thermoacoustic refrigerator uses the energy of sound or pressure waves to bring about refrigeration. The main components of a thermoacoustic refrigerator are the resonator, the stack, the acoustic driver and the heat exchangers as shown below in Figure 1.

OPERATING AND DESIGN PARAMETERS[6] The aim of the present analysis is to develop a design procedure for a thermoacoustic refrigerator and to determine the design parameters based on certain given operating parameters. Various operating and design parameters are given in Table 1. The properties of working gas like density, thermal conductivity, ratio of isobaric to isochoric specific heats and speed of sound in gas play an important part in this analysis. The working gas should have a low boiling point, high sound speed and should be inert to the components of the refrigerator. Table 1 Operating and Design parameters Operating parameters Design Parameters
resonator length, L stack length, Ls stack cross section area, A plate thickness, 2l plate spacing, 2y0 stack centre position, xs

Figure 1 Schematic of thermoacoustic refrigerator The standing pressure wave generated by acoustic driver in the resonator causes to-and-fro motion of the gas from one end of the stack to the other thereby providing a means of transporting heat. The gas near the pressure node gets cooled due to rarefaction and picks up heat from the stack making one end of stack cold. This gas moves towards pressure antinode and gets heated up due to compression. It loses heat to the stack and makes the other end of stack hot. Thus, a temperature gradient is set up along the stack length. With heat exchangers at the two ends of the stack, this temperature gradient can be used to produce refrigeration. A detailed description of thermoacoustic cooling cycle has been given by Arnott et. al[3]. Braun et. al[4] have developed a design optimization program based on the thermoacoustic simulation program known as DELTAE developed by Swift[1]. A design optimization procedure based on normalization of parameters has been given by Wetzel et. al[5] and Tijani et. al[6]

mean pressure, pm mean temperature, Tm frequency, f pressure amplitude, pa

ASSUMPTIONS[1,6] Following are the assumptions made in the present analysis 1. The length of the stack is much smaller as compared to the wavelength of the standing wave Ls<<λ (Short stack approximation). It can be assumed that the local pressure and velocity amplitude of the oscillating gas molecules is more or less same over the entire length of the stack. 2. The flow of the gas in the stack is assumed to be invicid. Friction at the inner walls of the resonator is also neglected.

3. The temperature difference across the stack ends is assumed to be small as compared to the mean temperature ?Tm<<Tm. It can be assumed that the thermo-physical properties of the gas do not vary significantly over the stack length and hence can be assumed constant. 4. The whole analysis is carried out at steady state operation of the refrigerator. The mean temperature of the gas is T m and the temperature gradient across the stack remains constant with time. 5. Conductivity neglected. of plate material is

where k=2π/λ. On account of the porosity ‘B’ of the stack, eqn.3 gets modified to,

u1 ?

pa sin(kx) ?m aB

(4)

The critical temperature gradient is an important parameter governing the refrigeration action of a thermoacoustic device. It is given by,

?Tcrit ?

(? ? 1)kB

?

cot(kx)

(5)

GOVERNING EQUATIONS[1] The equation governing the entire analysis is the general heat transfer equation expressed in terms of entropy transport. Neglecting viscosity of the gas, it can be written as

where β is the thermal expansion coefficient of the gas. For ideal gas, it is equal to the inverse of the absolute temperature. Of utmost concern in the design procedure are the cooling power produced by the stack, the work input needed to produce this cooling power and thus the COP. These are given by the following expressions :-

?T (

?s ? v ??s) ? ? ? ( K ?T ) ?t

(1)

? ? ? ? k ? Tm ? p1u1 (? ? 1) Q p 4

(6)

where s is specific entropy. The equation states that the rate of change of amount of heat at a certain point depends on convective flow of heat by means of entropy and conduction of heat. Various other expressions[1] relevant to present analysis such as those for the cooling power and the acoustic power can be derived from eq(1). Referring to Figure 2, the acoustic driver produces a region of maximum pressure variation at x=0 i.e a pressure antinode and velocity node. Subsequently, following equations can be written down for local pressure and velocity amplitudes in the resonator :-

? k ? LsTm ? 2? ( p1 ) 2 (? ? 1) ? Wp ? 4 ? mC p

(7)

COP ?

Tm ? a tan(kx) BLs (? ? 1)?

(8)

where, δk is the thermal penetration depth of the gas and is an important parameter governing the plate spacing of the stack. It is defined as the distance through which heat diffuses in the gas in time 1/ω and is given by

p1 ? pa cos(kx)

(2)

?k ?

u1 ?

pa sin(kx) ?m a

2K ? m C p?

(9)

(3)

Π is the perimeter of the cross section of the plate. Assuming the thickness of the plate to

be very small, the perimeter is equal to twice the plate width (Π=2w). Γ is the ratio of actual temperature gradient to the critical temperature gradient of the stack. DESIGN PROCEDURE This section describes the design procedure of a refrigerator producing a cooling power of 4 W and a cold end temperature of 210 K. The operating parameters chosen are shown in Table 2[6].

region within δk from the plate and hence the plate spacing is chosen to be 2δ k[8]. Hence, 2y0=0.212 mm.

Table 2 Operating and working gas parameters
Operating parameters pm=10 bar Tm=250 K f=400 Hz pa=0.2 bar Th=283 K Helium properties ρm=1.9244 kg/m3 a=935 m/s k=0.14 W/mK γ=1.67 Plate thickness The plate thickness can be found out from the expression Figure 3 Heat transport vs. distance from plate.

2l ? 2 y0 (

1 ? 1) B

(10)

Determining the design parameters Resonator length The smallest possible length of the resonator which will produce a standing wave is equal to quarter of the wavelength. A standing wave can be generated by keeping one end of resonator closed (driver end) and simulating the other end as ‘open’. This can be done by attaching a sufficiently large buffer[6,7] of gas at Tm at this end. For this most fundamental case we get L=58.8 cm. Plate spacing The thermal penetration depth of the gas at the given frequency is found out to be 0.106 mm. Figure 3 shows variation of heat flux with distance from the plate. It can be seen that almost all the transport of heat takes place in

where B is the porosity of the stack. The porosity is defined as the ratio of open area in the cross section to the total cross section area of the stack. The porosity is chosen such as not to disturb the acoustic standing wave significantly and is normally taken in the range 0.7-0.8[5,6]. Choosing B=0.75, we get 2l=0.07 mm. Stack length, centre position and area of cross section A graphical approach has been implemented to determine the stack length, centre position and the cross section area. It is assumed that the stack is made of ‘n’ parallel plates over one another and the cross section of the stack is a square of width ‘w’. In this case, the perimeter of the p late cross section becomes 4nl+4(n-1)y0. Knowing the fact that total cooling power produced by

the stack is ‘n’ times that produced by a single plate, it can be written

? ? ?n[n(l ? y ) ? y ] ?? p u [ ?Tm ? 1] (11) Q k 1 1 c 0 0 Ls ?Tcrit
Putting the known quantities in eq(11), an equation in three variables x, L s and n is obtained,

obtained. The cross section area of this stack comes out to be 7.06 cm2. This stack will produce a cooling power of 4 W at COP close to 2.

0.282n 2 ? 0.212n ?
509.54 ?0 tan(kx) sin(2kx)(0.221 ? 1) Ls
Similarly, from eq(8) we get, (12)

tan(kx) COP ? 0.75 Ls

(13)

Figure 5 Variation of number of plates with stack length. RESULTS AND DISCUSSIONS The design parameters obtained from the analysis are shown in Table 3 Table 3 Results Parameter L Ls A 2l 2y0 xs Value 588 mm 100 mm 706 mm2 0.07 mm 0.212 mm 100 mm

COP and n are plotted with Ls for different positions of stack centre in the resonator. Theoretically, every point on such a curve will give us a stack producing a cooling power of 4 W and a cold end temperature 210 K. Figure 4 and Figure 5 show variation of COP and n respectively with L s at different stack centre locations.

Figure 4 Variation of COP with stack length For instance, if Ls=0.1 m is chosen at xs=0.1 m, a stack with 95 parallel plates is

The stack length (Ls=10 cm) is very small as compared to the wavelength of the acoustic field (λ=235.2 cm). Hence, the assumption of short stack is perfectly valid. At the stack centre location (xs=10 cm) the critical temperature gradient is 1235.13 K/m. The actual temperature gradient across the stack ends is 750 K/m. Clearly, Γ=0.607<1 which is primary necessity[1] for the device to function as a refrigerator, is met. Another perfectly valid choice of stack producing the same cooling power at same

cold end temperature is L s=8.5 cm and xs=8cm. In this case, the stack will consists of 102 plates and a cross section of 8.2 cm2. The corresponding COP is 1.9. Tijani et al.[6] obtained a stack with cross section 11.8 cm2 and COP of 1.5 with same stack length and set of operation parameters. The difference in the cross section can be accounted for the authors in [6] had taken viscosity into account and hence needed more number of stack plates to produce 4 W power. Referring to Figure 4, high values of COP (~3) are expected from the analysis. But, it should be noted that this is COP of the stack and not of the refrigerator. Even when the viscosity is neglected, more power has to be input to account for losses at heat exchangers, driver, etc. (though not a part of this analysis). CONCLUSIONS A simple design procedure for a standing wave thermoacoustic refrigerator has been described. With a choice of operating parameters and helium as working gas, graphical approach has been used to determine the geometrical parameters of the stack. A refrigerator producing a cooling power of 4W at cold end temperature of 210K at steady state is designed theoretically using this design procedure. References
1. Swift G.W., Thermoacoustic engines, J Acoust Soc. Am, 84, (1998), pp 1146-1180 2. Swift G.W., Thermoacoustics- A unifying perspective of some engines and refrigerators, Acoustical Society of America Publications (2002). 3. Arnott W, Raspet R, Bass H.E, Thermoacoustic Engines, Ultrasonics Symposium, (1991), pp 995-1003. 4. Paek I.,Braun J.E., Mongeau L., Evaluation of standing-wave thermoacoustic cycles for

cooling applications, International Journal of Refrigeration, 32, (2007), pp 1059-1071. 5. Wetzel M., Herman C., Design Optimization of thermoacoustic refrigerators, International Journal of Refrigeration, 20, (1997), pp 3-21. 6. Tijani M.E.H., Zeegers J.C.H., deWaele A.T.A.M.., Design of thermoacoustic refrigerators, Cryognenics, 42, (2002), pp 4957 7. Tijani M.E.H., Zeegers J.C.H., deWaele A.T.A.M.., Construction and working of a thermoacoustic refrigerator, Cryognenics, 42, (2002), pp 59-66 8. Tijani M.E.H., Zeegers J.C.H., deWaele A.T.A.M., Optimal stack spacing for thermoacoustic Refrigeration, J Acoust Soc. Am, 112(1), (2002), pp 128-133


相关文章:
PROTEL Design设计菜单下RULES规则
Design 设计菜单下 RULES 规则 Routing 布线设计规则 1 Clearance Constraint-安全间距 Routing Corners-布线转角 Routing Layers-布线板层 Routing Priority-布线次序 ...
语言学中的五个design features
语言学中的五个design features_文学研究_人文社科_专业资料。今天我们来谈一谈语言学中的五个 design features. Chinese meaning 设计特点。 在讲五个设 计特点...
设计公司介绍
在 美国 BILKEY LLINAS DESIGN 设计公司的具体介绍(BLD 简称)大卫.罗克威尔(David Rockwell,以下简称 DR)是一位著名的美国建筑师和室内设计师,也是罗克威尔设计 ...
3DESIGN 和 JewelCAD 软件比较_ver2
3DESIGN 和 JewelCAD 软件比较_ver2_计算机硬件及网络_IT/计算机_专业资料。3DESIGN 和 JewelCAD 软件比较功能特征 曲面模型设计 实体模型设计 3DESIGN JewelCAD ...
RD方法(regression discontinuity design)
首先要说的是这篇文章使用的 RD 方法 (regression discontinuity design) 即断点回归方法。 , 断点回归方法是最近的政策评估中非常重要的一个方法, 他可以在没有...
lesson_5_design_and_culture
28. in tow directions 从两个方向 unacknowledged 不被承认的,未答复的 vital 极重要的,不可缺少的 the mass media 大众传媒 latter 后者的 guise design 外观...
GUI_Design_Studio使用手册
(链接 地址为 "C:\Documents and Settings\All Users\Documents\GUI Design Studio\Samples\Standard Components\Standard Components.gdp") 5) 6) 7) 8) 9) ...
Cadence Design entry HDL 教程_图文
Cadence Design entry HDL 教程_计算机软件及应用_IT/计算机_专业资料。Design entry HDL 教程 Concept 是 Cadence 公司自身开发的原理图输入工具,在业界拥有广泛的...
DesignBuilder操作手册
DesignBuilder操作手册_电脑基础知识_IT/计算机_专业资料。DesignBuilder 软件操作手册注意事项: 1.安装时需要注意的问题:安装目录及整个路径须为全英文,不能出现中文...
DHF-design history file
Design history file From Wikipedia, the free encyclopedia Jump to: navigation, search Design History File is a compilation of documentation that describes ...
更多相关标签: