pythonで簡単にロケットエンジンの性能計算：RocketCEAの使い方 (Windows編)

はじめに

【SpaceX】特集 <第1.5回> :ロケットの性能についてロケットの性能や打上能力は具体的にどういう指標なのか、どのように決まるのか、について解説します。...

みなさんこんにちは。前回記事の続きで、ロケットエンジンの性能を計算する方法を考えていた時たまたま発見した超便利な“RocketCEA”について、今回はその導入方法と使い方を簡単に説明します。

管理人も発見当初はサブ機のkubuntuで粘っていたのですが、ブログ執筆のためについにwindows版を導入しました。

CEAとは・RocketCEAとは

CEAとはNASAが開発した “Chemical Equilibrium with Applications” という化学平衡計算プログラムのことで、様々な化学種の燃焼計算や、計算結果に基づいたロケットエンジン燃焼器の性能計算ができます。fortran77で書かれておりフリーである反面、長年使い勝手が悪かった上、安保上の理由から最近は無償公開がされなくなってしまいました。

現在ではwebアプリとしてのみ使える状態です。

https://cearun.grc.nasa.gov/

しかし近年、消えてしまったと思われていたこのスタンドアローン版のCEAが有志の手で復活を遂げました。しかも使いやすいpythonのモジュールとして！公式ページはこちらから：

https://rocketcea.readthedocs.io/en/latest/

試した方法

本モジュールに最初に着目したqiita記事はこちらです。linux環境、mac環境はこの通りに導入すれば問題なく動きます。

インストール方法

実際に導入した手順です。.NET Frameworkが壊れてたり、おかしなバージョンのgccが入っていなければこの手順通りにやれば多分動くと思います！

トラブルがあった場合、可能な限りサポートいたします！

動作環境

OS

管理人がCADとクリスタ、他Steamなどゲームで使っているメインPCです。

ソフト環境

python: 3.10.1

その他競合しそうな環境：Git Bash, Anaconda

pythonの準備

実行に必要なライブラリです

https://rocketcea.readthedocs.io/en/latest/quickstart.html

numpyとmatplotlibを入れます/最新版に更新します。

pip install numpy
pip install matplotlib

pip install --upgrade numpy
pip install --upgrade matplotlib

gfortranの導入

rocketCEAを使うためにはfortranのコンパイラが必要です。ここでgcc gfortranを使います。linuxだと一瞬なのですが、windowsではここでいくつかコツがあります。順番通りにやって行きましょう。

一応の公式ガイド (これが原因で泥沼に)

https://rocketcea.readthedocs.io/en/latest/installgfortran.html#link-installgfortran

mingw64のダウンロード

下記のページにアクセスします。

https://sourceforge.net/projets/mingw-w64/

ここにある”ダウンロード”は　してはいけません。

動かないバージョンがdlされます。

Filesのタブから下にスクロールし、MinGW-W64 GCC-8.1.0から “x86_64-win32-seh”をダウンロードします。

ちなみにRocketCEA公式ページで紹介されている手順ではこのようなエラーが出てしまいます。

mingw64の解凍

無事にダウンロードができたら解凍します。

解凍フォルダmingw64の中身はこんな感じです。

解凍が完了したらmingw64フォルダごとCドライブ直下などわかりやすい場所に移動します。

環境変数の設定

下記にPATHを通す必要があります。コマンドプロンプトから下記の通り追加するかコントロールパネルを使います。

1. cmdで通す方法

set PATH=C:\MinGW\mingw64\bin;C:\MinGW\mingw64\lib;%PATH%

2. コントロールパネルで通す方法

コントロールパネルを開き、“環境変数”で検索。

“環境変数を編集”

環境変数の画面を開いたら、“システム環境変数”の“Path”を編集、下記を追加します。

C:\MinGW\mingw64\bin

C:\MinGW\mingw64\lib

念の為、私はC:\MinGW\mingw64も追加しました。

動作確認

以上で”インストール”は終了です。パスを通したことで直接コンパイラプログラム(gcc.exe)を叩くため特にほかの作業をしなくても動くはずです。システムを再起動し、gfortranが使えるようになっているかサンプルコードで確認します。

#test.c
#include <stdio.h>

main()
{
  printf("Hello World\n");
}

gccでコンパイラを実行してみます。

何やらエラーが出てきましたがとりあえず無事に動くことが確認できました。

RocketCEAの導入

gccが入ったらrocketcea本体を導入します。pythonのライブラリとして公開されています。

pip install rocketcea

以上です。

下記から直接ソースを落として使っても良いです。

公式ガイドの通り動作確認を行います。これで動けば導入は完了です。

python -c "from rocketcea.cea_obj import CEA_Obj; C=CEA_Obj(oxName='LOX', fuelName='LH2'); print(C.get_Isp())"
374.30361765576265

次ページ(使い方)へ

ABOUT ME

Trude (まきブロ管理人)

機械屋さん。多くの人に機械の面白さやエンジニアリングの楽しさを知ってもらうべく、解説や紹介記事を発信しています。

BLOG：https://makkiblog.com/