k0 omega/c "Adiabatic wave number"
Z0 rho0*c "Adiabatic characteristic impedance"
kv sqrt(-i*omega*rho0/mu) "Viscous wave number"
kth sqrt(-i*omega*rho0*Cp/kcond) "Thermal wave number"
Yv -besselj(2,kv*a)/besselj(0,kv*a) "Mean value of the scalar thermal field"
Yth -besselj(2,kth*a)/besselj(0,kth*a) "Mean value of the scalar viscous field"
Zc Z0/sqrt(Yv*(gamma-(gamma-1)*Yth)) "Mode characteristic impedance"
kc_sq k0^2*(gamma-(gamma-1)*Yth)/Yv "Square of mode wave number"
kc sqrt(kc_sq) "Bulk wave number"
Y_tl (i*S*sin(kc*L)/Zc) "Transmission line admittance"