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

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


相关文章:
123D Design实用详解(二)_图文
123D Design 实用详解(二)移动和合并 今天开始学习简单模型制作和各个工作栏的部分应用, 比如上一 次说到的基本几何体,它是一个非常方便的实体模型,方便我们在创...
Altium Design PCB拼板完整教程
Altium Design PCB拼板完整教程_计算机软件及应用_IT/计算机_专业资料。关于 PCB 拼板完整教程 一、为什么拼板 电路板设计完以后需要上 SMT 贴片流水线贴上元器件...
In Design 页面基本知识
In Design 页面基本知识_计算机软件及应用_IT/计算机_专业资料。10.1 页面基本操作在报纸、书籍和杂志的排版过程中,经常要用到页面的操作。在中文化 InDesign 中,...
DesignBuilder操作手册
DesignBuilder操作手册_电脑基础知识_IT/计算机_专业资料。DesignBuilder 软件操作手册注意事项: 1.安装时需要注意的问题:安装目录及整个路径须为全英文,不能出现中文...
Mechanism Design Theory
Mechanism Design Theory_经济/市场_经管营销_专业资料。背景知识:机制设计理论(Mechanism Design Theory) 一、什么是机制设计理论? 简单的讲,经济机制设计理论是研究...
In Design教学课件
In Design 备课第 1 课:工作界面的认识一:文档操作 1:新建页数不可以超过 9999 页; 2: “对页” :可以同时显示两个正在编辑的页面; 3:文件——新建:Ctrl+...
Design Explorer快速学习课件
Design Explorer快速学习课件_工学_高等教育_教育专区 暂无评价|0人阅读|0次下载|举报文档 Design Explorer快速学习课件_工学_高等教育_教育专区。ONEDAY 第一讲 1...
design compiler使用说明
60 用Design Analyzer设置特征参数: 在读入设计后,要说明设计的环境变量,这些环境变量包括了当信号到达端口(port)时, 他们的的驱动能力,或者输出端口的负载能力。 ...
GUI Design Studio使用手册
鼠标事件 鼠标事件 GUI Design Studio 提供 3 种鼠标事件: 鼠标经过事件——鼠标经过活动区域时触发该事件 常用于"Toggle Window"类型的连接,用于当鼠标在活动区域...
Design_Explorer快速学习课件
Design_Explorer快速学习课件_其它课程_高中教育_教育专区 暂无评价|0人阅读|0次下载|举报文档 Design_Explorer快速学习课件_其它课程_高中教育_教育专区。软件工具 ...
更多相关标签:
designer | designed | design by | design expert | android design | material design | design是什么牌子 | design for |