Spotlight #28 – Akustyka i muzyka: fale stojące, wzory Chladniego, mody pomieszczenia i fizyka dźwięku

Akustyka to najstarsza gałąź fizyki. Pitagoras badał stosunki długości strun. Galileusz upuszczał wahadła, by liczyć częstotliwości. Chladni posypywał piaskiem płyty skrzypiec i odkrywał kształty modów dźwiękowych. Sabine mierzył pogłos w salach wykładowych Harvardu i zapoczątkował akustykę architektoniczną. Dziś dziedzina ta obejmuje projektowanie sal koncertowych, sonar okrętów podwodnych, ultrasonografię medyczną, słuchawki z redukcją szumów i implanty ślimakowe. Wszystko to łączy równanie falowe — proste równanie różniczkowe cząstkowe opisujące propagację i interferencję zaburzeń ciśnienia.

Symulacje w tej fali eksplorują akustykę w sześciu skalach: jednowymiarowe fale stojące w rurze, dwuwymiarową geometrię węzłową na drgającej płycie, trójwymiarowe rezonanse pomieszczenia, czas pogłosu oraz mapowanie częstotliwości na miejsce wewnątrz ślimaka. Każda skala ujawnia coś innego o tym, jak fale zachowują się, napotykając granice.

1. Fale stojące w rurach — harmoniczne i warunki brzegowe

Gdy fala sinusoidalna przemieszcza się wzdłuż rury i odbija się od jej końca, fala biegnąca i odbita nakładają się, tworząc falę stojącą. Wzorzec przestrzenny jest zamrożony w miejscu: węzły (zerowe wychylenie) i strzałki (maksymalne wychylenie) występują w stałych pozycjach. Dopuszczalne częstotliwości są w całości wyznaczone przez długość rury i warunki brzegowe na każdym z końców.

Równanie falowe (1D, ciśnienie):
  ∂²p/∂t² = c² ∂²p/∂x²   gdzie c = prędkość dźwięku (343 m/s w powietrzu przy 20°C)

Warunki brzegowe (dla ciśnienia p i wychylenia u ≡ −∂p/∂x):
  Koniec otwarty:   węzeł ciśnienia      p = 0 w x = 0 lub L
  Koniec zamknięty: węzeł wychylenia     ∂p/∂x = 0 w x = 0 lub L

Częstotliwości harmoniczne:
  Otwarty–otwarty (oba końce otwarte):
    f_n = n · c / (2L),  n = 1, 2, 3, …  (obecne wszystkie harmoniczne)
    Strzałka na każdym końcu; podstawowa = c/2L

  Otwarty–zamknięty (jeden koniec otwarty, jeden zaślepiony):
    f_n = (2n−1) · c / (4L),  n = 1, 2, 3, …  (tylko nieparzyste harmoniczne)
    Podstawowa = c/4L, o oktawę niżej niż otwarty–otwarty tej samej długości
    Wyjaśnia, dlaczego zaślepiona piszczałka organowa brzmi o oktawę niżej

  Zamknięty–zamknięty (oba końce zaślepione):
    f_n = n · c / (2L),  n = 1, 2, 3, …  (te same częstotliwości co otwarty–otwarty,
    ale węzeł na każdym końcu zamiast strzałki)

Fala stojąca wychylenia (otwarty–otwarty, n-ty mod):
  u(x, t) = A · cos(nπx/L) · cos(2πf_n t)
  Węzły w x = L/(2n), 3L/(2n), …
  Strzałki w x = 0, L/n, 2L/n, …

Rzeczywiste instrumenty (korekta końcowa):
  Otwarty koniec nie jest idealnym węzłem ciśnienia; efektywna długość wynosi L + 0,6r
  gdzie r = promień rury (korekta końcowa Rayleigha).
  Uwzględnia częściowe promieniowanie na otwartym końcu.

Różnica między rurami otwarty–otwarty a otwarty–zamknięty wyjaśnia różnicę barwy między fletem a klarnetem. Flet jest otwarty na obu końcach i wytwarza wszystkie harmoniczne; klarnet jest efektywnie zamknięty przy stroiku i otwarty przy dzwonie, więc przedmuchuje się do trzeciej harmonicznej (kwinta i oktawa powyżej podstawowej), a nie do oktawy, a jego najniższy rejestr jest zdominowany przez podstawową i nieparzyste harmoniczne. Brakujące parzyste harmoniczne w widmie pojedynczego stroika nadają klarnetowi charakterystyczną „pustą" barwę.

2. Wzory Chladniego — kształt dźwięku

W 1787 roku Ernst Chladni posypał piaskiem kwadratową metalową płytę, pociągnął smyczkiem po jej krawędzi i odkrył, że piasek przemieszcza się z obszarów drgających i gromadzi wzdłuż linii zerowego ruchu — linii węzłowych. Każda częstotliwość rezonansowa wytwarzała inny, często zachwycająco piękny wzór geometryczny. Te figury, dziś nazywane wzorami Chladniego, są bezpośrednią wizualną mapą modów własnych drgającej płyty.

Równanie podstawowe (cienka płyta sprężysta, teoria Kirchhoffa–Love'a):
  D ∇&sup4;w + ρh ∂²w/∂t² = 0
  w(x,y,t) = ugięcie płyty
  D = Eh³/[12(1−ν²)] = sztywność zginania  (E = moduł Younga, h = grubość, ν = współczynnik Poissona)
  ρ = gęstość, h = grubość

Dla swobodnie drgającej płyty kwadratowej (bok a), rozdzielenie zmiennych w = W(x,y)e^(iωt):
  ∇&sup4;W = λ&sup4;W,   gdzie λ&sup4; = ρhω²/D

Płyta prostokątna swobodnie podparta (SS) — rozwiązanie ścisłe:
  W_mn(x,y) = sin(mπx/a) · sin(nπy/b)
  f_mn = (π/2) · (D/ρh)^(1/2) · [(m/a)² + (n/b)²]
  Linie węzłowe: pionowe w x = ka/m (k=1..m−1); poziome w y = lb/n (l=1..n−1)

Empiryczna reguła Chladniego (dla płyt kwadratowych, krawędzie swobodne):
  f_mn ∝ (m + 2n)²   ← przybliżenie dla swobodnych warunków brzegowych
  Rzeczywiste mody swobodnej płyty wymagają numerycznej analizy MES; wzory węzłowe
  są bardziej złożone niż w przypadku podparcia swobodnego.

Liczba linii węzłowych:
  Mod (m,n): m−1 linii węzłowych równoległych do Y, n−1 równoległych do X
  Lub m−1 + n−1 zakrzywionych linii węzłowych dla symetrii nieprostokątnej.

Degeneracja i mieszanie modów:
  Gdy dwa mody mają tę samą częstotliwość (f_mn = f_nm dla płyty kwadratowej),
  każda kombinacja liniowa jest również ważnym modem — rzeczywisty wzór zależy
  od warunków początkowych/brzegowych. To mieszanie tworzy obrócone i zakrzywione
  wzory węzłowe obserwowane w rzeczywistych eksperymentach Chladniego.

Wzory Chladniego są czymś więcej niż tylko piękne: lutnicy wykorzystują je do testowania jakości akustycznej wierzchniej i spodniej płyty skrzypiec. Wierzchnia płyta skrzypiec stradivariusa ma określone częstotliwości rezonansowe, które twórca kształtuje i stroi, usuwając drewno, aż wzory Chladniego odpowiadają pożądanym wzorcom. Współcześni projektanci głośników wykorzystują analizę metodą elementów skończonych wzorów węzłowych, by eliminować niepożądane rezonanse w membranach i stożkach.

3. Mody pomieszczenia — rezonanse w przestrzeniach prostokątnych

Sala koncertowa, studio nagraniowe czy salon to rezonator akustyczny. Fale dźwiękowe odbijające się od sześciu ścian tworzą trójwymiarowe fale stojące — mody pomieszczenia — przy określonych częstotliwościach wyznaczonych przez wymiary pomieszczenia. Przy tych częstotliwościach niektóre miejsca w pomieszczeniu mają bardzo wysokie ciśnienie (strzałki), a inne niemal zerowe (węzły). Wrażenia słuchowe zmieniają się drastycznie w zależności od tego, gdzie siedzisz.

Częstotliwości modalne (pomieszczenie prostokątne o sztywnych ścianach, L×W×H):
  f_lmn = (c/2) · √[(l/L)² + (m/W)² + (n/H)²]
  l, m, n = nieujemne liczby całkowite (nie wszystkie zerowe)
  c = prędkość dźwięku

Typy modów:
  Osiowe (jeden indeks niezerowy):     fale tylko wzdłuż jednej osi
    np. f_100 = c/(2L)         ← najprostsze, najsilniejsze mody
  Styczne (dwa indeksy):       fale wzdłuż dwóch osi, odbijające się od czterech ścian
  Skośne (trzy indeksy):        fale wzdłuż wszystkich trzech osi, odbijające się od wszystkich sześciu ścian

Gęstość modów (liczba modów poniżej częstotliwości f):
  N(f) ≈ (4πV/3c³)f³ + (πS/4c²)f² + (L_total/8c)f   (wzór Weyla)
  V = objętość pomieszczenia, S = całkowita powierzchnia ścian, L_total = suma długości wszystkich krawędzi
  Mody są rzadkie przy niskich częstotliwościach, gęste przy wysokich.

Częstotliwość Schroedera (przejście z reżimu modalnego do statystycznego):
  f_S ≈ 2000 · √(RT60/V)   Hz
  Poniżej f_S: dominują mody pomieszczenia; powyżej f_S: obowiązuje założenie pola rozproszonego.
  Typowe małe pomieszczenie (V=40 m³, RT60=0,4 s): f_S ≈ 200 Hz
  Sala koncertowa (V=20 000 m³, RT60=2 s): f_S ≈ 20 Hz

Kształty modów ciśnienia (sztywne ściany):
  p_lmn(x,y,z) = A · cos(lπx/L) · cos(mπy/W) · cos(nπz/H)
  Płaszczyzny węzłowe w pozycjach x = L(2k+1)/(2l), itd.

Problematyczne odstępy (mody zbyt oddalone od siebie):
  Duży odstęp między kolejnymi częstotliwościami modów → nierówna odpowiedź basowa.
  Zalecane proporcje pomieszczenia (kryterium Bonello):
    L:W:H ≈ 1,6:1,25:1 (lub 1,9:1,4:1) — minimalizuje pokrywające się częstotliwości modalne.

Mody pomieszczenia są głównym powodem, dla którego małe studia nagraniowe brzmią źle w zakresie basów. Poniżej około 200 Hz mody są szeroko rozstawione; jeśli siedzisz w strzałce ciśnienia dla modu osiowego 80 Hz, bas brzmi dudniąco; w węźle znika całkowicie. Projektanci studiów umieszczają absorbery, dyfuzory i pułapki basowe w narożnikach (gdzie ciśnienie jest najwyższe dla najniższych modów) i dobierają asymetryczne wymiary pomieszczenia, by rozłożyć częstotliwości modalne.

4. Pogłos sali koncertowej — RT60 i geometria pomieszczenia

Czas pogłosu RT60 to czas potrzebny, by dźwięk osłabł o 60 dB po zatrzymaniu źródła. Został zdefiniowany i po raz pierwszy zmierzony przez Wallace'a Clementa Sabine'a w 1900 roku podczas badania akustyki Fogg Art Museum w Cambridge. Jego empiryczny wzór, wyprowadzony ze starannych eksperymentów ze stoperem i piszczałką organową, pozostaje jednym z najbardziej użytecznych wzorów inżynierskich w akustyce architektonicznej.

Wzór Sabine'a (1900):
  RT60 = 0,161 · V / A
  V = objętość pomieszczenia (m³)
  A = całkowita absorpcja (m² sabin) = Σ α_i S_i
  α_i = współczynnik absorpcji powierzchni i (0 = odbicie doskonałe, 1 = absorpcja doskonała)
  Ważny dla niskiej średniej absorpcji (α_avg < 0,2)

Wzór Eyringa (1930, dokładniejszy przy wyższej absorpcji):
  RT60 = −0,161 · V / [S · ln(1 − α_avg)]
  S = całkowita powierzchnia
  Redukuje się do wzoru Sabine'a, gdy α_avg → 0

Średnia droga swobodna (średnia odległość między odbiciami):
  λ_mfp = 4V/S   (ważne dla pola rozproszonego)
  Liczba odbić na sekundę: c / λ_mfp = cS/(4V)
  Tempo zaniku energii: η = −cSα/(4V) → RT60 = 24 ln(10)/η ≈ 0,161 V/(Sα)

Optymalny RT60 według przeznaczenia pomieszczenia:
  Komora bezechowa:      < 0,1 s
  Studio nagraniowe:     0,2 – 0,4 s
  Mowa (sala lekcyjna):  0,4 – 0,8 s
  Opera:                 1,2 – 1,4 s
  Sala symfoniczna:      1,8 – 2,2 s
  Katedra:               4 – 10 s

EDT (Early Decay Time, wczesny czas zaniku):
  Zanik oceniany od 0 do −10 dB (ekstrapolowany do 60 dB).
  Koreluje ściślej z subiektywnym odczuciem „żywości" niż pełne RT60.
  Duży EDT względem RT60 sugeruje dominację późno docierających odbić.

Pierwszym zastosowaniem wzoru Sabine'a była renowacja Sanders Theatre w 1900 roku: dodając wyściełane siedzenia i draperie, podniósł absorpcję i obniżył RT60 z niezrozumiałych 5,5 sekundy do około 1,2 sekundy. Boston Symphony Hall, zaprojektowana z użyciem wzorów Sabine'a w 1900 roku, pozostaje jedną z akustycznie najlepszych sal koncertowych na świecie — bezpośrednią spuścizną zastosowania fizyki w architekturze.

5. Mechanika ślimaka — ucho wewnętrzne jako analizator częstotliwości

Ucho wewnętrzne wykonuje ciągłą transformatę Fouriera. Ślimak — wypełniona płynem, zwinięta rurka o długości około 35 mm — działa jak mechaniczny analizator częstotliwości: każda pozycja wzdłuż błony podstawnej rezonuje przy określonej częstotliwości, z wysokimi częstotliwościami blisko podstawy i niskimi blisko szczytu. Komórki włoskowate w każdej pozycji przekładają ruch mechaniczny na sygnały nerwowe, tworząc mapę tonotopową.

Tonotopowa mapa częstotliwości (funkcja Greenwooda, 1990):
  f(x) = A · (10^(a·x) − k)
  x = ułamkowa odległość od szczytu (0 = szczyt, 1 = podstawa)
  Ludzki ślimak: A = 165,4, a = 2,1, k = 0,88 (daje f w Hz)
  Szczyt:   ~20 Hz (niska częstotliwość)
  Podstawa: ~20 000 Hz (wysoka częstotliwość)
  ~3,5 mm na oktawę wzdłuż błony podstawnej

Fala biegnąca (model Georga von Békésy'ego):
  Gradient sztywności wzdłuż błony podstawnej (sztywna przy podstawie, wiotka przy szczycie)
  powoduje, że nadchodzące fale ciśnienia biegną od podstawy do szczytu.
  Każda fala narasta do maksymalnej amplitudy w miejscu charakterystycznym dla swojej częstotliwości,
  po czym szybko zanika — tworząc ostry pik kodujący częstotliwość.

Krzywa strojenia mechanicznego:
  Szerokość pasma rezonansu błony podstawnej: Q_mech ~ 3–10 (pasywna)
  Wzmocnienie aktywne (zewnętrzne komórki włoskowate, ruchliwość somatyczna OHC):
    OHC wzmacniają ruch, wyostrzając strojenie do Q ~ 30–100
    Ten „wzmacniacz ślimakowy" wyjaśnia 40 dB wzmocnienia i precyzyjną
    rozdzielczość częstotliwościową zdrowego ucha.

Supresja dwutonowa:
  Ton o częstotliwości f_1 może tłumić odpowiedź na f_2, jeśli są wystarczająco blisko siebie.
  Odzwierciedla kompresyjną nieliniowość wzmacniacza OHC.
  To źródło efektów maskowania wykorzystywanych w kompresji audio (MP3, AAC).

Pasmo krytyczne (skala Barka):
  Rozdzielczość częstotliwościowa ślimaka; tony w obrębie jednego pasma krytycznego maskują się nawzajem.
  Δf_CB ≈ 25 + 75[1 + 1,4(f/1000)²]^0,69 Hz (wzór Zwickera)
  Kodeki audio wykorzystują to zjawisko: odrzucają maskowane składowe poniżej progu maskowania.

Georg von Békésy otrzymał w 1961 roku Nagrodę Nobla w dziedzinie fizjologii lub medycyny za wykazanie fali biegnącej w wyciętych ludzkich ślimakach. Wzmacniacz ślimakowy — aktywna rola zewnętrznych komórek włoskowatych — został odkryty w latach 80. i wyjaśnia, jak zdrowe ucho osiąga zarówno skrajną czułość (wykrywanie drgań mniejszych niż średnica atomu wodoru), jak i precyzyjną selektywność częstotliwościową. Implanty ślimakowe omijają uszkodzone komórki włoskowate, stymulując bezpośrednio nerw słuchowy w wielu pozycjach elektrod wzdłuż mapy tonotopowej, wykorzystując ten sam przestrzenny kod częstotliwości, który naturalnie tworzy błona podstawna.

6. Dyfuzory akustyczne — rozpraszanie dźwięku w przestrzeni

Płaska ściana odbija dźwięk jak lustro: cała energia wraca w jednym kierunku zwierciadlanym. Dyfuzor reszt kwadratowych (QRD) — okresowy układ studzienek o różnych głębokościach — rozprasza padający dźwięk równomiernie na wiele kątów. Głębokości są dobierane z sekwencji reszt kwadratowych, konstrukcji teorioliczbowej odkrytej przez Manfreda Schroedera w 1975 roku. Rezultat: powierzchnia, która wydaje się akustycznie chropowata przy wszystkich częstotliwościach w swoim paśmie projektowym, eliminując olśnienie akustyczne i echo trzepoczące.

Sekwencja reszt kwadratowych (okres N, N = liczba pierwsza):
  s_n = n² mod N   (n = 0, 1, …, N−1)
  Przykład N=7: s = [0, 1, 4, 2, 2, 4, 1]
  Każda sekwencja jest symetryczna: s_n = s_{N-n}

Głębokość studzienki (maksymalna głębokość d_max dla najniższej częstotliwości projektowej f_min):
  d_n = s_n · λ_min / (2N)   lub   d_n = s_n · c / (2N · f_min)
  Całkowity zakres głębokości: d_max = (N−1) · c / (2N · f_min)
  Dla N=7, f_min=500 Hz: d_max = 6 · 343/(14·500) ≈ 294 mm

Rozproszone ciśnienie w polu dalekim (przybliżenie Fraunhofera, M okresów):
  p(θ) = Σ_{n=0}^{N−1} r_n · exp(iknd sinθ)   × sin(Mk d_x sinθ/2)/sin(k d_x sinθ/2)
  k = 2πf/c, d_x = szerokość studzienki, θ = kąt rozproszenia
  Jednorodne |p(θ)|² w funkcji θ jest celem projektowym (równa energia we wszystkich listkach dyfrakcyjnych)

Dlaczego to działa (teoria liczb):
  DFT współczynników odbicia QRD r_n = exp(iπs_n/N) ma stałą amplitudę:
  |DFT{r_n}| = stała dla wszystkich częstotliwości → płaskie widmo mocy kątowej.
  To połączenie sum Gaussa (teoria liczb) z dyfuzją akustyczną.

Współczynnik dyfuzji d:
  d = [(NΣE_i² − (ΣE_i)²) / ((N−1)(ΣE_i)²)]^(1/2)   (norma AES)
  d = 1 dla dyfuzji doskonałej, d = 0 dla odbicia zwierciadlanego.
  QRD zazwyczaj osiąga d > 0,8 w swoim paśmie projektowym.

Dyfuzory Schroedera są dziś standardem w studiach nagraniowych, salach koncertowych i obiektach nadawczych na całym świecie. Ponieważ ich konstrukcja jest zakorzeniona w teorii liczb, różne liczby pierwsze N dają różne szerokości okresów i zakresy częstotliwości. Rzeczywiste instalacje łączą kilka okresów obok siebie; listki dyfrakcyjne wynikające z okresowości sumują się konstruktywnie, dając silne rozproszone rozpraszanie. Związek między czystą matematyką (resztami kwadratowymi, sumami Gaussa) a akustyką pomieszczeń jest jednym z najbardziej eleganckich przykładów abstrakcyjnej matematyki znajdującej bezpośrednie zastosowanie inżynierskie.

Szerszy świat akustyki

Te sześć symulacji obejmuje akustykę od skali otworu rury po salę koncertową, dotykając jednej z najbardziej bezpośrednio postrzegalnych dziedzin fizyki na świecie. Dźwięk jest wszędzie, a rządzące nim równanie falowe powraca w całej fizyce: ta sama matematyka opisuje fale elektromagnetyczne, kwantowo-mechaniczne amplitudy prawdopodobieństwa, fale sejsmiczne w Ziemi i fale grawitacyjne w czasoprzestrzeni. Zrozumienie tego, jak granice, geometria i właściwości materiałowe kształtują fale akustyczne, buduje intuicję, która przenosi się bezpośrednio na wszystkie te dziedziny.