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

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

1:DNS厨
14/03/07 20:09:35.07 .net
CFD全般、圧縮・非圧縮、単相流・多相流何でも議論しましょう
ただし、OpenFOAMは専用スレかGoogleフォーラムで議論しましょう

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

2:DNS厨
14/03/07 20:10:06.56 .net
　　　　　　　,　- ―‐ -　､
　　　　　／　　　　　　　＼
　　　　 /　　　 ∧　∧　 , 　ヽ
　　　./　ｌ＼:/- ∨ -∨､!　,　',　　さあみんな集まってー！
　　　/ ハ.|/　　　　　　 ∨|,､ﾍ　　　Part3がはじまるよー！！
　 |ヽ' ヽ　　　 ●　　●　　　ﾉ! l
. 〈｢!ヽﾊ._　　　　＿＿　　　　_.lﾉ　|
　　く´ ＼.）　　ヽ. ノ　　　（.ﾉ￣
　　　＼　｀'ｰ-､＿＿_,_ - '´
　　　　｀　- ､ ||V V||　＼
　　　　　　　　| ||　　|| l＼ヽ

前スレ
自作・フリーソフト流体シミュレーションスレ
ｽﾚﾘﾝｸ(sim板)
自作・フリーソフト流体シミュレーションスレPart2
ｽﾚﾘﾝｸ(sim板)l50

3:DNS厨
14/03/07 20:12:10.71 .net
MAC法と乱流モデルの解説
URLﾘﾝｸ(www.caero.mech.tohoku.ac.jp)

有限要素法
URLﾘﾝｸ(www.colorado.edu)
URLﾘﾝｸ(www.colorado.edu)
URLﾘﾝｸ(www.colorado.edu)

LBM
URLﾘﾝｸ(serve.me.nus.edu.sg)
URLﾘﾝｸ(repository.kulib.kyoto-u.ac.jp)

スペクトル法
コロケーション法
URLﾘﾝｸ(citeseerx.ist.psu.edu)
ガラーキン法
URLﾘﾝｸ(www.made.gifu-u.ac.jp)
タウ法
URLﾘﾝｸ(eprints.lib.hokudai.ac.jp)

DNSアルゴリズム
Splitting法
URLﾘﾝｸ(www.dtic.mil)
Influence Matrix法
URLﾘﾝｸ(channelflow.org)
Kim, Moin & Moserの方法
URLﾘﾝｸ(citeseer.ist.psu.edu)

4:DNS厨
14/03/07 20:15:17.43 .net
フリーソフト
DNS：Channelflow
URLﾘﾝｸ(channelflow.org)

VOF
URLﾘﾝｸ(www.lmm.jussieu.fr)

OpenFVM
URLﾘﾝｸ(sourceforge.jp)

いろいろ：OpenFOAM
URLﾘﾝｸ(www.openfoam.com)

5:DNS厨
14/03/07 20:17:17.09 .net
初心者および勉強用リンク
地球惑星物理学演習:時間発展問題の数値解法
URLﾘﾝｸ(www-solid.eps.s.u-tokyo.ac.jp)
物理モデルの基礎と選択法　　乱流モデル
URLﾘﾝｸ(accc.riken.jp)
連立1次方程式の基礎知識　　反復解法の勉強
URLﾘﾝｸ(yebisu.cc.kyushu-u.ac.jp)

計算科学チュートリアル2008　　流体シミュレーションのベンチマーク
URLﾘﾝｸ(www2.cc.kyushu-u.ac.jp)
ベンチマークの一つ　二次元キャビティの論文
URLﾘﾝｸ(www.ann.jussieu.fr)
計算結果が載っているのでコードを組んで計算結果を比較してみよう

6:前スレ441
14/03/08 14:11:39.98 .net
>>1-5
禿乙です
突然前スレなくなっちゃったんでびっくりしてました
新スレでもどうぞよろしくお願いいたします

7:前スレ441
14/03/08 14:27:42.29 .net
おお、書き込めた
連投すみません
自宅のプロバイダが巻き込まれちゃったのでネカフェに出張なのですが
来た甲斐がありましたｗ

お勧めいただいた渦度-流れ関数法のキャビティ流れ、ようやくできました
分割メッシュ数１０×１０ですが（笑）、ちゃんと河村先生のプログラムと
同じ値を出してくれてます

で、調子にのってGHIAさんのデータと比べるべくメッシュ数とRe数上げてみたのですが……
とたんに収束遅くなるし油断すると発散するし……orz

なかなか奥が深いのですねｗ
とりあえず棚橋先生の「はじめて」の安定性のお話とか参考にして
このシンプルな解法でどこまでやれるのか見てみようと思ってます
長文失礼しました

8:DNS厨
14/03/08 21:01:45.96 .net
　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 ,...,ｨ-、
　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 /　j　 |´ﾞi
　　　　　　　　　　　　　＿　　＿　　　　　　　　　　　　　 /　,j　 i　ﾄ、
　　　　　　　　　 , -　´　　　　　　｀　-　iヽ　　　　　 rx,/　 j　 j　 j　|
　　　　　　　　／　　　　　∧　　　　　　　 !-、　　　 j　 ! 　 !　 !　 !　!
　　　　　 /　　　　/!　/　!　|、、　　　　 -ﾞ　　　ﾉ　 ! - 、　　ﾞ　j
　　　　　　 /　　　,__/ | /　　!. | ヽレｨ　　　 i　　 ,.ｲ　　ヽ　｀　、　ﾉ
　　　　　　j 　　 /ﾚ≧i、　　ヽレ≦ヽi、　　　i　 /　!　　 !ヽ　　　　j
　　　　　　i　　i　!,／ｨェﾞ　　　｀ｸケ,｀ヽヾ　　 | /　 i　　　　ヽ　　 j
　　　　　　|,、 iﾚﾞ i　fﾞUﾍ　　　 fﾞiUﾍ　| j j- , /|　　 ! 、　　　　　ノ
　　　　　　　Yｷ ! ﾞ弋夕　　　弋夕　´ ! !ﾐ/ﾞ　|　　ヽ｀　　 ,.-‐ ´!
　　　　　　　 !、!| ///　　　　　　　/// j j /　　{　　　ヽ_,／　　 /
　　　　　　　 | ﾞi !　　　 r----‐‐‐ 、　 jﾉ´　　　ヽ　　　＼　ノ/
　　　　　　　 |　V 、　　!　　　　　ﾉ,. ィ　´　　　ヽ　　　｀´ノ
　　　　　　　 |　　／｀iﾞﾞﾞi｀x―,xァ´ノ /||　　　　　　ヽ、＿_／´
　　　　　　　 | ／　　ﾉ　{,イ‐ﾍ/ く__ i　!i　　　　　　 _, - ´
　　　　　　　 !´ !　　｀ﾌヽ }-{/　／/　 |　　　_, -ｧﾞ´
　　　　　　 / 　 !|　　＼ヽ /／　　　　　　　／!
　　　　　 / 　　 !　　　　＼ﾞ´　　　　　　　 ´i　 |
┌──────────────┐
│>>6 >>7 いつの間にか落ちてたんで、はかせも吃驚したけど　　　.｜
│>とたんに収束遅くなるし油断すると発散するし……orz 　　　　　　.｜
│数値的な安定性、解の収束性、Reの変化による粘性の影響とか.｜
│体験できるので色々計算してみてね　　　　　　　　　　　　　　　　　 .｜
│メッシュ数 10x10、40x40、70x70、100x100、129x129 で計算して　.｜
│Ghia論文と比較すると解の収束が見られて面白いかもしれない　.｜
└──────────────┘

9:DNS厨
14/03/08 21:23:55.32 .net
　　　　　　　　　　　　／　　　　　　　｀丶、
　　　　　　　　　　｣∨　　　　　　　　　ヽ
　　　　　　　　　　7　　　　　　　　　　　　 '.
　　　　　　　　　 /　　　　ノ|　∧ ｀ト　　N|　　っ
.　　　　　　　 /N　　　イ/ |/　 ∨| ∧　　|　　っ
.　　　　　　　〃　Y⌒i　 Ⅳ∩　　　 ∩ 　}/V　　　　-=ﾆ二三三三三
　　　　　　　/ 　　{｛｛ |　 |　∪　　　 ∪　｛{-=ﾆ三／￣￣`'く>三三三
　　　　　　　　　／ゝ _|　 | ''　　　　　　 '' 人.三三|（‐～～‐）| |三三三
　　　　　 /　　厶-=ﾆ | :八‐-(_:::::::⌒)＜∧ﾉ三三|　ー--‐'　|_|三三三
　　 -=ﾆ二三三三三＼_,,,ﾌ三三三三三三三三|　　<・ )　｢三三三
二三三三三三三三三三三三三三三三三三ﾆ圦　　｛__>,,ﾉ三三三
三三三三三三三三三三三三三三三三三三三ﾆ≧===≦三三三
┌──────────────┐
│そういえば、Ghiaの論文で間違い発見？　　　　　　　　　　　　　　 .｜
│398ページのTABLE I、Re=3200、y=0.4531の時　　　　　　　　　　 .｜
│u=-0.86636 じゃなくて u=-0.086636 かもしれない　　　　　　　　　 .｜
│　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│>>7 流れ関数-渦度法である程度確認できたら次は、例えば　　 .｜
│MAC法やSMAC法でコード組むと3次元や他の流れへの適用が .｜
│可能になりますけど　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│↓SMAC法については初心者リンクで貼り忘れたけど　　　　　　 .｜
└──────────────┘
第４回：流体計算の時間進行法と計算アルゴリズム
URLﾘﾝｸ(fluid.web.nitech.ac.jp)

10:DNS厨
14/03/10 23:45:28.32 .net
　　　　　　　　　　　　-‐=‐- ミ
　　　　　　 .　 7ｱ⌒`＿_'⌒｀　、
　　　　　　 .　　　 ,.イ.: :7:. :.｀:: ..､Ｖ)__
　　　　　　ﾞ｀ｰ=ァ´.: /.: :/.: :. :.､:. :. ﾞ'く⌒`
　　　　 / 　／.:: .:/.: :/.: ハ:. :. :.＼:. :＼__)
.　　　　 {_..イ.:/.: .:/.: :/.: /　jﾄ､:. :. :. ＼:. :.〈
.　　　 /.: :l.: .: .: .: .: {／- ､　Χﾊ:. :. :ヽ:. :.
　　　 /.: .: |.:{.: .:.{.:: 仏=ミ　　ｨf〔7:. :. :. }:. :!
.　　　 {.: .: :八:. :. :. :. ｢Ｖ_ｿ　　弋ｿﾊ:.. :.八_{
　　　八: /.: .:Vﾊ:. :.｜``　　　' ｀八! :/
　　 /.: :イ.: .: :: :: :V 个ト .　ﾏﾌ ,.く
.　 /／"j|.: .: :: :: :/〉八　〈〕T爪　 :.
　/ﾞ　　 |.: .: :: ::/`ヽ〉:.＼__八__)　 :i
　　　　　|.:: ::. ::.{　　V〉:. ｢r-＜) ､ |
┌─────────────────┐
│流れ関数-渦度法、単純だけどそれなりに色々計算に使われてる　　　　　　　　 .｜
│みたいじゃないか　　　　　正直侮ってたじぇw　　　　　　　　　　　　　　　　　　　　 .｜
│カルマン渦シミュレーションモデル (Karman)　　　　　　　　　　　　　　　　　　　　　 .｜
│URLﾘﾝｸ(www.gfd-dennou.org) 　　　　　　　　　　 .｜
│渦度・流れ関数法によるマイクロポンプ内の流れの数値計算　　　　　　　　　　　.｜
│URLﾘﾝｸ(www2.nagare.or.jp) 　　　　　　　　　　.｜
│Volume Penalization法による翼果まわりの流れの数値解析　　　　　　　　　　　　.｜
│URLﾘﾝｸ(www.kurims.kyoto-u.ac.jp) .｜
└─────────────────┘

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

12:7
14/03/16 11:28:08.23 .net
>>8-10
おもしろそうな文献ありがとうございます
ふたたびネカフェから失礼しますｗ

プログラムは、取りあえず11×11、Re=100に戻し
そこから順にRe数、メッシュ数を上げて計算してみました

で、驚いたことが二つ
一つ目は11×11､Re=400で突如流速がほとんど0になったこと
41×41ではRe=3200､81×81では5000でおなじ現象が……
これが噂の「数値粘性」でしょうか？
まさかここまで酷いとは……
メッシュサイズは要注意なのですねｗ

二つ目は129×129､Re=7500
流速に時間的な振動が現れ、いっきに主流速がわからなくなるくらい優勢にｗ
最初はただの渦と思い「おお、これが乱流化なのか！？」とか思ってたのですが、
でもよく考えたら、移動壁付近の流速が壁と逆なんてありえなさすぎですｗ
河村先生の本見たらメッシュが粗いと危険とのことで
とりあえずメッシュを257×257にし、今はなんとかRe=10000まで安定して走り、
ちゃんとGHIAさんのと同じような流速分布がでるようになりましたｗ

自宅の古いノーパソでこんなことができるなんてめっちゃ感動です
次はまずパラビューで流線の可視化をやってみて、ついでＭＡＣ法も
書いてみようと思ってます
どうもいろいろとありがとうございました

13:DNS厨
14/03/17 07:13:24.74 .net
　　　　　　　　　　　　　　　　　　　　　／
　　　　　　　　　　　　　　　　　　　　ｰ―''⌒ヽ
　　　　　　　　　　　　　　　　　　　　　　　　 -┴―-　 ___
　　　　　　　　　　　　　　　　　　　　｛ｰ'´: : : : : : : : : : :{　}｀>
　　　　　　　　　　　　　　　　　　　　/ : : /: : |:..:.､|: : ｉ: :ヾ⌒:,　　　うぅ…管理局のウィルス兵器にやられた
　　　　　　　　　　　　　　　　　ー=彡:.: ｲ: : : |: : :.|: : |: : : |: : |
　　　　　　　　　　　　　　　　　＞{:／-{:､: 厶: :..|: : |: : : |: : |　　　コメントは直ってから
　　　　　　　　　　　,:⌒i___ 　　　 /::「￣|リーﾍ::|: : |_: :: |: : |　ﾆ=-
　　　　　　　　　　 {　　___}　　　 |::::{＿,|　艾亥|: : | }:::::|: : |　　　　　ﾆ=-
　　　　　　　　　　　ゝ.ノ　　　　　　ゞ:ﾍ　　｀¨¨¨¨|: :/ｨ{ノヾ:ノ⌒ヽ　　　　　　　　ﾆ=-　γ...ヽ
　　　　　　　　　　　　　　　 ,，　￣｀／i:≧号‐≦{イ:i:i:ﾊ ´:　:　:　: ｀¨¨`丶、　　　　ゝ_,.ｲ
　　　　　　　　　　　　 ,，　´￣￣／､:／¨¨¨`yi:i:i:i:i:i:i:i:ｉﾑ:　:ヾ(:　:　:　:　:　: ｀丶、　　 | 　 |
　　　　　　　　　　__／´￣｀ヽ　ー'　　ヽ／､:i:i:i:i:i:i:i:i:i:i:ﾑ:　:　:　:　:　:　:⌒⌒ :　:..ヽ　 | 　 |
　　　 ,，　--- ´ /: : : (⌒__): : :｀丶、　 (ヽ__＞ ´ ￣｀ヾ}:　:　:　:　:　:　:　:　:　: : ﾉ　 | 　 |
　　／ ) : :(ﾟωﾟ)/: : ,､_,､ : : : : : : : : __:＼　　　　　＼　　＼: : : : : : (　:　:　:　: ／　ヽ .| 　 |
　.::'゜ｰ'゜ : : : : 、 : :..(ﾟωﾟ) : : : :,､_,､(___ ヽ:＼　　　　　￣￣　｀ヾ: : : : ｝　:　:..／　 -‐　 | 　 |
　|: :r～y : : ､ﾉ ): : : : : : : : : : :(ﾟωﾟ): : : : : : : :ー―― 　____　＼￣￣~~　　　　　 _| 　 |
　|: :ー‐′: : :_´__ : : : :_,ﾉ ):_:_:_:_:_:_:_:_:_:_:_(⌒): : : (ﾟωﾟ): ﾉ ) : }　　 {　　　　 -‐　　　　 |
　{⌒ヽ: : : : :ゝ＿): : : : ´: : : _,ﾉ ): : : : :,､_,､ : : : : : : : : ,､_,､: : }　　　｝:-‐ 　　　ﾆ=-　| 　 |

14:名無しさん＠５周年
14/03/18 00:46:09.08 .net
>>11
固形物の変形とか、分裂は考える？
例えば、配管の中でウン◯チが崩壊しながら流れる混相流とか

15:名無しさん＠５周年
14/03/18 01:26:02.98 .net
>>14
固形物は球状粒等単純なもので変形はほとんどしない前提で考えています

16:DNS厨
14/03/18 02:11:48.50 .net
　　　　　　/- '/＼ {￣.＞､＼}　 .＼
　　　　 //}/,　 ,､ヽ　　　　＞､　　‘，
　　　　　i {/ /　/'　,ﾑ　ヽ　　　ヽ　　i
　　　　　Y　 i　 i / {', ＼　',　　　 ∧　 i
　　　　 /　　i　 i{―､ヽ -‐､}　　　　 i. ﾉ　　＿
　　　　 i　i　 i　.{代ﾘ　弋リ　　　　jﾏ　　　ヽヽ　,　-‐≦ ヽ
　　　　 ',. ',　 i　{ ､::　'_ _, :::j　　　/　＼　　ﾉ. '-´__／￣￣
　　　　　ﾏ.{､　　∨＞､ .／ﾉ　　／-､　.＼ {　f二_　〉
　　　　　　ヽ＼ r-､.ﾑ ,ﾋLﾉイ /　　 }　　〉ト ― j
　　　　　　　ムﾉ　Y/¨　 j '/ ､'i　　　 i ＼ i ＼≧彡
　　　　　 .／,' .ム／ﾏｬマ '　 , ミ　　　ヽ　/　　　 /j＞.､
　　　 .／/　 i ハニ -='　／､　∨　　ム'　＼＿ﾉj　　　＞
　r- ´　 ,.{　ヽ　　,イ＞.´　　j　ハ　　r-'　　　.ア＼　　　　　＞
┌───────────────┐
│いやはや、病院で薬貰って少し楽になったよ　　　　　　　　　　　　　　　.｜
│>>11 目的の流れを求めた後、比較実験する事になるだろうから　　　.｜
│先ずは、単相流で単純な円筒パイプ流れ=>曲がり管内流れシミュ　　 .｜
│レーションして既存の研究結果と比較してみようじゃないか　　　　　　　.｜
│曲がり管路内流れ　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│URLﾘﾝｸ(www.cc.utsunomiya-u.ac.jp) 　　　　　.｜
│　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│本題については固体粒子が均一に含まれていると仮定&流れの特性 .｜
│が知りたいなら、ざっくり非ニュートン流体として計算という手がある　 .｜
│非ニュートン流体の数値流動解析-流体内部構造からのアプローチ- .｜
│URLﾘﾝｸ(www.jsme-fed.org) 　　　　　　　　　.｜
│非ニュートン流体特性に着目した砂堆河床上の高粘性流れについて .｜
│URLﾘﾝｸ(library.jsce.or.jp) 　　　　.｜
│粒子の挙動が重要な場合は混相流+移動境界あたりで解くのかな？　.｜
└───────────────┘

17:DNS厨
14/03/18 02:48:46.23 .net
　　　　　　 /__ﾉV(__　_)I*=≦ｱ⌒'くi:i:i:i:i:＼
　　　　　 /⌒)__/⌒＼: : : ｀`丶、__´'くi:ｉ:ｉ:_〉
　　　　 .′/⌒{＼: : :.＼ : : : : : ＼｀`丶、
　　　　　　{ ./.: : :从: :丶: : ∨/: : : : : : : : : :＼
　　　　　　{/.: : :/.{　＼: : ＼∨/ : : : : : : : : : : ＼
　　　　___/.: : :,′{　　＼ : :｀∨＼: : : }}: : : : : : : ＼
　　　 ⌒,′: .:′､　　　＼:__|=- ≧=--‐: : : : :｜:.
　　　　 ′: : ｉ{___ ＼　　　 ´　i|≧=-: :.}}: : : : : : : |: :
　　 ___彡: : : :ｉ{≠ぅミ　　　ｋ抖竺=ミ}}: : : : : : : |: :
　　　　 }: : .:}八{　{/{ﾊ.　　 ⌒ {/{_い刈: : : : : : : |: :
　　　 ,: : : :乂__　Vrｸ　　　　　乂rクﾉ :／.: : : : ; : :
　　　 ′: : : }｝' ' '　　　　　　　　 ' ' 'ア: : : : : :_/_: :
　　　ｉ{: : :‘, 圦　　　　r--‐　vi　　／}.: : : :／　 }/
　　　八: : ∧: : : .　　　　　　ﾉ　　 __〉: : /　　　ん::
　　／　 ):/ : }}:. : : : : .　　　　　　 __〉: : :{　　　 {/ :
　　　　_/: : :.}} : : : : : :}＞=ぉi爪　 __ｱ. : :/} 　　 '{=-‐
　　　 ⌒ｱ: : :}}: : : :j{￣￣ ,ハ ____アﾞ.: :／ ,}i　　　 '{__
　　　 /: __.: :＞‐j{　　__/　 }_アﾞ.: :／__-=圦　 /⌒`
┌──────────────┐
│寝込んでいる間に新AAが増えてて草生えたじぇw　　　　　　　　　　 .｜
│>>12 数値粘性の現象といえば解がなまる、低Re計算結果と同様 .｜
│になる、幻影解が出る等だから一概には結論出せないかなぁ？　　.｜
│時間の刻み幅を十分小さくすると結果が変わるかもしれない　　　　.｜
│　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│他のコード自作するなり、出来合いのアプリ使うなりした時も　　　　.｜
│今回の経験が生きてくると思うので頑張ってね！　　　　　　　　　　　.｜
└──────────────┘
流体現象の数値シミュレーションにおける人工粘性のもつ意味について
URLﾘﾝｸ(www.kurims.kyoto-u.ac.jp)

18:DNS厨
14/03/18 02:58:33.59 .net
　　　　　　　　　　　　　 -=ﾆ　　　　　　　 ¨¨　 _
　　　　　　　　　　／　 ,ｨ仁ﾆ=-rく　　}ﾄ、_____ 八
　　　　　　　　　　　　　{才¨ 　　ﾞ^　　⌒V⌒ヽ二　＼
　　　　　　　 ′　　　　　　　 __　　　　　　　 }二]　　＼
　　　　　　　　」　　　　　-=ｧ7¨ 厂⌒＼_　　　 |ﾆ/
　　　　　_,　" !　／　　／　爪|　　　　∨∧.　ﾚ'
　　　　 "´　　〉　　. :　／}「 ].　　　 ∨∧
　∠,_ 　　　　 /　　 .:　／　 l!　ト、　　　 ∨∧　　　　　　　 }
./⌒　　　　 /　ｨ}/／　　　 l!　|,∧　　　　}/,∧　　　　　　.:::
{　　　　　　./__彡' /⌒丶　　}L_∨∧　　　.|//,∧　　　　 .::::::
　　ー=≦__　　　_}__　　　　 l!.　 ∨∧　　 ]　 /,∧　　　.::::::::|
∨ 　　　　7　　|i, 厂}心　　　　　 ∨∧　　　　//ﾊ　　.::::::::/
. ∨　　　 ]:./　八{　乂j_}　　　ｨ===zﾊ}　　　　 |i　.::::::／
.　∨　　　ﾚ　 |i .「　｀¨ '　　　　 }　r心　　　　　 ,ﾙ}::／
　　∨　　il　　|i从　　　　　　　　弋 _r癶　　　　 / ,ﾘ
.　　 ∨　.il　 .从込_　　 r　 __　　ー==彡　　　 /　八
　　　 ∨ ]!　 }/川/〕ト　乂__｀）ー=…7　＼ /　/　∧
.　　　　Ⅵ| 乂∧///ﾘ　ーr ===≦¨./　　　廴.　 //∧
　　　ー=彡/////,/　,/厂}____彡ヘ/　　　　/　　　//∧
.　　　　　　/////,/　,/ {　 ∧　　∠　　　　/＼　　　//∧
┌─────────────┐
│>>14 流体中の固体の分裂とかってどう計算するの？　　　　.｜
│固体粒子の変形、分裂、合体とか考えるだけで、うっ、頭がw.｜
│一応変形は↓とか見つけてきたじゃないか　　　　　　　　　　 .｜
└─────────────┘
二相系格子ボルツマン法を用いた流体中における粘弾性皮膜固体の挙動解析
URLﾘﾝｸ(www.jstage.jst.go.jp)

19:名無しさん＠５周年
14/03/18 15:21:32.68 .net
圧縮性流体コードをつかって衝撃波菅問題を解いているのですが，片側を真空(n=0,p=0)にしたいのですが，どうしたらよいでしょうか．
そのままでは発散してしまいます．

20:DNS厨
14/03/18 22:53:51.11 .net
　　　　　　＿＿　　　　　,'　厳
　　　 ,.. '"　 ,:　　￣ヽ、､!　密
　　 ,.'　, /. / l　|､　、　ﾞ'.|.　解
　　l ,' iｰｴ'、ヽ;ﾞｰ- l､:.　 |　で
　　l !//l:::::i　　l::::ヽヽ｀i :､!　試
　　lﾍj　､::ﾉ　　､::::ﾘ　| lﾅ|　行
　　〈 "　　　　　¨　" | ﾚﾞ|　し
　　　l＞...._ ,.‐-､　_,,.ｨi ⊂.　た
　　 |　:,イヽ｀ﾆｱ:/iヽ'　　|　？
　　 |./　ゝzヽ会't''　ヽ:i　ヽ.＿／
　　 /　　　　　　　　　ヽ :|
┌───────────┐
│>>19 一方が真空だと衝撃波上流側のマッハ数が　 .｜
│求められないんだけど？　つまり計算出来ません　　.｜
│思考鈍ってるので間違えてたらごめんなさいだけど　.｜
└───────────┘
衝撃波管
URLﾘﾝｸ(www.astro.phys.s.chiba-u.ac.jp)
Riemann Problem (Gas dynamics)　　Matlab Central
URLﾘﾝｸ(www.mathworks.com)
数値流体力学の基礎と宇宙気体力学
URLﾘﾝｸ(lyman.c.u-tokyo.ac.jp)
↑圧縮性流体勉強用

21:名無しさん＠５周年
14/03/18 22:55:11.61 .net
>>16
>>18
ありがとうございます。なんとなく筋が見えてきたのでやってみます

22:名無しさん＠５周年
14/03/19 01:56:08.42 .net
>>20
流体モデルでは真空中に広がっていくような計算は出来ないってことですかね？
そういう問題を解きたいのです．

23:DNS厨
14/03/20 07:00:43.14 .net
　　　　　　　　　　.．ｚ==ﾆﾍ￣　' ､
　　　　　　　　.:　　　´二＞=､ /V//,､
　　　　　　　 / 　／　　　孑z .ﾊヽY ;
　　　　　 /　 /　　　斗γ⌒ヽ　V
　　　　､　i,..イ　　 ,:'入_l{　　ﾘﾊ
　　　　 !ー　　| 　.〃{ﾟ::!ﾜ　　__.z:ｒ!　 !
　　　　 `ー=ｲiﾄ.　::i` ￣　　(ﾟ'ﾊ/i! i !
　　　　　　　　ﾘﾍ i !:､_　 _，r＝1ﾘ　人）
　　　　　　　　.::ﾊト.r!>― <ーイzイゝ-′
　　　　　　　/イ .:ﾘ　>.::∀:.L==ｲ.ﾘ
　　　　　 /.ノ .:/!　込/_Y_|　　 !1
　　　　 ..ノ/ｲｒ-ミ_＿_..｡!!__L__ /.:.､
　≧=z' (　(　 !:::::::..｀.＜::::::..,.´.:::::::イ_ヽ､_
　　　　　￣￣≧=-=:::::..＞ﾟ.:. ＜⌒¨´
　　　　　　　　　　　了三三ﾆ）
　　　　　　　　　　　　￣￣
┌───────────┐
│>>22 ふふん、流体衝撃波とは別物になる悪寒　　.｜
└───────────┘
無衝突衝撃波のミクロおよびマクロ構造と安定性
URLﾘﾝｸ(www.esst.kyushu-u.ac.jp)

24:12
14/03/20 19:34:42.20 .net
>>17
>今回の経験が生きてくると思うので頑張ってね！
はい、ありがとうございます！

時間の刻み幅のことはぜんぜん気づいてなかったので、
さっそく元のdt=10msのほかに、1msと0.1msも計算してみました
収束後（計算開始50秒後）小数点以下６桁まで表示して比べたところ、
この範囲内では流速の値は同じでした

それで、メッシュサイズの違いによる流速分布の変化をもう少し細かく見てみようと思い
追加で２１×２１、１２１×１２１、２４１×２４１を計算し、同じグラフに
プロットしたところ、どうもメッシュが粗いと解が鈍ってくるように思えてきました

　お目汚しの速度分布（恥ずいので数日で消すかもです）
　URLﾘﾝｸ(www.alicenetwork.net)

まだよくわからないのですが、いちおう今の時点では
数値粘性で鈍ってると考えてみようと思ってます

「幻影解」はﾃﾗｵｿﾛｼｽ……ですね
それにはまらないよう、気をつけたいですｗ
どうもありがとうございました

25:DNS厨
14/03/21 00:46:06.30 .net
　　　　　／ｌ, -‐. ´: : : : : : : : ￣: : ..｀. . ､
　　　　＿ｉ : :: : : : : : : : : : : : : : : : : : : : : : ヽ
　　　∠ : : : : : : : : : : : : : ; : : : ﾊ: : : :.: : : : : .ヽ
　　　　/: : : : : : : : : : : :':/ｌ : : / 丶 : : ﾄ; : : : : :丶
　　　 /: :/: :.' : : : : /: .:/　!: ./　　｀.:､: l !: : : : : : ｀、
　　　ｉ: :, : : ,' : : : -A-/-‐ﾋ/´ 　　｀乂_l:_!＿: : ､: i
　　　ｌ: : : : : : : ;ｨ.:/ ,ﾚ〒=､　　　,=〒=､ハ: ,.: .:. :!
　　│: : : : !: :i ﾚ,イ,.::::｀-'!　　　ｌ.,:::｀‐'lヽ V : : : l
　　　l: .:/ニｉ: :l　ｉ |ﾊ.::::.ﾊ}　　　　{ﾊ:::::.ﾊ! }/: : .: : l
　　　|: :{　 (ｌ: :ｌ　 ` ヽ.二.ﾉ　　　　ヽ.二.ﾉ ´ｉ: : .:/V
　　　|: :ヽ、 !: l 　 ``　　　　　　　　　　 ``│ｉV
　　　|: : : : ‐!: ｌ､　　　　 __　　＿_ 　　　 l: l
　　　|: : : : : :V: :＼　　 {　￣ヽ::::}　　　ノ:/
　　　|: : : : : : : : :: : ｀＞ゝ、＿___)/ .-‐.´:V:l
　　　|: : : : : : : : : ;イ　l　｀フ=く_}ヽ: : : : : : :.!
　　　|: : : : :. :: ／　|　ヽ／!:::::::l ｌ　ｌヽ : : : : ｌ
　　　| : : : : ／　　│　　ヽ〉:::〈/　 !　ヽ: : : l
┌───────────┐
│>>24 >プロットしたところ、どうもメッシュが　　　　　.｜
│　　　>粗いと解が鈍ってくるように思えてきました　.｜
│　　　　　　　　　逆に考えるんだ　　　　　　　　　　　 .｜
│　　「格子幅が小さくなると解が微分方程式の　　　.｜
│　　厳密解に近づく-数値解が収束する」　　　　　　.｜
│　　　　　　　　　と考えるんだ　　　　　　　　　　　　 .｜
│　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│お手本になる数値実験お疲れ様でしたけど　　　　.｜
│格子幅を変えて解の収束性を確認するのは　　　 .｜
│やっておくべき計算結果の検証方法の一つです　 .｜
└───────────┘
流れの数値解析における誤差評価と単調収束性
URLﾘﾝｸ(www.nagare.or.jp)
単調に収束しない場合もあるかもしれない

26:24
14/03/26 19:19:40.11 .net
>>25
ありがとうございます！

>逆に考えるんだ　　　　　　　　　　　
涼解です！！

>やっておくべき計算結果の検証方法の一つです
ぐはっ
なんか基本のきの字が抜けてたみたい
「あたしって、ほんとバカ」orz

>単調に収束しない場合も
むーん
「これが解だ」と思ってたら実は違ってた、ってことが
現実にありうる訳なんですね
ｶﾞｸﾌﾞﾙすぎｗｗｗ　気をつけたいです
どうもありがとうございました

27:▼▼▼経済板荒らし～元ナマポ受給者の知能▼▼▼
14/03/27 22:33:02.70 .net
下記スレ量産、２ちゃんねる経済板荒らし、自称坂本陽＠豊中市

「糞汁は経済板から消えるべき」
「産業歴史博物館設立実行委員会（エンジン博物館）」
※スレタイに「北斗の軍」と入ってるスレ立てもこのクソ野郎の仕業

(´-ω-`)y-~（禁煙マスター） ←　史上最悪の中年荒らし　（トリップを変えて別人になりすます技も覚えて来たようですｗ）

こいつのズリネタ　　　　　　　　　　URLﾘﾝｸ(ameblo.jp)

得意技：情報漏洩、ナマポ受給、糞ブログコピペ、ケムトレイル、人工地震、かまってちゃん、驚異の前言撤回ｗ、都合が悪いとすぐ話題脱線

ナマポで助けてもらい　　ｽﾚﾘﾝｸ(occult板:78番)
日本人滅びろ発言　　ｽﾚﾘﾝｸ(eco板:993番)

情報漏洩をきっかけに　　　ｽﾚﾘﾝｸ(eco板:646番)
　　　　　　　　　　　　　　　　　ｽﾚﾘﾝｸ(eco板:926番)
別れの挨拶１！　　　　　ｽﾚﾘﾝｸ(eco板:329番)
別れの挨拶２！　　　　ｽﾚﾘﾝｸ(occult板:303番)
別れの挨拶３！　　　　　ｽﾚﾘﾝｸ(eco板:54番)
そして告発の逆恨み！　　　　ｽﾚﾘﾝｸ(eco板:996番)
　　　　　　　　　　　　　　　　　　　URLﾘﾝｸ(ameblo.jp)
精神異常をきたした　　ｽﾚﾘﾝｸ(eco板:318番)
超低学力！　　　　　　　ｽﾚﾘﾝｸ(eco板:587番)

さらにロリコンも暴露！　ｽﾚﾘﾝｸ(eco板:580-629番)

何を言ってるのか？　　　ｽﾚﾘﾝｸ(eco板:597番)
無駄にスレを作ってる張本人ｗ　　ｽﾚﾘﾝｸ(occult板:256番)

これはひどい　　　　ｽﾚﾘﾝｸ(eco板:512番)

28:名無しさん＠５周年
14/03/29 16:12:06.25 .net
落ちてたのか、気付かなかった。

29:名無しさん＠５周年
14/03/30 05:00:56.18 .net
最近擬似圧縮性法って知ったのですが
これがあまりメジャーじゃないのはなぜなのですか？

30:DNS厨
14/03/30 23:55:47.55 .net
　　　　 '　／: ／: : : /: : ｨ≠＾!ﾐ: :: : : : : ＼__ノ〉
　　　　　,'　/: ://ﾊ: : : //´　　 |!　ヽ: : : : : : : ｲ
　　　　　|.::': ://: : : / }ﾂ　　　　 l!　　＼: : : : :ﾄヽ
　　　　　l:::!://: : : ,' /＿_　　　　|!　　　　 i: : : ::| ヾ!
.　　　　 ',:| /.: : : :|::| ＿__ヽ　　/　　- ―ﾊ: : : :|　!l
　　　、＿Y/: : : : !::衍 :i:iﾃx　/ 　ｘｔ≠ﾁ、}: : :!　!|
　　　　￣７: : : ﾊ个ｊ弋少　　　　 li::ij:ﾂﾉ/: : :/　ﾘ
　　　　　　|ﾊ: : |: :ヾ:l　　　　　　　｀￣ ./: : /、__ノ
　　　　　　|!:ﾚヽ!: : ヽ＞　　、 '__　　　/: //!￣
　　　　　　|!　　ヾ:ヾ个、　　- 　　イ: /ﾊﾉ
　　　　　　し　　　}: ::|ヾ人｀ ‐ｚ‐ '７ﾊ彡ｿ| ﾘ
　　　　　　　　　　|: ｣|: :|辷_7∧／　 |: : ::!￣￣￣ヽ
　　　　　　　　　　|: :|!: :| 　/ ∨__入_.!: : :ヽ、　　　）、
　　　　　　　　 .／!: :ｊ: ::ヽ/、/レ'￣三７ヾ: :ヽ__　　　＼
　　　　　　 .／　 |::/ 八:_::＞/ハ＼　〈　　T ￣ヽ.. 　　＼
　　　　　广ヽ　　i/::/　L ＿ﾊ!ヽｊ　 ≫'!:　　! ＼　＼! 　　＼
　　　　　 |　　　∨::/＼ﾊ　　ん＿ノ ./ |　　| 　ヽ　　',.　　　＼
┌──────────────┐
│>>28 まあ過疎スレだし多少はね？　　　　　　　　　　　　　　　　　　　.｜
│>>29 利点としては、圧力のポアソン方程式を収束させて解く必要が.｜
│無くなる事だろ？　stiffな問題とか圧力が振動するような問題には　.｜
│有効じゃないか？　要は使い所の問題で使ってる所は使ってるって .｜
│感じじゃないかなぁ？　　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│スペクトル法を使ったDNSはポアソン方程式一発で解けるので　　　.｜
│擬似圧縮性法は必要無いな　　　　　　　　　　　　　　　　　　　　　　　 .｜
└──────────────┘
せん断流の数値解析について
URLﾘﾝｸ(airex.tksc.jaxa.jp)
擬似圧縮性解法を用いた二次元円柱渦励振動の数値計算
URLﾘﾝｸ(www2.nagare.or.jp)

31:DNS厨
14/03/31 20:30:07.37 .net
　　　　　　　, ,.ｨ'´　　 ,　　｀ヽ、　　　　　　　　 :: :: :: :: :: :: :: :: :: ::
　　　　　　ｨ'"　　　, / , /.ハ､.ﾑ　　　　　た　　C: :: :: :: :: :: :: :: :::
　　　　　 ,'　　　　,アｴ/i/　二 : l　　　　 .の　　F:: :: :: :: :: :: :: :: :::
　　　　　 l　:　,!.i l i ftｵ　　fｫ.!ilﾘ　　　　し　　D　 :: :: :: :: :: :: :::
　　　　　 |　: l |.l |::: ﾋﾗ　　ﾋﾘ::lﾉ　　　　　い　　っ　　: :: :: :: :: :: :::
　　　　 |　: ｀:i : l　r―--tｯ　ﾉ　.　　　　よ　　て　 : :: :: :: :: :: :::
　　　　 |　.,ｨzi: :.lz.ﾞ..＿_,,..ｨ'l　　　　　　ね　　　　　: :: :: :: :: :: :::
　　　　　 l!/　　ヽﾊﾞｰt'j :i/,' ヽ.　　　　 ! !　　　　 :: :: :: :: :: :: ::.
　　　　　./li.､　　ﾄ.Vﾇfy'｣./　l'　　　　　　　　　　　　: :: :: :: :: :: ::
　　　 ./ :l! ﾞ､　　!<ﾞlﾐ/,｣ l　/!　　　　　　　　　　　　:: :: :: :: :: :: ::
　　　　,' : .|i　ヽ　|　ヽ!'　;l ,' :l　　　　　　　　　　　　:: :: :: :: :: :: ::
┌────────────────┐
│>>26 >なんか基本のきの字が抜けてたみたい　　　　　　　　　　　　　　　　　　.｜
│いえいえ、これが実は難しい問題で世に出てる本を見ても、個々の解法の　 .｜
│素晴らしさを説明して終わり～で具体的な検証はどうすればいいかが　　　　 .｜
│書かれてなかったりするから　　　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│V&Vの資料見ても概念的な説明だしね　　　　　　　　　　　　　　　　　　　　　　 .｜
│シミュレーションの品質保証とV&V　　　　　　　　　　　　　　　　　　　　　　　　　 .｜
│URLﾘﾝｸ(www.jsces.org) .｜
└────────────────┘
コンピュータによる流体力学
URLﾘﾝｸ(www.amazon.co.jp)コンピュータによる流体力学-J-H-ファーツィガー/dp/4621061658/
この本の2章とか具体的に何について気をつければいいか書いてあるけど

32:DNS厨
14/03/31 20:30:56.59 .net
　　　　　 - ─ -
　　-‐ｧ'´　　　　　　　　｀　、
. /　/　　　　　　　　　　　ヽ
.′/　　　　　 /!　∧　ﾄ､､　　ﾍ
|　i　　　　 / / j /　 ! ｊ　! }　　　',
|　l!　　　 / / /ｲ　　|/　ｊ八　}　 }
| 八　　 /}/ N　　　　　 N　}八乂
|　　Y^Y|　　|:|　　　　 |:|　ﾊﾉ
|　　｛（ ||　　U　　　　　 U　 {ﾉ　AA優先で文章を端折った
|　　乂||ｘ:x 　　　　　ｘ:x〉　今も反省してないけど
|　　｜ ||ト．丶. ＿ノ　 ,.ィ7
|　　ｊl＜{　｛｀マﾆﾆフ ´}ト､
|　 /　　＞＼/{　}V ＜　＼
┌───────────────┐
│V&Vとコンピュータによる流体力学両方の視点で書くと　　　　　　　　　　 .｜
│1.数学モデルの選定　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│非圧縮/圧縮、非粘性/粘性、2次元/3次元など解く方程式を選定します .｜
│2.計算モデルの選定と検証　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│離散化手法、計算格子、解法の選定をします　　　　　　　　　　　　　　　　.｜
│理論や解析解との比較やベンチマーク問題を解いてバグ出し、数値的　.｜
│安定性確認、格子幅による離散化誤差の見積もりと解の収束性の確認.｜
│3.計算結果と実験結果の定量的な比較と検証　　　　　　　　　　　　　　　 .｜
│2で数学的な問題を粗方潰しておいて、物理的な問題をここで検証します.｜
│　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 .｜
│って感じかなぁ？詳しくは勉強会やゼミで議論するべきかもしれない　　 .｜
└───────────────┘

33:名無しさん＠５周年
14/04/01 07:33:35.92 .net
法的にいいのかよ、これらってさ。

URLﾘﾝｸ(manta.blog.jp)

URLﾘﾝｸ(sisutore.blog.jp)

34:DNS厨
14/04/01 23:59:44.04 .net
　　　　　　　／　　ﾚ': : : : ―.＝７ｙー、　　ヽ
　　　　　 ./　　.／: ﾌ´: : : : : ／:￣＼　＼　∧
　　　　 '　／: ／: : : /: : ｨ≠＾!ﾐ: :: : : : : ＼__ノ〉
　　　　　,'　/: ://ﾊ: : : //´　　 |!　ヽ: : : : : : : ｲ
　　　　　|.::': ://: : : / }ﾂ　　　　 l!　　＼: : : : :ﾄヽ
　　　　　l:::!://: : : ,' /　　　　　　|!　　　　 i: : : ::| ヾ!
.　　　　 ',:| /.: : : :|::| 二二ゞ′/　ヽ- ―ﾊ: : : :|　!l
　　　、＿Y/: : : : !::衍 :i:iﾃx　/ 　ｘｔ≠ﾁ、}: : :!　!|
　　　　￣７: : : ﾊ个ｊ弋少　　　　 li::ij:ﾂﾉ/: : :/　ﾘ
　　　　　　|ﾊ: : |: :ヾ:l　　　　　　　｀￣ ./: : /、__ノ
　　　　　　|!:ﾚヽ!: : ヽ＞　　　 _'　　　　/: //!￣
　　　　　　|!　　ヾ:ヾ个、　　´='　　　イ: /ﾊﾉ
　　　　　　し　　　}: ::|ヾ人｀ ‐ｚ‐ '７ﾊ彡ｿ| ﾘ
　　　　　　　　　　|: ｣|: :|辷_7∧／　 |: : ::!￣￣￣ヽ
　　　　　　　　　　|: :|!: :| 　/ ∨__入_.!: : :ヽ、　　　）、
　　　　　　　　 .／!: :ｊ: ::ヽ/、/レ'￣三７ヾ: :ヽ__　　　＼
　　　　　　 .／　 |::/ 八:_::＞/ハ＼　〈　　T ￣ヽ.. 　　＼
　　　　　广ヽ　　i/::/　L ＿ﾊ!ヽｊ　 ≫'!:　　! ＼　＼! 　　＼
　　　　　 |　　　∨::/＼ﾊ　　ん＿ノ ./ |　　| 　ヽ　　',.　　　＼
　　　　　 |　　　 | ｊ'　　|＼＼_＿_＿_/ /!　 ,'　　　 ',　 ..i:　　　　　＼
　　　　　 |　　　 | |　　:!.　＼＼　＿_/　! /　　　　 |: : :l　　　　　　　ヽ
┌────────────────┐
│ふふん、4月1日がもう少しで終わるので未確定情報でも書こうじゃないかw.｜
│ある計算で二次元キャビティ流れで圧力の境界条件をノイマン境界条件　 .｜
│で解いたら128x128でもRe=1000でGhiaの論文の値からかなりずれたじぇ　.｜
│　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│Web上だとMac法、圧力をノイマン境界条件で解いてReが大きくなると　　　.｜
│流れ関数-渦度法と違ってくるとか報告してる資料があるが本当かもねw　 .｜
└────────────────┘

35:26
14/04/02 19:04:50.82 .net
>>31-32
㌧㌧㌧クスです！

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

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

36:名無しさん＠５周年
14/04/02 19:17:03.95 .net
多額な借金で炎上中。

URLﾘﾝｸ(manta.blog.jp)

URLﾘﾝｸ(sisutore.blog.jp)

URLﾘﾝｸ(livedeaichat.blog.fc2.com)

37:DNS厨
14/04/02 21:34:11.82 .net
　　　　　　　　　　　　　　＿_
　　　　　　　　　　　　 ,.　´　　　｀ヽ.
　　　　　　／￣｀ヽ　/　　, ‐ 、　, ‐'.
　　　　　　　/,　⌒ヽ　i ｜　　｀⌒　　⌒ｉ
　　　　　 { {　､__ノ ,ﾉ　|^v=｢￣｀H´￣}
　　　　　　い､　__ノ}:ｰｧ!　　ヽ＿_ハ__,ｲ‐.､
　　　　　　　,ゝ､　　,ｲ.:｜{`┐ ,. -‐～ｰ〈:.:. :.＼
　　　　　／. .:.:.｀ｉ＿}.: :ヽ｀ｰＶ : : :, へ.:.::',:.:. :. :.ヽ,
　　　／ . : .:.:./´ ＿}ヽ : ｀!7 : ／.:.::::,ﾊ.::}.:.:',:.ｌ.:. :.',
　　 / :, -‐ァ´| 　＿｣丿: : !{／ｉヽ.:.::/: :ｌ/.:.:.:l/.:. : :′
　　/ :/.:.:.:/　/ 、 ―j.:. .: .:i―┼‐∨.:. .:. .:.:.:l:. : : : :ｉ
　　|: : .:.:.:.i　｜　,ｌ￣l:.:. .: .:|＿_|＿_|:. :. :. .:.:.:|:. : : : :|
　　|: : .:.:.:.ｌ　ｌ／ | 　|.:. .: .:| 　 | 　 |:. :. :. .:.:.:|:. : : .:/
　　＼: .:.:.:＼＼,.| 　l.:. .: .:!＿_!_＿!:. :. :. .:.:.:!: : ／
┌───────────┐
│>>35　>具他例とかもあったらわかりやすいかも…….｜
│　　　　>２あたりは見当も付かないですorz　　　　　　.｜
│ほほっ、実はお前さんは既に実践しとるんじゃよ　　 .｜
│ベンチマーク問題=二次元キャビティ+Ghia論文　　　 .｜
│格子幅による解の収束性の確認=>>24の結果　　　 .｜
│　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 .｜
│今回の題材は>>32の1、2の練習にもってこい　　　　.｜
│　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 .｜
│孫悟空とクリリンがアルバイトをやってるうちに強く　 .｜
│なっていた亀仙流修行術みたいなもんじゃな　　　　　.｜
└───────────┘

38:名無しさん＠５周年
14/04/02 21:35:45.68 .net
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。

39:DNS厨
14/04/02 23:07:29.72 .net
､___ﾉヽ.＿,､_ﾉ､__,　　　　　　　　　　　　＿＿ﾊ＿＿
）　　　　　　　　ゝ,　　　　　　　　　　ヽ　　　　　　　/
ヽ　　　　　　　　 / ,,....:'"￣:＼　　　　ﾉ　　も　>>27ヽ._
/　|　　　　-|-iヽ〉.:.:.:.:.:.:.:.:.:.:.:.:.＼　　）　 .う　>>33　/
〉　ゝ‐'　　 /_ﾉ　（:;.-､:.:.:.:.:.:.:.:.:.:.:.:.:＼ヽ　　っ >>36 ヽ
､　ーｧ　　 /　　〉´ﾉ.:.:.:.:.;:ｨ'"｀ヽ、.:.:ｿ　　 ! ! >>38 （
ﾉ　（__　　 /^Lﾉ　､:.:.;.ｨ''､ /|　ﾊ .＼:.)　　　　　　　/
ヽ　ー､　　/＼・・（', /l Xヽl. /,__Ｘli ヽ}-､,.--v'⌒
〈　　‐'ﾞ　 /　　＼（./ /　 >l/　＜｀ lﾉ.lヽ.:l
　〉 | | |　　　/　　　〉〇‐'＿＿__ .ヽOﾉｲ/
〈　・.・.・　　/^Lﾉ　{ヽ_〈　　　｀ヾ! ,,.ｨ':／
. )　　　　　　　　　rﾞ'.:.:.｀''''===ｨ'.:.:／
　},..-､,.---､,..--､fﾞ'.:.:.:.:.:fi.:.:.:.<､了
　　　　　　/.:.:.:.:.:.:.:.:.:.:.:fjﾆ).:.:.:.:ヽ'ｿ
　　　　　./.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.､
┌───────────┐
│>>35　>シミュレーションの結果を他の人に伝える時 .｜
│　　　　>にも超重要っぽいですね　　　　　　　　　　　.｜
│そうですね　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│企業の技術者が会議でフルボッコされないための　 .｜
│方法論的な側面も強い気がするかもしれないw　　　 .｜
│　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 .｜
│>ファーツィガーさんの本をぽちってみました　　　　　.｜
│主に取り扱っているのが差分法と有限体積法です　 .｜
│3、4、6章で扱っている問題のソースコードがftpで　　.｜
│落とせるので無駄な買い物にはならないと思います　.｜
└───────────┘

40:名無しさん＠５周年
14/04/08 23:39:45.37 .net
ド素人が、実機流体の温度分布や速度分布、あるいは配管の熱応力分布を
可視化して表現するためにはどのようなプロセスを踏んだら宜しいか？

41:DNS厨
14/04/09 06:54:28.45 .net
　　　　　　　　　　　 -─‐ 厶　｝}、
　　　　　　　_　--　´　　　　　　　　｀ヽヘ
　　　　　｛｛｀7　　　 |　　　＼　　　　 ∧
　　　　　　　∨　'/　/│∧　{ ⌒＼　ﾘ r-‐-‐､
　　　　　　　,'　　　/⌒|/　＼x‐=ミ∨〈/ / /,ハ_
　　　　　　 |　i　i/ｘ=ミ.　　　 rしﾊ }!|　ゝつく　ﾉ∧　　温度境界層だにゃん
　　　　　　 _|　i　| i{ {しﾍ　　　Ｖ_ソ　|　|ハ￣　　 '､
　　　 _r'｢｛∨|i∧　 Vソ　　　　　''　|　|ノﾊ＿＿_ノﾍ
　　　｛_{└と从ﾊ　''　　　　,　　 |　　　|'.　　　　∧
　　　　| `┴''　| |ｉ八　　　ー　　　イ|/ー￢ﾍ　　　 -ﾍ
　　　　l.　　　 | |＼|｀　ーｚr‐≦／ │　　　 |＼
　　　ﾄ-―　ｉ_Ｌ. イ⌒´/ ∨:::∨　｣　　　　　　　　　　｝
　　　　　　　人　　　　　〉　〉〈　　 >　　　　　　＿＿ノ
　　　｛　　　　　＼ ∨　＼　∨ ／　　　　　ア｢
　　　　＼,,＿＿＿_､|　　　　｝　　　　　 /|│　　　　／⌒)
┌─────────────┐
│>>40　解きたい問題が、自然対流なのか、強制対流なのか?.｜
│流体の圧縮性は無視できるか?（数学モデルの選定）　　　　 .｜
│先行する研究があるかどうか?（計算モデルとベンチマーク） .｜
│あたりから調べましょう　　　　　　　　　　　　　　　　　　　　　　 .｜
│　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 .｜
│強制対流熱伝達の場合は速度場と温度場が独立なので　　.｜
│まずは温度場を無視して速度場がちゃんと解けるか調べ　　.｜
│ればいいと思うけど　　　　　　　　　　　　　　　　　　　　　　　　 .｜
└─────────────┘
平行平板間流れの層流熱伝達（理論解析と数値計算）
URLﾘﾝｸ(fluid.web.nitech.ac.jp)
プログラム付き強制対流熱伝達の資料です
これを読んで全然分からないなら勉強が必要かもしれない

42:DNS厨
14/04/09 23:55:09.65 .net
　　　　　 /⌒)__/⌒＼: : : ｀`丶、__´'くi:ｉ:ｉ:_〉
　　　　 .′/⌒{＼: : :.＼ : : : : : ＼｀`丶、
　　　　　　{ ./.: : :从: :丶: : ∨/: : : : : : : : : :＼
　　　　　　{/.: : :/.{　＼: : ＼∨/ : : : : : : : : : : ＼
　　　　___/.: : :,′{　　＼ : :｀∨＼: : : }}: : : : : : : ＼
　　　 ⌒,′: .:′､　　　＼:__|=- ≧=--‐: : : : :｜:
　　　　 ′: : ｉ{___ ＼　　　 ´　i|≧=-: :.}}: : : : : : : |:
　　 ___彡: : : :ｉ{≠ぅミ　　　ｋ抖竺=ミ}}: : : : : : : |:
　　　　 }: : .:}八{　{/{ﾊ.　　 ⌒ {/{_い刈: : : : : : : |:
　　　 ,: : : :乂__　Vrｸ　　　　　乂rクﾉ :／.: : : : ; :
　　　 ′: : : }｝' ' '　　　　　　　　 ' ' 'ア: : : : : :_/_:
　　　ｉ{: : :‘, 圦　　　　r--‐　vi　　／}.: : : :／　 }
　　　八: : ∧: : : .　　　　　　ﾉ　　 __〉: : /　　　ん
　　／　 ):/ : }}:. : : : : .　　　　　　 __〉: : :{　　　 {/
　　　　_/: : :.}} : : : : : :}＞=ぉi爪　 __ｱ. : :/} 　　 '{=
　　　 ⌒ｱ: : :}}: : : :j{￣￣ ,ハ ____アﾞ.: :／ ,}i　　　 '{_
　　　　 /: __.: :＞‐j{　　__/　 }_アﾞ.: :／__-=圦　 /⌒
　　　 /／___＼　j{ 　 .｢i:i:ヽアﾟ.: :／:ｉ:i:i:i:ｉ} ___} _/
　　　/　　￣￣} 　 _,ﾉi:i:i:/: : : /:i:i:i:i:i:i:i:ｉ|__ｱ_/　　　/
┌─────────────┐
│ふふん、どんどん逝ってみよう～　　　　　　　　　　　　　　　　 .｜
│流れが層流の場合、管内層流熱伝達のヌセルト数が解析　.｜
│的に求められるので計算モデルの検証に使えそうだな　　　.｜
│　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│乱流熱伝達は結構難しいと思う　　データをどうまとめれば　.｜
│いいかも良く考えてシミュレーションを行おうじゃないか　　　 .｜
└─────────────┘
熱交換器の伝熱評価
URLﾘﾝｸ(www.yama.tus.ac.jp)
局所加熱された円管内の強制対流乱流熱伝達管壁内熱伝導を考慮した数値解析
URLﾘﾝｸ(www.jstage.jst.go.jp)

43:DNS厨
14/04/13 19:54:15.68 .net
　　　　　　　　　　　　 .　 ´　￣　｀　　 .　　　　　　　　　 | j i/〉、
　　　　　　　　　　　　　／:::::::: /!::ﾄ、::::::::::::::＼　　　　　　　 |∩ 　/
　　　　　　　　　 /::::: ;ｨﾅ/ |::lxﾄ､‐ ､:::::::::: ヽ　　　　　　　　 |j　）ﾉ
　　　　　　　　　　　　　ｲi:l::ﾙ笊人l ﾉ笊ﾊi:::::::::::::: ヽ　　　／｀V　イ
　　　　　　　　　　　　　从i:| i.r' ｌi　　l.r' lﾄ |::::::::::::::::. ヽ　　／: ＼　Ｘﾉ
　　　　　　　　　　　　 /:::l:| |匕ﾘ　 l匕ﾘ .|:::ｌi::ｌ::ヽ::::ハ／: : : : : :＼／
　　 rf「ﾄ、　　　　　 /:::::l〈"" 　　 _ ""　|:::ｌi::|ﾙﾊ/／: : : : : : : :／
　　 |　　 i 　＿＿＿r‐ルﾙ＞　.「　　ﾉ .イﾚハj: ﾘ／: : : : : : : :／
　〈.ヽ/´￣/: : : : : : :| l|　　|l |〈　{_て .フ´　　/::/: : : : : : : : ／
　　ｰｌ　　 l: : : : : : : :| l|　　|l |ﾄﾄ、　xへ　／: : : : : : :: : :／
　　人＿人: : : : : : :! l|＿_|l lハ: ≧≠≦＝虫: : : : : :イ
　　　　　￣￣￣￣〈∧___|l |7∧: : : : : : : : : : : : : : : : ﾄ、
　　　　　　　　く＼∧fi|l // ∧: : : : : : : : : : : : : : : |: ＼＿
┌────────────┐
│>>40　どの程度の経験があるか分からなかったので　 .｜
│簡単な流体シミュレーション経験があるとして書いて　　.｜
│みました　　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　.｜
│シミュレーション経験無しなら二次元キャビティ計算　　 .｜
│あたりから始めてみては？　　　　Festina lente　　　　　.｜
└────────────┘

次ページ