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)