examples.phase.missOrientation.mesh1D

In this example a phase equation is solved in 1 dimension with a misorientation present. The phase equation is given by:

\[\tau_{\phi} \frac{\partial \phi}{\partial t} = \alpha^2 \nabla^2 \phi + \phi ( 1 - \phi ) m_1 ( \phi , T) - 2 s \phi | \nabla \theta | - \epsilon^2 \phi | \nabla \theta |^2\]

where

\[m_1(\phi, T) = \phi - \frac{1}{2} - T \phi ( 1 - \phi )\]

The initial conditions are:

\[\begin{split}\phi &= 1 \qquad \text{for $0 \le x \le L$} \\ \theta &= \begin{cases} 1 & \text{for $0 \le x \le L/2$} \\ 0 & \text{for $L/2 < x \le L$} \end{cases} \\ T &= 1 \qquad \text{for $0 \le x \le L$}\end{split}\]

and boundary conditions \(\phi = 1\) for \(x = 0\) and \(x = L\).

Here the phase equation is solved with an explicit technique.

The solution is allowed to evolve for steps = 100 time steps.

>>> from builtins import range
>>> for step in range(steps):
...     phaseEq.solve(phase, dt = timeStepDuration)

The solution is compared with test data. The test data was created with a FORTRAN code written by Ryo Kobayashi for phase field modeling. The following code opens the file mesh1D.gz extracts the data and compares it with the theta variable.

>>> import os
>>> from future.utils import text_to_native_str
>>> testData = numerix.loadtxt(os.path.splitext(__file__)[0] + text_to_native_str('.gz'))
>>> print(phase.allclose(testData))
1

Last updated on Feb 06, 2025. Created using Sphinx 7.1.2.