Zmic Zmic_real(freq)+i*Zmic_imag(freq) "Microphone impedance"
omega 2*pi*freq "Angular frequency"
mu mat1.def.eta(T0) "Dynamic viscosity"
gamma mat1.def.gamma "Ratio of specific heats"
Cp mat1.def.Cp(T0) "Heat capacity at constant pressure"
rho mat1.def.rho(pA,T0) Density
kcond mat1.def.k(T0) "Thermal conductivity"
c mat1.def.cs(T0) "Speed of sound"
k0 omega/c "Adiabatic wave number"
Z0 rho*c "Adiabatic characteristic impedance"
kv sqrt(-i*omega*rho/mu) "Viscous wave number"
kth sqrt(-i*omega*rho*Cp/kcond) "Thermal wave number"
Yv_in -besselj(2,kv*a_in)/besselj(0,kv*a_in) "Mean value of the scalar thermal field"
Yth_in -besselj(2,kth*a_in)/besselj(0,kth*a_in) "Mean value of the scalar viscous field"
Zc_in Z0/sqrt(Yv_in*(gamma-(gamma-1)*Yth_in)) "Mode characteristic impedance"
kc_sq_in k0^2*(gamma-(gamma-1)*Yth_in)/Yv_in "Square of mode wave number"
kc_in sqrt(kc_sq_in) "Mode wave number"
Yv_rt -besselj(2,kv*a_rt)/besselj(0,kv*a_rt) "Mean value of the scalar viscous field"
Yth_rt -besselj(2,kth*a_rt)/besselj(0,kth*a_rt) "Mean value of the scalar thermal field"
Zc_rt Z0/sqrt(Yv_rt*(gamma-(gamma-1)*Yth_rt)) "Mode characteristic impedance"
kc_sq_rt k0^2*(gamma-(gamma-1)*Yth_rt)/Yv_rt "Square of mode wave number"
kc_rt sqrt(kc_sq_rt) "Mode wave number"