  The formation of thermoacoustic shocks is revealed in a fluid complex plasma.
The thermoacoustic wave mode can be damped (or anti-damped) when the
contribution from the thermoacoustic interaction is lower (or higher) than that
due to the particle collision and/or the kinematic viscosity. In the nonlinear
regime, the thermoacoustic wave, propagating with the acoustic speed, can
evolve into small amplitude shocks whose dynamics is governed by the
Bateman-Burgers equation with nonlocal nonlinearity. The latter can cause the
shock fronts to be stable (or unstable) depending on the collision frequency
remains below (or above) a critical value and the thermal feedback is positive.
The existence of different kinds of shocks and their characteristics are
analyzed with the system parameters that characterize the thermal feedback,
thermal diffusion, heat capacity per fluid particle, the particle collision and
the fluid viscosity. A good agreement between the analytical and numerical
results are also noticed.
