fipy

0熱

1答えて

私は、2D Numpy配列の値を2D FiPyグリッドに転置して、システムを解決する必要があるのです。次のように私のコードは次のとおり tmpSourceとtmpSink前記numpyの配列で for x in range(0, size[0]): for y in range(0, size[1]): source.setValue(tmpSource.getGr

2熱

1答えて

Fipyで複数のPDEを解く

私はFipyを使って、Pythonで偏微分方程式のシステムを解こうとしています。私は少なくとも3人のPDEのシステムを持っています。このようなシステムを解決するにはどのアプローチが最善であるのだろうか？ Fipyは、PDEシステムのためにcoupled and uncoupledアプローチをサポートしています。私はFipyが2連以上の連立方程式をサポートしていても、もしそうでなければ（そうであれ

0熱

1答えて

Fipyの効率を上げる方法

私はfipyを使って、計算領域の分離された領域の拡散問題を解決しています。回路図を参照してください。 @Daniel Wheelerの助けを借りて、fipy.FaceVariableによって定義されたvaring coeffを持つFipyでモデル化されます。しかし、計算速度は私のコマンドを満たすことができませんでした。これは、有限差分法を使ったサイモンコードよりもはるかに遅いです。 fipy計算

0熱

1答えて

一般的な境界条件

-grad(y) + g(y) = 0ここで、gは不明な関数の一部ですyです。ここで私は仕事を得ることができないという単純な1次元の例だ： N=3 h=1./(float(N)-1.) mesh = Grid1D(nx=N, dx=h) c=CellVariable(mesh=mesh,value=0.5) ## Dirichlet boundary conditions #c.c

2熱

1答えて

PythonソースでPDEを解決する

私はBiologyに触発された問題に対処するためにFiPyを使用しています。本質的に、異なる点にソースとシンクがある2D平面を表現したいと思います。ソースは一定のレート（異なるソースは異なるレートを有することができる）で基板を放出し、シンクは固定レートで基板を消費する（異なるシンクは異なるレートを有することができる）。私のコード： import numpy.matlib from fipy i

2熱

1答えて

Fipyでの移流係数の操作：Advection Diffusion Equation

私はfipyを使ってPythonでの移流拡散方程式を解こうとしています。私はそれがドメインの中心を指すように対流係数を操作したいと思います。私のコードは、あなたが見ることができるように、時間に沿って、波が対流係数のベクトルによって確立された方向に移動 from fipy import * # Setting mesh and discretising space nx = 10 dx =