- ベストアンサー
差分法による電位の計算の精度について
ksugaharの回答
詳細な情報ありがとうございます。 境界条件が怪しいように思います。 ちなみに、点電荷が空間にぽつりとある状態は、 電磁界解析を行うのは、かなり大変です。 接地導体があると精度が出やすいのですが 無限遠がGNDである場合には精度がでません。
関連するQ&A
- SOR法による差分方程式の解法
プログラミングのカテゴリーか少し悩んだのですが、 物理関係で差分方程式を解いている人が多いかと思い、 こちらで質問させて頂きます。 解きたいのは、２Ｄポアソン方程式です。(1Ｄでも構いません) 微分方程式を差分化してSOR法を使用します。 SOR法は若干ながらメッシュサイズ、加速係数に 解（ポテンシャル形状）が依存します。 しかし、できるだけこういった条件に依存せずに、 正確な解を得たいです。 より安定に正しい（メッシュサイズなどに依存しない）解が得られる 方法など、SOR法について何かご存知の方、何でも構いませんので、 ご教授のほど宜しくお願いします。 又、複雑な構造（ダブルHEMT構造、局所的に高ドープなど）でも、計算を 収束し易くする方法など、何かありましたらアドバイスをお願いします。
- ベストアンサー
- 物理学
- 数値解析の手法（差分法）について
現在、とある2元の1階偏微分方程式（解はu,vでそれぞれ右と左に進行する波）を数値解析によって解こうと考えています。 数値解析の手段として、差分法がよく用いられると思いますが、 現在、私は、場所に関してはuは後退差分、vは前進差分を使い、 時間に関しては前進差分を使って解いています。 ネットでは場所に関しては中心差分、時間に関してはルンゲ=クッタやリープフロッグなどが 使われていることが多く、私もこの２つを用いて解いてみました。 偏微分方程式には線形項が含まれていたため、 線形問題に対する制約であるΔt/Δx << 1は最低満たすように刻み幅をいろいろ取り、 計算時間も辞さず計算機を動かしてみましたが、 ノイズが消えず、解析解に限りなく近づくには至りませんでした。。 Δt/Δx=0.0001なども試したのですが・・・ そこで、諦めて違う差分法を試し、 場所に関して、uは後退差分、vは前進差分を使い、 時間に関しては前進差分を使って見たところ、 Δt/Δx=0.01程度で解析解に近い、なかなか精度の良い数値解を得ることが出来ました。 2次の差分では上手くいかず、1次の差分だとわりかし上手くいく・・・ 精度的には中心差分やルンゲ=クッタなどの方がいいと思うのですが・・ 正直不思議でなりませんでした。。。 最初に試した差分法のコードミスかと思い、何回もコードを確認し直しましたが、 やはり解析解に近づくには至りませんでした。 こんなことってあるのでしょうか？？ 差分法でも場合によって使い分ける必要があるということでしょうか・・？ その場合分けするときの指標など、知っておられる方、教えて頂けると助かります。 問題によって
- 締切済み
- 数学・算数
- 点電荷が作る電位分布の求め方
点電荷が作る電位分布の求め方 お世話になります。 1次元の電位分布についての質問です。 高校の物理で習ったように、1[C]の点電荷（あるいは微小な大きさを持つ電荷）が原点にあるときの電位分布は、無限遠をゼロとして、 φ = (1/(4πε)) * (Q / r) ・・・(1) で表せますよね？ 同じ分布をポアソン方程式（div (grad φ) = -ρ/ε）から求めるにはどうすればよいでしょうか。 1次元の場合ポアソン方程式は単純な2階微分方程式になると思いますので、rで2階積分してみたのですが、原点以外ではρ= 0 のため φ が一次関数になってしまい、(1)のような反比例の関係にはなりそうにありません。 どこか考え方が間違ってるのだと思いますのでご指摘いただけると助かります。 よろしくお願い致します。
- ベストアンサー
- 物理学
- ポアソンの方程式を用いた電位の求め方
ポアソンの方程式を用いて原点に置かれた単位点電荷による電位を求めるという問題なのですが結果がガウスの法則によって求まる値と異なってしまいました。 一つ前の問題で電位がC/r(Cは定数)と表せることが求められたので、 ポアソンの方程式より ∇^2(C/r)=-ρ(r)/ε 両辺を半径rの円で体積積分すると、ガウスの発散定理より 左辺=∬∇(C/r)・dS =4πr^2/r^2×{∂(cr^2/r)/∂r} (極座標における発散より) =4πc 右辺=-1/ε∬∫ρ(r)dv =-1/ε (原点に単位点電荷があるのみなので) 右辺と左辺を比較することでc=-1/4πεとなり、電位は-1/4πεrとなるのですが、 ガウスの法則によって求めた電場を積分して求められる電位は1/4πεrとなります。（こちらは教科書に載っていたのであっていると思います） どこが間違っているのでしょうか。 基本的な問題ですがご存知の方、よろしくお願いいたします。
- ベストアンサー
- 物理学
- 熱伝導方程式を差分法で解くプログラム
前回熱伝導方程式の差分法の解き方について質問させて頂いたものです。 https://sp.okwave.jp/qa/q9363032.html 回答ありがとうございました。おかげさまで差分法で解くプログラムを書くことができました。 書けたのは良いのですが、間違っている可能性が大きいです。 私の考えでは、安定性条件を満たしていれば dt,dxを小さくすればするほど、収束温度は変わらず計算精度が高くなっていくと思ったのですが、 私の書いたプログラムではdxdtの値を変えると収束温度まで変わってしまいます。 プログラムに関しても初心者なので説明が難しいのですが、やはり収束温度が値によって変わるということはプログラムがどこか間違っているということでしょうか。 度々すみませんが、回答いただけるとありがたいです。よろしくお願いいたします。
- 締切済み
- 物理学
- ３次元の差分法について
よろしくお願いします。 現在、３次元の差分法について勉強しています。プログラムにて電位や電界の計算を行っているのですが、本やサイトでは簡単にするために２次元で扱っているものが多いです。 何かお勧めの参考書やサイトがあれば教えていただけますでしょうか。 よろしくお願いします。
- ベストアンサー
- 物理学
- 単精度計算と倍精度計算
Ｆ（ｘ）＝（１－ｘ・ｘ）１／２乗において、ｘの値を０．９、０．９９、…０．９９９９９と変化させた場合の単精度計算と倍精度計算の結果を求めるプログラムを教えてください。
- ベストアンサー
- C・C++・C#
- H2のエネルギー計算
H2のSchrodinger方程式は変分法を用いて，高精度の解が求められると思います．現在では100桁を超える高精度の解が求まっていると聞きました． ここで疑問に思ったことがあります．その100桁をも超える高精度の解といっても，本当に正しい値というのはどのようにして求まっているのでしょうか．実験でも誤差がでると思いますし，そもそも不確定性原理から100桁というのは正しい値といえるのでしょうか？
- ベストアンサー
- 化学
補足
返信が遅れて申し訳ありません。 回答していただきありがとうございます。 私も考えた結果、境界条件がよくないと感じました。 よって自由空間ではなく、設置されて0Vの４辺に囲まれた正方形内という条件で考えてみたいと思います。 電位の式が変わるので電気影像法で求めてみたいと思います。 また何かありましたらご教授よろしくお願いします。