A quasi-one-dimensional model of thermoacoustics in the presence of mean flow
Introduction
The interest in the interaction of acoustics and heat flow dates back more than 200 years [1], [2]. One of the investigated issues is the thermoacoustic (TA) effect: in narrow geometries thermo-viscous boundary layers interact with axial mean temperature gradients. This interaction causes a conversion of thermal energy to acoustic power or vice versa. Sondhauss [3] observed this effect first. He described the occurrence of spontaneous loud sound emissions in narrow pipes during the glass blowing process. Although Rayleigh [4] had already given the first qualitative explanations for this effect, only incomplete descriptions of the phenomenon were provided by different authors [5], [6]. Rott was the first author to completely describe the phenomenon analytically in a series of publications [7], [8], [9], [10], [11], [12]. Wheatley [13], his successor Swift [14] and finally in׳t panhuis [15] and in׳t panhuis et al. [16] improved the analytical understanding of the TA effect. Arnott et al. [17] enhanced the applicability of this approach for arbitrary pore shapes. Watanabe et al. [18] further improved these models by accounting for nonlinear effects such as drag or heat transfer. The TA effect has also been investigated in various computational fluid dynamics (CFD) simulations (e.g. [19], [20], [21], [22], [23]). Their analytical formulations facilitate the description of the axial propagation of acoustic velocity u1 and pressure p1 in terms of their cross-section averaged quantities. These transversally averaged parameters form a set of transport equations in the axial pore direction. Those transport equations are applied in a one-dimensional numerical tool [24] to predict the operating conditions in TA devices.
Previous publications have presented ideas for probable applications of TA devices [13], [14], and some authors claim that the technology is ripe for the market [25]. Nevertheless, almost 30 years after the first attempts of technical realizations of such apparatus, the technology is only applied in niche applications [26], [27], [28]. Thus, the main aim of researchers remains to expand the applicability of TA for cooling and power generation.
All cited authors restrict their descriptions to quiescent mean flow conditions. Until now, only a few publications have tried to deal with mean flow affected TA power conversion [29], [30], [31]. Bammann et al. [32] provide an overview over the published literature on flow-through TA refrigeration. One of the main reasons for the modest interest in their topic is the lack of proper modeling of the TA interaction effects under such conditions. All hitherto referenced publications consider the impact of mean flow only implicitly by adapting the axial mean temperature profile, or do not address it at all. Based on the assumption that the mean flow is due to streaming effects induced by acoustic quantities, a second-order enthalpy flux consideration leads to a distortion of the mean temperature profile [8], [33], [34], [35]. The present paper presents and analyzes an improved quasi-one-dimensional method, which explicitly takes the laminar mean flow into account in the perturbation equations. Using this approach enables the prediction of the performance of TA devices that exploit mean flow as a positive effect, e.g. Reid׳s refrigerator [30]. Due to the resulting one-dimensional form, the computational effort for the calculation of one configuration is comparatively small. Conceivable applications may, for example, arise in exhaust gas streams, where the waste heat enthalpy flux is directly transformed into acoustic power.
The derivation of the quasi-one-dimensional model is based on the method applied by in׳t panhuis et al. [16], [15], who derived TA transport equations that account for almost arbitrary stack pore geometry. Like all other proposed approaches, the model takes into account the so-called long-pore assumption, i.e. the hydraulic diameter is much smaller than the characteristic length in axial direction L. Further, in׳t panhuis [15] follows the common approach of neglecting the contribution from mean flow when linearizing the basic set of equations. In the present work, this assumption is not invoked and mean flow is explicitly incorporated.
Taking into account the contribution of mean flow leads to an increase in the number of terms in the considered equations. Hence, for the purpose of demonstrating the impact of explicitly accounting for mean flow, the solution is derived for a generic problem resulting in the least complex system of transport equations. The geometrical conditions and thermo-physical correlations are chosen to keep the equations short. Sticking to the general solution with only a few necessary assumptions is still feasible; however, this complicates the interpretation of the impact of mean flow due to the huge numbers of terms inside the system of equations. Thus, to simplify the model where possible, assumptions will be introduced. These assumptions due to compactness of equations are marked as such.
Apart from spatial distributions of the acoustic quantities for a single frequency, acoustic scattering matrices are used for validating the one-dimensional models against CFD results. The latter validation tool is mostly used to investigate the stability of the system in the linear regime by using 1D network models [36]. Although TA engines generally operate in nonlinear limit cycles, the identified growth rate at the onset of the acoustic instability is a quantitative measure for the achievable pressure amplitudes. Altogether, as the simplicity of the solution is a central request, this publication is restricted to the linear acoustic regime. If nonlinearities like the acoustic feedback into the mean temperature distribution or saturation effects have to be taken into account, all steps can be processed for the nonlinear situation.
The paper is organized as follows: In Section 2, the geometrical and mean flow conditions are presented. In Sections 3 and 4, the one-dimensional acoustic transport equations are derived for a parallel plate stack pore. The flow chart of the derivation is sketched in Fig. 1: First the compressible Navier–Stokes equations (NSE) are non-dimensionalized using narrow geometry assumptions. These dimensionless NSEs experience a simplified modal expansion. The mean flow field is solved analytically from the resulting zeroth-order set of equations. These solutions are used as input to the first-order set of equations, the linearized dimensionless NSEs. This set of equations for the acoustic field is then solved in segregated steps: First Green׳s function (GF) solution technique is applied on the axial momentum and energy equations. This yields formulations for their transversal shapes. Closure assumptions have to be made for two convective terms to obtain an analytical solution of these formulations. Integrating the resulting equations in the transversal direction creates a set of averaged perturbation equations. Finally, eliminating density and temperature perturbations leads to one-dimensional TA transport equations in terms of velocity and pressure oscillations. The resulting set of equations accounts for non-zero mean flow contributions, which appear explicitly in some terms. Thus, models based on this system of equations are named “explicit” in the rest of the paper. The existing model [37] is denoted as “implicit” (IMP), because in that model mean flow influences the system of equations only implicitly through its impact on the mean temperature profile, which is computed from a second-order energy equation that considers enthalpy transport by the mean flow [8]. For the explicit equations two alternative approaches for the necessary closure assumptions are investigated. Prediction accuracy of these two approaches, denoted by M I and M II, is tested for the spatial acoustic propagation for a single frequency and scattering behavior for different test cases against a full CFD simulation and the existing implicit one-dimensional solution method. Finally, an outlook to improved modeling approaches is presented.
Section snippets
Geometric settings
A channel of length L and constant hydraulic radius R is one of the most simple geometries to be considered for a TA system. Restricting the problem to a constant pore height simplifies the contact conditions of the fluid and solid regions in the problem. Choosing a channel geometry leads to a two-dimensional problem in space, which enormously reduces the size of the derived equations. The solid plates in the problem sketched in Fig. 2 are thick. The origin of the Cartesian coordinate
Non-dimensional, linearized Navier–Stokes equations of a slab pore
This section covers the first four steps of the sketch depicted in Fig. 1: At first the Navier–Stokes equations are non-dimensionalized. The resulting set of equations is then split into transport equations describing the mean flow field and the propagation of acoustic quantities. Here, a special modal expansion technique is applied.
As boundary layer effects are the cause of the TA effect, the Navier–Stokes equations are a suitable choice of equations describing the motion of the fluid under
Transversally averaged acoustic perturbation equations
The cross-sectional averaging of the LNSEs of a pore (Eqs. (26a), (26b), (26c), (26d), (26e), (26f)) leads to a system of one-dimensional, coupled differential equations. This set can then be evaluated numerically with little effort. Nevertheless, the averaged equations contain numerous terms. To avoid human errors in this derivation step, the process was carried out with a computational software capable of symbolic equation handling [42]. The derivation of the intermediate and final
Spatial propagation and scattering matrices
The accuracy of the closure assumptions chosen is validated against an existing one-dimensional solution denoted by IMP and a full CFD computation of one single pore. At first the x distributions for a single frequency are compared. Then, the frequency dependent scattering behavior is investigated in terms of scattering matrices.
Spatial propagation for a single frequency
The solid lines depicted in Fig. 4 represent the axial profiles of amplitude of the acoustic pressure (a) and the velocity fluctuation (b). Two cases are investigated: quiescent conditions (Ma0) and the mean flow affected case . The displayed phase Φ (dashed with empty symbols) is displayed relative to the reference phase of the acoustic pressure p1 located at the upstream end of the pore.
All boundary conditions are chosen to be equal except for BC of the mean upstream velocity.
Further improvements
The system of equations (58) derived in Section 4 used the two simplest modeling approaches, i.e. neglecting the convective terms for M I and inserting a spatial averaged formulation for quiescent mean conditions for M II. Furthermore, the mean parameters were restricted to be constant in the transversal direction. This procedure kept the expressions short, but both simplifications may be replaced by more accurate methods.
If laminar mean flow problems are considered, a parabolic velocity and
Conclusions
The one-dimensional explicit models derived in this study improve the accuracy of the description of the thermo-viscous interaction in narrow pores with mean flow. Earlier implicit models only incorporate the change in mean temperature distribution, which then interacts with the acoustic variables. The resulting scattering predictions of these implicit one-dimensional approaches are not affected by different, mean flow controlled temperature shapes. If mean flow is explicitly considered in the
References (65)
Vibrations of a gas column
Physica
(1949)Thermoacoustics
Archive of Applied Mechanics
(1980)- et al.
Experimental and computational visualization of the flow field in a thermoacoustic stack
Comptes Rendus Mécanique
(2003) - et al.
CFD simulation of thermoacoustic cooling
International Journal of Heat and Mass Transfer
(2010) - et al.
CFD simulation of a thermoacoustic engine with coiled resonator
International Communications in Heat and Mass Transfer
(2010) - et al.
A transient testing technique for the determination of matrix parameters of acoustic systems. Part Itheory and principles
Journal of Sound and Vibration
(1979) - et al.
Theory of a two source-location method for direct experimental evaluation of the four-pole parameters of an aeroacoustic element
Journal of Sound and Vibration
(1990) Time dependent boundary conditions for hyperbolic systems
Journal of Computational Physics
(1987)- et al.
A nonreflecting outflow boundary condition for subsonic Navier–Stokes calculations
Journal of Computational Physics
(1980) - et al.
Reconstruction of acoustic transfer matrices by instationary computational fluid dynamics
Journal of Sound and Vibration
(2001)
On the sound produced by a current of hydrogen gas passing through a tube
Journal of Natural Philosophy and Chemical Arts
Notiz über eine neue Art, die in einer an beiden Enden offenen Röhre enthaltene Luft in Schwingungen zu versetzen (note on a new method to cause an oscillation of air at both ends of an open pipe)
Annalen der Physik und Chemie
Über die Schallschwingungen der Luft in erhitzten Glasröhren und in gedeckten Pfeifen von ungleicher Weite on (acoustic oscillations of air in heated glas pipes and capped whistles of dissimilar width)
Annalen der Physik
Über den Einfluss der Wärmeleitung in einem Gase auf die Schallbewegung (on the influence of heat conduction in a gas on the acoustic movement)
Annalen der Physik
Damped and thermally driven acoustic oscillations in wide and narrow tubes
Zeitschrift für Angewandte Mathematik und Physik
The heating effect connected with non-linear oscillations in a resonance tube
Zeitschrift für Angewandte Mathematik und Physik
Thermally driven acoustic oscillations. Part IIIsecond-order heat flux
Zeitschrift für Angewandte Mathematik und Physik
Thermally driven acoustic oscillations. Part IVtubes with variable cross-section
Zeitschrift für Angewandte Mathematik und Physik
The description of simple waves by particle displacement
Zeitschrift für Angewandte Mathematik und Physik
Natural engines
Physics Today
Weakly nonlinear thermoacoustics for stacks with slowly varying pore cross-sections
Journal of Fluid Mechanics
General formulation of thermoacoustics for stacks having arbitrarily shaped pore cross sections
Journal of the Acoustical Society of America
A simplified model for linear and nonlinear processes in thermoacoustic prime movers. Part I. Model and linear theory
Journal of the Acoustical Society of America
Study of nonlinear processes of a large experimental thermoacoustic-stirling heat engine by using computational fluid dynamics
Journal of Applied Physics
Low mach number analysis of idealized thermoacoustic engines with numerical solution
Journal of the Acoustical Society of America
Le Froid Venu du Son (cooling by sound)
La Gazette du Vide
Cited by (2)
Dynamics of unsteady heat transfer in pulsating flow across a cylinder
2017, International Journal of Heat and Mass TransferCitation Excerpt :One novelty of the present paper is the combined use of CFD and SI to describe the dynamics of heat transfer in response to flow perturbations. Previously, the CFD/SI method was applied to problems in flame dynamics [51], aero-acoustics [52,53] and thermoacoustic stacks [54,55]. SI techniques were also used to identify approximate linear models for separated flows to compute feedback control laws, see e.g. [56,57].
Common features in the thermoacoustics of flames and engines
2018, International Journal of Spray and Combustion Dynamics