Comparing the equation for G2 for n_wheels==1 or n_wheels==2 in noise_landing_gear with the equation from the original source, I fell like the parenthesis are not correct. Instead of:
G2 = (13+np.log10(2.0(frequencyD/(velocity_fts(1-Mnp.math.cos(theta)))2.0)) \
(30+(frequencyD/(velocity_fts(1-Mnp.cos(theta))))8)-1(0.34H/D))* \
(np.math.sin(phi))2
I think it should read:
G2 = 13+np.log10(2.0(frequencyD/(velocity_fts(1-Mnp.math.cos(theta)))2.0) \
(30+(frequencyD/(velocity_fts(1-Mnp.cos(theta))))8)-1(0.34H/D))* \
(np.math.sin(phi))2
Similarly, for G1 for n_wheels = 4, instead of:
G1 = 12+np.log10(frequencyD/(velocity_fts(1-M*np.cos(theta))))2 \
(0.4+(frequencyD/(velocity_fts(1-Mnp.cos(theta))))2)*(-1.6)
I think it should read:
G1 = 12+np.log10(frequencyD/(velocity_fts(1-Mnp.cos(theta)))2 \
(0.4+(frequencyD/(velocity_fts(1-Mnp.cos(theta))))2)**(-1.6))
(although this is a bit trickier, because it looks like the log only applies to the Strouhal number squared. However, comparing with the other equations I think like this it'd be correct.) Would be nice to discuss this.
Furthermore, according to p. 28 in Fink, the case distinction for n_wheels as implemented in noise_landing_gear seems to be incorrect. Fink states:
Hence, there should be a distinction between n_wheels = 1 and n_wheels = 2, or am I wrong?
Would be nice to add the actual source of the equations/ method to the header: Martin R. Fink. Airframe Noise Prediction Method. Final Report, March 1977, US Department of Transportation.
G2 = (13+np.log10(2.0(frequencyD/(velocity_fts(1-Mnp.math.cos(theta)))2.0)) \ (30+(frequencyD/(velocity_fts(1-Mnp.cos(theta))))8)-1(0.34H/D))* \ (np.math.sin(phi))2 I think it should read: G2 = 13+np.log10(2.0(frequencyD/(velocity_fts(1-Mnp.math.cos(theta)))2.0) \ (30+(frequencyD/(velocity_fts(1-Mnp.cos(theta))))8)-1(0.34H/D))* \ (np.math.sin(phi))2
Similarly, for G1 for n_wheels = 4, instead of: G1 = 12+np.log10(frequencyD/(velocity_fts(1-M*np.cos(theta))))2 \ (0.4+(frequencyD/(velocity_fts(1-Mnp.cos(theta))))2)*(-1.6) I think it should read: G1 = 12+np.log10(frequencyD/(velocity_fts(1-Mnp.cos(theta)))2 \ (0.4+(frequencyD/(velocity_fts(1-Mnp.cos(theta))))2)**(-1.6)) (although this is a bit trickier, because it looks like the log only applies to the Strouhal number squared. However, comparing with the other equations I think like this it'd be correct.) Would be nice to discuss this.
Furthermore, according to p. 28 in Fink, the case distinction for n_wheels as implemented in noise_landing_gear seems to be incorrect. Fink states: Hence, there should be a distinction between n_wheels = 1 and n_wheels = 2, or am I wrong?
Would be nice to add the actual source of the equations/ method to the header: Martin R. Fink. Airframe Noise Prediction Method. Final Report, March 1977, US Department of Transportation.