現在「数値流体工学」著:荒川忠一 の本を参考に,差分法を用いた流体解析のプログラムを作成をしました.
そして,解析をしたところ計算の途中で差分値が0になり,止まってしまいます.
原因を探してみたのですが,予想としましては,
1.モデルの格子幅の大きさに依存(幅が大きいと早く止まる)
2.収束率を速めるために用いている緩和係数の大きさに依存
と考えましたが,実際に何をどうすれば解決するのか
という所までは考えが思い浮かばず,見当が付きません.
もし,このような経験をした方,もしくはこの本をご存知の方がいらっしゃいましたら,なにかアドバイスをいただきたいのです.
うまく問題事項が説明できませんので,ご質問とかありましたらお願いします.
大変困っております.
ぜひ,小さな情報でも構いません.よろしくお願いします.
No.2ベストアンサー
- 回答日時:
りゆうたいは計算したことがありません。
一般的数値計算の差ぶんとして答えます。計算機特有の問題.桁落ちしていませんか。
計算方法事体に間違いがないとして.
初期値に極端に変な値を使用している悪い問題の場合
が考えられます。
鞍部を挟んでさぶんしているなんて.やっていませんよね。解の巣に引っ掛かると.平気でこれやってしまいますから。
たんせいどと倍精度で答えが大幅に異なるのであれば.計算機特有の問題です。
ライブラリーがあれば4倍精度で計算してみてください。なければ.極端な場合には.自分で必要な桁の演算ルーチンを作ることになります。
帯域ライブラリーを使っている場合に.注意点としては.定義域とちいきの問題です。よく.範囲外の数値を入れで.0が帰ってきて(通常エラーの時に0を返す)そのまま使っていて.差分が0になった.なんて良くある話です。三角関数でよく発生します(使っているかどうかわかりませんけど)。
No.1
- 回答日時:
まずは、簡単な例(かつ、正解のわかっている問題)
で試されていると思いますが、
それが、正しい結果にならないわけですか?
今一度、本とプログラムを見比べ、バグをさがされては
いかがですか?
私も何度も教科書を疑ったことがありますが、
めったに教科書は間違っていないものです。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- その他(自然科学) 科学技術計算の仕事について 2 2023/02/04 18:09
- 統計学 Rを用いた「繰り返しがある直交表実験計画法」の分析方法 8 2023/08/01 17:58
- 大学・短大 女子大生(成人している大学2年以降の学生)の飲酒とリスクについて質問です。 5 2022/03/26 19:45
- 医学 女性の飲酒の是非について 2 2022/03/26 16:40
- 宇宙科学・天文学・天気 AIが答えた方程式 1 2023/02/20 00:12
- 生活習慣・嗜好品 お酒が強い体質の人はアルコール依存症になりやすいという話は事実ですか? 5 2022/03/26 18:35
- 物理学 時間を語るなら、(複数の時間の正体)を知る必要が有る。 1 2023/02/16 22:14
- 数学 原始関数の存在性の証明について 数学科の3回生です。院試の勉強でつまづいたので助けてほしいです。 R 6 2022/11/13 19:19
- その他(ソフトウェア) F-BASICで計算中の実行が中途で勝手に止まり、大変困っています。 2 2023/03/02 16:15
- 妊娠 出産経験のない女性は飲酒を全面禁止すべきではないでしょうか? 3 2022/04/12 17:27
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
65536は2の何乗なのでしょうか?
-
バッチファイルでウインドウを...
-
変化させるセルが変化しない
-
排他的論理和 BCC(水平パリテ...
-
骨折リスク評価のFRAXについて...
-
傾いた四角形内の範囲の条件式
-
タクシー料金の問題です
-
C++言語で、レナードジョーン...
-
計算式の文字列を解析して自分...
-
EXCELなどで「返す」という表現
-
べき乗の計算が遅い理由
-
数値計算の高速化 (cos, sin, exp)
-
ExcelのマクロVBA動作を軽くす...
-
prolog計算
-
CCurの引値を文字列にする時と...
-
VBAの再計算が反映されない件に...
-
Excel VBA TREND関数について
-
順列のプログラムについて(VB)
-
C言語初心者。静磁場の計算。台...
-
エクセルで特定のセルのみを任...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
65536は2の何乗なのでしょうか?
-
VBAの再計算が反映されない件に...
-
排他的論理和 BCC(水平パリテ...
-
EXCELなどで「返す」という表現
-
バッチファイルでウインドウを...
-
モジュラス103の計算とは何でし...
-
傾いた四角形内の範囲の条件式
-
Visual C++でdebugとreleaseで...
-
変化させるセルが変化しない
-
骨折リスク評価のFRAXについて...
-
C# 計算処理中に実行中ウィン...
-
VBAでの勤務時間計算
-
べき乗の計算が遅い理由
-
C言語についてです。 再帰を使...
-
Excel VBAにてFFT
-
数値計算の高速化 (cos, sin, exp)
-
VBとVBAの違い
-
VB6で正確なミリ秒を計測したい...
-
スレッド処理からダイアログを...
-
matlabで計算終了
おすすめ情報