自作・フリーソフト流体シミュレーションスレPart3 P35

**DNS厨** · 2014/03/07(金) 20:09:35.07

CFD全般、圧縮・非圧縮、単相流・多相流何でも議論しましょう
ただし、OpenFOAMは専用スレかGoogleフォーラムで議論しましょう

役に立つリンクは>>1から

**DNS厨** · 2014/03/07(金) 20:15:17.43

フリーソフト
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/

**DNS厨** · 2014/03/07(金) 20:17:17.09

初心者および勉強用リンク
地球惑星物理学演習:時間発展問題の数値解法
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
計算結果が載っているのでコードを組んで計算結果を比較してみよう

**前スレ441** · 2014/03/08(土) 14:11:39.98

>>1-5
禿乙です
突然前スレなくなっちゃったんでびっくりしてました
新スレでもどうぞよろしくお願いいたします

**名無しさん＠５周年** · 2014/03/16(日) 01:39:04.79

質問です
最終的に水道管をまげたときに内部を流れる固形物を含む液体の流れの
変化の解析を行いたいと考え、CFDに手を出したいのですが、
この分野の事前の知識が全くない状態で、どこから手をつけたらいいかわかりません
理論の習熟と、ソフトの使い方は最低限必要だと思うのですが、何を勉強し、何を使えばいいのかもわかりません。
どなたか、おおざっぱにどういう順番で勉強していけば１行目のことができるようになるか、教えていただけませんか？

**名無しさん＠５周年** · 2014/03/18(火) 00:46:09.08

>>5
固形物の変形とか、分裂は考える？
例えば、配管の中でウン◯チが崩壊しながら流れる混相流とか

**名無しさん＠５周年** · 2014/03/18(火) 01:26:02.98

>>6
固形物は球状粒等単純なもので変形はほとんどしない前提で考えています

**名無しさん＠５周年** · 2014/03/18(火) 15:21:32.68

圧縮性流体コードをつかって衝撃波菅問題を解いているのですが，片側を真空(n=0,p=0)にしたいのですが，どうしたらよいでしょうか．
そのままでは発散してしまいます．

**名無しさん＠５周年** · 2014/03/18(火) 22:55:11.61

>>1
>>18
ありがとうございます。なんとなく筋が見えてきたのでやってみます

**名無しさん＠５周年** · 2014/03/19(水) 01:56:08.42

>>1
流体モデルでは真空中に広がっていくような計算は出来ないってことですかね？
そういう問題を解きたいのです．

**名無しさん＠５周年** · 2014/03/29(土) 16:12:06.25

落ちてたのか、気付かなかった。

**名無しさん＠５周年** · 2014/03/30(日) 05:00:56.18

最近擬似圧縮性法って知ったのですが
これがあまりメジャーじゃないのはなぜなのですか？

**名無しさん＠５周年** · 2014/04/01(火) 07:33:35.92

法的にいいのかよ、これらってさ。

http://manta.blog.jp/

http://sisutore.blog.jp/archives/1000786211.html

26 · 2014/04/02(水) 19:04:50.82

>>31-32
㌧㌧㌧クスです！

V&Vの資料の考え方は、シミュレーションの結果を他の人に伝える時にも
超重要っぽいですね
ただ、お話はけっこうむずそうorz　具他例とかもあったらわかりやすいかも……
とりあえず素人な自分は、変な結果がでてきたらそれに気づけるよう、
ファーツィガーさんの本をぽちってみました
着いたら、気合入れて読んでみます！

>両方の視点で書くと
むーん……
この辺をよく理解してないと、恐ろしいことになっちゃうんだ
２あたりは見当も付かないですorz

**名無しさん＠５周年** · 2014/04/02(水) 19:17:03.95

多額な借金で炎上中。

http://manta.blog.jp/

http://sisutore.blog.jp/archives/1000786211.html

http://livedeaichat.blog.fc2.com/?no=1

**名無しさん＠５周年** · 2014/04/02(水) 21:35:45.68

ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。

**名無しさん＠５周年** · 2014/04/08(火) 23:39:45.37

ド素人が、実機流体の温度分布や速度分布、あるいは配管の熱応力分布を
可視化して表現するためにはどのようなプロセスを踏んだら宜しいか？

**名無しさん＠５周年** · 2014/04/16(水) 04:35:29.15

http://awabi.2ch.net/test/read.cgi/psycho/1393514224/5
↑　　↑　　　　↑　　↑　　　　↑　　↑

**名無しさん＠５周年** · 2014/04/16(水) 17:44:33.69

http://2ch.net/
otonaをつけたらいけるよー
↓
http://otona2ch.net/

**名無しさん＠５周年** · 2014/05/02(金) 04:26:26.21

非等方圧力で圧縮性流体(roeやhll)といてる先行研究とか知ってます？

忍法帖【Lv=40,xxxPT】(1+0：8) · 2014/05/04(日) 23:21:25.21

寧ろ逆に天国だろ
他に別に極楽だろ
当然正反対に理想郷だろ

**名無しさん＠５周年** · 2014/05/09(金) 02:28:23.74

>>1
ありがとうございます
でもこれらは非等方圧力ではないですよね？

**名無しさん＠５周年** · 2014/05/10(土) 01:11:47.05

>>1
うちの業界だけの言葉なんですかね？
圧力テンソルの対角成分が同じものでないことです。
例えば二次元なら
∇・P = (∂x Px) ex + (∂y Py) ey

**名無しさん＠５周年** · 2014/05/10(土) 10:19:13.53

>>1
独立な圧力がP一つではなくてPxとPyなんですよ

**名無しさん＠５周年** · 2014/05/10(土) 10:22:08.14

流体中での爆発とか？

**名無しさん＠５周年** · 2014/05/10(土) 10:25:42.51

それもあるかもしれないですね。
天体のジェットなんかを考える時はジェットの吹き出る方向と、それに垂直な方向では温度が違うんですよ

**名無しさん＠５周年** · 2014/05/10(土) 10:27:31.04

P = P⊥ diag(1, 1, 0) + P|| diag(0, 0, 1)
なんてかいたりしますね

**名無しさん＠５周年** · 2014/05/10(土) 10:29:02.72

通常の流体(等方圧力)では
P = P diag(1, 1, 1)なわけです

**名無しさん＠５周年** · 2014/05/10(土) 15:25:31.05

分子運動論的な言い方をすると、方向ごとに熱速度の大きさが違うので
面の方向ごとに圧力が違うってこと？
そういう世界もあるんだなあ。状態方程式も圧力の複数の成分が入ってるのかな

**名無しさん＠５周年** · 2014/05/11(日) 00:50:24.75

物理現象としては同じものを対象としてないって理解でOK?

**名無しさん＠５周年** · 2014/05/11(日) 01:49:01.72

>>30
同じものを対象にしていないとはどういう意味でしょうか。

http://www.physics.wisc.edu/grads/courses/726-f07/files/Section_14_Eq_General.pdf
みたいな感じですね。これは磁場が入ってますけど

**名無しさん＠５周年** · 2014/05/20(火) 01:03:01.31

初心者ですが、二流体モデルをやるにはどうすればよいんでしょうか？

**名無しさん＠５周年** · 2014/05/20(火) 04:53:37.90

>>32
collisionいれればいいんじゃない？

**名無しさん＠５周年** · 2014/05/21(水) 20:35:57.04

とりあえず、流体がシミュレーション通りに動くことを己の目で見たいのだけど
どこに行けばいい？

**名無しさん＠５周年** · 2014/05/22(木) 12:31:38.82

>>1
>>69
ひとまずもう少し勉強してみます。ちなみに液液です

**名無しさん＠５周年** · 2014/05/29(木) 14:41:12.86

fortranのmpi fftwのサンプルコードとか持ってる人いるかな．
fftwはドキュメントが不親切過ぎて．．．

**名無しさん＠５周年** · 2014/05/30(金) 18:38:11.09

>>1
そしたらみなさんFFTってMPIでどうやってるんですかね？
ライブラリ使わず？？
FFTW以外のライブラリを使う？？

**名無しさん＠５周年** · 2014/06/07(土) 16:57:23.82

前処理付き共役勾配法が今のところ最強ってことでいい？

**名無しさん＠５周年** · 2014/06/11(水) 02:52:05.71

なるほどね！今は無難に不完全コレスキーだから、ほかのも試してみようかしら

**名無しさん＠５周年** · 2014/06/11(水) 15:29:04.25

は日曜日以外毎日8時から18時過ぎまで長期間に渡って周囲に騒音をまき散らさないように配慮して工事する金儲けより公共の利益を重視するすばらしい会社です

**名無しさん＠５周年** · 2014/06/11(水) 16:40:17.83

朝日住宅は日曜日以外毎日8時から18時過ぎまで長期間に渡って周囲に騒音をまき散らさないように
配慮して工事する金儲けより公共の利益を重視するすばらしい会社です

**名無しさん＠５周年** · 2014/06/16(月) 00:27:44.16

Poisson方程式を周期境界条件で解きたいのですが，行列のrankが落ちてるので解けません．
どうしたらよいかご存知でしょうか．

**名無しさん＠５周年** · 2014/06/16(月) 23:58:13.74

残差切除法ってどうなんすかね？poisson方程式の解法として

**名無しさん＠５周年** · 2014/06/17(火) 00:33:56.75

>>42 計算格子のどこか一点(四隅のどこかとか)を0とかに固定すればおk。

**名無しさん＠５周年** · 2014/06/17(火) 22:59:47.33

ハカセじゃないと言葉が流暢だなぁ

なんだか面白いんでこのスレに張り付いてるんですが、
高卒程度の知識で楽しみながらやる勉強法ってありますか？

**名無しさん＠５周年** · 2014/06/17(火) 23:51:45.78

なんだ、今期はリドルにハマってんのか

**名無しさん＠５周年** · 2014/06/18(水) 03:21:14.43

>>1
すんません
出口側を切り落として行列を作り直すとはどういうことでしょうか

**名無しさん＠５周年** · 2014/06/19(木) 00:05:19.04

>>1
内点をf1~fnとするとf0=fn, fn+1 = f1という境界です．行列でかくと
2 -1 -1
-1 2 -1
....
-1 2 -1
-1 2 -1

右上と左下の-1が周期境界に対応してます

**名無しさん＠５周年** · 2014/06/19(木) 00:09:57.84

行列わかりづらいですね．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}
です

**名無しさん＠５周年** · 2014/06/20(金) 04:34:21.33

>>1
これは拡散方程式ですよね。しかも陽解法で解いてます。
ちょっと計算してみましたが仮に拡散方程式を陰解法で解いても行列はランク落ちしないです。

単純に考えてみると、周期境界のポアソン方程式は定数分の不定があるからランク落ちすると思うんですが、拡散方程式は方物型だから、コーシーデータが必要になるので、定数分の不定もなくなりますよね

**名無しさん＠５周年** · 2014/06/21(土) 00:16:21.99

>>1
フーリエで解けるのはわかるんです。0モードを考えないから

僕は差分でときたいんですねー

**名無しさん＠５周年** · 2014/06/22(日) 02:57:14.13

ひきこもり30代後半で童貞で低収入イラストレーターやってるバカ発見。
足立区に住んでいるそうだ
http://inumenken.blog.jp/archives/6580395.html

**名無しさん＠５周年** · 2014/06/25(水) 01:35:50.45

>>1
擬似逆行列ですか！勉強してみます

ちょっとおもったんですが、平均値がゼロという条件を課せばランク落ちしないんですが、帯行列でなくなるために随分ときにくくなって却下しました

**東日本大震災** · 2014/06/25(水) 16:46:41.07

あんなー目撃証言0とかでこういう事をするとなー東日本大震災とかになるだわ。
理解したー？

26才を超えてるのに精神病院がどうのこうのとか。
俺は小学生じゃないぞ。

これだけ年を取ってるのにお前はいったいなんなんだ？
友達じゃねーんだ。

**東日本大震災** · 2014/06/25(水) 16:50:34.60

あーはいそうですかー？
いちいちその時にそうやって祈っててどうもすいませんでした。
祈ってないのになー。

**天皇陛下** · 2014/06/25(水) 16:52:02.83

新聞配達員の最中に食べに行った刈谷市の辺のファミリーレストランを解らなくして東日本大震災

90 · 2014/06/26(木) 20:05:47.23

一点固定の方法でキャビティ流れをやろうとしましたが、うまくいきませんでした...。キャビティ流れでポアソン方程式を解く場合、行列のランク落ちをどうすればいいんでしょう？

**名無しさん＠５周年** · 2014/06/27(金) 00:56:16.81

>>57
あら、私(>>88)と同じ悩みを持つものが

**名無しさん＠５周年** · 2014/06/28(土) 09:53:19.28

>116
>>42ですが、僕がやって見た感じだと、sorのような収束が遅いスキームだと収束するんですが、CG法みたいに速いスキームのときだと収束しない場合が多いです

90 · 2014/06/28(土) 14:51:52.56

>116 解けました！速度の境界条件をミスしてました。

**名無しさん＠５周年** · 2014/07/02(水) 08:50:06.81

30代後半で童貞ひきこもり低収入イラストレーターを発見
足立区に住んでいるそうだ。
http://inumenken.blog.jp/archives/7002197.html

**名無しさん＠５周年** · 2014/07/13(日) 13:25:15.74

ネット版のホストみたいなの発見。

イケメンなら稼げるんだろうけど。
話が上手けりゃ稼げるかな。

誰かレポ頼む。
メンガでググると出てくる。

**名無しさん＠５周年** · 2014/07/21(月) 18:30:09.10

計算機の性能向上って頭打ちなのかなぁ…
ここ数年伸びが遅くなってる気がする

**名無しさん＠５周年** · 2014/07/23(水) 20:25:11.18

ここまで中身のある回答６２９のみｗ

**名無しさん＠５周年** · 2014/07/27(日) 17:14:53.40

qqqqqqaqaqaqaaqaqhttp://iup.2ch-library.com/i/i0289053-1303372232.jpg
http://iup.2ch-library.com/i/i0289054-1303372232.jpg
WqqAaqqQqqqqqqqQz いつの日か参考出品さaaqaqaqqれてた奴
AaaqQaaaaqaaaqqAqqaqaaqaqaqqqqaaqqqqaqaqqqaqqaqaqqqAaqqaqaaahqdvqqqqqqqqw

**名無しさん＠５周年** · 2014/07/27(日) 17:19:18.11

AaSAasAAqaaqq(qqqqa(AAaaaqaaaqqqq(ww(wwew))))

**名無しさん＠５周年** · 2014/07/28(月) 10:58:45.65

https://www.youtube.com/watch?v=LGlvDb3h4WA&;feature=youtu.be

**名無しさん＠５周年** · 2014/07/28(月) 12:58:59.54

>>1
民生品で一番進んでいるであろうゲームでさえも
マルチコア対応を堂々と謳ってるのって少ないし…
GPGPU…うーん…。。。
需要の問題か…？飽和気味なのかもね。

**名無しさん＠５周年** · 2014/08/03(日) 23:28:34.39

メニーコアCPU出ても、FDMとかFVMなら結局メモリバンドぶんしか
性能伸びない気がするんだけどなぁ。。。
OpenMPはスレッド数増やすと性能伸びないし。

**学会で変死** · 2014/08/06(水) 04:07:03.64

辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。

**クラークの超必殺** · 2014/08/06(水) 04:14:30.10

辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。
辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。辞めてストーカー。

**学会で変死** · 2014/08/06(水) 04:16:35.18

ふっふっふ時代の流れにインプット

マジレスを手伝った事にしたり首謀者のゴミにゴミに出された八つ当たりをして
無差別殺人か学会で変死

**DAY GAME** · 2014/08/06(水) 04:17:29.95

ふっふっふ時代の流れにインプット

首謀者のゴミのいじめを手伝って無差別殺人か学会で変死

**DAY GAME** · 2014/08/06(水) 04:18:43.14

ふっふっふ時代の流れにインプット

水槽の水を替える所から無かった事にして無差別殺人か学会で変死

**DAY GAME** · 2014/08/06(水) 04:35:32.27

ふっふっふ時代の流れにインプット

助けないと解らない事して無差別殺人か学会で変死

**名無しさん＠５周年** · 2014/08/10(日) 14:57:14.29

http://nonke.blog.jp/
載ってる

**名無しさん＠５周年** · 2014/08/14(木) 19:36:27.45

STAP騒動を見てて疑問に思ったんだけど、シミュレーションコード開発ってやろうと思えば実験より簡単にデータの捏造ができるとおもうんですよね。

複数人で同じコードを開発してれば、そいうのは防げるし、下手くそなコード書いてる人の指導とかもやりやすい
皆さんはコード開発は何人くらいで同じコードを書いてますか？

**名無しさん＠５周年** · 2014/08/16(土) 07:26:03.17

デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね
デブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ねデブ豚死ね

**名無しさん＠５周年** · 2014/08/16(土) 15:20:28.62

>>2 のフリーソフトのどれがお勧め？
どれ使ってやってるの？

**名無しさん＠５周年** · 2014/08/17(日) 08:26:17.20

>>148-151
詳しい情報ありがとう。
OpenFOAMじゃなく、Channelflowで遊ぶ理由はなに？
なにか優位性あるの？

いまのところ使う予定は未定。よっぽど暇なときにちょっと遊ぶかもって程度かな。
計算対象もいまのところ興味もってるもので無いし。

**名無しさん＠５周年** · 2014/08/17(日) 08:32:55.48

>>1
うちの場合みんなテーマがバラバラで、コード作りは各人任せなのよね。
それでほかの人からはノーチェックだから、仮になにか悪さしようと思ったらできる

まあ実験でも同じようにノーチェックなとこもあるのかもしれんけども

**名無しさん＠５周年** · 2014/08/17(日) 15:58:20.40

>>1
詳しい回答ありがとう。

**名無しさん＠５周年** · 2014/08/19(火) 18:58:22.09

プログラム言語の入門で最初によくでてくるのがHello World。
ワンチップマイコンだとLEDの点滅。

流体力学シミュレーションの最初の例として一般的なのはどんなの？

**名無しさん＠５周年** · 2014/08/20(水) 15:33:15.91

>>1
教えてくれてありがとう。
こんど試してみるときに、それのどれかをやってみることにします。

**名無しさん＠５周年** · 2014/08/21(木) 21:07:09.59

>>83
一次元拡散方程式(熱伝導方程式)を FTCSで解く、は如何でしょうか。
解けたあとはクーラン数の勉強もできます

**名無しさん＠５周年** · 2014/08/21(木) 22:43:52.50

>>85
ありがとうございます。
熱伝導方程式をとくんですか？
流体関係のサンプルじゃなくて？

まだソフトを調べてないんだけど、微分方程式とか偏微分方程式を
数値的に解くソフトであって流体力学専用ソフトじゃないの？

**名無しさん＠５周年** · 2014/09/08(月) 19:52:18.85

>>1
詳しい解説ありがとう。
流体力学受講しなかったんだよな。
受講してても使わないから忘れてたろうけど。

偏微分方程式を解くプログラムで
流体力学で解く問題に関係しそうなものが得意なようにチューニングされてるソフトって解釈したけど間違ってないかな？

前色々一般的な偏微分方程式を解くソフトを探してたときに
ＧＳＬライブラリかなにかに、
偏微分方程式を解くときはケースバイケースなのでGSLには含むことを考えなかった的な文書を読んだ記憶がある。

**名無しさん＠５周年** · 2014/09/08(月) 21:59:11.54

素人でもできる解析ソフトってあります？

**名無しさん＠５周年** · 2014/09/08(月) 22:12:24.70

自作ってスレについてるから
主旨に反するけど、自作しまくってたら
作ってる人がよっぽど作るの上手くないとバグだらけってことになりませんか？

作ったものを,githubとかにあげてもらって、多くの人が使いやすいようにして
バグを利用者にテストしてもらって発見しやすくした方が良いと思う。

バグ報告管理機能がgithubについてたかどうか知らないけど、それがあるところで公開したほうが。

それすると使い方とかを文書にする手間が発生するけど。
通常プログラム作る人は作って計算するのが主目的で
使い方は作った本人が良くわかってるから文書書くのは苦痛だと思いますが。

**名無しさん＠５周年** · 2014/09/11(木) 01:52:55.79

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]の．．．
の間違いです

**名無しさん＠５周年** · 2014/09/11(木) 08:02:41.78

波数で円を書いて2/3法内に入っている波数を使えばいいと思う。
√kx^2+ky^2である方向の波数の大きさわかるので。

**名無しさん＠５周年** · 2014/09/11(木) 08:13:47.46

>>91
なるほど！！！
ありがとうございます！

**名無しさん＠５周年** · 2014/09/12(金) 00:59:06.64

>>1
すみません！意味がわかりません！！

>>186
ちょいまち！これは>>183とはちがいますよね？>>183のやり方なら左右上下のカドは0埋め領域になるのでは
これはむしろ>>181で言うところの和集合になってるとおもいますが！

**名無しさん＠５周年** · 2014/09/12(金) 01:13:37.35

>>1
乱流になりきれないとは具体的にどのようになるの？

**183** · 2014/09/12(金) 07:18:54.28

混乱させてしまって申し訳ありません。
勝手に等方性を仮定してしまって、向きによってモード数が異なることに気持ち悪さを感じたので、ああいう書き方をしてしまいました。

Nx,Ｎyと書いてくれていたので等方である保証一切ないのに><

**183** · 2014/09/12(金) 07:27:55.48

ん。
質問本文はＮでした。

えっと2/3法を満たすという意味では、DNS厨さんの絵がそのものずばりで、そっから更に各方向のモード数が等しくなるように切り捨てると183になります。エイリアジングエラー除去という意味では過剰です。等方性を神様のように崇める人達以外は気にしなくて良いと思います。

**名無しさん＠５周年** · 2014/09/12(金) 16:51:43.87

　___　_
　ヽo,´-'─ 、　♪
　　 r, "~~~~"ヽ
　　 i. ,'ﾉﾚノﾚ!ﾚ〉　　　　☆ 日本のカクブソウは絶対に必須です ☆
　__ '!从.ﾟヮﾟﾉﾙ　　　総務省の『憲法改正国民投票法』のURLです。
　ゝﾝ〈(つY_i(つ　http://www.soumu.go.jp/senkyo/kokumin_touhyou/index.html
　｀,.く,§_,_,ゝ,
　　　~i_ﾝｲﾉ

**名無しさん＠５周年** · 2014/09/12(金) 21:59:51.31

皆さんご親切にありがとうございました

博士の図の様にやります！

**名無しさん＠５周年** · 2014/09/13(土) 12:13:35.08

なるほど。これは確かにええっすな

**名無しさん＠５周年** · 2014/09/16(火) 01:54:06.12

私もスペクトル法で質問があります
2D非圧縮NavierStokesを差分で解く時は非線形項をArakawa Jacobianで計算するのが一般的だと思うのですが、スペクトル法のときは気にせずそのまま計算していいのでしょうか

**名無しさん＠５周年** · 2014/09/16(火) 22:11:07.63

rotation form…はじめて聞きました

**名無しさん＠５周年** · 2014/09/22(月) 16:57:37.90

スペクトル法では、Rotational 型はだめで、Skew symmetric 型を使えということだったのでは？
ところで、-N/2, N/2 の領域になっているのは何か意味があるの？

**名無しさん＠５周年** · 2014/09/26(金) 00:49:08.11

FFTWのMPIって一次元分割なんかいっ！
ガッカリだよ

**名無しさん＠５周年** · 2014/09/30(火) 16:17:11.48

すみません、質問です
流体構造の連成解析などでは、弾性体の固有振動数が鍵になったりします

どうして固有振動数を調べる必要があるんでしょう？(´･ω･`)

**名無しさん＠５周年** · 2014/09/30(火) 17:21:53.69

自己解決しました、自励振動に関わってるのですね

**解析厨** · 2014/10/01(水) 05:12:47.66

ﾏﾙﾁｺｱ＆並行処理ﾌﾟﾛｸﾞﾗﾑ組める方望む
RSA Numbers 素因数分解に挑む・解析数式の考察
ttp://blog.livedoor.jp/superprojectx/archives/1009232905.html

**名無しさん＠５周年** · 2014/10/01(水) 12:45:41.28

>>1
円環の理ですね、わかりました(`･ω･´)ﾉｼ

**Q&A** · 2014/10/03(金) 16:30:59.24

ET原子
原子遊び
原子核抜き

水素がきめて　　2億年後
速い2000km/hタコメーター
車　電車など

水素崩壊　
★★『幸福の科学』統合スレッドpart761★★
http://wc2014.2ch.net/test/read.cgi/psy/1410957061/
Q&Aより

**名無しさん＠５周年** · 2014/10/07(火) 21:29:57.46

てんてー、相似則の適用範囲ってあるの？
なんで成り立つのか分かんないっす

**名無しさん＠５周年** · 2014/10/13(月) 23:13:22.82

天気予報とか台風の進路ってどうやって予測してるんですかねぇ
今回の台風みたいに九州で急に東に進みだすとかいう挙動も，
シミュレーションで再現できるんでしょうか・・・＞＜

**名無しさん＠５周年** · 2014/11/01(土) 11:25:27.55

シミュレーションていうけど人の死や経済・法律的に責任問題が生じうるような場面では
意思決定で使えないというのが原子力安全上の判断になったよね。

機械設計でシミュレーションは応用可能だろうが、危機管理の意思決定において少しでも不確実性が残る。
シミュレーション結果(ここではSPEEDI)を今までは判断基準としていたが、事故後の現在は参考資料に留めたよね。
これって結局、高レベル放射線は実測する必要があり事故の様子も確認しないと危機管理の意思決定には使えない
ということでしょ。
責任問題がいざ発生したとき、シミュレーション結果(確率的な誤差や曖昧さ含む)だけで避難指示と救援行動することは
危険だという事でしょ。
人が死んだり、天文学的な経済損失が発生して、一回きりで後戻りできない状況において
果たして理論上のシミュレーション結果だけで意思決定できるものなのかどうか。

**名無しさん＠５周年** · 2014/11/01(土) 12:53:59.22

今のシミュレーションはそこまで精度がないってことだろ。
条件がはっきりしてて仕組みもはっきりしてたらシミュレーションで十分って話になると思うよ。

例えばエミュレーターとPIC,PICの実機と比較すると
エミュレーターで遊んだ方がお金もかからないしほぼ一緒。

異なるのはチャタリングとか過電流がながれないとか。

原子力はもともと地震大国の日本で作るべきじゃなかったし、
継続決定は馬鹿すぎる選択だと考えてる。
頭悪すぎる非常に重要な選択を最悪の選択対応するのが得意。

そのつけはすぐにはっきりすることになるよ。

**名無しさん＠５周年** · 2014/11/01(土) 13:15:19.93

頭の悪い文章をここでひけらかすのやめてくれない？

**名無しさん＠５周年** · 2014/11/01(土) 13:18:58.20

大丈夫とか大したことになるはずないって
思いこんでるお前の方が馬鹿だよ。

薬ですら人体実験してから投与なのに。

俺らで今実験中って認識すらないだろ馬鹿すぎて。
どうなるかは先になってみないとわからない。

多分アウト。大量に漏れたら安全になるなんてことになるはずないだろ。

**名無しさん＠５周年** · 2014/11/01(土) 17:20:45.46

http://www.genanshin.jp/report/column/wave_20141029.html

**名無しさん＠５周年** · 2014/11/05(水) 23:37:36.77

円筒座標系で周方向に疑似スペクトル法（半径、軸方向は差分）を使って解こうとしているのですが、
時間進行方法や圧力の解き方がわかりません。
ご存知でしょうか。

**名無しさん＠５周年** · 2014/11/07(金) 07:29:22.09

軸中心付近のメッシュが細かくなりすぎるのとr=0近傍が特異点のように見えて鬱陶しいのはまた別のお話。。。

**名無しさん＠５周年** · 2014/11/09(日) 07:14:35.24

>>116です。
アルゴリズムはfractional step法を使って、圧力はsor法で解いて見ようと
思います。今回の計算対象は中心の流れがほとんど無いのでスリップ条件にしよう
と考えています。
貴重なご意見をありがとうございました。

**名無しさん＠５周年** · 2014/11/11(火) 04:37:11.94

http://monoist.atmarkit.co.jp/mn/articles/1410/31/news023.html
スパコン性能、2020年以降は伸びずに頭打ち!? ―その時CFDでは何ができるだろうか。東大生産研講演

NASAやBoeingは2030年に1ZFLOPS(1000Exa)級運用を目標にしてるらしいが、果たして現状を理解しているのか。
まだ性能向上が見込めるGPUクラスタでぶん回すのかねえ？w　死ぬ程効率悪いだろ。

**名無しさん＠５周年** · 2014/11/11(火) 17:58:54.92

いま世の中では粒子法と格子法どっちが主流なの？

**名無しさん＠５周年** · 2014/11/27(木) 01:29:30.31

2014年のアンシスカンファレンスは以下で行ないます。
http://www.youtube.com/watch?v=HJKHJEJgxgo
ユーザーの方もユーザーでない方もぜひご参加お待ちしております。

**名無しさん＠５周年** · 2014/12/07(日) 14:52:23.79

もうすぐ数値シンポですね。
ここの皆さんは参加されますか？

**名無しさん＠５周年** · 2014/12/07(日) 21:21:15.61

残念だが行けないorz

**学会員化した防犯ネットワーク** · 2014/12/20(土) 23:18:55.88

ビデオのテープの部分をおかしくして無差別殺人か学会で変死

**学会員化した防犯ネットワーク** · 2014/12/20(土) 23:20:11.42

うんだからー。
ドミーのトイレで田舎物とかやってこなくていいの。

違う違う違う違う。お前断られとるの。

**名無しさん＠５周年** · 2014/12/21(日) 22:59:48.96

原発推進法人年収９００万 / ４２歳ボーナス２３０万すべて国税

原子力研究開発機構原発推進国策法人
原子力研究開発機構原発推進国策法人

税金でボーナス２３０万原発推進東電天下り集団
税金でボーナス２３０万原発活用東電天下り集団

原子力研究開発機構東電天下り国税２３００億円じゃびじゃびゅ
原子力研究開発機構　東電天下り　国税２３００億円じゃｂじゃびゅ

**名無しさん＠５周年** · 2014/12/26(金) 07:00:32.85

PSEの良さってどんなところにあるの？

**名無しさん＠５周年** · 2014/12/27(土) 00:02:42.24

BSEは脳みそがスポンジになること

**名無しさん＠５周年** · 2014/12/27(土) 15:53:05.13

はかせー

http://www.isas.jaxa.jp/j/topics/topics/2008/0911.shtml

こういう流体の動きの模擬って、どの程度完璧なの？
それと、どうやったら、個人で同じようなこと出来ますか？

**名無しさん＠５周年** · 2014/12/27(土) 22:12:25.02

fringeしなくて良い予感。。。

**名無しさん＠５周年** · 2014/12/27(土) 23:51:47.83

>>129
うわー、ありえねー

**名無しさん＠５周年** · 2015/01/07(水) 23:31:40.77

遅くなったけど、教えてくれてありがとう、はかせ

**名無しさん＠５周年** · 2015/01/21(水) 11:21:33.71

俺はいじめられっこのゲイ野郎だぞ！AT和え

**名無しさん＠５周年** · 2015/01/22(木) 22:12:02.48

補間法について教えてください．

非構造格子（四面体）のセル中心で得られた流速を補間して節点での流速を求めたいのですが，
どのような方法がありますか？手法の名前や解説サイトでも助かります．

境界面の節点で滑りなしの速度境界条件を満たし，補間した値が二階微分可能であればなおうれしいです．

**名無しさん＠５周年** · 2015/01/24(土) 00:30:20.40

シミュレーションについて語り合う仲間がほしいならここがおすすめ。
http://blngs.com/
一回使ってみると、意外と使えるかと。

**264** · 2015/01/24(土) 21:46:26.98

>>1
ありがとうございました．
5章の最後がそれっぽいですね．ちゃんと読んでみようと思います．

**名無しさん＠５周年** · 2015/01/27(火) 17:04:12.59

スペクトル法に関して教えて下さい。
以下の論文では、位相シフト法をつかっていると記載していますが、
図1(b) では最大波数で大きな跳ねあがりが見られます。
これって普通なんですか？　位相シフト法の問題ですか？
https://www.jstage.jst.go.jp/article/kikaib/79/803/79_1281/_pdf

**名無しさん＠５周年** · 2015/01/28(水) 07:56:50.85

正カスケード(低波数から高波数)系の統計量なら十分な散逸領域が無くて最高波数のところに溜まっちゃうのでしょうね。

中間波数領域と違って散逸以外に逃がす先が無いから。

最高波数への供給が粘性散逸より大きいという話だと思います。

差分とかで起きなかったとしたら差分とったときに高波数を潰すような誤差を作り出すのだと思われる

**名無しさん＠５周年** · 2015/01/29(木) 19:32:14.91

http://jbbs.shitaraba.net/sports/42269/

**名無しさん＠５周年** · 2015/01/30(金) 04:56:59.94

☆☆☆☆☆
☆ 自民党、グッジョブですわ。 ☆
http://www.soumu.go.jp/senkyo/kokumin_touhyou/index.html

☆ 日本国民の皆様方、2016年7月の『第24回参議院選挙』で、改憲の参議院議員が
3分の2以上を超えると日本国憲法の改正です。皆様方、必ず投票に自ら足を運んでください。
そして、私たちの日本国憲法を絶対に改正しましょう。☆

**名無しさん＠５周年** · 2015/01/30(金) 09:40:14.75

>>1 >>270
速度場はきれいに落ちているので解像度の影響ということですね。

現在、1次元の位相シフト法ではうまくいくのですが、2次元の場合、高波数が跳ねてしまうので、
スキームの問題かと思いました。なにかミスがあるのだと思います。

ありがとうございました。

**名無しさん＠５周年** · 2015/01/30(金) 13:02:17.05

自作の差分法DNSプログラムでチャネル乱流計算をしています．
流れ方向の平均流速を出したところ，粘性低層は合うのですが，
対数則領域が過小に評価されてしまいます．
この原因は何だと考えられますでしょうか？(何らかの数値粘性で
チャネル中央部の流れが減衰してしまっているのでしょうか？？)
また，こういうことは，よく考えてプログラムを書かないと
しばしば起こるものなのでしょうか？？
なお，計算格子幅(壁座標)は壁方向が3で，流れ方向とスパン方向は30です．
また，対流項の離散化は2次の中心差分です．
どのようなアドバイスでもいただけたら大変ありがたいです．
よろしくお願いします．

**名無しさん＠５周年** · 2015/01/30(金) 14:11:20.01

>>142
流れ方向とスパン方向の解像度が粗すぎる。
平均量程度でも、dx+=18、dz+=9 ぐらいが必要。
乱流統計量まであわせたければ、さらにその倍は必要。

2次の中心差分で計算していれば数値粘性は入らない。

**名無しさん＠５周年** · 2015/01/30(金) 14:38:44.87

>>143さん
早速アドバイスをいただきありがとうございます．

格子数が足りな過ぎるのが原因だったのですね．
ご指摘の通りに細かくしてみます．
おっしゃるとおり中心差分には数値粘性が入っていませんので
私は相当とんちんかんなことを言っていますね．
有用なアドバイスを頂きありがとうございました．

今後LESもコーディングしたいと思いますが，その場合でも同程度の
解像度(dx+=18、dy+=3、dz+=9)が必要でしょうか？

**名無しさん＠５周年** · 2015/01/30(金) 15:02:04.64

>>144
森西先生の以下の論文が参考になるかと、
https://www.jstage.jst.go.jp/article/kikaib1979/66/647/66_647_1750/_pdf
解像度が粗いためにDNSでは散逸を捉えきれない影響を、SGS応力で補う観点からは、
論文中の図5に着目して、流量が減少する領域でLESの格子を設定する必要がある。
2次精度の場合はこの領域がかなり狭いということですね。

**名無しさん＠５周年** · 2015/01/30(金) 17:06:15.41

>>145さん
大変参考になる文献をご紹介いただきありがとうございました！
当方勉強不足なのでしっかりと勉強してみたいと思います．
図5によると，2次精度中心差分で精度を出すためには，
x+は20～40(z+は6～13)程度に選ぶ必要があるのですね．

当方の頭が悪くまだよく理解できていないのですが，
図5によると，格子が粗い場合の2次精度中心差分では
平均流速を過大に評価する傾向にあると思います．
当方の計算ではなぜ過小に評価されてしまったのでしょうか？？

**名無しさん＠５周年** · 2015/02/05(木) 12:53:46.07

>>146
図5において、2次精度の場合、dx+ < 40 ぐらいまでは、流量を過小評価している。
dx+ > 50 以上になると流量を過大評価している。
つまりあなたの計算結果はこれにあっていると言えると思います。

解像度を粗くして流量が減るのは、散逸が捉えられないため。
LESではGSで捉えられない、この散逸をSGSで補う。
従ってこの観点からは、格子分解能を流量が減少している領域で設定する必要がある。

**名無しさん＠５周年** · 2015/02/20(金) 11:34:38.75

か

◆cqxclbaQk/hG · 2015/02/21(土) 15:58:46.98

最初に無があった
無から有が生まれた
これが全ての真理

**名無しさん＠５周年** · 2015/03/07(土) 00:23:18.25

最近は擬似圧縮が流行っているのですか？

**名無しさん＠５周年** · 2015/03/09(月) 21:41:00.46

http://jbbs.shitaraba.net/sports/42269/

**名無しさん＠５周年** · 2015/03/14(土) 18:39:41.85

流体は関わりはないが、シミュレーションだけに偏るのはどうかと思う。
時々天気予報が大幅に外れることがあるけれど、あれはシミュレーション結果を過信
しすぎたためと思われるが。

TVの天気解説にもGPU使ったような風向きベクトル示した短時間のシミュレーション動画が
流れてるが実際の大気の流れはあのように単純ではないよなといつも思うのだがｗ

**名無しさん＠５周年** · 2015/03/15(日) 10:52:25.06

結局、全ての分子について運動方程式を解いたとしても精度は上がらないだろうから、所詮はシミュレーション。
春雨じゃ、濡れてまいろう。位の大雑把な気持ちでシミュレーションするのが病まなくて良い。

**名無しさん＠５周年** · 2015/03/16(月) 18:28:33.54

ビックカメ●札幌店の副店長
佐藤伸弦が暴行事件を起こしていた

**名無しさん＠５周年** · 2015/04/14(火) 15:19:13.85

http://i.imgur.com/V7P8RWS.jpg
やっぱうっすいから写真とるとかっこいいなぁ

**名無しさん＠５周年** · 2015/04/18(土) 10:12:47.84

また髪の毛の話してる

**名無しさん＠５周年** · 2015/04/21(火) 23:57:02.61

（2ch．scの方にも書きましたが、住人がいるか心配になったのでこっちにも書きます）

Sand management / Multiphase flows
http://www.lemma-ing.com/?page_id=976
http://www.lemma-ing.com/wp-content/uploads/2012/09/emulsion_shear_rate_comparison.png

例えばこの画像のように、せん断速度をスカラーの色コンターで表示するのをよく見かけるのですが、
自分の理解では、せん断速度は変形速度テンソル（行列）の要素だと理解していて、
どうして色コンターで表示できるのか分かりません。

自分の理解は間違っているのでしょうか？
それとも、この画像のshear rateは何か特別な処理をしてスカラー量に変換されたものなのでしょうか？

**名無しさん＠５周年** · 2015/04/22(水) 12:17:11.72

>>157
絶対値じゃないですかね

**294** · 2015/04/24(金) 00:11:02.30

御返事有難う御座います。

もし3次元の場合、見る方向によって同じ位置でもせん断速度の
大きさが変わるということですか？

**294** · 2015/04/28(火) 08:02:18.56

>それぞれの面に働くせん断速度が異なる
そうでした。ありがとうございました。当たり前でしたね。

**名無しさん＠５周年** · 2015/04/29(水) 23:07:00.78

有限差分法について質問があります。
セルの頂点に座標と全ての物理量（速度，圧力など）が定義されているタイプの離散化と，
セルの頂点にその座標を，また有限体積法のようにセルの中心で物理量が
局所値として定義されているタイプの離散化があることに気づきました。
（有限体積法では体積平均として定義されると思います。）

これらにはそれぞれメリット・デメリットなどあるのでしょうか？
後者は前者に比べ差分の計算の際に計算量が多くなってしまうように思いますが，
その分メリットなどあるのでしょうか？

**名無しさん＠５周年** · 2015/04/30(木) 07:34:34.55

スタガード格子のこと？
全ての物理量が同じ位置に置かれるのではなくて、格子点上に圧力やその他のスカラー、dx/2ずらした位置にu、dy/2ずらしてv、dz/2ずらしてwみたいになるはずだけど。
メリットは解が振動することに対する対策。

**名無しさん＠５周年** · 2015/05/01(金) 08:35:02.24

>>162
スタガード格子のように圧力と速度の定義位置をずらすのではなく，
すべての量がセル頂点にあるのか，セル中心にあるのかという違いのことです。
境界条件の与え方が違ってくると思います。
（たとえば壁面上にも格子があるのか，半セルずれた位置で与えるのか）
簡単なコードでどう違うか確認してみようにも，ゼロから作るのはちょっと大変かな・・・

**名無しさん＠５周年** · 2015/05/01(金) 09:39:56.61

それは例えば dxdy/dxdydz=1/dz みたいな形に約分されて有限差分法に見えてるだけで実際には有限体積法なのではなくて？

**303** · 2015/05/01(金) 18:19:44.88

そうか座標変換が絡むと物理量が同じようにセル中心にあってもFVMとFDMで式が厳密に同じにはならなくなるんだな

**名無しさん＠５周年** · 2015/05/04(月) 10:51:44.44

>>1 >>305
お返事ありがとうございます。
>>304で挙げされたのは非圧縮性ソルバの話，
>>305では非構造圧縮性ソルバ？の話で興味深く拝見しました。
非圧縮の場合はコロケートにしてかつ内挿・差分の扱いに気をつける必要があるみたいですね。
圧縮性非構造のほうでは検査体積の形や物理量再構築の際用いる量の取り方で精度が変わってくるということでしょうか。
まだ理解しきれていないので，じっくり読んでみたいと思います。
ちなみに私は構造格子圧縮性流のコードをコロケート風に書いています。
違う離散化で比較していないので，いまいちメリットがあるのか分からなくなり，質問させていただきました。

**名無しさん＠５周年** · 2015/06/11(木) 23:18:05.43

DualSPHysics使っている人いますか？
フリーの流体ソフトなんです。日本語の情報が無に等しいので日本では
まだ使ってる人いないのかな？
Windows64bit用をダウンロードしてみたけど使う以前に起動することすら苦労しそうです。

**名無しさん＠５周年** · 2015/06/12(金) 18:45:48.83

↑ごめんなさい
>>1 のまちがいでした。

**311** · 2015/06/18(木) 20:32:10.93

>>1 DualSPHysicsきっと楽しいですよ～
翻訳(ほとんど機械翻訳ですが)もある程度進んだので、
これからちょっとずつサンプル見ながらやっていこうと思います。

Blenderはそこそこ使えますがLBMは使いこなすまではまだ。。
少し手のこんだことすると途端にPCが激重になるので
シンプルなもの以外はあまり試せていないという状態です。
物理学としての流体自体の知識はさっぱりないです。

ところでAAかわいいですね。どうやって描くのか感心します。

**名無しさん＠５周年** · 2015/07/11(土) 11:49:38.01

近所に住んでいるホモ発見アプリ
iOS版
http://itunes.apple.com/jp/app/9monsters/id451345005?mt=8#

Android版
https://market.android.com/details?id=jp.co.applibros.alligatorxx

http://ninemonsters.com

**名無しさん＠５周年** · 2015/07/14(火) 00:55:27.64

LESについて質問させてください．
フィルタ化したナビエストークス方程式の粘性項は
つぎのようにあらわせるとあります．
∂/∂x_j {2 * (nu + nu_t) * D_{ij}}
例えばスマゴリンスキーモデルでは
乱流粘性係数nu_tの項は空間的に変化するのでしょうか？
いい換えると
nu_t * ∂/∂x_j {2 * D_{ij}}
とnu_tを外に出してはいけないのでしょうか？
よろしくおねがいします．

**名無しさん＠５周年** · 2015/07/15(水) 01:35:53.69

>>1　紹介頂きありがとうございます．勉強します．
　　　　LESのコード作りにためになる自作・フリーソフトをご存知でしたら教えてくれませんか．
　　　　OpenFOAMにはあまり良さを感じないので除外でお願いします．

**名無しさん＠５周年** · 2015/07/15(水) 02:51:17.56

http://i.imgur.com/dgw83fk.png
http://i.imgur.com/KXDcsNd.png
http://i.imgur.com/K9Y4HeD.png
http://i.imgur.com/0wZLGwy.png
http://i.imgur.com/yGz1U0r.png
http://i.imgur.com/HIfm1uU.png
http://i.imgur.com/p7Ux1yk.jpg
http://i.imgur.com/avOjOHe.jpg

**名無しさん＠５周年** · 2015/07/15(水) 03:16:15.24

病気だね

**名無しさん＠５周年** · 2015/07/15(水) 03:17:48.85

あれですよ　あれ　
家猫の夜中の大運動会　　　　　　　　　　　　地獄じゃなくて電波

**名無しさん＠５周年** · 2015/07/19(日) 21:44:42.91

ツイッターである画像が話題になっています。
https://twitter.com/oosakitakashi/status/621669244934041600

テレ東で中核派として登場した女性が
TBSに、一般市民としてインタビューを受けていたようですが

これは仕込みでしょうか？？？

TBSは中核派を一般市民として扱ってるのでしょうか？

**名無しさん＠５周年** · 2015/07/23(木) 09:24:56.02

赤羽の美容室ブルームで犯罪行為
ブルームコスタカレンデコラ社長中野隆一『なかのりゅういち』は
犯罪者である
お客様を赤羽警察署に偽証申告をして逮捕状をとる
共犯デコラGM 遠藤孝輔『えんどうこうすけ』
ホットペッパービューティー
川口池袋評判縮毛矯正カットパーマ格安カーラー美容院床屋散髪
美容室ヘアメイクネールサロンザマスターキースリー
今夜は眠れないお楽しみ様です。

**名無しさん＠５周年** · 2015/07/30(木) 02:42:29.81

世界教師マイトレーヤは大宣言の日に全人類にテレパシーで語りかけます。
この日から地球の歴史上初めて真の世界平和が始まります。

人類は世界資源の分かち合いに同意するでしょう。

**名無しさん＠５周年** · 2015/08/07(金) 18:44:02.98

安倍は許した
http://live.fc2.com/79752690/

**名無しさん＠５周年** · 2015/08/19(水) 17:26:27.13

>>1
壁面で速度0、でお前の頭の中では、速度勾配du/dy も0なわけだな。
ってことは、壁で速度ゼロ、速度勾配もゼロ、壁垂直方向に積分して速度0となｗｗ

驚異の新事実　壁面流れでは流体は流れることができなかった！

死ねよ、バカ

**名無しさん＠５周年** · 2015/08/20(木) 12:14:46.41

メディカルサイト歯科相模原日吉サンテラス歯科藤沢なのはな内科スマイル歯科アイ整形外科亀有リリオ歯科八王寺アクロスみなみの歯科足立ハート新宿くろさか歯科熊本ファミリー歯科桜台歯科森林公園滑川モール歯科小田原めぐみ歯科保土ヶ谷西谷富士見おとなこども歯科

**名無しさん＠５周年** · 2015/08/22(土) 19:38:25.39

壁で0であると仮定されただけのdu/dyを>>180先生はどうやって積分したのでしょうか

**名無しさん＠５周年** · 2015/08/22(土) 22:07:20.68

>>182
じゃあ、なんで摩擦速度が0になるんだ。バカなのアホなのｗｗｗ
指摘されても間違いさえ気づけないのｗｗ　死ねよ、クズ。

**名無しさん＠５周年** · 2015/08/23(日) 00:05:00.07

>>1 バカはやっぱり死ななきゃ治らないってかｗｗ
壁関数が使える壁面流れで、壁垂直方向の途中から速度勾配が発生する要因を挙げてみろｗｗ

死ねよｗｗｗ

**335** · 2015/08/23(日) 00:08:42.54

実際にそんな流れが存在しないことは分かってるよ
ただ>>1と>>331で提示された条件だけでは数学的に導けていないと指摘しただけ

**名無しさん＠５周年** · 2015/08/23(日) 00:13:18.28

>>185
バカは惨めだねｗｗ　途中から速度勾配が発生する（スリップ）ように
お前の頭の中もスリップしてるんだろｗｗｗｗ

**名無しさん＠５周年** · 2015/08/23(日) 05:18:09.63

リアルで話し相手いないんだろうなぁ

**名無しさん＠５周年** · 2015/08/23(日) 07:02:36.34

>>187
実在しないのに、数学的に導けてないだって。
それならなんだってありだわな。
ぼくのつくった境界条件ですかｗｗｗ

**名無しさん＠５周年** · 2015/08/23(日) 10:55:50.63

お前相手が一人だと思ってんだな

**名無しさん＠５周年** · 2015/08/23(日) 15:36:36.95

>>188
実際には存在しえない壁面摩擦速度0を仮定してる時点で何でもありだろ

**名無しさん＠５周年** · 2015/08/23(日) 18:26:54.57

無駄に煽るのはみっともないから止めんさい

**名無しさん＠５周年** · 2015/09/04(金) 00:35:06.64

[再]ためしてガッテントイレ問題大解決SP(c)2ch.net
http://nhk.2ch.net/test/read.cgi/livenhk/1441290593/

**名無しさん＠５周年** · 2015/09/05(土) 14:06:01.08

http://i.imgur.com/7aptG1x.jpg
http://i.imgur.com/jBI98ec.jpg
http://i.imgur.com/8q7xsUg.jpg
http://i.imgur.com/KTkwbu8.jpg
http://i.imgur.com/dgw83fk.png
http://i.imgur.com/ewKXFH6.jpg

**名無しさん＠５周年** · 2015/09/11(金) 04:14:12.39

>>190
空を自由に飛びたいな

**P35 ★** · NG

P35.cgi (Perl,SJIS) 2015/11/18 FOX.
>&ft;1 元のスレ
../sim/dat/1394190575.dat
P35 自作・フリーソフト流体シミュレーションスレPart3
311692 -> 53431 (バイト)
P35　テストに夢中。これ以降書いても消えちゃうかも

**名無しさん＠５周年** · 2015/12/12(土) 00:54:28.02

CFD格言　「ケチは損をする」
人間たるもの誰しも，高速に，素晴らしい解析結果を得たいと考えるものである。
しかしそのようなケチな心のせいで，時に大きな損をすることになる。
一時の焦りが粗い格子でのGoサインを出してしまうのである。
格子解像度が足りず境界層が出てこなかった，渦が消えてしまった，
格子生成を妥協した結果として衝撃波が異常に厚く，ギザついてしまった。
このような場合には計算結果は，全く使いものにならない。
一見つまらない解析にもふんだんに格子点数を使える余裕を持った心持ちを普段から維持していきたいものである。

**名無しさん＠５周年** · 2015/12/12(土) 01:05:08.56

CFD格言「ケチは損をする」
人間誰しも，よりはやく計算結果を得たいと考えてしまいがちである。しかし，焦り過ぎは禁物だ。
あなたの作った格子をよく見てみなさい。直交性がガタガタ，格子幅の拡大率は2倍近く，格子点数もかなりケチっちゃいましたね？
境界層速度分布はおかしくなり，渦は消え，衝撃波は異様に厚くなることでしょう。
このようなことにならないよう，リッチに格子点を使える余裕とおおらかな心を持ちたいものである。

**名無しさん＠５周年** · 2015/12/19(土) 21:27:40.92

むしろ限界まで格子粗くしたい

**名無しさん＠５周年** · 2016/03/09(水) 19:31:43.33

すいません　お聞きしたいのですが　テキストのメモ帳を開いたまま

ほかのサイトをクリックしても　メモ帳を表示画面に出しっぱなしにできる固定できる

フリーソフトとかないでしょうか

**名無しさん＠５周年** · 2016/05/19(木) 23:54:21.17

DNS厨兄貴は死んでしまったん？