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

**名無しさん** · NG

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

**名無しさん** · NG

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

**迷ったなら卍へ** · NG

http://homepage1.nifty.com/sads/index.html

**名無しさん** · NG

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

**2です** · NG

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

**>2さん** · NG

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

**2です** · NG

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

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

**2です** · NG

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

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

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

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

6 · NG

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

>8 · NG

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

**2です** · NG

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

**文無し** · NG

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

**2です** · NG

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

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

だぁ · NG

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

と · NG

テーラー展開の公式は
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 までを消すことができる。

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

と金 · NG

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

格子の数が11で、
２階微分を
{T(i-1)-2*T(i)+T(i+1)}/(⊿x^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

み · NG

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

**>16** · NG

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

**名無しさん** · NG

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

S · NG

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

み · NG

解析解が在るわけではないですが、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.

と · NG

#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);
}

**名無しさん** · NG

つづきをどうぞ＞２

**名無しさん** · NG

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

**名無しさん** · NG

http://www.yti.co.jp/kiryu/index2.html

**話追ってないけど > 22** · NG

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

**名無しさん** · NG

スタック？

**>27** · NG

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

**名無しさん** · NG

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

**無垢猫** · NG

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

**名無しさん** · NG

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

**名無しさん** · NG

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

**名無しさん** · NG

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

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

**名無しさん** · NG

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

**名無しさん＠１周年** · NG

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

**名無しさん＠１周年** · NG

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

**名無しさん＠１周年** · NG

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

**名無しさん＠１周年** · NG

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

**このスレ＠2周年** · NG

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

**名無しさん＠１周年** · NG

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

**名無しさん＠１周年** · NG

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

**名無しさん＠１周年** · NG

☆保証☆
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＄）

**名無しさん＠１周年** · NG

age

**名無しさん＠１周年** · NG

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

**名無しさん＠１周年** · NG

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

**名無しさん＠１周年** · NG

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

**名無しさん＠１周年** · NG

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

**名無しさん＠１周年** · NG

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

**名無しさん＠１周年** · NG

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

48 · NG

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

**名無しさん＠１周年** · NG

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

**名無しさん＠１周年** · NG

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

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

**名無しさん＠１周年** · NG

おいらあふぉでちたネ。

**山崎渉** · NG

（＾＾）

**初心者です** · NG

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

**忠告や** · NG

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

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

**初心者です** · NG

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

**名無しさん＠３周年** · NG

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

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

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

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

長レス、スマソ...

**名無しさん＠３周年** · NG

お金がないとか

**名無しさん＠３周年** · NG

まにあうんか？

**名無しさん＠３周年** · NG

いねえなあ

**名無しさん＠３周年** · NG

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

NG

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

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

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

**名無しさん＠３周年** · NG

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

**名無しさん＠３周年** · NG

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

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

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

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

**bloom** · NG

**bloom** · NG

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

**505** · NG

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

儲かる出会い系ビジネス

初心者でも簡単運営

写メール、画像対応

http://kgy999.net/open/

63 · NG

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

58 · NG

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

**名無しさん＠３周年** · NG

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

**548** · NG

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

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

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

http://kgy999.net/bbs/

**ピカレスク＠元々化学工学** · NG

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

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

**名無しさん＠３周年** · NG

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

**名無しさん＠３周年** · NG

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

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

**名無しさん＠３周年** · NG

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

**名無しさん＠３周年** · NG

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

**ピカレスク** · NG

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

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

**山崎渉** · NG

（＾＾）

75 · NG

>>78
ずばりそれ。

**名無しさん＠３周年** · NG

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

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

**ピット** · NG

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

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

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

**bloom** · NG

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

**名無しさん＠３周年** · NG

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

**名無しさん＠３周年** · NG

>>81
はは。

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

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

**名無しさん＠３周年** · NG

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

**bloom** · NG

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

**出会いNO1** · NG

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

**WIN** · NG

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

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

または、hiroaki_s_faith@hotmail.comに
件名を情報希望として送信してください。

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

**出会いNO1** · NG

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

**名無しさん＠３周年** · NG

**名無しさん＠３周年** · NG

保守

**bloom** · NG

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

**山崎渉** · NG

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

**名無しさん＠３周年** · NG

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

**bloom** · NG

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

**山崎渉** · NG

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

**名無しさん＠３周年** · NG

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

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

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

**名無しさん＠３周年** · NG

http://life.fam.cx/

NG

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

**名無しさん＠３周年** · NG

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

**名無しさん＠３周年** · NG

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

_ · NG

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

**bloom** · NG

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

**名無しさん＠３周年** · NG

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

**名無しさん＠３周年** · NG

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

**名無しさん＠３周年** · NG

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

**名無しさん＠３周年** · NG

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

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

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

**名無しさん＠３周年** · NG

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

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

**名無しさん＠３周年** · NG

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

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

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

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

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

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

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

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

**109** · NG

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

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

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

**名無しさん＠３周年** · NG

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

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

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

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

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

**109** · NG

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

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

**山崎渉** · NG

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

**名無しさん＠３周年** · NG

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

**100** · NG

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

_ · NG

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

**山崎渉** · NG

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

**名無しさん＠３周年** · NG

山崎渉うざ

**名無しさん＠３周年** · NG

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

**名無しさん＠３周年** · NG

すごいよ～♪

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

_ · NG

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

_ · NG

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

**あぼーん** · NG

あぼーん

**山崎渉** · NG

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

**ぼるじょあ** ◆yBEncckFOU · NG

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

**初心者a** · NG

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

**teisei** · NG

＞専門化

専門家

**山崎渉** · NG

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

**名無しさん＠３周年** · NG

あげ

**名無しさん＠３周年** · NG

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

**名無しさん＠３周年** · NG

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

**名無しさん＠３周年** · NG

(・с_・´)

**名無しさん＠３周年** · NG

あげ

**名無しさん＠３周年** · NG

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

**名無しさん＠３周年** · NG

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

美香 ◆RQ0Spv3q06 · NG

(´・_・)?

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

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

って何？

**名無しさん＠３周年** · NG

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

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

**名無しさん＠３周年** · NG

東北大学のスーパーコンピュータ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

**名無しさん＠３周年** · NG

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

**名無しさん＠３周年** · NG

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

**名無しさん＠３周年** · NG

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

何かニュースんなった？

**ぼるじょあ** ◆yBEncckFOU · NG

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

**名無しさん＠３周年** · 05/02/02 03:44:18

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

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

**名無しさん＠３周年** · 05/02/02 20:52:11

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

**名無しさん＠３周年** · 05/02/05 01:42:14

>>146
結構居るよ。

**テスト◇流星群** · 05/02/05 01:55:12

テスト

**名無しさん＠３周年** · 05/02/18 10:52:52

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

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

**名無しさん＠３周年** · 05/02/19 18:48:34

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

**ぼるじょあ** ◆yBEncckFOU · 05/02/20 08:45:24

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

**148** · 05/02/20 12:51:49

レスサンクス＞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
がちょいちょいあったので、こういう２次元シミュレーションでそれを再現できないかなあ、
って思ったんです。

**名無しさん＠３周年** · 05/02/21 21:06:26

OpenFOAMはどうよ?

**山.崎渉** · 05/02/22 18:58:34

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

**山.崎渉** · 05/02/22 19:12:04

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

**名無しさん＠３周年** · 05/03/12 18:23:21

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

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

**名無しさん＠３周年** · 05/03/12 20:17:00

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

**名無しさん＠３周年** · 05/03/12 23:03:31

メッシュ数による

**名無しさん＠３周年** · 05/03/13 03:00:55

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

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

**名無しさん＠３周年** · 05/03/13 10:04:04

>>159

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

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

**名無しさん＠３周年** · 05/03/13 13:07:44

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

**名無しさん＠３周年** · 05/03/13 19:49:54

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

**名無しさん＠３周年** · 05/03/15 01:20:02

とりあえず、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次元行ってみよう。ホットカーペット表面がくまなく暖まる時間が計算できるのか

**名無しさん＠３周年** · 05/03/15 20:07:00

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

**名無しさん＠３周年** · 05/03/15 22:42:45

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

**名無しさん＠３周年** · 05/03/16 03:38:59

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

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

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

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

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

**名無しさん＠３周年** · 05/03/16 04:07:57

コンソールアプリだけで頑張るなら、「エンジニアのための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などのフリーやシェアウェアの
開発･編集ツールをダウンロードしよう（ここでまたパスが必要だったりする）

**名無しさん＠３周年** · 05/03/16 04:31:49

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

**165** · 05/03/16 05:20:18

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

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

**名無しさん＠３周年** · 05/03/16 10:48:52

OpenFOAM 1.1 age

**名無しさん＠３周年** · 05/03/16 13:49:31

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

**名無しさん＠３周年** · 05/03/16 23:07:07

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

**名無しさん＠３周年** · 05/03/19 03:59:45

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

**名無しさん＠３周年** · 2005/03/22(火) 03:32:21

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++か

**名無しさん＠３周年** · 2005/03/22(火) 08:54:57

>>174

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

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

**名無しさん＠３周年** · 2005/03/22(火) 15:41:58

>>175　ｔｈｘ！

**170** · 2005/03/30(水) 01:04:29

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

**名無しさん＠３周年** · 2005/04/11(月) 19:47:48

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

**名無しさん＠３周年** · 2005/04/23(土) 03:36:34

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

**frontier7toiten** · 2005/04/25(月) 05:43:39

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

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

**frontier7toiten** · 2005/04/25(月) 06:06:40

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

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

**名無しさん＠３周年** · 2005/04/28(木) 20:24:49

＞コンピユター
ワロス

**名無しさん＠３周年** · 2005/04/28(木) 22:38:37

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

**名無しさん＠３周年** · 2005/04/29(金) 03:59:37

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

**名無しさん＠３周年** · 2005/04/29(金) 15:02:18

>>184

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

**名無しさん＠３周年** · 2005/05/02(月) 15:00:50

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

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

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

よろしくお願いします。

**名無しさん＠３周年** · 2005/05/03(火) 17:41:22

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

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

美香 ◆RQ0Spv3q06 · 2005/05/06(金) 08:08:12

(´・_・)?

**名無しさん＠３周年** · 2005/05/06(金) 20:56:19

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

**名無しさん＠３周年** · 2005/05/31(火) 20:46:34

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

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

**名無しさん＠３周年** · 2005/05/31(火) 22:56:31

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

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

**名無しさん＠３周年** · 2005/06/11(土) 02:39:27

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

**名無しさん＠３周年** · 2005/06/12(日) 00:31:02

>>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法の違いを教えてくれませんか？
今ひとつ違いがはっきりと判らないのです‥‥ σ(- -；

**190** · 2005/06/14(火) 15:49:23

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

**名無しさん＠３周年** · 2005/06/15(水) 20:05:03

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

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

**名無しさん＠３周年** · 2005/06/17(金) 21:34:48

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

**名無しさん＠３周年** · 2005/06/18(土) 01:35:40

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

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

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

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

**名無しさん＠３周年** · 2005/06/20(月) 02:14:48

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

**名無しさん＠３周年** · 2005/06/20(月) 04:19:12

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

**名無しさん＠３周年** · 2005/06/20(月) 20:17:15

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