Der Knallkrebs und die Physik der Kavitation

Biologie-Linse: Waffe, Signal – und eine Schere, die die Seite wechseln kann

Für den Krebs ist die Knallschere ein Mehrzweckwerkzeug. Die Stoßwelle betäubt oder tötet kleine Beutetiere 1; derselbe Knall dient zur Verteidigung und Verständigung 3. Auch ohne Geräusch trägt die Schere Information: In Experimenten von Hughes reagierten männliche Knallkrebse aggressiv auf einzelne Scheren, die in der geöffneten Drohhaltung fixiert waren – und wie stark, hing von der Größe der Schere im Verhältnis zur eigenen ab 5. Die geöffnete Schere wirkt als optisches Signal der Körpergröße, mit dem Rivalen einander einschätzen, bevor ein Kampf eskaliert.

Die Schere ist außerdem erstaunlich wandelbar. Bei Knallkrebsen der Gattung Alpheus können die große Knallschere und die kleine Greifschere die Rollen tauschen: Geht die Knallschere verloren, wird die Greifschere der Gegenseite im Lauf der folgenden Häutungen zu einer neuen Knallschere umgebaut. Mellon und Stephens zeigten 1978, dass sich Form und Funktion einer Schere sogar dann umwandeln, wenn nur ein Nerv auf der Gegenseite des Körpers durchtrennt wird – das Nervensystem steuert, welche Schere zur Waffe wird 6.

Knallkrebse prägen ganze Lebensräume akustisch. Sie gehören zu den wichtigsten Quellen biologischen Unterwasserlärms in flachen Buchten, Häfen und Küstengewässern. Für die Art Synalpheus parneomeris maßen Au und Banks Spitze-Spitze-Quellpegel von 183 bis 189 dB re 1 µPa in 1 m Abstand, mit Schallenergie von wenigen Kilohertz bis 200 kHz. Dieses Knistern kann Sonaranlagen einschränken und möglicherweise die Laute von Delfinen und Walen überlagern 7. Achtung: Unterwasser-Dezibel beziehen sich auf 1 µPa statt auf 20 µPa wie in Luft – die Werte sind daher nicht direkt mit Schallpegeln an Land vergleichbar.

Physik-Linse: Warum schnelles Wasser siedet

Normalerweise bringen wir Wasser zum Sieden, indem wir es erhitzen. Doch der Siedepunkt hängt genauso vom Druck ab wie von der Temperatur: Auf dem Gipfel des Mount Everest siedet Wasser schon bei etwa 70 °C. Senkt man den Druck weit genug, siedet Wasser sogar bei Zimmertemperatur. Der Knallkrebs senkt den Druck – mit Geschwindigkeit.

Die Bernoulli-Gleichung besagt: Entlang einer Stromlinie bleibt die Summe aus statischem Druck p und dynamischem Druck ½ρv² konstant. Je schneller das Wasser strömt, desto geringer ist sein statischer Druck:

Versluis und Kollegen schätzten die Geschwindigkeit des Strahls aus Hochgeschwindigkeitsaufnahmen der sich schließenden Schere: etwa 25 m/s 1 8. Mit p₀ = 101 kPa (Atmosphärendruck) und ρ = 998 kg/m³ für Wasser beträgt der dynamische Druck bei 25 m/s ½ · 998 · 25² ≈ 312 kPa – das Dreifache des Atmosphärendrucks. Schon lange vorher, bei etwa 14 m/s, erreicht der statische Druck den Dampfdruck von Wasser (2,3 kPa bei 20 °C). Das Wasser kann nicht mehr flüssig bleiben: Es kavitiert. Simulationen zeigen, dass es im Detail komplizierter ist: Der Unterdruck entsteht im Kern eines Wirbelrings, der sich um den Strahl aufrollt – der Hohlraum ist eher ringförmig als kugelförmig 8.

Ingenieure fassen das in der Kavitationszahl σ = (p_∞ − p_v) / (½ρv²) zusammen. Je kleiner σ, desto wahrscheinlicher ist Kavitation. Für einen Strahl mit 25 m/s ist σ ≈ 0,3. Die mechanische Nachbildung der Schere von Tang und Staack erzeugte Strahlen mit etwa 26 bis 37 m/s und σ = 0,14 bis 0,30 – vergleichbar mit dem echten Krebs 3.

Dann folgt der dramatische Teil. Ist der Strahl vorbei, wirkt der volle Umgebungsdruck auf die Blase, und sie fällt immer schneller in sich zusammen. Das bisschen Gas darin wird so heftig komprimiert, dass es sich auf mindestens 5.000 K erhitzt und für einen Augenblick leuchtet – Shrimpolumineszenz, nach Kenntnis der Autoren die erste Beobachtung dieser Art der Lichterzeugung bei einem Tier 2. Im selben Moment rast eine Stoßwelle mit Schallgeschwindigkeit durch das Wasser, rund 1.500 m/s 3.

Dieselbe Physik, die der Krebs als Waffe nutzt, ist für Ingenieure ein Ärgernis: Kavitationsblasen, die an Schiffsschrauben und Pumpenlaufrädern kollabieren, tragen mit der Zeit selbst gehärteten Stahl ab – man spricht von Kavitationserosion.

Mathematik-Linse: Wie lange dauert der Kollaps?

1917 berechnete Lord Rayleigh, wie ein leerer kugelförmiger Hohlraum in einer Flüssigkeit unter einem konstanten Druckunterschied kollabiert 4. Sein Ergebnis für die Kollapszeit ist erstaunlich einfach:

Nehmen wir eine Blase mit einem maximalen Radius von R_max = 3 mm an – eine Blase von wenigen Millimetern Größe. Mit ρ = 998 kg/m³ und p_∞ − p_v = 101.325 Pa − 2.339 Pa ≈ 99.000 Pa ergibt sich:

Aus der Formel folgen direkt zwei Dinge. Die Kollapszeit wächst proportional zur Blasengröße: Eine doppelt so große Blase braucht doppelt so lange. Und sie ist umgekehrt proportional zur Wurzel des Druckunterschieds: In tieferem Wasser, wo der Umgebungsdruck höher ist, kollabieren Blasen schneller.

Warum wird der Kollaps so heiß? Wird das Gas im Inneren schnell komprimiert, hat es keine Zeit, Wärme abzugeben – die Kompression ist adiabatisch. Für ein ideales Gas steigt die Temperatur dann wie

Dabei ist γ der Adiabatenexponent, das Verhältnis der Wärmekapazitäten c_p/c_V des Gases (1,4 für Luft). Eine Blase, die auf ein Zehntel ihres Radius schrumpft, erreicht bereits 293 K · 10^1,2 ≈ 4.600 K – dieselbe Größenordnung wie die mindestens 5.000 K, die aus dem Lichtblitz abgeleitet wurden 2. Die über die ganze Blase verteilte Energie wird in ein 10³ = 1.000-mal kleineres Volumen gebündelt. Echte Blasen enthalten Wasserdampf und verlieren etwas Wärme – die Abschätzung ist daher eine grobe Orientierung, keine genaue Vorhersage.

Informatik-Linse: Eine Blase simulieren

Rayleighs Formel setzt eine leere Blase voraus. Um eine echte Blase zu modellieren – mit Gas im Inneren, Oberflächenspannung und Viskosität –, nutzen Physiker die Rayleigh-Plesset-Gleichung, eine gewöhnliche Differentialgleichung für den Blasenradius R(t). Versluis und Kollegen reproduzierten mit einem Modell dieses Typs sowohl den zeitlichen Verlauf des Blasenradius als auch den Schall des Krebses 1:

Dabei ist p_B der Druck in der Flüssigkeit an der Blasenwand: Gasdruck plus Dampfdruck, abzüglich der Effekte von Oberflächenspannung und Viskosität.

Eine Größe in diesem Modell ist unbekannt: wie viel Gas die Blase des Krebses enthält. Gemessen hat das niemand. Das Programm unten macht deshalb, was Forschende mit einer unbekannten Eingangsgröße tun – es probiert mehrere Werte aus. Es simuliert den Kollaps mit 10, 100 und 1.000 Pa Gasdruck (beim größten Radius), gibt Kollapszeit, kleinsten Radius, höchste Wandgeschwindigkeit und die Temperaturabschätzung aus der Physik-Linse aus und zeichnet den Radius über der Zeit:

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp

# Wasser bei 20 °C
RHO, P_INF, P_V = 998.0, 101_325.0, 2_339.0   # Dichte (kg/m³), Umgebungs- und Dampfdruck (Pa)
SIGMA, MU = 0.0728, 1.0e-3                    # Oberflächenspannung (N/m), Viskosität (Pa·s)
T0 = 293.0                                    # Temperatur des Wassers (K)

R_MAX = 3.0e-3        # größter Blasenradius (m)
GAMMA = 1.4           # Polytropenexponent des Gases
C_WATER = 1482.0      # Schallgeschwindigkeit in Wasser (m/s)

def collapse(p_gas):
    """Simuliert den Kollaps ab R_MAX; p_gas = Gasdruck in der Blase bei R_MAX (Pa)."""
    def rayleigh_plesset(t, y):
        R, dR = y
        p_wall = p_gas * (R_MAX / R) ** (3 * GAMMA) + P_V - 2 * SIGMA / R - 4 * MU * dR / R
        ddR = ((p_wall - P_INF) / RHO - 1.5 * dR**2) / R
        return [dR, ddR]

    def turnaround(t, y):    # Ereignis: Die Blasenwand bleibt stehen und kehrt um
        return y[1]
    turnaround.terminal, turnaround.direction = True, 1

    return solve_ivp(rayleigh_plesset, (0, 2e-3), [R_MAX, 0.0], method="LSODA",
                     events=turnaround, rtol=1e-10, atol=1e-13, max_step=1e-6)

rayleigh = 0.915 * R_MAX * np.sqrt(RHO / (P_INF - P_V))
print(f"Rayleigh-Formel (leere Blase): {rayleigh * 1e6:.0f} µs\n")
print("Gas (Pa)   Kollaps (µs)   R_min (µm)   max. Geschw. (km/s)   T_max (1000 K)")

plt.rcParams["font.size"] = 13      # übereinander und mit größerer Schrift, damit die Grafik auch auf dem Handy lesbar bleibt
fig, (whole, end) = plt.subplots(2, 1, figsize=(6.4, 8.4), layout="constrained")
for p_gas in [10, 100, 1000]:     # wie viel Gas die Blase des Krebses enthält, hat niemand gemessen
    sol = collapse(p_gas)
    r_min, speed = sol.y[0, -1], np.abs(sol.y[1]).max()
    t_max = T0 * (R_MAX / r_min) ** (3 * (GAMMA - 1))    # adiabatische Erwärmung (Physik-Linse)
    print(f"{p_gas:8}   {sol.t[-1] * 1e6:12.0f}   {r_min * 1e6:10.1f}   {speed / 1e3:19.1f}   {t_max / 1e3:14.1f}")
    t_us, r_um = sol.t * 1e6, sol.y[0] * 1e6
    fast = np.abs(sol.y[1]) > C_WATER         # Wand schneller als der Schall: Das Modell gilt nicht mehr
    i = np.argmax(fast) if fast.any() else len(fast)
    for ax in (whole, end):
        line, = ax.plot(t_us[:i + 1], r_um[:i + 1])
        ax.plot(t_us[i:], r_um[i:], ":", color=line.get_color())
    end.annotate(f"{p_gas} Pa", (t_us[-1], r_um[-1]), xytext=(5, 0), textcoords="offset points",
                 va="center", backgroundcolor="white")

for ax in (whole, end):
    ax.axvline(rayleigh * 1e6, color="grey", linestyle="--", label="Rayleigh-Formel")
    ax.set_xlabel("Zeit (µs)")
whole.plot([], [], ":", color="grey", label="Wand schneller als Schall: Modell ungültig")
whole.set(ylabel="Blasenradius (µm)", title="Ganzer Kollaps: Die Kurven decken sich")
whole.legend(loc="lower left")
end.set(ylabel="Blasenradius (µm, logarithmisch)", xlim=(rayleigh * 1e6 - 4, rayleigh * 1e6 + 5),
        yscale="log", title="Letzte µs: Hier endet das Modell")
end.yaxis.set_major_formatter("{x:g}")
plt.show()

• In Colab öffnen
• .py herunterladen
• So führst du den Code aus

Ausgabe · das gibt der Code aus, bei jedem Build neu ausgeführt

Rayleigh-Formel (leere Blase): 276 µs

Gas (Pa)   Kollaps (µs)   R_min (µm)   max. Geschw. (km/s)   T_max (1000 K)
      10            275          3.0                  89.7           1160.6
     100            276         20.5                   5.1            116.3
    1000            279        137.0                   0.3             11.9

Was das Ergebnis zeigt:

• Die Kollapszeit ist robust. Egal wie viel Gas in der Blase ist: Sie kollabiert nach 275–279 µs, höchstens etwa 1 % neben Rayleighs Formel (1,1 % bei der größten Gasmenge). Einer Vorhersage, die kaum von einer unbekannten Eingangsgröße abhängt, kann man trauen.
• Der Endpunkt ist es nicht. Der kleinste Radius reicht von 137 µm bis hinunter zu 3 µm – ein Faktor von etwa 45. Weil die Temperatur mit (R_max/R_min)^{3(γ − 1)} wächst, schwankt die Abschätzung zwischen etwa 12.000 K und über einer Million Kelvin. Die Temperatur kann das Modell nicht festlegen.
• Das Modell zeigt seine eigenen Grenzen. Mit 10 oder 100 Pa Gas bewegte sich die Blasenwand, sobald die Blase kleiner als etwa 90 µm ist (gepunktet in Abb. 4), mit 5 bis 90 km/s – schneller als die Schallgeschwindigkeit in Wasser (etwa 1,5 km/s). Die Rayleigh-Plesset-Gleichung behandelt Wasser als inkompressibel und vernachlässigt Wärmeverluste; in dieser letzten Phase sind ihre Zahlen daher nicht mehr physikalisch. Echte Blasen werden durch Wasserdampf, Wärmeleitung und die Kompressibilität des Wassers abgebremst; der beim Krebs gemessene Lichtblitz deutet auf mindestens 5.000 K hin 2.

Auch einige Ideen aus der Programmierung lohnen einen zweiten Blick:

• Steifheit. Die Gleichung ist steif: Kurz vor dem Kollaps ändert sich der Radius extrem schnell. Ein adaptiver Löser wie LSODA verkleinert seine Schrittweite dort automatisch.
• Ereigniserkennung. Statt zu raten, wie lange simuliert werden muss, stoppt der Code genau dann, wenn die Blasenwand umkehrt – im Moment der stärksten Kompression.
• Validierung. Mit fast keinem Gas im Inneren muss die Simulation mit Rayleighs analytischem Ergebnis übereinstimmen – und das tut sie, auf weniger als eine Mikrosekunde genau. Numerik mit einem bekannten Grenzfall zu vergleichen, ist eine Gewohnheit, die sich lohnt.

Probier es selbst aus – jede Änderung ist eine Zeile:

• Setze R_MAX = 6.0e-3: Die Kollapszeit verdoppelt sich auf 551–557 µs und jeder Radius verdoppelt sich, Geschwindigkeiten und Temperaturen bleiben aber gleich – es kommt nur auf das Verhältnis R_max/R_min an.
• Setze P_INF = 201_325.0 (10 m Wassertiefe): Die Blase kollabiert schon nach 194–195 µs – und heftiger, alle Geschwindigkeiten und Temperaturen steigen.
• Ergänze 10_000 in der Liste der Gasdrücke: Das Gas federt den Kollaps ab, die Blase stoppt nach 308 µs bei etwa 800 µm und erwärmt sich nur auf etwa 1.400 K.

Rayleigh-Formel (leere Blase): 551 µs

Gas (Pa)   Kollaps (µs)   R_min (µm)   max. Geschw. (km/s)   T_max (1000 K)
      10            551          6.0                  89.7           1160.5
     100            552         41.0                   5.1            116.3
    1000            557        274.0                   0.3             11.9

Rayleigh-Formel (leere Blase): 194 µs

Gas (Pa)   Kollaps (µs)   R_min (µm)   max. Geschw. (km/s)   T_max (1000 K)
      10            194          1.7                 304.3           2332.3
     100            194         11.5                  17.1            233.5
    1000            195         77.4                   1.0             23.6

Rayleigh-Formel (leere Blase): 276 µs

Gas (Pa)   Kollaps (µs)   R_min (µm)   max. Geschw. (km/s)   T_max (1000 K)
      10            275          3.0                  89.7           1160.6
     100            276         20.5                   5.1            116.3
    1000            279        137.0                   0.3             11.9
   10000            308        799.8                   0.0              1.4