4

示されているように、3 つの偏微分方程式 (PDE) と変数の解析解があります。これらの方程式を使用して、\phi(x,y,t)、p(x,y,t)、C_{a}(x,y,t)、および C_{b}(x,y,t) について解きたいつまり、空間と時間の観点から。

pdepe( )Matlab には、1 次元の放物型楕円偏微分方程式の初期境界値問題を解決する関数があることを知っています。この関数またはMatlabの他の関数を使用して、以下で説明する2次元で結合された問題を解決する方法を知りたいです。

問題:

次の 2 つの方程式は、それぞれ 2 つの種 a と b の偏微分方程式を表します。 ここに画像の説明を入力

D_{h} と q は次のように与えられます。

ここに画像の説明を入力

ここに画像の説明を入力

ここで、R_{a}=R_{b}=R で、R は次のように与えられます。

ここに画像の説明を入力

最後に、最後の方程式は次のように与えられます。

ここに画像の説明を入力

初期および境界条件:

ここに画像の説明を入力

ドメインの合計サイズは 10 cm x 5 cm で、Y 字型のサブドメインの幅は 0.5 cm です。このサブドメインの初期の \phi は 0.50 ですが、周囲のマトリックスでは \phi= 0.26 です。1 Pa と 0 Pa の定数 p は、それぞれ境界 (1) と (2) で維持され、約 10^-3 mm^-1 の勾配に対応します。境界 (3) と (4) の p は、境界 (1) と (2) の間の線形勾配によって決定されます。C_{a} = 2 mol m^-3 および C_{b} = 0.2302 mol m^-3 の定数 C は境界 (3) で維持され、境界 (4) での濃度は C_{a} = に設定されます。 1 mol m^-3 および C_{b} = 0.4603 mol m^-3。境界 (1) での濃度は、境界 (3) と (4) の間の一定の勾配によって決定され、移流フラックス境界条件 $$(\frac{\partial C}{\partial x} = 0)$$ が設定されます。 (2) の出口で。

4

2 に答える 2

1

PDE ツールボックスはありますか?

はいの場合:pdetool進むべき道のようです (私はそれを持っていないので、検証したり実験したりすることはできません。自分で実験する必要があります)。

「いいえ」の場合:これまたはこれを調査する価値があるかもしれません。これらは基本的に、2D 波動方程式の FDM の実装です。それらのカーネルを取得して、連立方程式を解くための手段に変換できます。

おそらく簡単です:ここを見てください。これは、Matlab で使用できるかなりまともな FEM ツールキットです。

于 2012-10-28T16:27:02.830 に答える