reaktionsdiffusionsgleichung - uni ulm aktuellesΒ Β· 𝑛 =π·π‘›βˆ‡2π‘›βˆ’βˆ‡β‹… π‘˜1𝑛...

18

Upload: others

Post on 10-Oct-2019

1 views

Category:

Documents


0 download

TRANSCRIPT

Reaktionsdiffusionsgleichung

𝑒 = 𝑓 𝑒 + βˆ‡ β‹… π·βˆ‡π‘’

𝑒 = 𝑓 𝑒 βˆ’ βˆ‡ β‹… 𝑱

KontinuitΓ€tsgleichung

𝑱 = βˆ’π·βˆ‡π‘’

Fourierβ€˜sches Gesetz, Erstes Fickβ€˜sches Gesetz

β€žVon nichts kommt nichtsβ€œ

u: Dichte einer Grâße

f: Reaktionsterm (z. B. chemische Reaktion, Mitose, …)

falls 𝑓(𝑒) ≑ 0: u ist eine lokal erhaltene Grâße

(Masse, Impuls, Energie, Ladung, …)

J: Stromdichte von u

Divergenz von J negativ β†’ Senke β†’ Zufluß aus Umgebung

Divergenz von J positiv β†’ Quelle β†’ Abfluß in Umgebung

β€žFluß(-dichte) ist proportional zum

WΓ€rme-/Konzentrationsgradienten.β€œ

D: Diffusionstensor(-/koeffizient)

Reaktion Diffusion

Reaktionsdiffusionsgleichungen

𝑒 = 𝑒′′ + 𝑒 1 βˆ’ 𝑒

Fisher-Kolmogorov-Gleichung (1D)

Anwendung: Populationsdynamik

http://www.uni-muenster.de/Physik.AP/Purwins/RD/struktur-e.gif

𝑒 = 𝑓 𝑒, 𝑣 + π·π‘’βˆ‡2𝑒

𝑣 = 𝑔 𝑒, 𝑣 + π·π‘£βˆ‡2𝑣

Aktivator-Inhibitor-Systeme

𝒖 = 𝒇 𝒖 + βˆ‡ β‹… π·βˆ‡π’–

RD-System

https://homepages.warwick.ac.uk/

staff/D.Barkley/Research/

spiral_spectra/node1.html

Mit βˆ‡2𝑒 = Δ𝑒 = 𝛻 β‹… 𝛻𝑒 = π‘–πœ•2𝑒

πœ•π‘₯𝑖2 (Laplace-Operator)

𝑒 = 𝑒 1 βˆ’ 𝑒2 βˆ’ 𝑣 + πœ… + π·π‘’βˆ‡2𝑒

Z. B. Gierer-Meinhardt-Modell

𝜏 𝑣 = 𝑣 βˆ’ 𝑒 + π·π‘£βˆ‡2𝑣

Bewegung von Bakterien, Zellen etc.

Stylonychia

https://commons.wikimedia.org/wiki/File:Stylonychia_Side_view.jpg

ZufΓ€llig, β€žRandom Walkβ€œ: Diffusion Zielgerichtet: Taxis-Term

http://en.wikipedia.org/wiki/File:Random_walk_25000_not_animated.svg

Taxis

http://en.wikipedia.org/wiki/File:Chtxphenomen1.png

Außerdem

- Galvanotaxis

- Aerotaxis

- Anemotaxis

- Barotaxis

- Durotaxis

- Gravitaxis

- Hydrotaxis

- Magnetotaxis

- Phototaxis

- Rheotaxis

- Thermotaxis

etc.

(Chemo-)Taxis-Modell

𝑱chemo = πœ’(𝑒1, 𝑒2)βˆ‡π‘’2

𝑱diff = βˆ’π·βˆ‡π‘’1

𝑱 = 𝑱diff + 𝑱chemo 𝑒1 = 𝑓 𝑒1, 𝑒2 + βˆ‡ β‹… π·βˆ‡π‘’1 βˆ’ βˆ‡ β‹… πœ’ 𝑒1, 𝑒2 βˆ‡π‘’2

Reaktion Diffusion Chemotaxis

In KontinuitΓ€tsgleichung:

𝑒2 = 𝑔 𝑒1, 𝑒2 + βˆ‡ β‹… π·βˆ‡π‘’2

Reaktion Diffusion

Ziel: 𝑒1 (Zellen, Bakterien) soll sich entlang des Gradienten von 𝑒2 bewegen

mit

und

ChemosensitvitΓ€t

Keller-Segel-Modell (Keller & Segel 1970)

(Problem: Blow-Ups/SingularitΓ€ten:

Partikeldichte kann beliebig groß werden)

𝑱chemo = πœ’ 𝑒1, 𝑒2 βˆ‡π‘’2

πœ’ 𝑒1, 𝑒2 =π‘˜1𝑒1

π‘˜2 + 𝑒22

Lapidus & Schiller 1976

Woodward et al. 1995

(empirisch, β€žphΓ€nomenologischβ€œ)

πœ’ 𝑒1, 𝑒2 = π‘˜1𝑒1

Escherichia-coli-Modell

aus Mathematical Biology II, S. 260ff

Stewart et al. 2005

𝑛 = π·π‘›βˆ‡2𝑛 βˆ’ βˆ‡ β‹…

π‘˜1𝑛

π‘˜2 + 𝑐 2βˆ‡π‘ + π‘˜3𝑛 π‘˜4

𝑠2

π‘˜9 + 𝑠2βˆ’ 𝑛

Diffusion Chemotaxis Proliferation, Nekrose

𝑐 = π·π‘βˆ‡2𝑐 + π‘˜5𝑠

𝑛2

π‘˜6 + 𝑛2βˆ’ π‘˜7𝑛𝑐

Diffusion Produktion Verbrauch/Aufnahme

𝑠 = π·π‘ βˆ‡2𝑠 βˆ’ π‘˜8𝑛

𝑠2

π‘˜9 + 𝑠2

Diffusion Verbrauch

Budrene & Berg 1995

Death rate prop. to n

Birth rate ~ Logistic growth

Nahrungsaufnahme ~ Proliferation

Zelldichte

Chemo-

attraktor

Stimulans

Mathematical Biology II, S. 260ff

Vereinfachtes Modell

𝑏 = π·π‘βˆ‡2𝑏 βˆ’ π‘˜1βˆ‡ β‹…

𝑏

π‘˜6 + 𝑛 2βˆ‡π‘› + π‘˜2𝑏

π‘˜3𝑛2

π‘˜4 + 𝑛2βˆ’ 𝑏

Proliferation, NekroseDiffusion Chemotaxis

𝑛 = π·π‘›βˆ‡2𝑛 βˆ’ π‘˜5𝑏

𝑛2

π‘˜4 + 𝑛2

VerbrauchDiffusion

𝑏: Bakterienkonzentration (relativ, 𝑏0 ∈ 0, 1 mβˆ’2)

𝑛: NΓ€hrstoffkonzentration (relativ, 𝑛0 ∈ 0, 1 mβˆ’2)

Equation-BasedModeling