以下のコードで定義された「L 字型」ジオメトリで熱拡散の問題を解決したいと考えています。
from fipy import CellVariable, Gmsh2D, TransientTerm, DiffusionTerm
from fipy import Matplotlib2DViewer, Viewer
from fipy.tools import numerix
cellSize_1 = 0.05
mesh = Gmsh2D('''
cellSize = %(cellSize_1)g;
Point(1) = {0, 0, 0, cellSize};
Point(2) = {0, 5, 0, cellSize};
Point(3) = {1, 5, 0, cellSize};
Point(4) = {1, 1, 0, cellSize};
Point(5) = {4, 1, 0, cellSize};
Point(6) = {4, 0, 0, cellSize};
Line(7) = {1, 2};
Line(8) = {2, 3};
Line(9) = {3, 4};
Line(10) = {4, 5};
Line(11) = {5, 6};
Line(12) = {6, 1};
Physical Surface("Top") = {8};
Physical Surface("Inner") = {9, 10};
Physical Surface("Right") = {11};
Physical Surface("Inner") = {12, 7};
Line Loop(13) = {7, 8, 9, 10, 11, 12};
Plane Surface(14) = {13};
''' % locals()) # doctest: +GMSH
# Using this mesh, we can construct a solution variable
phi = CellVariable(name = "solution variable",
mesh = mesh,
value = 0.) # doctest: +GMSH
# We can now create a Viewer to see the mesh
viewer = None
from builtins import input
if __name__ == '__main__':
try:
viewer = Viewer(vars=phi, datamin=-1, datamax=3)
viewer.plotMesh()
except:
print("Unable to create a viewer for an irregular mesh (try "+
"Matplotlib2DViewer or MayaviViewer)")
ここで、「上」と「右」の境界を、定数値が 1 と 0 のディリクレ境界として定義したいと思います。また、「内側」と「右」の境界を、フラックスが 0 に等しいノイマン境界として定義したいと思います。制約の設定について? 私が使用した:
phi.constrain(1, where = mesh.Front)
phi.constrain(0, where = mesh.Top)
phi.faceGrad.constrain(0 * mesh.faceNormals, where = mesh.Inner)
phi.faceGrad.constrain(0 * mesh.faceNormals, where = mesh.Outter)
しかし、「AttributeError: 'Gmsh2D' オブジェクトに属性 'Front' がありません」というエラーだけが表示されました。どうすればこれを修正できますか? どうもありがとうございました。