問題タブ [fipy]

この質問は、おそらくPythonの知識が不足しているためですが、ここでは.

http://www.ctcms.nist.gov/fipy/examples/diffusion/generated/examples.diffusion.circle.html#module-examples.diffusion.circle

このような手順: ファイルを開く; data_var のコンテンツを読み取り/保存します。mesh = Gmsh2D(data_var) # gmsh でメッシュを生成

ありがとう

最近、FiPy パッケージをすべての依存関係と共に MacPorts を介して Macbook にインストールしました。Python で FiPy と NumPy をパッケージとして呼び出すのに問題はありません。

動作するようになったので、例を見ていきたいと思います。しかし、自分のコンピューターに「ベース ディレクトリ」または FiPy ディレクトリが見つかりません。

ベースディレクトリを見つけるにはどうすればよいですか? これらすべてを Macports 経由でインストールした場合、ベースディレクトリさえありますか?

注意として、私はPython27を使用しています。

助けてください！ありがとう。

Mac (OS X 10.10.2 および Python 2.7.6) に FiPy を新規インストールした後、ドキュメントの例をいくつか試しましたが、うまくいきませんでした。そこで、組み込みのテスト ルーチンを使用して FiPy をテストしました (Python シェルから "fipy.test()" と入力) 。 com/alfrenardi/3b1ad5ff814077b18882

LTS Ubuntuの新規インストールも試しましたが、同じ問題が発生しました。回避策はありますか？私の意見では、これはプラットフォームの問題ではなく、FiPy モジュールの問題のようです。

1D 拡散問題である fipy の例のリストから例を実装しようとしています。しかし、結果をプロットとして表示することはできません。例のコードで提案されているように、ビューアーを正しく定義しました。まだ役に立ちません。解ベクトルは正常に実行されます。しかし、ビューア機能を使用してプロットすることはできません。誰でも助けることができますか？ありがとう！

レベル セット 1D の問題に対して以下のコードを実行しようとしています (fipy Web ページの例)。このエラーが発生しています：

トレースバック (最後の最後の呼び出し): ファイル "C:/Users/sgowda/Documents/pde solver code/level set 1D.py"、20 行目、var.calcDistanceFunction() ファイル "C:\Users\sgowda\AppData\ Local\Continuum\Anaconda\lib\site-packages\fipy\variables\distanceVariable.py", line 335, in calcDistanceFunction raise Exception, "Neither not found on the $PATH" 例外: どちらlsmlibも$に見つからない道skfmmlsmlibskfmm

これを修正する方法を教えてください。distancefucntion() を調べてみましたが、このエラーの意味がわかりませんか?

これは簡単な質問かもしれませんが、勾配の外側にある空間的に変化する係数で PDE を解きたい場合、使用する正しい FiPy 構文は何ですか? これまでに見たすべての例は、勾配内の係数についてのみ説明しています。

例えば：

d/dt(シグマ) = (1/r) d/dr (r^0.5 d/dr(nu シグマ r^0.5))

（数値要因は無視しています）

Sigma(t,r) について解きたいと思います。d/dr の前で (1/r) を処理するにはどうすればよいですか?

この単純な方程式はマッサージできるので、勾配の外側にある空間的に変化する係数について心配する必要はありません (または単に係数を時間導関数項内に移動するだけです)。私が解決しようとしている実際の問題であり、そのトリックはもはや有効ではありません。たとえば、方程式が次のようになっている場合はどうすればよいですか。

d/dt (変数) = f(r) d^2/dr^2 (変数) + g(r) d/dr (変数)

どんな助けでも大歓迎です！

3D メッシュに変数があり、計画を切り出そうとしています。この質問が以前に尋ねられたことがないことに驚いています。簡単で一般的な問題に見えますが、良い方法が見つかりませんでした。アドバイスをいただければ幸いです。

3x3x5 の平行六面体があり、z 平面を抽出しようとしているとしましょう。

出力：

3D 変数は正しく塗りつぶされているように見えます。

まず、このリンクhttp://permalink.gmane.org/gmane.comp.python.fipy/1576で説明されている方法を使用しようとしました

CellVariableのcallメソッドは、 callメソッドを介して渡された一連の座標への補間を可能にします( callメソッドは、関数呼び出しのように括弧を使用してアクセスされます)。渡された各座標に対応する一連の値を返します。順序引数は、補間の順序を決定するだけです。

これが実際にどのように機能するかはわかりませんが、私が理解していることから、これは単一の平面を 0 次で補間する必要があるため、特定の平面で正確な値を抽出する必要があります。間違っている場合は修正してください。

奇妙なことに、それは機能しません。偶数のインデックスは問題ありませんが、奇数のインデックスは隣接するプレーンのコピーです。

どこかに愚かな間違いがあるに違いないと思いますが、手がかりがありません。何か案が？

Python で FiPy を使用していくつかの微分方程式を解こうとしていますが、初心者として、まだいくつかの問題があります。私がしていることは次のとおりです。セル変数を定義し、この変数の方程式を解き、それを更新します。繰り返しのたびにその値を保存したいと思います。次に例を示します。

私は自分の間違いに気づきました.「a」を更新するたびに「a_tt」の値が更新されるため、最後にすべて同じ要素を持つ配列があります。これを回避するには、代わりに何をすべきですか？

2 つ (またはそれ以上) のメッシュを連結することは可能ですか? 例：

同じことが 1D または 2D グリッドにも当てはまります。

中間の質問: オフセットを定義することは可能ですか?

私は、定常状態近似のための 2 次元熱伝導拡散モデルの開発にかなりの時間を費やしてきました。

簡単にするために、層流せん断膜、つまり底部での速度がゼロで、速度の増加が一定であると考えます。

熱容量は、温度に対して一定または直線的に増加する可能性があります。

境界条件は、一定の入口温度 (左) と一定の入力流束 (上) であり、すべての外面に勾配がないように強制されています。

こちらのコードを参照してください。

一定の熱容量を使用する場合、入力電力は出力電力に等しくなります。

非一定の熱容量を使用する場合、それらは大幅に異なります。熱容量が温度によって変化するほど、入力と出力の差が大きくなります。

可変速度係数 (ここでは v * c * rho) の導入は、fipy FAQ (拡散項の例のみを明示的に示しています) に記載されているとおりに行われました。グリッド解像度は出力電力を変更しません。したがって、グリッドの問題ではないと思います。また、過渡項を追加して、非常に高い時間ステップで解こうとしましたが、解は変わりません。

対流項を定義するときに何かひどく間違ったことをしたのではないかと心配していますが、エラーを見つけることができません。また、fipy がtheta(ランク = 0 のセル変数) をvelocity(ランク = 1 のセル変数) と混合し、それらを対流項に必要な面変数にキャストできるかどうかについても混乱しています。