粒子法(MPS法)について MPS法コードを作成しています。射流現象を解いているのですが、うまく行きません。圧力による斥力が強く働いてしまい、1つの粒子が計算領域全体に悪影響を及ぼし、一気に粒子が飛散してしまいます。MPS法の圧力振動の抑制についての論文通りに、1つ前の時間ステップの粒子数密度を用いて圧力を計算しなおす方法も試しましたが、やはりダメでした。 MPS法コードを作成されたことのある方で、コードの工夫した点や、成功・失敗談をよろしければお聞かせください。よろしくお願いします。
![粒子法(MPS法)について - MPS法コードを作成しています。射流現象を解いているのですが、うまく行きません。圧力による斥... - Yahoo!知恵袋](https://cdn-ak-scissors.b.st-hatena.com/image/square/1d07bee2b75b182ba712690f3a3464c29972e28b/height=288;version=1;width=512/https%3A%2F%2Fs.yimg.jp%2Fimages%2Fks%2Fclap%2Fimage%2Fogp%2Fogp.png)
Next: Contents 円周率の公式集 暫定版 Ver. 3.141 編集 : 松元隆二 December 26, 2000 暫定版です.書きかけの章があります. この公式集は,LATEXで書かれた原稿をLATEX2htmlで自動変換した ものです. 自動変換が正しく行われていない部分があり,一部の文章が欠けたり,式番号, 章番号が間違ってる部分があります.Adobe Acrobat版,ASCII pTEXDVI版で は正しく表示されます. Adobe Acrobat版,ASCII pTEXDVI版は,ここにあります. 数式は画像ファイルとして表示されます.したがって画像ファイル が いっぱいあります.ネットワークが細い人は見るのは止めたほうがよいでしょ う. Contents 独自定義の数学記号 円周率の公式 定義 微分積分学以前 Archimedes of Syracuse (
カシオの高精度計算サイトを試してみた。適当に「統計関数 > 超幾何分布」とたどって,例題の x = 14,n = 50,M = 60,N = 200 をそのまま試してみる。演算桁数は50桁まで選べるが,とりあえずデフォルトの演算桁数(22桁)では 確率密度 f(x,n,M,N) = 0.1342936640436631826 下側累積確率 P(x,n,M,N) = 0.4339792084682835697 と19桁まで出る。Rで検算(にならないか ^^;)してみたら,CPUによって結果が違う: > options(digits=16) > dhyper(14,60,140,50) [1] 0.1342936640436632 # Linux (Intel) [1] 0.1342936640436631 # Mac (PPC) > phyper(14,60,140,50) [1] 0.43
概要 CG法(Conjugate Gradient Methods)はM.R.HestenesとE.Stiefelによって1952年に提案された方法である[1]。 CG法は正定値対称行列に対して使われる連立一次方程式を反復法で解くための手法である。 行列の正値対称性 ベクトルの内積をのように書く。 実行列が正定値対称とは、 ということであり、が対称であるということは、 が成り立つということである。 CG法の基本原理 今、次のような線形同次方程式を解くとする。 CG法は回目の反復において、次のようにこの方程式の解や誤差を用いて定義される誤差のノルム (等号成立はのとき) を最小化するような近似解を部分空間の中から見つける方法である。但し、はクリロフ部分空間(Krylov Subspace)である。 つまりCG法は次のような連立一次方程式の近似解を探すための方法である。 このように部分空間の中
リリース、障害情報などのサービスのお知らせ
最新の人気エントリーの配信
処理を実行中です
j次のブックマーク
k前のブックマーク
lあとで読む
eコメント一覧を開く
oページを開く