ISO 9613-2 is the international standard for calculating outdoor sound propagation. It accounts for geometric divergence, atmospheric absorption, ground effects, barrier diffraction, and other attenuation factors. GeoNoise implements the core ISO 9613-2 propagation model (geometric divergence, atmospheric absorption, ground effect, barrier diffraction).

How does barrier diffraction work?

Barrier diffraction uses the Maekawa formula: A_bar = 10·log₁₀(3 + 20N) where N is the Fresnel number. Sound bends over and around barriers, with higher frequencies attenuated more than lower frequencies. GeoNoise calculates both over-top and side diffraction paths.

What is the ground dip phenomenon?

The ground dip is a frequency-dependent reduction in sound level caused by destructive interference between direct and ground-reflected sound waves. It typically occurs around 200-500 Hz depending on geometry. GeoNoise models this using coherent two-ray phasor summation.

What is the Delany-Bazley model?

The Delany-Bazley model calculates ground surface impedance from flow resistivity: Z = 1 + 9.08(f/σ)^(-0.75) - j·11.9(f/σ)^(-0.73). It's used to determine how much sound is reflected vs absorbed by the ground. GeoNoise supports both Delany-Bazley and the Miki extension.

GeoNoise - Environmental Noise Modeling & Acoustic Propagation Software

GeoNoise implements ISO 9613-based outdoor sound propagation with coherent phasor summation for multi-path interference modeling. All calculations are performed per octave band (63 Hz – 16 kHz).

L_p(f) = L_W(f) - A_{div} - A_{atm}(f) - A_{gr}(f) - A_{bar}(f)

geometric divergence (A_div) Spherical spreading: $A = 20\log_{10}(d) + 10\log_{10}(4\pi)$

▼

Point source spherical spreading models how sound energy spreads over an expanding sphere.

A_{div} = 20\log_{10}(d) + 10\log_{10}(4\pi)

where d = 3D distance (m), 10·log₁₀(4π) ≈ 10.99 dB for full spherical radiation.

Derivation: Power flux Φ = W/(4πr²), therefore L_p = L_W - 10·log₁₀(4πr²)

Cylindrical spreading (line sources): Uses 2π instead of 4π, as energy spreads over a cylinder's circumference rather than a sphere's surface area. Formula: A_div = 10·log₁₀(d) + 10·log₁₀(2π), where 10·log₁₀(2π) ≈ 7.98 dB.

atmospheric absorption (A_atm) ISO 9613-1: $A_{atm} = \alpha(f) \cdot d / 1000$

ISO 9613-1

▼

Air absorbs sound energy through viscous, thermal, and molecular relaxation processes.

A_{atm} = \alpha(f) \cdot \frac{d}{1000} \text{ dB}

ISO 9613-1: Calculates α from temperature, humidity, and pressure.

α = α_cl + (b_N·f²)/(f_r,N + f²/f_r,N) + (b_O·f²)/(f_r,O + f²/f_r,O)

The coefficient $\alpha$ combines classical losses with N₂ and O₂ molecular relaxation.

Simple: Uses the lookup table below at standard conditions (20°C, 70% RH).

63 Hz	125 Hz	250 Hz	500 Hz	1 kHz	2 kHz	4 kHz	8 kHz	16 kHz
0.1	0.4	1.0	1.9	3.7	9.7	33	117	392

Values in dB/km

ground effect (A_gr) $A_{gr} = A_s + A_r + A_m$ · Table 3 exponentials or two-ray phasor

ISO 9613-2

▼

Ground effect accounts for sound interaction with the ground surface. GeoNoise offers two models:

ISO 9613-2 Table 3 (Standard Compliant):

Divides the path into source, middle, and receiver regions. Each region uses height-dependent exponential functions per octave band.

  S ──────────────────────────────────────── R
  │←── As ──→│←────── Am ────────→│←── Ar ──→│
  │  source  │      middle        │ receiver │
  └──────────┴────────────────────┴──────────┘

A_{gr}(f) = A_s(f) + A_r(f) + A_m(f)

A_s, A_r per band (use h=h_s or h_r, d_p=min(30h, d)):

Band	Formula
63 Hz	$A = -1.5$ (fixed)
125 Hz	$A = -1.5 + G \cdot a'(h, d_p)$
250 Hz	$A = -1.5 + G \cdot b'(h, d_p)$
500 Hz	$A = -1.5 + G \cdot c'(h, d_p)$
1 kHz	$A = -1.5 + G \cdot d'(h, d_p)$
2–8 kHz	$A = -1.5(1-G)$

where the Table 3 functions are:

a'(h,d_p) = 1.5 + 3.0·e^-0.12(h-5)²·(1-e^-d_p/50) + 5.7·e^-0.09h²·(1-e^{-2.8×10⁻⁶d_p²})

b' = 1.5 + 8.6·e^-0.09h²·(1-e^-d_p/50) · c' = 1.5 + 14.0·e^-0.46h²·(…) · d' = 1.5 + 5.0·e^-0.9h²·(…)

A_m: 63 Hz → $-3q$ · 125–8 kHz → $-3q(1-G_m)$ · where $q = \max(0, 1 - 30(h_s+h_r)/d_p)$

Two-Ray Phasor Model (Physically Accurate):

Models coherent interference between direct and ground-reflected paths:

A_{gr} = -20\log_{10}\left|1 + \Gamma \cdot \frac{r_1}{r_2} \cdot e^{jk(r_2 - r_1)}\right|

Γ via Delany-Bazley impedance: Z_n = 1 + 9.08(f/σ)^-0.75 − j·11.9(f/σ)^-0.73. Captures constructive/destructive interference (ground dip).

Model Comparison:

Aspect	ISO 9613-2	Two-Ray Phasor
Standards	Compliant	Physics-based
Frequency	Exponential h/d functions	Wavelength-based phase
Interference	Empirical (Table 3)	Full wave interference
Ground dip	Not modeled	Naturally captured

ISO 9613-2:2024 Corrections:

The 2024 revision adds K_geo geometric correction (§7.3) to scale A_gr down at short distances where source and receiver regions overlap, and barrier ground cancellation (§7.4) to zero positive A_gr for barrier-blocked paths. See the ISO 9613-2:2024 revisions section below.

sommerfeld correction (probe) Spherical wave: $Q = \Gamma + (1-\Gamma)\,F(\nu)$

Weideman

▼

Corrects the plane-wave Fresnel coefficient for spherical wave spreading near the ground. Transitions from Q → 1 (near-field) to Q → Γ_plane (far-field).

Q = \Gamma_{\text{plane}} + (1 - \Gamma_{\text{plane}})\,F(\nu)

Numerical distance (complex):

\nu = \frac{1+i}{2}\sqrt{k\,r_2}\left(\cos\theta + \frac{1}{Z_n}\right)

Boundary loss factor via Faddeeva function:

F(\nu) = 1 + i\sqrt{\pi}\,\nu\;\text{wofz}(\nu) \qquad \text{wofz}(z) = e^{-z^2}\,\text{erfc}(-iz)

wofz computed via Weideman (1994) N=32 rational approximation — O(1), stable across the entire complex plane. No Taylor series or region switching.

Note: The grid engine uses a simplified asymptotic approximation (valid for |w| ≥ 4). The full Weideman N=32 computation is used in the probe engine only.

Limits: ν → 0 ⇒ F → 1 ⇒ Q → 1 (total reflection) · |ν| → ∞ ⇒ F → 0 ⇒ Q → Γ_plane · Passive surfaces: |Q| ≤ 1 enforced.

barrier diffraction Maekawa: $A_{bar} = 10\log_{10}(3 + 20N)$

Maekawa

▼

Sound diffracts over and around barriers. Attenuation depends on the path difference δ.

Over-Top Diffraction (vertical):

    S                    R
     \                  /
      \   δ = A+B-d   /
       A ──→ ● ←── B
             │
      ═══════╧══════   Barrier

\delta = A + B - d \quad | \quad N = \frac{2\delta f}{c}

Insertion loss:

A_{bar} = 10\log_{10}(3 + 20N) \quad \text{capped at 20-25 dB}

Side Diffraction (horizontal):

Sound can also bend around the ends of finite-length barriers:

      ════════════════  Barrier (top view)
     ●                  ●
   Left               Right
   Edge               Edge
     ↑                  ↑
  S ─┘                  └─ R
(around left)      (around right)

Path difference for side: $\delta = |S \to Edge| + |Edge \to R| - |S \to R|$

The path with minimum δ (over-top or around sides) determines the effective attenuation.

Side Diffraction Toggle:

Mode	Behavior
Off	ISO 9613-2 infinite barrier assumption. Only over-top diffraction.
Auto	Enable side diffraction for barriers < 50m (recommended default).
On	Always compute side diffraction for all barriers.

Low frequencies diffract easily → less attenuation. High frequencies are blocked effectively.

ISO 9613-2:2024 Alternative:

The 2024 revision replaces the Maekawa formula with D_z = −10·log₁₀(1 + 2 + C₂/λ·C₃·δ·K_met), which includes explicit meteorological correction for wind/temperature gradient effects on barrier effectiveness. See the ISO 9613-2:2024 revisions section below.

ISO 9613-2:2024 revisions K_geo correction, barrier cancellation, revised D_z + K_met

ISO 9613-2:2024

▼

The 2024 revision (ISO 9613-2:2024 — "Engineering method") supersedes ISO 9613-2:1996 with three corrections to the ground effect and barrier diffraction models. All three are optional in GeoNoise — select the ISO 9613-2:2024 calculation profile or enable individually.

1. K_geo Geometric Correction (§7.3)

At short source–receiver distances, the three-region ground model (A_s + A_r + A_m) becomes geometrically invalid because source and receiver regions overlap. K_geo reduces A_gr in this regime:

K_{geo} = \frac{d_p^2 + (h_s - h_r)^2}{d_p^2 + (h_s + h_r)^2}

where d_p = horizontal distance [m], h_s = source height [m], h_r = receiver height [m]. Applied as: A_gr,2024 = K_geo × A_gr.

Behavior: K_geo → 1.0 for d_p >> h_s + h_r (1996 and 2024 agree at distance). Significant only below ~20 m. Dimensionless: [m²]/[m²].

d_p	h_s	h_r	K_geo	Effect
100 m	1 m	1.5 m	0.999	<0.1% — negligible
20 m	1 m	1.5 m	0.985	~2% reduction
5 m	1 m	1.5 m	0.808	~19% reduction

2. Barrier Ground-Effect Cancellation (§7.4)

When a barrier blocks the direct path, the ground reflection geometry is disrupted. The 2024 revision cancels positive ground attenuation for blocked paths:

A_{gr} = 0 \quad \text{when } A_{gr} > 0 \text{ and path is barrier-blocked}

This is the conservative choice — it predicts higher noise levels at the receiver. Ground boosts (A_gr < 0 from constructive interference) are preserved.

Algebraically: the standard writes A_bar = D_z − A_gr. In the total sum A_div + A_atm + A_gr + A_bar, the A_gr terms cancel. GeoNoise implements this equivalently by setting A_gr = 0.

3. Revised Barrier D_z + K_met (§7.4)

Replaces the 1996 Maekawa formula with a revised insertion loss that includes explicit meteorological correction:

D_z = -10\log_{10}\!\left(1 + 2 + \frac{C_2}{\lambda}\,C_3\,\delta\,K_{met}\right)

Condition: δ > z_min. If δ ≤ z_min, D_z = 0.

z_{min} = -\frac{2\lambda}{C_2 \cdot C_3}

where C₂ = 20 (thin/single-edge) or 40 (thick/building), C₃ = 1 (single diffraction), λ = c/f [m], δ = path difference [m].

Meteorological correction:

K_{met} = \exp\!\left(-\frac{1}{2000}\sqrt{\frac{(\max(d_{ss},d_{sr})+e)\cdot\min(d_{ss},d_{sr})\cdot d}{2(\delta - z_{min})}}\right)

where d_ss, d_sr = source/receiver to diffraction edge [m], e = edge separation [m] (0 for thin barriers), d = total distance [m]. For lateral (horizontal) diffraction: K_met = 1.

Why revised: The 1996 K_met over-corrected at large distances with low barriers, predicting near-zero barrier attenuation that contradicted field measurements. The 2024 formula uses (δ − z_min) in the denominator and max/min ordering to prevent this.

D_z produces ~2–3 dB less attenuation than 1996 Maekawa for the same geometry, consistent with the revision intent. Caps: 20 dB (thin), 25 dB (thick).

Dimensional Verification:

K_geo: [m²]/[m²] = dimensionless ✓
z_min: [m]/dimensionless = [m] ✓
K_met exponent: [1/m]·√([m³]/[m]) = [1/m]·[m] = dimensionless ✓
D_z inner: dimensionless + [1/m]·[m]·dimensionless = dimensionless ✓

wall reflections Image source method: $|\Gamma| \approx 0.9$

▼

Reflections from building walls computed using virtual image sources:

         Wall
           │
   S ──────●────── R
    \      │      /
     \     │     /
      S'   │   (S' = mirror of S across wall)

|\Gamma| \approx 0.9 \text{ per reflection (10% absorption)}

Reflections are validated for line-of-sight and checked against other buildings.

coherent phasor summation $p_{total} = \sum_i p_i \cdot e^{j\phi_i}$

▼

Multiple paths from the same source are summed with phase (capturing interference):

p_{total} = \sum_i p_i \cdot e^{j\phi_i} \quad \text{where } \phi_i = -k \cdot d_i + \phi_{reflection}

This captures:

Constructive interference (paths in-phase → +6 dB)
Destructive interference (paths out-of-phase → deep nulls)
Ground dip phenomenon at specific frequencies
Comb filtering from wall reflections

Multiple sources are summed energetically (incoherent):

L_{p,total} = 10\log_{10}\left[ \sum_j 10^{L_{p,j}/10} \right]

source templates CNOSSOS-EU

▼

GeoNoise includes a library of 40 pre-defined noise source templates based on standardized emission data from ASHRAE Handbook, BS 5228, and FHWA RCNM.

Categories:

Industrial (15) — compressors, pumps, generators, fans, transformers
HVAC (10) — rooftop units, chillers, AHU, heat pumps, VRF
Construction (10) — excavators, pile drivers, jackhammers, bulldozers
Generic (5) — white/pink noise, speech, tonal alarm, subwoofer

Data sourced from ASHRAE Handbook, BS 5228, and FHWA RCNM.

accuracy & limitations Expected precision and model assumptions

▼

Accuracy pending measurements and validation

Known Limitations:

Flat terrain assumed (no topography effects)
Neutral atmosphere (no temperature/wind refraction)
First-order wall reflections only (2nd/3rd order settings in UI are reserved for future implementation)
2.5D geometry (vertical plane diffraction)
Homogeneous ground (single type per scene)
All sources are omnidirectional (no directivity pattern support yet)
Single diffracting edge per path (multiple barrier diffraction not yet implemented)
Moving source animation caps: 3 entities (CPU), 16 entities (GPU)
Grid engine Sommerfeld correction uses asymptotic approximation (|w| ≥ 4 only)

ISO 9613-2:2024 Support:

K_geo geometric correction (§7.3) — optional
Barrier ground-effect cancellation (§7.4) — optional
Revised D_z + K_met barrier formula (§7.4) — optional
All default OFF for backward compatibility with 1996

Numerical Guards:

Distance floor: d_min = 0.1m to prevent division by zero
Angle limits: cos θ clamped for grazing incidence
Impedance cap: |Z| bounded to avoid runaway coefficients
Level floor: -100 dB (acoustic silence threshold)

Map Integration:

Real-world map backgrounds via Mapbox GL JS with synchronized pan/zoom

references Standards and academic sources

▼

ISO 9613-1:1993 — Acoustics — Attenuation of sound during propagation outdoors — Part 1: Calculation of absorption by atmosphere
ISO 9613-2:1996 — Acoustics — Attenuation of sound during propagation outdoors — Part 2: General method of calculation
Maekawa, Z. (1968). "Noise reduction by screens." Applied Acoustics, 1(3), 157-173
Delany, M.E. & Bazley, E.N. (1970). "Acoustical properties of fibrous absorbent materials." Applied Acoustics, 3(2), 105-116
Pierce, A.D. (1981). "Acoustics: An Introduction to Its Physical Principles and Applications." McGraw-Hill
Chien, C.F. & Soroka, W.W. (1975). "Sound propagation along an impedance plane." JASA, 58(6). Defines Q = Γ + (1−Γ)F(w) and the boundary-loss factor
Kearns, J.A. (1989). NASA report — asymptotic spherical-wave over impedance plane, numerical distance, grazing-incidence limits
Weideman, J.A.C. (1994). "Computation of the Complex Error Function." SIAM J. Numer. Anal., 31(5). N=32 rational approximation for wofz(z)
Poppe, G.P.M. & Wijers, C.M.J. (1990). "Algorithm 680: Evaluation of the Complex Error Function." ACM TOMS. Stable wofz(z) reference implementation
Taraldsen, G. (2005). "A note on reflection of spherical waves." JASA. Chien & Soroka lineage, Sommerfeld F context
Al Azah, M. & Chandler-Wilde, S.N. (2020). "Computation of the Complex Error Function using Modified Trapezoidal Rules." Integral representations + symmetry
Gautschi, W. (1969). NASA report — early wofz(z) computational reference for complex arguments
Miki, Y. (1990). "Acoustical properties of porous materials — Modifications of Delany-Bazley models." J. Acoust. Soc. Jpn., 11(1)
Shaw, E.A.G. (1974). "Transformation of sound pressure level from the free field to the eardrum in the horizontal plane." JASA, 56(6), 1848-1861. Measured head-related transfer functions showing 10-15 dB ILD at 6-8 kHz
Kuhn, G.F. (1977). "Model for the interaural time differences in the azimuthal plane." JASA, 62(1), 157-167. ILD increases linearly 1-4 kHz, plateaus at 10 dB
Blauert, J. (1997). "Spatial Hearing: The Psychophysics of Human Sound Localization." MIT Press. Standard reference for ILD of 8-10 dB at 4-8 kHz for rigid sphere model matching average human head
Woodworth, R.S. (1938). "Experimental Psychology." Holt. ITD formula for rigid sphere: ITD = (a/c)(theta + sin(theta)), basis for analytical HRTF model
Algazi, V.R. et al. (2001). "The CIPIC HRTF Database." IEEE WASPAA. Measured HRTF sets at 25 azimuths x 50 elevations for 45 subjects
ISO 9613-2:2024 — Acoustics — Attenuation of sound during propagation outdoors — Part 2: Engineering method for the prediction of sound pressure levels. Supersedes ISO 9613-2:1996. K_geo geometric correction (§7.3), barrier ground cancellation (§7.4), revised D_z + K_met (§7.4)

CORE LEVELS

Overall level at receiver:

L_p = L_W - A_div - A_atm - A_gr - A_bar

Energetic summation of multiple sources:

L_sum = 10 log₁₀( Σ 10^L_i/10 )

Sound power conversions:

L_W = 10 log₁₀(P / 10^-12)

P = 10^-12 · 10^L_W/10

GEOMETRY

Plan distance:

d_2D = √((x₂-x₁)² + (y₂-y₁)²)

3D distance:

d_3D = √((x₂-x₁)² + (y₂-y₁)² + (z₂-z₁)²)

Direct path used for diffraction reference:

d_direct = √(d_2D² + (h_s - h_r)²)

SPREADING

Point source (spherical):

A_div = 20 log₁₀(d) + 10 log₁₀(4π)

Line source (cylindrical):

A_div = 10 log₁₀(d) + 10 log₁₀(2π)

ATMOSPHERIC ABSORPTION

Path loss:

A_atm = α · d / 1000

Where α is the atmospheric attenuation coefficient in dB/km and d is the distance in m

Speed of sound and wavelength:

c = 331.3 + 0.606 T_C

λ = c / f

ISO 9613-1 absorption coefficient:

α = 8.686 f² [ (1.84e-11 (p/p₀)^-1 (T/T₀)^1/2)

+ (T/T₀)^-2.5 (0.01275 e^-2239.1/T (f_r,O/(f_r,O² + f²))

+ 0.1068 e^-3352/T (f_r,N/(f_r,N² + f²))) ]

C = -6.8346 (T₀₁/T)^1.261 + 4.6151

p_sat = p₀ · 10^C, h = RH · p_sat / p

f_r,O = (p/p₀) [24 + 4.04e4 h ((0.02 + h)/(0.391 + h))]

f_r,N = (p/p₀) (T/T₀)^-1/2 [9 + 280 h e^{-4.17((T/T₀)^-1/3 - 1)}]

GROUND EFFECTS

ISO 9613-2 Tables 3-4 (per-band, primary model):

A_gr = A_s + A_r + A_m

A_s, A_r per band (Table 3 exponential functions):

63 Hz: A = −1.5 (fixed)

125 Hz: A = −1.5 + G·a'(h, d_p)

250 Hz: A = −1.5 + G·b'(h, d_p)

500 Hz: A = −1.5 + G·c'(h, d_p)

1 kHz: A = −1.5 + G·d'(h, d_p)

2–8 kHz: A = −1.5(1 − G)

a' = 1.5 + 3.0·e^{−0.12(h−5)²}·(1−e^−d_p/50) + 5.7·e^−0.09h²·(1−e^{−2.8×10⁻⁶d_p²})

A_m (63 Hz) = −3q

A_m (125–8 kHz) = −3q(1 − G_m)

q = max(0, 1 − 30(h_s+h_r)/d_p)

G = ground factor (0=hard, 0.5=mixed, 1=soft)

Two-ray phasor geometry:

r₁ = √(d² + (h_s - h_r)²), r₂ = √(d² + (h_s + h_r)²)

Normalized impedance — Delany-Bazley (0.01 < f/σ < 1.0):

ζ = 1 + 9.08 (f/σ)^-0.75 - j 11.9 (f/σ)^-0.73

Miki extension (f/σ > 1.0):

ζ = 1 + 5.50 (f/σ)^-0.632 - j 8.43 (f/σ)^-0.632

Reflection coefficient (locally reacting):

Γ = (ζ cosθ - 1) / (ζ cosθ + 1)

Two-ray ratio and attenuation:

ratio = 1 + Γ (r₁/r₂) e^{-jk(r₂ - r₁)}

A_gr = -20 log₁₀(max(|ratio|, 10^-6))

BARRIER DIFFRACTION

Thin barrier (single edge, cap 20 dB):

A = √(d(S,I)² + (h_b - h_s)²)

B = √(d(I,R)² + (h_b - h_r)²)

δ = A + B - d_direct

Thick barrier / building (double edge, cap 25 dB):

A = √(d(S,I_in)² + (h_b - h_s)²)

B = d(I_in, I_out) (roof width)

C = √(d(I_out,R)² + (h_b - h_r)²)

δ = A + B + C - d_direct

Maekawa-style insertion loss:

N = 2δ / λ

Thin: A_bar = 10 log₁₀(3 + 20N) cap 20 dB

Thick: A_bar = 10 log₁₀(3 + 40N) cap 25 dB

Blocked path total (ISO 9613-2 §7.4 additive):

A_total = A_div + A_atm(d_actual) + A_bar + A_gr

ISO 9613-2:2024 revised barrier (§7.4):

D_z = −10 log₁₀(1 + 2 + C₂/λ · C₃ · δ · K_met)

z_min = −2λ / (C₂ · C₃) — below this, D_z = 0

C₂ = 20 (thin), 40 (thick) C₃ = 1 (single diffraction)

K_met meteorological correction:

K_met = exp(−(1/2000) · √((max(d_ss,d_sr)+e) · min(d_ss,d_sr) · d / (2(δ−z_min))))

Lateral diffraction: K_met = 1

ISO 9613-2:2024 ground corrections (§7.3):

K_geo = (d_p² + (h_s−h_r)²) / (d_p² + (h_s+h_r)²)

Barrier cancellation: A_gr = 0 when A_gr > 0 and barrier-blocked

WEIGHTING

Octave-band weighting:

L_A = L_band + A(f)

L_C = L_band + C(f)

MOVING SOURCES

Waypoint-based animation for sources and head probes with real-time noise map recompute.

Animation System:

Waypoint tool (W key): Click-to-place waypoints on canvas
Catmull-Rom spline: Smooth interpolation between waypoints
Dual-rate loop: 60 Hz visual rendering, adaptive physics recompute (20 Hz single entity, 5 Hz multi-entity CPU, 60 Hz GPU)
Transport controls: Play/pause, scrub, speed (m/s), loop, ping-pong, reverse

Capacity:

CPU backend: max 3 concurrent animated entities
GPU backend: max 16 concurrent animated entities

Noise map, receivers, and probes recompute in real-time as entities move along their paths. Incremental recompute for single-source movement; full recompute for multi-source.

BINAURAL HEAD PROBE

Virtual binaural listener with L/R ear separation and head-related transfer function modeling.

Head Geometry:

L/R ears: 17.5 cm separation, rotatable head azimuth
Dual probes: Independent coherent ray-tracing per ear
Placement: L key to place, drag to move, lollipop to rotate

HRTF Model (Analytical):

Per-path processing: Arrival angle computed for each traced path
Head shadow: Frequency-dependent attenuation — low frequencies pass through, high frequencies (>2 kHz) attenuated up to 10× on shadow side
ILD/ITD readout: Interaural Level Difference and Interaural Time Difference displayed in inspector

References: Shaw (1974), Kuhn (1977), Woodworth & Schlosberg (1954)

GeoNoise uses two calculation engines optimized for different use cases. Understanding when each is applied helps interpret results correctly.

GRID ENGINE

Used for: Mic Grids, Receivers, Noise Maps

Path Finding:

Single direct path per source-receiver pair
Checks for barrier occlusion → computes diffraction
Ground effect via A_gr model (ISO or Two-Ray)

Ground Effect Models:

ISO 9613-2: Empirical per-band coefficients (no interference)
Two-Ray: Phasor interference between direct + ground reflection

Summation:

L = 10·log₁₀(Σ 10^Lᵢ/10)

Energy (power) summation across sources.

Trade-off: Fast computation for heat maps. With Two-Ray enabled, captures ground dip. No wall reflections or multi-path beyond ground.

The grid engine supports WebGPU acceleration for real-time noise map computation at 60 Hz during animation. GPU↔CPU fallback is automatic — receivers and panels compute on CPU while grid maps leverage GPU parallel processing.

PROBE ENGINE

Used for: Probe microphones only

Path Finding (Ray Tracing):

Direct path
Ground-reflected path (specular reflection)
Wall reflections (first-order image sources)
Side diffraction around barrier edges
Over-top diffraction

Ground Effect Models:

Physical Impedance: Complex reflection coefficient from surface impedance
ISO 9613-2: Tabulated absorption coefficients (0.01–0.60)

Impedance Models (when Physical Impedance selected):

Delany-Bazley:

$Z_n = 1 + 9.08(f/\sigma)^{-0.75} - j \cdot 11.9(f/\sigma)^{-0.73}$

Delany-Bazley-Miki:

$Z_n = 1 + 5.5(f/\sigma)^{-0.632} - j \cdot 8.43(f/\sigma)^{-0.632}$

Reflection Coefficient:

$\Gamma = \frac{Z_n \cos\theta - 1}{Z_n \cos\theta + 1}$

Summation:

p = Σ pᵢ·e^jφᵢ where φ = -k·d

Phasor (coherent) summation with phase from path length.

Trade-off: Full multi-path interference (ground dip, wall comb filtering, diffraction patterns). Slower, used only for point analysis.

CALCULATION ARCHITECTURE

┌─────────────────────────────────────────────────────────────────┐
│                      GeoNoise Calculation Flow                  │
├─────────────────────────────────────────────────────────────────┤
│                                                                 │
│  ┌──────────────┐              ┌──────────────┐                 │
│  │   Sources    │              │  Receivers   │                 │
│  │ (with Lw/f)  │              │    Grids     │                 │
│  └──────┬───────┘              │    Probes    │                 │
│         │                      └──────┬───────┘                 │
│         ▼                             │                         │
│  ┌──────────────────────────────────────────────────────────────┤
│  │                    PATH FINDING                              │
│  ├──────────────────────────────────────────────────────────────┤
│  │                                                              │
│  │   ┌─────────────────────┐    ┌─────────────────────────────┐ │
│  │   │   GRID ENGINE       │    │      PROBE ENGINE           │ │
│  │   │                     │    │                             │ │
│  │   │  • Single path      │    │  • Ray tracing (multi-path) │ │
│  │   │  • Direct + best    │    │  • Ground reflection rays   │ │
│  │   │    diffraction      │    │  • Wall reflection rays     │ │
│  │   │  • A_gr handles     │    │  • Side diffraction paths   │ │
│  │   │    ground effect    │    │  • Used for: Probes only    │ │
│  │   └─────────┬───────────┘    └─────────────┬───────────────┘ │
│  │             │                              │                 │
│  └─────────────┼──────────────────────────────┼─────────────────┤
│                ▼                              ▼                 │
│  ┌──────────────────────────────────────────────────────────────┤
│  │                 ATTENUATION (ISO 9613-2)                     │
│  ├──────────────────────────────────────────────────────────────┤
│  │   For each path, apply per-band:                             │
│  │   • A_div  (geometric spreading)                             │
│  │   • A_atm  (atmospheric absorption)                          │
│  │   • A_gr   (ground effect — ISO or Two-Ray phasor)           │
│  │   • A_bar  (barrier diffraction if blocked)                  │
│  └──────────────────────────────────────────────────────────────┤
│                              │                                  │
│                              ▼                                  │
│  ┌──────────────────────────────────────────────────────────────┤
│  │                    SUMMATION                                 │
│  ├──────────────────────────────────────────────────────────────┤
│  │                                                              │
│  │   ┌─────────────────────┐    ┌─────────────────────────────┐ │
│  │   │   INCOHERENT        │    │      COHERENT (Phasor)      │ │
│  │   │                     │    │                             │ │
│  │   │  L = 10·log₁₀(Σ10^  │    │  p = Σ pᵢ·e^(jφᵢ)          │ │
│  │   │          (Lᵢ/10))   │    │  L = 20·log₁₀|p|            │ │
│  │   │                     │    │                             │ │
│  │   │  • Power summation  │    │  • Phase-aware interference │ │
│  │   │  • Fast, stable     │    │  • Captures nulls & peaks   │ │
│  │   │  • Used by: Grids   │    │  • Used by: Probes          │ │
│  │   └─────────────────────┘    └─────────────────────────────┘ │
│  │                                                              │
│  └──────────────────────────────────────────────────────────────┘
└─────────────────────────────────────────────────────────────────┘

COMPARISON SUMMARY

Aspect	Grid Engine	Probe Engine
Path count	1 per source	Multiple (ray traced)
Ground effect	ISO or Two-Ray phasor	Physical Impedance or ISO 9613-2
Wall reflections	No	Yes (first-order)
Summation	Incoherent (power)	Coherent (phasor)
Ground interference	With Two-Ray only	Always (if enabled)
Multi-path interference	No (walls, etc.)	Yes (walls, diffraction)
Ground dip	With Two-Ray enabled	Naturally emerges
Performance	Fast (grid-parallel)	Slower (per-probe)
Output	Heatmap, single dB value	Frequency response curve

WHY TWO ENGINES?

Grids/Receivers need to compute thousands of points quickly for heat maps. Using multi-path ray tracing everywhere would be prohibitively slow. The ISO 9613-2 single-path approach is designed for this: it's an engineering standard optimized for outdoor prediction across large areas with acceptable accuracy.

Probes are point measurements where you want detailed frequency response. Ray tracing + coherent summation reveals interference effects invisible to incoherent methods: the ground dip phenomenon, wall reflection comb filtering, and barrier edge diffraction patterns.

Both engines share the same ISO 9613-2 attenuation functions for A_div, A_atm, A_gr, and A_bar. The difference is in how many paths are traced and how they're combined at the receiver.