自作・フリーソフト流体シミュレーションスレPart3 P35
CFD全般、圧縮・非圧縮、単相流・多相流何でも議論しましょう
ただし、OpenFOAMは専用スレかGoogleフォーラムで議論しましょう
役に立つリンクは>>1から フリーソフト
DNS:Channelflow
ttp://channelflow.org/
VOF
ttp://www.lmm.jussieu.fr/~zaleski/codes/
OpenFVM
ttp://sourceforge.jp/projects/sfnet_openfvm/releases/
いろいろ:OpenFOAM
ttp://www.openfoam.com/ 初心者および勉強用リンク
地球惑星物理学演習:時間発展問題の数値解法
ttp://www-solid.eps.s.u-tokyo.ac.jp/~ataru/edu/ensyu08.pdf
物理モデルの基礎と選択法 乱流モデル
ttp://accc.riken.jp/secure/4650/080121_yamamoto.pdf
連立1次方程式の基礎知識 反復解法の勉強
ttp://yebisu.cc.kyushu-u.ac.jp/~watanabe/RESERCH/MANUSCRIPT/TUTORIAL/leq.pdf
計算科学チュートリアル2008 流体シミュレーションのベンチマーク
ttp://www2.cc.kyushu-u.ac.jp/scp/users/forum/tutorial20080501-02/3-3-RI2T_tutorial080502_matsumoto_2-2.pdf
ベンチマークの一つ 二次元キャビティの論文
ttp://www.ann.jussieu.fr/~frey/papers/Navier-Stokes/Ghia%20V.,%20High-Resolutions%20for%20incompressible%20flow%20using%20the%20Navier-Stokes%20equations%20and%20a%20multi-grid%20method.pdf
計算結果が載っているのでコードを組んで計算結果を比較してみよう >>1-5
禿乙です
突然前スレなくなっちゃったんでびっくりしてました
新スレでもどうぞよろしくお願いいたします 質問です
最終的に水道管をまげたときに内部を流れる固形物を含む液体の流れの
変化の解析を行いたいと考え、CFDに手を出したいのですが、
この分野の事前の知識が全くない状態で、どこから手をつけたらいいかわかりません
理論の習熟と、ソフトの使い方は最低限必要だと思うのですが、何を勉強し、何を使えばいいのかもわかりません。
どなたか、おおざっぱにどういう順番で勉強していけば1行目のことができるようになるか、教えていただけませんか? >>5
固形物の変形とか、分裂は考える?
例えば、配管の中でウン◯チが崩壊しながら流れる混相流とか >>6
固形物は球状粒等単純なもので変形はほとんどしない前提で考えています 圧縮性流体コードをつかって衝撃波菅問題を解いているのですが,片側を真空(n=0,p=0)にしたいのですが,どうしたらよいでしょうか.
そのままでは発散してしまいます. >>1
>>18
ありがとうございます。なんとなく筋が見えてきたのでやってみます >>1
流体モデルでは真空中に広がっていくような計算は出来ないってことですかね?
そういう問題を解きたいのです. 最近擬似圧縮性法って知ったのですが
これがあまりメジャーじゃないのはなぜなのですか? >>31-32
dddクスです!
V&Vの資料の考え方は、シミュレーションの結果を他の人に伝える時にも
超重要っぽいですね
ただ、お話はけっこうむずそうorz 具他例とかもあったらわかりやすいかも……
とりあえず素人な自分は、変な結果がでてきたらそれに気づけるよう、
ファーツィガーさんの本をぽちってみました
着いたら、気合入れて読んでみます!
>両方の視点で書くと
むーん……
この辺をよく理解してないと、恐ろしいことになっちゃうんだ
2あたりは見当も付かないですorz ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。 ド素人が、実機流体の温度分布や速度分布、あるいは配管の熱応力分布を
可視化して表現するためにはどのようなプロセスを踏んだら宜しいか? 非等方圧力で圧縮性流体(roeやhll)といてる先行研究とか知ってます? 寧ろ逆に天国だろ
他に別に極楽だろ
当然正反対に理想郷だろ >>1
ありがとうございます
でもこれらは非等方圧力ではないですよね? >>1
うちの業界だけの言葉なんですかね?
圧力テンソルの対角成分が同じものでないことです。
例えば二次元なら
∇・P = (∂x Px) ex + (∂y Py) ey >>1
独立な圧力がP一つではなくてPxとPyなんですよ それもあるかもしれないですね。
天体のジェットなんかを考える時はジェットの吹き出る方向と、それに垂直な方向では温度が違うんですよ P = P⊥ diag(1, 1, 0) + P|| diag(0, 0, 1)
なんてかいたりしますね 通常の流体(等方圧力)では
P = P diag(1, 1, 1)なわけです 分子運動論的な言い方をすると、方向ごとに熱速度の大きさが違うので
面の方向ごとに圧力が違うってこと?
そういう世界もあるんだなあ。状態方程式も圧力の複数の成分が入ってるのかな 物理現象としては同じものを対象としてないって理解でOK? 初心者ですが、二流体モデルをやるにはどうすればよいんでしょうか? >>32
collisionいれればいいんじゃない? とりあえず、流体がシミュレーション通りに動くことを己の目で見たいのだけど
どこに行けばいい? >>1
>>69
ひとまずもう少し勉強してみます。ちなみに液液です fortranのmpi fftwのサンプルコードとか持ってる人いるかな.
fftwはドキュメントが不親切過ぎて... >>1
そしたらみなさんFFTってMPIでどうやってるんですかね?
ライブラリ使わず??
FFTW以外のライブラリを使う?? 前処理付き共役勾配法が今のところ最強ってことでいい? なるほどね!今は無難に不完全コレスキーだから、ほかのも試してみようかしら は日曜日以外毎日8時から18時過ぎまで長期間に渡って周囲に騒音をまき散らさないように 配慮して工事する金儲けより公共の利益を重視するすばらしい会社です 朝日住宅は日曜日以外毎日8時から18時過ぎまで長期間に渡って周囲に騒音をまき散らさないように
配慮して工事する金儲けより公共の利益を重視するすばらしい会社です Poisson方程式を周期境界条件で解きたいのですが,行列のrankが落ちてるので解けません.
どうしたらよいかご存知でしょうか. 残差切除法ってどうなんすかね?poisson方程式の解法として >>42 計算格子のどこか一点(四隅のどこかとか)を0とかに固定すればおk。 ハカセじゃないと言葉が流暢だなぁ
なんだか面白いんでこのスレに張り付いてるんですが、
高卒程度の知識で楽しみながらやる勉強法ってありますか? >>1
すんません
出口側を切り落として行列を作り直すとはどういうことでしょうか >>1
内点をf1~fnとするとf0=fn, fn+1 = f1という境界です.行列でかくと
2 -1 -1
-1 2 -1
....
-1 2 -1
-1 2 -1
右上と左下の-1が周期境界に対応してます 行列わかりづらいですね.texで書きますと,
\begin{array}{ccccc}
2 & -1 & 0 & 0 & -1 \\
-1 & 2 & -1 & 0 & 0 \\
0 & -1 & 2 & -1 & 0 \\
0 & 0 & -1 & 2 & -1 \\
-1 & 0 & 0 & -1 & 2 \\
\end{array}
です >>1
これは拡散方程式ですよね。しかも陽解法で解いてます。
ちょっと計算してみましたが仮に拡散方程式を陰解法で解いても行列はランク落ちしないです。
単純に考えてみると、周期境界のポアソン方程式は定数分の不定があるからランク落ちすると思うんですが、拡散方程式は方物型だから、コーシーデータが必要になるので、定数分の不定もなくなりますよね >>1
フーリエで解けるのはわかるんです。0モードを考えないから
僕は差分でときたいんですねー ひきこもり30代後半で童貞で低収入イラストレーターやってるバカ発見。
足立区に住んでいるそうだ
http://inumenken.blog.jp/archives/6580395.html >>1
擬似逆行列ですか!勉強してみます
ちょっとおもったんですが、平均値がゼロという条件を課せばランク落ちしないんですが、帯行列でなくなるために随分ときにくくなって却下しました あんなー目撃証言0とかでこういう事をするとなー東日本大震災とかになるだわ。
理解したー?
26才を超えてるのに精神病院がどうのこうのとか。
俺は小学生じゃないぞ。
これだけ年を取ってるのにお前はいったいなんなんだ?
友達じゃねーんだ。 あーはいそうですかー?
いちいちその時にそうやって祈っててどうもすいませんでした。
祈ってないのになー。 新聞配達員の最中に食べに行った刈谷市の辺のファミリーレストランを解らなくして東日本大震災 一点固定の方法でキャビティ流れをやろうとしましたが、うまくいきませんでした...。キャビティ流れでポアソン方程式を解く場合、行列のランク落ちをどうすればいいんでしょう? >>57
あら、私(>>88)と同じ悩みを持つものが >116
>>42ですが、僕がやって見た感じだと、sorのような収束が遅いスキームだと収束するんですが、CG法みたいに速いスキームのときだと収束しない場合が多いです >116 解けました! 速度の境界条件をミスしてました。 ネット版のホストみたいなの発見。
イケメンなら稼げるんだろうけど。
話が上手けりゃ稼げるかな。
誰かレポ頼む。
メンガでググると出てくる。 計算機の性能向上って頭打ちなのかなぁ…
ここ数年伸びが遅くなってる気がする qqqqqqaqaqaqaaqaqhttp://iup.2ch-library.com/i/i0289053-1303372232.jpg
http://iup.2ch-library.com/i/i0289054-1303372232.jpg
WqqAaqqQqqqqqqqQz いつの日か参考出品さaaqaqaqqれてた奴
AaaqQaaaaqaaaqqAqqaqaaqaqaqqqqaaqqqqaqaqqqaqqaqaqqqAaqqaqaaahqdvqqqqqqqqw AaSAasAAqaaqq(qqqqa(AAaaaqaaaqqqq(ww(wwew)))) >>1
民生品で一番進んでいるであろうゲームでさえも
マルチコア対応を堂々と謳ってるのって少ないし…
GPGPU…うーん…。。。
需要の問題か…?飽和気味なのかもね。 メニーコアCPU出ても、FDMとかFVMなら結局メモリバンドぶんしか
性能伸びない気がするんだけどなぁ。。。
OpenMPはスレッド数増やすと性能伸びないし。 辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。 辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。 ふっふっふ時代の流れにインプット
マジレスを手伝った事にしたり首謀者のゴミにゴミに出された八つ当たりをして
無差別殺人か学会で変死 ふっふっふ時代の流れにインプット
首謀者のゴミのいじめを手伝って無差別殺人か学会で変死 ふっふっふ時代の流れにインプット
水槽の水を替える所から無かった事にして無差別殺人か学会で変死 ふっふっふ時代の流れにインプット
助けないと解らない事して無差別殺人か学会で変死 STAP騒動を見てて疑問に思ったんだけど、シミュレーションコード開発ってやろうと思えば実験より簡単にデータの捏造ができるとおもうんですよね。
複数人で同じコードを開発してれば、そいうのは防げるし、下手くそなコード書いてる人の指導とかもやりやすい
皆さんはコード開発は何人くらいで同じコードを書いてますか? デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね >>2 のフリーソフトのどれがお勧め?
どれ使ってやってるの? >>148-151
詳しい情報ありがとう。
OpenFOAMじゃなく、Channelflowで遊ぶ理由はなに?
なにか優位性あるの?
いまのところ使う予定は未定。よっぽど暇なときにちょっと遊ぶかもって程度かな。
計算対象もいまのところ興味もってるもので無いし。 >>1
うちの場合みんなテーマがバラバラで、コード作りは各人任せなのよね。
それでほかの人からはノーチェックだから、仮になにか悪さしようと思ったらできる
まあ実験でも同じようにノーチェックなとこもあるのかもしれんけども プログラム言語の入門で最初によくでてくるのがHello World。
ワンチップマイコンだとLEDの点滅。
流体力学シミュレーションの最初の例として一般的なのはどんなの? >>1
教えてくれてありがとう。
こんど試してみるときに、それのどれかをやってみることにします。 >>83
一次元拡散方程式(熱伝導方程式)を FTCSで解く、は如何でしょうか。
解けたあとはクーラン数の勉強もできます >>85
ありがとうございます。
熱伝導方程式をとくんですか?
流体関係のサンプルじゃなくて?
まだソフトを調べてないんだけど、微分方程式とか偏微分方程式を
数値的に解くソフトであって流体力学専用ソフトじゃないの? >>1
詳しい解説ありがとう。
流体力学受講しなかったんだよな。
受講してても使わないから忘れてたろうけど。
偏微分方程式を解くプログラムで
流体力学で解く問題に関係しそうなものが得意なようにチューニングされてるソフトって解釈したけど間違ってないかな?
前色々一般的な偏微分方程式を解くソフトを探してたときに
GSLライブラリかなにかに、
偏微分方程式を解くときはケースバイケースなのでGSLには含むことを考えなかった的な文書を読んだ記憶がある。 自作ってスレについてるから
主旨に反するけど、自作しまくってたら
作ってる人がよっぽど作るの上手くないとバグだらけってことになりませんか?
作ったものを,githubとかにあげてもらって、多くの人が使いやすいようにして
バグを利用者にテストしてもらって発見しやすくした方が良いと思う。
バグ報告管理機能がgithubについてたかどうか知らないけど、それがあるところで公開したほうが。
それすると使い方とかを文書にする手間が発生するけど。
通常プログラム作る人は作って計算するのが主目的で
使い方は作った本人が良くわかってるから文書書くのは苦痛だと思いますが。 x方向の[-N/2, -2/3(N/2)]および[2/3(N/2), N/2]とy方向の[-N/2, -2/3(N/2)]および[2/3(N/2), N/2]の...
の間違いです 波数で円を書いて2/3法内に入っている波数を使えばいいと思う。
√kx^2+ky^2である方向の波数の大きさわかるので。 >>1
すみません!意味がわかりません!!
>>186
ちょいまち!これは>>183とはちがいますよね?>>183のやり方なら左右上下のカドは0埋め領域になるのでは
これはむしろ>>181で言うところの和集合になってるとおもいますが! >>1
乱流になりきれないとは具体的にどのようになるの? 混乱させてしまって申し訳ありません。
勝手に等方性を仮定してしまって、向きによってモード数が異なることに気持ち悪さを感じたので、ああいう書き方をしてしまいました。
Nx,Nyと書いてくれていたので等方である保証一切ないのに>< ん。
質問本文はNでした。
えっと2/3法を満たすという意味では、DNS厨さんの絵がそのものずばりで、そっから更に各方向のモード数が等しくなるように切り捨てると183になります。エイリアジングエラー除去という意味では過剰です。等方性を神様のように崇める人達以外は気にしなくて良いと思います。 ___ _
ヽo,´-'─ 、 ♪
r, "~~~~"ヽ
i. ,'ノレノレ!レ〉 ☆ 日本のカクブソウは絶対に必須です ☆
__ '!从.゚ ヮ゚ノル 総務省の『憲法改正国民投票法』のURLです。
ゝン〈(つY_i(つ http://www.soumu.go.jp/senkyo/kokumin_touhyou/index.html
`,.く,§_,_,ゝ,
~i_ンイノ 皆さんご親切にありがとうございました
博士の図の様にやります! 私もスペクトル法で質問があります
2D非圧縮NavierStokesを差分で解く時は非線形項をArakawa Jacobianで計算するのが一般的だと思うのですが、スペクトル法のときは気にせずそのまま計算していいのでしょうか