P0 5e6[Pa] Source pressure amplitude v0 P0/(rho0*c0) Source velocity amplitude f0 1e5[Hz] Source frequency T0 1/f0 Source period omega0 2*pi*f0 Angular frequency beta 10 Coefficient of nonlinearity c0 1481.44[m/s] Speed of sound rho0 999.6[kg/m^3] Density mu0 1.0016e-3[Pa*s] Viscosity d_diff 4/3*mu0/rho0 Diffusivity of sound Gamma 2*beta*P0/(omega0*rho0*d_diff) Goldberg number x_sh rho0*c0^3/(beta*P0*omega0) Shock formation distance L 4.5*x_sh Computational domain length Nt ceil(L/(c0*T0)) Number of cycles N0 8 Number of harmonics to resolve