Allen-Cahn-Gleichung

Die Allen-Cahn-Gleichung (manchmal auch zeit-abhängige Ginzburg-Landau-Gleichung) ist eine semilineare parabolische partielle Differentialgleichung und Reaktionsdiffusionsgleichung. Sie wird unter anderem verwendet, um den Phasenübergang von binären Legierungen zu beschreiben. Des Weiteren wird sie auch zur Modellierung der Kristallzüchtung verwendet.[1]

Die Gleichung ist nach John W. Cahn und seinem Doktoranden Sam Allen benannt.[2]

Allen-Cahn-Gleichung

Sei

  • Ω R n {\displaystyle \Omega \subset \mathbb {R} ^{n}} eine offene Menge,
  • v 0 ( x ) L 2 ( Ω ) {\displaystyle v_{0}(x)\in L^{2}(\Omega )} eine beliebige Initialfunktion,
  • ε > 0 {\displaystyle \varepsilon >0} und T > 0 {\displaystyle T>0} zwei Konstanten.

Gesucht ist eine Lösungsfunktion v ( x , t ) : Ω × [ 0 , T ] R {\displaystyle v(x,t):\Omega \times [0,T]\to \mathbb {R} } , welche die Allen-Cahn-Gleichung[3]

t v Δ x v = 1 ε 2 f ( v ) , Ω × [ 0 , T ] {\displaystyle \partial _{t}v-\Delta _{x}v={\frac {1}{\varepsilon ^{2}}}f(v),\quad \Omega \times [0,T]}

löst, wobei

  • Δ x {\displaystyle \Delta _{x}} der Laplace-Operator nach x {\displaystyle x} ist,
  • f ( v ) = F ( v ) {\displaystyle f(v)=F'(v)} die Ableitung für ein nicht-negatives F C 1 ( R ) {\displaystyle F\in C^{1}(\mathbb {R} )} mit zwei Minima F ( ± 1 ) = 0 {\displaystyle F(\pm 1)=0} ist.

Häufig verwendet man folgende Initialbedingung mit der Neumann-Randbedingung

{ v ( x , 0 ) = v 0 ( x ) , Ω × { 0 } n v = 0 , Ω × [ 0 , T ] {\displaystyle {\begin{cases}v(x,0)=v_{0}(x),&\Omega \times \{0\}\\\partial _{n}v=0,&\partial \Omega \times [0,T]\\\end{cases}}}

wobei n v {\displaystyle \partial _{n}v} die Normalenableitung n v = v n {\displaystyle \partial _{n}v=\nabla v\cdot n} (und n {\displaystyle n} die äußere Normale) ist.

F {\displaystyle F} ist ein Energiepotential, häufig wählt man dafür die Funktion

F ( v ) = ( v 2 1 ) 2 4 , f ( v ) = v 3 v {\displaystyle F(v)={\frac {(v^{2}-1)^{2}}{4}},\qquad f(v)=v^{3}-v}

dann schreibt sich die Allen-Cahn-Gleichung als

t v Δ x v = 1 ε 2 ( v 3 v ) , Ω × [ 0 , T ] {\displaystyle \partial _{t}v-\Delta _{x}v={\frac {1}{\varepsilon ^{2}}}\left(v^{3}-v\right),\quad \Omega \times [0,T]}

Herleitung

Betrachte das Freie-Energie-Funktional

E ε [ v ] = Ω ( 1 2 | v | 2 + 1 ε 2 F ( v ) ) d x , {\displaystyle E_{\varepsilon }[v]=\int _{\Omega }\left({\frac {1}{2}}|\nabla v|^{2}+{\frac {1}{\varepsilon ^{2}}}F(v)\right)dx,}

dann erhält man die Allen-Cahn-Gleichung, wenn man den L 2 {\displaystyle L^{2}} -Gradientenfluss des Funktionals berechnet, das bedeutet man berechnet die Gleichung

t v = δ E ε [ v ] δ v , {\displaystyle \partial _{t}v=-{\frac {\delta E_{\varepsilon }[v]}{\delta v}},}

wobei wir rechts die Funktionalableitung genommen haben.[3]

Literatur

  • Sören Bartels: Numerical Methods for Nonlinear Partial Differential Equations. Hrsg.: Springer International Publishing. Deutschland 2015. 

Einzelnachweise

  1. T. Philippe, H. Henry, M. Plapp: A regularized phase-field model for faceting in a kinetically controlled crystal growth. In: Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences. Band 476, Nr. 2241, September 2020, doi:10.1098/rspa.2020.0227, PMID 33071578, PMC 7544347 (freier Volltext). 
  2. S. M. Allen und J. W. Cahn: A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening. In: Acta. Metal. Band 27, 1979, S. 1084–1095. 
  3. a b Sören Bartels: Numerical Methods for Nonlinear Partial Differential Equations. Hrsg.: Springer International Publishing. Deutschland 2015, S. 153. 
Normdaten (Sachbegriff): GND: 1292845384 (lobid, OGND, AKS)