Python/PIVの検証用に管内の粒子流れ動画を作ってみた

  • このエントリーをはてなブックマークに追加

粒子をトラッキングして流れの速度場を計測するPIVというのを学びたく、まずはサンプルの管内流れ動画をPythonで作ってみました。これで良いのかどうかはわかりませんが、Pythonによる流れ動画の作り方を備忘録として残します。

こんにちは。wat(@watlablog)です。本日はPIVコードのテストに使う粒子流れ動画を生成する案を記録します

PIVコード検証用のサンプル動画

PIVとは?

PIV(Particle Image Velocimetry)とは、粒子画像流速測定法を意味し、流れ場にトレーサ粒子を配置させ、粒子の速度を計測する技術です。

トレーサ粒子は液体であればまさに粒子、気体であれば人体に無害な煙状の物体を使うことが一般的です。

以下にPIVで行われている事をざっくり図で示します。
PIVは高速度カメラ等で動画(複数の画像)を撮影します。各画像間で何秒経っているかは撮影時の条件から自明なので、各画像内のトレーサ粒子の移動量を知る事ができれば、速度を算出する事が可能です。

PIVのイメージ図

簡単そうに見えますが、実際の撮影画像から精度良くトレーサ粒子の移動量を抽出するためには様々なノウハウが必要です。餅は餅屋という事で以下にPIV専門の会社様のリンクを記載しておきます。

外部リンク:カトウ光研株式会社:PIVとは

流体分野の実験について

機械工学の世界では航空機、自動車、その他産業機械は屋外や水中等で高速で移動したり、様々な流体の中で使用される機械を扱います。これら流体力を受ける機械は、場合によっては本来の性能を出せないことが多々あります。

そのため製品形状に対して数値シミュレーションを使って設計を行ったりするのですが、流体の解析は適切な壁の設定、乱流モデル、境界条件、メッシュ形状…と様々な要因で結果が変わるので、検証目的としても実験による流れの把握は重要です。

また、流体シミュレーションを精度良くやろうと思った時には、場合によっては億単位のセル数を必要とし、スパコンで何日も計算することもあります。実機があるなら実験の方がはやいというケースもあるのがこの分野の特徴ですね。

このように、計算における科学的発展が目覚ましい現代においても、流体分野は(他の分野もそうですが)理論流体力学や数値流体力学と共に、実験流体力学も重要な役割を担っています。

サンプル動画を作ってみたきっかけ

この記事を書いた理由は、勉強のために自分でもPIVのコードを作って使ってみたかったからです。そのためにはコードを試す撮影画像が必要なのですが…個人で実験設備を揃えるのは至難のワザなので、粒子の流れ画像をなんとなく作ってしまおうと思いました。

この記事で作成するサンプル動画

この記事では以下のようなGIF動画を作成します。

サンプルの粒子流れ動画

黒い点の一つ一つがトレーサ粒子を模擬しており、初期配置は完全ランダムにしてみました。
また、粒子の大きさも時刻によってある範囲で変わるようにしており、数値調整次第では途中で消えたり出てきたりといった挙動もできるようにしてあります。

図の\(y\)軸中心の流速が速く、\(y=\pm 2\)では遅いという管内流れの特徴も模擬しています。

流体は粘性を持ちますので、壁と流体の摩擦から壁付近では流速が遅くなります。実際は管径や粘性、レイノルズ数によってこの速度分布が変わりますが、今回のプログラムではここまで表現していません。ご興味のある方はハーゲンポアズイユ流れとかでGoogle検索してみて下さい。

以下に速度分布を同時プロットした動画を示しますが、今回はscipyのnorm.pdf(確率密度関数)を使って速度分布を模擬しました。実際は管のシミュレーションとかして求めるのが良いのでしょうけど、まだ未熟なためそれは後日。

速度分布を図示

norm.pdfは「Pythonでガウス分布を持つノイズの作り方と調整方法」でも使っていますので、是非参考にしてみて下さい。

Pythonで粒子流れ動画を作るコード例

Advertisements

全コード

全コードを以下に示します。GIF動画を作るcreate_gif関数と粒子画像を作るgenerate_particle_movie関数で構成されています。詳細はコード内コメントに記載しましたのでご確認下さい。

generate_particle_movie関数はwx(横幅), wy(高さ), sf(粒子密度のスケールファクター), frame(フレーム数), center(流れ中心), sigma(速度分布の標準偏差), vsf(速度のスケールファクター)を引数にしています。

是非色々な流れ場を作って遊んでみて下さい。

色々な流れ場を作ってみる

デフォルトの流れ

デフォルトは以下のコードです(関数実行時に何も指定しなくても動くようにしました)。

デフォルトの流れ

一様流

速度分布のsigmaとvsf(速度に対するスケールファクター)を変更して一様流を作ってみます。

 

一様流

グラフの上下端で速度変化が少なくなりました。流体シミュレーションで言えばスリップ境界になったようなイメージ。

センターオフセット流れ

流速の速い部分を\(y=1\)とオフセットさせたセンターオフセット流れ(ネーミングは一般的ではないと思いますが)を以下に示します。

コード変更点は以下。centerの値を0以外にします。ここでは流れの速さをわかりやすくするためにvsfも変えています。

センターオフセット流れ

軸を消して撮影画像のようにしてみる

画像解析で使うために、グラフの情報(軸の値、ラベル等)は要らないかも知れません。以下のコードのようにラベル定義をコメントアウトし、作ったtickの設定をFalseにする事で消す事ができます。

こんな感じ。

軸を消してみた結果

枠線すら要らない、という場合はもっと簡単で、.axis('off')を書けば良いだけです。

こんな感じ。これならそのまま解析できそうですね。

まとめ

PIVという流体の実験計測、分析技術をちょっと学んでみたく、まずはサンプルの流れ場の動画を作成する事をしてみました。

あまり需要の無い記事だとは思うので完全に自分用の備忘録のようなものです。

本記事ではPIVの概要と流れを模擬して動画を生成するPythonコードを紹介しました。

パラメータ変更で色々な流れ方をするようにできたと思います。

この動画で本当にPIV解析ができるのかどうかはまだ謎ですが、やってみて不具合があったらこのページを修正してみようと思います。

本当にPIV解析ができるようになるのか…は今後の勉強次第!
Twitterでも関連情報をつぶやいているので、wat(@watlablog)のフォローお待ちしています!

  • このエントリーをはてなブックマークに追加

SNSでもご購読できます。

コメント

コメントを残す

*