自作・フリーソフト流体シミュレーションスレPart3
at SIM
[1からを表示]
50:名無しさん@5周年
14/05/02 04:26:26.21 .net
非等方圧力で圧縮性流体(roeやhll)といてる先行研究とか知ってます?
51:DNS厨
14/05/04 01:16:59.53 .net
,. -─-、
/ ̄`ヽ { ヽ
{ ゝ-─-ゝ ノ
`7 ハヽ ハヽ <
/ しノ しノ ヽ
│ ,.-‐(^>─-、 |
│ { _人__ \ ノ
'、 くマ{ V7 ̄ノ /
\ \ゝー ´_/ /、
/ ヘ._ ̄ ̄___, イ| \
. / ヽ__\ ̄ /__」/ / ̄`ヽ
/\ <.__\/_」 { { }
. ( \,.ィ  ̄ ̄ \ゝ.___,ノ
`ー一'′| o `|
| o |
┌────────────────┐
│>>50 前スレ259の人かな? .|
│A New High Resolution Scheme for Compressible Viscous Flow with Shocks .|
│URLリンク(aero-comlab.stanford.edu) .|
│圧縮性流れのアルゴリズム(II) 離散化方程式の解法 .|
│URLリンク(www.ifs.tohoku.ac.jp) .|
│↑こっちは先行研究とかじゃないけど国内の研究 .|
└────────────────┘
52: 忍法帖【Lv=40,xxxPT】(1+0:8)
14/05/04 23:21:25.21 .net
寧ろ逆に天国だろ
他に別に極楽だろ
当然正反対に理想郷だろ
53:名無しさん@5周年
14/05/09 02:28:23.74 .net
>>51
ありがとうございます
でもこれらは非等方圧力ではないですよね?
54:DNS厨
14/05/09 23:48:36.15 .net
__ ,' な
,.. '" ,:  ̄ ヽ、、! に
,.' , /. / l |、 、 ゙'.|. 言
l ,' iーエ'、 ヽ;゙ー- l、:. | っ
l !//l:::::i l::::ヽヽ`i :、! て
lヘj 、::ノ 、::::リ | lナ| ん
〈 " ¨ " | レ゙| の
l>...._ ,.‐-、 _,,.ィi ⊂. ?
| :,イヽ`ニア:/iヽ' |
|./ ゝzヽ会't'' ヽ:i ヽ._/
/ ヽ :|
┌───────────┐
│>>53 非等方性流れ≒非等方圧力 じゃないの? .|
│ .|
│非等方圧力をどういう意味で使ってるのかと .|
│どういう流れ場を想定しているか教えて欲しいけど.|
└───────────┘
55:名無しさん@5周年
14/05/10 01:11:47.05 .net
>>54
うちの業界だけの言葉なんですかね?
圧力テンソルの対角成分が同じものでないことです。
例えば二次元なら
∇・P = (∂x Px) ex + (∂y Py) ey
56:DNS厨
14/05/10 06:56:21.82 .net
. - ─ - .
/ ` .ノ}
. . ' ト、
/ / ∧. ∧ 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:名無しさん@5周年
14/05/10 10:19:13.53 .net
>>56
独立な圧力がP一つではなくてPxとPyなんですよ
58:名無しさん@5周年
14/05/10 10:22:08.14 .net
流体中での爆発とか?
59:名無しさん@5周年
14/05/10 10:25:42.51 .net
それもあるかもしれないですね。
天体のジェットなんかを考える時はジェットの吹き出る方向と、それに垂直な方向では温度が違うんですよ
60:名無しさん@5周年
14/05/10 10:27:31.04 .net
P = P⊥ diag(1, 1, 0) + P|| diag(0, 0, 1)
なんてかいたりしますね
61:名無しさん@5周年
14/05/10 10:29:02.72 .net
通常の流体(等方圧力)では
P = P diag(1, 1, 1)なわけです
62:名無しさん@5周年
14/05/10 15:25:31.05 .net
分子運動論的な言い方をすると、方向ごとに熱速度の大きさが違うので
面の方向ごとに圧力が違うってこと?
そういう世界もあるんだなあ。状態方程式も圧力の複数の成分が入ってるのかな
63:名無しさん@5周年
14/05/11 00:50:24.75 .net
物理現象としては同じものを対象としてないって理解でOK?
64:名無しさん@5周年
14/05/11 01:49:01.72 .net
>>63
同じものを対象にしていないとはどういう意味でしょうか。
URLリンク(www.physics.wisc.edu)
みたいな感じですね。これは磁場が入ってますけど
65:DNS厨
14/05/11 22:53:30.30 .net
,、- '''"~~~~゙゙゙ ''‐ 、
,、.'" __,,,,,_ ゙' 、 / ヽヽ / /
,.':: ,、- ''"´ `゙゙゙''‐ 、 ゙' , ,,、、、 /_ \ / / /
,.'::: .,. " ゙'.、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プラズマの圧力非等方度同定手法の研究
URLリンク(ir.soken.ac.jp)
Stability of anisotropic-pressure tokamak equilibria to ideal ballooning modes
URLリンク(research.microsoft.com)
Magnetic hydrodynamics of two-dimensional flows in a plasma with an anisotropic pressure
URLリンク(link.springer.com)
2.ミラープラズマの安定性
URLリンク(jasosx.ils.uec.ac.jp)
>>64 電磁流体力学(MHD)のCGL理論では非等方圧力テンソルは磁場に垂直な圧力と
磁場に平行な圧力で与えられる...ちぃ覚えた
66:DNS厨
14/05/11 23:07:24.99 .net
' /: /: : : /: : ィ≠^!ミ: :: : : : : \__ノ〉
,' /: ://ハ: : : //´ |! ヽ: : : : : : : イ
|.::': ://: : : / }ツ 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 キーワードで絞り込めば逝けるんじゃないかい? .|
└─────────────┘
MAGNETIC FIELD AMPLIFICATION AND EVOLUTION IN TURBULENT
COLLISIONLESS MHD:AN APPLICATION TO THE ICM
URLリンク(arxiv.org)
Numerical Prediction of Three-Dimensional Non-Equilibrium Flows
Using the Regularized Gaussian Moment Closure
URLリンク(arrow.utias.utoronto.ca)
67:名無しさん@5周年
14/05/20 01:03:01.31 .net
初心者ですが、二流体モデルをやるにはどうすればよいんでしょうか?
68:名無しさん@5周年
14/05/20 04:53:37.90 .net
>>67
collisionいれればいいんじゃない?
69:DNS厨
14/05/21 00:56:07.65 .net
-‐=‐- ミ
. 7ア⌒`__'⌒` 、
. ,.イ.: :7:. :.`:: ..、 V)__
゙`ー=ァ´.: /.: :/.: :. :.、:. :. ゙'く⌒`
/ /.:: .:/.: :/.: ハ:. :. :.\:. :\__)
. {_..イ.:/.: .:/.: :/.: / jト、:. :. :. \:. :.〈
. /.: :l.: .: .: .: .: {/- 、 Χハ:. :. :ヽ:. :.
/.: .: |.:{.: .:.{.:: 仏=ミ ィf〔7:. :. :. }:. :!
. {.: .: :八:. :. :. :. 「V_ソ 弋ソ ハ:.. :.八_{
八: /.: .:Vハ:. :.|`` ' `八! :/
/.: :イ.: .: :: :: :V 个ト . マフ ,.く
. //"j|.: .: :: :: :/〉八 〈〕T爪 :.
/゙ |.: .: :: ::/`ヽ〉:.\__八__) :i
|.:: ::. ::.{ V〉:. 「r-<) 、 |
┌──────────────────┐
│>>67 流体と粒子の混合物の場合の論文 2章あたりを読んでみてね .|
│2流体モデルによる流動層の理論解析と数値実験 .|
│URLリンク(hdl.handle.net) .|
│URLリンク(repo.biblio.tuat.ac.jp) .|
│モデルによって必要なものが変わる上に実験や経験式が必要な印象なのだけど... .|
│使ってる人の意見を求めてみようじゃないか! .|
└──────────────────┘
70:DNS厨
14/05/21 01:07:36.98 .net
.. 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:名無しさん@5周年
14/05/21 20:35:57.04 .net
とりあえず、流体がシミュレーション通りに動くことを己の目で見たいのだけど
どこに行けばいい?
72:DNS厨
14/05/22 06:54:19.10 .net
/ ` ヽ、
(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:名無しさん@5周年
14/05/22 12:31:38.82 .net
>>68
>>68
ひとまずもう少し勉強してみます。ちなみに液液です
74:DNS厨
14/05/27 01:01:25.75 .net
', / 単 二
固 気 .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 .|
│URLリンク(www.eng.tau.ac.il) .|
│これでもダメなら .|
│"two liquid phases" "Two Fluid Model" とか .|
│"liquid-liquid" "Two Fluid Model" でググるとか? .|
└───────────┘
75:名無しさん@5周年
14/05/29 14:41:12.86 .net
fortranのmpi fftwのサンプルコードとか持ってる人いるかな.
fftwはドキュメントが不親切過ぎて...
76:DNS厨
14/05/30 06:23:08.28 .net
(⌒ヽ (⌒ヽ / ∧ ⌒ヽ、
` ` ´ / /! / ! ヽ ヽ ̄  ̄_)  ̄_)
/ / / { / 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 .|
│URLリンク(www.fftw.org) .|
│これは確かに厳しいかもしれない .|
│FFTW MPI 使ってるフリーソフトがあれば参考になるんだけど .|
└─────────────┘
77:名無しさん@5周年
14/05/30 18:38:11.09 .net
>>76
そしたらみなさんFFTってMPIでどうやってるんですかね?
ライブラリ使わず??
FFTW以外のライブラリを使う??
78:名無しさん@5周年
14/05/31 01:28:36.60 .net
そもそも今どきMPIを書くのか
大筋は↓で C から F への読み替えで一回ややこしくなっただけじゃないか、と大雑把に
変換後の値は各プロセスにばらまかれてるから集めるなら MPI API コールかね。
マルチプロセスでスタート、MPI と FFTW の INIT
↓
横に何人いるかは分かってます。
全体でLxMのFFTやりたいんで自分の割り当て量ください。
↓
各自、自分のメモリ取ります。
↓
MPI-FFTプラン立てます。
↓
計算します。
79:DNS厨
14/05/31 09:06:08.18 .net
/ ,:イ>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 .|
│URLリンク(www.fftw.org) .|
│マルチスレッドFFTWもCのサンプルから調査したのすっかり忘れてた.|
│ .|
│FFTのテストは三角関数辺りから始めて複雑なものはOctaveや .|
│Scilabなどの組み込み系と比較したな .|
└───────────────┘
80:名無しさん@5周年
14/06/07 16:57:23.82 .net
前処理付き共役勾配法が今のところ最強ってことでいい?
81:DNS厨
14/06/10 00:19:35.39 .net
/ ` ヽ、
(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法の最近の前処理のロバスト性と効率化について
URLリンク(www.kurims.kyoto-u.ac.jp)
並列非線型共役勾配法アルゴリズムとその性能評価
URLリンク(www.ssisc.org)
数値計算の基礎 共役勾配法の勉強
URLリンク(www7.ocn.ne.jp)
82:DNS厨
14/06/10 21:51:06.20 .net
/⌒)__/⌒\: : : ``丶、__´'くi:i:i:_〉
. .′/⌒{\: : :.\ : : : : : \``丶、
{ ./.: : :从: :丶: : ∨/: : : : : : : : : :\
{/.: : :/.{ \: : \∨/ : : : : : : : : : : \
. ___/.: : :,′{ \ : :`∨\: : : }}: : : : : : : \
. ⌒,′: .:′、 \:__|=- ≧=--‐: : : : :|:.\
′: : i{___ \ ´ i|≧=-: :.}}: : : : : : : |: : :
___彡: : : :i{≠ぅミ k抖竺=ミ}}: : : : : : : |: : :
}: : .:}八{ {/{ハ. ⌒ {/{_い刈: : : : : : : |: : :
,: : : :乂__ Vrク 乂rク ノ :/.: : : : ; : : :
. ′: : : }}' ' ' ' ' 'ア: : : : : :_/_: : :
. i{: : :‘, 圦 r--‐ vi /}.: : : :/ }/ :
八: : ∧: : : . ノ __〉: : / ん: :
/ ):/ : }}:. : : : : . __〉: : :{ {/ : :
_/: : :.}} : : : : : :}>=ぉi爪 __ア. : :/} '{=-‐
⌒ア: : :}}: : : :j{ ̄ ̄ ,ハ ____ア゙.: :/ ,}i '{___
/: __.: :>‐j{ __/ }_ア゙.: :/__-=圦 /⌒`
┌──────────────┐
│>>80 情報提供ありがとう! .|
│ふふん、CG法の基礎だけじゃなくてソースもあった方が良いよね?.|
└──────────────┘
北大 センター開発・ライブラリプログラム
URLリンク(www.hucc.hokudai.ac.jp)
ICCG法によるポアッソン方程式の数値解析コード-p.150
URLリンク(www.hucc.hokudai.ac.jp)
URLリンク(www.hucc.hokudai.ac.jp)
Matlab Central Conjugate Gradient Method
URLリンク(www.mathworks.com)
83:名無しさん@5周年
14/06/11 02:52:05.71 .net
なるほどね!今は無難に不完全コレスキーだから、ほかのも試してみようかしら
84:名無しさん@5周年
14/06/11 15:29:04.25 .net
は日曜日以外毎日8時から18時過ぎまで長期間に渡って周囲に騒音をまき散らさないように 配慮して工事する金儲けより公共の利益を重視するすばらしい会社です
85:名無しさん@5周年
14/06/11 16:40:17.83 .net
朝日住宅は日曜日以外毎日8時から18時過ぎまで長期間に渡って周囲に騒音をまき散らさないように
配慮して工事する金儲けより公共の利益を重視するすばらしい会社です
86:DNS厨
14/06/13 00:49:15.60 .net
__ ,' な>>84
,.. '" ,:  ̄ ヽ、、! に>>85
,.' , /. / l |、 、 ゙'.|. 言
l ,' iーエ'、 ヽ;゙ー- l、:. | っ
l !//l:::::i l::::ヽヽ`i :、! て
lヘj 、::ノ 、::::リ | lナ| ん
〈 " ¨ " | レ゙| の
l>...._ ,.‐-、 _,,.ィi ⊂. ?
| :,イヽ`ニア:/iヽ' |
|./ ゝzヽ会't'' ヽ:i ヽ._/
/ ヽ :|
┌──────────────┐
│ちなみに>>81の最後のリンクと>>82のリンクは初心者用です .|
│はかせの他のレスも初心者リンクがごっちゃになってるんでよろしく.|
│あと、国内大型計算機センターで北大みたいにソースを公開してる.|
│所は他には見つかりませんでした 残念 .|
└──────────────┘
87:名無しさん@5周年
14/06/16 00:27:44.16 .net
Poisson方程式を周期境界条件で解きたいのですが,行列のrankが落ちてるので解けません.
どうしたらよいかご存知でしょうか.
88:名無しさん@5周年
14/06/16 23:58:13.74 .net
残差切除法ってどうなんすかね?poisson方程式の解法として
89:名無しさん@5周年
14/06/17 00:33:56.75 .net
>>87 計算格子のどこか一点(四隅のどこかとか)を0とかに固定すればおk。
90:DNS厨
14/06/17 19:31:54.49 .net
‖ 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:DNS厨
14/06/17 19:43:45.09 .net
/ / / / /// / } } 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:名無しさん@5周年
14/06/17 22:59:47.33 .net
ハカセじゃないと言葉が流暢だなぁ
なんだか面白いんでこのスレに張り付いてるんですが、
高卒程度の知識で楽しみながらやる勉強法ってありますか?
93:名無しさん@5周年
14/06/17 23:51:45.78 .net
なんだ、今期はリドルにハマってんのか
94:DNS厨
14/06/17 23:58:15.17 .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 \≧彡
┌─────────────────┐
│>>87 提唱者のHPがあるね 非圧縮性流体も解けるみたいだしCG法と比較して .|
│みたいじゃないか .|
│残差切除法による計算サービスサイト .|
│URLリンク(www.asiri.co.jp) .|
│だ円形境界値問題の数値解法―残差切除法について .|
│URLリンク(www.asiri.co.jp) .|
│非圧縮流れによる残差切除法を用いた高精度流れ解析法の検証 .|
│URLリンク(www2.nagare.or.jp) .|
│検証というならGhiaの論文と比較しないといかんでしょ? .|
│ .|
│>>87 どうにもならない場合はムーア・ペンローズの擬似逆行列という手も .|
│多軸自在継手を用いた可変構造の形態解析および一般逆行列に基づく流体解析.|
│URLリンク(dspace.jaist.ac.jp) .|
│Matlab/Octaveのpinv関数だね 特異値分解できれば求められるじぇ .|
└─────────────────┘
95:DNS厨
14/06/18 00:47:01.84 .net
,':. : : : : : : :.|: ; :,ィ: :/ !: :/ !: ハ: :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:DNS厨
14/06/18 00:55:37.83 .net
/:.:./:.:./:.:.:.:.:.:.:.:.ィ:.∧:∧:.:.:.:.:.:.ヽ:.:.:.:.:.:.:.:.ヽ
. /:.:.:.:.:.:.:/:.:.:.:.:.:.:.:/:.:/ 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:名無しさん@5周年
14/06/18 03:21:14.43 .net
>>90
すんません
出口側を切り落として行列を作り直すとはどういうことでしょうか
98:DNS厨
14/06/18 07:04:46.92 .net
-─-
'" \ノ|
/ └- 、
/ /| ⌒ヽ
ー=彡 , / | i\ ,
/ / / 八 { \ ′
| / \>─\ト- i
| |´ハ| \ ,ハ |
| |/' 下示ミx | ト、 :|
|ハ _,,|__ ___ 乂):::リ | |ヽ} |
У _〉_〉_〉_〉 ̄ ,,,, | |ノノ |
ノ /´) :| }': |
厂{ _/ __ -‐ 、 ノ /::: |
/ \ /| ト __∨___ノ < ∨‐- 、 |
. / \|〃 /⌒ ⌒\/ , |
∧ ノ__ _\ i ,
. / \_//::::::/ l |::::::\ヽ | ′
{ / /::::::::/ ,__ , |:::::::::::| | | i
. \ :i /::::::::/ / /`ー─;l| |:::::::::::| | | |
ー┤  ̄ ̄ / | /´|:  ̄ ̄ | } |
廴___/`ー| | 乂____ノ /' |
┌────────────┐
│>>97 おはようございます .|
│周期境界条件はどうやって満たしてるの? .|
│例えば差分法で端と端で周期境界条件を満たした場合.|
│そのままだとrankが落ちますけど これは分かる? .|
└────────────┘
99:名無しさん@5周年
14/06/19 00:05:19.04 .net
>>98
内点をf1~fnとするとf0=fn, fn+1 = f1という境界です.行列でかくと
2 -1 -1
-1 2 -1
....
-1 2 -1
-1 2 -1
右上と左下の-1が周期境界に対応してます
100:名無しさん@5周年
14/06/19 00:09:57.84 .net
行列わかりづらいですね.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}
です
101:DNS厨
14/06/19 07:05:47.62 .net
___
/ `丶、
」∨ ヽ
7 '.
/ ノ| ∧ `ト N| っ
. /N イ/ |/ ∨| ∧ | っ
. 〃 Y⌒i W∩ ∩ }/V -=ニ二三三三三
/ {{{ | | ∪ ∪ {{-=ニ三/ ̄ ̄`'く>三三三
/ゝ _| | '' '' 人.三三|(‐〜〜‐)| |三三三
/ 厶-=ニ | :八‐-(_:::::::⌒)<∧ノ三三| ー--‐' |_|三三三
-=ニ二三三三三\_,,,フ三三三三三三三三| <・ ) 「 三三三
二三三三三三三三三三三三三三三三三三ニ圦 {__>,,ノ三三三
三三三三三三三三三三三三三三三三三三三ニ≧===≦三三三
┌────────────┐
│>>100 ごめんなさい 昨日は寝不足で寝ちゃった>< .|
│1次元拡散方程式で周期境界条件、ソース付き .|
│URLリンク(rnumata.org) .|
│参考になると思うけど .|
│ .|
│rankが落ちるということは線形独立じゃない方程式が .|
│存在するって事なので式の数減らせば?って事だっ .|
│たんだけど そんな事しなくても解けるねw .|
└────────────┘
102:名無しさん@5周年
14/06/20 04:34:21.33 .net
>>101
これは拡散方程式ですよね。しかも陽解法で解いてます。
ちょっと計算してみましたが仮に拡散方程式を陰解法で解いても行列はランク落ちしないです。
単純に考えてみると、周期境界のポアソン方程式は定数分の不定があるからランク落ちすると思うんですが、拡散方程式は方物型だから、コーシーデータが必要になるので、定数分の不定もなくなりますよね
103:DNS厨
14/06/20 20:54:31.52 .net
./ ./: フ´: : : : : /: ̄\ \ ∧
' /: /: : : /: : ィ≠^!ミ: :: : : : : \__ノ〉
,' /: ://ハ: : : //´ |! ヽ: : : : : : : イ
|.::': ://: : : / }ツ l! \: : : : :トヽ
l:::!://: : : ,' / |! i: : : ::| ヾ!
. ',:| /.: : : :|::| 二二ゞ′/ ヽ- ―ハ: : : :| !l
、_Y/: : : : !::| ./ }: : :! !|
 ̄ 7: : : ハ个j 〓〓 .〓〓 /: : :/ リ
|ハ: : |: :ヾ:l /: : /、__ノ
|!:レヽ!: : ヽ> ' /: //! ̄
|! ヾ:ヾ个 、 .⊂⊃ イ: /ハノ
し }: ::|ヾ人` ‐z‐ '7ハ彡ソ| リ
|: 」|: :|辷_7∧/ |: : ::! ̄ ̄ ̄ヽ
|: :|!: :| / ∨__入_.!: : :ヽ、 )、
./!: :j: ::ヽ/、/レ' ̄三 7ヾ: :ヽ__ \
./ |::/ 八:_::>/ハ \ 〈 T  ̄ヽ.. \
广ヽ i/::/ L _ハ!ヽj ≫'!: ! \ \! \
| ∨::/\ ハ ん_ノ ./ | | ヽ ',. \
┌──────────────┐
│ふふん↓ようやくゆっくりレスできるじゃないか .|
│>>100 放物型と楕円型も確認しないで>>100のようなレスを .|
│するとは己も焼きが回ったな〜w .|
│周期境界条件のポアソン方程式の解法については答えはもう .|
│出ている >>91のフーリエ・ガラーキン法なら簡単に解けるじぇ . |
│問題があるとすると、二次元の時は(0,0)モード、三次元の時は .|
│(0,0,0)モードに解の空間平均値を代入することかな? .|
│ .|
│それはそれとして>>100のリンクは>>95の例として勉強になるね .|
└──────────────┘
104:名無し
14/06/21 00:03:48.18 l1sgvTEJj
初めまして、流体の勉強してるんですけど、ここで質問いいですかね。
円管内層流の熱流体のシミュレーションを行いたいんですけど、
円管外の熱が時変の場合のシミュレーションって複雑になるんですか?
出力は円管内の流体の温度です
勉強し始めたばっかでよくわからないんですけど、この場合のプログラムってどうなるんですかね。
MATLABやSCILABで行う予定です
105:名無しさん@5周年
14/06/21 00:16:21.99 .net
>>103
フーリエで解けるのはわかるんです。0モードを考えないから
僕は差分でときたいんですねー
106:DNS厨
14/06/21 11:12:50.42 .net
__ -―- _
x <´ ,、 .,、_, >=====、
/ ,r、__ノ::Y:::::|、,、 //‐-ヽ, ヽ
/ ,f-―v-'⌒`/i-.ヽ∧:::} }
/l ,r.-< ̄./ ./ \_フ ./
/ i // //,イ /, } .}\
{ ,.イ { .// /.ヽ/ iリ ト,} .ト, ヽ
{ / _,l .{ {,/ ヽ .} } __, .} lヽ .}
\ / fr } iリ 〈fテ't、 /ノ´リ} .ト,}
__\// ヽ_} l `¨´ ,rチ }//}
_// ̄`ヽ<_,ノl ', '_ ,¨./},イリ
/ ./ .} / l l ト,.〈〜ニ= ,イ,イ.i/
x <´, / } /、 l l i >- ,r ' } .i }リ'
x <´ l .\{.l l >,\リ l,/ }リ
/ 、 _ // `゙フ `i!,/<::::i>`ヽ ./,|/i
./ .`i } ー->ー-r一’ , ヽ, /::::::::::\ .} .′
ヽ i/´¨/. / ./ .Y::::::::::::::::::\
┌─────────────┐
│差分法で解きたいというならそれももう出てる .|
│ムーア・ペンローズの擬似逆行列使えばいいじゃないか? .|
│ .|
│周期境界条件で0モードも0でいいからFFTを使う方法= .|
│フーリエ・ガラーキンで無問題で差分法にこだわる必要無い.|
│んじゃないかい? .|
│発動条件は厳しいが、これ精度も良いし計算自体は単純 .|
│だしで代替の必要感じないんだがなぁw .|
└─────────────┘
107:名無しさん@5周年
14/06/22 02:57:14.13 .net
ひきこもり30代後半で童貞で低収入イラストレーターやってるバカ発見。
足立区に住んでいるそうだ
URLリンク(inumenken.blog.jp)
108:名無しさん@5周年
14/06/25 01:35:50.45 .net
>>106
擬似逆行列ですか!勉強してみます
ちょっとおもったんですが、平均値がゼロという条件を課せばランク落ちしないんですが、帯行列でなくなるために随分ときにくくなって却下しました
109:名無しさん@5周年
14/06/25 11:34:27.19 .net
◎2chスレッド勢いランキングサイトリスト◎
★+ニュース板
・ 2NN (推奨サイト)
・ 2chTimes
★+ニュース板新着
・ 2NN新着
・ Headline BBY
・ unker Headline
★+ニュース板その他
・ Desktop2ch
・ 記者別一覧
★全板
・ 全板縦断勢いランキング (推奨サイト)
・ スレッドランキング総合ランキング
・ ログ速
★全板実況込み
・ 2勢 (推奨サイト)
・ READ2CH
・ i-ikioi
※ 要タイトル検索
※ 2chブラウザ併用推奨
110:東日本大震災
14/06/25 16:46:41.07 .net
あんなー目撃証言0とかでこういう事をするとなー東日本大震災とかになるだわ。
理解したー?
26才を超えてるのに精神病院がどうのこうのとか。
俺は小学生じゃないぞ。
これだけ年を取ってるのにお前はいったいなんなんだ?
友達じゃねーんだ。
111:東日本大震災
14/06/25 16:50:34.60 .net
あーはいそうですかー?
いちいちその時にそうやって祈っててどうもすいませんでした。
祈ってないのになー。
112:天皇陛下
14/06/25 16:52:02.83 .net
新聞配達員の最中に食べに行った刈谷市の辺のファミリーレストランを解らなくして東日本大震災
113:DNS厨
14/06/26 07:08:10.74 .net
- ─ -
-‐ァ'´ ` 、
. / / ヽ
.′/ /! ∧ ト、 、 ヘ
| i / / j / ! j ! } ',
| l! / / /イ |/ j八 } }
| 八 /}/ N N }八乂
| Y^Y| |:| |:| ハノ
| {( || U U {ノ
| 乂||x:x x:x〉
| | ||ト.丶. _ ノ ,.ィ7
| jl<{ {` マニニフ ´}ト、
| / > \/{ }V < \
┌───────────┐
│>>108 無数にある解の内の一つを何らかの条件を .|
│付けて決定すれば良いので... .|
│>>89 みたいな方法でもいけるかもしれない .|
│これって全境界をノイマン境界条件で解いた時に .|
│やってたのを見たことがあるけど .|
└───────────┘
次ページ最新レス表示スレッドの検索類似スレ一覧話題のニュースおまかせリスト▼オプションを表示暇つぶし2ch
1808日前に更新/315 KB
担当:undef