TOP カテ一覧 スレ一覧 100〜終まで 2ch元 削除依頼
nAsAnASAメーデー
【最適化】Isight/OPTIMUS/modeFRONTIER
nAsAnASAメーデー
☆有限要素法(Finite Element Method)★
☆有限要素法(Finite Element Method)★
SNS狂鬼ムラカミ
【最適化】Isight/OPTIMUS/modeFRONTIER
☆有限要素法(Finite Element Method)★
【最適化】Isight/OPTIMUS/modeFRONTIER
偏った標本調査・世論調査・統計を吊し上げるスレ
自作・フリーソフト流体シミュレーションスレPart3
- 1 :2014/03/07 〜 最終レス :2018/08/08
- CFD全般、圧縮・非圧縮、単相流・多相流何でも議論しましょう
ただし、OpenFOAMは専用スレかGoogleフォーラムで議論しましょう
役に立つリンクは>>3から
- 2 :
- , - ―‐ - 、
/ \
/ ∧ ∧ , ヽ
./ l\:/- ∨ -∨、! , ', さあみんな集まってー!
/ ハ.|/ ∨|,、ヘ Part3がはじまるよー!!
|ヽ' ヽ ● ● ノ! l
. 〈「!ヽハ._ __ _.lノ |
く´ \.) ヽ. ノ (.ノ  ̄
\ `'ー-、 ___,_ - '´
` - 、 ||V V|| \
| || || l\ ヽ
前スレ
自作・フリーソフト流体シミュレーションスレ
ttp://kamome.2ch.sc/test/read.cgi/sim/1303299935
自作・フリーソフト流体シミュレーションスレPart2
ttp://ikura.2ch.sc/test/read.cgi/sim/1335112085/l50
- 3 :
- MAC法と乱流モデルの解説
ttp://www.caero.mech.tohoku.ac.jp/publicData/Daiguji/
有限要素法
ttp://www.colorado.edu/engineering/cas/courses.d/IFEM.d/
ttp://www.colorado.edu/engineering/CAS/courses.d/AFEM.d/
ttp://www.colorado.edu/engineering/CAS/courses.d/NFEM.d/
LBM
ttp://serve.me.nus.edu.sg/shuchang/Teaching/ME5361.pdf
ttp://repository.kulib.kyoto-u.ac.jp/dspace/bitstream/2433/97133/1/KJ00004711538.pdf
スペクトル法
コロケーション法
ttp://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.36.98&rank=1
ガラーキン法
ttp://www.made.gifu-u.ac.jp/~tanaka/LectureNote/numerical_analysis.pdf
タウ法
ttp://eprints.lib.hokudai.ac.jp/dspace/bitstream/2115/20118/1/29.%E9%BB%92%E7%94%B0%E6%98%8E%E6%85%88.pdf
DNSアルゴリズム
Splitting法
ttp://www.dtic.mil/cgi-bin/GetTRDoc?Location=U2&doc=GetTRDoc.pdf&AD=ADA205883
Influence Matrix法
ttp://channelflow.org/dokuwiki/lib/exe/fetch.php?media=docs:chflowguide.pdf
Kim, Moin & Moserの方法
http://citeseer.ist.psu.edu/viewdoc/summary?doi=10.1.1.147.3090
- 4 :
- フリーソフト
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/
- 5 :
- 初心者および勉強用リンク
地球惑星物理学演習:時間発展問題の数値解法
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
計算結果が載っているのでコードを組んで計算結果を比較してみよう
- 6 :
- >>1-5
禿乙です
突然前スレなくなっちゃったんでびっくりしてました
新スレでもどうぞよろしくお願いいたします
- 7 :
- おお、書き込めた
連投すみません
自宅のプロバイダが巻き込まれちゃったのでネカフェに出張なのですが
来た甲斐がありましたw
お勧めいただいた渦度-流れ関数法のキャビティ流れ、ようやくできました
分割メッシュ数10×10ですが(笑)、ちゃんと河村先生のプログラムと
同じ値を出してくれてます
で、調子にのってGHIAさんのデータと比べるべくメッシュ数とRe数上げてみたのですが……
とたんに収束遅くなるし油断すると発散するし……orz
なかなか奥が深いのですねw
とりあえず棚橋先生の「はじめて」の安定性のお話とか参考にして
このシンプルな解法でどこまでやれるのか見てみようと思ってます
長文失礼しました
- 8 :
- ,...,ィ-、
/ 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 :
- / `丶、
」∨ ヽ
7 '.
/ ノ| ∧ `ト N| っ
. /N イ/ |/ ∨| ∧ | っ
. 〃 Y⌒i W∩ ∩ }/V -=ニ二三三三三
/ {{{ | | ∪ ∪ {{-=ニ三/ ̄ ̄`'く>三三三
/ゝ _| | '' '' 人.三三|(‐〜〜‐)| |三三三
/ 厶-=ニ | :八‐-(_:::::::⌒)<∧ノ三三| ー--‐' |_|三三三
-=ニ二三三三三\_,,,フ三三三三三三三三| <・ ) 「 三三三
二三三三三三三三三三三三三三三三三三ニ圦 {__>,,ノ三三三
三三三三三三三三三三三三三三三三三三三ニ≧===≦三三三
┌───────────────────────────┐
│そういえば、Ghiaの論文で間違い発見? .|
│398ページのTABLE I、Re=3200、y=0.4531の時 .|
│u=-0.86636 じゃなくて u=-0.086636 かもしれない .|
│ .|
│>>7 流れ関数-渦度法である程度確認できたら次は、例えば .|
│MAC法やSMAC法でコード組むと3次元や他の流れへの適用が .|
│可能になりますけど .|
│↓SMAC法については初心者リンクで貼り忘れたけど .|
└───────────────────────────┘
第4回:流体計算の時間進行法と計算アルゴリズム
ttp://fluid.web.nitech.ac.jp/Gotoh_Home_page/Edu/Public_course/Text/Text_4th.pdf
- 10 :
- -‐=‐- ミ
. 7ア⌒`__'⌒` 、
. ,.イ.: :7:. :.`:: ..、 V)__
゙`ー=ァ´.: /.: :/.: :. :.、:. :. ゙'く⌒`
/ /.:: .:/.: :/.: ハ:. :. :.\:. :\__)
. {_..イ.:/.: .:/.: :/.: / jト、:. :. :. \:. :.〈
. /.: :l.: .: .: .: .: {/- 、 Χハ:. :. :ヽ:. :.
/.: .: |.:{.: .:.{.:: 仏=ミ ィf〔7:. :. :. }:. :!
. {.: .: :八:. :. :. :. 「V_ソ 弋ソ ハ:.. :.八_{
八: /.: .:Vハ:. :.|`` ' `八! :/
/.: :イ.: .: :: :: :V 个ト . マフ ,.く
. //"j|.: .: :: :: :/〉八 〈〕T爪 :.
/゙ |.: .: :: ::/`ヽ〉:.\__八__) :i
|.:: ::. ::.{ V〉:. 「r-<) 、 |
┌─────────────────────────────────┐
│流れ関数-渦度法、単純だけどそれなりに色々計算に使われてる .|
│みたいじゃないか 正直侮ってたじぇw .|
│カルマン渦シミュレーションモデル (Karman) .|
│ttp://www.gfd-dennou.org/arch/davis/stpk/demo/Karman/ .|
│渦度・流れ関数法によるマイクロポンプ内の流れの数値計算 .|
│ttp://www2.nagare.or.jp/jscfd/cfds15/papers/B04/B04-3.pdf .|
│Volume Penalization法による翼果まわりの流れの数値解析 .|
│ttp://www.kurims.kyoto-u.ac.jp/~kyodo/kokyuroku/contents/pdf/1808-10.pdf .|
└─────────────────────────────────┘
- 11 :
- 質問です
最終的に水道管をまげたときに内部を流れる固形物を含む液体の流れの
変化の解析を行いたいと考え、CFDに手を出したいのですが、
この分野の事前の知識が全くない状態で、どこから手をつけたらいいかわかりません
理論の習熟と、ソフトの使い方は最低限必要だと思うのですが、何を勉強し、何を使えばいいのかもわかりません。
どなたか、おおざっぱにどういう順番で勉強していけば1行目のことができるようになるか、教えていただけませんか?
- 12 :
- >>8-10
おもしろそうな文献ありがとうございます
ふたたびネカフェから失礼しますw
プログラムは、取りあえず11×11、Re=100に戻し
そこから順にRe数、メッシュ数を上げて計算してみました
で、驚いたことが二つ
一つ目は11×11、Re=400で突如流速がほとんど0になったこと
41×41ではRe=3200、81×81では5000でおなじ現象が……
これが噂の「数値粘性」でしょうか?
まさかここまで酷いとは……
メッシュサイズは要注意なのですねw
二つ目は129×129、Re=7500
流速に時間的な振動が現れ、いっきに主流速がわからなくなるくらい優勢にw
最初はただの渦と思い「おお、これが乱流化なのか!?」とか思ってたのですが、
でもよく考えたら、移動壁付近の流速が壁と逆なんてありえなさすぎですw
河村先生の本見たらメッシュが粗いと危険とのことで
とりあえずメッシュを257×257にし、今はなんとかRe=10000まで安定して走り、
ちゃんとGHIAさんのと同じような流速分布がでるようになりましたw
自宅の古いノーパソでこんなことができるなんてめっちゃ感動です
次はまずパラビューで流線の可視化をやってみて、ついでMAC法も
書いてみようと思ってます
どうもいろいろとありがとうございました
- 13 :
- /
ー―''⌒ヽ
-┴―- ___
{ー'´: : : : : : : : : : :{ }`>
/ : : /: : |:..:.、|: : i: :ヾ⌒:, うぅ…管理局のウィルス兵器にやられた
ー=彡:.: イ: : : |: : :.|: : |: : : |: : |
>{:/-{:、: 厶: :..|: : |: : : |: : | コメントは直ってから
,:⌒i___ /::「 ̄|リーヘ::|: : |_: :: |: : | ニ=-
{ ___} |::::{_,| 艾亥|: : | }:::::|: : | ニ=-
ゝ.ノ ゞ:ヘ `¨¨¨¨|: :/ィ{ノヾ:ノ⌒ヽ ニ=- γ...ヽ
,,  ̄`/i:≧号‐≦{イ:i:i:ハ ´: : : : `¨¨`丶、 ゝ_,.イ
,, ´ ̄ ̄/、:/¨¨¨`yi:i:i:i:i:i:i:i:iム: :ヾ(: : : : : : ` 丶 、 | |
__/´ ̄ `ヽ ー' ヽ/、:i:i:i:i:i:i:i:i:i:i:ム: : : : : : :⌒⌒ : :..ヽ | |
,, --- ´ /: : : (⌒__): : :` 丶、 (ヽ__> ´  ̄ `ヾ}: : : : : : : : : : : ノ | |
/ ) : :(゚ω゚)/: : ,、_,、 : : : : : : : : __:\ \ \: : : : : : ( : : : : / ヽ .| |
.::'゜ー'゜ : : : : 、 : :..(゚ω゚) : : : :,、_,、(___ ヽ:\  ̄ ̄ `ヾ: : : : } : :../ -‐ | |
|: :r〜y : : 、ノ ): : : : : : : : : : :(゚ω゚): : : : : : : :ー――― ____ \ ̄ ̄~~ _| |
|: :ー‐′: : :_´__ : : : :_,ノ ):_:_:_:_:_:_:_:_:_:_:_(⌒): : : (゚ω゚): ノ ) : } { -‐ |
{⌒ヽ: : : : :ゝ _): : : : ´: : : _,ノ ): : : : :,、_,、 : : : : : : : : ,、_,、: : } }:-‐ ニ=- | |
- 14 :
- >>11
固形物の変形とか、分裂は考える?
例えば、配管の中でウン◯チが崩壊しながら流れる混相流とか
- 15 :
- >>14
固形物は球状粒等単純なもので変形はほとんどしない前提で考えています
- 16 :
- /- '/\ { ̄.> 、\} .\
//}/, ,、ヽ >、 ‘,
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 目的の流れを求めた後、比較実験する事になるだろうから .|
│先ずは、単相流で単純な円筒パイプ流れ=>曲がり管内流れシミュ .|
│レーションして既存の研究結果と比較してみようじゃないか .|
│曲がり管路内流れ .|
│ttp://www.cc.utsunomiya-u.ac.jp/~sugiyama/Paper/a6_0.pdf .|
│ .|
│本題については固体粒子が均一に含まれていると仮定&流れの特性 .|
│が知りたいなら、ざっくり非ニュートン流体として計算という手がある .|
│非ニュートン流体の数値流動解析-流体内部構造からのアプローチ- .|
│ttp://www.jsme-fed.org/newsletters/2005_9/no5.html .|
│非ニュートン流体特性に着目した砂堆河床上の高粘性流れについて .|
│ttp://library.jsce.or.jp/jsce/open/00561/2005/7-12_7012f.pdf .|
│粒子の挙動が重要な場合は混相流+移動境界あたりで解くのかな? .|
└─────────────────────────────┘
- 17 :
- /__ノV(__ _)I*=≦ア⌒'くi:i:i:i:i:\
/⌒)__/⌒\: : : ``丶、__´'くi:i:i:_〉
.′/⌒{\: : :.\ : : : : : \``丶、
{ ./.: : :从: :丶: : ∨/: : : : : : : : : :\
{/.: : :/.{ \: : \∨/ : : : : : : : : : : \
___/.: : :,′{ \ : :`∨\: : : }}: : : : : : : \
⌒,′: .:′、 \:__|=- ≧=--‐: : : : :|:.
′: : i{___ \ ´ i|≧=-: :.}}: : : : : : : |: :
___彡: : : :i{≠ぅミ k抖竺=ミ}}: : : : : : : |: :
}: : .:}八{ {/{ハ. ⌒ {/{_い刈: : : : : : : |: :
,: : : :乂__ Vrク 乂rク ノ :/.: : : : ; : :
′: : : }}' ' ' ' ' 'ア: : : : : :_/_: :
i{: : :‘, 圦 r--‐ vi /}.: : : :/ }/
八: : ∧: : : . ノ __〉: : / ん::
/ ):/ : }}:. : : : : . __〉: : :{ {/ :
_/: : :.}} : : : : : :}>=ぉi爪 __ア. : :/} '{=-‐
⌒ア: : :}}: : : :j{ ̄ ̄ ,ハ ____ア゙.: :/ ,}i '{__
/: __.: :>‐j{ __/ }_ア゙.: :/__-=圦 /⌒`
┌────────────────────────────┐
│寝込んでいる間に新AAが増えてて草生えたじぇw .|
│>>12 数値粘性の現象といえば解がなまる、低Re計算結果と同様 .|
│になる、幻影解が出る等だから一概には結論出せないかなぁ? .|
│時間の刻み幅を十分小さくすると結果が変わるかもしれない .|
│ .|
│他のコード自作するなり、出来合いのアプリ使うなりした時も .|
│今回の経験が生きてくると思うので頑張ってね! .|
└────────────────────────────┘
流体現象の数値シミュレーションにおける人工粘性のもつ意味について
ttp://www.kurims.kyoto-u.ac.jp/~kyodo/kokyuroku/contents/pdf/1288-1.pdf
- 18 :
- -=ニ ¨¨ _
/ ,ィ仁ニ=-rく }ト、_____ 八
{才¨ ゙^ ⌒V⌒ヽ二 \
′ __ }二] \
」 -=ァ7¨ 厂⌒\_ |ニ/
_, " ! / / 爪| ∨∧. レ'
"´ 〉 . : /}「 ]. ∨∧
∠,_ / .: / l! ト、 ∨∧ }
./⌒ / ィ}// l! |,∧ }/,∧ .:::
{ ./__彡' /⌒丶 }L_∨∧ .|//,∧ .::::::
ー=≦__ _}__ l!. ∨∧ ] /,∧ .::::::::|
∨ 7 |i, 厂}心 ∨∧ //ハ .::::::::/
. ∨ ]:./ 八{ 乂j_} ィ===zハ} |i .::::::/
. ∨ レ |i .「 `¨ ' } r心 ,ル}::/
∨ il |i从 弋 _r癶 / ,リ
. ∨ .il .从込_ r __ ー==彡 / 八
∨ ]! }/川/〕ト 乂__`) ー=…7 \ / / ∧
. Y| 乂∧///リ ーr ===≦¨./ 廴. //∧
ー=彡/////,/ ,/厂}____彡ヘ/ / //∧
. /////,/ ,/ { ∧ ∠ /\ //∧
┌─────────────────────────┐
│>>14 流体中の固体の分裂とかってどう計算するの? .|
│固体粒子の変形、分裂、合体とか考えるだけで、うっ、頭がw.|
│一応変形は↓とか見つけてきたじゃないか .|
└─────────────────────────┘
二相系格子ボルツマン法を用いた流体中における粘弾性皮膜固体の挙動解析
ttps://www.jstage.jst.go.jp/article/kikaib/77/773/77_773_4/_pdf
- 19 :
- 圧縮性流体コードをつかって衝撃波菅問題を解いているのですが,片側を真空(n=0,p=0)にしたいのですが,どうしたらよいでしょうか.
そのままでは発散してしまいます.
- 20 :
- __ ,' 厳
,.. '" ,:  ̄ ヽ、、! 密
,.' , /. / l |、 、 ゙'.|. 解
l ,' iーエ'、 ヽ;゙ー- l、:. | で
l !//l:::::i l::::ヽヽ`i :、! 試
lヘj 、::ノ 、::::リ | lナ| 行
〈 " ¨ " | レ゙| し
l>...._ ,.‐-、 _,,.ィi ⊂. た
| :,イヽ`ニア:/iヽ' | ?
|./ ゝzヽ会't'' ヽ:i ヽ._/
/ ヽ :|
┌──────────────────────┐
│>>19 一方が真空だと衝撃波上流側のマッハ数が .|
│求められないんだけど? つまり計算出来ません .|
│思考鈍ってるので間違えてたらごめんなさいだけど .|
└──────────────────────┘
衝撃波管
ttp://www.astro.phys.s.chiba-u.ac.jp/netlab/cans/cans1d/md_shktb/Readme.pdf
Riemann Problem (Gas dynamics) Matlab Central
ttp://www.mathworks.com/matlabcentral/fileexchange/45666-riemann-problem-gas-dynamics
数値流体力学の基礎と宇宙気体力学
ttp://lyman.c.u-tokyo.ac.jp/~hachisu/lecture/hydro/num2012sum.pdf
↑圧縮性流体勉強用
- 21 :
- >>16
>>18
ありがとうございます。なんとなく筋が見えてきたのでやってみます
- 22 :
- >>20
流体モデルでは真空中に広がっていくような計算は出来ないってことですかね?
そういう問題を解きたいのです.
- 23 :
- .. z==ニヘ ̄ ' 、
.: ´二>=、 /V//,、
/ / 孑z .ハ ヽY ;
/ / 斗γ⌒ヽ V
、 i,..イ ,:'入_l{ リハ
!ー | .〃{゚::!ワ __.z:r! !
`ー=イiト. ::i`  ̄ (゚'ハ/i! i !
リヘ i !:、_ _,r=1リ 人 )
.::ハト.r!>― <ーイzイゝ-′
/イ .:リ >.::∀:.L==イ.リ
/.ノ .:/! 込/_Y_| !1
..ノ/イr-ミ___..。!!__L__ /.:.、
≧=z' ( ( !:::::::..`.<::::::..,.´.:::::::イ_ヽ、_
 ̄ ̄≧=-=:::::..>゚.:. <⌒¨´
了三三ニ)
 ̄ ̄
┌─────────────────────┐
│>>22 ふふん、流体衝撃波とは別物になる悪寒 .|
└─────────────────────┘
無衝突衝撃波のミクロおよびマクロ構造と安定性
ttp://www.esst.kyushu-u.ac.jp/~space/research/shock_j.html
- 24 :
- >>17
>今回の経験が生きてくると思うので頑張ってね!
はい、ありがとうございます!
時間の刻み幅のことはぜんぜん気づいてなかったので、
さっそく元のdt=10msのほかに、1msと0.1msも計算してみました
収束後(計算開始50秒後)小数点以下6桁まで表示して比べたところ、
この範囲内では流速の値は同じでした
それで、メッシュサイズの違いによる流速分布の変化をもう少し細かく見てみようと思い
追加で21×21、121×121、241×241を計算し、同じグラフに
プロットしたところ、どうもメッシュが粗いと解が鈍ってくるように思えてきました
お目汚しの速度分布(恥ずいので数日で消すかもです)
ttp://www.alicenetwork.net/picture/view/alp5108
まだよくわからないのですが、いちおう今の時点では
数値粘性で鈍ってると考えてみようと思ってます
「幻影解」はテラオソロシス……ですね
それにはまらないよう、気をつけたいですw
どうもありがとうございました
- 25 :
- /l, -‐. ´: : : : : : : :  ̄: : ..`. . 、
_i : :: : : : : : : : : : : : : : : : : : : : : : ヽ
∠ : : : : : : : : : : : : : ; : : : ハ: : : :.: : : : : .ヽ
/: : : : : : : : : : : :':/l : : / 丶 : : ト; : : : : :丶
/: :/: :.' : : : : /: .:/ !: ./ `.:、: l !: : : : : : `、
i: :, : : ,' : : : -A-/-‐ヒ/´ `乂_l:_!_: : 、: i
l: : : : : : : ;ィ.:/ ,レ〒=、 ,=〒=、ハ: ,.: .:. :!
│: : : : !: :i レ,イ,.::::`-'! l.,:::`‐'lヽ V : : : l
l: .:/ニi: :l i |ハ.::::.ハ} {ハ:::::.ハ! }/: : .: : l
|: :{ (l: :l ` ヽ.二.ノ ヽ.二.ノ ´i: : .:/V
|: :ヽ、 !: l `` ``│iV
|: : : : ‐!: l、 __ __ l: l
|: : : : : :V: :\ {  ̄ ヽ::::} ノ:/
|: : : : : : : : :: : `> ゝ、____)/ .-‐.´:V:l
|: : : : : : : : : ;イ l ` フ=く_}ヽ: : : : : : :.!
|: : : : :. :: / | ヽ/!:::::::l l lヽ : : : : l
| : : : : / │ ヽ 〉:::〈/ ! ヽ: : : l
┌─────────────────────┐
│>>24 >プロットしたところ、どうもメッシュが .|
│ >粗いと解が鈍ってくるように思えてきました .|
│ 逆に考えるんだ .|
│ 「格子幅が小さくなると解が微分方程式の .|
│ 厳密解に近づく-数値解が収束する」 .|
│ と考えるんだ .|
│ .|
│お手本になる数値実験お疲れ様でしたけど .|
│格子幅を変えて解の収束性を確認するのは .|
│やっておくべき計算結果の検証方法の一つです .|
└─────────────────────┘
流れの数値解析における誤差評価と単調収束性
ttp://www.nagare.or.jp/download/noauth.html?d=21-1-t01.pdf&dir=17
単調に収束しない場合もあるかもしれない
- 26 :
- >>25
ありがとうございます!
>逆に考えるんだ
涼解です!!
>やっておくべき計算結果の検証方法の一つです
ぐはっ
なんか基本のきの字が抜けてたみたい
「あたしって、ほんとバカ」orz
>単調に収束しない場合も
むーん
「これが解だ」と思ってたら実は違ってた、ってことが
現実にありうる訳なんですね
ガクブルすぎwww 気をつけたいです
どうもありがとうございました
- 27 :
- 下記スレ量産、2ちゃんねる経済板荒らし、自称坂本陽@豊中市
「糞汁は経済板から消えるべき」
「産業歴史博物館設立実行委員会(エンジン博物館)」
※スレタイに「北斗の軍」と入ってるスレ立てもこのクソ野郎の仕業
(´-ω-`)y-~(禁煙マスター) ← 史上最悪の中年荒らし (トリップを変えて別人になりすます技も覚えて来たようですw)
こいつのズリネタ http://ameblo.jp/garo-farann/image-11515076726-12507503444.html
得意技:情報漏洩、ナマポ受給、糞ブログコピペ、ケムトレイル、人工地震、かまってちゃん、驚異の前言撤回w、都合が悪いとすぐ話題脱線
ナマポで助けてもらい http://toro.2ch.sc/test/read.cgi/occult/1354090242/78
日本人滅びろ発言 http://kohada.2ch.sc/test/read.cgi/eco/1342954026/993
情報漏洩をきっかけに http://kohada.2ch.sc/test/read.cgi/eco/1326591504/646
http://kohada.2ch.sc/test/read.cgi/eco/1342954026/926
別れの挨拶1! http://kohada.2ch.sc/test/read.cgi/eco/1356935230/329
別れの挨拶2! http://toro.2ch.sc/test/read.cgi/occult/1354090242/303
別れの挨拶3! http://kohada.2ch.sc/test/read.cgi/eco/1356855992/54
そしてKの逆恨み! http://kohada.2ch.sc/test/read.cgi/eco/1342954026/996
http://ameblo.jp/garo-farann/entry-11462157964.html
精神異常をきたした http://kohada.2ch.sc/test/read.cgi/eco/1374582405/318
超低学力! http://kohada.2ch.sc/test/read.cgi/eco/1369571885/587
さらにロリコンも暴露! http://kohada.2ch.sc/test/read.cgi/eco/1376273900/580-629
何を言ってるのか? http://kohada.2ch.sc/test/read.cgi/eco/1375448326/597
無駄にスレを作ってる張本人w http://toro.2ch.sc/test/read.cgi/occult/1349421194/256
これはひどい http://kohada.2ch.sc/test/read.cgi/eco/1372265052/512
- 28 :
- 落ちてたのか、気付かなかった。
- 29 :
- 最近擬似圧縮性法って知ったのですが
これがあまりメジャーじゃないのはなぜなのですか?
- 30 :
- ' /: /: : : /: : ィ≠^!ミ: :: : : : : \__ノ〉
,' /: ://ハ: : : //´ |! ヽ: : : : : : : イ
|.::': ://: : : / }ツ l! \: : : : :トヽ
l:::!://: : : ,' /__ |! i: : : ::| ヾ!
. ',:| /.: : : :|::| ___ヽ / - ―ハ: : : :| !l
、_Y/: : : : !::衍 :i:iテx / xt≠チ、}: : :! !|
 ̄ 7: : : ハ个j弋少 li::ij:ツノ/: : :/ リ
|ハ: : |: :ヾ:l ` ̄ ./: : /、__ノ
|!:レヽ!: : ヽ> 、 '__ /: //! ̄
|! ヾ:ヾ个 、 - イ: /ハノ
し }: ::|ヾ人` ‐z‐ '7ハ彡ソ| リ
|: 」|: :|辷_7∧/ |: : ::! ̄ ̄ ̄ヽ
|: :|!: :| / ∨__入_.!: : :ヽ、 )、
./!: :j: ::ヽ/、/レ' ̄三 7ヾ: :ヽ__ \
./ |::/ 八:_::>/ハ \ 〈 T  ̄ヽ.. \
广ヽ i/::/ L _ハ!ヽj ≫'!: ! \ \! \
| ∨::/\ ハ ん_ノ ./ | | ヽ ',. \
┌────────────────────────────┐
│>>28 まあ過疎スレだし多少はね? .|
│>>29 利点としては、圧力のポアソン方程式を収束させて解く必要が.|
│無くなる事だろ? stiffな問題とか圧力が振動するような問題には .|
│有効じゃないか? 要は使い所の問題で使ってる所は使ってるって .|
│感じじゃないかなぁ? .|
│スペクトル法を使ったDNSはポアソン方程式一発で解けるので .|
│擬似圧縮性法は必要無いな .|
└────────────────────────────┘
せん断流の数値解析について
ttp://airex.tksc.jaxa.jp/dr/prc/japan/contents/AA0004168026/nalsp0019026.pdf
擬似圧縮性解法を用いた二次元円柱渦励振動の数値計算
ttp://www2.nagare.or.jp/jscfd/cfds14/pdf/c02-2.pdf
- 31 :
- , ,.ィ'´ , `ヽ、 :: :: :: :: :: :: :: :: :: ::
ィ'" , / , /.ハ、.ム た 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 .|
│ttp://www.jsces.org/research/hqc/gallery_file/13_jsme_sangaku_120910.pdf .|
└────────────────────────────────┘
コンピュータによる流体力学
ttp://www.amazon.co.jp/コンピュータによる流体力学-J-H-ファーツィガー/dp/4621061658/
この本の2章とか具体的に何について気をつければいいか書いてあるけど
- 32 :
- - ── -
-‐ァ'´ ` 、
. / / ヽ
.′/ /! ∧ ト、 、 ヘ
| i / / j / ! j ! } ',
| l! / / /イ |/ j八 } }
| 八 /}/ N N }八乂
| Y^Y| |:| |:| ハノ
| {( || U U {ノ AA優先で文章を端折った
| 乂||x:x x:x〉 今も反省してないけど
| | ||ト.丶. _ ノ ,.ィ7
| jl<{ {` マニニフ ´}ト、
| / > \/{ }V < \
┌──────────────────────────────┐
│V&Vとコンピュータによる流体力学両方の視点で書くと .|
│1.数学モデルの選定 .|
│非圧縮/圧縮、非粘性/粘性、2次元/3次元など解く方程式を選定します .|
│2.計算モデルの選定と検証 .|
│離散化手法、計算格子、解法の選定をします .|
│理論や解析解との比較やベンチマーク問題を解いてバグ出し、数値的 .|
│安定性確認、格子幅による離散化誤差の見積もりと解の収束性の確認.|
│3.計算結果と実験結果の定量的な比較と検証 .|
│2で数学的な問題を粗方潰しておいて、物理的な問題をここで検証します.|
│ .|
│って感じかなぁ?詳しくは勉強会やゼミで議論するべきかもしれない .|
└──────────────────────────────┘
- 33 :
- 法的にいいのかよ、これらってさ。
http://manta.blog.jp/
http://sisutore.blog.jp/archives/1000786211.html
- 34 :
- / レ': : : : ―.=7yー 、 ヽ
./ ./: フ´: : : : : /: ̄\ \ ∧
' /: /: : : /: : ィ≠^!ミ: :: : : : : \__ノ〉
,' /: ://ハ: : : //´ |! ヽ: : : : : : : イ
|.::': ://: : : / }ツ l! \: : : : :トヽ
l:::!://: : : ,' / |! i: : : ::| ヾ!
. ',:| /.: : : :|::| 二二ゞ′/ ヽ- ―ハ: : : :| !l
、_Y/: : : : !::衍 :i:iテx / xt≠チ、}: : :! !|
 ̄ 7: : : ハ个j弋少 li::ij:ツノ/: : :/ リ
|ハ: : |: :ヾ:l ` ̄ ./: : /、__ノ
|!:レヽ!: : ヽ> _' /: //! ̄
|! ヾ:ヾ个 、 ´=' イ: /ハノ
し }: ::|ヾ人` ‐z‐ '7ハ彡ソ| リ
|: 」|: :|辷_7∧/ |: : ::! ̄ ̄ ̄ヽ
|: :|!: :| / ∨__入_.!: : :ヽ、 )、
./!: :j: ::ヽ/、/レ' ̄三 7ヾ: :ヽ__ \
./ |::/ 八:_::>/ハ \ 〈 T  ̄ヽ.. \
广ヽ i/::/ L _ハ!ヽj ≫'!: ! \ \! \
| ∨::/\ ハ ん_ノ ./ | | ヽ ',. \
| | j' |\\_____/ /! ,' ', ..i: \
| | | :!. \\ __/ ! / |: : :l ヽ
┌───────────────────────────────┐
│ふふん、4月1日がもう少しで終わるので未確定情報でも書こうじゃないかw.|
│ある計算で二次元キャビティ流れで圧力の境界条件をノイマン境界条件 .|
│で解いたら128x128でもRe=1000でGhiaの論文の値からかなりずれたじぇ .|
│ .|
│Web上だとMac法、圧力をノイマン境界条件で解いてReが大きくなると .|
│流れ関数-渦度法と違ってくるとか報告してる資料があるが本当かもねw .|
└───────────────────────────────┘
- 35 :
- >>31-32
dddクスです!
V&Vの資料の考え方は、シミュレーションの結果を他の人に伝える時にも
超重要っぽいですね
ただ、お話はけっこうむずそうorz 具他例とかもあったらわかりやすいかも……
とりあえず素人な自分は、変な結果がでてきたらそれに気づけるよう、
ファーツィガーさんの本をぽちってみました
着いたら、気合入れて読んでみます!
>両方の視点で書くと
むーん……
この辺をよく理解してないと、恐ろしいことになっちゃうんだ
2あたりは見当も付かないですorz
- 36 :
- 多額な借金で炎上中。
http://manta.blog.jp/
http://sisutore.blog.jp/archives/1000786211.html
http://livedeaichat.blog.fc2.com/?no=1
- 37 :
- __
,. ´ `ヽ.
/ ̄`ヽ / , ‐ 、 , ‐'.
/, ⌒ヽ i | `⌒ ⌒i
{ { 、__ノ ,ノ |^v=「  ̄`H´ ̄}
い、 __ノ}:ーァ! ヽ__ハ__,イ‐.、
,ゝ、 ,イ.:|{`┐ ,. -‐〜ー〈:.:. :.\
/. .:.:.`i_}.: :ヽ`ーV : : :, へ.:.::',:.:. :. :.ヽ,
/ . : .:.:./´ _}ヽ : `!7 : /.:.::::,ハ.::}.:.:',:.l.:. :.',
/ :, -‐ァ´| _」丿: : !{/iヽ.:.::/: :l/.:.:.:l/.:. : :′
/ :/.:.:.:/ / 、 ―j.:. .: .:i―┼‐∨.:. .:. .:.:.:l:. : : : :i
|: : .:.:.:.i | ,l ̄l:.:. .: .:|__|__|:. :. :. .:.:.:|:. : : : :|
|: : .:.:.:.l l/ | |.:. .: .:| | |:. :. :. .:.:.:|:. : : .:/
\: .:.:.:\ \,.| l.:. .: .:!__!__!:. :. :. .:.:.:!: : /
┌──────────────────────┐
│>>35 >具他例とかもあったらわかりやすいかも…….|
│ >2あたりは見当も付かないですorz .|
│ほほっ、実はお前さんは既に実践しとるんじゃよ .|
│ベンチマーク問題=二次元キャビティ+Ghia論文 .|
│格子幅による解の収束性の確認=>>24の結果 .|
│ .|
│今回の題材は>>32の1、2の練習にもってこい .|
│ .|
│孫悟空とクリリンがアルバイトをやってるうちに強く .|
│なっていた亀仙流修行術みたいなもんじゃな .|
└──────────────────────┘
- 38 :
- ひんがら目気色悪すぎこっち見んなR。ひんがら目気色悪すぎこっち見んなR。ひんがら目気色悪すぎこっち見んなR。
ひんがら目気色悪すぎこっち見んなR。ひんがら目気色悪すぎこっち見んなR。ひんがら目気色悪すぎこっち見んなR。
ひんがら目気色悪すぎこっち見んなR。ひんがら目気色悪すぎこっち見んなR。ひんがら目気色悪すぎこっち見んなR。
ひんがら目気色悪すぎこっち見んなR。ひんがら目気色悪すぎこっち見んなR。ひんがら目気色悪すぎこっち見んなR。
ひんがら目気色悪すぎこっち見んなR。ひんがら目気色悪すぎこっち見んなR。ひんがら目気色悪すぎこっち見んなR。
ひんがら目気色悪すぎこっち見んなR。ひんがら目気色悪すぎこっち見んなR。ひんがら目気色悪すぎこっち見んなR。
- 39 :
- 、___ノヽ._,、_ノ、__, __ハ__
) ゝ, ヽ /
ヽ / ,,....:'" ̄:\ ノ も >>27ヽ._
/ | -|-iヽ 〉.:.:.:.:.:.:.:.:.:.:.:.:.\ ) .う >>33 /
〉 ゝ‐' /_ノ (:;.-、:.:.:.:.:.:.:.:.:.:.:.:.:\ヽ っ >>36 ヽ
、 ーァ / 〉´ノ.:.:.:.:.;:ィ'"`ヽ、.:.:ソ ! ! >>38 (
ノ (__ /^Lノ 、:.:.;.ィ''、 /| ハ .\:.) /
ヽ ー、 /\・・ (', /l Xヽl. /,__Xli ヽ}-、,.--v'⌒
〈 ‐'゙ / \ (./ / >l/ <` lノ.lヽ.:l
〉 | | | / 〉〇‐'____ .ヽOノイ/
〈 ・.・.・ /^Lノ {ヽ_〈 `ヾ! ,,.ィ':/
. ) r゙'.:.:.`''''===ィ'.:.:/
},..-、,.---、,..--、f゙'.:.:.:.:.:fi.:.:.:.<、了
/.:.:.:.:.:.:.:.:.:.:.:fjニ).:.:.:.:ヽ'ソ
./.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.、
┌──────────────────────┐
│>>35 >シミュレーションの結果を他の人に伝える時 .|
│ >にも超重要っぽいですね .|
│そうですね .|
│企業の技術者が会議でフルボッコされないための .|
│方法論的な側面も強い気がするかもしれないw .|
│ .|
│>ファーツィガーさんの本をぽちってみました .|
│主に取り扱っているのが差分法と有限体積法です .|
│3、4、6章で扱っている問題のソースコードがftpで .|
│落とせるので無駄な買い物にはならないと思います .|
└──────────────────────┘
- 40 :
- ド素人が、実機流体の温度分布や速度分布、あるいは配管の熱応力分布を
可視化して表現するためにはどのようなプロセスを踏んだら宜しいか?
- 41 :
- -─‐ 厶 }}、
_ -- ´ `ヽヘ
{{`7 | \ ∧
∨ '/ /│∧ { ⌒\ リ r-‐-‐、
,' /⌒|/ \x‐=ミ∨〈/ / /,ハ_
| i i/x=ミ. rしハ }!| ゝつく ノ∧ 温度境界層だにゃん
_| i | i{ {しヘ V_ソ | |ハ  ̄ '、
_r'「{∨|i∧ Vソ '' | |ノハ___ノヘ
{_{└と 从 ハ '' , | |'. ∧
| `┴'' | |i八 ー イ|/ー¬ヘ -ヘ
l. | |\|` ーzr‐≦/ │ |\
ト-― i_L. イ⌒´/ ∨:::∨ 」 }
人 〉 〉〈 > __ノ
{ \ ∨ \ ∨ / ア「
\,,____、| } /|│ /⌒)
┌─────────────────────────┐
│>>40 解きたい問題が、自然対流なのか、強制対流なのか?.|
│流体の圧縮性は無視できるか?(数学モデルの選定) .|
│先行する研究があるかどうか?(計算モデルとベンチマーク) .|
│あたりから調べましょう .|
│ .|
│強制対流熱伝達の場合は速度場と温度場が独立なので .|
│まずは温度場を無視して速度場がちゃんと解けるか調べ .|
│ればいいと思うけど .|
└─────────────────────────┘
平行平板間流れの層流熱伝達(理論解析と数値計算)
ttp://fluid.web.nitech.ac.jp/Gotoh_Home_page/Edu/Public_course/Text/Text_6th_pm.pdf
プログラム付き強制対流熱伝達の資料です
これを読んで全然分からないなら勉強が必要かもしれない
- 42 :
- /⌒)__/⌒\: : : ``丶、__´'くi:i:i:_〉
.′/⌒{\: : :.\ : : : : : \``丶、
{ ./.: : :从: :丶: : ∨/: : : : : : : : : :\
{/.: : :/.{ \: : \∨/ : : : : : : : : : : \
___/.: : :,′{ \ : :`∨\: : : }}: : : : : : : \
⌒,′: .:′、 \:__|=- ≧=--‐: : : : :|:
′: : i{___ \ ´ i|≧=-: :.}}: : : : : : : |:
___彡: : : :i{≠ぅミ k抖竺=ミ}}: : : : : : : |:
}: : .:}八{ {/{ハ. ⌒ {/{_い刈: : : : : : : |:
,: : : :乂__ Vrク 乂rク ノ :/.: : : : ; :
′: : : }}' ' ' ' ' 'ア: : : : : :_/_:
i{: : :‘, 圦 r--‐ vi /}.: : : :/ }
八: : ∧: : : . ノ __〉: : / ん
/ ):/ : }}:. : : : : . __〉: : :{ {/
_/: : :.}} : : : : : :}>=ぉi爪 __ア. : :/} '{=
⌒ア: : :}}: : : :j{ ̄ ̄ ,ハ ____ア゙.: :/ ,}i '{_
/: __.: :>‐j{ __/ }_ア゙.: :/__-=圦 /⌒
//___\ j{ .「i:i:ヽア゚.: :/:i:i:i:i:i} ___} _/
/  ̄ ̄} _,ノi:i:i:/: : : /:i:i:i:i:i:i:i:i|__ア_/ /
┌─────────────────────────┐
│ふふん、どんどん逝ってみよう〜 .|
│流れが層流の場合、管内層流熱伝達のヌセルト数が解析 .|
│的に求められるので計算モデルの検証に使えそうだな .|
│ .|
│乱流熱伝達は結構難しいと思う データをどうまとめれば .|
│いいかも良く考えてシミュレーションを行おうじゃないか .|
└─────────────────────────┘
熱交換器の伝熱評価
ttp://www.yama.tus.ac.jp/kikai/tempfile/6-heat.pdf
局所加熱された円管内の強制対流乱流熱伝達 管壁内熱伝導を考慮した数値解析
ttps://www.jstage.jst.go.jp/article/kikaib1979/53/492/53_492_2528/_pdf
- 43 :
- . ´  ̄ ` . | j i/〉、
/:::::::: /!::ト、::::::::::::::\ |∩ /
/::::: ;ィナ/ |::lxト、‐ 、:::::::::: ヽ |j ) ノ
イi:l::ル笊 人l ノ笊ハi:::::::::::::: ヽ /`V イ
从i:| i.r' li l.r' lト |::::::::::::::::. ヽ /: \ Xノ
/:::l:| |匕リ l匕リ .|:::li::l::ヽ::::ハ/: : : : : :\/
rf「ト、 /:::::l〈"" _ "" |:::li::|ルハ//: : : : : : : :/
| i ___r‐ルル> .「 ノ .イレハj: リ/: : : : : : : :/
〈.ヽ/´ ̄/: : : : : : :| l| |l |〈 {_て .フ´ /::/: : : : : : : : /
ーl l: : : : : : : :| l| |l |トト、 xへ /: : : : : : :: : :/
人_人: : : : : : :! l|__|l lハ: ≧≠≦=虫: : : : : :イ
 ̄ ̄ ̄ ̄〈∧___|l |7∧: : : : : : : : : : : : : : : : ト、
く\∧fi|l // ∧: : : : : : : : : : : : : : : |: \ _
┌───────────────────────┐
│>>40 どの程度の経験があるか分からなかったので .|
│簡単な流体シミュレーション経験があるとして書いて .|
│みました .|
│ .|
│シミュレーション経験無しなら二次元キャビティ計算 .|
│あたりから始めてみては? Festina lente .|
└───────────────────────┘
- 44 :
- http://2ch.sc/
otonaをつけたらいけるよー
↓
http://otona2ch.sc/
- 45 :
- ,r-―――-- 、
./ -=ヽ.o,r-\、
/ ,/⌒`ー--<´ヽ}
.{ ./ ,/V,、 .ト, i、__ ヽ
.{ /i { ● ` リ●`i .}.__
ヽ{∧ リ (^⌒) ,} ,}{_ノ
{ ヽ i`ー-ニ--</ /´
.} /`ー-'Y:::ヽ,__Yシ
/ {ー-/ _,`ー' ) {
/ .,r┴=< .: ヽ .{
.{ ゞ/´`ーt'^'i´ヽヽ,
.{ , >t-<,_,,,,,,_.〉 }
`ヽト./:::/iNヽ/ ヽ:::V,リ
. /、/`^ ′ ir-ヽ
ー’ `ー ′
┌───────────────────────┐
│着た北来た喜多、乗ってきた〜 .|
│XP難民の救世主Ubuntu 14.04 LTSついにリリース! .|
│なんと5年サポートだぞっ .|
│ttp://blog.livedoor.jp/itsoku/archives/38345124.html .|
│という事で人柱逝きます! Unityだけは簡便なw .|
└───────────────────────┘
- 46 :
- ====-、 ゙ミt、 ,、}ミi lミi
====-、\ Vl}.//.l/ {/
\ .}ニ/
lミ/ -==-、
ヾtzzzチ' ノツ
___
,.'" ̄ヽ ,..-、 / ,_,ィ ,..ィ ,. `ヽ、
,' l ./ V /ィ_-Z_l/ {/ヽ、! . ヽ.
,' l._ ,' .l /ゝ' ̄ ヽ‐-Vl ゙f'i
! .ノーZ! l ,ハ / ̄ ̄`ヽ`ひ 〉,、 〈
、 ., .l 'ヾt、 ,..-‐' ゝ\ ノ .// ノ /}
`7"::`ヽ、 ヽ.__ノ l `' .、/ィ_,,...-‐ァ''アミr‐_,,..ィr''"~" /:::::
,'::::::::::::::::`::ー'_ `ヽ l .._ `< ン/zチ=リ.‐''"/\ / ,':::::::
.!::::::::::::::::::::::::::::ヽ.__. V \ .Vュヽニ二ノ`'ー'',..- l l:::::
l::::::::::::::::::::::::::::::::;l :l! ヽ ムf"¨>‐、 ノ .l : l::::
ヽ::::::::::::::::::::::::::::;'、_./ ゙; ムー/ / \_ / ./:::
ヽ::::::::::::::::::;;::-‐''゙ ',. ム! i ` , ', ./::
 ̄ ̄ ̄ ̄ヾミヽ ヽ `ヽ、 .,.ィ'ノl/::
::ヽ.._ 、 /:`゙''':::::
:::::::: ̄`::.....___/::::::
┌───────────────────────┐
│1.初っ端、GRUBをMBRにインストールするのを忘れて .|
│起動せず=>1回目ちんだ〜 .|
│2.Unity から GNOME compiz へ移行、仮想画面使うと .|
│制御不能=>2回目ちんだ〜 .|
│3./home だけフォーマットしなかったからか?起動時に .|
│毎回、「エラー: malformed file」 どゆこと〜 .|
│ttps://forums.ubuntulinux.jp/viewtopic.php?pid=103111 .|
└───────────────────────┘
- 47 :
- -‐=‐- ミ
. 7ア⌒`__'⌒` 、
. ,.イ.: :7:. :.`:: ..、 V)__
゙`ー=ァ´.: /.: :/.: :. :.、:. :. ゙'く⌒`
/ /.:: .:/.: :/.: ハ:. :. :.\:. :\__)
. {_..イ.:/.: .:/.: :/.: / jト、:. :. :. \:. :.〈
. /.: :l.: .: .: .: .: {/- 、 Χハ:. :. :ヽ:. :.
/.: .: |.:{.: .:.{.:: 仏=ミ ィf〔7:. :. :. }:. :!
. {.: .: :八:. :. :. :. 「V_ソ 弋ソ ハ:.. :.八_{
八: /.: .:Vハ:. :.|`` ' `八! :/
/.: :イ.: .: :: :: :V 个ト . マフ ,.く
. //"j|.: .: :: :: :/〉八 〈〕T爪 :.
/゙ |.: .: :: ::/`ヽ〉:.\__八__) :i
┌──────────────────────────┐
│ふふん、取り敢えずUbuntu14、エラーが出なくなったので .|
│早速>>4の幾つかを入れて逝こうじゃないか .|
│ .|
│まずは、channelflow のインストール .|
│最新1.4.2をconfigure〜makeすると案の定失敗したorz .|
│どうも、octave の oct.h がかなり内容が変わったみたいだな .|
│己はウィザード級のハッカーじゃないからw ソース改変無理 .|
│ここは、octaveツールを捨ててインストールしかないじぇ .|
│ .|
│対応は簡単で configure の後にできる Makefile を開いて .|
│am__append_1 = programs-octave をコメントアウト .|
│これで、programs-octave フォルダに入らないのでmakeが .|
│無事終了して、インストールできる .|
└──────────────────────────┘
ちなみに、octaveツールは乱流の不安定周期解?を得られるというもの
ttp://www.math.wisc.edu/~waleffe/ECS/ECS2.pdf
↑Waleffeの解と同様のものが求められるみたい
- 48 :
- /  ̄`ー、 __,,.--‐‐――‐-、
, '´ : ヽ._,,.-‐¬ ̄ : . \
/: . . :.::{::.:. : : :{:.:.. . : : : . . : : : . . /^
.′:. . _ 、 :、.:.:::i;;:.:. : :|::::::.:.__;;.-―‐- 、::_:.:. : /
|;;:;::.:. . . | ヽ{  ̄ヽ;_}:_,.-‐ァ'´::.::.::.::./ /::.ヽ`ー:'ヘ
l;;;;:;:;:::.:. : . . ! j!,.-‐ Vハゝ'´::.__::.::.::.::.::.^^::.::.::.::\::.::.:{
二i;;;:;:;:::.:.:.:. : :// ‐' !::ハ::.::.{ __ヽ.::.::.::.::.::.::.::.::.::.`、::.:ヽ
._-三_,};;;:;:;:::.:.:.:.:/レ' ! }ノ ヽ::ヽ二ノ::.::.::.::.::.::.::.::.::.::.:〉::.::.:
-;/.:.::t;;:;:;::::.:.:.:j -‐' _,,. rく´ \::.::.::.::.::.::.::.::.::.::.::.:/::.::.::.:
/ .:.:::ヽい从;{ __,,,.. -、___,{イ `ヽ、::._::.::.::._;;.-'´ノ::./
!::.:....  ̄:./ ヽ_ヘヽ、 ___ __,二 -‐=-‐'´
ヽ:_:::::.:.:.:.... /___ 、 v'´ .:.::_>'´ ̄~ ̄ ̄
 ̄`ー‐''´三`ー-r一''¨´ ̄
 ̄―二三三f´,r'´ . . :.:._;.-‐‐―‐‐―‐--- ::.__:.:.:._z、__;;..
 ̄三^ヽ:.__;.-‐' 三二― ̄  ̄ 二―  ̄
┌──────────────────────────┐
│>>4 の2番目の surfer もコンパイルエラーが出たりしたけど .|
│動くところまで確認できたよ .|
│手順を書かなきゃRないんだけど... .|
│ .|
│パトラッシュ…疲れたろう 僕も疲れたんだ… .|
│なんだかとっても眠いんだ… .|
└──────────────────────────┘
- 49 :
- / `丶、
」∨ ヽ
7 '.
/ ノ| ∧ `ト N| っ
. /N イ/ |/ ∨| ∧ | っ
. 〃 Y⌒i W∩ ∩ }/V -=ニ二三三三三
/ {{{ | | ∪ ∪ {{-=ニ三/ ̄ ̄`'く>三三三
/ゝ _| | '' '' 人.三三|(‐〜〜‐)| |三三三
/ 厶-=ニ | :八‐-(_:::::::⌒)<∧ノ三三| ー--‐' |_|三三三
-=ニ二三三三三\_,,,フ三三三三三三三三| <・ ) 「 三三三
二三三三三三三三三三三三三三三三三三ニ圦 {__>,,ノ三三三
三三三三三三三三三三三三三三三三三三三ニ≧===≦三三三
┌──────────────────────────┐
│surfer インストール方法は初代スレ 164〜330 参照ですけど .|
│ttp://ikura.2ch.sc/test/read.cgi/sim/1303299935 .|
│Ubuntu14 は mpich のインクルードファイルのパスが今までと .|
│変わったので GetpartC のコンパイルは簡単になりました .|
│ .|
│surfer のインストールは問題無いかと思ったら .|
│./Drop < falling.par でエラーになっちゃった .|
│surfer-3.1.7/source/Drop.f の32行目に rewind(11) .|
│surfer-3.1.7/source/time_step.f の85行目に rewind(66) を .|
│追加して make し直せばエラーが出なくなります .|
│ .|
│ttp://www.alicenetwork.net/picture/view/alp5143 .|
│./Drop < falling.par してgnuplot.3000を描画〜良いかもしれない.|
│という事で土筆のはかまを取る作業に戻ります .|
└──────────────────────────┘
- 50 :
- 非等方圧力で圧縮性流体(roeやhll)といてる先行研究とか知ってます?
- 51 :
- ,. -─-、
/ ̄`ヽ { ヽ
{ ゝ-──-ゝ ノ
`7 ハヽ ハヽ <
/ しノ しノ ヽ
│ ,.-‐(^>─-、 |
│ { _人__ \ ノ
'、 くマ{ V7 ̄ノ /
\ \ゝー ´_/ /、
/ ヘ._ ̄ ̄___, イ| \
. / ヽ__\ ̄ /__」/ / ̄`ヽ
/\ <.__\/_」 { { }
. ( \,.ィ  ̄ ̄ \ゝ.___,ノ
`ー一'′| o `|
| o |
┌────────────────────────────────┐
│>>50 前スレ259の人かな? .|
│A New High Resolution Scheme for Compressible Viscous Flow with Shocks .|
│ttp://aero-comlab.stanford.edu/Papers/AIAA-1995-466-282.pdf .|
│圧縮性流れのアルゴリズム(II) 離散化方程式の解法 .|
│ttp://www.ifs.tohoku.ac.jp/edge/publications/implicitmethod.pdf .|
│↑こっちは先行研究とかじゃないけど国内の研究 .|
└────────────────────────────────┘
- 52 :
- 寧ろ逆に天国だろ
他に別に極楽だろ
当然正反対に理想郷だろ
- 53 :
- >>51
ありがとうございます
でもこれらは非等方圧力ではないですよね?
- 54 :
- __ ,' な
,.. '" ,:  ̄ ヽ、、! に
,.' , /. / l |、 、 ゙'.|. 言
l ,' iーエ'、 ヽ;゙ー- l、:. | っ
l !//l:::::i l::::ヽヽ`i :、! て
lヘj 、::ノ 、::::リ | lナ| ん
〈 " ¨ " | レ゙| の
l>...._ ,.‐-、 _,,.ィi ⊂. ?
| :,イヽ`ニア:/iヽ' |
|./ ゝzヽ会't'' ヽ:i ヽ._/
/ ヽ :|
┌─────────────────────┐
│>>53 非等方性流れ≒非等方圧力 じゃないの? .|
│ .|
│非等方圧力をどういう意味で使ってるのかと .|
│どういう流れ場を想定しているか教えて欲しいけど.|
└─────────────────────┘
- 55 :
- >>54
うちの業界だけの言葉なんですかね?
圧力テンソルの対角成分が同じものでないことです。
例えば二次元なら
∇・P = (∂x Px) ex + (∂y Py) ey
- 56 :
- . - ── - .
/ ` .ノ}
. . ' ト、
/ / ∧. ∧ Vハ
' / . 斗ト / ∨ ,ィ十 、 Y
. i / /レ__レ' ∨ _ハ! ハ l
. | |i /〃云k .ィ云`ヾリ | ?
从 iル' ! |l匕l| |l匕l| }l| ト、 |
. 人ハ! .t仆ソ t仆ソ ノ | | i. !
|N|"""" "" | レ' !
| 〈 「 ̄ ̄ } | | |
| i〕ト . ー‐ ' . イl| | |
| li i i _∧ ̄ ノ^ト、 l! ル' !
| ll l /〈 X´ リ ト、 l
| ll/|_ ∨}__ ∨ .| \ |
/ / V ∨ ─┘ ヽ リ
/ i | \. ∨/ > i | ∨
┌──────────────────┐
│>>55 うん、だから非等方性流れは .|
│∂P/∂x≠∂P/∂y .|
│じゃないの?ってことなんだけど... .|
└──────────────────┘
- 57 :
- >>56
独立な圧力がP一つではなくてPxとPyなんですよ
- 58 :
- 流体中での爆発とか?
- 59 :
- それもあるかもしれないですね。
天体のジェットなんかを考える時はジェットの吹き出る方向と、それに垂直な方向では温度が違うんですよ
- 60 :
- P = P⊥ diag(1, 1, 0) + P|| diag(0, 0, 1)
なんてかいたりしますね
- 61 :
- 通常の流体(等方圧力)では
P = P diag(1, 1, 1)なわけです
- 62 :
- 分子運動論的な言い方をすると、方向ごとに熱速度の大きさが違うので
面の方向ごとに圧力が違うってこと?
そういう世界もあるんだなあ。状態方程式も圧力の複数の成分が入ってるのかな
- 63 :
- 物理現象としては同じものを対象としてないって理解でOK?
- 64 :
- >>63
同じものを対象にしていないとはどういう意味でしょうか。
http://www.physics.wisc.edu/grads/courses/726-f07/files/Section_14_Eq_General.pdf
みたいな感じですね。これは磁場が入ってますけど
- 65 :
- ,、- '''"~~~~゙゙゙ ''‐ 、
,、.'" __,,,,,_ ゙' 、 / ヽヽ / /
,.':: ,、- ''"´ `゙゙゙''‐ 、 ゙' , ,,、、、 /_ \ / / /
,.'::: .,. " ゙'.、r'" ゙'、 / . ̄ ── / ' '
,.':::: , " , ., , ,, 、, , 、 、 ,r" ', / _/ / /
_,i、 ,.' .i ,','l !l .!l l. ! ',, l .l、/ ........ !
゙,~ ,.ri. .{ .!.!,i,',l ,',!,! l', l',l,,,!/ .::::::::::::::!
゙、:::::,.' .::! l゙、'l"l,' !,'!'!'.レ !' .'''.,' .:::::::::::::::::l
,>' .:r"l ',.l,、‐''''i- ''i;' .:::::::::::::::::::,'
〈-、::::゙'‐!. ,. '、ゝ;;;:j /. .::::::::::::::::::::,'
゙、 ゙''‐-l ', .ヽ`.´ ' ,.' .::::::::::::::::::::/ 非等方圧力
`,i‐、..}i '.;:. ゙ヽ ヾソ' .:::::::::::::::::::,.' ちぃ覚えた!
/ ;. ゙,゙, ';: ノ,,、-''" .::::::::::::::::::,r'
非等方圧力MHD平衡解析に基づくLHDプラズマの圧力非等方度同定手法の研究
ttp://ir.soken.ac.jp/index.php?active_action=repository_view_main_item_detail&page_id=29&block_id=155&item_id=4053&item_no=1
Stability of anisotropic-pressure tokamak equilibria to ideal ballooning modes
ttp://research.microsoft.com/en-us/um/people/cmbishop/downloads/bishop%20-%20sapteibm%20-%20nf85.pdf
Magnetic hydrodynamics of two-dimensional flows in a plasma with an anisotropic pressure
ttp://link.springer.com/article/10.1007%2FBF01015087
2.ミラープラズマの安定性
ttp://jasosx.ils.uec.ac.jp/JSPF/JSPF_TEXT/jspf1999/jspf1999_06/jspf1999_06-695.pdf
>>64 電磁流体力学(MHD)のCGL理論では非等方圧力テンソルは磁場に垂直な圧力と
磁場に平行な圧力で与えられる...ちぃ覚えた
- 66 :
- ' /: /: : : /: : ィ≠^!ミ: :: : : : : \__ノ〉
,' /: ://ハ: : : //´ |! ヽ: : : : : : : イ
|.::': ://: : : / }ツ l! \: : : : :トヽ
l:::!://: : : ,' / |! i: : : ::| ヾ!
',:| /.: : : :|::| 二二ゞ′/ ヽ- ―ハ: : : :| !l
、_Y/: : : : !::衍 :i:iテx / xt≠チ、}: : :! !|
 ̄ 7: : : ハ个j弋少 li::ij:ツノ/: : :/ リ
|ハ: : |: :ヾ:l ` ̄ ./: : /、__ノ
|!:レヽ!: : ヽ> _' /: //! ̄
|! ヾ:ヾ个 、 ´=' イ: /ハノ
し }: ::|ヾ人` ‐z‐ '7ハ彡ソ| リ
|: 」|: :|辷_7∧/ |: : ::! ̄ ̄ ̄ヽ
|: :|!: :| / ∨__入_.!: : :ヽ、 )、
./!: :j: ::ヽ/、/レ' ̄三 7ヾ: :ヽ__ \
./ |::/ 八:_::>/ハ \ 〈 T  ̄ヽ.. \
广ヽ i/::/ L _ハ!ヽj ≫'!: ! \ \! \
| ∨::/\ ハ ん_ノ ./ | | ヽ ',. \
| | j' |\\_____/ /! ,' ', ..i: \
| | | :!. \\ __/ ! / |: : :l ヽ
| |::| |:. \\ .// .l!:∨ ∧
┌──────────────────────────┐
│ふふん、ということで>>50の条件に合う論文検索してみたじぇ .|
│"anisotropic pressure" compressible Riemann でググっただけ.|
│だけどw キーワードで絞り込めばRるんじゃないかい? .|
└──────────────────────────┘
MAGNETIC FIELD AMPLIFICATION AND EVOLUTION IN TURBULENT
COLLISIONLESS MHD:AN APPLICATION TO THE ICM
ttp://arxiv.org/pdf/1305.5654.pdf
Numerical Prediction of Three-Dimensional Non-Equilibrium Flows
Using the Regularized Gaussian Moment Closure
ttp://arrow.utias.utoronto.ca/~groth/publications/AIAA-2014-0229-lam.pdf
- 67 :
- 初心者ですが、二流体モデルをやるにはどうすればよいんでしょうか?
- 68 :
- >>67
collisionいれればいいんじゃない?
- 69 :
- -‐=‐- ミ
. 7ア⌒`__'⌒` 、
. ,.イ.: :7:. :.`:: ..、 V)__
゙`ー=ァ´.: /.: :/.: :. :.、:. :. ゙'く⌒`
/ /.:: .:/.: :/.: ハ:. :. :.\:. :\__)
. {_..イ.:/.: .:/.: :/.: / jト、:. :. :. \:. :.〈
. /.: :l.: .: .: .: .: {/- 、 Χハ:. :. :ヽ:. :.
/.: .: |.:{.: .:.{.:: 仏=ミ ィf〔7:. :. :. }:. :!
. {.: .: :八:. :. :. :. 「V_ソ 弋ソ ハ:.. :.八_{
八: /.: .:Vハ:. :.|`` ' `八! :/
/.: :イ.: .: :: :: :V 个ト . マフ ,.く
. //"j|.: .: :: :: :/〉八 〈〕T爪 :.
/゙ |.: .: :: ::/`ヽ〉:.\__八__) :i
|.:: ::. ::.{ V〉:. 「r-<) 、 |
┌────────────────────────────────────┐
│>>67 流体と粒子の混合物の場合の論文 2章あたりを読んでみてね .|
│2流体モデルによる流動層の理論解析と数値実験 .|
│ttp://hdl.handle.net/10636/297 .|
│ttp://repo.biblio.tuat.ac.jp/bitstream/10636/297/1/199703HirayamaOsamu_F.pdf .|
│モデルによって必要なものが変わる上に実験や経験式が必要な印象なのだけど... .|
│使ってる人の意見を求めてみようじゃないか! .|
└────────────────────────────────────┘
- 70 :
- .. z==ニヘ ̄ ' 、
.: ´二>=、 /V//,、
/ / 孑z .ハ ヽY ;
/ / 斗γ⌒ヽ V
、 i,..イ ,:'入_l{ リハ
!ー | .〃{゚::!ワ __.z:r! !
`ー=イiト. ::i`  ̄ (゚'ハ/i! i !
リヘ i !:、_ _,r=1リ 人 )
.::ハト.r!>― <ーイzイゝ-′
/イ .:リ >.::∀:.L==イ.リ
/.ノ .:/! 込/_Y_| !1
..ノ/イr-ミ___..。!!__L__ /.:.、
≧=z' ( ( !:::::::..`.<::::::..,.´.:::::::イ_ヽ、_
 ̄ ̄≧=-=:::::..>゚.:. <⌒¨´
了三三ニ)
 ̄ ̄
┌───────────────────┐
│う〜ん、>>54、>>56はアホな事書いてるなw .|
│流動する流体でも圧力の等方性が成り立つ .|
│圧力勾配を与えると流れが生じるわけだしw .|
└───────────────────┘
- 71 :
- とりあえず、流体がシミュレーション通りに動くことを己の目で見たいのだけど
どこに行けばいい?
- 72 :
- / ` ヽ、
(V ∧ \
/ ハ / ヽ ハ ヽ
⌒ ⌒/// ヽi⌒ヽ ヽ
/ /|/__レ V__ヽ } '
| //rァヽ 'rァ 、V ハ }
| | | { {:r´:} {:r´:} l l ハ
| / ! ! 弋ソ 弋ソ レ イ
| | |"" "" lレ !
| `'| |、 /  ̄ ̄ ヽ ノj |
| V ` ゝ__ ノィ´ V ′
| /〈 \/ヽ\ l
l // ∨、:.:.∧〉 l\ ,′
| /  ̄7 ∨:.:.∨ < ヽl
| / \ ヽ:./ / ヽ
l l l l `ヽ´ ´ ハ
V l ,′ | l
/ V l l
┌─────────────────────┐
│>>71 自分でコードを書いた? .|
│リアルタイムに結果を見たい? .|
│情報が少なくてどうすればいいか分からないけど .|
└─────────────────────┘
- 73 :
- >>68
>>68
ひとまずもう少し勉強してみます。ちなみに液液です
- 74 :
- ', / 単 二
固 気 .l . ・ ・ . . ,' .純 相
液 液 | ..’ ,.-‐ '', ̄, ̄`ヽ、 、 l .に 流
/ / .| , / /,' /| .i `ヾ'| 考 の
固 液 | :./ ,'-/.,'._/ l _l_l__i_ | え 組
気 液 l ,' l./--./ 、l__ヽli::. .| て 合
| / / : l , l/ fミjミi fZミiヽヾ .', も せ
♪ 気 ./ .ヽ、ic ヒzリ ヒzリ l | iY \ ♪ を
気 / ゙, リXx _¨ xX.|.lノ .l//\__/
___,.-、二l ヽ / `i ,.|ノ | :
‘,ヽ>..__ヽ_,,..-r'" l | ,
・ |! ltXjXメ l .| .
: |: .l ', ::|
’ l: | ∨ ::|
┌─────────────────────┐
│なので、重要な情報は早めにね .|
│液液の二流体モデル .|
│Modelling and Control of Two-Phase Phenomena:.|
│Liquid-Liquid Two-Phase Flow Systems .|
│ttp://www.eng.tau.ac.il/~brauner/brauner.pdf .|
│これでもダメなら .|
│"two liquid phases" "Two Fluid Model" とか .|
│"liquid-liquid" "Two Fluid Model" でググるとか? .|
└─────────────────────┘
- 75 :
- fortranのmpi fftwのサンプルコードとか持ってる人いるかな.
fftwはドキュメントが不親切過ぎて...
- 76 :
- (⌒ヽ (⌒ヽ / ∧ ⌒ヽ、
` ` ´ / /! / ! ヽ ヽ ̄  ̄_)  ̄_)
/ / / { / l ヽ ヽ、 ヽ ヽ
/ ,′ / | / ヽ lヽ ヽヽ ヽ ',
′ l / |/ ヽl_ \! \ ヽ i !
l l ハ / ノ >─<´ヽヽl }
| l l T7 ̄_ミ /てヽ ヽ l l_ } \
/ { { ハ l { ん(..ハ !///ハ } | | ヽ ! 、_)
(_/ Vヽ ヘ ! V/し:} {/:し/リ ′ ! l ヽ } │ 、
Vヽ ', V_:ソ ゞ=' l ,′ノ ノ ! 、_)
し ヽ i "" "" | l__/ |
__ i } U /! / |
ヾ:\ ', ゝ u _ -‐ ⌒ヽ ィ レ !
}´⌒ヽ ヽ |`ー- _└‐_ニ´‐/ /\ |
{ ', ∨ /l \/ / | \ |
ゝ、__ノヘ l / レイ: . :ヘ / ,′ \ l
\ \ ヽ / / { 〉: 〈 7 」 ヽ |
┌──────────────────────────┐
│マルチスレッドFFTWなら多少は分かるんだけど .|
│6.13 FFTW MPI Fortran Interface .|
│ttp://www.fftw.org/doc/FFTW-MPI-Fortran-Interface.html .|
│これは確かに厳しいかもしれない .|
│FFTW MPI 使ってるフリーソフトがあれば参考になるんだけど .|
└──────────────────────────┘
- 77 :
- >>76
そしたらみなさんFFTってMPIでどうやってるんですかね?
ライブラリ使わず??
FFTW以外のライブラリを使う??
- 78 :
- そもそも今どきMPIを書くのか
大筋は↓で C から F への読み替えで一回ややこしくなっただけじゃないか、と大雑把に
変換後の値は各プロセスにばらまかれてるから集めるなら MPI API コールかね。
マルチプロセスでスタート、MPI と FFTW の INIT
↓
横に何人いるかは分かってます。
全体でLxMのFFTやりたいんで自分の割り当て量ください。
↓
各自、自分のメモリ取ります。
↓
MPI-FFTプラン立てます。
↓
計算します。
- 79 :
- / ,:イ>i: : : : : ヽ: : :\ ヽ
i /: :l.: :jヘ、: : : : : ヽ : : :\ i
l ,/: : :l.:/ l \: : : : : :ヽ: : : :ヽ、 l
レ: : : :':/ i \: : : : :`!: : : :∨
i:.; : :ル'─ヘ ─ヘ; : : :l: : : : l
l:, : : :行i:i:iテ xi:i:iチ:.: :l: : : : :l
l;.: : : l 弋り 弋乂j: : :l: : : : :l
l: : : : l、"" "" l: : :l: : : : :l
l:.lヾ.: :` 、 ヘ ,: : :l: : : : :l
l:l: :ヽj: : : ;>t-<彡从イ: : :イ
jiゞ,x=≠=、i´fヘツ x‐' i: : :l;.:.l
/r'{二ニ{ili}ー-:、,ゝ、l: : ソ`´`ー、
,x'"7=ミノ iヽ、__,/ ,ノ.~l: :/ i、
, ・´ ,' ゞ{_,ノ-‐'笠ニ-―l:/ ヘ
; i x// ̄´ i:!、 \
i l .i// ≧ーi、__ \
ゝ、_,乂〈i ,//::::::ll:::ヽ、 \
l @``・、 _,・'/:::::::::ツ.l:::::::∧ ゝ
┌─────────────────────────────┐
│>>78 >大筋は↓で C から F への読み替えで... .|
│考えてみたらFFTW MPIもCのサンプルがたくさん出てくるんだから .|
│CからFへ変換すれば良かったんじゃないかw ありがとうだぞ! .|
│6.3 2d MPI example と>>75でどうにかなりそう? >>75 .|
│ttp://www.fftw.org/fftw3_doc/2d-MPI-example.html .|
│マルチスレッドFFTWもCのサンプルから調査したのすっかり忘れてた.|
│ .|
│FFTのテストは三角関数辺りから始めて複雑なものはOctaveや .|
│Scilabなどの組み込み系と比較したな .|
└─────────────────────────────┘
- 80 :
- 前処理付き共役勾配法が今のところ最強ってことでいい?
- 81 :
- / ` ヽ、
(V ∧ \
/ ハ / ヽ ハ ヽ
⌒ ⌒/// ヽi⌒ヽ ヽ
/ /|/__レ V__ヽ } '
| //rァヽ 'rァ 、V ハ }
| | | { {:r´:} {:r´:} l l ハ
| / ! ! 弋ソ 弋ソ レ イ
| | |"" "" lレ !
| `'| |、 /  ̄ ̄ ヽ ノj |
| V ` ゝ__ ノィ´ V ′
| /〈 \/ヽ\ l
l // ∨、:.:.∧〉 l\ ,′
| /  ̄7 ∨:.:.∨ < ヽl
| / \ ヽ:./ / ヽ
l l l l `ヽ´ ´ ハ
┌────────────────────────────┐
│>>80 前処理にも色々あって最適なもの探すのが大変そうだけど .|
└────────────────────────────┘
CG法の最近の前処理のロバスト性と効率化について
ttp://www.kurims.kyoto-u.ac.jp/~kyodo/kokyuroku/contents/pdf/1362-3.pdf
並列非線型共役勾配法アルゴリズムとその性能評価
ttp://www.ssisc.org/swopp09/paper.pdf
数値計算の基礎 共役勾配法の勉強
ttp://www7.ocn.ne.jp/~kawa1/numeric.pdf
- 82 :
- /⌒)__/⌒\: : : ``丶、__´'くi:i:i:_〉
. .′/⌒{\: : :.\ : : : : : \``丶、
{ ./.: : :从: :丶: : ∨/: : : : : : : : : :\
{/.: : :/.{ \: : \∨/ : : : : : : : : : : \
. ___/.: : :,′{ \ : :`∨\: : : }}: : : : : : : \
. ⌒,′: .:′、 \:__|=- ≧=--‐: : : : :|:.\
′: : i{___ \ ´ i|≧=-: :.}}: : : : : : : |: : :
___彡: : : :i{≠ぅミ k抖竺=ミ}}: : : : : : : |: : :
}: : .:}八{ {/{ハ. ⌒ {/{_い刈: : : : : : : |: : :
,: : : :乂__ Vrク 乂rク ノ :/.: : : : ; : : :
. ′: : : }}' ' ' ' ' 'ア: : : : : :_/_: : :
. i{: : :‘, 圦 r--‐ vi /}.: : : :/ }/ :
八: : ∧: : : . ノ __〉: : / ん: :
/ ):/ : }}:. : : : : . __〉: : :{ {/ : :
_/: : :.}} : : : : : :}>=ぉi爪 __ア. : :/} '{=-‐
⌒ア: : :}}: : : :j{ ̄ ̄ ,ハ ____ア゙.: :/ ,}i '{___
/: __.: :>‐j{ __/ }_ア゙.: :/__-=圦 /⌒`
┌────────────────────────────┐
│>>80 情報提供ありがとう! .|
│ふふん、CG法の基礎だけじゃなくてソースもあった方が良いよね?.|
└────────────────────────────┘
北大 センター開発・ライブラリプログラム
ttp://www.hucc.hokudai.ac.jp/library.html
ICCG法によるポアッソン方程式の数値解析コード-p.150
ttp://www.hucc.hokudai.ac.jp/pdf/library/kaihatsu_hokoku/k14.pdf
ttp://www.hucc.hokudai.ac.jp/library/psiccg.lzh
Matlab Central Conjugate Gradient Method
ttp://www.mathworks.com/matlabcentral/fileexchange/22494-conjugate-gradient-method
- 83 :
- なるほどね!今は無難に不完全コレスキーだから、ほかのも試してみようかしら
- 84 :
- は日曜日以外毎日8時から18時過ぎまで長期間に渡って周囲に騒音をまき散らさないように 配慮して工事する金儲けより公共の利益を重視するすばらしい会社です
- 85 :
- 朝日住宅は日曜日以外毎日8時から18時過ぎまで長期間に渡って周囲に騒音をまき散らさないように
配慮して工事する金儲けより公共の利益を重視するすばらしい会社です
- 86 :
- __ ,' な>>84
,.. '" ,:  ̄ ヽ、、! に>>85
,.' , /. / l |、 、 ゙'.|. 言
l ,' iーエ'、 ヽ;゙ー- l、:. | っ
l !//l:::::i l::::ヽヽ`i :、! て
lヘj 、::ノ 、::::リ | lナ| ん
〈 " ¨ " | レ゙| の
l>...._ ,.‐-、 _,,.ィi ⊂. ?
| :,イヽ`ニア:/iヽ' |
|./ ゝzヽ会't'' ヽ:i ヽ._/
/ ヽ :|
┌────────────────────────────┐
│ちなみに>>81の最後のリンクと>>82のリンクは初心者用です .|
│はかせの他のレスも初心者リンクがごっちゃになってるんでよろしく.|
│あと、国内大型計算機センターで北大みたいにソースを公開してる.|
│所は他には見つかりませんでした 残念 .|
└────────────────────────────┘
- 87 :
- Poisson方程式を周期境界条件で解きたいのですが,行列のrankが落ちてるので解けません.
どうしたらよいかご存知でしょうか.
- 88 :
- 残差切除法ってどうなんすかね?poisson方程式の解法として
- 89 :
- >>87 計算格子のどこか一点(四隅のどこかとか)を0とかに固定すればおk。
- 90 :
- ‖ l/ ノ / | / / / / l/l l l l l ト、
// / // / / l // / / / / /l ,,,xz斥 j j} l ; \
/ / / / / / /l // / / /l / / ,,x芹 ‖ j ‖ l / \
/ // イ {/ // / |___/ / / /// / ,;爻::::ノ /l/ / // / \
/ イ / イ / / l / ̄/ー----<_,,,// 癶;;゚/ メ`/ / / / } ヽ
/ /八  ̄ ̄ ゙、 イい -‐气宀ir====云=ミz / `~'~ / / l / / }
./ { \ ゙、 ヽ \ l/\ \:::::゚:::ノ / l:.. / / / / / /
{ 八 \ Y \ \ `'''冖''''" }ノ  ̄ / / / // ノ} /
\ \ i `ト...._三ニ= __ノ / / ///::://ーj/ _
\ \ ヽ l } l l `ト `ー-=xz-x〜スノ / /::l// l:::{/:::/::::::::::{::::::::::
`〜、 \_ Y j / l l \ `ト 、  ̄ ̄ / /:::::::l/:::l l:::::::::::::::::::::::::::::::::::::::
\ヽ l\ } j / ノj l \ ー- _/{八 /::::::::::{::::::l l::::::::::::::::::::::::::::::::::::::::::
┌────────────────────────────────────────┐
│>>87 周期境界なんだから入り口と出口の物理量が同じになるんでrankが落ちるのは当り前っス .|
│出口側を一旦切り棄てて行列を作り直してポアソン方程式を解いた後、出口側に入り口の値を .|
│挿入しなおせばいいんスよ 差分で解いてるんでしょ? .|
└────────────────────────────────────────┘
- 91 :
- / / / / /// / } } l l
/// / / / } / j ‖ / 八 } l
/ //l / / ノ,,,__ l | / / ∧ ゙、 j l
〜'" / l /l / /  ̄l\l / _ l\\ / l l
/ ̄ ̄ ̄l-┐/ / ,,,__ `トl l /'" /  ̄Tヽ / l l
ノ / j / { `宀〜ミx,,ヽ | l /l// l / / j }
(_/ / /癶 \ ヾ、 l /ノ彡示气,,_} / / } ハl j
〈_ イ `ー--― ''" ̄ ̄ ̄フ }/ んj!:::} )》=‐ / j / jノ
-――-、 ____ ---―'''" ノ:. 弋;;//`/ / / / /
 ̄ ̄ \い 、 ~"~ ////イ/ /\
_ -―――--┘ ト `ー-‐x,,,, ー=彡/ / \
 ̄ ̄ l l \{ ヽj \ ;;;; ̄ / ノ ノ 〕
\ l\ l r―l \ -‐''" / // / /
)ノト ノ ヽ从ハ |  ̄~''〜、_ー'"j┐ / //イ _彡イ/( /
)/ ノ }ノl ハ 斤 l 从{/ 八( __ \(
┌────────────────────────────────┐
│えー、フーリエ・スペクトル・コロケーション法の場合は .|
│周期境界条件は常に満たすとして出口側は常に切り捨てている状態っスね .|
│感覚的には端点の値は同じなのでどちらか不要、数値計算的には行列が .|
│非正則になるのでどちらかを切り捨てておく必要があるって所っスか? .|
│ .|
│フーリエ・スペクトル・ガラーキン法は単純にフーリエ係数との加減乗除になる.|
│んで、全方向周期境界条件の場合は本当に単純な計算になるっスよ .|
└────────────────────────────────┘
- 92 :
- ハカセじゃないと言葉が流暢だなぁ
なんだか面白いんでこのスレに張り付いてるんですが、
高卒程度の知識で楽しみながらやる勉強法ってありますか?
- 93 :
- なんだ、今期はリドルにハマってんのか
- 94 :
- /- '/\ { ̄.> 、\} .\
//}/, ,、ヽ >、 ‘,
i {/ / /' ,ム ヽ ヽ i
Y i i / {', \ ', ∧ i
/ i i{―、ヽ -‐、} i. ノ _
i i i .{代リ 弋リ jマ ヽ ヽ , -‐≦ ヽ
',. ', i { 、:: '_ _, :::j / \ ノ. '-´__/ ̄ ̄
マ.{、 ∨> 、 ./ノ /-、 .\ { f二_ 〉
ヽ \ r-、.ム ,ヒLノ イ / } 〉 ト ― j
ムノ Y/¨ j '/ 、'i i \ i \≧彡
┌──────────────────────────────────┐
│>>87 提唱者のHPがあるね 非圧縮性流体も解けるみたいだしCG法と比較して .|
│みたいじゃないか .|
│残差切除法による計算サービスサイト .|
│ttp://www.asiri.co.jp/product/rc_method/index.html .|
│だ円形境界値問題の数値解法―残差切除法について .|
│ttp://www.asiri.co.jp/product/rc_method/pdf/rc_ronbun.pdf .|
│非圧縮流れによる残差切除法を用いた高精度流れ解析法の検証 .|
│ttp://www2.nagare.or.jp/cfd/cfd26/paper/B10-1.pdf .|
│検証というならGhiaの論文と比較しないといかんでしょ? .|
│ .|
│>>87 どうにもならない場合はムーア・ペンローズの擬似逆行列という手も .|
│多軸自在継手を用いた可変構造の形態解析および一般逆行列に基づく流体解析.|
│ttps://dspace.jaist.ac.jp/dspace/bitstream/10119/11444/1/paper.pdf .|
│Matlab/Octaveのpinv関数だね 特異値分解できれば求められるじぇ .|
└──────────────────────────────────┘
- 95 :
- ,':. : : : : : : :.|: ; :,ィ: :/ !: :/ !: ハ: :l : : `、
l: : : : : : : : :|/!/ レ_≧、l∠ レ' _l;_l: : : : l
|: : : : : : : : :| //::::::「 テ=、:l: :j: : :!
|: : : : : : : : :l 〈 ハ:::::d /:::`:i } V:. : .}
l: : : : : : : : :l ,,, ヽ:ニ:ノ {.し:ソ ! ,イ: :N
│: : : : : i: : | `ー'゚ ,,/: :.N
!:: :. : : : l: : | ' |:: : |
! : : : : : l: : | ,. ‐‐- 、 l: : :|
, ‐‐-、__レヘ.:.、: : !: :.lト、 ` −‐' /:: : :.l
/: : _: : ∨ `' W V!:ノ、 丶 __ . . .ィ:´:: :: :: : :|
/: /:::::::〉: / \ |ニ‐‐、::;lヘ: : ;、 : 丿
{: :{:::::::/: | !,へ ヾ二_:ゝ l `´ `~
l: :l::::::!/\ \ / ヽつ |\ ,‐- 、
l: :ヽ:/ . . . . \ \‐-、 / ─孑ヽ.!. .⊂/,─ }
┌────────────────────────────┐
│>>92 >>5の初心者リンクの最初なんかは学部生が勉強する .|
│内容なのであまり関係無いですよ 計算手法を覚えて結果を .|
│どう評価するかの世界ですしね .|
│流体シミュレーションをやる場合、対象が非線形偏微分方程式の .|
│Navier-Stokes方程式になります いきなり非線形偏微分方程式 .|
│解くのは大変なので、>>5の最初のリンクで勉強しながら .|
│簡単な線形常微分方程式を解いて解析解と比較 .|
│簡単な非線形常微分方程式を解いて計算例と比較 .|
│簡単な線形偏微分方程式を解いて解析解と比較 .|
│簡単な非線形偏微分方程式を解いて計算例と比較 .|
│第一ゴール、流体シミュとしてキャビティ流れを解いて論文と比較 .|
│みたいな感じでやるのがいいかなと思ってます .|
└────────────────────────────┘
- 96 :
- /:.:./:.:./:.:.:.:.:.:.:.:.ィ:.∧:∧:.:.:.:.:.:.ヽ:.:.:.:.:.:.:.:.ヽ
. /:.:.:.:.:.:.:/:.:.:.:.:.:.:.:/:.:/ V:∧:.:.:.:.:.:∧:.:.:.:.:.:.:.:∧
,:.:.:./:.:.:./:.:.:.:.:.:./:.:/ V:∧:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.::::.
. .:.:./:./:./:.:.:.:.:/:;/ ∨:ハ:.:.:.:.:.:.:.V:.V:.:.:.:.:.:i
l:.:.:.:/:./:.ム彡 ' ゝー≠-弋:ハ:.:.:.:.:.:.:}:.:.}:.:.:.:.:.:.!
|:.:.:仏ァ -ぇ ぇ=z- ヽ:.ハ:.:.:.:.:.ハ:.ハ:.:.:.:.:.:|
|:.:.:.:イ,代歹 伐ツ 圦:リ:.:.:.:,:.:.∧:.:.:.l:!:.:|
l:.:/:/| ー彡 - 彡' |:.:.:.:/:./rハ:.:.:.リ:.:l
;.:.:.:.:} l:.:.:/:./ム ノ:.:.:.:.:::
V:.:.:; |:./:./r.イ:.:/:.:.:/
ヽ:.l ヽ レi;//:.:/:./:.:.; '
从 /:.:/:./:/
^ ⌒ イ:.:.:.:./'′
> < ト、:.イ 長道と光合成
ー┬一 リ V
「、| _ -=彡 l したい
| >zz< |
____ 入 || |_ =彡^ >
r// / ヽ || -一 >zz-
┌────────────────────────────┐
│>>93 誰が為に我は行く〜今期はシドニアの騎士とJOJOと蟲師 .|
│リドルとキャプテン・アース今のところ外れなんですがw .|
│ .|
│プログラミングでもツールでも色んな問題解くのが上達の早道 .|
│まずはサンプル的な問題から .|
│はっきりわかんだね .|
└────────────────────────────┘
- 97 :
- >>90
すんません
出口側を切り落として行列を作り直すとはどういうことでしょうか
- 98 :
- -──-
'" \ノ|
/ └- 、
/ /| ⌒ヽ
ー=彡 , / | i\ ,
/ / / 八 { \ ′
| / \>─\ト- i
| |´ハ| \ ,ハ |
| |/' 下示ミx | ト、 :|
|ハ _,,|__ ___ 乂):::リ | |ヽ} |
У _〉_〉_〉_〉 ̄ ,,,, | |ノノ |
ノ /´) :| }': |
厂{ _/ __ -‐ 、 ノ /::: |
/ \ /| ト __∨___ノ < ∨‐- 、 |
. / \|〃 /⌒ ⌒\/ , |
∧ ノ__ _\ i ,
. / \_//::::::/ l |::::::\ヽ | ′
{ / /::::::::/ ,__ , |:::::::::::| | | i
. \ :i /::::::::/ / /`ー──;l| |:::::::::::| | | |
ー┤  ̄ ̄ / | /´|:  ̄ ̄ | } |
廴___/`ー| | 乂____ノ /' |
┌───────────────────────┐
│>>97 おはようございます .|
│周期境界条件はどうやって満たしてるの? .|
│例えば差分法で端と端で周期境界条件を満たした場合.|
│そのままだとrankが落ちますけど これは分かる? .|
└───────────────────────┘
- 99 :
- >>98
内点をf1~fnとするとf0=fn, fn+1 = f1という境界です.行列でかくと
2 -1 -1
-1 2 -1
....
-1 2 -1
-1 2 -1
右上と左下の-1が周期境界に対応してます
- 100 :
- 行列わかりづらいですね.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}
です
100〜のスレッドの続きを読む
ANSYS
NullPointerExceptionを「ぬるぽ」と呼ぶスレ2
俺がandroidから書き込むときにAAを持ってくるスレ
移民を1000万人も受け入れたら日本はどうなるのか?
オープンソースCFDツールボックス OpenFOAM part3
CAE技術者の苦悩を愚痴れ
SNS狂鬼ムラカミ
俺がandroidから書き込むときにAAを持ってくるスレ
NullPointerExceptionを「ぬるぽ」と呼ぶスレ2
商用CFDについて語るスレ
--------------------
【殺人団体革マル枝野幸男】衆参同日選について「(可能性は)99%だと思っている。あるつもりで準備している」
[集会所] NMB48変態スレ☆2
快盗戦隊ルパンレンジャーVS警察戦隊パトレンジャー アンチスレ27
自殺を考えている既婚男性3
【AKB48】入山杏奈ちゃん応援スレ☆86【あんにん】
【ゲスト】合同誌にまつわる話81【アンソロ】
スタートの練習したら?
☆☆☆ 西内まりや総合本スレ ☆☆☆ Part.15
【スギ花粉】殺伐広場へR162周山街道【注意】49周目
多読のすゝめ 8冊目
【オリンピア】不二子 TYPE A+★9【獲得枚数変動】
シャロン・オズボーンの一日
【外科矯正】顎変形症・上顎前突・下顎後退【classII】 Part.2
夫の浮気を疑う中年女性の心理
レートは高いのにタイトル童貞のインチキ棋士豊島
Jリーガーブログ総合スレ Part2
【生物板ID導入議論スレ】part2
オリジナルサピエンスとネアンデルタールと現代人2
【農家と膿家は】農家の暗部 part151【違うんです】
信長の野望で過大、過小評価されている武将スレ19
TOP カテ一覧 スレ一覧 100〜終まで 2ch元 削除依頼