流体力学のシミュレーションソフト

誰でも使えていろいろ遊べる安いソフトってありませんか？
渦巻きを作ってその減衰を見てみたいんです。
よろしくお願いします。

自分で作ればよいのでは?
製品は高いですよ

僕も有限体積法を使った流体解析をやりたいのですが、
なかなか難しくて作れません。
周りに手本にできそうな人もいないし、
いいホームページありませんか。

ん〜どうでしょう、ちょっとお聞きしたいのですが
定常1次元の熱伝導を解ける位の知識はお持ちなんでしょうか?
それとも、プログラムに対するアレルギーでも?
まず簡単なこのプログラムでも組んでみてはどうでしょうか

簡単な課題を出してもらえませんか？３時間以内程度で
できる（２さんみたいなプロなら１０分くらい？）やつ。
ちょっとウォーミングアップをしてから話をしたい。。
できれば明確な問題設定とともにお願いします。m(==)m

プロではありませんが･･
それではこういうのはどうでしょう
熱伝導率k=1の棒があります棒の一端は100℃ もう一方は0℃
定常状態の温度分布はどうなるでしょう?

基礎方程式は(Tは温度,xは棒の長手方向の座標,y,z方向無視)
k∂＾2T/∂x＾2=0
境界条件は
T(1)=100, T(10)=0
で格子の数は10個と決めます(とりあえず)

続き (上の条件で組めるなら見る必要なしです)
これを、中心差分で差分すると、ある任意の座標点iに対して
(i=2〜10 1と11はもう与えてしまったから)
(ごめんなさい格子数は11にすべきでした、修正お願いします)

{T(i-1)+2*T(i)+T(i+1)}/(�凅^2)=0
�凅は例えば 1/(11-1) =0.1 としましょうか
この11は格子の数です

すると、11元連立方程式のできあがりです。
後はこれを解くだけ。
ここまで良いですか?

ちなみに理論解は1次関数です(お分かりだと思いますが)
中心差分について説明が必要なら言ってください
できる限りお答えします

むー、ありがとうございます。では手元のC++で取り組んでみます。
不慣れなもので何から手をつけたらいいのかわかんないですが、
ゆっくり考えてみます。では後ほど（笑）。

差分の３点公式というやつでは？たしか寺寛の応用編で見たことある。

>10
そうなんですか? 私は微分の定義を2回適用して
limitを外した物とおもってたんですが
御存知なら詳しい話を教えてください。

寺寛の応用編ではhを十分小さな正数f(x)を微分可能な
関数で有れば微分商の定義は以下のように定義できる
[f(x+h)-f(x)]/h=f'(x)
これを前進差分商と後退差分商とをTaylorの公式を適用する
ことで上記の式を求めています

ごめんなさい!!!
式おもいっきり間違ってます
差分式は
↓ここです
{T(i-1)-2*T(i)+T(i+1)}/(�凅^2)=0
です。

6さん他 コードを組んだ方 本当に申し訳ありませんでした

差分方での熱流体解析ソルバー、ノード数限定であれば、フリーで入手できます。
宇宙業界のソフトで　SINDA/FLUINTっていうやつ。
二相流も扱えるらしい。

テーラー展開の公式は
f(x+d) = f(x) + d f'(x) + d^2/2 f''(x) + d^3/3! f'''(x) + d^4/4! f''''(x) +

これを利用して
f(x+2d) = f + 2d f' + 2d^2 f'' + (4/3)d^3 f''' + (2/3)d^4 f'''' - …
f(x+d) = f + d f' + d^2/2 f'' + d^3/6 f''' + d^4/24 f'''' + …
f(x) = f
f(x-d) = f - d f' + d^2/2 f'' - d^3/6 f''' + d^4/24 f'''' + …
f(x-2d) = f - 2d f' + 2d^2 f'' - (4/3)d^3 f''' + (2/3)d^4 f'''' - …

さて３点公式は
{f(x+d) - 2 f(x) + f(x-d)} = d^2 f'' + d^4/12 f'''' + d^6/360 f'''''' + …
となる。

f(x-2d)〜f(x+2d) を使って、f''''の係数も消してしまうことを考える。
そのためには
{f(x+2d) - 16 f(x+d) + 30 f(x) - 16 f(x-d) + f(x-2d)}
= - 12 d^2 f'' + (2/3)d^6 f'''''' + …
これは５点公式とよばれ f''についての３点公式よりも精密な近似となる。

同様にしてf(x-3d)〜f(x+3d)を使う７点公式を作れば d^6f/dx^6 までを消すことができ
f(x-4d)〜f(x+4d)を使う９点公式を作れば d^8f/dx^8 までを消すことができる。

こんな雰囲気でよかったと思います。計算間違いはあるかもしれません。

T''=0
T(1)=100
T(11)=0
理論解は
T = 110 - 10 x
ですね。

格子の数が11で、
２階微分を
{T(i-1)-2*T(i)+T(i+1)}/(�凅^2)
で近似する。

点2について、i = 2を代入して
{T(1)-2*T(2)+T(3)}/0.01 = 0
整理して
T(1) - 2 T(2) + T(3) = 0

結局、解くべき連立方程式は
T(1) = 100
T(1) - 2 T(2) + T(3) = 0
T(2) - 2 T(3) + T(4) = 0
T(3) - 2 T(4) + T(5) = 0
T(4) - 2 T(5) + T(6) = 0
T(5) - 2 T(6) + T(7) = 0
T(6) - 2 T(7) + T(8) = 0
T(7) - 2 T(8) + T(9) = 0
T(8) - 2 T(9) + T(10) = 0
T(9) - 2 T(10) + T(11) = 0
T(11) = 0

ちょっと覗いてみては？
http://www.cfd-online.com/Resources/

行列の対角化ですね。どうするんだっけ？(爆)

NS方程式による流体解析プログラムを作ったのですが、数値的に正確な解が得られてるのか確認したいのですが、ポアズイユ流れ以外に何かありませんか？
どこかにベンチマークになるような資料はないでしょうか？

CIP法というのがあります。
http://ciprus2.es.titech.ac.jp/index.html

解析解が在るわけではないですが、GhiaらによるCavity flowが
よく用いられます。＞NSソルヴァーのベンチマーク
U. Ghia, K.N. Ghia, C.T. Shin, High-re solution for incompressible
flow 〜, J. Comput. Phys. 48 (1982) 387-411.

#include <stdio.h>
#define N 11

void main() {
double ma[N][N], mb[N][N];
int n=N;
double t, s;
int i, j, k;

/// 係数行列の初期値を設定 ///
for (i=0; i<n; i++)
for (j=0; j<n; j++)
ma[i][j] = mb[i][j] = 0;
for (i=0; i<n; i++) mb[i][i] = 1;
ma[0][0] = ma[n-1][n-1] = 1;
for (i=1; i<n-1; i++) {
ma[i][i-1] = ma[i][i+1] = 1;
ma[i][i] = -2;
}

/// maの逆行列がmbに入る ///
for (j=0; j<n; j++) {
t = ma[j][j];
for (k=0; k<n; k++) {
ma[j][k] /= t;
mb[j][k] /= t;
}
for (i=0; i<n; i++) {
if (i!=j) {
s = ma[i][j];
for (k=0; k<n; k++) {
ma[i][k] -= ma[j][k] * s;
mb[i][k] -= mb[j][k] * s;
}}}
}

for (i=0; i<n; i++) printf("%f\n", mb[i][0] * 100.0);
}

つづきをどうぞ＞２

熱伝導方程式っぽいのは？

double ma[N][N], mb[N][N];
こんなのは、普通スタックには置かんだろ。

スタック？

mainの上に放り出すか、
double x[][]; ---&gt; static double x[][];
くらいにはするだろうっちゅうこった。

どーゆー利点があるんだい？

というか、結局は本体をmainでなく外部関数にしろってことでしょ。
まあ良いんじゃない。簡単の為に上のようにしたんだろうから。

ちがうでしょ。効率のもんだいでしょ。

class matrix {..};
matrix a = new matrix(n);
　　:
a.inverse();
みたいなのをもういじる必要のない形で
ライブラリ化しておきたいんですが、
参考意見お願いします。並列化なども。

昔、fjに流れた。
もしくは、octaveの中にwrapperが入っていたはず。

自作したら、めちゃめちゃクラスができた。MatrixとかSymMatrixとかBaseMatrix
とかEigenとか
%

言ってることがわからない。octaveって？wrapperって？

http://www.mediafoundry.gr.jp/~crg/katalog/gazou/miko.jpg

　　　　　　　　おにぎりワッショイ！！
　　　　　＼＼　 おにぎりワッショイ！！　／／
　+　　　+　＼＼　おにぎりワッショイ！！／+
　　　　　　　　　　　　　　　　　　　　　　　　　　　　+
.　　　+　　 ／■＼　　／■＼　 ／■＼　　+
　　　　　　（　~∀~∩（~∀~∩）（　~ー~）
　+　　（(　（つ　　　ノ（つ　　丿（つ　　つ　)）　　+
　　　　　　 ヽ　 （　ノ　（　ヽノ　　）　）　）
　　　　　　　（＿）し'　し（＿）　（＿）＿）

　　　　　　　　おにぎりワッショイ！！
　　　　　＼＼　 おにぎりワッショイ！！　／／
　+　　　+　＼＼　おにぎりワッショイ！！／+
　　　　　　　　　　　　　　　　　　　　　　　　　　　　+
.　　　+　　 ／■＼　
　　　　　　（　~∀~∩
　+　　（(　（つ　　　ノ　+
　　　　　　 ヽ　 （　ノ
　　　　　　　（＿）し'

　　　　　おにぎりワッショイ！！
　　＼＼　 おにぎりワッショイ！！　／／
+　+　＼＼　おにぎりワッショイ！！／+
　　　　+
.　　　+　　 ／■＼　　　　+
　　　　　　（　~∀~∩
　+　　（(　（つ　　　ノ　　+
　　　　　　 ヽ　 （　ノ
　　　　　　　（＿）し'

　　　　　祝2周年ワッショイ！！
　　＼＼　 祝2周年ワッショイ！！　／／
+　+　＼＼　祝2周年ワッショイ！！／+
　　　　+
.　　　+　　 ／■＼　　　　+
　　　　　　（　~∀~∩
　+　　（(　（つ　　　ノ　　+
　　　　　　 ヽ　 （　ノ
　　　　　　　（＿）し'

メッシュ切りが簡単で計算も比較的速いなんて都合の良いお勧めのソフトはある？

ウリの攻撃ツールでチョパーリどもが
大慌てニダﾞ!いい気味ニダ!
＿＿　＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
　　　.∨
　　　　　　　ｶﾀｶﾀ＿＿＿_
　　　∧＿∧　....|｜＼　ＬＧ＼
　　 <丶｀∀´>　 |｜ 　|￣￣ |
　┌（　　つ／￣|｜／ ￣￣/
　 |└ ヽ　|二二二」二二二二二二二二」
　　￣］|_＿７_７ .|　||　　　　　　　 |　||
　　／￣＼　 ／ . || 　 　　　　／　 ||
　 ◎　　◎..［＿___||　　　 　 .［＿＿||

☆保証☆
1ヶ月で1通のメールが登録メールアドレスに届きます。
このメールを確認（メールに記された広告主のURLをクリックする事が必要です）するだけで小切手1300円が家に直接配達されます。
通帳等の口座番号は必要有りません。
日本のサイトでもメール1通を確認すれば１〜５円貰えますが、
Net先進国アメリカでは10ドル貰えます。（国としてのスケールの違いですね。）
このサイトはアメリカのサイトで、2002年現在、日本ではほとんど加入者がいない状態です。
（もう世界的にはアメリカ、韓国等で200万人以上の会員が加入済みで実績が有ります）
http://www.MintMail.com/?m=2177421
このURLをクリックして入力すれば登録完了です。
（上記でからでないと獲得不可）
☆経験談☆
私もこの広報文を見てすぐ加入し、あちらこちらの掲示板10箇所ぐらいにこのメモを書いておきました。
先月、家にUSAから小切手が送金され,
これを地元の大和銀行で円に変えたみたら驚くことに60200円になりました。
おっ、と思いサイトに入って私の会員数を確認してみると、
私の直接会員が17人で、その会員が集めた会員が48人だったんです。
こんなに早く会員が増えるなんて、と驚きました。（ちゃんと確認できるところが安心します。）
でもこんなに早く私にもできたのは、まだ日本では会員が少ないからだと思います。
正直早めに会員に登録したほうがベストです。
早ければ早いほど上位会員になって獲得金額が多くなります。
☆収入の内容☆
ミントメールでは月10＄分のメールが届きます。
読めば（実は英語ですが確認だけでいいです）10＄が積立できます。
そして集めた会員がメールを確認してくれればその会員を抽選したあなたには6＄が積立できます。
そしてその会員がまた会員を集めたとしたらあなたには上位の上位の会員として4＄をもらいます。
コツコツが基本ですね!
もしあなたが会員を20人を集めてまたその会員が20人ずつ集めたとしたら
10＄⇒１３００円（あなたの積立金）
120＄⇒１５６００円（1段階下位会員よりの積立金：20人×6＄）
1600＄⇒２０８０００円（2段階下位会員よりの積立金：20人× 20人×4＄）

age

>>42
いいから速く俺に金送れよ。

http://www.windperfect.co.jp/

>>34がスレを潰したようです。

パソコンなどの箱入り電子機器をうまく冷却したいけど、シミュレーションを使うと事前にどの程度のことが分かるの？

>>47
筐体内の温度分布とか出るよ。
因みに47さんはリーマンですか？それとも単なる自作マニア？

>>48
リーマンだよ。
正確な温度を素早く得られるの？

>>49
素早くとはどの程度のオーダーですか？
因みに精度を上げれば、もちろん計算時間がかかるよ。
まぁ、これ以上は２ｃｈで聞くより、CFDソフト売ってるところに
ベンチマークを頼んだほうが手っ取り早いと思います。

>>47
CPUやDISK等のパーツの温度なら±5℃くらいの誤差まで分かるかな。
PCくらいの箱だったら計算時間は数時間。箱をモデリングするのに
数日かかるかな。

>>50
３０分以内。無茶ですな、やっぱり。
思った経路に風が流れるか通るか分かれば、とりあえずＯＫと思ってました。
>>51
±５℃なら思っていたより精度が良く、数時間なら待てるかも。
でも、モデリングに数日とは少し萎えました。

とにかく助言ありがとうございました。ベンチマークも活用してみます。

おいらあふぉでちたネ。

（＾＾）

現在、流体計算をすべて自作プログラムで（メッシュも）計算して研究してます。
しかしなぜ汎用プログラムを使わないのかと言われ困ってます。
自作プログラムの利点として自分の好きなようにプログラム改正ができることと
物理的な理解が深まることがあると自分では思ってますが、汎用プログラムを
使うことの欠点があるとすればなにかあるでしょうか？
知っている人がいれば誰かお願いいたします。

>>55
ならば君のソフトで日構造講師は対応できるのかい？
不連続メッシュ、テトラメッシュは？
移動協会は？複写は？使える乱流モデルは？
悪いけど、アマチュアの核プログラムなんて、工学的には
つかえないものばかり。党内噴射、ターボ機器、
熱交換器、などなど説いたことあるのかい？
制度検証してみたかい？

逆に使えるもので自身があるなら、自分で会社つくりなさい。

私は学生なんで卒業するためにある流体問題を解いてるだけです。
一応移動境界値問題です。
ただ、卒業審査で
「なんで自分で作ってんの？なんで汎用ソフト使わないの？」
っていう質問がよくあるのでその対策で教えていただきたいのです。
簡単にでもいいのでどうか汎用ソフトの欠点および自作プログラムの利点
があれば（教授を納得させることができるぐらいの）教えて頂けないでしょうか？
宜しくお願いいたします。

>>56
汎用ソフト作ることが>>55の目的じゃないだろーよ )ｳﾞｫｹ

>>55 >>57
学生さんなら自分で作るってゆーのはいいよ。（全部は無理かも知れんけど）
中身の理解ができるしね。自分の研究に訳の分からんモノ使ってても不安でしょ。

まあ、CFDで何かの現象を調べて考察することが主目的ならソフト作ることに時間
とられてちゃ本末転倒かもしれんけどね。

移動境界のコード書けるんなら立派なもんぢゃ。「なんで汎用使わないの」とかいう
DQN教官には「「学生なんだから！」」と答えとき。

長レス、スマソ...

お金がないとか

まにあうんか？

いねえなあ

親切に教えていただける方がいらっしゃる気配がなさそうだったのでしばらく
板を見てませんでした．
ありがとうございます．参考になりましたっていうか審査の不安が少し消えました．
それでは審査がんばります．
ありがとうございました．

>>62
移動境界できるし、学生なんだから精度わるくても
計算できりゃいいという安直な考えならば審査うける□はありません。

まずそのコードの信頼性を証明しなさい（それがテーマなんだろうけど）。
要は汎用コードとの差は実績が決定的。

>>58みたいなヴぁか野郎を信じるとなんでも
いいから計算できた。よかたよかたになる。
それじゃまずいと思われ。学生なんだからで修士/博士が取れたら
だれも苦労しないよ。

確かに「要は汎用コードとの差は実績が決定的」ですが、企業が何億もつぎこんで
開発した汎用プログラムと学生の自作とを比べたら実績（精度）で勝つのは無理だと思います。
学生だから・・・という甘えがかなりありますが。
ただ「開発」ではなく「大学研究」という観点から自作を推し進めてきましたが、
やはり審査の質問はきついです・・・
まあ　がんばります。
>>58さんも>>63さんもありがとうございました。

お前らたのむから聞いてやんなさい

ーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーー
是非とも御協力ください。
下のＨＰをクリックして、現在２位の小玉に票を入れてください。
いずれ大仁田をぬいて１位にしたいのです。
できれば何回もおながいします。

http://kenji03.hp.infoseek.co.jp/cgi-bin/kage/votec.cgi

ーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーー

v

http://www.agemasukudasai.com/bloom/

■■出会い系サイト運営システムレンタル■■

儲かる出会い系ビジネス

初心者でも簡単運営

写メール、画像対応

http://kgy999.net/open/

>>64
とりあえず2次元の
簡単な例題でもやってみてください。
健闘祈ります。

>>64
がんがれー。作った過程も無駄じゃないからね。

>>69
でもそれじゃ、あくまで課題演習。
間違っても土建屋さんや電中研みたいに研究発表と称して、
発表することのないようにね。

■■無料レンタル掲示板■■

どんどんレンタルして下さい

ランキングありジャンルも豊富です


http://kgy999.net/bbs/

>>55
　欠点。
・高い。
・下手をすれば、ブラックボックス化。

　一からプログラミングをすることは、素晴しいことだと思います。正直、羨ま
しい。企業に勤めてから、業務でＣＦＤをやり始めると、正直そんな時間はあり
ません。とはいえ、>>63さんの言うように、「現実との」照らし合わせは重要で
すね。で、>>63の書き込みにも疑問はあります。「汎用コードの実績って、そん
なに信頼に値するものかな？」と。汎用コードの計算結果と現実の照らし合わせ
の発表会があるくらいですから・・・・。今、家だから、その発表会の資料が手
許にないけど、自動車表面からの乱流剥離点の解析結果なんか、結構ばらばらだ
ったと記憶している。

mathematicaで流体力学のシミュレーションできますか？

>>55
欠点
･融通が利かない。

2次元問題や、軸対象流れであるのに、
3次元でメッシュを切らなければいけなかったりする場合がある。
無駄なメモリの消費と、計算時間の増大なんかが出ちゃいます。
ただ、自作コードがうまく動いているかなどのチェックに当てられる時間に対して、
どっちが余計かかるかなぁというところではありますが。

>>75
>2次元問題や、軸対象流れであるのに、
>3次元でメッシュを切らなければいけなかったりする場合がある。
今時そんな制限のある汎用ソルバーなんてほとんどないだろ？

つーか
２次元問題って必要あるの？

>>75
　Ｃ○Ｘのことでつか？「３次元でメッシュを・・・」というわけで、>>76さん、
あるんですわ。まあ、ＣＦ○、高性能とは思いますが。

>>77
　計算資源節約のためには、できるならば２次元で行ないます。特に、非定常
計算は。

（＾＾）

>>78
ずばりそれ。

汎用コードは精度が高い、ってのが前提になってるようだが、
どっからそんな法螺が出てきたんだか…。

スキームやモデルの研究をしてる人間が聞いたら失笑するよ。
そもそも「汎用」コードなんて未だ有り得ないものだ。

多様化するＰＣビジネスの中で、確実に急成長している堅実ビジネスです。
私自身、様々なネットビジネスを試してきましたが、
今回初めて自信を持ってご紹介できます。
初期資金はたったの1000円！
参考までに他では、1万円〜２万円が相場です。
ばらつきはありますが、平均日商8000円位になっています。
現在、土日を除く毎日が給料日です。毎日数時間、頑張って掲示板に
書き込むだけで、口座に毎日入金されるようになります。

詳しい内容をメールにて返信いたします。件名「情報希望」の上、
[email protected]まで、お問い合わせ下さい。

名前ランのピットをクリック→件名に情報希望→送信でメールを送れます。
私と同じヒマな人は試しにメールしてみてください。

http://www.agemasukudasai.com/bloom/

↓↓　ちょっとエッチなDVDが今ならなんと!一枚900円！　↓↓
　　　　http://www.net-de-dvd.com/
↓↓　特売品は、なんと一枚600円!!　↓↓
　　　　http://www.net-de-dvd.com/
↓↓　数に限りあり!今だけ！早いもの勝ち！急げ　↓ ↓
　　　　http://www.net-de-dvd.com/

>>81
はは。

汎用コードってのは、ま、使い方次第では役に立つこともある、
という感じで使われているんじゃないだろうか。

ていうか、そもそも汎用だろうがなんだろうが流体解析のコード
で他人に「失笑」されないコードってあるの？　あるなら、みんな
でそれ使や終わっちまう話なんだけど。

少なくとも、多くの汎用ｺｰﾄﾞは多くの研究者のｵﾅﾆｰｺｰﾄﾞよりは遥かに実用的。

http://www.agemasukudasai.com/bloom/

出会えるサイトNO1
　　　女性に大人気女性の方も訪問してね！！
http://asamade.net/cgi-bin/pc_i_j_ez-index.cgi
　　　　　　　　代理店募集中
　　　　　　貴方も月100万円稼げる！
　　　　　　http://asamade.net/web/
　

多様化するＰＣビジネスの中で、確実に急成長している堅実ビジネスです。
私自身、様々なネットビジネスを試してきましたが、
今回初めて自信を持ってご紹介できます。
資金は後にも先にもたったの1000円のみです。
参考までに他では、1万円〜２万円が相場です。
ばらつきはありますが、平均日商8000円位になっています。
現在、土日を除く毎日が給料日です。掲示板に
書き込むだけで、口座に毎日入金されるようになります。
同時に２００以上の裏情報も入手できます。

興味のある方、ヒマな方は名前ランにあるWINをクリックして
件名に情報希望を入力し送信してください。

または、[email protected]に
件名を情報希望として送信してください。

送っていただいても説明以外のメールは送信しませんので
安心してください。（１通のみ送信します）

出会えるサイトNO1
　　　女性に大人気女性の方も訪問してね！！
http://asamade.net/cgi-bin/pc_i_j_ez-index.cgi
　　　　　　　　代理店募集中
　　　　　　貴方も月100万円稼げる！
　　　　　　http://asamade.net/web/

@

保守

http://www2.leverage.jp/start/

　　 ∧＿∧
　　（　　＾＾ ）＜ ぬるぽ（＾＾）

∧＿∧
　　（　　＾＾ ）＜ ぬるぽ（＾＾）

http://homepage.mac.com/ayaya16/

　　　 　∧＿∧
ﾋﾟｭ.ｰ　(　　＾＾ ） ＜これからも僕を応援して下さいね（＾＾）。
　　＝〔~∪￣￣〕
　　＝ ◎――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉

大学が、税金使って、イギリスや米国製の汎用コード何故導入すんのか？
非常に怒りを覚える。

大学とは、超スーパーソフトを生み出す人材を育成するとこだろ。
自ら作成するか、汎用コード導入するにしても国産に限定すべき。

国民の税金で、国立大学がヘボな欧米製ソフトを導入し、
そのマーケ活動に加担し、国産ソフト駆逐に寄与する。
その腐りきった理念が、さっぱ理解でけん。

http://life.fam.cx/

格子で空間をいくつかに分けて流体をシミュレートするときって
重力加速度を扱うときどんな感じでやったらいいんですか？
空間が格子で分けられてるわけだから
どの微小空間にどれだけの流体が含まれてるかってことは
わかるかもしれないけど、
もし流体が空間内を落下する時、流体が今空間のどこにあるのかって言う
正確な情報を扱うことができなくなると思うんですが・・・
俺流体に関してはまったく知識ない度素人なんでよくわからないんですが
俺にもわかるように誰か簡単に教えてください。

☆☆☆ネット広告掲載で副収入！☆☆☆
http://look2.info/kurikuri/
PCでも携帯サイトでもOKです！いますぐ無料で作成できます。
よくあるあの「うざい広告」を貼り付けておくだけで
知らず知らずに収入になります。
PCの方は画面表示にそってお進みください。広告の説明をいたします。
⇒かわいいこのあそこ
⇒一番上の広告をクリックして出会い広告へ進んでください！
出会い広告をそのままクリックすれば広告説明の画面に移動します。
↑こちらのサイトは一般者からのＨ画像投稿掲示板です。
あわせてお楽しみください。
http://look3.info/kurikuri/
☆☆☆出会い系サイトは「無料の時代」です！☆☆☆
当サイトは広告料金収入で運営されておりますから、登録・ご利用はすべて「無料」です。
アドレス等も非通知で、送信できますから安心してご利用いただけます。

みてね♪
http://s-rf9.free-city.net/page003.html

http://homepage.mac.com/hiroyuki43/jaz10.html

http://homepage.mac.com/ayaya16/

ここってなかなかいいと思いません？安いし
　
http://www.dvd-yuis.com/

航空宇宙技術研究所のUPACS参加しろ

◎大人になってから◎
http://endou.kir.jp/betu/linkvp2/linkvp.html

【クレジットカードのショッピング枠を現金化！】

クレジットカードで買い物をすると、商品代金の８０〜９０％のキャッシュバック！！
１０万円の買い物をすれば８〜９万円の現金がお手元に！！

■■朝日カード■■
http://www2.pekori.to/~ryu/asahi/

>100
漏れも流体に関しては度素人なのでキミがよくわからない理解で
よくわからない悩み方をしているのはよくわかる。

漏れ流素人的な闇〜
一様流じゃない場合、時間ステップごとに隣から流れてきた流体を
中に残っている流体と単純に平均しちゃったら「拡散」とか「熱流束」
とか「粘性」とかが勝手に発生してしまいますよね。どうやってそれを
ゼロにしてるんでしょうか？

>>109
言ってる意味が>>100と同様さっぱりわからんのは私も素人だからなのか･･･
どうせまともに答える人なんて1年に一人くらいしか来ないんだから
暇つぶしと思って私の相手をしてくださいYO

>時間ステップごとに隣から流れてきた流体を中に残っている流体と単純に平均しちゃったら

隣の計算点orセルから流入があるてことですか？
単純に平均しちゃうとは流体が混合するという意味ですか？

>「拡散」とか「熱流束」とか「粘性」とかが勝手に発生してしまいますよね。

勝手にってどういう意味でしょうか？できれば詳しく説明してほしいんですけど。

>どうやってそれをゼロにしてるんでしょうか？

ゼロにしなきゃいけないんですか？
拡散も粘性もゼロって完全流体か何かでしょうか？

なんかｱﾌｫなレスでｽﾏｿ･･･

いえ、動的問題になると途端に私の頭が重力多体問題とかのレベルまで
退化してしまうもので。

勝手にっていうのはどういう意味なんでしょうね（w　1セルの状態を
変数1組であらわすにはセルの中を一様になるまで混合しないと
いけないような気がしたわけですよ。しかもそうやってできる見かけの
粘性って言うのは流れの方向によって変わってくる（層流境界層が
きちんと格子に平行に流れたらそういう混合が起こらない、とか）
……やり方が根本的に間違ってるような気がします。

ちょっと話を変えて、完全流体の定常流れを流れ関数のポアソン方程式？
で解く場合、粗い数値解（たとえば差分法に限定して）は厳密解に比べて
どの方向におかしくなるのでしょうか。
1. 連続の式が近似的にしか満たされない。
2. 微妙に渦あり流れになる。
3. 差分の意味では厳密に正則なのでまったく問題ない。（そんなのあり？）

重力多体問題へ退化って、それは退化なのですか？（ぉ
流体計算の粒子法とか格子ボルツマン法とかはそっち方面に近いのですかね？
私はあんまり知りませんけど。

えーっと、話が変わったほうから（わら
流れ関数のポアソン方程式を残差がなくなるように繰り返し計算するわけですよね。
完全流体なので上の作業は連続の式の残差を０にするのと等価になります。要は
計算が終われば連続の式は満たされたことになります。ポテンシャル流の便利な
ところですね。
つーことで結局何が問題なのかということですが、多分格子依存性が出てくるという
ことではないでしょうか？流線を書いてみるとガタガタした感じになると。そんだけ。

などと書いてますがあんまり専門ではないので大嘘書いてるかもしれません。ただ
少なくとも楕円型方程式の数値計算の場合に主要な問題となるのは計算時間であって
数値計算の教科書にはそれ以外のことはあんまり書いてないような気がします。

んで元の話に戻りますが、その見かけの粘性って何のことだろうと考えていたの
ですが
１．実はレイノルズ応力のことを言っていた。
２．粗い格子によって生じる余分な粘性効果のことを言っていた。
のどっちかではないかなーと思ったのですがどうでしょうか？１０９さんが物理や数学
には詳しいけど流体力学独特の考え方に慣れてないと１のように疑問に思うかな〜？
あるいはもっと詳しい人なら２のことをいっているのかな〜？と。２だったら格子を
めちゃくちゃ細かく切ったらそういう問題は起こりませんよと言えるかと思いますが。

こんなレスでいいんでしょうか＞偉い人
つーかなげえよ＞俺

見かけの粘性って２．です。見かけじゃなくて人工的な項、でした。スマソ。
δu/δt - uδu/δxを素直に計算しろ、というのが結論になるのですね。

ところで100はどこに行った〜

　__∧＿∧_
　|（　　＾＾ ）|　＜寝るぽ（＾＾）
　|＼⌒⌒⌒＼
　＼ |⌒⌒⌒~|　　　　　　　　　山崎渉
　　 ~￣￣￣￣

>>100
いやあ、悩みますね。界面を問題にする計算なら、その辺はっきりさせるんだけど・・・。
数値拡散の問題でもあるのかな。

>>115
やっぱ難しいんですね・・・。
みんな質問の意味がよくわからないって言ってたから
そういうの関係なく計算できちゃうのかと思ってた。

http://homepage.mac.com/hiroyuki44/

　__∧＿∧_
　|（　　＾＾ ）|　＜寝るぽ（＾＾）
　|＼⌒⌒⌒＼
　＼ |⌒⌒⌒~|　　　　　　　　　山崎渉
　　 ~￣￣￣￣

山崎渉うざ

☆頑張ってまーす！！☆女の子が作ったサイトです☆
　　　　　　　☆見て見て！！
http://yahooo.s2.x-beat.com/linkvp2/linkvp2.html

すごいよ〜♪

http://angely.h.fc2.com/page001.html

http://homepage.mac.com/hiroyuki44/

http://homepage.mac.com/hiroyuki44/

あぼーん

　　 ∧＿∧
　　（　　＾＾ ）＜ ぬるぽ（＾＾）

　　　 　∧＿∧　 ∧＿∧
ﾋﾟｭ.ｰ　（　 ・３・） (　　＾＾ ） ＜これからも僕たちを応援して下さいね（＾＾）。
　　＝〔~∪￣￣￣∪￣￣〕
　　＝ ◎――――――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉&ぼるじょあ

ある回転している円盤の中に液体を入れる区画を作り
その回転している円盤の区画に、中心点からそれぞれ液体を入れた時の、遠心力だ
内部の力関係を分析したいと思うんだけど、そういうのってシミュレーションで出来るかね？？
専門化に頼んだ方がいいかな？
専門化に頼んだ場合費用はどれくらい必要でしょうか？？

＞専門化

専門家

　　　 (⌒V⌒)
　　　│ ＾ ＾ │＜これからも僕を応援して下さいね（＾＾）。
　　⊂|　　　　|つ
　　　（＿）（＿）　　　　　　　　　　　　　　　　　　　　　　山崎パン

あげ

>>127
一日１０万円として、バクッと４日仕事だな。だから、４０万円。

Mathematicaを使って流体のシミュレーションを行いたいのですが，
参考になるサイトは無いでしょうか？
ご教示願います．

(・с_・´)

あげ

>>131
そういう仕事やる方は個人で給料いくらぐらいもらえんですか？

http://www.fsis.iis.u-tokyo.ac.jp/result/software/index.html

(´・_・)?

アーミールーバ＝仙道会＝仙界の古老

って何？

OpenFOAMってのが GPLで出たみたいだけどどうなのかな？
これまでの FOAMとの違いとか，他の CFDソフトと比較とか．
>>136との比較も興味ある．

http://www.opencfd.co.uk/openfoam/index.html#openfoam

東北大学のスーパーコンピュータSX-7が、HPCチャレンジベンチマークにおいて世界最高速記録を達成!!

　NECはこのたび、東北大学情報シナジーセンター(根元義章センター長http://www.isc.tohoku.ac.jp/)に納入し、2003年の1月から稼動しておりますスーパーコンピュータSX-7http://www.cc.tohoku.ac.jp/service/super/super-2.htmlにおいて、HPC(高性能計算)分野での
新しいベンチマークテストであるHPCチャレンジベンチマーク(http://icl.cs.utk.edu/hpcc/)で、28項目中16項目で世界最高速を達成いたしました。

　HPCチャレンジベンチマークは、米国政府の援助のもと、テネシー大学のJ.ドンガラ(J.Dongarra)博士を中心に米国や日本のHPC関係者が参画して、作成されたものであります。
従来、高性能計算機の性能番付として主流であったTOP500ベンチマークで使用されてきたLinpack(連立方程式の計算)を補完し、高性能計算機の性能を多面的な観点から評価する目
的で開発されたベンチマークであります。同ベンチマークはHPCシステムで重要となる複数の性能テストから構成されており、コンピュータの演算性能だけではなく、メモリアクセ
ス性能、並列処理のためのライブラリであるMPI(Message Passing Interface)の通信性能など、HPCシステムにおける様々な機能の性能を総合的に評価することが出来ます。

　今回のベンチマークテストでは、NECのSXシリーズのようなベクトル型スーパーコンピュータだけではなく、汎用プロセッサを使用したスカラー型スーパーコンピュータの性
能も数多く登録されており、今後のHPCベンチマークテストの主流になるものと見込まれております。

　このたび、SX-7が世界最高速を達成した、HPCチャレンジベンチマークは、Linpackベンチマークも含んだ7カテゴリー28項目のベンチマークから構成されております。

プレスリリース全文はこちらです。
日本電気株式会社プレスリリース　2004年12月21日
http://www.nec.co.jp/press/ja/0412/2102.html

>>139
それちょっと遊んでみたいかも。

ν速から来たお
記念カキコ！

>>142
何でそんなとっから来たんだ？

何かニュースんなった？

　　　 　∧＿∧　 ∧＿∧
ﾋﾟｭ.ｰ　（　 ・３・） (　　＾＾ ） ＜これからも僕たちを応援して下さいね（＾＾）。
　　＝〔~∪￣￣￣∪￣￣〕
　　＝ ◎――――――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉&ぼるじょあ

リアルタイム流体シュミュレーション
ttp://download.developer.nvidia.com/developer/SDK/Individual_Samples/samples.html
のVertex Texture Fetch Water ってところからダウンロード
ソースもある

ただしGeforce FX6800とかSM3.0か使えるグラボだけな
いいかげんCPUなんか使ってないでGPUにうつろうや

絵がキレイならOKという人は少ないからね。

>>146
結構居るよ。

テスト

この前、味噌汁作ってるときに鍋の中の対流って火加減によって結構違って見えて面白かった。
で、なんかでシミュレーションできないかなあ、って思ったんだけど、
こんな感じの2次元モデルがサクサク作れるようなソフトでおすすめってある？

家でいじりたいので、高いソフトは無理・・・

それは三次元モデルじゃないと無理っぽいよ。
でも、解析して何が明らかになるのかさっぱりわからない。

　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　　 　∧＿∧　 ∧＿∧　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
ﾋﾟｭ.ｰ　（　 ・３・） (　　＾＾ ） ＜これからも僕たちを応援して下さいね（＾＾）。
　　＝〔~∪￣￣￣∪￣￣〕 　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　＝ ◎――――――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉&ぼるじょあ
　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　

レスサンクス＞150

そっかー。
なんかね、強火だと鍋の外側から沸きあがるような対流になってて、
弱火だともうちょっと内側から湧き上がるような対流みたいに見えたんだ。
これは加熱位置か加熱温度の関係で対流が変わるってことかな？って思った。

不思議になってGoogleで調べてたら、加熱位置で対流が逆になる話とか
ｈttp://panasonic.co.jp/ism/dougama/html/02.html
２次元のシミュレーションの話とか
http://www.nesi.co.jp/01service/service2_2.html
http://www.eri.u-tokyo.ac.jp/hiroki/convection/2dconv/2dsuihei.html
がちょいちょいあったので、こういう２次元シミュレーションでそれを再現できないかなあ、
って思ったんです。

OpenFOAMはどうよ?

　　　
　　
━―━―━―━―━―━―━―━―━[JR山崎駅（＾＾）]━―━―━―━―━―━―━―━―━―
　　　　　　　　　
　　　 　∧＿∧
ﾋﾟｭ.ｰ　(　　＾＾ ） ＜これからも僕を応援して下さいね（＾＾）。　　　　　　　　　　　　　　　　　　　　　　　　　
　　＝〔~∪￣￣〕　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　＝ ◎――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉　　　　　　　　　　　　　　　　　　　　　　　
　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　__∧＿∧_　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　|（　　＾＾ ）|　＜寝るぽ（＾＾）　　　　　　　　　　　　　　　　　　　　　　　　
　|＼⌒⌒⌒＼　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　＼ |⌒⌒⌒~|　　　　　　　　　山崎渉　　　　　　　　　　　　　
　　 ~￣￣￣￣　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　 ∧＿∧　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　（　　＾＾ ）＜ ぬるぽ（＾＾）　　　　　　
　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　　 (⌒V⌒)　　　　　　　　　　　　　　　　　　　　
　　　│ ＾ ＾ │＜これからも僕を応援して下さいね（＾＾）。　　　
　　⊂|　　　　|つ　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　　（＿）（＿）　　　　　　　　　　　　　　　　　　　　　　山崎パン　
　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　　 　∧＿∧　 ∧＿∧
ﾋﾟｭ.ｰ　（　 ・３・） (　　＾＾ ） ＜これからも僕たちを応援して下さいね（＾＾）。
　　＝〔~∪￣￣￣∪￣￣〕　　　　　　　　　　　　　　　　　　　　　　　　　　
　　＝ ◎――――――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉&ぼるじょあ

...これからも僕を応援して下さいね（＾＾）。　　　
　　
━―━―━―━―━―━―━―━―━[JR山崎駅（＾＾）]━―━―━―━―━―━―━―━―━―
　　　　　　　　　
　　　 　∧＿∧
ﾋﾟｭ.ｰ　(　　＾＾ ） ＜これからも僕を応援して下さいね（＾＾）。　　　　　　　　　　　　　　　　　　　　　　　　　
　　＝〔~∪￣￣〕　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　＝ ◎――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉　　　　　　　　　　　　　　　　　　　　　　　
　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　__∧＿∧_　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　|（　　＾＾ ）|　＜寝るぽ（＾＾）　　　　　　
　|＼⌒⌒⌒＼　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　＼ |⌒⌒⌒~|　　　　　　　　　山崎渉　　　　　　　　　　　　　
　　 ~￣￣￣￣　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　 ∧＿∧　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　（　　＾＾ ）＜ ぬるぽ（＾＾）　　　　　　
　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　　 (⌒V⌒)　　　　　　　　　　　　　　　　　　　　
　　　│ ＾ ＾ │＜これからも僕を応援して下さいね（＾＾）。　　　
　　⊂|　　　　|つ　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　　（＿）（＿）　　　　　　　　　　　　　　　　　　　　　　山崎パン　
　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　　 　∧＿∧　 ∧＿∧
ﾋﾟｭ.ｰ　（　 ・３・） (　　＾＾ ） ＜これからも僕たちを応援して下さいね（＾＾）。
　　＝〔~∪￣￣￣∪￣￣〕　　　　　　　　　　　　　　　　　　　　　　　　　　
　　＝ ◎――――――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉&ぼるじょあ

3次元でやる場合、すぐメモリオーバー？

2次元の行列だけで計算しておいて、次の部分を計算したら前の部分を
ファイルに書き込んだ後すぐ領域捨てて次の部分、という風にやらないと
メモリオーバーになっちゃうのかな？

>>156
自作のプログラムか？

メッシュ数による

で、Cで書いてみようとしたが、そのバッファ部分のスワップが出来ないぽ
値渡しで無く参照渡しにすべきだろうが、さて配列ポインタってどうしよう？

単なるa,bのスワップならポインタとtemp使うんだけど、配列の場合どうしたものか
a[i]の先頭アドレスがaとだけ知っている

>>159

素人なら、行列計算の配列をグローバル変数にしてしまえ。

配列の渡し方など悩まなくて済む。

それで結局3次元に挑戦は永久に出来ない、っと

ムーアの法則に期待するしかない！

とりあえず、1次元熱伝導ならmalloc()とfree()使ってどうにか動いた
2次元の奴を動かすのに2次元配列作ったり、3次元作るのに4次元作ってたら
ムダも多いし
double a[10]とdouble *s作って、*sのためにmalloc場所確保、そこに*s=a[i]の式書いて
printf("%6.2lf",:*(s+i))で出力、その後a[i]=*(s+i)で値戻して後にfree(s)で開放

while(count<10)で連続計算、その間printf("%p %p",a,s)などとやってメモリ位置が無駄になっていないの確認

流体で圧力も関わるとかの場合、これが構造体やクラス（C++の場合）になるのか？
次は2次元行ってみよう。ホットカーペット表面がくまなく暖まる時間が計算できるのか

境界条件を設定するのが大変そうだね。
＞ホットカーペット

>>163
皆さん、シミュレーションの勉強には何使ってます？
今年から始める予定なんですが、おすすめの本があれば教えてください。
一応今のところ数値計算の基本書（CとかForranで微分とか書いてある）だけ
多少勉強しました。

下手に数値計算法の本だけ買ってきて、しかも1995年刊とかだと
Cで書かれたコードが、()と{}の間に変数名を定義する奴だったりして
今の奴だと動かなかったり

「Visual C++による理工学問題の解法（日刊工業新聞社刊）」という本が最近出ている。
現在のwindows中心のパソコン環境で、何か格好良いしC++でやろうかな
と言って形だけやりたい人ならこれが良い目。
Windowsの登場でgraphics.hが消えてしまった現在、画像出力やグラフ出力はgnuplotに
頼るのがベストだが、上記の本だとこれだとシングルダイアログベースのアプリで絵を描いてくれる
しかし、VC++によるコンソールアプリケーションの作り方の手順や
MFCAppwizardEXEの作り方の手順が分からなかったり、
いざ出来てもメンバ関数追加の方法が分からなかったりする場合ちと辛い。

「VC++ あんしん入門」（内堀高士著、技術評論社刊）辺りを読めば、
VC++の基礎的な使い方がかなり纏まっている。
CやC++を使おうとして、いきなりVC++買って「どの項目選択したら
俺が作りたいもの出来るのかわかんねぇー！」とならないように。
しかしこの本は文法に関しての記述は少ない。
各大学のホームページを漁ると、案外詳しい説明がある事も

文法抑えようと、東京電機大から出ている「学生のためのC&C++」やっても良いが
この本だけだとかなり抜けている部分も多い。C++だと参照やconstの使い方が
かなり難しいが、これについては載っていない。
実はクラスの作り方、演算子定義とかのややこしい部分は、かなり込み入った本を
やらないと載っていない。高橋マナの「やさしいC」「やさしいC++」は結構分厚いが
これですらクラスメンバの初期化について載っていなかったりする

そんな面倒なもの要らない場合が大半だろうが、これが分からないと
「オブジェクト指向計算力学入門（培風館）」はこれが分からないと理解できない

コンソールアプリだけで頑張るなら、「エンジニアのためのCプログラミング
（日刊工業新聞社刊）」が結構良い。

東京電機大から出ている本だと、「学生のための基礎C++builder」
「学生のための応用C++builder」も分かりやすく、結果が図示されてくれて面白い
基礎〜〜からやらないと、C++builder扱うの初めてならこれも混乱してしまう。
C++builderはストリンググリッドというエクセルの枠みたいなもので
入出力が出来るので、VBみたいに素人にも扱いやすい。
VC++にこんなコンポーネントは無い。しかしこの本は微分方程式は扱っていない。
ただ、C++Builderはコンソールアプリケーションを作るのが不可能ではないが
多少手間を要する（出来た雛形の一部を消したりする）ので、初期の文法の
勉強はちょっと厄介。何より、C++builderでCやC++の文法の勉強をする事を
前提にしてある本は多分存在しない。

数値計算だけに拘るなら、やらん方が良い本の方が多い
「VisualC++6.0プログラミング入門（アスキー出版局）」は、文法やMFCの構造その他が
相当分かっていないと使えないし、数値計算に関係する内容は殆ど無い。
しかしこれがVisualC++買って来ると一緒についてくるので「まあ最初は基礎からやるか」
と思って始めると挫折する
「10日でおぼえるVisualC++6.0入門教室（SHOEISHA）」も、途中からはWin32APIの説明ばかりなので
これも挫折する。10日でおぼえるの文句に引かれて買わないように

文法覚えるためにコンソールアプリだけやるなら、BCC落してきてパソコンにインストールした方が
良いかもしれない。しかし最初にパス設定などやるので、パソコン自体の初心者が
いきなりやるとウィンドウズ自体ぶっ壊れる恐れも、完全に0ではない。
またDOS窓からコマンド書き込んだりして面倒なので、CPADなどのフリーやシェアウェアの
開発･編集ツールをダウンロードしよう（ここでまたパスが必要だったりする）

CやC++で流体や熱伝導や波動をシミュレーションする本は、
殆ど見かけたことが無いしあってもやたら難しい。
JavaやVB、VBAあたりの方が分かりやすいが、
これでも3次元の流れを扱うものは存在しない。
また2次元流であっても、ちょっとややこしくなると出力に2分とか
3分とか掛かり始める。自分の作ったものが悪いのか、
パソコンが止まってしまったのか分かりにくい場合もあるが焦らないように

VBAだと「エクセル・シミュレーション入門（山海堂）」
Javaだと「流体・熱流動の数値シミュレーション（森北出版）」などが
流体力学を扱っている。しかしJavaは出力まで滅法遅い。

VBで問題になるのは、VBにはCのようにコンソールアプリで出力させるものが
現在だとフリーソフトでぼちぼちとあるだけで、またそこからいきなり
構造化BASICベースのVBをやった場合の落差が大きい事。まあ文法が易しいので
殆ど問題にならないが、VBでも難しい本で勉強した場合文法さっぱりになる
これも東京電機大の「学生のためのVisual Basic」あたりが無難か。

別スレで、「味噌汁煮てて対流で出来る模様のパターンが面白かったけど、
シミュレーションとか簡単に出来ますか？」と書き込まれていたが、
素人に3次元流体シミュレーションが出来るようになるにはあと20年は掛かるか？
C#にちょっと期待している（殆ど触った事が無いが）

・・・・と、あれこれやった挙句、自分の丁度良いレベルはVBAだなと納得している罠orz

>>166-168
夜分遅くにレスありがとうございます。参考にさせていただきます。
私自身、プログラミング経験がほとんどなく、言語に対する好き嫌いはありませんので、
いくつか触ってみてしっくりくるのを使ってみようと思います。特に挫折しなければCで行きたいなと。
出来れば基本的にコンソールだけですませたいので、CはBCCを使ってます。
（VCはワークスペーストかいろいろあって面倒だったので）

流体数値計算の本は、あっても各種方程式のモデルと説明に大半のページが割かれていて、
初心者にはハードルが高いですよね。（当然知っていないといけない内容なんでしょうけど）
ということでひとまず、単純なところから手をつけていこうと思います。

OpenFOAM 1.1 age

>>168 delphiやフォートランは？
>>170　やっぱlinuxか

www.asahi-net.or.jp/~uc3k-ymd/Lesson/index01.html
このページはポインタ廻りが纏まっているか。

>172参考に、2次元熱伝導も成功
しかしC++のnew delete [] は上手く行ったが、mallocが上手く行かん
どうも取得メモリバイト数違う

Cのmallocも成功
double** s;
s=(double**)malloc(n*double*);と、外側で2次行列のメモリ領域確保の後
for文で
for(i=0;i<m;i++)
{
s[i]=double(*)malloc(k*double);
}
と内側でも確保。後にfree(s)は別のループ抜ける時1回きりで良いようだ（？）

しかしメモリ周りﾏﾝﾄﾞｸｾ。これが出来なきゃCやC++の意味が無い。
流体でやるなら少なくとも構造体の2〜3次配列だが、中に関数含むならクラスの出来る
C++か

>>174

流体計算など反復計算をやる場合の配列は、一つの連続したメモリ空間を使用した方が高速化には好ましい。
配列の次元ごとにメモリを確保すると、メモリ空間はトビトビに作成される可能性がある。


ちなみに、貴殿のやり方だとメモリ解放は
for(i=0;i<m;i++){
free(s[i]);
}
free(s);
となる。

>>175　ｔｈｘ！

OpenFOAMいい!
ところで、これは、階層型メッシュは使えるのかな？

船舶工学の本に、3次元流体シミュのフォートランコードが載っていた

VC#で流体シミュやる本があったら、結構行けるかも知れん
ただC#ってまだマイナー言語だし、マイナーなまま消えかねん様相だし

体力学のコンピユター解析によって証明されておりますが、
専門当方は、入力＜出力と言う、流体機械装置を創って有ります。
流家である皆様は、コンピユター解析のデータを信用なさいますか。
コンピュター解析のデータと実際のデータとは、一致するとお思いですか。
この様な事を書いたのは、入力＜出力と書いて有る為に、皆様は絶対に信用しないだろうと思ったからです。
当方では、実験データと、コンピュタ解析データがよく一致しております。
何方か、当方の流体機械装置を、シュミレーションして、入力＜出力が真実である事を証明して頂けませんか。
興味の有る方は、御返事下さい。
詳しいデータは，ＨＰに載せて有ります。不明の所は、御説明申し上げます。　

　　　　　ホームページアドレス http://www.sutv.zaq.ne.jp/ckbil303　　　　　

流体力学のコンピユター解析
当方は、入力＜出力と言う、流体機械装置を創って有ります。
流体力学のコンピユター解析によって証明されておりますが、
専門家である皆様は、コンピユター解析のデータを信用なさいますか。
コンピュター解析のデータと実際のデータとは、一致するとお思いですか。
この様な事を書いたのは、入力＜出力と書いて有る為に、皆様は絶対に信用しないだろうと思ったからです。
当方では、実験データと、コンピュター解析データがよく一致しております。
何方か、当方の流体機械装置を、シュミレーションして、入力＜出力が真実である事を証明して頂けませんか。
興味の有る方は、御返事下さい。
詳しいデータは，ＨＰに載せて有ります。不明の所は、御説明申し上げます。　

　　　　　ホームページアドレス http://www.sutv.zaq.ne.jp/ckbil303　　　　　

＞コンピユター
ワロス

http://yuuki502.hp.infoseek.co.jp/battleload.mp3

東北大の流体研ってどうなん？

>>184

差分法が大好きなところ。

突然ですが、初歩の質問です。

ナビエ・ストークスの方程式の、発見の経緯
（実験や検証など）の記述のある本（もしくはサイト)を
探しています。

方程式を見ればそのままじゃないかと思われ
るかもしれませんが、発見者がどういう思考で
辿りついたのかに興味があります。

よろしくお願いします。

位置と流速の関数だろから、それを引数に含む関数仮定しておいて
それを微小位置、時間動かしてから、テイラー級数展開して
引き算してまあいっかなーーー、？？？？？

これは導き方の説明であって、経緯の説明じゃないな

(´・_・)?

１階：移流
２階：拡散
３階：散逸
４階：振動

シミュレートといったら有限要素法くらいしか知らない俺が
流体と剛体を同時に扱えるシミュレータを作ろうとか言ってるんですが
何を勉強したらいいですか？
もし可能なら剛体の破壊とかもシミュレートできるようにしたいんですが。

プログラミングはそこそこできるつもり。

>何を勉強したらいいですか？

何を勉強したらいいかわからないようでは無理だと思う

>>190 粒子法（SPH,MPS）
>>191 連続体力学(参考書多数あり)

>>190
> 流体と剛体を同時に扱えるシミュレータを作ろうとか言ってるんですが

東京工業大の矢部先生が開発したCIP法はいかが？
（と言っても俺もよくは知らないんだけど）

矢部研究室ホームページのCIP法によるシミュレーションCG
http://www.mech.titech.ac.jp/~ryuutai/graphics.html （）

(株)アドバンストテクノロジー の
CIP法による流体・固体・気体解析ソフトウェア PARCIPHAL
http://www.ad-tech.co.jp/08_15.html

↑の入門セミナー＆講習会
http://www.ad-tech.co.jp/10_0_PARCIPHAL.html

書籍「CIP法―原子から宇宙までを解くマルチスケール解法」
http://www.amazon.co.jp/exec/obidos/ASIN/4627918313/249-7725382-5173960

CIP法から派生したIDO法というのもある（東京工業大学の青木先生が開発）
http://www.nr.titech.ac.jp/~taoki/Activity/ido.html



つか，誰か，CIP法とIDO法の違いを教えてくれませんか？
今ひとつ違いがはっきりと判らないのです‥‥ σ(- -；

>>191-193
よーし、パパ頑張っちゃうぞ。

球面での流体シミュレーションを緯度経度座標と
有限差分法でやってみようかと思っているのですが，
（１）極が特異点になる
（２）高緯度で東西方向の格子間隔がとても小さくなり，CFL条件が厳しくなる
の二つをどうすればいいか判らず困っています。
特に（１）ではどんなことが問題でどんなマズいことが
起こるのかもよく理解できていません。
どなたか，緯度経度座標系＋有限差分法についての
良い解説資料やサンプルコードをご存じないでしょうか？

※ 陰陽格子のようなキメラ格子を使えば両方とも解決なのですが，
　 それを使う前に緯度経度座標でプログラミングして，
　 球面での自分の差分法のコーディングが間違っていないか
　 確かめたいので，まずは緯度経度座標で，と考えています。

>>170
使えます．
ブロックごとにfineにしたりとかできます．
サポート無しだとつらいかも．

素人なりに3次元熱伝導方程式シミュ？のC版、どうにか作って動かせた。
3重ポインタ気味悪いと思いつつも使えた
double ***sとか
*(*(*(s+k)+j)+i)とか。シンタックスシュガー使うと却って訳分からなくなるので
丁寧に書いた。

しかしまともに動くようになるまで、計算のおかしい部分が多くて
なかなか直せなかった。
どこからmallocして、どのタイミングでどっからfreeするんだ？とか思いつつ
for文で括った中でなんとか作って行った
ちゃんと中に伝わるようになり、使うアドレスも一定である事も表示できた。

しかし5×5×5で、周辺1、中0の状態から計算するだけで
1時間以上もかけて直すようじゃ、まだポインタの本質分かってないんだろうな
ポインタ関連の本3冊も読んだのに。

次は構造体にして流体にして圧力項も含んでみようかな。

197>>
動的配列は、C++のベクタ(std::vector)を使用してみたら？
すごーくらくちん。
ポインタは諸刃のやいばが本質のようです。
（　ポインタ不要論　で検索してみて）

かくしてポインタあやふやなまま、STLの海に投げ出されたのであった

200×200×200の配列
メモリ埋まる？

構造体のポインタ配列開放についてちょっとあやふやな点がある。
とにかく実行後、出力も終わってからエラーが出る

　　　　　　　　　　　　, -'"´ ￣｀丶､_
　　　　　　　　 　 ,.∩　　　　　　　　 ｀ヽ
　　　　　　　　　〃∪'´￣｀二二人＼ 　ヽ
　　　　　　　　　| ツ´￣￣￣￣´　ヾ ヽ.　',
　　　　　　　　　|ﾊ　,ニ、 　　,. - ､　| | | ｌ |
　　　　　　　　　| ﾊ ｨハ 　 　 ,二ヽ. | | | | |　同じ板にコピペするとそのままだけど、
　　　　　　　　　| | |　じ'　　　|トJ〉　 /)} ｌ | 違う板にコピペすると鬼のような怖い顔
　　　　　　　　　|　ﾊ　 ､'＿,　 ￣,,　厶ｲ川|　に変わる摩訶不思議な佳子様コピペ。
　　　　　　　　　l　l /＼　　　 ..　イV＼川 |
　　　　　　　　 ,'　l l　,ｲ ｀ｌ￣´ ／　　 /ヽl l
　　　　　　　　 ｌ　| ｌ ハ 　`メ､　　　 〃　 ヽヽ､__ノ
　　　　　　　　 l　 ∨ └‐ｲ「ト--ｧ'´　　 　 ﾊヽ__ノ
　　　　　　　　　ヽ/　　}　 l」」　/　　　　　/ }`ー
　　　　　　　　　 〈_ｎ| 八　　　/ /　　 　 /ノ
　　　　　　　　　 〈二二人 ｃ /＼/ /　, ｲ
　　　　　　　 　　　/　 /厂 ／＼__>＜ ｛＿
　　　　　　　 　　 /　 /　／　 ／ハ ＼＼　ト--- ､
　　　　　　　 　 _/　 /／　 ／　　ハ　 ＼＼　　 　 ＼
　　 ,　-‐ ''"´　/　 ／ 　／￣＼　｀丶、　　 ＼ ＼　　ヽ
　 /　 '´￣￣ / ／ 　／ ＼ 　 　＼　　＼　 　 ヽ　＼　ハ
　 |　　　　, -,--'　 ／　　　 ＼　 　 ヽ　 　＼　　 ヽ　　∨
　 ヽ､　　//　　 ／　　　　　　　ヽ　　 ＿＿__..--- ､　　ハ
　　 　 ￣〈| |　/　　　´￣￣￣￣￣/ /　　　　　 　 ￣`ヽ｝
　　　　　　しｰﾍ　　　　　　　　　 　 {__{＿_.... 、　 　 　 /　|
　　　　　　　　　 ` ー───　''"´　　　 　 　 ｀ヽ､　　|　/
　　　　　　　　　　　　　　　　　　　　　　　　　　　　 ` ┴′

↑
全然変わらないじゃんか

>>195
球座標はやったことないけど，円筒座標でならいつもやってます．
ってことで以下は全て円筒での話．
球もこの発展形でいけると思う．

(1) とりあえず以下を見てみて．
K. Fukagata and N. Kasagi, JCP 181, 478-498 (2002)

(2)は，周方向成分に陰解法を使えばよし．TDMAが定番かな．

>>204
情報ありがとうございます！ えと，
・TDMA⇒3重対角行列解法(orアルゴリズム)
ですか？

『 諸君 私は流体が好きだ 』
http://www.geocities.co.jp/Technopolis/4404/i_love_fluid.html

元ネタは勿論，あの 「諸君 私は戦争が好きだ」

非等方性乱流って，どれくらいのレイノルズ数から発生するものなのでしょうか？

>>206 の中で 「非等方性乱流の限りを尽くし‥‥ Re = O（10＾9） の様な流体計算」 って
あるけど，この 10＾9ってオーダーのレイノルズ数って，非等方性乱流の起こりうる値として
妥当なのですか？

>>207-208
つ [ http://science3.2ch.net/test/read.cgi/kikai/1092193978/144-146 ]

>>205
はいそうです．
ちなみにちゃんとやるなら普通の TDMAでなく，
その cyclic版を使った方がいいかもです．

>>210
重ねてご指導ありがとうございました m(_ _)m"

今のお金持ちの人はこの国がどうなってもいいのだろう。自分達の財産さえ残せれば。
構造改革とは国民に「他人を妬む心を植え付けること」 そして「他人を陥れることを推奨すること」 です。
経団連の奥田会長は「お金持ちの税金を安くし、消費税を１８％にすればいい」などと発言している。
本来なら革命でも起こすべきですが、一般大衆は無知ですし、牙を抜かれてしまい
何の気力もない。豚に餌を与えるが如くか。政府認定御用学者はマスコミでヘチョなことばか言って
私どもアホな国民を誘導する始末。新生銀行の社外取締役にデヴィッドロックフェラー�V世の名前が．．この人こそ親玉ですかね。
少子高齢化はそれほど問題でないとの認識である。それより若者の能力の質が落ちてる方がよほど大問題なのだ。
昔の社長は自分の手取りを少なくして、社員の給料を上げてたそのため一般の人達は裕福になり、
中流階級が生まれました。上・下流の二極化が進んでいる、そして中流階級は崩壊し
金持ちは海外の高級ブランドを貧乏人は中国などの安い製品を買う。中流崩壊は自国の製品の客をなくすのだ。
アメリカがいい例だ。ゼネラルエレクトリックは金融業でなんとかもってる。かくして日本には中流階級がいなくなり、
少数のお金持ちはガツガツと強欲にお金を稼ぎ、大多数の国民は保険には入れず病気や怪我などした時の生活は相当きついはずである。
ニューオリンズ市を見よ。まともな教育も受けれず、避難することもできずに苦しむアメリカ市民を、あれが日本の１０年後の姿かなぁ。

バブルリング（水中で息吐くとできる空気の輪っか。専門用語では渦輪）のシミュレーションをしたいんだけどどういうソフト使ったらいいのかな？

数値拡散のないコード
渦法とか粒子法とか

移流後に反応か、反応後に移流かで
出力結果が異なるのか

平衡考えると後者が良いのだろうが
計算ﾏﾝﾄﾞｸｻ

息をするのもﾏﾝﾄﾞｸｻ

じゃ氏ね

流れ方向と一致したメッシュを作成すれば数値拡散はかなり減らせるはずなので
メッシュも同時に作成するようにすればいいんでない？

すればぁｗ

>>198のいうようにvectorを見てみたら
下手するとかえってややこしくなるか？

しかし構造体で配列作った場合、関数を通したら
struct str
{
double a;
double b;
};
で、bを関数に通したらaにせっかく値が入っていたのに
何故か0になっちゃってたりする。
怖くて外部関数使いにくいなあ

動的配列使っているけど、こりゃ中も
double* a;
とかにした方が良いのかな？

VIPから来ました（　＾ω＾）

記念明子

そっとしとけクズ

同サイズ動的配列2枚使ったら、とにかく計算してくれた
最初に配列間計算して動配1つ目に代入、
次いでメンバ同士計算で動配2つ目に代入、
元の配列に戻す。
返り値のみアドレス型、引数は実値でﾖｶﾀ

でも「問題のあるポインタ代入」と警告されたorz

外装関数の内部変数追加して、アドレス返りを纏めて行ったら
エラーが消えてくれた。ﾖｶﾀ

構造体配列をアドレス受けで引数にして
void返しで中処理だけの関数作ったら
問題なく動いた

　　　　　　　　 ﾊ,,ﾊ
　　　　　 　 ('(ﾟ∀ﾟ∩_　おいらをどこかのスレに送って！
　　　　　　／ヽ 　 〈／＼　お別れの時にはお土産を持たせてね！
　　　　 ／|￣￣￣|.＼／
　　　 　 　| .ﾓﾂ煮..|／
　　　　　　 ￣￣￣
現在の所持品：たばこ・ライター・コーヒー・ブラックブラック・枕・ケータイ電話
睡眠薬・聖教新聞 ・ダッチワイフ・外付けSCSI340MHDD・ネットランナー4月号
TYG02・小嶋進社長・ゲルトモ・ディプスファンタジア・抵抗1kΩ(黄紫橙金)・媚薬
バーボン・リップクリーム・セガサターン バールのような物・バイブ・50インチ液晶TV
ボンカレー（甘口）・ジローラモ・エースヘルム・ＭＡＲ全巻セット・切断済みの亀頭
東京タワーのミニチュアキーホルダー ・無償の愛・おっぱいプリン・ｲｵﾅ・電話代の請求書・ウナコーワ
NHKの請求書・新庄のグラブ ・山田オルタナティブ・クイッククエンチガム・小隕石・こち亀71巻・こち亀78巻
ベゲタミン錠A・ニートの息子･マウス・炭酸ボンベ

　　　　　　　 ﾊ,,ﾊ
　　　　　 　 ('(ﾟ∀ﾟ∩_　おいらをどこかのスレに送って！
　　　　　　／ヽ 　 〈／＼　お別れの時にはお土産を持たせてね！
　　　　 ／|￣￣￣|.＼／
　　　 　 　| .ﾓﾂ煮..|／
　　　　　　 ￣￣￣
現在の所持品：DJかわちン（ｒｙ･J('ｰ`)し（嫁）・(`Д´)（そうた）
たばこ・ライター・コーヒー・ブラックブラック・枕・ケータイ電話
睡眠薬・聖教新聞 ・ダッチワイフ・外付けSCSI340MHDD・ネットランナー4月号
TYG02・小嶋進社長・ゲルトモ・ディプスファンタジア・抵抗1kΩ(黄紫橙金)・媚薬
バーボン・リップクリーム・セガサターン・バールのような物・バイブ・50インチ液晶TV
ボンカレー（甘口）・ジローラモ・エースヘルム・ＭＡＲ全巻セット・ハムとレタス
ゾフィーソフビ人形・キングカズ・外連味なく弁当2・イワシの缶詰・4番佐伯
アルビレックス新潟カレー・レディケット・ ストロング金剛・ディープインパクト
モルガ ・ニラ茶・タカラトミーの株・ＰＳ２ ・塩辛 ・エビマヨ・蟹味噌・納豆
加藤鷹・シナモンロール・猿岩石・クエン酸・ＵＮＯ・わら人形・いくら丼・おでん
遠近両用眼鏡・のどごし生・バタピー・ツチノコ・りぼん５月号付録・青汁・リポＤ

>>228
http://science4.2ch.net/test/read.cgi/material/949004149/96

>>198のいうベクタコンテナを調べてみた

これちゃんとメモリ開放してくれるんだろうな？ｗ

C#でも作れた。VBやJAVAより速いな

ここらで質量保存の式（連続の式に圧力追加）
も入れた計算をやってみよう
その後温度も。

にょろ？

関数化するならC、C++よりC#が、「勉強し易い」ｗ
参照も意外と早いなあ

流体力学のシミュレーションソフト

流体力学のシミュレーションソフト

流体力学のシミュレーションソフト

流体力学のシミュレーションソフト

流体力学のシミュレーションソフト

流休カ学のシュミレーショソンフト

流休カ学のシュミレーショソンフ

流休カ学のシュミレーショソンフト

流休カ学のシュミレーショソンフト

C#で完成！

面白いなあ
壁面を速度0にしてρを大きくしたら、次の層が
より外と比べて逆に動いたりしている

300行弱で書けたが、グラフィック表示考えたらこの倍以上になるな
現在出来るのはとりあえずファイル入出力だけ。データを流体の式に従い
計算してファイル出力する
初期条件と境界条件に汎用性持たせるのとかもどうやれば良いかなあ？

自分はぶきっちょゆえ「風を読め！」とか言われてもちっとも分からない訳だが
これならいろいろ観察が出来て大変面白い。

>>260
公開きぼん

シミュレーションというより趣味レーションだなｗ
まだ2次元だけだし次は3次元でやってみるか。

とりあえずは下に回転子のある洗濯機みたいな奴で
（洗濯物がある場合のシミュは自分には到底無理だ
移動境界とかまださっぱりだしｗｗｗ）

参考文献：河村本を3冊
エクセルVBAシミュの奴をC#の配列に書き直したｗ

有限要素法も将来やるとしたら、行列の各種計算も練習するか

>>260
スキームは何？
Mac法とか？

河村本だから勿論それ
動的配列取り回して収束させた

　∧＿∧　月曜日？ボコボコにしてやんよ
　(　･ω･)＝つ≡つ 　　(日曜 )
　(っ　≡つ＝つ 　　　 ∩∩　)　月曜から守ってください
　/　　　) ﾊﾞﾊﾞﾊﾞﾊﾞ 　　 / /　)
　( /￣∪ 　　　　　　　 LLニ⊃

　 　　　　　　　　　・,' ,., '・,';，.　ｵﾗｵﾗ平日様のお通りだぞ
　（ 現実）＝つ≡つ);;)日(;;(⊂≡⊂＝（月曜 ）
　(っ　≡つ＝つ　　(っ　 ⊂)　　⊂＝⊂≡　⊂)
　/　　　) ﾊﾞﾊﾞﾊﾞﾊﾞ ｜　 　｜ ﾊﾞﾊﾞﾊﾞﾊﾞ (　　　＼
　( /￣∪　　　　 　 ∪￣ ∪　　　　　　 ∪￣＼ )

　∧＿∧
　(　･ω･)　　そのとき、僕は見てるしかなかった・・・。
　U　　U　　　僕はなんて無力なんだろう・・・。
　|　　　|
　U￣∪

来年はオープンフォームがもっとはやりますように!

今年はOpenFOAMの年になりますように!w

オープンソースCFDツールボックス OpenFOAM part2
http://science5.2ch.net/test/read.cgi/sim/1159261169/

何か出力変だと思ったら()の位置間違ってたorz

VOF法＠26周年：2007/01/30(火) 00:06:14 ID:rdyYl0PV
VOF法のサンプルコードを探しているのですが，書籍でもホームページでも
いいので，どなたか在処をご存じでないでしょうか？

※ できたらVOF法についての良い日本語文献についても，ご存じの方が
　 おられたら教えて頂けませんか。

エロスレで純愛デートしましょう18日目
http://sakura03.bbspink.com/test/read.cgi/sureh/1168875512
このﾚｽを見た人は4日後にあなたの大切な人がいなくなるでしょう････
それが嫌ならこのﾚｽを5つの板にｺﾋﾟﾍﾟしてください｡
信じるか信じないかはあなた次第です｡

OpenFOAM age

千の風の挙動を流体シミュレーション

OpenFOAM
FrontFlow/Red,Blue
Adventure

このやたらsageているのは何が目的なの？

>>442　グリッド数10の3乗で1000個ｗ

Open Open OpenFOAM!

これスクリプト？
なんでこのスレなんだろうねぇ

知らない

なんという良スレ・・・
スレタイを見ただけでワクワクしてしまった
このスレは間違いなく伸びる
　　　　　　　
　　 ／￣＼
　　|　＾o＾　|　　
　　 ＼＿／

OpenFOAM イイネ!

埋まっちゃうよｗ

あげ

このスレに気に入らない情報があるんだろうなｗ

歌手 aiko がモロ盗作･･･　↓まとめ
ttp://a.weep.jp/gossip/?aiko%2F%C5%F0%BA%EE%B5%BF%CF%C7%CC%E4%C2%EA%2F%B2%E8%C1%FC

国分が愛想つかした理由がよくわかりました