L 0.5[m] Domain length f 1e5[Hz] Source frequency P0 1e6[Pa] Source amplitude T 1/f Source period beta 10 Coefficient of nonlinearity omega 2*pi*f Angular frequency c 1481.44[m/s] Speed of sound k omega/c Wave number rho 999.6[kg/m^3] Density NT L/c/T Number of cycles Xsh (rho*c^3)/(beta*P0*omega) Shock formation distance Tend (NT-1)*T Total traveling time to simulate Nh 4 Number of harmonics to resolve dx0 c/f/6 Mesh size resolving actuation frequency using quartic elements dx dx0/Nh Mesh size resolving up to Nh'th harmonic dt 0.2*dx/c CFL condition on time step (quadratic elements)