工学解析って何？！　日本オンリーワン株式会社

キャビティフローでの騒音発生問題についての考察

図のようなキャビティと気体の流れが発生させる騒音について以下のように考察する。

問題を解析で解く事と実験的に気体粒子の軌跡を考えてみると想像ではあるが以下の図のようになると考えられる。

気体は、キャビティのエッジからバックステップフローによるウエイクを生じ、流線の巻き込み現象が起こり、一部はキャビティ内で旋回する。この時に生じる圧力変動は、線形音響理論から次の2つの周波数の式で表される値に近似される。

騒音もこの圧力変動によって生じる気体の混合から生じることが予想される。

よって周波数の応答は、圧力変動の周波数応答に準じるものと成るであろう。また微小な渦の生成・崩壊による高周波成分が観測されるであろう。

東京大学生産技術研究所で行われている新幹線車両つなぎ目で生じる騒音の解析について、このことを提案したい。

OpenCascade

OpenCascadeとSALOMEは、フランスのOpenCascade　Technology社のリリースするGNUオープンソースの構造・流体解析のツールである。オープンソースであり高機能であるため、最低限の投資でトップパフォーマンスをたたき出す。

OpenCascadeTechnology社は、多くのサポートベンダーとカスタマーを抱え、日本では矢崎総業やMITSUTOYOがサポートベンダーであり開発サポート体制が敷かれている。

商用ライセンスを含み、OpenCascadeとSALOMEは、全てオープンでフリーで誰でも入手可能であり5万ダウンロードの実績がある。

OpenCascadeは、３Dモデリングと可視化部分のPRE・POSTでC++のライブラリーに加え、多くのフリーソフトたとえばVTK、Qｔなどで構成されている。３Dモデリング機能は、IGES/STEPのデータ読み込み機能を備え、UG社のCADなどの対応している。

マルチフィジカル・マルチスケール・連成シミュレーションや統合されたシミュレーションを実現する。 またSALOMEは、OpenCascade上で動作するシミュレーション統合プラットフォームであり、IGES/STEPの形状データ読み込み後、テトラ・プリズム・ヘキサ・混合のメッシュ生成や物理データ貼り付け、境界条件の貼り付けを行い、解析計算と可視化までサポートしている。CADと解析をドラマティックに融合リンクするツールである。

OpenCascade Technology社カスタマーサポートマネージャSergey ZertchaninovにパーミッションをいただいたのでOpenCascade本体の画像を載せる。

Maxima

GNUのフリーソフトMaximaを紹介する。
Maximaのホームページ

Maximaは、MITのMacsymaSystemをWilliam　F.SchelterがCommonLispに実装したものです。

(%i3) solve(x*x*x+1=0,x);
(%o3) [x=-(sqrt(3)*%i-1)/2,x=(sqrt(3)*%i+1)/2,x=-1]

など簡単に数式を処理いたします。
(%i3)の入力で(%o3)を返します。


図形描画機能もあります。
(%i4)plot3d(x^4+y^4,[x,-1,1],[y,-1,1]);

使ってみてください。

OPENFOAM

OpenFoamは、ロンドン大学インペリアルカレッジ・ゴスマン研究室に在籍していたHenry Weller(現在英国：OpenCFD Ltd.代表)らによって開発が進められ2004年にGNUのオープンソースとしてリリースされている流体解析のツールである。

Javaで駆動されるPre部分とC++ライブラリーによる各種ソルバー(PDEソルバー)とユーザーアプリの部分から構成されており、解析結果の可視化部分は、別のオープンソースのParaviewで行う様にデータ変換のツールが用意されている。Javaで駆動されるPre部分は管理画面のケースブラウザとメッシャーツール（ブロックメッシャー）とからなる。

ソルバーは、PDEソルバーといい、非圧縮性・圧縮性、乱流（多数のLESモデル）、2層流、自由表面、燃焼問題を扱える。またユーテリティーには、メッシュ操作（並列計算用分割、ダイナミックメッシュ）、LESモデル用POST、渦度、揚力、抗力算出が可能であり、ジョブ管理からデバックレベル管理などの管理画面が用意され、並列処理（LAM/MPI）も可能である。

海外では勿論のことで、オープンソースのCFDツール・高機能であるため日本でも学会等が取り上げ始めている。日本機械学会・日本建築学会・日本音響学会などで、OpenFoamを用いて解析がなされ論文が出ている。

Wikki.co.ukのHrvoje JasakにパーミッションをいただいたのでUSNavyのLarge Eddy Simulationの解析例を載せる。

ローレンツ方程式(その２）

前回に、余計な項が含まれていたので、解きなおした。

>function xdot=f(x,t)
xdot(1)=-10*x(1)+10*x(2)
xdot(2)=-x(1)*x(3)+28*x(1)-x(2)
xdot(3)=x(1)*x(2)-8/3*x(3)
endfunction
>x0=[1;1;1];
>t=linspace(0,500,1000);
>x=lsode("f",x0.t);
>plot(t,x);

Oh,No!

アメリカの金星探査ロケットで起こったことである。それまで順調に飛行していたのが、軌道修正時に行方不明になった。何回か軌道修正して順調に飛行を続けていたが、金星目前で軌道修正した途端、予定されている軌道をはずれ、地球から所在が把握できなくなったのである。

巨額の投資をして開発した高価なロケットが行方不明となり、徹底的な調査が行われた。

原因が判明した。なんとたった一つのコンマがピリオドに化けていたことがその重大な事態を導いたことが分かる。

DO 100 I=1.10と言う記述が見つかりプログラムがDO100I=1.10と解釈して変数DO100Iに1.10を代入しただけとなっていることが判明した。

Oh,No!

ローレンツ方程式

流体のカオス的ふるまいを表す方程式としてローレンツ方程式がある。1963年Lorenz,E.N.が”Deterministic Nonperiodic Flow",Journal of Atmospheric Sciencesに発表した決定論的非周期流れの方程式である。

dx/dt=-px+py
dy/dt=-xz+rx-y
dz/dt=xy-bz

ローレンツがこの論文で与えたのは、
p=10
r=28
b=8/3
であり、この元に解かれるX,Y,Zの軌跡（アトラクタ）が示された。

この方程式を私も解いてみた。OctaveというGNUのソフトを使い解く。

>function xdot=f(x,t)
xdot=zeros(1,3)
xdot(1)=-10*x(1)+10*x(2)
xdot(2)=-x(1)*x(3)+28*x(1)-x(2)
xdot(3)=x(1)*x(2)-8/3*x(3)
endfunction
>x0=[1;1;1];
>t=linspace(0,500,1000);
>x=lsode("f",x0.t);
>plot(t,x);

以上の設定で解けるはずである。

アトラクタは示せないが、確かにカオス的解が得られる。このことが流体のふるまいを示していると言うのいである。アトラクタを示せばよりわかり易いであろうが、たとえば小川の中の石の周りの流れがこのことを表している。同じようなパターンの流れを示したかと思えば異なるパターンにぱっと変わり、また元に戻りこれが非周期的に起こる。川の中の石の周りの流れをじっと見て余生を過ごしたレオナルド・ダ・ビンチもこういったことが連立常微分方程式の初期値問題でこんなことが表せると気づいたらさぞかし喜んだであろう。

地球シミュレータとTSUBAME

2002年3月に登場したES（Earth Simulator)は、ベクトルプロセッサを5,120個並列につなぐことによって40TFLOPSの性能を示し当時コンピュートニックショックと言うほど世界を圧倒してきた。しかしESの登場によって触発されたアメリカは、巻き返しのため躍起で高性能のスパコンを開発し始める。今やBlueJeanという米国IBMのスパコンがESの10倍ほど高速で世界一の高速演算を誇る。

現在日本で一番高速演算するスパコンは、東工大のTSUBAMEであり、世界8位である。X86のCPUを並列につなげて高速演算させるのでX86対応のソフトなら何でも走る。また省電力化にも工夫がなされ空調に工夫がある。ダウンフローとサイドフローの空調で温度調節がうまくなされ、ESよりかなり省電力で運転できる。
　　　<

職人と日本が捨ててはならぬもの

大阪下町の職人。
東京大田区の町工場の職人。
に代表される人たちは、こう思うであろう。
”今俺が作ってやるから、見てろ！！！”

時には図面を見ずに、時には地べたに蝋で図面を書いて作り始める。私は、東京電力㈱や中部電力㈱の発電所サイトで登録され働く職人を見てきた。図面を完全に頭に描いて作っている。手間が早い。段取りに時間をかけ、段取り７割りというまで入念に仕事の準備をする。仕事始めから終わりまでをきっちりこなし養生して次の仕事につなげる。

挨拶は、”ご安全に！”である。
無精ひげを生やしている人もいたり、首に手ぬぐいを巻いている人もいて、地べたに座って昼食をとったりする。新日鉄㈱設計会社の東海設計㈱の工務部の職人は、首に手ぬぐいを巻いたりしてSUN端末のCADで図面を書いたりしている。

今や職人たちがCAD図面や官庁申請書を書いている時代であり、そうできないと工事できない時代である。

元々勘が良く、運動神経抜群なのであるから物覚えも速いし、物づくりを毎日やっているのであるからCADなんか簡単に操作する。出来ないわけがない。

職人は、”俺は学校が嫌いだった”といい、今では人生の中で職聖を見出し立派に毎日働いている。こういった人の中に超ベテランがボウシンとして存在して、勤めている。

この人たちの物づくりは、経験から来る理屈と自信に満ちており筋が通っている。時には理論解や解析解を凌駕するほど精度が高い。

日本は、第二次世界大戦の時、零式戦闘機や戦艦大和や１万メートル射程の大砲を作り出した脅威の物づくり大国である。戦争兵器だけでなく今日は職人技によって開発されたものに溢れている。ありとあらゆる建物は全部職人が作ったものである。

これらのことを職人に感謝して、日本が絶対に捨ててはならぬものとして継承してゆきたい。このような分野で解析や分析の仕事をしてゆこうと思う。