Zaco intop1(p)/intop1(pabe.iomega*(w+eps))/(rho0*c0) "Specific acoustic impedance"
pext_1 at3_spatial(0,0,-1[m],pabe.p_t,'minc') "Acoustic pressure at 1 m"
prms sqrt(0.5*pext_1*conj(pext_1)) "RMS pressure at 1 m"
TVR 20*log10(prms/Vrms/1[uPa/V]) "Transmitting Voltage Response (TVR)"
pext_Zeval at3_spatial(0,0,Zeval,pabe.p_t,'minc') "Acoustic pressure at Zeval"
Ifront 0.5*pext_Zeval*conj(pext_Zeval)/(rho0*c0) "On-axis intensity at Zeval"
Ptot -intop1(pabe.I_bndx*pabe.nx+pabe.I_bndy*pabe.ny+pabe.I_bndz*pabe.nz) "Total radiated power"
Iave Ptot/(4*pi*Zeval^2) "Average intensity of monopole source at Zeval"
Di Ifront/Iave "Intensity directivity"
DI 10*log10(Di) "Directivity index of tonpilz transducer"
k0 2*pi*freq/c0 "Wave number"
pzt_stress F_pre/intop2(1) "Nominal compressive stress at the PZT"
a sqrt(intop3(1)/pi) "Radius of the equivalent flanged piston"
DI_fl_pist 10*log10((k0*a)^2/(1-2*besselj(1,2*k0*a)/(2*k0*a))) "Directivity index of flanged piston"