57877615f5
Former-commit-id:2b462d8665[formerly133dc97f67] [formerlya02aeb236c] [formerlya02aeb236c[formerly9f19e3f712]] [formerly2b462d8665[formerly133dc97f67] [formerlya02aeb236c] [formerlya02aeb236c[formerly9f19e3f712]] [formerly06a8b51d6d[formerlya02aeb236c[formerly9f19e3f712] [formerly06a8b51d6d[formerly 64fa9254b946eae7e61bbc3f513b7c3696c4f54f]]]]] Former-commit-id:06a8b51d6dFormer-commit-id:2c3569dd39[formerly9bb8decbcf] [formerly8e80217e59] [formerlye2ecdcfe33[formerly377dcd10b9] [formerly8e80217e59[formerly3360eb6c5f]]] Former-commit-id:e2ecdcfe33[formerly377dcd10b9] Former-commit-id:e2ecdcfe33Former-commit-id:7dbd17a5aa
63 lines
2.7 KiB
Python
Executable file
63 lines
2.7 KiB
Python
Executable file
import matplotlib
|
|
matplotlib.rc('text', usetex = True)
|
|
import pylab
|
|
import numpy as np
|
|
|
|
## interface tracking profiles
|
|
N = 500
|
|
delta = 0.6
|
|
X = -1 + 2. * np.arange(N) / (N - 1)
|
|
pylab.plot(X, (1 - np.tanh(4. * X / delta)) / 2, ## phase field tanh profiles
|
|
X, (X + 1) / 2, ## level set distance function
|
|
X, (1.4 + np.tanh(4. * X / delta)) / 4, ## composition profile
|
|
X, X < 0, 'k--', ## sharp interface
|
|
linewidth = 5)
|
|
|
|
## legend
|
|
pylab.legend((r'phase field', r'level set', r'composition', r'sharp interface'), shadow = True, loc = (0.01, 0.55))
|
|
|
|
ltext = pylab.gca().get_legend().get_texts()
|
|
pylab.setp(ltext[0], fontsize = 20, color = 'b')
|
|
pylab.setp(ltext[1], fontsize = 20, color = 'g')
|
|
pylab.setp(ltext[2], fontsize = 20, color = 'r')
|
|
pylab.setp(ltext[3], fontsize = 20, color = 'k')
|
|
|
|
## the arrow
|
|
height = 0.1
|
|
offset = 0.02
|
|
pylab.plot((-delta / 2., delta / 2), (height, height), 'k', linewidth = 2)
|
|
pylab.plot((-delta / 2, -delta / 2 + offset * 2), (height, height - offset), 'k', linewidth = 2)
|
|
pylab.plot((-delta / 2, -delta / 2 + offset * 2), (height, height + offset), 'k', linewidth = 2)
|
|
pylab.plot((delta / 2, delta / 2 - offset * 2), (height, height - offset), 'k', linewidth = 2)
|
|
pylab.plot((delta / 2, delta / 2 - offset * 2), (height, height + offset), 'k', linewidth = 2)
|
|
pylab.text(-0.06, height - 0.06, r'$\delta$', {'color' : 'k', 'fontsize' : 24})
|
|
|
|
## X-axis label
|
|
pylab.xticks((-1, 0, 1), ('-1', '0', '1'), color = 'k', size = 20)
|
|
|
|
## Left Y-axis labels
|
|
pylab.ylabel(r'\bf{phase field} $\phi$', {'color' : 'b',
|
|
'fontsize' : 20 })
|
|
pylab.yticks((0, 0.5, 1), ('0', '.5', '1'), color = 'k', size = 20)
|
|
|
|
## Right Y-axis labels
|
|
pylab.text(1.05, 0.5, r"\bf{level set} $\phi$", {'color' : 'g', 'fontsize' : 20},
|
|
horizontalalignment = 'left',
|
|
verticalalignment = 'center',
|
|
rotation = 90,
|
|
clip_on = False)
|
|
pylab.text(1.01, -0.02, "-1", {'color' : 'k', 'fontsize' : 20})
|
|
pylab.text(1.01, 0.98, "1", {'color' : 'k', 'fontsize' : 20})
|
|
pylab.text(1.01, 0.48, "0", {'color' : 'k', 'fontsize' : 20})
|
|
|
|
## level set equations
|
|
pylab.text(0.1, 0.85, r'$|\nabla\phi| = 1,$ \newline $ \frac{\partial \phi}{\partial t} + U|\nabla \phi| = 0$', {'color' : 'g', 'fontsize' : 20})
|
|
|
|
## phase field equations
|
|
pylab.text(0.2, 0.15, r'$\mathcal{F} = \int f\left( \phi, c \right) dV,$ \newline $ \frac{ \partial \phi } { \partial t } = -M_{ \phi } \frac{ \delta \mathcal{F} } { \delta \phi }$', {'color' : 'b', 'fontsize' : 20})
|
|
|
|
## these went wrong in pdf in a previous version
|
|
pylab.text(-.9,.42,r'gamma: $\gamma$', {'color': 'r', 'fontsize': 20})
|
|
pylab.text(-.9,.36,r'Omega: $\Omega$', {'color': 'b', 'fontsize': 20})
|
|
|
|
pylab.show()
|