Dear fipy users and developers,
I'm beginning my adventure with computer modelling in FiPy (which is a
great software), and I have some questions considering a modelling of
the superlattice.
I have a paraelectric/ferroelectric (PE/FE) superlattice, consist of
ferroelectric and paraelectric layers. The aim of the simulation is to
calculate the polarization and electrostatic potential distribution
versus temperature.
Because I have different sets of equations for P and \phi in PE and FE
I've decided to make two meshes, one with FE layer and second with PE layer.
The problem I have is to incorporate boundary conditions. I have a
boundary conditions in a form (at the interface between FE and PE
layers) {latex}:
\begin{eqnarray}
\frac{\partial \phi^f}{\partial z} - \epsilon_p \frac{\partial
\phi^p}{\partial z} & = & 4 \pi P, \label{pefe1}\\
\phi^f & = & \phi^p, \label{pefe2}\\
\frac{\partial P}{\partial z} & = & \lambda P \label{pefe3}
\end{eqnarray}
The question I would like to ask is: how to incorporate those conditions
and join these two meshes together.
Any help would be greatly appreciated.
Respectfully yours,
Dariusz Kedzierski
