• 締切済み

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

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

みんなの回答

  • metzner
  • ベストアンサー率60% (69/114)
回答No.1

順次積分するだけです。まず0を積分して、次にその答えをrを掛けて。。。

関連するQ&A

  • 流体力学 二重円筒

    次の問題を自分で解いてみたのですが正解であるか不安です。正解かどうか教えてください。お願いいたします。 問題   図のような同軸の二重円筒に流体を満たし内側円筒のみを回転させ、軸トルクを測定することで 流体の粘度を測定しする装置について考える。内側円筒と外側円筒の半径をそれぞれR_i=20mm,  R_o=30mmとし、回転角速度をΩとする。以下の問いに答えよ。 (1)回転開始から時間が経過し、流れが定常に発達した状況を想定する。代表長さを流体層厚さd  =R_o-R_i,代表速さを内側円筒壁面の速度とするレイノズル数Reがおおよそ40以下の場合、流れ は周方向速度のみをもつ軸対称クエット流れとなり、軸トルクの計測から粘度が求まる。測定可能 な最低粘度を0.01Pa・s とするとき、設定可能な最大の円筒回転角速度Ω_maxを求めなさい。ただ し、流体の密度はρ=10^3kg/m^3とする。 (2)(1)の条件が満たされるとき、流れの周方向速度成分u_θは以下の常微分方程式に従う。これ  を適切な境界条件のもとで解き、u_θの半径方向速度分布を求め図示しなさい。          d^2(u_θ)/dr^2+1/r(d(u_θ)/dr)-u_θ/r^2=0  なおω=u_θ/rなる変数変換を用いてよい。 (3)回転角速度Ω=1rad/sに設定し、粘度μ=0.05Pa・sの液体を計測する。トルクの計測において  円筒上下の境界や機械的な摩擦の影響がないものとして、(2)で得られた速度分布をもとに粘性 せん断応力から得られる軸トルクTの値を求めよ。なお円筒の軸方向長さL=50mmとし、せん断応 力は         τ_rθ=μ(d(u_θ)/dr-u_θ/r) で与えられるものとする。 (自分の答え) (1)     Re=u_θd/ν=40より     u_θmax=40×10^-5/(0.03-0.02)=0.04 u_θmax=R_iΩより       Ωmax=2rad/s (2) 常微分方程式を次のように変形する        r^2u"+ru'-u=0 ここでu=r^λとおいて代入して得られる特性方程式は        λ^2=1         したがって一般解は        u_θ=C_1r_^-1+C_2r       境界条件r=R_iのときu=ΩmaxR_i,r=R_oのときu=0より        C_2=Ω_maxR_i^2/(R_i^2-R_o^2)=-1.6 C_1=-R_o^2C_2=1.44×10^-3 u_θ=(1.44×10^-3)r^-1-1.6r (3)     u'=(-1.44×10^-3)r^-2-1.6 ニュートンの粘性法則より          τ_rθ=μu'     よってトルクTは      T=2πR_iL×(τ_rθ|r=R_i)×R_i=-9.04×10^-6[N・M]     

  • 流体力学の問題

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

  • 円管内の流体の温度分布について

    円管内に粘度の高い樹脂を流す場合についての質問です。 円管の管壁を均一に加熱しているとし、また、円管に流入する樹脂の流量も一定とします。念のため、円管に流入する樹脂の温度<円管の管壁温度です。 このとき、中を通る流体の円管断面での流速分布として、 円管中心と、管壁付近ではどちらの流速が大きいのでしょうか? また、円管内に流入する樹脂の流量が微妙に変動した場合 円管中心と管壁では流速の変化の割合が異なるでしょうか。 もし円管断面での温度分布が無ければ、せん断の影響で中心の方が流速が大きいと考えますが、一方で円管断面の中心の温度が低いので、粘度が上がって端部のほうが流速が大きくなることも考えられるため、判断に迷っております。

  • 流体力学の問題がわかりません><

    外半径がa,内半径がbの水平な同心円管のあいだを流れる流体(粘性率η)について、管の長さをl,両端の圧力差をΔpとすると、単位時間に流れる流体の体積は次式で与えられることを示せ。 πΔp/8ηl{a^4–b^4–(a^2–b^2)^2/log(a/b)} どうかよろしくお願いいたします。

  • 流体力学の質問

    質問1 以下はネットから引っ張ってきたんですが、流体の粘性応力は考えなくていいんですか? 「直径Dの円管について長さLの区間を考えます。 この区間の圧力損失をΔpとします。 また、円管壁面での壁面せん断応力をτwとします。 この時、運動量保存式を考えると πDLτw=πD^2/4Δp となる。」 質問2 管摩擦係数は壁面の粗さにもよるはずなのに、ブラジウスの近似式では、レイノルズ数しか出てこないのはなぜでしょうか? 質問3 ある円管で、壁面でのせん断応力は、μdu/dy(μは流体の粘性係数、uは流体速度、yは壁面からの距離)の、y=0での値で正しいですか? 多くてすみませんが、宜しくお願いします。

  • 流体力学についての質問なのですが・・・

    流体力学についての質問なのですが・・・ 水平に設置された円管内を非圧縮性流体が一定の流量で流れていて、円管の途中で直径が半分になっていたとき 管内平均流速・最大流速・壁面せん断応力・レイノルズ数は何倍になっているか、という問題なのですが 円管直径をdとすると、ハーゲンポアズイユの式から 平均流速 Um(R=d)=-d^2/32*dp/dx um(R=d/2)=-d^2/32*1/4*dp/dx=1/4*Um と解答を導いて1/4倍という答えを求めたのですが、よく考えてみると円管を細くすると流速は上がると思うのですが計算で導けなくて困っています。 この解答では間違っているのでしょうか? せん断応力Τ(R=d)=d/4*ΔP/L τ(R=d/2)=d/8*ΔP/L レイノルズ数 Re(R=d)=Um・L/v Re(R=2/d)=1/4*Um・L/v (自分の答え) 平均流速 1/4倍 最大流速 1/4倍 せん断応力 1/2倍 レイノルズ数 1/4倍

  • 流体力学/円管~円筒隙間の流速

    内径r1の円管の中に、外径r2の円筒が入っており、 これらは同軸で、且つ軸は鉛直方向を向いている状態です。 この円管~円筒の隙間に密度ρ、粘度ηの流体を上端から注入した時、 この流体の鉛直方向の流速(単位時間当たりの体積)を求めたいのですが、 どのように算出すれば良いでしょうか。 円管と円筒の上端/下端位置は同一。また、上端/下端は大気開放です。 また、注入時の初速は与えず、イメージとしては、上端に広がった水平面に 流体が広がっていて、そこから自然に隙間に流れ込むイメージです。 ハーゲンポアズイユの式の変形で行けるような気がしてますが・・・。

  • 非ニュートン(指数則)流体の圧力損失の求め方

    現在、プラスチック用金型の流路設計を行っております。 その際に、流路の圧力損失により金型の形状を決めるのですが、その圧力損失の計算に苦慮しております。 流路を流れる流体が水のようなニュートン流体の場合の圧力損失の計算式は、 (円管) ?P=8ηLQ/πR^4 (平行平板) ?P=12ηLQ/wh^3 ?P:圧力損失(Pa) η:見かけ粘度(Pa・s) L:円管および平行平板長(mm) Q:流量(mm^3/sec.) R:円管半径(mm) w:平行平板幅(mm) h:平行平板間隙(mm) 以上の計算式で求めることができると思います。 プラスチック、すなわち樹脂は非ニュートン流体ですので、基本計算式はこの通りだと思いますが、せん断速度により見かけ粘度が変化するので、この計算式をそのまま使ってはいけないのではないかと考えております。 どなたか、非ニュートン流体の圧力損失の導き方をご存知でしたらご教示下さい。 宜しくお願い致します。

  • 同心二重円筒間の流れが理解できません・・・

    初投稿させて頂きます。 ここ3、4日ずっと考えていますが全く分からず投稿しました。 宜しく御願いします。 内円筒(半径r1)と外円筒(半径r2)の間に、非ニュートン流体が満たされています。 外円筒を静止させて(角速度ω2=0)、内円筒だけω1の角速度で回転させました。 このとき、mおよびnがべき乗(指数)則定数であることを示せという問題なのですが…、 答えとして、 γ={(r2/r)^(2/n)}・【(2・ω1)/n[{(1/kr)^(2/n)}-1]】 となると書かれているのですが、なぜこのような形になるのか全く理解できません。 γは記入の関係上、ドット(・)が付けることが出来ませんでした。 せん断率もしくはひずみ速度と御考えください。 (以下同様です。) 詳細を書きます。          せん断ひずみ率がγ=γ(rθ)で表わせます。【( )内の文字は下付文字です。】  べき乗(指数)則流体は単純せん断流れと仮定します。               ここで、べき乗則流体はτ=mγと表わせます。 mは流体に関しての物質定数です。 宜しく御願いします。

  • 流体力学の問題で…

    非圧縮性流体の中を球(半径R)が速度V(ベクトル)で移動しているとき流体の流れはどうなるのかという問題で、ポテンシャル流だとしてそのポテンシャルφはラプラス方程式(球座標)を満足し、rとVの関数になり、無限遠で0になるという条件を満足しなければならないんですが、参考書では(ランダウさんの流体力学1、非圧縮性流体の章)φはCv・grad(1/r)とおいてるんです。ラプラス方程式の解だからφを1/r^a(aは自然数)とするならわかるんですが。grad(1/r)はなぜなんでしょう?