ただし、OpenFOAMは専用スレかGoogleフォーラムで議論しましょう
役に立つリンクはから
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 計算結果が載っているのでコードを組んで計算結果を比較してみよう
禿乙です
突然前スレなくなっちゃったんでびっくりしてました
新スレでもどうぞよろしくお願いいたします
最終的に水道管をまげたときに内部を流れる固形物を含む液体の流れの
変化の解析を行いたいと考え、CFDに手を出したいのですが、
この分野の事前の知識が全くない状態で、どこから手をつけたらいいかわかりません
理論の習熟と、ソフトの使い方は最低限必要だと思うのですが、何を勉強し、何を使えばいいのかもわかりません。
どなたか、おおざっぱにどういう順番で勉強していけば１行目のことができるようになるか、教えていただけませんか？
固形物の変形とか、分裂は考える？
例えば、配管の中でウン◯チが崩壊しながら流れる混相流とか
固形物は球状粒等単純なもので変形はほとんどしない前提で考えています
そのままでは発散してしまいます．
ありがとうございます。なんとなく筋が見えてきたのでやってみます
流体モデルでは真空中に広がっていくような計算は出来ないってことですかね？
そういう問題を解きたいのです．
これがあまりメジャーじゃないのはなぜなのですか？
ttp://http://manta.blog.jp/ ttp://http://sisutore.blog.jp/archives/1000786211.html ㌧㌧㌧クスです！
V&Vの資料の考え方は、シミュレーションの結果を他の人に伝える時にも
超重要っぽいですね
ただ、お話はけっこうむずそうorz 具他例とかもあったらわかりやすいかも……
とりあえず素人な自分は、変な結果がでてきたらそれに気づけるよう、
ファーツィガーさんの本をぽちってみました
着いたら、気合入れて読んでみます！
>両方の視点で書くと
むーん……
この辺をよく理解してないと、恐ろしいことになっちゃうんだ
２あたりは見当も付かないですorz
ttp://http://manta.blog.jp/ ttp://http://sisutore.blog.jp/archives/1000786211.html ttp://http://livedeaichat.blog.fc2.com/?no=1 ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
可視化して表現するためにはどのようなプロセスを踏んだら宜しいか？
ttp://awabi.2ch.net/test/read.cgi/psycho/1393514224/5 ↑ ↑ ↑ ↑ ↑ ↑
ttp://http://2ch.net/ otonaをつけたらいけるよー
↓
ttp://http://otona2ch.net/ 他に別に極楽だろ
当然正反対に理想郷だろ
ありがとうございます
でもこれらは非等方圧力ではないですよね？
うちの業界だけの言葉なんですかね？
圧力テンソルの対角成分が同じものでないことです。
例えば二次元なら
∇・P = (∂x Px) ex + (∂y Py) ey
独立な圧力がP一つではなくてPxとPyなんですよ
天体のジェットなんかを考える時はジェットの吹き出る方向と、それに垂直な方向では温度が違うんですよ
なんてかいたりしますね
P = P diag(1, 1, 1)なわけです
面の方向ごとに圧力が違うってこと？
そういう世界もあるんだなあ。状態方程式も圧力の複数の成分が入ってるのかな
同じものを対象にしていないとはどういう意味でしょうか。
ttp://http://www.physics.wisc.edu/grads/courses/726-f07/files/Section_14_Eq_General.pdf みたいな感じですね。これは磁場が入ってますけど
collisionいれればいいんじゃない？
どこに行けばいい？
ひとまずもう少し勉強してみます。ちなみに液液です
fftwはドキュメントが不親切過ぎて．．．
そしたらみなさんFFTってMPIでどうやってるんですかね？
ライブラリ使わず？？
FFTW以外のライブラリを使う？？
配慮して工事する金儲けより公共の利益を重視するすばらしい会社です
どうしたらよいかご存知でしょうか．
なんだか面白いんでこのスレに張り付いてるんですが、
高卒程度の知識で楽しみながらやる勉強法ってありますか？
すんません
出口側を切り落として行列を作り直すとはどういうことでしょうか
内点をf1~fnとするとf0=fn, fn+1 = f1という境界です．行列でかくと
2 -1 -1
-1 2 -1
....
-1 2 -1
-1 2 -1
右上と左下の-1が周期境界に対応してます
\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}
です
これは拡散方程式ですよね。しかも陽解法で解いてます。
ちょっと計算してみましたが仮に拡散方程式を陰解法で解いても行列はランク落ちしないです。
単純に考えてみると、周期境界のポアソン方程式は定数分の不定があるからランク落ちすると思うんですが、拡散方程式は方物型だから、コーシーデータが必要になるので、定数分の不定もなくなりますよね
フーリエで解けるのはわかるんです。0モードを考えないから
僕は差分でときたいんですねー
足立区に住んでいるそうだ
ttp://http://inumenken.blog.jp/archives/6580395.html 擬似逆行列ですか！勉強してみます
ちょっとおもったんですが、平均値がゼロという条件を課せばランク落ちしないんですが、帯行列でなくなるために随分ときにくくなって却下しました
理解したー？
26才を超えてるのに精神病院がどうのこうのとか。
俺は小学生じゃないぞ。
これだけ年を取ってるのにお前はいったいなんなんだ？
友達じゃねーんだ。
いちいちその時にそうやって祈っててどうもすいませんでした。
祈ってないのになー。
あら、私()と同じ悩みを持つものが
ですが、僕がやって見た感じだと、sorのような収束が遅いスキームだと収束するんですが、CG法みたいに速いスキームのときだと収束しない場合が多いです
足立区に住んでいるそうだ。
ttp://http://inumenken.blog.jp/archives/7002197.html イケメンなら稼げるんだろうけど。
話が上手けりゃ稼げるかな。
誰かレポ頼む。
メンガでググると出てくる。
ここ数年伸びが遅くなってる気がする
ttp://iup.2ch-library.com/i/i0289053-1303372232.jpg 
ttp://iup.2ch-library.com/i/i0289054-1303372232.jpg WqqAaqqQqqqqqqqQz いつの日か参考出品さaaqaqaqqれてた奴
AaaqQaaaaqaaaqqAqqaqaaqaqaqqqqaaqqqqaqaqqqaqqaqaqqqAaqqaqaaahqdvqqqqqqqqw
ttp://https://www.youtube.com/watch?v=LGlvDb3h4WA&feature=youtu.be 民生品で一番進んでいるであろうゲームでさえも
マルチコア対応を堂々と謳ってるのって少ないし…
GPGPU…うーん…。。。
需要の問題か…？飽和気味なのかもね。
性能伸びない気がするんだけどなぁ。。。
OpenMPはスレッド数増やすと性能伸びないし。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
マジレスを手伝った事にしたり首謀者のゴミにゴミに出された八つ当たりをして
無差別殺人か学会で変死
首謀者のゴミのいじめを手伝って無差別殺人か学会で変死
水槽の水を替える所から無かった事にして無差別殺人か学会で変死
助けないと解らない事して無差別殺人か学会で変死
ttp://http://nonke.blog.jp/ 載ってる
複数人で同じコードを開発してれば、そいうのは防げるし、下手くそなコード書いてる人の指導とかもやりやすい
皆さんはコード開発は何人くらいで同じコードを書いてますか？
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
どれ使ってやってるの？
詳しい情報ありがとう。
OpenFOAMじゃなく、Channelflowで遊ぶ理由はなに？
なにか優位性あるの？
いまのところ使う予定は未定。よっぽど暇なときにちょっと遊ぶかもって程度かな。
計算対象もいまのところ興味もってるもので無いし。
うちの場合みんなテーマがバラバラで、コード作りは各人任せなのよね。
それでほかの人からはノーチェックだから、仮になにか悪さしようと思ったらできる
まあ実験でも同じようにノーチェックなとこもあるのかもしれんけども
詳しい回答ありがとう。
ワンチップマイコンだとLEDの点滅。
流体力学シミュレーションの最初の例として一般的なのはどんなの？
教えてくれてありがとう。
こんど試してみるときに、それのどれかをやってみることにします。
一次元拡散方程式(熱伝導方程式)を FTCSで解く、は如何でしょうか。
解けたあとはクーラン数の勉強もできます
ありがとうございます。
熱伝導方程式をとくんですか？
流体関係のサンプルじゃなくて？
まだソフトを調べてないんだけど、微分方程式とか偏微分方程式を
数値的に解くソフトであって流体力学専用ソフトじゃないの？
詳しい解説ありがとう。
流体力学受講しなかったんだよな。
受講してても使わないから忘れてたろうけど。
偏微分方程式を解くプログラムで
流体力学で解く問題に関係しそうなものが得意なようにチューニングされてるソフトって解釈したけど間違ってないかな？
前色々一般的な偏微分方程式を解くソフトを探してたときに
ＧＳＬライブラリかなにかに、
偏微分方程式を解くときはケースバイケースなのでGSLには含むことを考えなかった的な文書を読んだ記憶がある。
主旨に反するけど、自作しまくってたら
作ってる人がよっぽど作るの上手くないとバグだらけってことになりませんか？
作ったものを,githubとかにあげてもらって、多くの人が使いやすいようにして
バグを利用者にテストしてもらって発見しやすくした方が良いと思う。
バグ報告管理機能がgithubについてたかどうか知らないけど、それがあるところで公開したほうが。
それすると使い方とかを文書にする手間が発生するけど。
通常プログラム作る人は作って計算するのが主目的で
使い方は作った本人が良くわかってるから文書書くのは苦痛だと思いますが。
の間違いです
√kx^2+ky^2である方向の波数の大きさわかるので。
なるほど！！！
ありがとうございます！
すみません！意味がわかりません！！
ちょいまち！これはとはちがいますよね？のやり方なら左右上下のカドは0埋め領域になるのでは
これはむしろで言うところの和集合になってるとおもいますが！
乱流になりきれないとは具体的にどのようになるの？
勝手に等方性を仮定してしまって、向きによってモード数が異なることに気持ち悪さを感じたので、ああいう書き方をしてしまいました。
Nx,Ｎyと書いてくれていたので等方である保証一切ないのに><
質問本文はＮでした。
えっと2/3法を満たすという意味では、DNS厨さんの絵がそのものずばりで、そっから更に各方向のモード数が等しくなるように切り捨てると183になります。エイリアジングエラー除去という意味では過剰です。等方性を神様のように崇める人達以外は気にしなくて良いと思います。
ヽo,´-'─ 、 ♪
r, "~~~~"ヽ
i. ,'ﾉﾚノﾚ!ﾚ〉 ☆ 日本のカクブソウは絶対に必須です ☆
__ '!从.ﾟ ヮﾟﾉﾙ 総務省の『憲法改正国民投票法』のURLです。
ゝﾝ〈(つY_i(つ
ttp://http://www.soumu.go.jp/senkyo/kokumin_touhyou/index.html ｀,.く,§_,_,ゝ,
~i_ﾝｲﾉ
博士の図の様にやります！
2D非圧縮NavierStokesを差分で解く時は非線形項をArakawa Jacobianで計算するのが一般的だと思うのですが、スペクトル法のときは気にせずそのまま計算していいのでしょうか
ところで、-N/2, N/2 の領域になっているのは何か意味があるの？
ガッカリだよ
流体構造の連成解析などでは、弾性体の固有振動数が鍵になったりします
どうして固有振動数を調べる必要があるんでしょう？(´･ω･`)
RSA Numbers 素因数分解に挑む・解析数式の考察
ttp://blog.livedoor.jp/superprojectx/archives/1009232905.html 円環の理ですね、わかりました(`･ω･´)ﾉｼ
原子遊び
原子核抜き
水素がきめて 2億年後
速い2000km/hタコメーター
車 電車など
水素崩壊
★★『幸福の科学』統合スレッドpart761★★
ttp://wc2014.2ch.net/test/read.cgi/psy/1410957061/ Q&Aより
なんで成り立つのか分かんないっす
今回の台風みたいに九州で急に東に進みだすとかいう挙動も，
シミュレーションで再現できるんでしょうか・・・＞＜
意思決定で使えないというのが原子力安全上の判断になったよね。
機械設計でシミュレーションは応用可能だろうが、危機管理の意思決定において少しでも不確実性が残る。
シミュレーション結果(ここではSPEEDI)を今までは判断基準としていたが、事故後の現在は参考資料に留めたよね。
これって結局、高レベル放射線は実測する必要があり事故の様子も確認しないと危機管理の意思決定には使えない
ということでしょ。
責任問題がいざ発生したとき、シミュレーション結果(確率的な誤差や曖昧さ含む)だけで避難指示と救援行動することは
危険だという事でしょ。
人が死んだり、天文学的な経済損失が発生して、一回きりで後戻りできない状況において
果たして理論上のシミュレーション結果だけで意思決定できるものなのかどうか。
条件がはっきりしてて仕組みもはっきりしてたらシミュレーションで十分って話になると思うよ。
例えばエミュレーターとPIC,PICの実機と比較すると
エミュレーターで遊んだ方がお金もかからないしほぼ一緒。
異なるのはチャタリングとか過電流がながれないとか。
原子力はもともと地震大国の日本で作るべきじゃなかったし、
継続決定は馬鹿すぎる選択だと考えてる。
頭悪すぎる非常に重要な選択を最悪の選択対応するのが得意。
そのつけはすぐにはっきりすることになるよ。
思いこんでるお前の方が馬鹿だよ。
薬ですら人体実験してから投与なのに。
俺らで今実験中って認識すらないだろ馬鹿すぎて。
どうなるかは先になってみないとわからない。
多分アウト。大量に漏れたら安全になるなんてことになるはずないだろ。
ttp://http://www.genanshin.jp/report/column/wave_20141029.html 時間進行方法や圧力の解き方がわかりません。
ご存知でしょうか。
アルゴリズムはfractional step法を使って、圧力はsor法で解いて見ようと
思います。今回の計算対象は中心の流れがほとんど無いのでスリップ条件にしよう
と考えています。
貴重なご意見をありがとうございました。
ttp://http://monoist.atmarkit.co.jp/mn/articles/1410/31/news023.html スパコン性能、2020年以降は伸びずに頭打ち!? ―その時CFDでは何ができるだろうか。東大生産研講演
NASAやBoeingは2030年に1ZFLOPS(1000Exa)級運用を目標にしてるらしいが、果たして現状を理解しているのか。
まだ性能向上が見込めるGPUクラスタでぶん回すのかねえ？w 死ぬ程効率悪いだろ。
ttp://http://www.youtube.com/watch?v=HJKHJEJgxgo ユーザーの方もユーザーでない方もぜひご参加お待ちしております。
ここの皆さんは参加されますか？
ドミーのトイレで田舎物とかやってこなくていいの。
違う違う違う違う。お前断られとるの。
原子力研究開発機構 原発推進 国策法人
原子力研究開発機構 原発推進 国策法人
税金でボーナス２３０万 原発推進 東電天下り集団
税金でボーナス２３０万 原発活用 東電天下り集団
原子力研究開発機構 東電天下り 国税２３００億円じゃびじゃびゅ
原子力研究開発機構 東電天下り 国税２３００億円じゃｂじゃびゅ
ttp://http://www.isas.jaxa.jp/j/topics/topics/2008/0911.shtml こういう流体の動きの模擬って、どの程度完璧なの？
それと、どうやったら、個人で同じようなこと出来ますか？
うわー、ありえねー
非構造格子（四面体）のセル中心で得られた流速を補間して節点での流速を求めたいのですが，
どのような方法がありますか？手法の名前や解説サイトでも助かります．
境界面の節点で滑りなしの速度境界条件を満たし，補間した値が二階微分可能であればなおうれしいです．
ttp://http://blngs.com/ 一回使ってみると、意外と使えるかと。
ありがとうございました．
5章の最後がそれっぽいですね．ちゃんと読んでみようと思います．
以下の論文では、位相シフト法をつかっていると記載していますが、
図1(b) では最大波数で大きな跳ねあがりが見られます。
これって普通なんですか？ 位相シフト法の問題ですか？
ttp://https://www.jstage.jst.go.jp/article/kikaib/79/803/79_1281/_pdf 中間波数領域と違って散逸以外に逃がす先が無いから。
最高波数への供給が粘性散逸より大きいという話だと思います。
差分とかで起きなかったとしたら差分とったときに高波数を潰すような誤差を作り出すのだと思われる
ttp://http://jbbs.shitaraba.net/sports/42269/ ☆ 自民党、グッジョブですわ。 ☆
ttp://http://www.soumu.go.jp/senkyo/kokumin_touhyou/index.html ☆ 日本国民の皆様方、2016年7月の『第24回 参議院選挙』で、改憲の参議院議員が
3分の2以上を超えると日本国憲法の改正です。皆様方、必ず投票に自ら足を運んでください。
そして、私たちの日本国憲法を絶対に改正しましょう。☆
速度場はきれいに落ちているので解像度の影響ということですね。
現在、1次元の位相シフト法ではうまくいくのですが、2次元の場合、高波数が跳ねてしまうので、
スキームの問題かと思いました。なにかミスがあるのだと思います。
ありがとうございました。
流れ方向の平均流速を出したところ，粘性低層は合うのですが，
対数則領域が過小に評価されてしまいます．
この原因は何だと考えられますでしょうか？(何らかの数値粘性で
チャネル中央部の流れが減衰してしまっているのでしょうか？？)
また，こういうことは，よく考えてプログラムを書かないと
しばしば起こるものなのでしょうか？？
なお，計算格子幅(壁座標)は壁方向が3で，流れ方向とスパン方向は30です．
また，対流項の離散化は2次の中心差分です．
どのようなアドバイスでもいただけたら大変ありがたいです．
よろしくお願いします．
流れ方向とスパン方向の解像度が粗すぎる。
平均量程度でも、dx+=18、dz+=9 ぐらいが必要。
乱流統計量まであわせたければ、さらにその倍は必要。
2次の中心差分で計算していれば数値粘性は入らない。
早速アドバイスをいただきありがとうございます．
格子数が足りな過ぎるのが原因だったのですね．
ご指摘の通りに細かくしてみます．
おっしゃるとおり中心差分には数値粘性が入っていませんので
私は相当とんちんかんなことを言っていますね．
有用なアドバイスを頂きありがとうございました．
今後LESもコーディングしたいと思いますが，その場合でも同程度の
解像度(dx+=18、dy+=3、dz+=9)が必要でしょうか？
森西先生の以下の論文が参考になるかと、
ttp://https://www.jstage.jst.go.jp/article/kikaib1979/66/647/66_647_1750/_pdf 解像度が粗いためにDNSでは散逸を捉えきれない影響を、SGS応力で補う観点からは、
論文中の図5に着目して、流量が減少する領域でLESの格子を設定する必要がある。
2次精度の場合はこの領域がかなり狭いということですね。
大変参考になる文献をご紹介いただきありがとうございました！
当方勉強不足なのでしっかりと勉強してみたいと思います．
図5によると，2次精度中心差分で精度を出すためには，
x+は20～40(z+は6～13)程度に選ぶ必要があるのですね．
当方の頭が悪くまだよく理解できていないのですが，
図5によると，格子が粗い場合の2次精度中心差分では
平均流速を過大に評価する傾向にあると思います．
当方の計算ではなぜ過小に評価されてしまったのでしょうか？？
図5において、2次精度の場合、dx+ < 40 ぐらいまでは、流量を過小評価している。
dx+ > 50 以上になると流量を過大評価している。
つまりあなたの計算結果はこれにあっていると言えると思います。
解像度を粗くして流量が減るのは、散逸が捉えられないため。
LESではGSで捉えられない、この散逸をSGSで補う。
従ってこの観点からは、格子分解能を流量が減少している領域で設定する必要がある。
無から有が生まれた
これが全ての真理
ttp://http://jbbs.shitaraba.net/sports/42269/ 時々天気予報が大幅に外れることがあるけれど、あれはシミュレーション結果を過信
しすぎたためと思われるが。
TVの天気解説にもGPU使ったような風向きベクトル示した短時間のシミュレーション動画が
流れてるが実際の大気の流れはあのように単純ではないよなといつも思うのだがｗ
春雨じゃ、濡れてまいろう。位の大雑把な気持ちでシミュレーションするのが病まなくて良い。
佐藤伸弦が暴行事件を起こしていた
ttp://i.imgur.com/V7P8RWS.jpg やっぱうっすいから写真とるとかっこいいなぁ
Sand management / Multiphase flows
ttp://http://www.lemma-ing.com/?page_id=976 
ttp://www.lemma-ing.com/wp-content/uploads/2012/09/emulsion_shear_rate_comparison.png 例えばこの画像のように、せん断速度をスカラーの色コンターで表示するのをよく見かけるのですが、
自分の理解では、せん断速度は変形速度テンソル（行列）の要素だと理解していて、
どうして色コンターで表示できるのか分かりません。
自分の理解は間違っているのでしょうか？
それとも、この画像のshear rateは何か特別な処理をしてスカラー量に変換されたものなのでしょうか？
絶対値じゃないですかね
もし3次元の場合、見る方向によって同じ位置でもせん断速度の
大きさが変わるということですか？
そうでした。ありがとうございました。当たり前でしたね。
セルの頂点に座標と全ての物理量（速度，圧力など）が定義されているタイプの離散化と，
セルの頂点にその座標を，また有限体積法のようにセルの中心で物理量が
局所値として定義されているタイプの離散化があることに気づきました。
（有限体積法では体積平均として定義されると思います。）
これらにはそれぞれメリット・デメリットなどあるのでしょうか？
後者は前者に比べ差分の計算の際に計算量が多くなってしまうように思いますが，
その分メリットなどあるのでしょうか？
全ての物理量が同じ位置に置かれるのではなくて、格子点上に圧力やその他のスカラー、dx/2ずらした位置にu、dy/2ずらしてv、dz/2ずらしてwみたいになるはずだけど。
メリットは解が振動することに対する対策。
スタガード格子のように圧力と速度の定義位置をずらすのではなく，
すべての量がセル頂点にあるのか，セル中心にあるのかという違いのことです。
境界条件の与え方が違ってくると思います。
（たとえば壁面上にも格子があるのか，半セルずれた位置で与えるのか）
簡単なコードでどう違うか確認してみようにも，ゼロから作るのはちょっと大変かな・・・
お返事ありがとうございます。
で挙げされたのは非圧縮性ソルバの話，
では非構造圧縮性ソルバ？の話で興味深く拝見しました。
非圧縮の場合はコロケートにしてかつ内挿・差分の扱いに気をつける必要があるみたいですね。
圧縮性非構造のほうでは検査体積の形や物理量再構築の際用いる量の取り方で精度が変わってくるということでしょうか。
まだ理解しきれていないので，じっくり読んでみたいと思います。
ちなみに私は構造格子圧縮性流のコードをコロケート風に書いています。
違う離散化で比較していないので，いまいちメリットがあるのか分からなくなり，質問させていただきました。
フリーの流体ソフトなんです。日本語の情報が無に等しいので日本では
まだ使ってる人いないのかな？
Windows64bit用をダウンロードしてみたけど使う以前に起動することすら苦労しそうです。
のまちがいでした。
翻訳(ほとんど機械翻訳ですが)もある程度進んだので、
これからちょっとずつサンプル見ながらやっていこうと思います。
Blenderはそこそこ使えますがLBMは使いこなすまではまだ。。
少し手のこんだことすると途端にPCが激重になるので
シンプルなもの以外はあまり試せていないという状態です。
物理学としての流体自体の知識はさっぱりないです。
ところでAAかわいいですね。どうやって描くのか感心します。
iOS版
ttp://http://itunes.apple.com/jp/app/9monsters/id451345005?mt=8# Android版
ttp://https://market.android.com/details?id=jp.co.applibros.alligatorxx ttp://http://ninemonsters.com フィルタ化したナビエストークス方程式の粘性項は
つぎのようにあらわせるとあります．
∂/∂x_j {2 * (nu + nu_t) * D_{ij}}
例えばスマゴリンスキーモデルでは
乱流粘性係数nu_tの項は空間的に変化するのでしょうか？
いい換えると
nu_t * ∂/∂x_j {2 * D_{ij}}
とnu_tを外に出してはいけないのでしょうか？
よろしくおねがいします．
LESのコード作りにためになる自作・フリーソフトをご存知でしたら教えてくれませんか．
OpenFOAMにはあまり良さを感じないので除外でお願いします．
ttp://i.imgur.com/dgw83fk.png 
ttp://i.imgur.com/KXDcsNd.png 
ttp://i.imgur.com/K9Y4HeD.png 
ttp://i.imgur.com/0wZLGwy.png 
ttp://i.imgur.com/yGz1U0r.png 
ttp://i.imgur.com/HIfm1uU.png 
ttp://i.imgur.com/p7Ux1yk.jpg 
ttp://i.imgur.com/avOjOHe.jpg 家猫の夜中の大運動会 地獄じゃなくて電波
ttp://https://twitter.com/oosakitakashi/status/621669244934041600 テレ東で中核派として登場した女性が
TBSに、一般市民としてインタビューを受けていたようですが
これは仕込みでしょうか？？？
TBSは中核派を一般市民として扱ってるのでしょうか？
ブルーム コスタ カレン デコラ 社長 中野隆一『なかの りゅういち』は
犯罪者である
お客様を赤羽警察署に偽証申告をして逮捕状をとる
共犯 デコラGM 遠藤孝輔『えんどうこうすけ』
ホットペッパービューティー
川口 池袋 評判 縮毛矯正 カット パーマ 格安 カーラー 美容院 床屋 散髪
美容室 ヘア メイク ネールサロン ザマスターキースリー
今夜は眠れない お楽しみ様です。
この日から地球の歴史上初めて真の世界平和が始まります。
人類は世界資源の分かち合いに同意するでしょう。
ttp://http://live.fc2.com/79752690/ 壁面で速度0、でお前の頭の中では、速度勾配du/dy も0なわけだな。
ってことは、壁で速度ゼロ、速度勾配もゼロ、壁垂直方向に積分して速度0となｗｗ
驚異の新事実 壁面流れでは流体は流れることができなかった！
死ねよ、バカ
じゃあ、なんで摩擦速度が0になるんだ。バカなのアホなのｗｗｗ
指摘されても間違いさえ気づけないのｗｗ 死ねよ、クズ。
壁関数が使える壁面流れで、壁垂直方向の途中から速度勾配が発生する要因を挙げてみろｗｗ
死ねよｗｗｗ
ただとで提示された条件だけでは数学的に導けていないと指摘しただけ
バカは惨めだねｗｗ 途中から速度勾配が発生する（スリップ）ように
お前の頭の中もスリップしてるんだろｗｗｗｗ
実在しないのに、数学的に導けてないだって。
それならなんだってありだわな。
ぼくのつくった境界条件ですかｗｗｗ
実際には存在しえない壁面摩擦速度0を仮定してる時点で何でもありだろ
ttp://nhk.2ch.net/test/read.cgi/livenhk/1441290593/ 
ttp://i.imgur.com/7aptG1x.jpg 
ttp://i.imgur.com/jBI98ec.jpg 
ttp://i.imgur.com/8q7xsUg.jpg 
ttp://i.imgur.com/KTkwbu8.jpg ttp://i.imgur.com/dgw83fk.png 
ttp://i.imgur.com/ewKXFH6.jpg 空を自由に飛びたいな
>&ft;1 元のスレ
../sim/dat/1394190575.dat
P35 自作・フリーソフト流体シミュレーションスレPart3
311692 -> 53431 (バイト)
P35 テストに夢中。これ以降書いても消えちゃうかも
人間たるもの誰しも，高速に，素晴らしい解析結果を得たいと考えるものである。
しかしそのようなケチな心のせいで，時に大きな損をすることになる。
一時の焦りが粗い格子でのGoサインを出してしまうのである。
格子解像度が足りず境界層が出てこなかった，渦が消えてしまった，
格子生成を妥協した結果として衝撃波が異常に厚く，ギザついてしまった。
このような場合には計算結果は，全く使いものにならない。
一見つまらない解析にもふんだんに格子点数を使える余裕を持った心持ちを普段から維持していきたいものである。
人間誰しも，よりはやく計算結果を得たいと考えてしまいがちである。しかし，焦り過ぎは禁物だ。
あなたの作った格子をよく見てみなさい。直交性がガタガタ，格子幅の拡大率は2倍近く，格子点数もかなりケチっちゃいましたね？
境界層速度分布はおかしくなり，渦は消え，衝撃波は異様に厚くなることでしょう。
このようなことにならないよう，リッチに格子点を使える余裕とおおらかな心を持ちたいものである。
ほかのサイトをクリックしても メモ帳を表示画面に出しっぱなしにできる固定できる
フリーソフトとかないでしょうか
運営のP35爆撃のせいで色々へこんだのでスレ建て直すわ
そっちでも運営がやりたい放題ならおーぷんに亡命するw
上田泰己氏 理研の税金無駄使い、954万円高級家具カッシーナ・イクスシーの指定購入も大問題
ttp://http://1000nichi.blog73.fc2.com/blog-entry-7696.html ttp://https://www.youtube.com/watch?v=quIHgwuF6r4 ttp://nhk2.2ch.net/test/read.cgi/livenhk/1478337570/ ttp://https://www.youtube.com/watch?v=_SOMgpdUc98 ttp://http://d.hatena.ne.jp/kouhou999/ ttp://https://youtu.be/7EekMD3GGHQ ttp://https://youtu.be/uH-WOOcNZ0s ttp://https://youtu.be/xHx5MbIGEoY 上田泰己氏 理研の税金無駄使い、954万円高級家具カッシーナ・イクスシーの指定購入も大問題
ttp://http://1000nichi.blog73.fc2.com/blog-entry-7696.html ttp://https://youtu.be/7EekMD3GGHQ ttp://https://youtu.be/uH-WOOcNZ0s ttp://https://youtu.be/xHx5MbIGEoY ttp://https://youtu.be/7EekMD3GGHQ ttp://https://youtu.be/uH-WOOcNZ0s ttp://https://youtu.be/xHx5MbIGEoY ttp://http://www.nikkansports.com/general/nikkan/news/1750500.html ボードゲームのオリジナルオーダー制作
ttp://http://www.logygames.com/logy/ordermade.html カードゲームを自作する1 【自宅でカード印刷】
ttp://http://tanishi.org/?p=801 500円ボードゲームのすすめ
ttp://http://kazuma.yaekumo.com/c0011_20161206advent.html 100円ショップでボードゲームを自作しよう
ttp://https://sites.google.com/site/jun1sboardgames/blog/makeyourbg ノーアイデアでボードゲームを作ろう第1回「100円ショップで物を買う」
ttp://http://boardgamelove.com/archives/boardgame-make-1/ ボードゲーム市場がクラウドファンディングの出現で急成長を遂げ市場規模を拡大中
ttp://http://gigazine.net/news/20150820-board-game-crowdfunding/ アナログゲームは実際アツい！ 識者が語る国内アナログゲーム市場の現況とゲームメカニクス
ttp://http://www.famitsu.com/news/201608/24114021.html 実際のところ、自作ボードゲームってどれぐらい売れるもんなの？
ttp://http://roy.hatenablog.com/entry/2016/11/04/225009 カフェも急増 ボードゲームにアラサーがハマる理由
ttp://http://style.nikkei.com/article/DGXMZO10921930R21C16A2000000?channel=DF260120166491 ボードゲーム印刷製造リンク集
ttp://http://www.tgiw.info/2014/09/link-production.html ttp://https://youtu.be/quIHgwuF6r4 ttp://https://youtu.be/quIHgwuF6r4 ここで何か質問する事があるかもしれません
よろしくお願いします
先に閉店かー！安部マリーオ！
ttp://https://www.youtube.com/watch?v=dpizDazz9Fg ttp://https://www.youtube.com/watch?v=HpX_uJUX92U ttp://https://www.youtube.com/watch?v=XHV9lUNP6rA ttp://https://www.youtube.com/watch?v=4UPAWP730a0 デジモンテイマーズは神アニメだよ
デジモンテイマーズは話題作だよ
デジモンテイマーズは意欲作だよ
デジモンテイマーズは超大作だよ
デジモンテイマーズは良作だよ
デジモンテイマーズは佳作だよ
デジモンテイマーズは上作だよ
デジモンテイマーズは秀作だよ
デジモンテイマーズは傑作だよ
デジモンテイマーズは名作だよ
ttp://https://blogs.yahoo.co.jp/sqjqs742/15340852.html を奪うつもりだ。 ← 民主主義の崩壊
ttp://http://www.data-max.co.jp/280113_ymh_02/ ↑ マスコミは 9条しか報道しないが 自民案
の真の恐怖は21条など言論の自由を奪うこと
ttp://http://blog.goo.ne.jp/kimito39/e/ec37220f64a8e1d6ed732dd0ab95cbf0 ↑超危険な緊急事態条項で人権無視の内閣独裁に！
ttp://https://www.youtube.com/watch?v=h9x2n5CKhn8 上のビデオで 自民党は 国民に基本的人権
は必要ないと怖いことを平気で言う。
ttp://http://xn--nyqy26a13k.jp/archives/31687 ↑ 都民ファーストも安倍と同じく 憲法改正で 人権
無視の大日本帝国憲法に戻すつもりだから
絶対に投票してはだめだ。 民主主義が崩壊する
ttp://http://blog.goo.ne.jp/ngc2497/e/8899f65988fe0f35496934dc972e2489 ↑ ネトウヨ= 安倍サポーター工作員はネットで国民を騙す。
ttp://https://dot.asahi.com/aera/2016071100108.html?page=3 ttp://http://blog.goo.ne.jp/kimito39/e/c0dd73d58121b6446cf4165c96ebb674 ↑ 安倍自民を操るカルト右翼「日本会議」は国民主権否定。
国民投票や選挙では自民党、維新、小池新党に絶対に入れるな。
市場に必要なことを妄想してみた。6年間の来場者数推移（2016年4月時点調べ）
ttp://http://bodoge.hoobby.net/columns/00001 ボードゲーム市場がクラウドファンディングの出現で急成長を遂げ市場規模を拡大中
ttp://http://gigazine.net/news/20150820-board-game-crowdfunding/ 中世っぽいデザインの金属サイコロ＆ダイスカップ「Rhythm Metal Gaming Dice」
ttp://http://gigazine.net/news/20140207-rhythm-metal-dice/ ファンタジー世界っぽい15種類のデザインのコインセット「Fantasy Coins」
ttp://http://gigazine.net/news/20140415-fantasy-coins/ デザイン戦略やタイポグラフィの歴史などがゲームしつつ理解できるトランプ「The Design Deck」レビュー
ttp://http://gigazine.net/news/20161027-design-deck-review/ ボードゲームのオリジナルオーダー制作
ttp://http://www.logygames.com/logy/ordermade.html 500円ボードゲームのすすめ
ttp://http://kazuma.yaekumo.com/c0011_20161206advent.html はじめてのボードゲーム制作記
ttp://http://nrmgoraku.hateblo.jp/entry/2017/05/22/221258 ゼロからボドゲを作って販売して分かった、3つのこと
ttp://http://begin-boardgames.seesaa.net/category/23191867-1.html 靴箱でテーブルサッカーゲームの作り方
ttp://http://www.handful.jp/curation/4207 ゲームマーケット2017春お疲れ様でした！！来場者は1000人増の13000人！
ttp://http://boardgamenumber.tokyo/2017/05/16/post-416/ 2017年開催のボードゲームイベント一覧
ttp://http://nicobodo.com/archives/19131088.html げんｋになってしのうとおもったんだけどな
じゃまおおいな
それを使って何をするかというのが全く決まらない
アイデア下さい！
ttp://http://www.youtube.com/watch?v=AP4JfKlSdFg ttp://https://ja.m.wikipedia.org/wiki/%E7%B7%B4%E9%A6%AC%E5%8C%BA%E6%8A%95%E7%A5%A8%E6%89%80%E8%A5%B2%E6%92%83%E4%BA%8B%E4%BB%B6 並列化ができるようになりましたが、いざ計算を行うとOUTPUT_FORMATに設定しておいた
PARAVIEWが適用されていないようで、.vtkファイルが1種類も出力されていません
まったく原因がわからないのですが、どうしたらいいでしょうか？
osはwindowsです
ttp://http://www.youtube.com/watch?v=43NKlCb1mOk ttp://http://www.youtube.com/watch?v=VWqe1L-8hSY ttp://http://www.youtube.com/watch?v=RYvuuNKZkZA ttp://http://www.youtube.com/watch?v=AP4JfKlSdFg 
ttps://i.imgur.com/ZvCasWc.jpg 