Pythonで球体周りの抗力係数のパラメタスタディ

こんにちは('ω')ノ 本日はOpenFOAMでの流体周りの抗力係数を確認してみました。 黒実線が実験値の近似曲線で、青丸が今回のシミュレーションの結果です。 横軸にレイノルズ数と縦軸に抗力係数の常用対数をとったものです。 抗力係数：$${C_{d}=\frac{F_{D}}{\frac{1}{2}\rho U^2A}}$$ レイノルズ数：$${Re=\frac{UL}{\nu}}$$ そこそこ試験結果と近しい値をとりました。 全てのレイノルズ数領域で定常解析かつ乱流