圧縮・非圧縮性流体の統一解法

1 名前：み [2000/04/01(土) 18:42] 圧縮性流と非圧縮性流を統一的に扱うスキームは構成可能か？ 2 名前：名無しさん [2000/04/02(日) 12:13] もうすでに作ってしまっている人を私は知っています。 あと1年ぐらいで世の中に出回るでしょう！！！ 3 名前：名無しさん [2000/04/03(月) 22:18] どんな？まだ内緒なのかな。 4 名前：名無しさん [2000/04/06(木) 10:01] 圧縮性流と非圧縮性流はどんなスキーム？ 5 名前：み [2000/04/06(木) 21:04] 例えばTVD系の保存型スキーム＞圧縮 例えばMAC系のProjectionスキーム＞非圧縮 6 名前：名無しさん [2000/04/07(金) 11:08] 用語が全然わからない… 7 名前：名無し [2000/04/30(日) 16:26] もう売ってるって。 8 名前：nanashi-san > 7 [2000/04/30(日) 17:57] もしかしてCIPのことか？ あれは信用されてないぞ（わら 9 名前：＞８ [2000/05/05(金) 19:44] そんなのしらん。 宇宙業界では、何でもできる流体解析ソフトが３０年の改良を経て使われている。 一般人は、知らないだけ。 10 名前：奈々資産 [2000/05/05(金) 20:56] 宇宙業界?? 30年前ってことはICEかな。 11 名前：名無しさん [2000/05/13(土) 08:46] >8 CIP法(東工大の教授のアレですよね？)の場合は、まだ適用範囲が狭いからそう思われているのでは？ 個人的にはあの解法はロジック的に単純で好きです。 12 名前：名無しさん [2000/05/13(土) 12:25] いや、数値計算や流体の重要な部分を落としているから単純になってい るだけ。実際、問題が多い。信用されないのにはそういう訳がある。 ちなみに、某氏の計算には嘘がある。 13 名前：12＠補 [2000/05/13(土) 12:26] ＞ CIP 14 名前：通りすがり [2000/05/14(日) 01:28] そーなんですか。どんなとこ省いてるんですか？ 15 名前：名無しさん [2000/05/15(月) 17:34] CIPって何の略なん？ 16 名前：通りすがり２＞１４ [2000/05/15(月) 18:58] 　例えば古典的な多項式を用いて補間を行っている事。 良く知られているように、相界面やcontact discontinuityなどの領域 では物理量の値や微係数が不連続になる。すなわち、その様な領域では テーラー展開は収束し難いか、あるいは収束しない。よって、テーラー 展開から得られる古典的な多項式は、流体の空間プロファイルの補間に は相応しくない。これは昔から分かっている事で、NeumannやLax、比 較的最近ではHartenやOsherら、多くの研究者が「人口粘性」を導入す るというアプローチでこの問題に挑んできた。CIPではこの問題につい てはほとんど考慮されておらず、非常に古典的な地平にまで逆戻りして しまっている。 17 名前：通りすがり３＞１４ [2000/05/15(月) 19:10] 　例えば音波を陰的に扱ったというだけで統一解法だとしてしまっている事。 CIPでは陰的に離散化された拡散方程式のような恰好の方程式で音波を扱 うが、良く知られているように、後退オイラーによる拡散方程式の数値 解は拡散数が大きい領域で非常に（数値的に）拡散的になってしまう。 つまり、音速のCFLが（約）1以上の領域では音波は高精度には記述され ないことになる。このため、例えば音速が流速の数倍〜数十倍程度の条 件の時に音速のCFLを無視して流速のCFLのみでタイムステップを決定し てしまうと、しっかり音速のCFLまでを考慮してタイムステップを決定 した場合とは異なる解、つまり間違った解が供給されてしまう。 18 名前：通りすがり４ [2000/05/15(月) 19:16] これはCIPだからというわけではないが、CIPを用いた計算例には何故か、 渦が発生する問題であるにもかかわらず実粘性項を加えずに計算を行っ ているものがある。実粘性無し、つまりレイノルズ数無限大の渦流は有 限個の格子では解像する事が不可能。したがって数値粘性からくる、物 理的には意味の無い架空のレイノルズ数を用いた結果を得る事になる。 19 名前：例の名無しさん [2000/05/16(火) 11:00] ちょっと調べた。なるほど、ＣＩＰってHermiteを使ってるんやな？ そら確かに解に不連続があればどないもならんスキームやろうけど、 滑らかであればええんちゃうの？ まぁ、ただ多次元への拡張が 難しそうやけど（大概のスキームに言えることやけどね）。 （難しいってのは、技術的ではなく物理的に） 20 名前：通りすがり [2000/05/17(水) 13:49] 14です。丁寧な解答どうも＞16,17,18 数値粘性に無頓着なのは困ったもんですね。 ところで、もともとCFLは最も時間スケールが 短い(特性勾配がキツイ）現象を基準にするのでは？＞17 21 名前：>20 [2000/05/17(水) 16:51] しかしそれだと、流速が音速よりずっと小さく、かつ流速の時間スケー ルで起こる現象を解析したい、なんていう場合に凄く非効率的になるよ ね？＞ 最も時間スケールが短い現象を基準 そこで、音波を陰解法で扱って効率を上げてしまおうってわけ。 22 名前：>19 [2000/05/17(水) 16:52] 特性法とか、そっち方面の話かな。＞ 技術的ではなく物理的に 23 名前：通りすがり [2000/05/17(水) 18:07] ＞20 数十倍も違えば、たしかに。 最近、コンピュータパワーに頼り過ぎてる、おれを実感。 ウソよりいいかと。 らじゃ、です。 24 名前：すまん [2000/05/17(水) 18:17] ありゃ、23は＞21。 25 名前：もっと [2000/05/20(土) 23:07] 統一解に近づいている解法は無いのかな？ 26 名前：口 [2000/05/23(火) 14:09] R.Klein, "Semi-Implicit Extension of a Godunov-Type〜", J.Comput.Phys., 121, (1995) 213-237 なんて良いかも。CIPの論文よりずっと色々な議論がされていて面白い。 27 名前：ＨＨＨ [2000/06/05(月) 00:56] ＣＩＰも結局，圧縮性の問題を解くときには，数値粘性を入れています． 28 名前：なっち [2000/06/05(月) 15:12] でもそれはショックフロントにしか効かないね。＞数値粘性 29 名前：なるほど [2000/06/16(金) 20:12] 音波は後退オイラーで解いているのか。これでは数値粘性がキツイねえ。 30 名前：>15 mailto:今更 [2000/07/04(火) 12:19] Cubic Interpolated Propagation 31 名前：助教授 mailto:sage [2000/07/14(金) 10:14] TVDと異なり，CIPはリーマン問題をクリアしていません． 32 名前：ke [2000/07/31(月) 21:53] Ｔ大の？＞助教授 33 名前：> [2000/09/20(水) 14:17] 最近は色々と研究されてますね。 「〜Low Mach Number Flow〜」なるタイトルの論文は皆そうでしょう。 34 名前：> [2000/09/20(水) 14:19] あ、わたしもCIPはアラが多すぎると思います。 話が感覚的過ぎるんですよね。 35 名前：名無しさんZ80 mailto:sage [2000/09/22(金) 09:20] SPHとCIPの併用じゃ駄目？ ってゆーか、いまそんなの作ってる。 カオス挙動はΔtをやたら小さくしてなんとか。 やっぱり擬似粘性はいってるけど。 36 名前：hage mailto:hage [2000/09/22(金) 14:33] どの部分をSPHで、どの部分をCIP？ 37 名前：＞３４ [2000/10/04(水) 06:43] ＣＩＰは，直感的で分かりやすいから使いやすく， 発展しているのではないでしょうか．．． 38 名前：＞３７ [2000/10/04(水) 13:36] なるほどね。しかし実際は直感で分りやすいだけでは駄目なはず。 例えば非常に単純で分かり易い手法を提案したとしても、保存則 を満足するとは限らないといった手法は信頼できません。保存則を 離散的に満たさないスキームは衝撃波の位置を正確に捕らえられない と言う事実は基本中の基本です。スキームは単純（simple）であるべ きです、しかし単純であれば何でもいいと言う訳ではない。 ＣＩＰはCONSERVATIVEでない。 これが大きな問題ではないでしょうか？ 39 名前：名無しさん＠そうだ選挙にいこう [2000/10/04(水) 16:53] >>37 人間の直感だけで済んでれば苦労は要らないんだけどね。 40 名前：名無しさん＠そうだ選挙にいこう [2000/10/04(水) 16:59] >>38 ＞ＣＩＰはCONSERVATIVEでない。これが大きな問題ではないでしょうか？ 他にも問題は在るけど、１番大きいのはコレかもね。 一応、ショックフロントでの保存性の崩れを修正する為の人口粘性が 使われてはいるけど、その効かせ具合を調節する方法がまったく示さ れていない。弱い衝撃波の場合には人口粘性係数は0.5程度でも大丈夫 なようだけど、強いショックになると、なんと3〜5なんていう大きい 値が必要となることもある。こんなに広いレンジのある係数を決定す るための手段が示されていないんだから、恐ろしくて衝撃波問題には 使えない。 41 名前：名無しさん＠そうだ選挙にいこう [2000/10/05(木) 03:54] ここ、何気に専門家が居る気がする（わ 42 名前：ななし [2000/10/05(木) 15:36] 時に直感で理解できなくても良い結果がでる手法があります。 一番いい例はＡＵＳＭ Ｓｃｈｅｍｅ。はっきり言ってでっち上げ というか, cook up されたものです。しかし単純で結構うまくいく という経験的事実から広く使われているようです。(有限体積法で 使われるので保存則は満足してます。) 43 名前：名無しさん＠そうだ選挙にいこう [2000/10/05(木) 22:12] なんの略？＞AUSM 44 名前：ななし [2000/10/06(金) 01:57] AUSM = Advection Upstream Splitting Method. cf. Liou, M.S., and Steffen, C.J.: "A New Flux Splitting Scheme", Journal of Computational Physics, 107, pp23-39. 45 名前：＞３８ [2000/10/09(月) 21:22] ＞保存則を満足するとは限らないといった手法は信頼できません。 いや〜，それは，あなたが実際にＣＩＰで計算したことがないから分からない だけでしょう．．． さらに，最近のＣＩＰの論文をきちんと読んでいないからそう言った考えに なってしまったのでしょう．．． 確かに，当初の論文は，非保存系でのフォーミュレーションから スキームを構築されて発展してきました．（そのおかげで，圧縮，非圧縮 の統一解法への道も開かれました） その後，いろいろな人がＣＩＰを使うようになり，別のグループが ＣＩＰの保存性について議論している論文も２，３年前から出ています． 非常に保存性の良い，スキームであることが分かりました． ＃　ＣＩＰは，保存性が悪い！と言う論文が出ているのでしょうか？？ ＃＃　「直感的で分かりやすく，良いスキームだから」 ＃＃　私があえて言わなくとも，これだけ大勢の人が使っていれば， ＃＃　言う必要も無いでしょう．．． 46 名前：名無しさん＠そうだ選挙にいこう [2000/10/10(火) 10:54] ＞別のグループが どのグループ？ 47 名前：名無しさん＠そうだ選挙にいこう [2000/10/10(火) 11:03] ＞＃　ＣＩＰは，保存性が悪い！と言う論文が出ているのでしょうか？？ 保存型スキームより保存性が悪いのは自明だよね。 ＞＃＃　私があえて言わなくとも，これだけ大勢の人が使っていれば， 衝撃波問題への利用からは撤退したグループが幾つか在るね。 48 名前：名無しさん＠１周年 [2000/10/10(火) 12:01] ＞非常に保存性の良い，スキームであることが分かりました． もちろんオイラーをフルに解いた場合について解析してるんだよね？ まさか、一定速度で対流させた場合だけ、なんていうんじゃないよね？ 49 名前：名無 [2000/10/10(火) 17:28] ＞＃　ＣＩＰは，保存性が悪い！と言う論文が出ているのでしょうか？？ 衝撃波の解法としてはTVD系スキームより劣る、っていう論文なら在るよね。 Full paperには成ってないかもしれないけど。 50 名前：>38 [2000/10/12(木) 13:51] ＞非常に保存性の良い，スキームであることが分かりました． これは理論的に？ それとも実験的に？ 理論的にであれば OK, それはGOOD NEWSですよ！ しかし単にいくつかの限られた 問題に対して実験的に示されただけではやっぱり信用できない ですね。 ＞＃ＣＩＰは，保存性が悪い！と言う論文が出ているのでしょうか？？ 問題はCIPは保存則を満たすと理論的に証明されているのかどうか だと思います。それが保証されていないのなら、例え実際には問題 はなくとも、使おうという気にはなれませんよね。なぜなら 良い保存系スキームなど他に沢山存在するであろうから。 51 名前：プラズマン＝ランダウ [2000/10/12(木) 15:39] なんとなく、45は衝撃波解析と保存性の関係が理解できて いないような気がする。 52 名前：＞51 [2000/10/14(土) 00:51] そう言った考えを否定した論文もCIPの論文に書いてあります． ＃　具体的には「数値流体力学会誌の解説」です． ＃　何号かは失念しましたが．．．（一年以内のものです） 53 名前：＞50 [2000/10/14(土) 01:02] 「良い保存系スキームなど他に沢山存在するであろうから。 」 specialケースの問題を解いて，納得してください．．． ＃　私は，工学的な立場で考えますので， ＃「より単純で，汎用性があって，みんなに使ってもらえるスキーム」 ＃　を考えています．．． 54 名前：a [2000/10/14(土) 01:20] >52 どう言った考えを否定しているのかもう少し詳しく 書いて頂けますか？ 55 名前：猫八百 [2000/10/14(土) 01:26] >>52,>>53 言ってる事に全く具体性が無い。上に幾つも出ている「CIPの 欠点」についての指摘に具体的に答えてください。 56 名前：名無しさん＠１周年 [2000/10/14(土) 01:33] で、保存性の良し悪しはどういう例題で議論したのよ？>>52-53 57 名前：ちょっと質問 [2000/10/14(土) 01:42] >>52-53 §あなたは衝撃波問題を扱ったことがありますか？ §CIPにはノイマン型（もしくは、その改良版）の人口粘性が入って いますが、何故そのようなものが必要なのでしょうか？ 58 名前：名無しさん＠１周年 [2000/10/14(土) 02:24] >みんなに使ってもらえるスキーム 少なくとも、衝撃波問題に限って言えばTVD、ENO系スキームの方が 圧倒的に広く使われてるね。 59 名前：若葉名無しさん [2000/10/15(日) 11:58] kame.tadaima.com/2ch/ 60 名前：名無しサンシティ [2000/10/15(日) 12:36] >>59 なんでよりによってこのスレに？？ 61 名前：名無しさん＠１周年 [2000/10/17(火) 11:06] 戻ってこないのかな？ 62 名前：名無しさん [2000/10/17(火) 13:45] で結局のところ、圧縮・非圧縮性流体の統一解法は存在すんのかいな？ 63 名前：名無しさん＠１周年 [2000/10/18(水) 17:17] 圧縮性スキームで、time step をめちゃくちゃ細かくすれば何でも解けます 64 名前：>63 [2000/10/19(木) 14:32] それって本当？ それってＣＦＬナンバーを小さくするって 事だよね？ だけどそうすると一般的には数値粘性効果が絶大 になってろくでもない解を得ることになるのでは？ （注意： 人工粘性じゃなくて数値粘性だよ） 65 名前：名無しさん＠そうだ選挙にいこう [2000/10/19(木) 16:13] それはスキームによりけりだね>>64 66 名前：名無しさん＠ [2000/10/20(金) 08:42] ＞６５ へぇ〜。64の言ってることはおいらの知ってるスキームの全てに 言えることだからそうだろうなと思ってたけど。じゃあ、そうなら ないスキームって例えばどんなスキームか例挙げてもらえますか？ 67 名前：名無しさん＠ＣＩＰ [2000/10/25(水) 09:22] Prof. CIP www.mech.titech.ac.jp/~ryuutai/yabejp.html CIP SOFT CO. LTD. www.iijnet.or.jp/cipsoft/index.html 68 名前：名無しさん＠１周年 [2000/10/25(水) 10:58] 上のウェブでこんなん見つけた。 ”非保存形式における保存スキームの開発 （第１３回数値流体シンポジウム　1999年１２月）” www.mech.titech.ac.jp/~ryuutai/publist.html ＣＩＰを保存則を満たすように修正したと主張しているが。 69 名前：名無しさん＠１周年 [2000/10/25(水) 19:38] >>68 一定速度の対流とバーガーズだけ、ね。しかも１次元。 衝撃波の計算にどれだけ影響が有るものか。 まったく無いわけではないだろうけど。 メモリ使用量も増やし、これだけややこしい事をしつつ、 １次元onlyではね。多次元では更にややこしく成りそうだ。 70 名前：名無しさん＠１周年 mailto:sage [2000/10/25(水) 23:28] この方法は元々数学的仮定に不備があるのを，あとから無理矢理 こじつけようとしている感があるね． 71 名前：名無しさん＠１周年 [2000/10/26(木) 13:20] 数学的仮定の不備？ 72 名前：名無しさん＠１周年 [2000/10/26(木) 14:22] でもCIP SOFTなんて会社が成立してるぐらいだから CIPも捨てたもんでもないんじゃないの？ よく知らないけど。 ところで私はCIPの原理を勉強したことがないのですが、 誰か簡単に説明できますか？ 簡単にできなければいいですけど。 ちょっと興味あるな。 73 名前：名無しさん＠そうだ選挙にいこう [2000/10/26(木) 19:17] >>72 ttp://www.sciencedirect.com/science?_ob=IssueURL&_tockey=%23TOC%235301%231999%23998809997%23130395%23FLT%23Volume_119,_Issues_2-3,_Pages_99-316_(2_June_1999)&_auth=y&_acct=C000010898&_version=1&_urlVersion=0&_userid=132065&md5=bbc97c63330c1cf18940fab0d1b8265e ここの8番目の論文で引用されてる [1,　2, 17] を読むべし。 しかしこれ、人口粘性の最新の論文らしいけど、ちょっとヒドイね。 74 名前：名無しさん＠１周年 [2000/10/27(金) 22:48] >>72 親族系企業のようだから、少し差っ引いて考えた方が良いかもよ（笑 75 名前：名無しさん＠１周年 [2000/10/29(日) 00:14] わしどっちかっていうと、ソルバを用いた解析が主体なんだけど、 CIPは感覚的に捕らえやすい。ソルバとしてそれなりにタフだけど、系が閉じていないので場合によっては解の正確性に欠ける。 TVDは系が閉じてるけど、条件によっては簡単に解が出てこないときがある(衝撃波頭とか)と大まかに捕らえてるんですが、これ じゃあかんですか？ 76 名前：名無しさん＠１周年 [2000/10/29(日) 01:49] 意味わかんない 77 名前：>75 [2000/10/29(日) 11:09] ＞CIPは感覚的に捕らえやすい このセリフよく聞きますが、具体的に説明してもらえますか？ 感覚的にどう捕らえているのですか？ よろしくお願いします。 また、「系が閉じてる」の系とは何の系のことでしょうか？ もう少し具体的によろしく。 78 名前：＞77 [2000/11/07(火) 13:34] 物理量だけでなく，その勾配量も未知量として，解く． 79 名前：名無しさん＠そうだ選挙にいこう [2000/11/07(火) 14:21] それだけならエルミート法って言って昔から使われてる方法だよ。 80 名前：名無しさん＠１周年 [2000/11/07(火) 15:58] えっ！ CIPってHermiteなの？ なぁ〜んだ、ただのパクリじゃん！ 勾配も未知量ってことなら３次元では大変な数の未知量を扱う ことになるね。 81 名前：名無しさん＠そうだ選挙にいこう mailto:SAGE [2000/11/07(火) 17:11] >>80 上の方に有った「保存なんたらCIP」ってのは、さらに 未知量を増やすらしい… 82 名前：名無しさん＠１周年 [2000/11/10(金) 06:55] 一度，液滴変形のシミュレーションで使ってみたことがあるんですが， クーラン数をやたらと小さく取らないとうまくいかず，時間がかかりす ぎるという印象があります． 83 名前：名無しさん＠１周年 [2000/11/10(金) 13:42] >>82 他の方法では試した？ 84 名前：名無しさん＠１周年 [2000/11/11(土) 06:37] >>83 私の場合，Lax-Wendroff, Godunov's first order upwind, Roe's first order upwindと 比較したんですが，CIPは非常に小さいクーラン数を用いて長い計算時間をかけない と安定しませんでした． しかし，CIPのMLでは，他のスキームと比較して議論するのは好ましく無いというような ことをおっしゃられていた方がいますので，CIPを正しく評価するための検証法は他にある のかもしれません． 85 名前：83 [2000/11/11(土) 17:40] >>84 ちなみに、流体は非圧縮性ですか？ 表面張力のモデリングは何を用いました？ ＞他のスキームと比較して議論するのは好ましく無い どういう理屈でしょうね？？ 86 名前：名無しさん＠１周年 [2000/11/14(火) 11:49] >ちなみに、流体は非圧縮性ですか？ 液相，気相とも圧縮性で解いてますが，密度比を大きく取ってないので 若干ごまかしたような感じになってます． 衝撃波が液滴を通過する際の液滴変形挙動を追ってみようかと思って計算してみました． >表面張力のモデリングは何を用いました？ これに関しては厳密に考慮してません．界面における法線方向の境界条件に表面張力の 生データが入っているだけです． こんなとこなんですが，CIPはクーラン数の制約が大きい印象を受けたんで採用しませんでした． 87 名前：ここは黒人音楽板です mailto:ここは黒人音楽板です [2000/11/14(火) 18:09] >>86 界面の捕獲はどのように？「生データ」というのは？ 88 名前：名無しさん＠キティ立入禁止 [2000/11/24(金) 15:22] 結局CIPって大したスキームじゃないのね。論文を一つ見たけど、 （上にあった保存なんたらのやつね）、あまりの稚拙さに情けなくなっちゃった。 東工大ってこの程度なの？ いや、正確には、そのY教授ってそんな 程度なの？ って思っちゃった。日本人としては、日本発のスキーム ってことで期待したいんだけど、この人には期待できないです。 89 名前：あはは、 [2000/11/24(金) 20:17] 辛らつですな>>88 どのあたりが稚拙だと思いました？ 90 名前：Y [2000/11/28(火) 23:25] 早く答えろ>>88！ 91 名前：名無しさん＠そうだ選挙にいこう [2000/12/05(火) 22:55] Osherのグループに居た人が今度、新しいのを発表するよ。 92 名前：名無しさん＠キティ立入禁止 [2000/12/07(木) 16:34] >91 A NEW Central Schemeだろ？ Riemann問題を回避するってやつ。 ちがう？ 93 名前：名無しさん＠そうだ選挙にいこう [2000/12/08(金) 20:36] >>92 ちがう 94 名前：名無しさん＠キティ立入禁止 [2000/12/09(土) 10:39] ＞９３ じゃ、何よ？ どうせ何も知らないんだろ。 95 名前：名無しさん＠そうだ選挙にいこう [2000/12/09(土) 10:54] >>94　君、つまんないから出てって良いよ 96 名前：あぼーん mailto:あぼーん [あぼーん] あぼーん 97 名前：↑↑↑↑↑↑↑↑↑↑↑↑↑↑↑↑↑↑↑↑↑↑↑↑↑↑↑ mailto:sage [2000/12/09(土) 13:35] あ〜あ、やっちゃったね。このスレ。 せっかく盛り上がってたスレッドなのに、あんたらのレスで台無しだよ。 なんでここで、そんなレスしかできないわけ？ 空気読めないの？ だからあんたらは駄目な奴だって言われるんだよ。 わかってるの？ それにしても、もったいないなあ、せっかくここまで育ったスレッドなのに。 ここまででおしまいかよ。 まあ、しかしやっちゃったものはしょうがない。 これからはもうちょっとマシなレスするように心がけろよ。 98 名前：名無しさん＠１周年 [2000/12/10(日) 02:45] もう終わりなの？ もちょっと続けてよ。 99 名前：名無しさん＠そうだ選挙にいこう [2000/12/15(金) 21:18] Fedkiw 100 名前：名無しさん＠１周年 [2000/12/15(金) 22:59] 今月の数値流体シンポに出てくるのでは？＜矢部先生のグループ それに期待しませう 101 名前：名無しさん＠１周年 [2000/12/16(土) 01:38] ＞９９ Fedkiwか。会ったことある。お化け流体法だろ。 102 名前：名無しさん＠そうだ選挙にいこう mailto:sage [2000/12/17(日) 20:33] >>100　毎年出てるだろ（ﾜﾗ 103 名前：名無しさん＠１周年 [2000/12/26(火) 18:36] で、どうでした？今年の数値流体シンポは。 104 名前：名無しさん＠１周年 [2001/01/04(木) 16:03] ＞１０３ 大したことなかった。いつものようにただ集まっただけよ。 105 名前：遅すぎるレス [2001/01/18(木) 21:18] >>64 ＞それって本当？ それってＣＦＬナンバーを小さくするって ＞事だよね？ だけどそうすると一般的には数値粘性効果が絶大 ＞になってろくでもない解を得ることになるのでは？ これに関する文献、誰か知らないかな？ 106 名前：名無しさん＠１周年 [2001/01/19(金) 12:29] Advectionに関しては普通はそうなる。これに関しては CFDの本ならなんでも載ってるはずだよ。 あと、大学でＣＦＤを学んだのならやってるはずだよ。 そう、フーリエ分解して複素数でなんじゃかんじゃするやつ。 Dissipation Errorだよね。これって(Dispersionもそうだけど) ＣＦＬナンバーに依存するんだよね。 107 名前：名無しさん＠そうだ選挙にいこう [2001/01/19(金) 17:52] おかしいな、僕の使ってる対流スキームはCFLを小さくしても 精度が落ちない…。比較的普通のスキームだと思うけど…。 108 名前：107 [2001/01/19(金) 18:12] CFL＝0.1で2,000step計算させた結果と CFL＝0.001で200,000step計算させた結 果が変化ほとんど無し。 109 名前：名無しさん＠１周年 [2001/01/27(土) 15:26] >107 それって具体的にどんなスキーム？ まさか線形じゃないよね？ 106の言ってるのは線形の スキームだと思うよ。例えば、Lax-Friedrichsとか FTCSとか。なんて思うがね。 110 名前：名無しさん＠そうだ選挙にいこう [2001/01/29(月) 20:39] 意味わからん 111 名前：名無こ [2001/01/30(火) 00:00] >>109 どういう意味で線形？ 112 名前：名無しさん [2001/01/30(火) 15:14] ＞１１０，１１１ ってことは、Godunovの定理とかvon Neumannの安定解析 とか知らないんだね。勉強するべし。 ちなみに、Godunovの定理とは「線形でOscillation-freeな スキームはせいぜい１次精度でしかない」。つまり、高精度の non-oscillatoryスキームは非線形でなければならない。 そこで、流束制限関数なんかが登場するわけ。 （線形ってのは、ある格子点での数値解がその前の時間ステップの 解の線形結合であるという意味。そうでなければvon Neumannの安 定解析が使えない。それを克服するためにTVDという非線形の 安定解析が開発された。） 113 名前：名無しさん＠１周年 [2001/01/30(火) 20:58] は？線形でも非線型でも、CFL依存性が少ない スキームは在るけど？ 114 名前：名無しさん＠１周年 [2001/01/30(火) 20:59] 線形だとCFLが小さい時に拡散的になって、 非線型だとそうはならない？？ なんか話が変だと思うな。 115 名前：名無しさん [2001/01/31(水) 06:21] ＞線形でも非線型でも、CFL依存性が少ないスキームは在るけど？ 誰も「線形だからCFL依存性がある」とは言ってないですね。 線形だからその依存性が解析的に調べられるってことですね。 もちろん依存性の少ないものもあるでしょう（LPEスキームとか）。 基本的なスキームの解析はどんなCFDの本にも書いてます。 本をご覧ください。 ＞線形だとCFLが小さい時に拡散的になって、非線型だとそうは ＞ならない？？ そんなこと誰も言ってないですね。非線形だと解析が難しい からCFL依存性を解析的に調べるのが困難だから、わからない と言ってるだけですね。（ちなみにTVDを用いてもdissipationや dispersionはわからない。単に安定性しかわからない。） 要は非線形だとやってみないとわからない、ってことですね。 まず本を読んで勉強してくださいね。例えばLax-Friedrichs スキームのdissipation, dispersion errorのCFL依存性を まず安定解析で導出してください（大抵の本には書いてますが）。 それから実際にscalar advectionを解いて見てください、 いろんなCFLで。CFL＝０．０１程度でどんな解もほぼペチャンコ になります。 とにかくまず本を読みましょう！ それじゃ！ 116 名前：名無しさん＠そうだ選挙にいこう [2001/02/01(木) 18:23] じゃ結局>>106はLax-Friedrichsの事なんすね。 117 名前：名無しさん＠１周年 [2001/02/01(木) 19:48] >>106の「普通は」と>>109の「まさか」の意味が 何度読んでも分からないね。 118 名前：名無しさん＠１周年 [2001/02/02(金) 00:50] そればかりは幾ら本を読んでも載っていまい… 119 名前：名無しさん＠キティ立入禁止 [2001/02/02(金) 10:25] >116,117,118 君達、つまんないから出てって良いよ 120 名前：あぼーん mailto:あぼーん [あぼーん] あぼーん 121 名前：名無しさん [2001/02/02(金) 10:33] あれっ、とうとう終わっちゃったね、このスレ。 前は結構楽しんで読んでたんだけど。 まぁ、しょうがないか。 内容のレベルもかなり下がってきてたしね。 また、誰か何か書いてね、時々はチェックするからさ。 122 名前：あぼーん mailto:あぼーん [あぼーん] あぼーん 123 名前：あぼーん mailto:あぼーん [あぼーん] あぼーん 124 名前：あぼーん mailto:あぼーん [あぼーん] あぼーん 125 名前：あぼーん mailto:あぼーん [あぼーん] あぼーん 126 名前：あぼーん mailto:あぼーん [あぼーん] あぼーん 127 名前：あぼーん mailto:あぼーん [あぼーん] あぼーん 128 名前：あぼーん mailto:あぼーん [あぼーん] あぼーん 129 名前：名無しさん＠１周年 [2001/02/02(金) 12:02] age 130 名前：57 [2001/02/05(月) 14:13] いつになったら答えてくれるのだろうか… >>57 131 名前：非公開＠個人情報保護のため [2001/02/05(月) 17:04] 2次元CIPCSLプログラムを構築し，回転流動場に適用したけど，うまくいかん． 132 名前：名無しさん＠１周年 [2001/02/06(火) 00:15] もっと情報を提供しなければ まともなレスは無いと思いなさい… 133 名前：名無しさん＠１周年 [2001/02/06(火) 16:43] 132さん、まぁまぁ、そう言わずに。 131さん、うまくいかないのはコードにバグがあるからだよ。 または、CIPCSLとやら自体が駄目なんだよ。どっちかなんだよ。 取り敢えずどっちの問題なのか調べてみてね。 134 名前：131 [2001/02/07(水) 13:28] かなり複雑なコードを必要とするんだよ，2DのCIPCSLって． だから普及しないとおもう．3Dなんて全くもって無理． 結局理論はいいんだろうけど， 実用厳しいんじゃダメなんじゃないか？ 135 名前：名無しさん [2001/02/07(水) 13:47] なるほどね。でも理論的背景が本当にしっかりしてるんだったら 将来はポピュラーになる可能性はあるね。圧縮性のほうのＭＵＳＣＬ とやらも昔は厄介と思われてたらしいけど、やっぱ理論がしっかり しているということと、コンピューターの速度の増加によって 標準になったもんね。 136 名前：名無しさん＠１周年 [2001/02/07(水) 13:59] どんな計算条件を使ったのかが全然分かんない。 137 名前：名無しさん [2001/02/07(水) 13:59] ＞結局理論はいいんだろうけど どのあたりが？ 138 名前：134 [2001/02/07(水) 15:41] >>137 物理量の保存を全領域(グリッド点だけでなくて，グリッド内も)に わたって，満たすようにしようという点。 ただ，適用する方程式が限られてる気がする。 運動方程式には無理ちゃうか。 139 名前：名無しさん [2001/02/07(水) 16:07] >>138 ＞グリッド点だけでなくて，グリッド内も 保存CIPの場合はグリッド（っていうかセル）内だけでしょ？ でもそれ（「グリッド点だけでなくて，グリッド内も」）、 べつに特別じゃないでしょ。 例えば、簡単に線形補間でプロファイルを見積もったとしよう。 この時、セル内でその線形関数を積分し、それを全セルに渡って 足し合わせても結果は �杷_{i}Δx だよね。たぶん高次関数で 補間した場合もこうなるよ。 140 名前：名無しさん [2001/02/07(水) 16:12] 上のほうに在った保存CIPの論文（日本語）の中で 「セル内においてfが一定であるという仮定を利用」 しないと �杷_{i}Δx にはならない、みたいな記述が 在るけど、>>139のような理由でそれは嘘です。 141 名前：138 [2001/02/07(水) 16:58] うん，139は間違ってるね。 グリッド内でf＝一定にはしないです。 高次補間関数を積分しているので。 142 名前：名無しさん＠そうだ選挙にいこう [2001/02/07(水) 17:36] 結局CIPCSLとやらには特別な点は見当たらないという事ですか。 143 名前：名無しさん＠１周年 [2001/02/07(水) 17:38] 一応，保存性を満足しているのでは？←ｱｲﾏｲ 144 名前：名無しさん＠そうだ選挙にいこう [2001/02/07(水) 18:06] 保存型スキームを使えば普通じゃん、それ 145 名前：139 [2001/02/07(水) 18:08] ん？意味分かってる？>>141 146 名前：141 [2001/02/07(水) 19:13] >>145 理解してますよ． たしかに「グリッド点での物理量」保存しません。訂正します。 >高次関数で補間した場合も・・・ というのは，意味がわかりませんが。 CIPCSLはグリッド点上の物理量fの全グリッド点和が 保存されるのではなく， セル内のfの総量∫fdxの全領域和が保存されるわけです． すなわち，全領域での保存性が成り立つと言うことです． (日本語おかしいですか？) で，�杷_{i}Δx ではないです． 147 名前：名無しさん＠１周年 [2001/02/07(水) 19:23] >>144 保存型スキームってのは収束するまで反復計算するやつのことを いっているのかな？ 148 名前：139 [2001/02/07(水) 19:55] >>146 残念ながら、理解してないよ。もう一度だけ繰り返そう。 離散空間{x_{1}, x_{2},...,x_{n}}上に f_{i} (i=1,...,n) が定義されているとする。区間[x_{i}, x_{i+1}]内で f を 線形関数で補間し（一定とするわけではない）、その区間での 積分をとれば ρ_{i+1/2} = (f_{i}+f_{i+1})Δx/2。 ここでΔx≡(x_{i+1}-x_{i})とした。 ρ を全領域で総和すれば �買ﾏ = �杷_{i}Δx。 高次関数を使って同じ計算をしても同様の結果を得る。 つまり、普通のスキームでも∫fdx は保存されている。 上で示したように f のみを用いたスキームでは ∫fdx = �杷_{i}Δx なんだから。 149 名前：139 [2001/02/07(水) 19:58] これ↑常識。 ∫fdx が保存されるのはCIPCSLだけ、というのも、セル内で f を一定とした場合にだけ∫fdx = �杷_{i}Δx となる、という のも嘘。 150 名前：139 [2001/02/07(水) 20:01] Δx = 一定という仮定が落ちてたな>>148 151 名前：146 [2001/02/07(水) 20:25] 「Σ∫fdx」と「 �杷_{i}Δx 」の違いは分かってますよね？失礼。 結局CIPCSLってのはセル内の総物理量の保存を考えるってことで， セルが大きくなっても精度よく計算できると考えることができると 思うんですが・・． さらに，1Dはまだしも２Dの場合には(方法にもよりますが)， グリッド内の物理量の分布形状を精度良く(あまり使いたくない言葉) 表現することができるのですよ． 線積分とか面積分とかを用いて． 152 名前：＞151 [2001/02/07(水) 20:35] たしかに君の方が分かっていない。 153 名前：＞151 [2001/02/07(水) 20:36] 理系だからといって言語力が低くて良いというわけではありません。 154 名前：151 [2001/02/07(水) 21:01] あれ，伝わりませんでしたか・・ ちなみにずっと間違えていましたが， ∫fdxではなく，∫Fdxでした． Fはグリッド内の補間関数です。 155 名前：139 [2001/02/07(水) 21:16] >>151,>>154 君は他人の発言の意図を読むのが苦手なようだね。 >>148をもう一度良く読もう。ちゃんと補間関数の 積分について議論し、答えを出しているから。 なんだったら2次や3次関数を用いた場合について の証明を自分でしてみれば良い。 >>152-153 気持ちは良く分かる。 156 名前：151 [2001/02/07(水) 21:30] ご迷惑をおかけしましたね。 自分の中で解決しました(最悪のパターン？) CIP法を発展させて(退化してる気もする)， 保存性を向上させたものがCIPCSLということ？ (ん？まだわかってない？) ちなみに，CIPかCIPCSLを実際にプログラムに 組みこんでるひといます？ 157 名前：＞１３９ [2001/02/08(木) 14:34] ＞ρ を全領域で総和すれば �買ﾏ = �杷_{i}Δx。 ＞高次関数を使って同じ計算をしても同様の結果を得る。 これは間違い。例えばHermiteを使うとスロープが 入ってくるから式は変わる。「高次関数」ってどんな のを想定しているのかしら？ 「２次や３次関数」って どうやって作るの？ 確かに区間[x_{i}, x_{i+1}]内で f_{i},f{i+1},f_{i+2}を使えば２次関数が作れる が、その積分を全ての区間で足しても�杷_{i}Δx にはならないよ。（これは実際に確かめた。） ＞つまり、普通のスキームでも∫fdx は保存されている。 これも間違い。保存されているかどうかは、Fluxが telescopingかどうかであって、ある物理量の全総和 が�杷_{i}Δxで表されることとは何の関係もない。 ある時間ｔで�杷_{i}Δx＝１だったとして、 次のタイムステップで�杷_{i}Δx＝１のままかどうか が問題なんでしょ。 158 名前：139 [2001/02/08(木) 21:34] >>157 ＞例えばHermiteを使うとスロープが〜 知ってるよ。だから、CIPは定速度場上でも完全には 保存しない。上では保存型CIPと、fのみを使った保存 型スキームとの比較をしているつもり。 ＞が、その積分を全ての区間で足しても�杷_{i}Δx〜 なるけど？どういう計算をしたの？？ ＞＞つまり、普通のスキームでも∫fdx は保存されている。 これも上同様、保存型スキームについての話。だって、 保存型スキームの話をしてるんだから。 159 名前：139 [2001/02/08(木) 21:39] >>151 ＞結局CIPCSLってのはセル内の総物理量の保存を考えるってことで， ＞セルが大きくなっても精度よく計算できると考えることができると っていうか、格子上に置いてる変数の数が増えてるん だから、精度が良くなって当然なんじゃ？だって、f のみを使ったスキームに置き換えて考えればつまり、 格子数を増やしたってことに当たるんだから。 それ以上の効果が有るんなら面白いけど。 160 名前：139 [2001/02/08(木) 21:43] >>156 いや、もともと保存しないスキームを保存するよう にしたんだから一応 発展 なんじゃない？ でも、こうする事で、これまでCIPの利点とされて きた特性までもが消えてしまいそうな気もするが…。 161 名前：>158 [2001/02/09(金) 03:47] >＞が、その積分を全ての区間で足しても�杷_{i}Δx〜 >なるけど？どういう計算をしたの？？ 単に積分して足しただけだよ。境界から十分離れていれば 確かに�杷_{i}Δxと書けるが、どうしても境界付近の 値がf_{i}Δxにならない。ｃf_{i}Δxのように１でない 定数cが掛かってくるよ。つまり�� c_i f_{i}Δxのように は書ける。 162 名前：157ではないです [2001/02/09(金) 05:02] >>161 境界がどうのって言ってるんじゃないと思うよ. >単に積分して足しただけだよ。 この「積分」ってのは解析的に積分したもののことをいってる？ それとも数値的に積分したもののこと？ 163 名前：139 [2001/02/09(金) 13:52] >>162 ＞境界がどうのって言ってるんじゃないと思うよ. あたり。 >>161 普通、保存してるかどうかを確かめる時には周期境界を 仮定するんだよ。だって、境界から余分なフラックスが 流れ込んできてたら、どんなスキームだって保存しない じゃない。 周期境界を仮定すれば、片方の境界から流れ出たフラッ クスと同じ量のものがもう片方から流れ込む。 164 名前：> [2001/02/09(金) 16:42] ＞この「積分」ってのは解析的に積分したもののことをいってる？ いってる。積分と言えば普通、解析的なものを指す、と思うが。 どうでしょう。 なるほど、139の言ってることを理解した。 普通のスキームってのは discrete conservationを満たす スキームってことで。 言い換えれば�杷_{i}Δxを保存するスキーム。 で、高次関数を使ってセル内で解を再構築してそれを積分和し たものも�杷_{i}Δxとなると。（セル内での積分値はもちろん 異なるが足せば同じになると）故に、�杷_{i}Δxを保存する スキームは高次関数で再構築されたfの積分和も保存しているこ とになる。高次関数を使ったのだからその積分和のほうが精度が いいと思ってしまいがちだが、そうではないと。と言う訳で、 �杷_{i}Δxを保存すれば∫（高次関数補間）dxも保存していると。 165 名前：>>> [2001/02/09(金) 16:47] ＞＞境界がどうのって言ってるんじゃないと思うよ. ＞あたり あなたがそういうことを言ってるとか言ってないとかの問題じゃな くて。境界があればquadratic reconstructionでは積分和が �杷_{i}Δxにならないってことは事実。あなたが境界がどうのこ うのと言ってないだけに、私が指摘しただけですよ。それだけですよ。 ＞普通、保存してるかどうかを確かめる時には周期境界を仮定するんだよ。 保存してるかどうかはスキームを見ればわかると思うが、まぁ計算 してみないと保存してるかどうかわからないスキームもあるのだろう。 なるほど、スキームは明らかに 保存ではないけど、実際計算してみて大体保存していると後で チェックするってことはあるだろうな。 ＞境界から余分なフラックスが流れ込んできてたら、どんな ＞スキームだって保存しないじゃない。 保存というのは一定という意味ではなくて、保存量の変化が 境界の値にのみ依存するという意味です。それが保存則。境界から 余分なフラックスが入ってきたり出ていったりすれば当然その差し 引きの分だけ増加又は減少するはずです（これが保存則）。保存型 スキームはそのような変化も数値的に再現できなければならない。 だからフラックスのtelescopingが必要なんですね。で、要は境界 があれば164にあるようなことは厳密には正しくないってことです。 （ところで、境界を考慮することは重要ですよ。境界のない問題なん て、Model Problems以外じゃほとんど皆無です。） 166 名前：139 [2001/02/09(金) 17:43] >>165 よくよく思い返してみたら、ちょっとした疑問が。 >>157 ＞ある時間ｔで�杷_{i}Δx＝１だったとして、 ＞次のタイムステップで�杷_{i}Δx＝１のままかどうか ＞が問題なんでしょ。 これで「境界の影響は無視する」ってことで了解が 取れていると思ったんだが。 167 名前：> [2001/02/09(金) 17:46] なるほど。 168 名前：139 [2001/02/09(金) 17:46] >>165 ＞＞どんなスキームだって保存しないじゃない。 これは言葉足らずだった。意図したのは、 　　　　どんなスキームだって�杷_{i}Δxは保存しないじゃない。 ということ。 169 名前：> [2001/02/09(金) 17:50] なるほど。 170 名前：> [2001/02/09(金) 18:02] という訳で、160にあるように、保存ＣＩＰは別に他の保存型 スキームよりも保存性が格別に良いとは言えないが、しかし、 保存型ではなかったＣＩＰが保存性を装備したというこ とは一つの発展である、ってことですかね。 しかしなぜ保存性がそんなに大切なのでしょうかね、皆さん。 171 名前：156など [2001/02/09(金) 18:24] 移流項の計算にCIP法を用いたCCUP法などで，DAMBREAKの計算を行なうと， 水柱の隅角部で変な現象が起きます．（やったことある人には分かる） しかし人工粘性を用いるとこの現象が消えていい結果が得られるのです. しかし，CIPCSLは「厳密な保存性を保証するため」人工粘性がなくても， この現象が起きないと言う触れ込みだったと思います． 172 名前：156など [2001/02/09(金) 18:58] >この現象・・・ 「こういう」に訂正します. 173 名前：名無しさん＠１周年 [2001/02/12(月) 19:58] ＞人工粘性を用いるとこの現象が消えていい結果が得られるのです. ＞CIPCSLは「厳密な保存性を保証するため」人工粘性がなくても， ＞この現象が起きない ということは人口粘性でも保存性を保てるということかな？ 違うだろうけど、保存性と人工粘性の関係がよくわからないな。 174 名前：amigo [2001/02/13(火) 16:46] CIPの利点、簡単なのでどんな方程式にも使いやすい。でも、理論的根拠 は、ほとんどない。なぜなら、音速をオイラーで陰的に解いている。これでは、 ほとんど精度がでない。計算の精度は、もっとも低い精度の部分で決定されるはず。 また、時間の刻みを u だけで決めている。ということは、u+cs,u-csの部分は 完全に捕らえられていないはず。非線形の場合は、陰解法でも、max(｜u｜, ｜u+cs｜,｜u-cs｜)からdtを決定する必要がある。線形方程式の安定性理論と 非線形方程式の安定理論は、違う。 175 名前：名無しさん＠そうだ選挙にいこう [2001/02/13(火) 18:28] ＞理論的根拠は、ほとんどない。 ＞なぜなら、音速をオイラーで陰的に解いている。 このあたりの話の繋がり方が理解できません 176 名前：名無しさん＠そうだ選挙にいこう [2001/02/13(火) 18:29] ＞計算の精度は、もっとも低い精度の部分で決定される ここも意味不明です 177 名前：名無しさん＠そうだ選挙にいこう [2001/02/13(火) 18:30] >>171 ＞人工粘性を用いるとこの現象が消えて どの人口粘性？？ 178 名前：amigo [2001/02/13(火) 19:36] CIPでは、advection , non-advectionに分けるでしょう。 前者は、時間と空間の両方３次精度、後者は時間１次精度。精度の高い計算と 悪い計算を組み合わせて使っても、結局、悪部分で決定される。つまり、時間は １次精度なわけ。 A=(真の値）＋（１次の誤差） B=(真の値）＋（３次の誤差） A+B=(真の値）＋（１次の誤差）＋（３次の誤差） 　　=(真の値）＋（１次の誤差） A*B=(真の値）＋（１次の誤差） 179 名前：名無しさん＠そうだ選挙にいこう [2001/02/13(火) 19:38] >時間と空間の両方３次精度 時間は1次 180 名前：名無しさん＠そうだ選挙にいこう [2001/02/13(火) 19:41] >A+B、A*B この+と*がまた意味不明なんだが 181 名前：名無しさん＠１周年 [2001/02/13(火) 20:40] ＞計算の精度は、もっとも低い精度の部分で決定される ここはわかる。エラー解析の基本。 だけど、それ以外のところは私にもわかりませんです。 182 名前：非公開＠個人情報保護のため [2001/02/13(火) 21:08] >>177 よくCIPの論文に載ってるヤツ >>178 初期のCIPの論文見たことある？ 183 名前：名無しさん＠そうだ選挙にいこう [2001/02/13(火) 21:15] なんでランキン・ユゴニオから求められる、ショック・フロントに 必要な人口粘性が、ダム・ブレイク程度の問題でも必要になるのか。 謎だ。 184 名前：名無しさん＠そうだ選挙にいこう [2001/02/13(火) 21:16] ＞しかし，CIPCSLは「厳密な保存性を保証するため」人工粘性がなくても， ＞この現象が起きないと言う触れ込みだったと思います． もうこんなアプリケーションを行なってるんだろうか。 見たこと無いが。 185 名前：名無しさん＠そうだ選挙にいこう [2001/02/13(火) 21:17] ＞計算の精度は、もっとも低い精度の部分で決定される 空間に関してはそう単純でもないけどね 186 名前：名無しさん＠そうだ選挙にいこう [2001/02/19(月) 19:50] なんか中途半端なトコで止まったなぁ。 187 名前：名無しさん＠１周年 [2001/02/20(火) 13:35] ＞空間に関してはそう単純でもないけどね 単純だよ。なんで単純じゃないの？ なにか 実例があればよろしく。 188 名前：名無しさん＠１周年 [2001/02/20(火) 14:31] >187 隣接2点で3次精度なのか， 隣接4点で3次精度なのか， 189 名前：名無しさん＠１周年 [2001/02/20(火) 18:10] >>187 対流項を3次、拡散項を2次； 対流項を2次、拡散項を2次： この両者による結果の圧倒的な差を無視して、 両者を「2次精度」とひとくくりにして良いのか 190 名前：187 [2001/02/20(火) 18:42] >188 それは離散化の仕方によって変わるでしょう。 つまり、２点でもＨｅｒｍｉｔｅのように 傾きまで考慮すれば４次精度になるし、 ４点使えば傾き無しで３次精度にできるでしょう。 要は、その後のはなしでしょう。離散化を選んで、 さて、一番低い精度が支配的になるか、と。 >189 その「圧倒的な差」というものを知らないので 無視するつもりはないですけど。具体的にそれはどんな 差なのでしょうか？ うーんと、ちょっと考えても、前者は いくら対流項が３次でも、拡散が２次なら、、、 やっぱ結果として２次だと思うなぁ〜。（もちろん 解が特異点などをふくむ場合は理屈通りにはいかない のはわかるけど・・・） もちろん、前者のほうが後者 に比べて誤差が小さいであろうことは想像できますが。 やっぱりその誤差は２次のオーダーで落ちていくのでは ないでしょうか？ とりあえず、その「圧倒的な差」について 説明して頂けますか？ よろしくお願いします。 191 名前：名無しさん＠そうだ選挙にいこう [2001/02/22(木) 19:32] >>その「圧倒的な差」について説明して頂けますか？ 少し自分の手を動かそう。 192 名前：名無しさん＠１周年 [2001/02/22(木) 19:47] >>その「圧倒的な差」について説明して頂けますか？ ＞少し自分の手を動かそう。 なるほど、これじゃ議論にならないわけだ。 もういいよ、このスレ。ほんとつまんなくなっちゃったよ。 193 名前：名無しさん＠１周年 [2001/02/22(木) 19:53] >>190 ＞つまり、２点でもＨｅｒｍｉｔｅのように ＞傾きまで考慮すれば４次精度になるし、 ３次。 194 名前：名無しさん＠１周年 [2001/02/22(木) 19:54] >>192 ＞もういいよ、このスレ。ほんとつまんなくなっちゃったよ。 自分の頭で考えずに何でも聞こうとする人が居るから。 195 名前：名無しさん＠１周年 [2001/02/23(金) 08:34] >>193 ４次精度だよ。少し自分の手を動かそう。 196 名前：名無しさん＠１周年 [2001/02/23(金) 08:39] さぁ、みんなで手を動かそう！ 197 名前：名無しさん＠１周年 [2001/02/23(金) 12:32] 手、手、手、 198 名前：名無しさん＠１周年 [2001/02/23(金) 15:02] よーし、それじゃぁ俺も動かすとするか。それっ！ シャカ シャカ シャカ シャカ シャカ シャカ シャカ シャカ 199 名前：名無しさん＠１周年 [2001/02/23(金) 20:32] >>195 やっぱりそうか。君、対流方程式を偶数次精度の スキームで解くとどういう事が起こるのか知らな いんだね、要するに。 200 名前：名無しさん＠１周年 [2001/02/24(土) 11:50] ほんともういいよ、このスレ。ほんとつまんない。 特に199君, 君がいろいろよく知っているであろうこと は分ったからさぁ〜。なんだかんだと根拠のかけらも示さず言っておきながら、何ですかって聞かれたら 「自分で考えろ」だって。「圧倒的な差」とか「偶数精度の・・・」 とか、何が言いたいのかはっきりしないよ。圧倒的な何の差だ？ とりあえず、議論をしようとは思ってないだろ？ そうじゃなきゃそんな書き方はしないよな。「自分の頭で考えずに何でも聞こうとする人が居るから。」だって？ 190を見るかぎり少しは考えてるように見えるけどね、少なくとも俺には。 ちなみに195は明らかに皮肉ってるだけだろ。 こんなこと書いたらまた何かあやふやな表現で返してくるんだろ。 お好きにどうぞ。もう誰も君に質問してこないと思うよ、安心しな。なんか一昔前 の大学教授を思わせるね、君は。「そんなこともしらんのか！」「自分で考えろ！」とかね。実は自分もよく知らないってのが大半だったね。後は知ってるけどただでは教えない嫌〜なタイプ。 201 名前：名無しさん＠１周年 [2001/02/24(土) 13:49] >なんか一昔前の大学教授を思わせるね まだまだいますって、そういう教授・・・鬱 202 名前：名無しさん＠１周年 [2001/02/26(月) 17:15] >>200 つまんないなら来なくて良いよマジで 203 名前：名無しさん＠１周年 [2001/03/01(木) 21:34] なんかJCPに載ったね、圧縮流と非圧縮流を一緒に解く とかっていうの 204 名前：名無しさん＠１周年 [2001/03/05(月) 14:18] >>190 ＞つまり、２点でもＨｅｒｍｉｔｅのように ＞傾きまで考慮すれば４次精度になるし、 ＞４点使えば傾き無しで３次精度にできるでしょう。 比較のし方が間違ってるね。 205 名前：名無しさん＠１周年 [2001/03/06(火) 11:50] >204 いや、間違ってないね。 206 名前：名無しさん＠１周年 [2001/03/06(火) 12:43] 3次関数使って4次になるわけないじゃん。 どうせ f_{n+1} = f_{n} + O(dx_{4}) で、 やったー4次だ！なんて思ってるんだろう（笑 207 名前：名無しさん＠１周年 [2001/03/06(火) 12:49] >>170 ＞しかしなぜ保存性がそんなに大切なのでしょうかね、皆さん。 衝撃波。凝縮や化学反応等による質量転移。非常に多くのタイム ステップを費やす必要のある問題各種。 208 名前：名無しさん＠１周年 [2001/03/06(火) 15:04] >206 何笑ってるの？ f_{n+1} = f_{n} + O(dx_{4})だったら４次精度じゃん。 209 名前：名無しさん＠１周年 [2001/03/06(火) 15:06] ＞２０７ それらの問題において、”なぜ”保存性が大切なのでしょうかね？ 皆さん。 210 名前：名無しさん＠１周年 [2001/03/06(火) 15:07] ふわぁ〜〜〜〜あ。 すまん。 あくびした。 211 名前：名無しさん＠１周年 [2001/03/06(火) 15:11] 数値流体関係の掲示板みたいなのどっかにあったら誰か教えてください。 212 名前：名無しさん＠１周年 [2001/03/06(火) 15:11] >>208 大笑い。ネタならネタと言ってくれよ、相手せずに済むから。 213 名前：名無しさん＠１周年 [2001/03/06(火) 15:26] 実質初心者板だね 214 名前：名無しさん＠１周年 [2001/03/06(火) 16:51] 確かに質は落ちているね。いや、最初からこんなもんだったと言うべきかな。 215 名前：208 [2001/03/06(火) 16:58] ネタじゃないですよ。私は純粋にそう思うんですけど。あきれて物が言えないってのなら、何も言わなくても結構です。ごめんなさい。 216 名前：仕様書無しさん [2001/03/07(水) 13:47] >>211　日本じゃやってる人間少ないから　無いんじゃない？ メーリングリスト探したほうがいいかも 217 名前：名無しさん＠そうだ選挙にいこう [2001/03/07(水) 15:53] 欧米ならCFDとかCEMの学科が在ったりするのにな〜 218 名前：流体 mailto:sage [2001/03/08(木) 06:04] >>211 www67.tcup.com/6730/netu.html 219 名前：名無しさん＠１周年 [2001/03/20(火) 23:24] >>217 CFDに絞ったcurriculumがあるって事？欧州は知らないけど、北米じゃ 聞いたことないよ。 220 名前：名無しさん＠京都板じゃないよ [2001/03/24(土) 01:17] ついに登場、声で会話できる2ch風のチャットだＹＯー。 ニチャンネラの生ボイスをきかせてＹＯ！ がんがん、みんな語っちゃオウＹＯ。 まぢに、ニチャンネラ同士で、朝まで討論会してみないかいぃ？ 最大15人くらいまで同時会話できるっぽい！ どんどん。きてくれＹＯ！ 2ch-voice.e412.com 221 名前：age mailto:sage [2001/03/27(火) 02:40] もっと勉強しましょう。このスレもっと早くくれば良かったな。 終わってやんの。 222 名前：名無しさん＠１周年 [2001/03/29(木) 17:02] ageなきゃ終ったままだよ（ｗ 223 名前：age [2001/03/30(金) 11:25] age age age age age ago ago ago ago?! 224 名前：名無しさん＠１周年 [2001/03/31(土) 18:03] >>218 なんか書き込み辛い雰囲気だなぁ〜 寺子屋のまん中で踊れっていわれてるようだ。 225 名前：名無しさん＠１周年 [2001/05/10(木) 17:12] フェドさんの論文読んだよ。なるほどね。 226 名前：日本アメリカ化計画 [2001/05/13(日) 11:53] Ａ型は徒党を組んで国民を操ろうとする。注意せよ！ 227 名前：名無しさん＠１周年 [2001/05/15(火) 19:36] >>221 で、何が言いたかったの？ 228 名前：名無しさん＠１周年 mailto:sage [2001/05/16(水) 02:13] みなさん、はじめまして、わたしは１の妹です。 おにいちゃんが今朝自殺しました 庭の物置で首を吊っていました 数日前からブツブツと何か独り言を言ったり、 更生がどうとか、天から選ばれたとか、ご飯も食べないし 家族中が心配していたんです。 そしたら、こんなことになっちゃって…しんじられない 遺書のようなものがあって、そこにはこう書かれていたんです。 俺の渾身の投稿も空振りに終わった…、最後に命燃やし尽くすことは出来なかった…と そしてブックマーク、一番上と書いてありました。 その意味が何なのか知りたくてさっきネットに繋いでみたらここだったんです。 まさかおにいちゃんがこんなことで自殺したなんて…わたし… 遺書はわたしが見つけて、そのことを警察の方にも両親にも知らせてませんでした。 こんな事実を知った後となっては、なおさら知らせるわけにはいきません。 おにいちゃんはなんでこんなことしたんだろう…わたしは妹としてそれを知りたい。 だから、わたしはしばらくこの２ちゃんねる掲示板を観察しようと思います。 なぜ、おにいちゃんはここまで追い詰められたのか…その答えを探しにいきます。 おにいちゃん、これでいいよね…、おにいちゃんが１2歳の誕生日にくれた ふでばこ今も使ってるんだよ…おにいちゃん、本当にありがとう 229 名前：名無しさん＠１周年 [2001/05/24(木) 13:09] 某計算法のドンは基礎知識に問題があるね 230 名前：名無しさん脚 [2001/05/24(木) 17:33] >>229 日吉方面ですか？ 231 名前：名無しさん＠１周年 [2001/05/24(木) 19:52] ちがうだろ 232 名前：名無しさん＠１周年 mailto:sage [2001/05/26(土) 06:35] １の主治医です。 この度、このようなスレッドを１が立てるに至ったことは、 主治医として、大変残念な事であり、また、治療の効果が まだまだ現れていないことを証明しているため、そろそろ 最終的な決断を下す必要があるようです。 みなさんお聞きになったことがあるかもしれませんが、 必ずしも心の病は、特殊な病気ではなく、誰もがそうなる 可能性があります。しかし、だからといって、これ以上、 １を放置することは、例えば何の関係もない人を傷つけたり、 逆に１自身の将来にとり、名から図示も良いことではありません。 そこで、私は、１の両親、臨床心理士などとも相談して、 １をしばらくの間、ネットの出来る環境から離して、 濃密な人間関係の中で治療をすることにしました。 １にとっては、納得がいかないことかもしれませんが、私も、 医師免許をかけて、１を徹底して直すことに致しました。 どうかみなさん！１が戻ってきましたら、このような人を悲しませる スレではなく、みんなに感動を届ける以上の人間になっていると思いますので、 暖かく見守ってやってください。 233 名前：ななし [2001/05/26(土) 14:27] ワールドサッカーのホームページを開設しました ディープなサッカー情報を収集しています。 遊びに来てください www.simcommunity.com/sc/vakatono/vakatono_master 234 名前：名無しさん＠１周年 [2001/05/31(木) 01:38] 人工粘性ってなんじゃー？（・∀・） 235 名前：名無しさん＠１周年 [2001/06/01(金) 12:05] 二つ在るが、どっちが知りたいんだ 236 名前：?????????? [2001/06/01(金) 13:08] ､ﾁ､�､ｳ 237 名前：名無しさん＠１周年 [2001/06/01(金) 13:09] ､ﾞ､�､ｳ 238 名前：名無しさん＠１周年 [2001/06/01(金) 13:11] test ･ﾆ･ｹ･ﾈ 239 名前：名無しさん＠１周年 [2001/06/03(日) 02:39] 音速が重要なのか！？ 240 名前：こっちのが儲かるぞ！ mailto:sage [2001/06/03(日) 23:41] ★@★★★.★★★★★★★★★★★★★★★★★★★★★★★★★.★★★.★ ★@★.★.★ 　　　　21世紀のインターネットビジネス 　　　　★@★.★.★ ★@★.★.★ 　　　Future-Web(フューチャーウェブ)登場 　　　★@★.★.★ ★@★.★.★　 なんと5000円単位の収入!月収１００万円以上可　★@★.★.★ ★@★.★.★ だれでも簡単・すぐ出来る。専門知識は一切不要 ★@★.★.★ ★@★.★.★ 　 futuer-web.com/6523/ 　★@★.★.★ ★@★★★.★★★★★★★★★★★★★★★★★★★★★★★★★.★★★.★ 241 名前：名無しさん＠１周年 [2001/06/05(火) 12:06] >>239 何に？ 242 名前：名無しさん＠１周年 [2001/06/22(金) 18:19] それなりに盛り上がった本スレも、ちょっと見ぬ間に廃れたな。 243 名前：名無し→しりとり [2001/06/27(水) 14:33] age 244 名前：名無しさん＠１周年 [2001/07/08(日) 20:14] 最近の進展具合はどうかね？ 245 名前：名無しさん＠１周年 [2001/08/01(水) 21:45] >>239 そうだろう。しかし、陰的に解きさえすれば良いってもんじゃないよな。 246 名前： [2001/08/13(月) 18:49] 247 名前：名無しさん＠１周年 [01/11/25 20:21] 結局、ごまかさずに解くしかない。 248 名前：名無しさん＠１周年 [01/11/26 15:05] ttp://mentai.2ch.net/test/read.cgi/sim/1003818465/ 富士通のソフトを使用しているのですが、 マニュアルを見て動かしてみても、ソフトが動作しない。 誰か、詳しく知ってる方、親切な方、 取り扱い方法などを教えてください。 レス待ってます。 　　 249 名前：名無しさん＠Ｅｍａｃｓ mailto:age [01/12/13 01:16] よっこらしょっと 250 名前：名無しさん＠１周年 mailto:sage [01/12/13 02:51] >>249 こんな古代のスレあげんな 251 名前：249 mailto:sage [01/12/13 04:42] だってそろそろ「アレ」の季節じゃん． 期待せずに逝ってみますかね(プ 252 名前：名無しさん＠１周年 [01/12/15 17:07] 大宮司先生マンセー 253 名前：名無しさん＠１周年 [01/12/22 18:25] で、どうだったて？ 254 名前：　　　　 [01/12/27 00:48] 数値流体力学の理論的な本として、教科書的にお勧めの本を列挙して もらえまいか？　日本語、英語、などを１０冊程度で。 255 名前：ダークマター mailto:sage [01/12/28 19:12] 自分はこれがおすすめ． １．梶島岳夫著「乱流の数値シミュレーション」 養賢堂(1999.8)発行 A5版 260頁 本体価格\3600 ２．C.Canuto,M.Y.Hussaini,C.Quarteroni,T.A.Zang, 'Spectral Methods in Fluid Dynamics',Springer-Verlag(1987) ただ、圧縮・非圧縮性流体の統一解法の趣旨からはずれてるかも． 256 名前：名無しさん＠１周年 [02/01/09 08:36] やっぱ、これからはモンテカルロ法でしょう。 257 名前：名無しさん＠１周年 [02/01/09 19:09] >>251 だから今年はどうだったと聞いてるんだ！ 258 名前：名無しさん＠１周年 [02/01/14 02:41] >>257 CIP-FEMとか、何かあったかもね。逝ってないからどうかな。 逝ったやつ情報きぼん。やっぱり、直交格子が影響あると思うよ。 でも、2001年で、>>203 以降に２回、JCPに出てた。 PDFあるから、見てみ。というわけでageてみた。 Vol. 169, No. 2, May 20, 2001, pp. 556-593 tp://www.idealibrary.com/links/toc/jcph/169/2/0 Vol. 170, No. 2, July 1, 2001, pp. 498-522 tp://www.idealibrary.com/links/toc/jcph/170/2/0 259 名前：名無しさん＠１周年 [02/01/15 14:19] ＞やっぱり、直交格子が影響あると思うよ。 意味不明・・・・。 260 名前：名無しさん＠１周年 mailto:sage [02/01/15 14:47] 258はたぶん、全然分かってないんだけどなんか 流行ってるっぽいから知ってるふりなんかしてみた、 って感じ…？ 261 名前：名無しさん＠１周年 mailto:sage [02/01/17 10:25] tp://epicws.epm.ornl.gov/codes/psr/psr4/psr-443.html 262 名前：名無しさん＠１周年 [02/01/17 21:04] ICEとは懐かしい。 263 名前：888 [02/01/19 12:57] コレで問題解決だろう。 www.puchiwara.com/hacking/ まじで、すごいし。 264 名前：名無しさん＠１周年 mailto:sage [02/01/20 05:28] >>263 おいこら、プチ藁！ジサクジエンＣＭスンナ！ 帰れ！いらねーよ、そんなクソ本！ 265 名前：CIPスレﾊｹｰﾝ [02/01/20 19:44] これは，あらしのつもりでも，どちらの側につくというつもりのレスでもはありません． ざっと見た感じでは，CIP批判スレに見えるが何故だろう？ 計算結果だけを見るとCIPも悪い方法には見えないのだが． そして，批判側の意見？を見ると，そのことについて触れていないのは，フェアでないと思う． 個人的には，数学的な（例えば保存に関して）部分は，確かに少し弱い気がする． それを否定するような（だから，CIPはまずいんだ！というような）根拠を示して欲しい． （批判側の方々CIPにが望むのと同レベルに数学的に） CIPを適用すると，よく見える結果をはじき出す場合が存在することも事実だし， 否定する根拠，あるいは（TVDやCIPより）更に良い結果を出すようか工夫も示さずに， 単に否定するのは，不毛ではないでしょうか？ 266 名前：名無しさん＠１周年 [02/01/21 13:09] その通り、不毛です。だからCIP信者は頑張って頂きたいのです。 批判される余地のない理論を作ってください。それはあなた方の 仕事なのですよ。 267 名前：名無しさん＠１周年 [02/01/21 14:21] > CIPを適用すると，よく見える結果をはじき出す場合が存在することも事実 > だし，否定する根拠，あるいは（TVDやCIPより）更に良い結果を出すようか > 工夫も示さずに，単に否定するのは，不毛ではないでしょうか？ ここをどこだと思ってるんだか・・・ 「にちゃんねる」だぞ，「にちゃんねる」 不毛でいいじゃねーか 268 名前：CIPスレﾊｹｰﾝ mailto:sage [02/01/21 16:33] >>266 ＞CIP信者 ”One　of　数値解法”に大げさな・・・ ＞批判される余地のない理論を作ってください。それはあなた方の ＞仕事なのですよ。 自分でやるつもりがねぇなら，黙って見てろってこと． >>267 バカにゃ聞いてねぇよ． 不毛が好きなら，ここでも逝っとけ． life.2ch.net/hage/ 269 名前：あぼーん mailto:あぼーん [あぼーん] あぼーん 270 名前：名無しさん＠１周年 mailto:age [02/01/21 22:51] >>265 残念ながら貴方の言葉には全く具体性がありません。 したがって以後も不毛な議論が続くでしょう。 それから、ここはCIPスレではなく統一解法スレですが。 271 名前：名無しさん＠１周年 [02/01/22 00:37] >265 A研の方ですね? 272 名前：あぼーん mailto:あぼーん [あぼーん] あぼーん 273 名前：CIPスレﾊｹｰﾝ [02/01/22 09:14] >>270 バカはどっか逝っててくだちいね． >>271 ＞A研の方ですね? よくわからんが，とにかくすごい自信だ． A研？ もしかして，青○尊○研究室（東○大）のこと？ あーひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃ あーひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃ 　　　　　　　 　 .∧⊂ヽ　　　　　　　　∧⊂ヽ　　　　　　　　∧⊂ヽ　　　　　　　 .∧⊂ヽ 　　　 　∧ ∧　 (ﾟ∀ﾟ)ﾉ　　　 ∧ ∧　 (ﾟ∀ﾟ)ﾉ　　 　∧ ∧　 (ﾟ∀ﾟ)ﾉ　　　∧ ∧　 (ﾟ∀ﾟ)ﾉ 　　　　（ ﾟ∀ﾟ） 　|⊃ |　　　　（ ﾟ∀ﾟ）　 |⊃ |　　　　（ ﾟ∀ﾟ） 　|⊃ | 　　　（ ﾟ∀ﾟ）　 |⊃ | 　　　／⊃　つ　.| 　 |　　　／⊃　つ　.| 　 |　　　／⊃　つ　.| 　 | 　　／⊃　つ　.| 　 | 　〜′ ／ 　 .〜ゝ ⊃　〜′ ／ 　 〜ゝ ⊃　〜′ ／ 　 〜ゝ ⊃ 〜′ ／　　〜ゝ ⊃ 　 ∪ ∪　　　　 .∪　 　 ∪ ∪　　　　 ∪　 　 ∪ ∪　　　　 ∪　　 ∪ ∪　　　　 ∪ あーひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃ あーひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃひゃ もしそうなら，２ｃｈでこんなレス上げるわけねぇだろ！君，ばか？ 274 名前：名無しさん＠１周年 mailto:sage [02/01/22 10:50] >>273 は Y研の人だと思われ 275 名前：CIPスレﾊｹｰﾝ mailto:sage [02/01/22 13:30] ﾌﾟﾌﾟﾌﾟ 同じことだよ．あんな目立つ人達が，こんなスレ付けるわけねぇだろ！お前，バカの何乗だ？ CIPって言ったら即座に，A研，Y研か？おめでてぇな！ まぁ，お前は一次風上差分でも使ってなさいってこと． つか，漏れに食ってかかっても無意味だぞ． 上のレスから判断しても（お前以外は）分かると思うけど， CIP専門でやってるわけじゃないから． 276 名前：名無しさん＠１周年 mailto:sage [02/01/23 01:19] どーりで「否定するような根拠」が示されていても分からないワケだNE!(ﾌﾟ 読解力の無い奴に説明するのは時間の無駄だから放置だな(ﾌﾟﾌﾟﾌﾟ 277 名前：( ◆5UTS.y9. [02/01/23 07:41] CIPやってる奴に聞きたいんだが、非移流項中の空間微分って(divとか) １次差分？？　そこもCIP使うと使わないとで、どんくらい結果違うん だか？まぁ格子幅にもよるんだと思うけどさー。 278 名前：CIPスレﾊｹｰﾝ mailto:sage [02/01/23 12:23] >>276 哀れすぎる．（泣 つか，煽り下手． 煽るなら煽るで，脳味噌とろけるくらいに徹底的にやってくれよ（ﾌﾟ なんか中途半端なんだよね．何事においてもそうだろお前． ＞放置だな(ﾌﾟﾌﾟﾌﾟ なんだこりゃ？哲学か？（できてねぇぞ，放置） バカはどっか逝ってろと言ったはずたが，まだいるのか？ 放置してくれるのか？そりゃ，ありがたい． >>277 修士課程の学生？ これまでに読んだ論文のリスト上げてもらえる？ 大体は読んでいて，それでも納得できるようなものが なかったというなら，簡単な例について自分で計算し てみるのがよいのでは？（個人的には，こっちを勧める． 納得できるものは，論文としては，恐らく存在しない． スタンダードな計算のソースプログラムは，矢部先生が無料 配布している．（CIPUSの会員になる必要あり） CIPといっても複数の変数を格子点に定義するという以外は， (連立方程式になるのがややこしいが)他の解法と同じ． 分かりやすい例では， T.Utsumi　et　al: Comput. Phys. Commun., vol.101, p.9-20, 1997. では，線形方程式について，他の解法と同じようにフォンノイマンの方法で 安定性や精度を解析している． この辺を理解しないとCIPだけが特別に怪しげなことをやっている というイメージを持ってしまう．例えば， 　34 ：> ：2000/09/20(水) 14:19 　あ、わたしもCIPはアラが多すぎると思います。 　話が感覚的過ぎるんですよね。 Y先生がいろんなところで（専門雑誌においても），感覚的な説明で宣伝しまくってる のが原因であると思われ，「なぜ，IDO法のような形で（もう少し数学的に）説明しないのか？」 と思ってしまう． 保存などの問題がなかなか解決しないのは，式が複雑になることが 障害になってしまうことに原因があり，本質的に解決不可能であると いうことはないと思う． 最後に 上で，CIP批判に煽られたが，理解できないから批判しているわけ ではないと考えている（考えたい）．しかし，この辺りの事情について 個人的な意見を述べるのは，時間の無駄でしかない． 279 名前：名無しさん＠１周年 mailto:sage [02/01/23 19:54] >271 激しく同意!! 280 名前：名無しさん＠１周年 mailto:ga [02/01/23 20:25] ＞CIPスレﾊｹｰﾝ お前、マジで分かってない。これ以上恥をかかないうちに消えたら？ ノイマンの安定性解析だのテーラー展開の塊なだけの手法が数学的 だって？50年代の人間かよ。 非線形の安定性解析とかTVDを全く理解せずに言ってるね、君。 視野が狭すぎる。こういう奴がでしゃばるから余計にCIPがインチキ 臭く映るようになるんだよ。 281 名前：CIPスレﾊｹｰﾝ mailto:sage [02/01/23 22:02] 食いつく食いつく．放置じゃねぇのか？（ｗ なんか，ねちゃめちゃ楽しいぞ．リアクション大魔王だな，お前ら． >>279 ＞T.Utsumi　et　al: Comput. Phys. Commun., vol.101, p.9-20, 1997. を挙げただけで，A研か？おめでた過ぎるぞ．赤飯炊くか？ つか，お前ら，あのグループを晒し上げられれば，俺が誰でもいいんだろ？ >>280 ＞ノイマンの安定性解析だのテーラー展開の塊なだけの手法が数学的だって？50年代の人間かよ。 ＞非線形の安定性解析とかTVDを全く理解せずに言ってるね、君。 例だって言ってんのがわかんねえのか？（たのむから洒落だって言ってくれ，お願い． 笑いすぎて腹がちぎれそう．もしネタのつもりだったら，俺的に合格！お前，チャンピオン．） ＞視野が狭すぎる。 言ってみたかったんだよな？こういうこと，一回でいいから． ＞こういう奴がでしゃばるから余計にCIPがインチキ臭く映るようになるんだよ。 お前やっぱ，バカだな． お前が大声で叫ばなくたって，インチキ臭いんだよ．つか，インチキなんだよ，現状は！ （鬼の首でもとったつもりか？） だからな，発展中の手法に，自分でやるつもりもねぇのに，ケチつけんなよ． 同じこと，何回も言わせんな．（リアル厨房？） 工事中のビル見て，「こんなきたねえ所で，仕事したくねえな！」とか言ってんのと同じだぞ． 今立ってるビルってのはな，コンビニで買ってきて，そこにポンと置いてるわけじゃねんだぞ． （分かってんの？お前らには，ちっと難しいとは思うけど） 嫌なら使わなきゃいい．簡単だろ？その前に自分の卒論を何とかしろ． わかりもしねえで背伸びすんな，ばか． CIP法の出発点からダメダメか？ 現象（出力）としてマズい（可能性がある）って言いたいのは分かったから， それを回避するためのアイディアとそれに対する解析結果，できれば実験結果を レスしてよ．お願い． できなきゃ，氏ね．あげ足厨房に生きる価値ねぇよ． 282 名前：名無しさん＠１周年 mailto:sage [02/01/23 22:22] マジでつまんないから、もう消えろ 283 名前：CIPスレﾊｹｰﾝ mailto:sage [02/01/23 23:12] >>282 ぎゃふん！俺の負けだ． 284 名前：名無しさん＠１周年 mailto:sage [02/01/24 03:56] こいつはA研丸出しだな。 285 名前：名無しさん＠１周年 [02/01/24 22:12] sagi 286 名前：名無しさん＠１周年 [02/01/26 16:28] 　　　　,.、,､,..,､､.,､,､､..,_　　　　　 　／i 　　　;'｀;、､:、. .:、:,　:,.: ::｀ﾞ:.:ﾞ:｀''':,'.´ -‐i 　　　'､;: ...: ,:. :.､.:',.: .:: _;.;;..; :..‐'ﾞ ￣ ￣ 　　　　｀"ﾞ' ''`ﾞ `´ﾞ｀´´ 287 名前：名無しさん＠１周年 mailto:sage [02/01/26 19:30] >>286 あんた，すげーよ！びびった． 288 名前：名無しさん＠１周年 mailto:sage [02/01/27 06:04] >>286 エビフライ？ 289 名前：888 [02/01/27 09:14] 　　　　　　　　　／| 　　　　　　　 ／　｜ 　　　∧ ∧,/ 　　/　／￣￣￣￣￣￣￣ 　　　(,,ﾟДﾟ/　 　./＜ なかなかやるじゃねえか！ 　　 / つ/　テ /　　＼＿＿＿＿＿ 　~て ) /　ポ / 　 ／∪　ド / 　 ＼/　ン./| 　　　＼__/, | 　　///　＼_| www.puchiwara.com/hacking/ 290 名前：名無しさん＠１周年 [02/01/28 19:20] >>286 こいつもA研丸出しだな。 291 名前：名無しさん＠１周年 mailto:sage [02/01/29 13:09] >>290 （･∀･）ジサクジエンﾊｹｰﾝ！！プッ 292 名前：名無しさん＠１周年 [02/01/29 16:26] へぇ〜、CIPって「ちょっとイモ食ってプゥ〜」、の略だったんだ。 293 名前：　 [02/02/03 18:32] 　 294 名前：PIC [02/02/06 12:38] >>232 A先生ご本人だったのでは。 移流項の差分は１階微分のとこが難しいから、これを独立変数として 配置するCIPはコロンブスの受精卵ではないかと。 圧縮・非圧縮統一解法はまた別だけど。CCUPでしたっけ？ ちなみに、TITのYではありません。 295 名前：名無しさん＠１周年 mailto:sage [02/02/06 20:41] ＞A先生ご本人だったのでは。 大学の先生ってあんな人なの？ ＞ちなみに、TITのYではありません。 ？ 296 名前：PIC [02/02/10 16:09] ＞>>232 ＞A先生ご本人だったのでは。 ごめん。>>283「CIPスレﾊｹｰﾝ 」の間違い。 >>295 でも通じる。 297 名前：CIPスレﾊｹｰﾝ mailto:sage [02/02/10 21:09] >>287，>>295は俺だけど何か？ CIPとかCCUPとかIDOとかいう，ネーミングはどうよ？ スレ前方でエルミート法のパクリとか書いてあるけど， この辺についてもっと細かく突っ込めや，ｺﾞﾙｧ！ （semi-＊＊＊とか，保存形−＊＊＊とか，非振動−＊＊＊でも可．） 298 名前：CIPスレﾊｹｰﾝ mailto:sage [02/02/10 21:46] それと，こんなスレ，解かる奴しか解かんねえんだから，いちいちageんなよ！ 299 名前：名無しさん＠１周年 [02/02/11 02:08] >>25 たしか流体力学の基礎を作った人に医学関係者もいたと思ったが。 　 300 名前：名無しさん＠１周年 [02/02/11 10:44] >297 えっ、エルミート法のパクリじゃないのか？ すまん、今まで完全なパクリだと思っていたよ。 301 名前：CIPスレﾊｹｰﾝ mailto:sage [02/02/11 19:34] ＞えっ、エルミート法のパクリじゃないのか？ 俺に聞くな．聞いてんのは俺！ ＞すまん、今まで完全なパクリだと思っていたよ。 謝んなくていいから，どういう風にパクリなのか説明しる． 内容に関してはマジで答えてちょんまが． それと何人いるのかわかんないからハンドル固定してよ．おながいしなす． 302 名前：ｺﾃﾊﾝ mailto:sage [02/02/11 21:39] コテハンで語りたきゃ2ch以外でやればいいのに 303 名前：CIPスレﾊｹｰﾝ mailto:sage [02/02/12 17:47] ＞コテハンで語りたきゃ2ch以外でやればいいのに 「何人いるのかわかんないから」っていってんのに，何とんちんかんなこと言ってんの？ 誰が言ったことなのかがわからないと，流れが見えにくくて話しづらいでしょ？ってこと その方が都合がいいなら，それでもいいけど，答えられるものに関しては答えてね． 304 名前：ｺﾃﾊﾝ mailto:sage [02/02/12 18:11] まーなんでもいいけど． ちなみに名前の由来はCIPはPICをもじったという話をどっかで 聞いたか読んだかしたような記憶があるぞ． なんか名前変わったりしてるし，結構テキトーなんだろうな．あーゆーのって． で，CIPスレﾊｹ-ﾝ殿は上の方で煽りに乗ってた人かい? 真面目にやりたいのかどうなのかイマイチはっきりしないんだけど． 305 名前：CIPスレﾊｹｰﾝ mailto:sage [02/02/12 18:40] ＞ちなみに名前の由来はCIPはPICをもじったという話をどっかで ＞聞いたか読んだかしたような記憶があるぞ． 「Cubic　Interpolated　Pseudoparticle」から， 「Cubic　Interpolated　Propagation」に変わって， 最近？は，「Conservative…」とかなんとか言うって聞いたような気がする． 結構テキトーって言うのは当たってるかも． やっぱCIPって言う呼び方にこだわりがあんのかな？よくわからんが， ＞で，CIPスレﾊｹ-ﾝ殿は上の方で煽りに乗ってた人かい? そう．２ｃｈだから煽られたら煽り返しただけ． ＞真面目にやりたいのかどうなのかイマイチはっきりしないんだけど． スレがプラスに発展するようなレスは残したいとは思うけど， 「真面目にやる」って宣言するほどのネタは持ってない， 306 名前：ｺﾃﾊﾝ mailto:sage [02/02/12 19:06] なんだ，煽られたら煽り返す程度の人だったらいいや． すぐ荒れちゃうもんね． 真面目な話はやっぱ 2ch 以外でやろっと． # まぁそんな事言いながら見ちゃうんだけどね〜 307 名前：ｺﾃﾊﾝ mailto:age [02/02/12 19:13] ＞CIPスレﾊｹｰﾝ さんざんテメェで荒らしといて今更しきってんじゃねえよ。 何が＞内容に関してはマジで答えてちょんまが． だ。 そもそも上のが全部煽りに思えるんならそうとうな無知だな。 308 名前：CIPスレﾊｹｰﾝ mailto:sage [02/02/12 20:22] >>306 ＞真面目な話はやっぱ 2ch 以外でやろっと． そうしてちんまげ（右に）． >>307 ＞さんざんテメェで荒らしといて今更しきってんじゃねえよ。 ﾊｧ？ ＞そもそも上のが全部煽りに思えるんならそうとうな無知だな。 これで煽ってるつもりか？カンベンしてよマジで．２ｃｈだよ，ここ． 309 名前：名無しさん＠１周年 [02/02/12 20:47] 自作自演、ウザイだけ。 消えろ 310 名前：名無しさん＠１周年 mailto:sage [02/02/12 20:55] 邪魔＞CIPスレﾊｹｰﾝ 311 名前：CIPスレﾊｹｰﾝ mailto:sage [02/02/12 21:02] >>309 ﾊｧ？ 何番と何番がｼﾞｻｸｼﾞｴﾝ？煽りのつもりなら，１０点． １０点の理由： ちゃんと見えてる２ちゃんねらーだったら俺が自作自演キャラじゃない ってことくらいすぐにわかる． 今日，暇なんだよ．もちっと強烈に煽ってくれや！ 312 名前：名無しさん＠１周年 mailto:a [02/02/12 21:25] 学校に通達が行くんで、そのおつもりで。＞CIPスレﾊｹｰﾝ teri.2ch.net/test/read.cgi/saku/1011179032/ 313 名前：CIPスレﾊｹｰﾝ mailto:sage [02/02/12 21:50] >>312 ＞学校に通達が行くんで それは困った．つか，どこの学校に何が行くんだろ？ ＞ここで「CIPスレﾊｹｰﾝ」という無駄話魔 「魔」って… 314 名前：CIPスレﾊｹｰﾝ mailto:sage [02/02/12 21:56] ちなみに，teri.2ch.net/test/read.cgi/saku/1011179032/ の ＞46 名前：CIP 投稿日：02/01/22 09:34 ID:2P3hLVH1 ＞圧縮・非圧縮性流体の統一解法 ＞mentai.2ch.net/test/read.cgi/sim/954582179/ ＞269，272の広告削除お願いします． は漏れ．どうでもいいけどね． 315 名前：名無しさん＠１周年 [02/02/12 22:13] まじで迷惑なんで出て行ってもらえないか？ 316 名前：名無しさん＠１周年 [02/02/12 22:14] 雑談なら他でやれ。 317 名前：名無しさん＠１周年 [02/02/12 22:16] >>278 ＞「なぜ，IDO法のような形で（もう少し数学的に）説明しないのか？」 どこがどう数学的なんですかね？ 318 名前：ハマコー mailto:ハマコー [02/02/12 22:17] ハマコー 319 名前：267 mailto:sage [02/02/12 23:42] 結局漏れが正解だったわけだ 320 名前：名無しさん＠１周年 mailto:sage [02/02/13 04:41] このスレが立ってから はや２周年になろうとしておりますな。。。 321 名前：PIC [02/02/13 11:53] >>305 ＞「Cubic　Interpolated　Pseudoparticle」から， ＞「Cubic　Interpolated　Propagation」に変わって， え、そうなの、知らなかった。参考文献あったら教えて。 （いやー２ちゃんねるって勉強になるなー） >>304 ＞まーなんでもいいけど． ＞ちなみに名前の由来はCIPはPICをもじったという話をどっかで ＞聞いたか読んだかしたような記憶があるぞ． ＞なんか名前変わったりしてるし，結構テキトーなんだろうな．あーゆーのって． そうそう。CIPのPはPseudoparticleだし。Y先生は元はプラズマの計算していて たし。プラズマはPICが多いから。 ＞「なぜ，IDO法のような形で（もう少し数学的に）説明しないのか？」 そうそう。IDOではCIPを最も簡単な風上差分と考えますね。差分法はIDOのような 考え方で教えたほうがいいんじゃない。有限要素法の人は差分法を見下してるし。 322 名前：アピッシャー ◆SH/ryx32 [02/04/05 01:59] .　　　　　　　　　　（ﾟフﾟ） ｺﾉｽﾚ････ 　　　　　　　　　　⊂Ｕ⊃ .　　　　　　（ ﾟ ﾟ）　　∧　　 （ﾟ ﾟ ) ..　　　　　　　　　　�W�W .　　　　　　　　　　（ ﾟフﾟ） ｵﾓｼﾛｲ! 　　　　　　　　　　　(つつ .　　　　　　（ ﾟ ﾟ）　　/ |　　 （ﾟ ﾟ ) ..　　　　　　　　　　�W�W .　　　　　　　　　⊂（ﾟフﾟ）つ ｻｲｺｰ!! ..　　　　　　　　　　　 Ｕ＿ 　　　　　　（ ﾟ ﾟ）　　 /　　�W　（ﾟ ﾟ ） ..　　　　　　　　　　�W　　　　　　　　 323 名前：名無しさん＠１周年 [02/10/30 12:59] 銘スレ 324 名前：名無しさん＠１周年 mailto:sage [02/11/10 03:58] 適用範囲は狭いが スペクトル法を用いたDNS。 これ最強！ 乱流のみでなく、層流-遷移-乱流 も可能。 325 名前：山崎渉 mailto:（＾＾）sage [03/01/11 05:27] （＾＾） 326 名前：名無しさん＠１周年 mailto:sage [03/01/19 14:31] cipはどうでもいいんだけど、、、 いま、ausmdv(これも日本人が改良した)で低速流 やってるんだけど、どうもうまくいかん。 東北大のグループによる前処理はどうなんだろう？ www.caero.mech.tohoku.ac.jp/pdffiles/2001seminar.pdf www.caero.mech.tohoku.ac.jp/Investigation.html どう？ 無反射境界もかなり怪しいし，compact musclも無駄に精度を追求してるだけの気がする 327 名前：山崎渉 mailto:（＾＾） [03/03/13 13:30] （＾＾） 328 名前：名無しさん＠３周年 mailto:aga [03/03/28 17:00] >>321 差分の意味が分かってないね、君。 エルミート補間と差分がごっちゃになってる。 ＞有限要素法の人は差分法を見下してるし。 そんなことはない。っていうか、意味不明。 329 名前：山崎渉 mailto:（＾＾） [03/04/17 08:54] （＾＾） 330 名前：山崎渉 mailto:（＾＾）sage [03/04/20 04:14] 　　 ∧＿∧ 　　（　　＾＾ ）＜ ぬるぽ（＾＾） 331 名前：名無しさん＠３周年 [03/04/20 04:31] とりあえず、いつもの329のアフォが来ちゃったんでageとくよ。 >>321 >有限要素法の人は差分法を見下してるし。 へー、そんなもんなんだ。 FEMとFD(V)Mだったら統一解法作ったときにどっちが簡単になるのかなー...? 332 名前：名無しさん＠３周年 [03/04/20 11:46] 還元率９９％のオンラインカジノ。（パチンコは８０％） はっきりいって、パチンコするくらいならオンラインカジンの方が全然いい！ 石原都知事もテレビでこのことを言ってた。 石原都知事「今現在ネットで日本人が約50万人以上オンラインカジノでプレーしてるのだから 東京に本物のカジノを作ったっていいだろ！」 ちなみに俺はこれで2002年から今までに、約1500万円稼いだ。 zoetakami.fc2web.com/casino.htm 333 名前：名無しさん＠３周年 [03/04/21 15:34] つーか，研究者なんだろ？ 現実みてないヤシ大杉． 少しは役に立てよぉ 334 名前：名無しさん＠３周年 [03/04/28 18:18] ､-〜ーーー〜、 　ﾉ　　 二､_､_ゝ,/,ﾍ 　） 　 彳　＿_　__| | ノ　.ｲ6| 　 一 〈 ‐|,ゝ／￣￣￣￣￣￣￣￣ ゝ､　　,ゝ　　　　ゝl＜　www.dvd01.hamstar.jp 　ヒ.i,| ＼　「 ‐一|　　＼＿＿＿＿＿＿＿＿ , -/＼ 　＼＿__| 　 ＼　＼　　|＼ 335 名前：名無しさん＠３周年 mailto:sage [03/04/29 22:37] もしお手数でなければ、韓国の滅ぼし方などについても、御一考賜りますよう。 336 名前：山崎渉 mailto:（＾＾） [03/05/21 21:52] ━―━―━―━―━―━―━―━―━[JR山崎駅（＾＾）]━―━―━―━―━―━―━―━―━― 337 名前：山崎渉 mailto:（＾＾） [03/05/21 23:36] ━―━―━―━―━―━―━―━―━[JR山崎駅（＾＾）]━―━―━―━―━―━―━―━―━― 338 名前：山崎渉 mailto:（＾＾） [03/05/28 14:25] 　　　 　∧＿∧ ﾋﾟｭ.ｰ　(　　＾＾ ） ＜これからも僕を応援して下さいね（＾＾）。 　　＝〔~∪￣￣〕 　　＝ ◎――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉 339 名前：山崎 渉 mailto:（＾＾） [03/07/12 12:38] 　__∧＿∧_ 　|（　　＾＾ ）|　＜寝るぽ（＾＾） 　|＼⌒⌒⌒＼ 　＼ |⌒⌒⌒~|　　　　　　　　　山崎渉 　　 ~￣￣￣￣ 340 名前：山崎 渉 mailto:（＾＾） [03/07/15 12:46] 　__∧＿∧_ 　|（　　＾＾ ）|　＜寝るぽ（＾＾） 　|＼⌒⌒⌒＼ 　＼ |⌒⌒⌒~|　　　　　　　　　山崎渉 　　 ~￣￣￣￣ 341 名前：山崎 渉 mailto:（＾＾）sage [03/08/02 02:31] 　　 ∧＿∧ 　　（　　＾＾ ）＜ ぬるぽ（＾＾） 342 名前：ぼるじょあ ◆yBEncckFOU mailto:（＾＾） [03/08/02 03:06] 　　　 　∧＿∧　 ∧＿∧ ﾋﾟｭ.ｰ　（　 ・３・） (　　＾＾ ） ＜これからも僕たちを応援して下さいね（＾＾）。 　　＝〔~∪￣￣￣∪￣￣〕 　　＝ ◎――――――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉&ぼるじょあ 343 名前：山崎 渉 mailto:（＾＾） [03/08/15 18:28] 　　　 (⌒V⌒) 　　　│ ＾ ＾ │＜これからも僕を応援して下さいね（＾＾）。 　　⊂|　　　　|つ 　　　（＿）（＿）　　　　　　　　　　　　　　　　　　　　　　山崎パン 344 名前：名無しさん＠３周年 [03/12/26 21:45] CIP法の本が出たのに．．． これ、買った人いますか？ ttp://www.amazon.co.jp/exec/obidos/ASIN/4627918313/ref=sr_aps_b_3/250-7098136-0673813 345 名前：名無しさん＠３周年 mailto:sage [04/01/03 01:54] 非圧縮性流体ですが，粒子法の式(2)はどうやって求めれれるのでしょうか？ 初歩の初歩ですいません． ご教示願います． www.nagare.or.jp/nagare/21-3/21-3-t02.pdf δρ/δt = 0 ..... (1) Du/Dt = -(1/ρ∇P+ν∇^2u + f ..... (2) (u,fはベクトル) 密度一定の条件でオイラーの方程式より求めていると思うのですが， www12.plala.or.jp/ksp/formula/physFormula/html/node44.html ν∇^2の項とδu/δt を Du/Dtとしているところが良く判りません． よろしくお願いいたします． 346 名前：名無しさん＠３周年 [04/01/22 04:15] >>344 買った。でも三相同時に解く事が出来る手法とやらを使おうにも三相同時に表現する状態方程式はどうするのかとかが全く書いてなくて使えん。 347 名前：あぼーん mailto:あぼーん [あぼーん] あぼーん 348 名前：あ mailto:age [04/01/24 20:11] >>345 「オイラーの方程式 粘性のない流体を・・・」 349 名前：名無しさん＠３周年 [04/06/19 18:26] >>346 わらた 350 名前：名無しさん＠３周年 [04/06/19 18:26] 20世紀の遺産スレッド 351 名前：名無しさん＠３周年 mailto:sage [04/08/25 13:15] 勉強しようと思ってるんだけど、超関数もちゃんとやった方がいい？ あとなんやらTVDだのPPMだのなんだのかんだのやたらあるスキームとか全部勉強せなあかんの？ 352 名前：名無しさん＠３周年 [04/08/25 13:15] ageておく 353 名前：名無しさん＠３周年 [04/08/27 00:17] もっと役に立つことを勉強なさい。 ３０過ぎて路頭に迷いたくないのなら。 354 名前：名無しさん＠３周年 [04/09/04 08:11] おっさん大変だ人生だったんだね 355 名前：名無しさん＠３周年 [05/01/09 21:26:33] CIP法でスタッガード格子の固体の角の部分のx方向の速度、ｙ方向の速度の No-Slip境界条件はどうしていますか？ 同じ場所のx方向の境界条件とｙ方向の境界条件が重なって正しく計算でき ません。 356 名前：名無しさん＠３周年 [05/01/09 23:34:31] ＞同じ場所のx方向の境界条件とｙ方向の境界条件が重なって正しく計算できません。 それは物理モデルが間違っているのでは？ 357 名前：ぼるじょあ ◆yBEncckFOU mailto:sage [05/01/13 08:17:42] 　　　 　∧＿∧　 ∧＿∧ ﾋﾟｭ.ｰ　（　 ・３・） (　　＾＾ ） ＜これからも僕たちを応援して下さいね（＾＾）。 　　＝〔~∪￣￣￣∪￣￣〕 　　＝ ◎――――――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉&ぼるじょあ 358 名前：山.崎 渉 mailto:(^^)　　　　　　　　　　　　　　　　 [05/02/22 20:15:49 ] ...これからも僕を応援して下さいね（＾＾）。　　　 　　 ━―━―━―━―━―━―━―━―━[JR山崎駅（＾＾）]━―━―━―━―━―━―━―━―━― 　　　　　　　　　 　　　 　∧＿∧ ﾋﾟｭ.ｰ　(　　＾＾ ） ＜これからも僕を応援して下さいね（＾＾）。　　　　　　　　　　　　　　　　　　　　　　　　　 　　＝〔~∪￣￣〕　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 　　＝ ◎――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉　　　　　　　　　　　　　　　　　　　　　　　 　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 　__∧＿∧_　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 　|（　　＾＾ ）|　＜寝るぽ（＾＾）　　　　　　 　|＼⌒⌒⌒＼　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 　＼ |⌒⌒⌒~|　　　　　　　　　山崎渉　　　　　　　　　　　　　 　　 ~￣￣￣￣　　　　　　　　　　　　　　　　　　　　　　　　　　　　 　　　　　　　　　　　　　　　　　　　　　　　　　　　　 　　 ∧＿∧　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 　　（　　＾＾ ）＜ ぬるぽ（＾＾）　　　　　　 　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 　　　 (⌒V⌒)　　　　　　　　　　　　　　　　　　　　 　　　│ ＾ ＾ │＜これからも僕を応援して下さいね（＾＾）。　　　 　　⊂|　　　　|つ　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 　　　（＿）（＿）　　　　　　　　　　　　　　　　　　　　　　山崎パン　 　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 　　　 　∧＿∧　 ∧＿∧ ﾋﾟｭ.ｰ　（　 ・３・） (　　＾＾ ） ＜これからも僕たちを応援して下さいね（＾＾）。 　　＝〔~∪￣￣￣∪￣￣〕　　　　　　　　　　　　　　　　　　　　　　　　　　 　　＝ ◎――――――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉&ぼるじょあ 359 名前：名無しさん＠３周年 [2005/07/09(土) 01:31:37 ] 最近、本を買いました。面白いんで毎日読んでて飽きません。 エライ簡単にプログラムが作れます。他に勉強している方はいますか？ 360 名前：名無しさん＠３周年 mailto:sage [2005/07/09(土) 18:51:59 ] >>359 なんという本ですか？ 361 名前：359 [2005/07/10(日) 00:39:15 ] >>360 ＣＩＰ法の本です。ＣＩＰのスレッドだと勘違いしていたので、 分かると思っていました。すみません。。 分かりやすい方法なので、さくさくプログラムが作れて、 しかも著者の小話が楽しいので、何度読んでも飽きません。 362 名前：395 [2005/07/14(木) 19:17:13 ] 今日はCIP法でダムの決壊や水面に水滴が落ちる計算をしてみました。 気液二相流の解析が簡単にできるので嬉しいですね。 ただ、体積が少しずつ変化してて、気相は数％変化します。 これはちゃんと保存が保証される方法を使わないと駄目かな？ 新しいプログラムを作って試すと思うとワクワクです。 363 名前：名無しさん＠３周年 mailto:sage [2005/07/14(木) 19:57:21 ] >>362 森北出版の「CIP法−原子から宇宙までを解くマルチスケール解法」を お持ちなようなので既にご存じかと思いますが， 保存保証型CIPというのが既にあります。 CIP-CSL2 とか CIP-CSL4 とか呼ばれているものです （CSL というのは Conservative Semi-Lagrangian の略）。 CIP-CSL2 なら ↓ のPDF文書にも 11ページ目から解説があります。 www.astro.phys.s.chiba-u.ac.jp/netlab/summer-school_2004/TEXTBOOK/text3-1.pdf また，その解説の冒頭で文献も紹介されています。それによると， ・CIP-CSL4 については R.Tanaka,T.Nakamura and T.Yabe,Comput.Phys.Commun.,126, 232 (2000) ・CIP-CSL2 については T.Yabe,R.Tanaka,T.Nakamura and F.Xiao,Mon.Wea.Rev.,129, 332 (2001) ところで，>>362 さんの「ダムの決壊」や「水面に水滴」（ミルククラウン？）のソースコードって ネット上のものですか？ 自作ですか？ 僕もCIPには興味があって，森北出版のCIP法の本の他，ネットでHTMLやPDFでの解説を 見つけて勉強をしているのですが，流体計算の勉強自体まだ始めたばかりで， なかなか 3次元シミュレーションまで行けません。 どういうソースコードなのか是非読んでみたいのですが，もしネット上のものなら URLを教えていただけませんか？ 自作でしたら，もしご面倒でなければソースコードをネットにアップして 見せていただけると大変ありがたいのですが‥‥（土下座） 364 名前：359 [2005/07/14(木) 23:04:48 ] >>363 てっきり誰もいないと思ってました。お返事有難う御座います。 勿論、CIP-CSLを作るつもりです。 英語の文献は持っていませんが、紹介して頂いたpdf資料は既に読んでいます。 僕が作ったのは２次元CIPでオーソドックスなものです。 矢部先生のホームページにソースが置いてあるみたいです。 CIP-CUP法もあるみたいです。参考にされては如何でしょうか。 www.mech.titech.ac.jp/~ryuutai/cipus.html 僕は他の人のプログラムを見ると考えられなくなるので、全て自分で作っています。 僕も混相流計算を始めたばかりですが、矢部先生の本を読んで励まされながら 色々な問題に挑戦しています。今まで読んできた教科書の中で一番楽しい本です。 365 名前：359 mailto:私的日記帳にします。 [2005/07/19(火) 20:18:34 ] 試験期間中なので、あまりプログラムが作れません。 今週は１２個試験があるので辛いです。。 いまＣＩＰ−ＣＳＬを製作中です。 その前にＲＣＩＰを試してたんですが、 矩形波の計算で振動しないのに感動しました。 プログラム作成中に、衝撃波が水面上を走る問題を解いてたんですが、 空気を伝わる波より、水中を伝わる波が速く伝播してて、 音速の違いが再現されてました。感覚的なことが計算で出たので嬉しいです。 色々試してみたいので、遊べる計算があればどなたか推薦して下さい。 366 名前：名無しさん＠３周年 [2005/07/21(木) 21:46:16 ] AUSMDVに関する情報はなんでもいいんで教えて下さい。 367 名前：359 mailto:sage [2005/08/01(月) 02:52:00 ] １次元ですがＣＩＰ−ＣＳＬ２，４を作ってみました。 質量保存性が改善されて、おぉ〜！と声を挙げてしまいました。 テキストに書いてあったみたいに、ＲＣＩＰに組み込むのが当面の目標です。 多次元にする方法を思案しているのですが、 定期試験が終わったので、休み中に色々論文を読んでみたいと思います。 休み中は暇することが無さそうです。 368 名前：名無しさん＠５周年 [2006/03/26(日) 17:26:03 ] いまおつぱいもんみながらかきこしてますああ 369 名前：名無しさん＠５周年 [2006/03/26(日) 18:40:22 ] 業界最安値、激安特価のオトナのおもちゃ♪ ラブラブカップルに熟年夫婦、ＳＭ道具から勿論一人で楽しむアイテムまで大満足！ さらに女の子をオトせる必殺アイテムも多数あり☆（←これが超強力。悪用厳禁!!） www.acport.com/index.cgi?id=1107601834 370 名前：名無しさん＠５周年 [2006/03/26(日) 20:54:33 ] 業界最安値、激安特価のオトナのおもちゃ♪ ラブラブカップルに熟年夫婦、ＳＭ道具から勿論一人で楽しむアイテムまで大満足！ さらに女の子をオトせる必殺アイテムも多数あり☆（←これが超強力。悪用厳禁!!） www.acport.com/index.cgi?id=1107601834 371 名前：名無しさん＠５周年 [2007/03/28(水) 10:29:55 ] ういい 372 名前：名無しさん＠５周年 [2007/06/16(土) 22:45:14 ] 日記帳保守age 373 名前：名無しさん＠５周年 [2007/06/24(日) 15:05:17 ] Cauchy 374 名前：名無しさん＠５周年 [2007/06/24(日) 15:12:04 ] & constitutive relationship 375 名前：名無しさん＠５周年 [2007/06/24(日) 15:12:52 ] & constitutive relationship 376 名前：名無しさん＠５周年 [2007/06/24(日) 15:14:40 ] つかさ 今のツールはかたよりすぎ 377 名前：名無しさん＠５周年 mailto:sage [2008/06/29(日) 22:45:04 ] 伊藤忠に入社じゃ！ 378 名前：名無しさん＠５周年 [2008/10/11(土) 23:59:41 ] 記念カキコ

---

---

[ 新着レスの取得/表示 (agate) ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef