Subfield: Optics

Thermoacoustic Engine

A thermoacoustic engine is a device that converts heat into acoustic power, or sound waves - a form of mechanical work. Like many other heat machines, it can be operated in reverse to become a refrigerator, using sound to pump heat from a cold to a hot reservoir. The high operating frequencies reduce heat conduction and eliminate the need for any working chamber confinement. Unlike many other engine types, the thermoacoustic engine has no moving parts except the working fluid itself.

The efficiencies of thermoacoustic machines are typically lower than other engine types, but they have advantages in set up and maintenance costs. This creates opportunities for renewable energy applications, such as solar-thermal power plants and utilization of waste heat. Our analysis will focus on the creation of acoustic energy within the system, ignoring the extraction or conversion for powering external devices.

Part A: Sound wave in a closed tube

Consider a thermally insulating tube of length $L$ and cross-sectional area $S$, whose axis lies along the $x$ direction. The two ends of the tube are located at $x=0$ and $x=L$. The tube is filled with an ideal gas and is sealed on both ends. At equilibrium, the gas has temperature $T_{0}$, pressure $p_{0}$ and mass density $\rho_{0}$. Assume that viscosity can be ignored and that the gas motion is only in the $x$ direction. The gas properties are uniform in the perpendicular $y$ and $z$ directions.

<img_4558>

Figure 1
Context question:
A.1 When a standing sound wave forms, the gas elements oscillate in the $x$ direction with angular frequency $\omega$. The amplitude of the oscillations depends on each element's equilibrium position $x$ along the tube. The longitudinal displacement of each gas element from its equilibrium position $x$ is given by

$$
u(x, t)=a \sin (k x) \cos (\omega t)=u_{1}(x) \cos (\omega t)
\tag{1}
$$

(please note the $u$ here describes the displacement of a gas element)

where $a \ll L$ is a positive constant, $k=2 \pi / \lambda$ is the wavenumber and $\lambda$ is the wavelength. What is the maximum possible wavelength $\lambda_{\max }$ in this system?
Context answer:
\boxed{$\lambda_{\max }=2 L$}


Extra Supplementary Reading Materials:

We will assume throughout the question an oscillation mode of $\lambda=\lambda_{\max }$.

Now, consider a narrow parcel of gas, located at rest between $x$ and $x+\Delta x(\Delta x \ll L)$. As a result of the displacement wave of Task A.1, the parcel oscillates along the $x$ axis and undergoes a change in volume and other thermodynamic properties.

Throughout the following tasks assume all these changes to the thermodynamic properties to be small compared to the unperturbed values.
Context question:
A.2 The parcel volume $V(x, t)$ oscillates around the equilibrium value of $V_{0}=S \Delta x$ and has the form

$$
V(x, t)=V_{0}+V_{1}(x) \cos (\omega t) .
\tag{2}
$$

Obtain an expression for $V_{1}(x)$ in terms of $V_{0}, a, k$ and $x$.
Context answer:
\boxed{$V_{1}(x)=a k V_{0} \cos (k x)$}


Context question:
A.3 Assume that the total pressure of the gas, as a result of the sound wave, takes the approximate form

$$
p(x, t)=p_{0}-p_{1}(x) \cos (\omega t) .
\tag{3}
$$

Considering the forces acting on the parcel of gas, compute the amplitude $p_{1}(x)$ of the pressure oscillation to leading order, in terms of the position $x$, the equilibrium density $\rho_{0}$, the displacement amplitude $a$ and the wave parameters $k$ and $\omega$.
Context answer:
\boxed{$p_{1}(x)=a \frac{\omega^{2}}{k} \rho_{0} \cos (k x)$}


Extra Supplementary Reading Materials:

At acoustic frequencies, the thermal conductivity of the gas can be neglected. We will treat the expansion and contraction of gas parcels as purely adiabatic, satisfying the relation $p V^{\gamma}=$ const, where $\gamma$ is the adiabatic constant.
Context question:
A.4 Use the relation above and the results of the previous tasks to obtain an expression for the speed of sound waves $c=\omega / k$ in the tube, to first order. Express your answer in terms of $p_{0}, \rho_{0}$ and the adiabatic constant $\gamma$.
Context answer:
\boxed{$c=\sqrt{\frac{\gamma p_{0}}{\rho_{0}}}$}


Context question:
A.5 The change in the gas temperature due to the adiabatic expansion and contraction, as a result of the sound wave, takes the form:

$$
T(x, t)=T_{0}-T_{1}(x) \cos (\omega t)
\tag{4}
$$

Compute the amplitude $T_{1}(x)$ of the temperature oscillations in terms of $T_{0}, \gamma$, $a, k$ and $x$.
Context answer:
\boxed{$a k(\gamma-1) T_{0} \cos (k x)$}


Context question:
A.6 For the purpose of this task only, we assume a weak thermal interaction between the tube and the gas. As a result, the standing sound wave remains almost unchanged, but the gas can exchange a small amount of heat with the tube. The heating due to viscosity can be neglected.

For each of the points in Figure 2 (A, C at the edges of the tube, B at the center) state whether the temperature of the tube at that point will increase, decrease or remain the same over a long time.

<img_4517>

Figure 2
Context answer:
\boxed{证明题}


Extra Supplementary Reading Materials:

Part B: Sound wave amplification induced by external thermal contact

A stack of thin well-spaced solid plates is placed inside the tube. The plates of the stack are aligned in parallel to the tube axis, so as not to obstruct the flow of gas along the tube. The center of the stack is positioned at $x_{0}=L / 4$, and spans a width of $\ell \ll L$ along the tube axis, filling its entire cross section. The right and left edges of the stack are held at temperature difference $\tau$. The left edge of the stack, at $x_{H}=x_{0}-\ell / 2$, is held by an external thermal reservoir at temperature $T_{H}=T_{0}+\tau / 2$, and at the same time, its right edge, at $x_{C}=x_{0}+\ell / 2$, is held at a temperature $T_{C}=T_{0}-\tau / 2$.

The plate stack allows a slight longitudinal heat flow to maintain a constant temperature gradient between its edges, such that $T_{\text {plate }}(x)=T_{0}-\frac{x-x_{0}}{\ell} \tau$.

<img_4412>

Figure 3. A sketch of the system. (A) and (B) denote the hot and cold heat reservoirs respectively. (D) denotes the stack.

To analyze the effect of the thermal contact between the plate stack and the gas on the sound waves in the tube, make the following assumptions:

- As in the previous part, all changes to the thermodynamic properties are small compared to the unperturbed values.
- The system operates in the fundamental standing-wave mode of the longest possible wavelength. It is only slightly modified by the presence of the plate stack.
- The stack is much shorter than the wavelength $\ell \ll \lambda_{\max }$, and can be positioned far enough from both displacement and pressure nodes, so that the displacement $u(x, t) \approx u\left(x_{0}, t\right)$ and the pressure $p(x, t) \approx p\left(x_{0}, t\right)$ may be considered uniform over the entire length of the stack.
- We may neglect any edge effects, caused by the parcels moving in and out of the stack.
- The temperature difference between the ends of the plate stack, i.e. between the hot and the cold reservoirs, is small compared to the absolute temperature: $\tau \ll T_{0}$.
- Heat conduction through the stack, through the gas, and along the tube are all negligible. The only significant sources of heat transfer are convection due to the motion of the gas and conduction between the gas and the stack.
Context question:
B.1 Consider a specific parcel of gas in the region of the stack, originally at $x_{0}=L / 4 .$ As the parcel moves within the stack, the local temperature of the nearby part of the stack changes as follows:

$$
T_{\mathrm{env}}(t)=T_{0}-T_{\mathrm{st}} \cos (\omega t)
\tag{5}
$$

Express $T_{\mathrm{st}}$ in terms of $a, \tau$ and $\ell$.
Context answer:
\boxed{$\frac{a \tau}{\ell \sqrt{2}}$}


Context question:
B.2 Above which critical temperature difference $\tau_{\text {cr }}$ will the gas be conveying heat from the hot reservoir to the cold one? Express $\tau_{\mathrm{cr}}$ in terms of $T_{0}, \gamma, k$ and $\ell$.
Context answer:
\boxed{$\tau_{\mathrm{cr}}=k \ell(\gamma-1) T_{0}$}


Context question:
B.3 Obtain the general approximate expression for the heat flow $\frac{d Q}{d t}$ into a small parcel of gas as a linear function of its volume and pressure change rates. Express your answer in terms of the rate of volume change $\frac{d V}{d t}$, the rate of pressure change $\frac{d p}{d t}$, the unperturbed equilibrium values of parcel pressure and volume $p_{0}, V_{0}$ and the adiabatic index $\gamma$. (You may use the expression for the molar heat capacity at constant volume $c_{v}=\frac{R}{\gamma-1}$, where $R$ is the gas constant.)
Context answer:
\boxed{$\frac{1}{\gamma-1} V_{0} \frac{d p}{d t}+\frac{\gamma}{\gamma-1} p_{0} \frac{d V}{d t}$}


Extra Supplementary Reading Materials:

The limited heat flow rate between the parcel and the stack causes a phase difference between the pressure and volume oscillations of the parcel. We will see how this generates work.

Let the heat flux into the parcel from the stack be proportional to the temperature difference between the parcel and the neighboring element of the stack, given approximately by $\frac{d Q}{d t}=-\beta V_{0}\left(T_{\mathrm{st}}-T_{1}\right) \cos (\omega t)$. Here $T_{1}$ and $T_{\mathrm{st}}$ are the temperature oscillation amplitudes of the gas parcel and the neighbouring stack from Tasks A. 5 and B. 1 , respectively, and $\beta>0$ is a constant. Assume that at the machine's operating frequencies, the change in gas temperature as a result of this heat flow is insignificant compared to both $T_{1}$ and $T_{\mathrm{st}}$.
Context question:
B.4 In order to calculate work, we will consider a change to the volume of the moving parcel as a result of the thermal contact with the stack. Let us write the pressure and the volume of the parcel under the stack's influence in the form:

$$
\begin{gathered}
p=p_{0}+p_{a} \sin (\omega t)-p_{b} \cos (\omega t), \\
V=V_{0}+V_{a} \sin (\omega t)+V_{b} \cos (\omega t) .
\end{gathered}
\tag{6}
$$

Given $p_{a}$ and $p_{b}$, find the coefficients $V_{a}$ and $V_{b}$. Express your answer in terms of $p_{a}, p_{b}, p_{0}, V_{0}, \gamma, \tau, \tau_{\text {cr }}, \beta, \omega, a$ and $\ell$.
Context answer:
\boxed{$V_{b}=\frac{1}{\gamma} p_{b} \cdot \frac{V_{0}}{p_{0}}$ , $V_{a}=(-\frac{1}{\gamma} p_{a}-\frac{\gamma-1}{\gamma} \frac{\beta}{\omega} \frac{a}{\ell \sqrt{2}}(\tau-\tau_{\mathrm{cr}})) \cdot \frac{V_{0}}{p_{0}}$}


Context question:
B.5 Obtain an approximate expression for the acoustic work per unit volume $w$ produced by the gas parcel over one cycle. Integrate over the volume of the stack to obtain the total work $W_{\text {tot }}$ generated by the gas over one cycle. Express $W_{\text {tot }}$ in terms of $\gamma, \tau, \tau_{\text {cr }}, \beta, \omega, a, k$ and $S$.
Context answer:
\boxed{$\frac{\pi}{2 \omega}(\gamma-1) \beta(\tau-\tau_{\mathrm{cr}}) k a^{2} S$}


Question: B.6 Obtain an approximate expression for the heat $Q_{\text {tot }}$ transported from the left side of the plane $x=x_{0}$ to the right, over a cycle. Express your answer in terms of $\tau, \tau_{\text {cr }}, \beta, \omega, a, S, \ell$

Solution:
We want to compute the amount of heat convection over one cycle. This means that we need to take the amount of heat moving in or out of the parcel, and weigh it by the position of the parcel at that time. Thus, the total heat conveyed by the parcel, integrated along a cycle, is:

$$
Q_{\text {tot }}=\frac{1}{\Delta x} \int \frac{d Q}{d t} u \cdot d t
$$

This expression can be computed to leading order using $\frac{d Q}{d t}=\beta V_{0}\left(T_{\mathrm{st}}-T_{1}\right) \cdot \cos (\omega t)$ and the unperturbed displacement $u\left(x_{0}, t\right)=\frac{a}{\sqrt{2}} \cos (\omega t)$. This gives

$$
Q_{\mathrm{tot}}=\frac{\pi}{\omega} \beta V_{0}\left(T_{\mathrm{st}}-T_{1}\right) \frac{a}{\sqrt{2}}=\frac{\pi}{\omega} \beta V_{0} \cdot \frac{a}{\ell \sqrt{2}}\left(\tau-\tau_{\mathrm{cr}}\right) \cdot \frac{a}{\sqrt{2}}=\frac{\pi}{2 \omega} \beta\left(\tau-\tau_{\mathrm{cr}}\right) \frac{a^{2} S}{\ell} .
$$

Final answer: $\frac{\pi}{2 \omega} \beta(\tau-\tau_{\mathrm{cr}}) \frac{a^{2} S}{\ell}$