The present research analyzes random porous thermoacoustic stack systems analytically, experimentally, and numerically with a primary objective to develop a comprehensive analytical porous media modeling for random porous (such as Reticulated Vitreous Carbon (RVC) foams) environment. Mathematical models are developed for flow, thermal, and energy fields within the random porous medium stack. The Darcy and Brinkman-Forchheimer-extended Darcy models are used for modeling the momentum equation and local thermal equilibrium assumption between the porous matrix and trapped fluid in the void space for energy equation. The expressions of temperature, energy flux density, and acoustic work absorbed or produced by a thermoacoustic device are compared with existing literature and observed good agreements. After obtaining the flow and thermal fields’ information, the present study examines the entropy generation distribution within the stack. One important item revealed in this study is that entropy generation inside the porous medium completely follows the trend of the imaginary part of Rott’s first function profile. Another major contribution of this research is to identify the location of maximum entropy generation which is identical to the location of maximum thermoacoustic heat and work transport. The expression of Nusselt number for steady flow cannot be used in oscillatory random porous medium because of the phase difference between the temperature gradient at the wall and the temperature difference between the wall and the space averaged temperature. The present research experimentally examines novel stack configuration by considering “alternating conducting and insulating materials” as stack in thermoacoustic devices. The objective of considering such stack arrangement is to reduce the conduction heat transfer loss from the hot end of the stack to the cold end, thereby increasing the performance of the stack. Eight different heterogeneous stack arrangements are studied in this research. The performance of the heterogeneous stack arrangement is compared with the typical homogeneous stacks. This research shows that heterogeneous stacks can be used in thermoacoustic devices particularly in small (millimeter) scale thermoacoustic devices. Numerically the present study investigates the influence of working fluid, geometric, and operating conditions on stack performance by solving the full Navier-Stokes, mass, energy equation, and equation of state.

The present research analyzes random porous thermoacoustic stack systems analytically, experimentally, and numerically with a primary objective to develop a comprehensive analytical porous media modeling for random porous (such as Reticulated Vitreous Carbon (RVC) foams) environment. Mathematical models are developed for flow, thermal, and energy fields within the random porous medium stack. The Darcy and Brinkman-Forchheimer-extended Darcy models are used for modeling the momentum equation and local thermal equilibrium assumption between the porous matrix and trapped fluid in the void space for energy equation. The expressions of temperature, energy flux density, and acoustic work absorbed or produced by a thermoacoustic device are compared with existing literature and observed good agreements. After obtaining the flow and thermal fields’ information, the present study examines the entropy generation distribution within the stack. One important item revealed in this study is that entropy generation inside the porous medium completely follows the trend of the imaginary part of Rott’s first function profile. Another major contribution of this research is to identify the location of maximum entropy generation which is identical to the location of maximum thermoacoustic heat and work transport. The expression of Nusselt number for steady flow cannot be used in oscillatory random porous medium because of the phase difference between the temperature gradient at the wall and the temperature difference between the wall and the space averaged temperature. The present research experimentally examines novel stack configuration by considering “alternating conducting and insulating materials” as stack in thermoacoustic devices. The objective of considering such stack arrangement is to reduce the conduction heat transfer loss from the hot end of the stack to the cold end, thereby increasing the performance of the stack. Eight different heterogeneous stack arrangements are studied in this research. The performance of the heterogeneous stack arrangement is compared with the typical homogeneous stacks. This research shows that heterogeneous stacks can be used in thermoacoustic devices particularly in small (millimeter) scale thermoacoustic devices. Numerically the present study investigates the influence of working fluid, geometric, and operating conditions on stack performance by solving the full Navier-Stokes, mass, energy equation, and equation of state.

Thermoacoustic devices are a group of systems that make use of the thermoacoustic effect to achieve an energy conversion between thermal and acoustic energy. The thermoacoustic effect occurs when a solid boundary is introduced into an acoustic field, and a non-zero net heat transportation takes place while the net mass transfer remains null. Thermoacoustic technologies are gaining an increasing research interest because of their potential applications for building alternative prime movers or heat pumps which do not use working fluids causing environmental damage and require very little maintenance due to their lack of moving part. However, the operation of this type of system is yet to be fully understood: fluid flow and heat transfer processes within the system components such as thermoacoustic stacks and heat exchangers still require a lot of attention. The performance of the system working with relatively low amplitude acoustic wave can be predicted by the linear thermoacoustic theory, which is already well developed. However, a high amplitude acoustic wave is usually required in order to achieve high power density or high power output. Unfortunately, the performance of such systems can be seriously degraded due to nonlinear effects, such as turbulence, minor loss or high proportion of harmonics. The lack of understanding of these effects impedes the design and construction of high efficiency systems. The work described in this thesis is focused on the study of flow phenomena taking place around parallel plate stack placed in a standing wave thermoacoustic resonator, by using advanced flow diagnostics techniques such as particle image velocimetry (PIV) and hot wire anemometry (HWA). In order to carry out the experimental study, a standing wave thermoacoustic device working at relatively low frequency of 13.1Hz was designed, commissioned and tested. The frequency response of this device was carefully investigated and compared with the analytical results using linear acoustic equations and a linear model of the loudspeaker. A further comparison with the analytical results obtained with the modelling tool DeltaEC (Design Environment for Low-amplitude Thermoacoustic Energy Conversion) was also presented. The resonator was driven from low to large pressure amplitudes with drive ratios up to 10%. A good agreement is obtained for small amplitudes, but the discrepancies become larger when the driving amplitude is increased. The analysis reveals that the large discrepancy at high amplitude can be attributed to minor losses. Following the above preliminary work, a more comprehensive study of the flow field around parallel-plate stacks was conducted by means of PIV and HWA. It was shown that the flow around the two studied parallel-plate stacks exhibits rather complicated flow features when the amplitude of the acoustic oscillation varies. Symmetrical and asymmetrical vortex shedding phenomena are observed and two distinct modes of generating 'vortex streets' are identified. It shown that a velocity related parameter such as the Reynolds number, defined on the plate thickness and the velocity amplitude at the entrance to the stack, and a geometrical parameter are not sufficient to define the flow characteristics in this type of flow problem. It is also proposed to introduce an extra frequency related parameter such as the Keulegan-Carpenter number (KC) and to carry out a similarity analysis in order to understand better the physics behind the flow phenomena and their controlling parameters. Typical ensemble-averaged velocity fields are used in the analysis above. However, the detailed flow features obtained from the ensemble averaged flow fields and the instantaneous flow fields could be different in a substantial way. The flow behaviour, its kinematics, dynamics and scales of turbulence, therefore are further investigated by using the classical Reynolds decomposition to separate the instantaneous velocity fields into ensemble-averaged mean velocity fields and fluctuations in a set of predetermined phases within an oscillation cycle. The mean velocity field and the fluctuation intensity distributions are investigated over the acoustic oscillation cycle. By using fast Fourier transform (FFT) spatial filtering techniques, the velocity fluctuation is further divided into large- and small-scale fluctuations, and their physical significance is discussed. The physics behind the flow phenomena are further studied by carrying out an analysis of the wake flow during the ejection part of the flow cycle, where either closed re-circulating vortices or alternating vortex shedding can be observed. A similarity analysis of the governing Navier-Stokes equations is then undertaken in order to derive the similarity criteria governing the wake flow behaviour. Similarity numbers including two types of Reynolds number, the KC number and a non-dimensional stack configuration parameter are considered. The influence of these parameters on the flow behaviour is discussed by investigating the experimental data obtained, along with additional data from literature.