[論文レビュー] Dynamics of a Tube Conveying Fluid

Shuhei Shima, Tamotsu MIZUGUCHI|ArXiv.org|May 16, 2001

Vibration and Dynamic Analysis被引用数 5

ひとこと要約

本稿では、一般の境界条件における流体を運ぶチューブの動的挙動をモデル化する非線形積分微分方程式を構築した。その結果、固定端のチューブはホップ分岐に起因するフラッターを示し、ピン支持のチューブはピッチフォーク分岐に起因する回転不安定性を示すことが明らかになった。数値シミュレーションにより、これらの不安定性が確認され、系が連続体極限に近づくにつれて理論的臨界値に収束することが示された。

ABSTRACT

A tube conveying a large amount of fluid with a free outlet does not sit still. We construct and analyze a nonlinear evolution equation describing such phenomena. Two types of boundary conditions at the inlet are considered, one for which it is clamped and one for which it is hinged. Analyzing the linear stability of the trivial solution, we find that with the former boundary conditions, it exhibits a ``flutter'' instability, while with the latter boundary conditions, it exhibits a ``rotation'' instability. These instabilities and the nonlinear behaviors of the system are also studied numerically.

研究の動機と目的

• 一般の入口境界条件における流体を運ぶチューブの大きな振幅運動を最小限で物理的に現実的であるモデルを構築すること。
• 自明解（まっすぐな状態）の線形安定性を解析し、境界条件に応じて不安定性の種類（フラッターまたは回転）を特定すること。
• 線形不安定性とその後の非線形動的挙動を捉える非線形発展方程式を導出し、数値的に検証すること。
• 離散モデルの臨界値が N → ∞ に近づくにつれて理論的連続体極限に収束することを確認すること。

提案手法

• N 個のジョイントと N 個のセグメントを持つ離散的アーチレートドモデルを構築し、角度依存の力によって弾性力、流体の流れ、粘性力をモデル化する。
• 固定されたセグメント長を満たすためにラグランジュ乗数を用い、角度 θn に対して N 個の独立した発展方程式を導出する。
• 連続体極限（N → ∞）を取ることで、チューブの形状と動的挙動を記述する非線形積分微分方程式を導出する。
• 自明解の線形安定性解析を実施し、固有値問題を解いて不安定性の臨界パラメータを特定する。
• 有限 N に対する発展方程式の数値的シミュレーションを実施し、分岐の種類と理論的極限への収束を観察する。
• 数値的臨界値（εc(N), Tc(N), εh(N)）を関数 a/(N−b)+c にフィットさせ、N→∞ における外挿を行い、理論的予測の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 入口が固定（位置と方向が固定）されている場合、流体を運ぶチューブにどのような不安定性が生じるか？
• RQ2 入口がピン支持（位置のみ固定）されている場合、系はどのように振る舞い、運動の始まりを支配する分岐の種類は何か？
• RQ3 離散モデルから導出された非線形発展方程式は、不安定化の発生を正確に予測できるか？また、連続体極限において理論的値に収束するか？
• RQ4 セグメント数 N が増加するにつれて、臨界流速パラメータ（εc, εh）および振動周波数（ωc）はどのように変化するか？
• RQ5 慣性効果を含むモデルと含まないモデルの間で、非線形動的挙動にどのような定性的な差が生じるか？

主な発見

• 固定端チューブは εc ≈ 37.69 でホップ分岐を経て、持続的な振動的フラッター運動を示す。
• ピン支持チューブは、ゴルドストーンモードを伴うピッチフォーク分岐を経て εh ≈ 20.16 で、自発的な回転不安定性を示す。
• 数値的シミュレーションにより、臨界値 εc(N)、Tc(N)、εh(N) が N → ∞ に近づくにつれて理論的極限に収束することが確認され、εc(N) → 38.0 ± 0.9、εh(N) → 20.3 ± 0.6 となる。
• 分岐は上昇的であることが示され、臨界点を超えて ε が増加するにつれて変形の大きさが単調に増加した。
• 慣性効果によりより複雑な動的挙動が生じる：固定端チューブではカオスへの周期倍増加ルートが観察され、ピン支持チューブでは2回目のホップ分岐とその後の周期倍増加が発生した。
• モデルは吸引と吐き出しの両方のチューブに適用可能であり、流体力は各点における運動量の変化に起因するためである。