シミュレイト趣味

　
　　＜＜水路に発生した孤立波の運動＞＞

　三角関数のような連続した波ではなく、ひと山づつ独立した安定な波が孤立波です。
　下記のプログラムで、水路の中の孤立波を体験してみて下さい。

　　　　プログラム名：　USFlow.exe　プログラム倉庫へ

（注）「Ｍｉｃｒｏｓｏｆｔ　Ｖｉｓｕａｌ　Ｂａｓｉｃ　２００５」による実行形式での配布のため、「．ＮＥＴ　ＦＲＡＭＥＷＯＲＫ」が機能している環境以外では実行できません。

＜＜１．現象模擬空間の構成＞＞

　水路といっても、幅は２００ｍにしてあるので河川に近い想定です。
　底の勾配は１万分の１、定常状態での水深と流量は２ｍと２８５立方ｍ／秒です。
　距離２００ｍを単位区間として、２０区間で４ｋｍの模擬演算を行います。

　空間表現は、分布を扱う次元で言えば、水路の流れに沿う距離方向の１次元です。
　水路幅方向には流れの一様性が仮定されていますし、高さ方向も水面の水位や水深
　に纏めて扱うことで水圧の垂直方向分布などの取り扱いは避けられています。

　演算上の時間刻みは１回３０秒にしてあり、６０回分＝３０分間の模擬となります。
　ただし、画面のグラフ更新は１秒単位ごとなので実現象の３０倍早い上映です。

＜＜２．模擬された孤立波の波速と波長＞＞

　掲載のプログラム実行によって画面に現れる孤立波は、波速で約４．５ｍ／秒、波長は１２００ｍくらいになるでしょうか。
　定常状態の流速である０．８ｍ／秒をはるかに上回る「長波」の波速になっていますので、この模擬演算は検証の関門をひとつクリアしています。
　なお、波形が徐々に低くなるのは、水路底面での摩擦により運動エネルギーが損失を受けることによるものです。

＜＜３．滑らかな現象表現と、それを可能とする背景＞＞

　ここで注目を頂きたいのが、たかだか６区間程度で孤立波１個を滑らかな流量分布として表現している点です。

　じつは、この点が独創的な手法であり、距離方向の微分演算を単純な差分ではなく自然スプラインを応用することにより近似曲線次数の高い（誤差の少ない）ものとし、波動現象を扱ううえで大切な、計算の安定性も獲得しています。
（差分の誤差が無意味な波動を生み、収拾がつかなくなる傾向が回避されています。）

　このことを直感いただくために、流量分布曲線は自然スプラインが与える各区間ごとの３次曲線で描き、近似曲線が滑らかにつながることも示しています。
　なお、水面の高さである水位の分布グラフは、位置のエネルギーに流速のエネルギーも加えた一括のエネルギー分布しか扱っていないため、水位だけを分離して近似曲線で描くことが難しく、折れ線グラフでの表示となっています。

＜＜４．実現象との対比による意味づけ＞＞

　ところで、実際にこのような現象が起こりうるでしょうか。
　答えはＹＥＳです。
　大きな川では、河口部から数ｋｍ～数十ｋｍ程度上流に潮止め堰が設けられることがありますが、その堰のゲートを故意に急に動かして一定時間放流を増やし、再び急に戻すと、こんなふうな津波とも呼べそうな孤立波が発生してしまいます。

　もちろん、ダムや堰の管理を預かる方々は、こんな現象を生じないよう厳重に注意をしていますので、通常は心配する必要はありません。

　そのほかでは、狭い運河などで何かの拍子にできた孤立波が、かなり遠くまで伝わる現象が知られています。

　孤立波は、英語でｓｏｌｉｔｏｎ（ソリトン）という物質的な名前が与えられていますが、光ファイバーの中を伝わる赤外線のごく短いパルスもソリトンと呼べる場合があるそうです。
　光と水面の波とでは細かい性質は色々違いますが、波動の形式を持つエネルギーの固まりが遠くまで伝わる点で、ソリトンとしての共通性があると言えそうです。

＜＜５．端部での反射波について＞＞

　最後に、流れの最下端で波動エネルギーを綺麗に全部消滅させられない関係で、若干の反射波が上流に戻っていく様子まで模擬されています。　これはご愛嬌です。
　現実には、最下端に何か障害物が無い限りこんな反射波は起きないと思います。

　
　　＜＜レイトレーシングによる３次元グラフィック＞＞

　２０年ほど前の月刊アスキーに掲載され、下記のとおり単行本にもなったプログラム（Ｎ８８ＢＡＳＩＣ版）のうち、特にレイトレーシングの部分に惚れ込み、ＱｕｉｃｋＢａｓｉｃ時代も引き続き楽しんで来ましたが、最新のＶｉｓｕａｌＢａｓｉｃ文法の習得を兼ねて、移植と若干の機能追加をしたので実行形式プログラムを紹介します。

［２００７年５月上旬　追加］
　半ば本業に関連する「ダム湖の荒廃湖岸に緑を～」という論文に連動した３Ｄ画像を作成する機会がごく最近あり、自己の将来活動への種まきとしてＰＲします。
　追加配布するのは、RayTracing5.exeと関連データファイルで、RayTracing3との相違は次のとおりです。
　プログラム：
　　１）図形データ数の増大に対応したデータ配列サイズの倍増
　　２）水面下など透明体内部の透過において距離減衰（距離6.0ごと半減）を追加
　　３）上記距離減衰分の補填として、その４割を透明体自身の色で補う処理を追加
　　４）透過や反射による枝分かれ判定処理における打切り光強度を精密化
　　５）画像表示枠サイズが宣言値より縦・横４ドットずつ削られていた分を補正
　データファイル：　（２２種から２３種に増加）
　　１）巨大プランター（PLANTER.txt）を１番目に追加し、その他を繰り下げ
　　２）正１２面体の２種のデータにおいて、一部の自発光強度が強過ぎたのを修正

　結果的に、巨大プランターでは６０余りの立体図形を扱うことになり、処理時間が他のデータよりも長く掛かりますので、その点はご理解を願います。
　また全般的に透明体を透過した色は以前よりも若干ですが、くすみがちとなります。
　透過での距離減衰を追加した目的は、水深による見え方の差を現実に近づけるためであり、一定の効果が得られていることは右の画像でも判断頂けるかと思います。

[追加終わり]

　　プログラム名：　RayTracing3.exe　又は　RayTracing5.exe プログラム倉庫へ

（注）「Ｍｉｃｒｏｓｏｆｔ　Ｖｉｓｕａｌ　Ｂａｓｉｃ　２００５」による実行形式での配布のため、「．ＮＥＴ　ＦＲＡＭＥＷＯＲＫ」が機能している環境以外では実行できません。

　データファイルを多数含むＺＩＰ形式の圧縮ファイルですので、全て同じフォルダに復元してからご使用下さい。

　図形入力が簡略化されていなかったり、ピクチャーマッピング機能も無いなど、最新のレイトレーシング技法とは比べようもありませんが、基本的な機能の把握には十分に役立つと思いますし、何よりも手作り的な懐かしさがここにはあります。

◇書籍名：PC‐9801 3次元グラフィックス入門アスキーブックス
　　　　　　　　著者　守川穣（元　福井高等専門学校　教官）
　　　　Ａｍａｚｏｎなどに掲載されており、入手は可能と思われます。
　　　　図形データや計算方法などの詳細解説は、本書を利用下さい。

　なおプログラムの移植や機能追加などの改変については、著者の守川先生からご了承を頂いています。

　ソースプログラムは、初めてのＶｉｓｕａｌＢａｓｉｃ挑戦で、改変部分はお見せできるような代物でないですし、アスキー社さんへの配慮もあって掲載しない方針です。
　ただ、どんなめちゃくちゃソースでも気にしないという方は、ご相談下さい。

＜＜改変の概要＞＞
　以下にオリジナルプログラムからの変更点を記載しておきます。

１）Ｎ８８ＢＡＳＩＣから（ＱｕｉｃｋＢａｓｉｃ経由で）ＶｉｓｕａｌＢａｓｉｃへ
　言語仕様に合わせた変更。（型宣言、ＩＦなどの構文、関数、描画コマンドなど）

２）イベント駆動型の新しい大枠構成を造り、ボタンクリック等のイベント検出により
　元のプログラム全体をグラフィック画面の再描画サブルーチンとして呼び出す形式に
　組み替え。

３）最初のクラス宣言のところで、大部分の変数を前もって定義しておく形式とした。
　（旧Ｂａｓｉｃでは、宣言なしでサブルーチン枠を越えていたことへの補償。）

４）フルカラー化の試験を兼ねてグラディエーションのあるタイトル画面を追加。

５）初期のＰＣ－９８０１が８色どまりだったため、縦２×横２計４ドットで１画素を
　表していたが、フルカラー化に合わせて１画素１ドットに精細化を図った。

６）プログラム内のＤＡＴＡ文や、それを読むＲＥＡＤ文が使えなくなったことに対応
　して、立体図形データを複数のテキストデータファイルに独立させ、リストボックス
　で選択されたデータファイルを対象として処理する形式に変更。
　（実行形式プログラムのユーザーは、自己作成の図形データを読ませたい場合には、
　　データファイルのどれかを名前を変えずに書き換える方法が採用可能。）

７）視点を変える機能として、視野参照点を通る垂直軸回りの回転と、垂直方向の移動
　をテキストボックスに入力したうえで再実行できるように追加した。

＜＜移植での苦労話＞＞
　イベント駆動型ということで、グラフィック画面のリフレッシュサブルーチンを抜け出さないと画面が表示されず、旧来のように順次表示ドットが書き変わる課程が見られなくなったことが一番のダメージでした。
　その違和感に悩まされ、遠回りの苦労を強いられました。

＜＜データファイルについて＞＞
　データファイルで紹介する立体図形の組み合わせ中には、著者の守川先生がお造りになられたデータもありますが、それを大きく翻案したものや新たにデータを作ったものも加えています。
　とくに正１２面体は、根拠資料は既に無くしましたがデータを作るのに難儀した思い出深い立体です・・・。