Hydrogen, helium and thermo-acoustic refrigerators

In this work the design and analysis of 1 kW thermo-acoustic refrigerators with hydrogen and helium for the temperature difference of 38 K is discussed. Helium is the best for thermoacoustic refrigerators compared to the other competent gases. But hydrogen is chosen since it is less expensive and better thermophysical properties compared to helium. The best parallel plates geometry with 15% blockage is chosen for the stack and heat exchangers. The effect of resonance frequency of hydrogen and helium varying from 400–600 Hz on the theoretical performance is discussed. The coefficient of performance and the power density of 1.65 and 40.3 kW/m3 for hydrogen, and 1.58 and 19.2 kW/m3 for helium is reported for the optimized designs, respectively. The theoretical results are compared with the DeltaEC software results, shows the cooling power and coefficient of performance of 590 W and 1.11 for hydrogen, and 687 W and 1.25 for helium, respectively.


Introduction and review of literature
The stack-based standing wave thermoacoustic refrigerators are the low cost, eco-friendly and promising technology and which can be built indigenously with only moving component as magnetic coil for room temperature refrigeration.The thermoacoustic coolers are expected to replace the existing vapour compressor refrigeration systems for household and commercial applications.In the recent years, thermoacoustic cooling receives the attentions of the researchers worldwide to optimize the coolers both theoretically and practically, ranging from 3-500 W at a reasonable and promising coefficient of performance of 0.5-1.8[1][2][3][4].Thermoacoustic coolers may be driven by thermoacoustic engine or loudspeaker.The loudspeaker that drives the thermoacoustic cooler is the simple and readily available in the commercial market, and which can be modified easily to improve efficiency compared to the thermoacoustic engines [5][6][7].The pictorial representation of the thermo-acoustic refrigerator powered by commercial loudspeaker is depicted in Fig. 1.The commercial loudspeaker firmly fixed with the duct left to the ambient heat exchanger, followed by the parallel plate stack and cold heat-exchanger.The taper section fastened right to the cold heat-exchanger.The compact quarter wave length divergent-hemispherical resonator fastened to the taper section (Fig. 1) is the efficient design over the small diameter tube design as discussed elsewhere [3,4,8].
The loudspeaker supplies acoustic energy because of its vibrating diaphragm to the filled in gas of a resonator.Upon electrical energy to the loudspeaker, the gas starts oscillating through the porous stack-heat-exchangers, the taper and the divergent section.The methodology of choosing the geometrical specifications and the positions of the stack, heat-exchangers, resonator and the loudspeaker of the cooler is found elsewhere [9][10][11].The gas experiences maximum compression in the ambient heatexchanger region and maximum rarefaction (expansion) in the cold heat exchanger region during the front and back motion of the coil, respectively.The oscillating gas is considered as the adiabatic parcel since the whole outer surface of the device is assumed to be insulated perfectly with no heat leak losses.
The four thermodynamic processes (Fig. 2) simar to the Bell Coleman's air refrigeration cycle (reverse Brayton or Joule cycle) are: the process (1-2) and (3-4) represents adiabatic gas compression and expansion, respectively and (4-1) and (2-3) represents isobaric heat absorption and rejection, respectively.The P-v, T-s and P-h thermodynamic diagrams with the variation in the temperature, enthalpy and size of the oscillating gas parcel is given in Fig. 2. The gas oscillations gradually decrease the gas temperature in the cold heat-exchanger region, absorbing heat of the cold heat-exchanger plates connected to cooling space.Contrary the gas temperature increases in the ambient heat-exchanger region, reject heat to the atmosphere via the ambient heat-exchanger plates at room temperature using the water-cooling system [1,3].
Likewise, the compression and expansion of gas caused by the loudspeaker creates difference in temperature across the ambient and cold heat-exchangers through the stack for heat pumping from the cooling space to the surrounding environment (Fig. 1).In the literature it is found that the working gases for the thermoacoustic refrigerators is chosen based on the better thermophysical properties (high sound velocity and thermal conductivity) and low cost [1,10].And the other criterion for the selection of the best gas is the Prandtl number (σ), is the ratio of viscous to thermal diffusion.From the literature [12] it is found that the cooler performance increases with decrease in the Prandtl number of the oscillating gas.Hence the pure or the gas mixture with low Prandtl number is the best for the thermoacoustic coolers.Prandtl number is proportional to dynamic viscosity (μ) and isobaric specific heat (C p ) and varies inversely with thermal conductivity (k) of the oscillating gas.In this research work hydrogen and helium is chosen since both the gases meets the required conditions as the working gases (Table 1).The hydrogen has the higher isobaric specific heat, leading to the marginal increase in the Prandtl Fig. 1 The pictorial representation showing the thermo-acoustic refrigerator powered by commercial loudspeaker Fig. 2 The P-v, T-s, and P-h diagrams of a thermo-acoustic cooler with circles around the state points representing the size of a gas parcel undergoing thermodynamic process number compared to the helium.In this research work the thermoacoustic cooler with 1 kW capacity is theoretically optimized with hydrogen and helium using the parallel plates stack and heat exchangers system.Therefore, the primary objective is to study the effect of resonance frequency of hydrogen and helium operating from 400-600 Hz on the theoretical performance of the cooler is discussed.
The resonance frequency for both the gases is decided based on their better cooling capacity (Q), performance (COP) and power-density of the cooler for the compact design.Hydrogen is cheaper compared to helium and hence an attempt is made to use hydrogen as the competent working gas over helium for the possibility to decrease the cost of the coolers towards commercialization.But the hydrogen is highly inflammable and with proper safety measures the use of hydrogen in thermoacoustic coolers is justifiable.The second objective is to study the effect of increasing the length of the cold heatexchanger (l 4 ) from two times the gas displacement (2x 1 ) to four times (4x 1 ) on the theoretical performance is discussed in Sects. 2 and 3.In this research work the length of the cold heat-exchanger equal to four times the gas displacement (l 4 = 4x 1 ) is considered to avoid the space constraints associated in mounting the thermocouples with the small length at two times the gas displacement (l 4 = 2x 1 ) as discussed elsewhere [3].The theoretically optimized design results with hydrogen and helium are validated with the DeltaEC software simulation results [13], and the conclusions are drawn by considering the design assumptions and limitations.

Design and analysis of parallel plates stack and heat-exchangers
The design and analysis of the parallel plates stack and heat-exchangers system at 85% porosity 'ε' is discussed.The cross-sectional geometry for the parallel plates stack, cold and ambient heat-exchangers are the same but the length varies.The porosity represents the percentage of the cross-sectional area of the stack and heat-exchangers open for the gas to oscillate.In other words, it is also called as the blockage ratio 'B', represents the percentage of the solid area offers blockage for the gas to oscillate along its length.The porosity and blockage ratio are given by where 'l' and 'y' are referred to the parallel plates halfthickness and spacing, respectively of the stack-heat exchangers system (Fig. 3).The stack is usually made of mylar material, whereas the cold and ambient (1) Fig. 3 A pictorial view of the parallel plates stack and heat exchanger heat-exchangers is of copper, both made precisely using the popular 3D printing technology.The thermophysical properties of the mylar stack and copper parallel plates at the mean gas temperature T mg 283 K are: k s = 0.158 Wm −1 K −1 , ρ s = 1350.5kgm −3 and C ps = 1047.1 Jkg −1 K −1 for the stack; and k c = 398.96Wm −1 K −1 , ρ c = 9000 kgm −3 and C pc = 420 Jkg −1 K −1 for the cold and ambient heatexchangers.The 1 kW cooler design is optimized with hydrogen and helium for the temperature difference (θ) of 38 K, assuming the ambient heat-exchanger temperature (T ax ) at 301 K using the circulating cooling water system as found in the published literature [2,3].The cooler is operating with 10 bar pressure at 3% drive ratio 'D' is the ratio of dynamic pressure amplitude of oscillating gas 'P a ' to gas pressure 'P' .
In the literature, it is found that the drive ratio should be less than or equal to 3% to avoid the acoustic non-linearities of the system [14].The viscous penetration depth 'δ v ' and the thermal penetration depths 'δ k ' of the oscillating hydrogen and helium are the very thin layers around the parallel plates stack, are given by [9] The oscillating gas should have comparatively a lower value of 'δ v ' which resist the gas movement compared to 'δ k ' through which the thermo-acoustic heat pumping occurs.The design parameters and the variation (3) in the values of the penetration depths ('δ v ' and 'δ k ') of the hydrogen and helium versus the resonance frequency 'f' are given in Table 2.The penetration depths are inversely proportional to the resonance frequency of the oscillating gas.
The amplitude of velocity and pressure of the oscillating gas measured at a point 'x' away from the speaker surface is given by [5] The gas parcel excursion (displacement amplitude, x 1 ) is expressed as [5] where 'u 1 'and 'ω' are the amplitude of gas velocity and angular velocity at a particular location, respectively.For a typical engine or refrigerator, the penetration depths are smaller than the oscillating gas parcel excursion x 1 is given by [15] From the literature it is found that to minimize the altering of the acoustic standing wave in the stack and heat exchangers system, it is recommended to use the parallel plates oscillating gas spacing (2y) of 2δ k to 4δ k [16][17][18].Hence it is chosen to use a gas spacing of 4δ k to avoid the acoustic disturbances.The stack and heat exchangers plates thickness (2 l) at 85% porosity is calculated using Eq. ( 1).The stack and heat-exchangers plates spacing (2y) and thickness (2 l) varies inversely with the resonance frequency for both hydrogen and helium as given in Table 3.The derived parameters, namely the angular velocity (ω), wave number (c), and wavelength (λ) of the hydrogen and helium versus the resonance frequency are presented in Table 4.For the stack-based thermoacoustic cooler or engine, the ( 5)

Design parameters
Viscous and thermal penetration depths Design requirements:  parameter which determines the cooler performance is the stack critical temperature (θ cr ).
It is the state of the thermo-acoustic stack at which the stack does not pump heat with no temperature difference (θ) exist, is given by [15] The design parameters are normalized to the 11 dimensionless parameters to simplify the design process and the effect of resonance frequency on the normalized parameters for hydrogen and helium are given in Table 5.If θ < θ cr then the normalized critical temperature 'Γ' (ratio of temperature difference to critical temperature) will be less than one and is true for both hydrogen and helium (Table 5).At this condition the stack plates transfer heat from the cold heat-exchanger to the ambient heat-exchanger to produce refrigeration effect in the cooling space upon the acoustic power input by the loudspeaker (Fig. 1) [19,20].The hydrogen has the lower value of the normalized critical temperature (0.4) compared to helium (0.5).Hence the hydrogen cooler can produce better refrigeration effect compared to helium at (9) the required temperature difference (θ) of 38 K (Table 2).The normalized stack-length 'l 3n ' and centre-position 'X n ' values are optimized between the stack diameter (d 1 ) and the stack coefficient of performance (COP s ).They are chosen based on the balance between 'd 1 ' and 'COP s ' .The thermal conductivity 'k' of the parallel plates mylar stack affects its performance because of the reverse conduction of heat from the ambient to cold heat-exchanger and decreases the 'COP s ' as discussed elsewhere [1,6,21].The normalized cooling-power 'Q n ' and acoustic-power 'W n ' equations are derived by considering the heat-capacity ratio and the thermal conductivity of both the gas and the stack using the Rott's thermoacoustic theory is given by [9,22,23] where ε s , Λ and Λ 1 in Eqs.(10) and (11) are defined as the dimensionless stack heat capacity correction factor, and (10)  the gas and the stack heat conduction correction factors, respectively, are given by Using the data given in Tables 1, 2, 3 and 4 in Eqs. ( 10) -( 11), the 'Q n ' and 'W n ' for the optimized parallel plates stack are determined and presented in Table 5.The coefficient of performance of the stack 'COP s ' is given by The stack performance versus the resonance frequency for hydrogen and helium is calculated using Eq. ( 15) (Table 6).The performance of the stack varies inversely with the resonance frequency because of the increase in the fluid turbulence at higher frequency.The hydrogen has lower performance of the stack compared to helium since the hydrogen has the larger stack length (l 3 ) which absorb more acoustic power (W s ), leading to decrease in the 'COP s ' .Using the data given in Tables 1, 2 and 5 in the equation Q n = Q/(PuA 1 ), the variation of the stack diameter 'd 1 ' versus the resonance frequency 'f' is determined (Table 6).The 'd 1 ' is proportional to 'f' for both the hydrogen and helium.This is because the increase in the resonance frequency increases the gas turbulence which decreases the thermal penetration depth (Table 2) and hence the stack diameter increases with increase in the resonance frequency (Table 6).Using the normalized equations and the optimized values given in Table 5, the stack-length (l 3 ) and the stack-centre point (X) are determined.By using Eq. ( 7), the oscillating gas displacement amplitude (x 1 ) at the right-side to the cold end of the (12) parallel plates stack is determined.The length of the cold heat-exchanger (l 4 ) at the cold-side of the stack equal to the length 2x 1 and 4x 1 is determined (Table 6).
The cold heat-exchanger length equal to the ambient heat-exchanger is taken as: l 4 = 4x 1 = l 2 (Table 6).Further, the length between the speaker and the ambient heat-exchanger 'l 1 ' is determined.Using the stack normalized acoustic-power equation W s /(PuA 1 ) = W n and substituting the known values, the acoustic-power dissipation in to the stack 'W s ' is determined.Similarly, the acoustic-power dissipation in the cold heat-exchanger 'W cx ' and the ambient heat-exchanger 'W ax ' are determined using Eq. ( 11) by setting the normalized critical temperature (Γ) equal to zero.The values of the normalized equations l 4n , X 4n , l 2n and X 2n versus the resonance frequency are determined (Table 5).The frequency of 600 Hz for hydrogen and 400 Hz for helium are chosen to be optimum by considering the balance between the optimized design and performance parameters values, namely the COP s , d 1 , W s , W cx , and W ax as given in Table 6.

Design and analysis of divergent resonator
The quarter wave-length divergent hollow resonator is a rigid structure made of thin stainless-steel sheets.It connects the loudspeaker and the stack and heat exchangers system and the right-side geometry is the taper-divergenthemispherical (TDH) design (Fig. 1).The compact TDH design is chosen over the other geometries because of its lower acoustic dissipation losses and improved cooler power-density as discussed in the published literature [1,24,25].The taper length (L 2 = 20 mm), the diameter ratio (d r = 0.15) and the hemispherical radius (r h = 0.5d 1 ) are chosen for the final shape for both hydrogen and helium (Fig. 4).At this geometry of the resonator the acoustic power dissipation losses are minimal and which improves cooling power, coefficient of performance and powerdensity of the cooler as discussed elsewhere [4,8].The divergent section of the resonator is the buffer volume terminated with the hemispherical design (Fig. 4), simulates Table 6 The optimized specifications and performance parameters of the parallel plates stack and heat exchangers system at 2y = 4δ k versus working gas and frequency open end for the oscillating gas where the velocity is maximum [1,9,26].The stainless-steel resonator sheet with a thickness of 0.5 mm is self-sufficient to sustain the dynamic pressure amplitude of 10 bar gas pressure.The stainlesssteel sheet is the best choice compared to aluminium and Delrin materials to minimise the acoustic dissipation losses and gas leakage and threading problems associated during experimentation as discussed elsewhere [2,3].The stainless-steel material has the following properties: thermal conductivity (k) 14.06 Wm −1 K −1 , density (ρ) 7924.2 kgm −3 and specific heat (C p ) 445.91 Jkg −1 K −1 .The speaker modified with the gas back volume (called driver) which match the speaker frequency with the resonance frequency to increase the efficiency from 3% to about 50-90% possible [6,[27][28][29].
The gas back volume (Fig. 4) which match the loudspeaker frequency and the resonance frequency.It is of 600 Hz for hydrogen and 400 Hz for helium.The moving coil spring constant 's' of the loudspeaker [5] is where 'm' is the mass of the vibrating material.The gas back volume for the driver [5] is (16) where 'γ' is the index of compression or expansion of gas, 'P' is the oscillating gas pressure, 'A b ' is the area of cross section of back volume.In the literature it is found that the driver has maximum efficiency when the area of cross section of back volume (A b ), the vibrating diaphragm of driver (A d ) and the cross section of stack (A 1 ) are equal [30][31][32][33][34].In Table 7 the design and electro-mechanical parameters of the speaker with gas back volume for hydrogen and helium are presented.Substituting the data of Tables 1, 2 and 7 in Eqs. ( 16) and ( 17), the spring constant 's' and gas back volume 'V b ' and its length 'l b ' are determined for both hydrogen and helium (Table 7).The resonator full length of quarter-wavelength design 'L t ' measured from the driver surface and the hemispherical end for both the hydrogen and helium is given by For the taper-divergent-hemispherical design (Fig. 4), the angle of diverging section 'θ 2 ' varies inversely with its diverging length 'L 3 ' . Because of the better properties of hydrogen compared to helium (Tables 1, 2, 3 and 4), the ( 17)   hydrogen at 600 Hz has the smaller length 'L 3 ' and the angle 'θ 2 ' , and the smaller resonator length 'L t ' and the volume 'V T ' .This in turn leading to 2.1 times better power-density over helium at 400 Hz (refer Table 8).Thus, the compact hydrogen resonator has only the 0.6 times the acoustic-power absorption 'W r ' in to the internal surfaces of the resonator compared to the helium resonator (Table 8).
Using the values of Table 1 and 2 in Eq. ( 5) and ( 6) and the acoustic-power dissipation in the resonator internal surface area 'w r ' between the driver the and ambient heat-exchanger, and in the taper-divergent-hemispherical portion of the resonator with hydrogen and helium is given by [5] Hydrogen has the lower acoustic power dissipation compared to helium because of its compact design.The acoustic-power dissipation into the individual resonator component is determined by multiplying Eq. ( 19) with its internal surface area.The total acoustic-power input to the resonator 'W r ' is the sum of all the individual resonator component acoustic power dissipation.Excluding the heatexchangers (Refer Fig. 4), the stainless-steel sheet resonator and the stack length is assumed to be super insulated or placed under the vacuum vessel to avoid the heat leak losses as discussed elsewhere [1,3,35].Using the values given in Tables 6 and 8, the total acoustic-power dissipation 'W t ' into the parallel plates stack, heat-exchangers, and the resonator is given by The heat rejection rate of the ambient heat-exchanger 'Q r ' , performance of the cooler 'COP' , Carnot's 'COPC' and relative 'COPR' with hydrogen and helium for the 1 kW cooler is calculated using Eqs.( 21) - (24) given below.(19) w r = 0.25ρu 2  1 δ v ω + 0.25 The hydrogen offers minimum acoustic resistance since it has good thermo-physical properties over helium (refer Table 1).The hydrogen resonator has the smaller heat rejection rate 'Q r ' since it has the smaller total acoustic-power dissipation 'W t ' in the cooler components (Table 8) compared to helium.Thus, the hydrogen cooler resulting with the better coefficient of performance and power-density over helium because of the decrease in the acoustic-power dissipation in the heat exchangers (Table 6).The low-cost hydrogen, which is about 2.5 times less expensive compared to helium.Hence hydrogen may be considered as one of the best working gasses for thermoacoustic coolers with the necessary safety precautions to avoid the danger of explosion.

DeltaEC results and discussion
In this section the theoretically optimized quarter wavelength taper-divergent-hemispherical design resonator (Fig. 4) with hydrogen and helium is tested in the 'Design environment for low-amplitude thermo-acoustic Energy Conversion-DeltaEC' software [13] for validation.The software anticipates the physical performance of the thermoacoustic devices (coolers and engines) using the onedimensional (1D) computer program for linear analysis.The results are analysed by considering the assumptions, constraints and the limitations of the DeltaEC software as well as the theoretical design.The DeltaEC integrates the 1D Nikolaus Rott's governing thermoacoustic wave, heat, and acoustic-power equations of the oscillating gas for the model proposed by a designer in a computer programme [15,22,23].The software uses the geometric sequential segments of the thermoacoustic cooler in a programme, namely the loudspeaker, duct, ambient heat-exchanger, stack, cold heat-exchanger, taper, cones, and the compliance.A designer develops the Del-taEC program using the statements, namely the BEGIN, Insert, Kill, Append, VESPEAKER, Duct, HX, STKSLAB, CONE, COMPLIANCE, HARDEND, etc.These segments do not have the heat leak losses with the surrounding environment since they are assumed to be insulated perfectly by default including the driver using the INSU-LATE statement.The driver electric power input 'E' with (24) COPR = COP COPC the gas back volume is equal to the sum of the acousticpower output 'W t ' (Eq.( 20)) and its energy conversion loss.Therefore, in DeltaEC software the ambient heatexchanger rejects heat 'Q R ' which is equal to the sum 'E' and 'Q' (cooling-power) (refer Table 9).The effect of the temperature difference of the ambient and cold heatexchangers 'θ' at 38 K and 28 K by setting the target of 301 K for the ambient heat-exchanger and the effect of the cold heat-exchanger length 'l 4 ' equal to 2x 1 and 4x 1 is discussed.For both hydrogen and helium, the frequency of gas 'f' varies inversely with the cold heat-exchanger length 'l 4 ' irrespective of the temperature difference 'θ' .The average gas temperature 'T' , mean temperature of the gas across the heat exchangers 'T mg ' , and mean temperature of the heat exchangers 'T mx ' are presented in Table 9.
It is found that the 'T' > 'T mg ' and 'T mg ' > 'T mx ' for both hydrogen and helium which is given by T > T mg > T mx .
For the temperature difference of 38 K, the DeltaEC predicts 283 K for 'T mg ' and 282 K for 'T mx ' at the constant ambient heat exchanger temperature of 301 K for both hydrogen and helium.It jusifies the assumptions of taking the 'T mg ' one kelvin greater than the 'T mx ' that is T mg = (T mx + 1 K) in the theoretical design (Table 2).The driver electric power input 'E' and its efficiency 'η ea ' is proportional to the cold heat-exchanger length 'l 4 ' , and inversely proportional to the temperature difference 'θ' for both hydrogen and helium (Table 9).This is because the increase in the length of 'l 4 ' has the higher acoustic-power absorption 'W cx ' (Table 6) which in turn increases the driver electric power input 'E' .The driver efficiency 'η ea ' increases as the oscillating gas offers smaller resistance on the vibrating diapharagm of the driver.Also the increase in the length of 'l 4 ' decreases the cooling power 'Q' and the performance of the cooler (Eq.[22][23][24].It is because of the increase in the acoustic-power absorption in the cold heat-exchanger 'W cx ' .The cooling-power and performance of the cooler (Eq.22) decreases with increase in the temperature difference for both hydrogen and helium because of the decrease in the cold heat-exchanger temperature.But the relative coefficient of performance (Eq.24) increases with increase in the temperature difference 'θ' because of the decrease in the cold heat-exchanger temperature 'T cx ' for both hydrogen and helium.The dynamic acoustic pressure amplitude 'P a ' and the resonator components solid material temperature 'T s ' as a function of the distance 'x' measured from the driver surface along the resonator length for hydrogen and helium for the optimized TDH resonator (Fig. 4 and Table 8) is shown in Figs.5-8.
Referring to Figs. 5 and 6, the gas dynamic acoustic pressure amplitude 'P a ' at the driver surface where 'x' equal to zero is found to be 0.3 bar for both hydrogen and helium, which is same as that of the theoretical design (Table 2).The gas dynamic acoustic pressure amplitude remains the same along the length of the duct (100% porosity) and decreases in the parallel plates stack and heat-exchangers system (85% porosity) for both hydrogen and helium.It decreases inordinately in the small length 20 mm taper section because of the sudden decrease in the diameter ratio of 15% (Fig. 4), and goes negative in the remaining portion of the resonator (Figs. 5 and 6).The variation of the resonator components solid material temperature 'T s ' as a function of the distance 'x' measured from the driver surface along the resonator length for both hydrogen and helium are shown in Figs.7 and  8.There is a small drop in 'T s ' along the length of the 'l 2 ' since it is targeted to a constant temperature of 'T ax ' at 301 K.It decreases drastically along the stack length since it pumps heat in the stack region and remains constant in the remaining length as shown in Figs.7 and 8.

Conclusion
The effect of resonance frequency of hydrogen and helium on the stack and heat-exchangers parallel plates spacing and thickness, the wavelength and the performance of the stack is discussed in Sect. 2. The increase in resonance frequency decreases these four parameters (Tables 3, 4 and 6) for both the gases.It is resulting with the very thin stack-heat exchangers plate spacing and thickness (Table 3), and decrease the resonator length.The decrease in resonator length improves the power density of the cooler (Table 8).Referring to Table 6, with the increase in resonance frequency, hydrogen shows small increase in the stack diameter over helium since it has better thermo-physical properties for a cooler (Table 1).Hence the cooler is optimized to operate with hydrogen at 600 Hz and helium at 400 Hz.The choice of operating the cooler at these two frequencies is justified based on its higher performance of the parallel plates stack and the cooler, and the power-density (Tables 6 and 8).The hydrogen shows 4.4% increase in the cooler performance and 110% increase in the powerdensity compared to helium.
For the taper-divergent-hemispherical resonator design (Fig. 4), the DeltaEC predicts the small decrease

Fig. 4 A
Fig. 4 A schematic representation of the resonator with taper-divergent-hemispherical design

Fig. 5 Fig. 6
Fig. 5 Dynamic acoustic pressure amplitude of hydrogen as a function of gas position

Fig. 7
Fig. 7 Resonator solid material temperature of hydrogen as a function of gas position

Table 1
Thermo-physical properties of hydrogen and helium at 10 bar and 283 K

Table 3
Stack-heat exchangers plate spacing and thickness versus working gas and resonance frequency

Table 4
Derived parameters used in the design process

Table 5
Effect of the resonance frequency on Dimensionless (Normalized) design parameters optimized for hydrogen and helium

Table 7
The driver electro-mechanical parameters with hydrogen and helium

Table 8
The optimized quarter-wavelength resonator parameters for hydrogen and helium

Table 9
Simulation results of DeltaEC for the taper-divergent-hemispherical resonator design versus the working gas, temperature difference and cold heat-exchanger length