• ベストアンサー

化学工学 円管内濃度分布

半径Rの管の内部にz方向に速さu(r,z)の流れが存在する。 半径rおよびr+Δr、長さΔzの円管で囲まれた円環部のΔtにおける物質収支を考えて以下の濃度分布を与える基礎式を求めよ。 ただし、定常状態とする。 ∂c/∂t + ∂(uc)/∂z = D(1/r)∂{r(∂c/∂r)}/∂r + D ∂^2c/∂z^2 ちなみに、z軸は管に平行で、Dは拡散係数です。 これを解きたいのですが、何をやっていいのかさっぱり分かりません。 考慮するのは、流れと拡散と…他に何かあるのでしょうか。 解き方が分かる方、宜しくお願いします。

  • 化学
  • 回答数1
  • ありがとう数1

質問者が選んだベストアンサー

  • ベストアンサー
  • tareo
  • ベストアンサー率58% (7/12)
回答No.1

 仰る通り、流れと拡散による物質収支を考える事が基本です。後フィックの法則ですかね。具体的に言うと、まず、流れが発生していないとしたときの、微小体積Δv、ある時間Δt内の拡散による物質収支を考えます。そして、次にz軸方向の流れによる、微小体積Δv、ある時間Δt内の物質収支を考えます。この場合、物質流はucで表される事に注意してください。  例をあげるとz方向の拡散による物質収支は拡散流束をJ(r,z,t)とすると、 {J(r,z,t,)-j(r,z+Δz,t)}・Δt={J(r,z,t,)ーJ(r,z,t,)+∂J/∂z・Δz}・Δt=∂J/∂z・Δz・Δtの様な感じです。  このときテイラー展開を使います。そして、最終的に上で求めた、物質収支のすべての和がΔt、Δvにおける、濃度の変化量(∂c/∂t・Δt・Δv)に等しいとおいて整理すれば求める式が得られると思います。

meines20
質問者

お礼

お礼が遅くなってしまってすいません。 おかげさまで答えを導くことができました。 この手の問題に慣れていないもので、円筒座標系に変換するのに手間取ってしまいましたが… ご解答ありがとうございました。

関連するQ&A

  • 円管内速度分布

    全くわかりません。 何から始めてよいのかわからないので困っています。 よろしくお願いします。 半径R、長さLの管の内部流れについて、半径rおよびr+Δrの円筒で囲まれた円環部の力のバランスを考えて以下の速度分布を与える基礎式を求めよ。ただし、定常状態とする。 円柱を横にした図があり、左端にz=0、Po、右端にz=L、Plが書いてあり管内のrにおける速度をu、r+Δrにおける速度をu+Δuと書いてあります。 基礎式↓ 0=μ/r×δ/δr×{r(δu/δr)}+(Po-Pl)/L あるいは、 0=μ/r×δ/δr×{r(δu/δr)}-δP/δz ここで、μは流体の粘度である。

  • 蒸発後の濃度分布の計算について

    蒸発後の濃度分布の計算について質問です。 10L程の鉛直にした円筒容器底部に数ml程の液体を入れて自然蒸発をさせた実験を行っています。 そこで、円筒容器内に気化した蒸気の濃度分布をエクセルで計算しています。 時間(t)を行のセルに、位置(x)を列のセルに入力して、その間のセルに一次元非定常拡散方程式(半無限領域での拡散)C=C0(1-erf(x/(Dt)^0.5)を用いて濃度を算出しました。 一次元非定常拡散方程式(半無限領域での拡散)C=C0(1-erf(x/(Dt)^0.5)では、円筒容器底部(液体の侵入口)の濃度は常に飽和濃度となってしまいます。 蒸発後の時間(t)のセルにD・d^2c/dc^2・⊿tを組んでみたのですがうまくいきませんでした。 蒸発後の濃度分布の計算はどのように計算式を組めばよいのでしょうか。 よろしくお願いします。

  • 層流速度分布について

    流体の問題なのですが… 非圧縮性粘性流体が水平円管内を定常流れで層流になって流れるとするとき、 円管内の層流速度分布を求めたい。 参考書などを調べたところ 管内速度Vz(半径r0の円)は Vz=-(1/4μ)(dp/dz)(r0^2-r^2) とありました。(多分、この式だと思います) このことから円管内の速度分布は回転放物面であると 考えられる、とあるのですが、 なぜこの式が出てきて、この式から回転放物面と分かるのでしょうか? また、速度分布を求めよと言われた場合 これは式で表すのでしょうか?それともグラフにするのでしょうか? 流体力学の知識がほとんどないため、全く分かりません…。 うまく伝えられないのですが、どなたか解説お願いします。

  • 伝熱工学 移動現象 熱工学に関する問題です。

    問題集に略解しか掲載されておらず、勉強が滞ってしまいました。以下の問題の詳しい解答が知りたいです。しっかりと理解を深めたいのでよろしくお願いします。 内径 r[mm] 外径 R[mm] 長さ h[m]の金属円管の中に 円管と同じ半径をもつ半導体が埋め込まれている。金属円管の熱伝導度はk[j/(m・s・k)]である。この発熱体の発熱量はQ[w]であり、発熱体の上下端は完全に断熱されており、高さ方向には均一に発熱するとする。 (1)金属円管内で熱移動が定常であるとき、金属円管内での熱流速qを円管中心からの距離rの関数として表しなさい (2)金属円管の内壁の温度がt[k]であるとき金属円管外壁の温度Tを求めなさい。ただし定常熱伝導であるので一次元のフーリエの熱伝導の式が適用できるとする。 以上です。よろしくおねがいします。

  • 層流の速度分布について

        R ーーーーーーーーーーーーーーー    r|→U \  U=Um(1-rの2乗/Rの2乗)     |    \ ーーーーーーー→|ーーーー    0|  Um/     |   / ーーーーーーーーーーーーーーー (円管内の層流の速度分布) 半径Rで真ん中の速度Umのとき、半径rで速度Uの速度はU=Um(1-rの2乗/Rの2乗)になるのですが、この式の求め方がわからないのです 詳しく教えてください

  • 流体力学の問題-同心二重円管

    今流体力学の勉強をしており 同心二重円管の問題につまづいております 教えてください。 問 長さhの同心二重円管の円管間の隙間に   はニュートン流体で満たされている。外側の(円管半径R1)は   角速度ωで回転している。内側の円管(半径R2)を静止させるのに   必要なトルクMはいくらか?粘度η 自分で分かっているのは 速度分布を求めてせん断応力から解くということが分かっています 速度分布は円柱座標のナビエストークスの式のθ方向から d/dr*{1/r*(d/dr*(Vθ*r))}=0 を解けば速度分布が出ると思うのですが この微分方程式の解き方が分からなくつまづいております よろしくお願いします

  • 円環状線電荷分布の問題をガウス則を用いて解きたい

    次のような円環状線電荷分布のつくる電場を求める問題をガウス則を用いて解こうと思った場合には具体的にどのようにすればよいでしょうか? (問) xy平面上の半径aの円環に線密度λで電荷が一様に分布している場合のz軸上の電場を求めよ この手の対称性の良くないモデルにおいて電場を求めるにはガウス則は適しておらず、クーロン則と重ね合わせの理を用いて解くべきであるというのは承知した上での質問です。 また、この問以外の対称性の良くないモデルにガウス則を適用するための汎用的な手法があればご教示願います。

  • 流体力学の問題

    図のように垂直に置かれた半径Rの円管内を定常的に流れ落ちる、密度ρ、粘性定数乳μのニュートン流体について考える。この円管のz軸にそって距離Lだけ離れた点における圧力は各々、P0、PL(P0>PL)であり、この圧力勾配は円管内で常に一定である。流れはは層流で十分発達したハーゲンポアズイユ流れであるものとして次の問いに答えなさい。  なお、解答にあたっては重力加速度はgとし、問題に定義していない物理量を使用する場合はその定義を明らかにした上で使用すること、また、答えのみでなくその導出過程も示すこと。 という問題で、 (1)円柱座標系r、θ、zにおいて微笑流体要素の体積はdr、dθ、dzおよび半径rを用いてどのようにあらわせられるか。 (2)流体要素についてz軸方向の力のつり合いを考える事により、z方向の速度vzに関する微分方程式を導出せよ。 なのですが、 (1)はrdrdθdzだと思うのですが、 (2)についてはどういうアプローチで解いていくのかが分かりません。 教えていただけるとありがたいです。

  • 慣性モーメント

    水平にした円環の慣性モーメント;M(r_1^2+r_2^2)/2 垂直にした円環の慣性モーメント;M((r_1^2+r_2^2)/4+d^2/12)   の二つの式を求め方がさっぱりわかりません。 薄い円板の慣性モーメント;Mr^2/2とすると水平にした円環の慣性モーメントのほうはドーナツ型をしていて、厚さは無視できるから外半径で求めたモーメントとから内半径で求めたモーメントをひけばいいのではないかと思ったんですけどなんかちがうみたいで、I_bはさっぱりわかりません。 I_a;水平にした円環の慣性モーメント I_b;垂直にした円環の慣性モーメント M;全体の質量 r_1;円環の内半径 r_2;円環の外半径 d;円環の厚さ r;薄い円板の半径

  • ノイマン境界条件の球の拡散方程式の解析解について

    こんばんは。ノイマン境界条件の球の拡散方程式の解析解が必要なのですが、 WEBではなかなか見つかりません。 ∂c/∂t=D/r^2*∂/∂r(r^2*∂c/∂r) に対して、初期値c=c0、r=Rのとき∂c/∂r=a(定数)での 非定常のものなのです。 ヒントを含めましてよろしくお願いいたします。