解析手法は、回転物体周りのストークス流れ場を解き、表面せん断応力を計算し、その応力を積分してトルクを求めることです。 半径 (R) の浸漬球が、粘度 (\mu) のニュートン流体中で角速度 (\omega) で回転する場合、低レイノルズ数の解から抵抗トルクは次のように求められます。

[ G=-8\pi\mu\omega R^3 ]

ここで、負号は流体力学的トルクが回転に逆らうことを示します。

クリーピング流れでは、トルクは完全に粘性応力によって決まります。すなわち、回転する表面の境界条件を用いて線形ストークス方程式を解き、表面における接線方向応力を評価し、そのモーメントを物体全体にわたって積分します。

ストークス流れ法が適用できる理由

流体慣性が無視できる

この方法は、粘性項と比較して慣性項を無視できるほどレイノルズ数が十分に小さい場合に適用できます。支配的な運動量方程式は次のようになります。

[ -\nabla p+\mu\nabla^2\mathbf{v}=0, \qquad \nabla\cdot\mathbf{v}=0 ]

静止した流体中で球が回転する場合、流れは軸対称で、純粋な方位角方向の流れとなります。

トルクと速度の線形関係

ストークス方程式は線形であるため、トルクは角速度に正比例します。また、流体粘度および球半径の3乗にも比例します。

[ |G|=8\pi\mu\omega R^3 ]

この比例関係は、並進する球に対する線形抗力則の回転版です。

回転流れ場を解く

境界条件を適用する

球の中心を原点とする球座標を使用します。軸周りに回転する球について、速度条件は次のとおりです。

[ v_r=0, \qquad v_\theta=0, \qquad v_\phi(R,\theta)=\omega R\sin\theta ]

最後の条件は、球表面の接線方向速度を表します。球から十分遠方では、流体は静止しています。

[ v_\phi(r,\theta)\rightarrow 0 \quad \text{as} \quad r\rightarrow\infty ]

方位角方向速度を求める

これらの条件のもとでストークス方程式を解くと、次式が得られます。

[ v_\phi(r,\theta)

\frac{\omega R^3\sin\theta}{r^2} ]

速度は (1/r^2) に比例して減衰し、クリーピング流れにおける回転擾乱が局所的であることを示します。

流れの分布を解釈する

流体速度は (\sin\theta=1) となる球の赤道で最大となり、回転軸上ではゼロになります。流れ場全体は、回転する表面と周囲流体との粘性結合によって生じます。

表面せん断応力の計算

接線方向応力を評価する

軸対称な方位角方向流れでは、関連する球座標系の応力成分は次のとおりです。

[ \tau_{r\phi}

\mu \left( \frac{\partial v_\phi}{\partial r}

\frac{v_\phi}{r} \right) ]

速度場を代入し、(r=R) で評価すると、次式が得られます。

[ \tau_{r\phi}(R,\theta)

-3\mu\omega\sin\theta ]

負号は、流体応力が球の回転に逆らうことを示します。

トルクアームを考慮する

極角 (\theta) にある表面要素から回転軸までの垂直距離は次のとおりです。

[ R\sin\theta ]

したがって、微小トルクは、表面要素に作用する接線方向力にこのモーメントアームを乗じたものとなります。

トルクの積分

球全体にわたって積分する

表面積要素は次のとおりです。

[ dA=R^2\sin\theta,d\theta,d\phi ]

したがって、全トルクは次式で表されます。

[ G

\int_0^{2\pi} \int_0^\pi \tau_{r\phi}(R,\theta) (R\sin\theta) R^2\sin\theta ,d\theta,d\phi ]

せん断応力を代入すると、次のようになります。

[ G

-3\mu\omega R^3 \int_0^{2\pi}d\phi \int_0^\pi\sin^3\theta,d\theta ]

次式を用いると、

[ \int_0^{2\pi}d\phi=2\pi, \qquad \int_0^\pi\sin^3\theta,d\theta=\frac{4}{3} ]

次の結果が得られます。

[ \boxed{G=-8\pi\mu\omega R^3} ]

必要な印加トルクは、同じ大きさで符号が反対になります。

[ \boxed{G_{\text{applied}}=8\pi\mu\omega R^3} ]

トルクと動力の関係

回転を維持するために必要な粘性動力は次のとおりです。

[ P=|G|\omega ]

回転する球については、

[ \boxed{P=8\pi\mu\omega^2R^3} ]

これにより、流体中で散逸する機械動力を直接見積もることができます。

単位操作への方法の拡張

回転抵抗係数を用いる

その他の浸漬物体についても、同じ手順を回転抵抗係数を用いて次のように表せます。

[ G=-\zeta_r\omega ]

球の場合は、

[ \zeta_r=8\pi\mu R^3 ]

この係数は、物体形状、拘束条件、壁面への近接、および他の移動表面の存在によって変化します。

攪拌機に適用する

低速で高粘度の単位操作では、インペラまたは回転攪拌機に必要なトルクを、クリーピング流れ場を解くこと、または適切な抵抗係数を求めることによって算出できます。得られたトルクは、せん断応力、動力消費、およびモーター要件の計算に利用できます。

局所量と全体量を区別する

応力分布は、表面全体にわたる局所的な流体荷重を表します。積分されたトルクは、装置のサイジングおよび動力計算に必要な全体量です。

トレードオフを理解する

低レイノルズ数の仮定には制約がある

解析結果は、流体慣性が無視できることを仮定しています。レイノルズ数がもはや小さくない場合、対流慣性によって速度場と応力場が変化し、(G\propto\mu\omega R^3) の関係は適切でなくなる可能性があります。

球の結果は普遍的ではない

係数 (8\pi) は、無限に広がるニュートン流体中で孤立した球が回転する場合に特有のものです。この係数を、インペラ、円柱、壁面近傍の粒子、または拘束容器に直接適用してはなりません。

ニュートン流体の挙動を仮定している

導出では、一定粘度およびニュートン流体の構成則を仮定しています。高分子溶液、懸濁液、生体流体、その他の非ニュートン流体では、修正された応力モデルと異なるトルク・速度関係が必要になる場合があります。

境界効果は重要になる可能性がある

壁面や近傍表面は速度場を変化させ、流体力学的抵抗を増加させます。そのため、パイロットプラント設備では、壁面補正係数、数値的なストークス流れ解、または実験的キャリブレーションが必要になる場合があります。

目的に適した選択をする

球に対する解析解は、基準および実用的な初期計算として有用ですが、設備に適用する前に形状と流動領域を確認する必要があります。

主な目的が孤立した回転球のトルク推定である場合： (G=-8\pi\mu\omega R^3) を使用します。符号は、加えられた回転に対する抵抗を表します。
主な目的が低速攪拌機のサイジングである場合： 球の係数を直接適用するのではなく、攪拌機の形状に固有の回転抵抗を求めます。
主な目的が局所せん断の予測である場合： 全トルクを積分して求める前に、速度場と表面応力分布を使用します。
主な目的がクリーピング流れの範囲外での運転である場合： 流体慣性を考慮し、適切な高レイノルズ数相関式、数値モデル、または実験を使用します。

基本原理は単純です。粘性流れ場を解き、表面応力を評価し、そのモーメントを積分して流体力学的トルクを求めます。

まとめ表：

主な項目	結果または適用
支配モデル	線形クリーピング流れのストークス方程式
球のトルク	`G = -8πμωR³`
回転抵抗	孤立した球では `ζᵣ = 8πμR³`
必要動力	`P = 8πμω²R³`
計算方法	流れを解き、表面応力を評価し、そのモーメントを積分する
主な制約	形状、壁面、慣性、および非ニュートン流体には補正または異なるモデルが必要

お客様の用途に適した信頼性の高い低速混合・流体プロセス設備をお探しですか？ LABPARKにお問い合わせください。LABPARKは、化学工学、バイオプロセス・バイオテクノロジー、環境・水処理向けの教育・職業訓練用単位操作パイロットプラントを提供しています。当社のソリューションは、実践的なトレーニング、研究、プロセス開発、スケールアップ向けに設計された用途重視のパイロットプラントシステムにより、大学、研究機関、企業がトルク、せん断、消費電力、および設備性能を評価することを支援します。