数値流体力学はもう死んでる

日本の数値流体力学の第一人者だった桑原先生は亡くなられたらしい。

最近は、中国にも追いつかれてきたからね。

>>200
シミュレーション系の研究室に入るか
シミュレーション系研究室の日本人の学生と上辺だけでも友達になれば
コード集め放題だもんな

勉強したいんでコード見せて
の一言で簡単に手に入るw

優秀な奴はアメリカに留学しているだろｗ

>>202
本当に優秀だったら本国には帰らないぞ
中国が行け行けの時は本国に帰ってこない
優秀な奴を呼び戻すのに必死だったのを知らんのか？

逆だろ。軍事機密まで持ち出す中国系に、アメリカは盗まれないように必死だぞ。

差金決済取引

>>196
即レスだが、OpenFOAMのVOFでContactAngleを設定すると
それっぽい絵が出来るようだ

http://www.youtube.com/watch?v=5XZqxPJm5w4

シミュレーションて何が楽しいの?
自分も大学(学生)で数値シミュレーションの研究をやってるけど
所詮道具じゃん。
全くやりがいを感じないんだけど…

まったくだね。。でも大学の研究ならまだいいじゃん。

ベンダーに就職してシミュに関わったら人生終わりだよ。。
メーカーに転職出来ないので、この業界を転々とするしかない。。

>>207
シミュレーションは、同時に、シミュレーションと実測が一致することをやらないと
確かに面白くもなんともありません。
　普通は、たとえば、水槽ゆらして波作って、
シミュレーションと実測がだいたい合っている/合っている　なんてことをあるはず。

>>209
＞＞シミュレーションは、同時に、シミュレーションと実測が一致することをやらないと
確かに面白くもなんともありません。

実測も行っていますが、一致しようとしまいと私にとっては関心が持てません

あ

い

やー

にーはおぅ

このスレも死んでいる

その通り！！

>>210
自分は面白いけどな〜
DNSとかエミュレーターで遊んでる感覚だよ

自分にとっては
自作プログラム=ゲーム
論文=ゲームの説明書や攻略本
みたいなもので

テストプログラムとして
シミュレーション系の論文読んでパラメータ合わせて計算して結果比較
本チャンのゲームとして
実験系の論文読んでパラメータ化して自作プログラム動かして遊んでる

いいね。
遊びの感覚で、傍目無駄としか思えないパラメータ空間で処理を
走らせると・・・まあたいていは無駄なんだけどｗたまに面白いことが起こる。

取説つくる感覚、はとても正しい。
というか猿が気がつかない事に気付く・思いつく、
猿でも理解できるまで噛み砕いて説明する、のがリサーチャのお仕事だからね。
なんかしらんけど面白いでしょ？な「報告」だけでも意味はあるけどね
できればもう一歩踏み込むと、たのしい。
まあ反響は数年間で全世界で１０人くらい、な事が多いけどな・・・

いいね！そういうの。

>>12
乱流現象では渦層(vortex sheet)の組織構造や
渦線(vortex line)の時系列・空間的な変化が
DNSやLESで空間・時間的に一括でデータが手に入るようになって
詳細がわかってきた

この辺、つまり乱流の空間・時間的にまとまったデータなんて実験で得るのは難しい
熱線流速計じゃまずダメでLDVやPIVのこれからに期待というところなんだけど
乱流構造については、理論やらモデルが先行してた所にシミュレーションの結果が入ってきて
かなり影響を与えている

DNSは枯れた道具になってしまったけど
道具は道具なりに使い方次第って所だろうね

>>12
それと、層流乱流遷移なんかはDNSで一気に発展したね

・微小な撹乱のみしか扱えない線形安定性理論
・有限な振幅も扱えるが体系化されていない弱非線形性理論
に加えて1980年以降
線形から強非線形まで扱えるDNSが実験との差を埋める役割を果たした

例えば、平面クエット流や円管内流れなどは線形安定性では絶対安定なのだけれど
実験では乱流が観測される

この辺りの差異に対する回答の一つとして、DNSを使って
有限振幅の撹乱の非線形発達が層流乱流遷移を引き起こす可能性があることを
DNSのパイオニアたるOrszag等が示している

Orszag さんの名前を２chで見るとは！
最近 orz に似てなくも無いと思うようになってしまった。

10年前のに返答するって素敵w

朝っぱらからちゃちゃ入れる君も素敵w

朝っぱらから・・何この２人w

>>223
10年前のレスに返答したというより
>>12
のような考えを持っていてこのスレを読んでる人に向けたレス

乱流の研究ってこの先どういう方向に向かうのかイマイチ分からないんだよなぁ
俺が理解できないだけなのかもしれんけど

>>226
aliasing error はどの手法で除いてる？　またFFTは何を使ってるの？

>>228
aliasing error除去
単純だし、今はメモリが安いので

3/2 rule

他に
padding法
phase shift法
あえて取らない（擬スペクトル法のまま）
があるんだっけ？


FFT
自作よりも早いし使い勝手が良いので

fftw3

ただ、OpenMP非対応なので
AMD ACMLかIntel MKLのFFTに移行するかも

DNSってDirect Numerical Simulationだけど、これだと何の直接計算か分からないし
Direct Navier Stokesのほうが適切じゃない？
そもそも他の分野では前者の意味でDNSって言葉を使ったりはしないのかな？

aliasing error って（具体的には）なに？

ともあれ、英字３文字で何かの略称にする習慣はそろそろ止めてもらいたいのが本音だ。
わけわからんよw

>>231
>aliasing error って（具体的には）なに？

格子解像度が粗いと高周波の成分が低周波成分として捉えられてしまう問題
離散信号のフーリエ変換に付きまとう問題なので他分野でもお目にかかる

aliasing errorの影響について書いてる論文が少ないのだけど、除去しない場合
自分の経験上（点数が少ない場合）、二次不安定性から乱流への過程でエネルギーのカスケードが
上手くいかないようで乱流にならなかった
具体的に書くと平均速度分布が対数則に乗らない

理論上は解像度を十分に上げれば気にしなくてもいいはず

10年前のレスというか
10年間でこの分野が
進歩した

その分を、昔の人に報告してやったって感じかな

そういう処置をしてもDNSと言ってもいいもんなのかな
DNSの本来の意味って、陽的なモデルを使ってないってことなのかね？だったらいいんだけど…

>>232
さんくす。
Corseな格子系が短波長成分を取りこぼしてしまう、その部分が計算モデル自体に保存則を
持たせているせいで完全に消えてなくなるわけではなく、他の形(低周波成分)に化けて現れてしまう、
という事かな。非線型計算でも対処が必要な共通の悩みだね。

>>229
ありがとうございます。私はある理由があって、位相シフトを使っています。
>>234
モデルというわけでなく誤差ですね。差分法ではほぼ間違いなくこの誤差が含まれていますが、
波数空間で計算するスペクトル法は除くことができるということです。

スペクトル法でのaliasing errorも中心差分でのwiggleも、
解像度の粗さが原因という点で本質的には同じだよね？

つーかこんなスレタイでこんな議論をしてるとはシュールというか・・・ｗ

やはり１０年のうちに進むべきものはちゃんと進んだ、ということか
もっとも多少古いスキームでも格子数ガンガンにとればモダンとされる
スキームに匹敵するからね。偉い先生が若い頃に作ったのを
走らせて、なんじゃこりゃ、結構いけてる！と本人びっくりという話も聞く。
メーカー様さま、という要素もある。
やっぱり一番じゃないとね・・・・・そのあたりは生まれついての
日本人じゃない上に理系じゃない政治家さんには理解できないだろうが。

擬似スペクトルって何？
普通のスペクトル法とは違うの？

>>237
なんか線形や非圧縮と圧縮で用語が違うし・・・コード自体がまるで違うから
自分は圧縮で作って使ってるけど、非圧縮のことはかなり知らないし、
そういう人は多いのでは無いだろうか、と自分の勉強不足を楽観してみるテスト

>>197
同意。他と比べて少しも良くない。
過剰宣伝だよ、あれ。

まさに「どうしてこうなった・・・」って感じだなｗ

LANLの96年のスライド．
http://www.lanl.gov/orgs/t/t3/secondlevel/history/viewgraphs.pdf

最後のページ(Great Lies of CFD Codes)マジ面白いわｗｗｗ
今でも通じる話ばかりｗｗｗ

CIPはまだマシよ．過剰宣伝はしているが，少なくともアルゴリズムはちゃんと公開されているし．
最近の本では(バグはちょっと入っているが)動くサンプルも入っているし．

そこらへんにはアルゴリズムの一部だけ述べて解析例ばっかり書いたのを
論文に出して(そもそもこんなの論文になるのか?しかも結構まともな学会)
市販ソフト作って商売までしている奴もいるぜ．

ビジネスとしてはアリなんだろうけどね・・・
学会がある程度リーダーシップ持って不届きな連中を排除するような仕組みを作らんと駄目だな

245>>
偉い奴が率先してそうやっているから更に問題．
どんだけ金に植えているんだよ，全く．

ヨナもそうおもう。

>>247
何だそれは？