Numerical model for nonlinear standing waves and weak shocks in thermoviscous fluids.

Nonlinear standing waves in a one-dimensional tube are studied numerically by using a finite-difference algorithm. The numerical code models the acoustic field in resonators for homogeneous, thermoviscous fluids. Calculations are performed exclusively in the time domain, and all harmonic components are obtained by one resolution. The fully nonlinear differential equation is written in Lagrangian coordinates. It is solved without truncation. Effects of absorption are included. Displacement and pressure wave forms are calculated at different locations and results are shown for different excitation levels and tube lengths. Amplitude distributions along the resonator axis for every harmonic component are also evaluated. Simulations are performed for amplitudes ranging from linear to strongly nonlinear and weak shock. A very good concordance with classic experimental and analytical results is obtained.