戻すのを忘れないでください。中には、コロナホール由来の高速風領域とそうでない低速風領域が混在しています。観測されることが知られています。シミュレーションで示しました。1次元プロファイルです。その真ん中に磁場の弱い領域ができていることがわかります。相対論的速度でジェットが吹き出していると考えられています。よって報告されました。もとで現れる断熱加速です。バルク速度を見ていることを指摘しています。渦問題を収録しています。標準問題として広く使われています。コメントアウトして、計算してみてください。非物理的な効果が累積して途中で計算が破綻するはずです。速度差がある2つの領域の境界層では、大きな渦が成長することが知られています。いう流体力学の不安定です。渦の成長を見てみましょう。)の中では、以下のような設定を用いています。関数で設定しています。下側密度に相当するパラメーターαを使って、上下に密度差を入れています。倍の密度差ができます。磁場は面に垂直方向を向いています。以上に設定してください。不安定)が出来ています。成長の仕方が大きく変わります。)に変えてみましょう。このとき、横向き磁場の張力が上下方向の運動にブレーキをかけようとします。渦の成長がどう変わるのか、計算して確かめてみましょう。磁力線が激しくつなぎ変わる「磁気リコネクション」が起きます。磁気リコネクションを取り扱うことができます。周囲のプラズマにぶつかっていきます。速度程度であることがわかります。簡単のため、上流プラズマの音速〜典型的な音速とします。以下の関係が成り立ちます。リコネクションジェットのマッハ数が1を超えてしまうことがわかります。周辺に生成される衝撃波構造を議論しました。使ったコードの改良版をそのまま収録しています。以下のような磁場反転構造を用いています。プラズマベータ(磁気圧に対するプラズマ圧力の比)を上流のインフロー領域で評価したものです。簡単のため、系全体で温度($T=p/\rho$)を一定にしています。0)付近では磁気圧が働かないため、代わりにプラズマの圧力を増やして全体構造を支えています。このような構造を「プラズマシート」あるいは「電流シート」と呼びます。程度の弱い磁場を追加して、原点近傍でリコネクションを起こしています。基本的なケースでは、対称性を利用して計算量を節約することができます。第一象限のみを解く並列コードです。上の図は、デフォルトの設定で計算した磁気リコネクションのアウトフロー速度です。リコネクション点から細いジェットが吹き出しています。下流側に衝撃波や不連続面を作っているのがわかります。した物理に興味のある方は「生存圏研究」の解説記事をご覧ください。これらの結果をまとめたものが、下の2016年版「プラズモイド・ダイアグラム」です。未解明の基本問題が残っています。著者の腕次第で研究成果を出すことは可能です。コードが皆様の研究をアシストできれば幸いです。何箇所かコードを修正する必要があります。画面が真っ黒になってしまうことがあります。ソフトウェアレンダリングモードで動作させる必要があります。いうオプションを追加すると良いでしょう。このオプションは、付属プログラムのコメント部分に記述しています。見つからない、というエラーが出る。あれば、以下のライブラリを追加インストールしてください。場合は、西田さんのメモが参考になると思います。いうエラーのため、ウィンドウが表示されない。こちらかこちらを参考に対処してください。ノートブックで表示されるグラフを大きくしたい。documentation)。グラフの大きさを任意に変更できるようになります。並列版のデータ形式は?シリアル版と同様、1つのファイルにデータを書き出します。並列プロセス毎にファイルを書き出します。各ファイルに1グリッドずつののりしろが含まれていることを考慮してください。汎用言語ですから、いろいろな解析ができるはずです。データを読み込んだ後は、ぜひ自分の手でプログラミングしてみてください。のりしろは1グリッドで足りるのでしょうか?はい。補間した基本量をノード間転送で調整しています。この過程で2グリッド先の情報が伝わっています。mainp2.f90)の先頭で定義しています。...)で定義しています。専用のサブルーチン(set_eta.f90)等で定義しています。物理的な意味は特にありません。mpibc*.f90(並列版)で設定します。サブルーチンを用意していますので、サンプルコードを参考にしてください。コードを書く必要はありません。境界条件には特に注意してください。エラーを累積した後で結果を変えるので厄介です。複数の境界設定ファイルを利用していました。ノード内通信を共有メモリのデータコピーで代用します。これは実験的な機能です。(experimental)3次元計算はできるのでしょうか?いまのところ、できません。配列の次元を増やすだけで、そのまま3次元計算できるはずです。腕に覚えのある方は、コードを書き換えて3次元化してみて下さい。特性変数を補間した方が良いでしょう。基本量の補間で代用しています。基本量の違いはあまり気にならないと思います。実質精度を上げることを目指しています。リミターを採用しています。下線のものをお使いください。広がる外側境界を考えます。法の4状態の仮定が成立しなくなることがあります。問題に関係しています。ソルバーは非線形ソルバーなので、ガリレイ不変ではありません。最後の引数は、ただの整数では?はい。コードを読みやすくするために使われています。27分(HLL)〜33分(HLLD)かかりました。性能に応じて解像度を調整すれば良いでしょう。このコードは実際の研究に使われているのでしょうか?はい。国内外のいくつかの研究グループに使っていただいています。開発者とは独立に進められた研究の成果です。宛に問い合わせてくだされば、常識的な範囲内で対応します。サポートの必要度合いによっては、共同研究という形で進めた方が良い場合もあります。その場合は、共同研究を申し込んでくだされば幸いです。紹介してくださると嬉しいです。文面の例を以下に示します。言及しても良いと思います。紹介していただければと思います。CFL数を小さめ(0.5以下)に設定しているのはなぜでしょうか?作者の好みです。両側のセル境界からやってきた波が真ん中でぶつかり、直感に合わないと思うのです。良いそうですが、私はどうも納得がいきません。データファイルの番号にあわせて変更してください。途中データのファイル番号を自動検索します。新しいバージョンに更新した方が良いでしょうか?毎年コードの改良を重ねています。更新することをお勧めします。フラックスの部分です。並列版コードが同梱されています。変数で並列プロセス数を指定してください。コマンドで並列プログラムをコンパイルできます。ap2.out)ができるはずです。プログラムの実行方法は、例えば以下のようにします。オプションは並列プロセス数を意味しています。変数のプロセス数と一致していなければなりません。コマンドはプログラムの(標準)出力を画面に表示しつつ、ファイル(output.log)に記録します。並列計算が思ったほど速くない。並列が思いがけず有効になっていることがあります。コンパイラのオプションをよく確認してください。オプションを追加してください。プログラムが高速に動作するようです。使って並列スレッド数を設定します。./a.out全体の環境変数を変えたくなければ、以下でも良いでしょう。オンになってしまうことがあります。オフにすることができます。プログラムが異常終了する。シェルのスタックサイズの制限を解除してください。使うことができるそうです。使うと計算結果がおかしい/プログラムが異常終了する。利用時に問題が多いです。別のコンパイラを使ってください。タイムに注目してください。自分で検証してみると良いと思います。オプションを追加します。次のようなジョブスクリプト(test.sh)を用意します。ジョブを投入してください。慣れていくと良いと思います。全ての言語に対応しています。モジュールを追加インストールすれば大丈夫です。途中画面を出力したままログアウト(デタッチ)して、後で再び接続(アタッチ)することができます。デフォルトでは少し使いにくいので、以下のサイトを参考にカスタマイズしても良いと思います。ファイルを編集して行います。同じぐらい紛らわしいのは承知しています。最後まで目を通してくださり、ありがとうございます。インストールすることができます。類似の方法でインストールできると思います。インストールされます。ライブラリのインストール」を参考にインストールしてください。最初からインストールされています。インストールすれば良いと思います。コマンドを使えば良いようです。参考になると思います。開発ツールが一通り入っています。インストールしてください。使ったプログラムのコンパイル・実行に影響が出ることがあります。できるだけ最新系列のものをインストールするようにして下さい。tools)を追加インストールする必要があります。次のコマンドを試してみてください。ダイアログの指示に従ってください。それ以上の操作は必要ありません。コンパイラは有償製品です。日本ではエクセルソフトが販売代理店になっています。Editionを無償でダウンロードして利用することができます。パッケージ集でインストールすることができます。以下のどちらかのサイトからファイルを入手することをお勧めします。標準のインストーラーパッケージを配布しています。長年、コンパイラのバイナリファイルを配布してくださっているサイト。バイナリをダウンロードできます。後者を使う場合は、アーカイブファイルを以下のように手動で展開します。以上の最適化オプションを使うようにしてください。プログラムが速く動作するようになります。並列計算を行うことができます。併用したハイブリッド並列が有効です。スパコンと同じ並列条件でプログラムを開発・実行することができます。コンパイル・インストールする必要があります。自分でビルド・インストールすることをお勧めします。それぞれ別のディレクトリにインストールすることを考えます。オプションを指定しなくても結構です。ファイルがインストールされます。インストールするには以下のようにします。処理には時間がかかりますが、気長に待ってください。並列コンパイルに対応しています。オプションを使ってください。インストールした手順です。オプションを調整してください。インストールするには以下のようにします。組み合わせてインストールした方が良いでしょう。インストールした手順です。コンパイラと組み合わせる必要があると思います。確認することができます。ディストリビューションを使って、インストールすることができます。使う方法を紹介します。パッケージをダウンロードして、インストールします。--「Download」タブからファイルをダウンロードできます。ライブラリをインストールしていきます。インストールしてみましょう。インストールすることをお勧めします。便利に作業できるようになります。インストールしてください。オプションを使います。シェルを立ち上げてみましょう。確認してみてください。Doドキュメント利用事例、スパコンのジョブスクリプトの例を紹介。初期設定方法(特に並列版)。問題のベンチマーク結果。よるアクセラレーター対応。集団隣接通信(MPI_Neighbor_*)のテスト。ハイブリッド並列性能の向上。法の修正理論と論文化。2020-xx-yy2019-12-27コメント・ドキュメントの微修正。ノートブックを追加しました。よるインタラクティブ操作にも対応。(和田智秀さんに感謝)ムービー作成用のサンプルスクリプトを追加しました。ファイル番号に加えてファイル名を指定できるようにしました。2019-01-25途中データからのリスタート方法を変更しました。とき、ファイルを自動検索します。計算領域を多次元分割するようにしました。2013)を選択可能に。ファイル構成を整理しました。よるスレッド並列性能を改良しました。通信にノンブロッキング通信を採用。ソルバーのバグを修正しました。その他、多くの微修正。2017-07-28相対論ジェットの境界層問題を収録しました。コード(垂直伝搬専用)を収録。可視化ルーチンを改良。2017-04-28データのファイル形式を変更しました。使わないようにしました。(近藤光志さん、堀田英之さん、松本倫明さんに感謝)収録課題の解説をひととおり用意しました。可視化ルーチンを改良。サンプルノートブックを一時的に削除。可視化ルーチンを改良。サンプルノートブックを追加。ディレクティブを追加。ソルバーを統合・高速化。