![](http://oshiete.xgoo.jp/images/v2/pc/qa/question_title.png?5a7ff87)
直線電流周辺での荷電粒子の運動
(0,0,z)方向に一様に電流Iが流れているとします。ビオサバールの法則からその周辺にはB = μI/2πr (r=√{x^2+y^2})の磁場が発生するはずです。
ここで、ある座標(r_0,0,0)から速度(0,v_0,0)を持たせて粒子を打ち出した場合の軌跡は計算可能でしょうか。自分はma = q(v×B)を書き下し、dz/dt = k{ln(r/r_0)} , k≡qμI/2mπ 以上計算が出来ませんでした…。
数値計算でシミュレーションして見たところ、z方向にはほぼ等速で移動しつつ、電流周りに周期的に回転しながらr方向に振動する粒子の運動が検出されました。
数値計算の結果の妥当性を調べたいのですが、上の式は解析的に解けるのでしょうか。。そしてそもそもこの系は、全体としてエネルギーは増加するのでしょうか…。
No.3ベストアンサー
- 回答日時:
>無理やり毎ステップにエネルギー保存則を適用して速度を収束させていますが
保存則の式や収束のアルゴリズムがわからないので何ともいいがたいです。v=一定なので,方向さえ変えてやればいいわけですが…。
>これはxyz座標系の限界、というべきものなのでしょうか。…円筒座標ですべきなのでしょうか。
直線電流の場合,ローレンツ力はφ成分を持たないため,φ方向の運動について角運動量保存が成立します。この対称性があるためにこの問題では,円筒座標が計算ステップが少なくて有利であるように思われます。円電流もその対称性から,円筒座標または球座標の方がエレガントになるような気がしますが,実際やってみないとわかりません。
![「直線電流周辺での荷電粒子の運動」の回答画像3](http://oshiete.xgoo.jp/_/bucket/oshietegoo/images/media/0/1289033_5497e4c39bced/M.jpg)
円筒座標は外積計算が不安だったため保留にしていましたが、やはりそちらの方が良さそうですね。。
御意見を参考にやり直してみることにします、ご丁寧にどうもありがとうございました!
No.2
- 回答日時:
大切なことを忘れました。
粒子が受ける力は必ず速度と垂直な方向ですから,仕事をされることがありません。したがって速さは一定でエネルギーが増加することはないはずです。
この回答への補足
ご丁寧に画像まで付けて下さりありがとうございます。
以下、問いとは少し離れた質問になってしまうため無理にお答えして頂く必要はございませんが気になったので…お時間がありましたら返答をお願いいたします。
xyz座標系でRunge-Kutta法を用いて計算した結果、無理やり毎ステップにエネルギー保存則を適用して速度を収束させていますがどうしてもrが徐々に増大してしまいます。また、時間刻みによって周期が変わってしまいます。
これはxyz座標系の限界、というべきものなのでしょうか。いずれこれを直線電流から円形電流に変えて粒子シミュレーションを行いたいと考えていますが、円筒座標ですべきなのでしょうか。
No.1
- 回答日時:
円柱座標(r,φ,z)を用いると少し見通しよくなりますが,解析的には解けないかもしれません。
以下「'」で時間微分を表します。r成分:r''-rφ'^2 = -kz'/r …(1)
φ成分:1/r・d/dt(r^2φ') = 0 …(2)
z成分:z'' = kr'/r …(3)
(2)より,r^2φ' = h = r0v0(一定)
∴ φ' = h/r^2
(3)より,z' = k ln(r/r0)
以上を(1)に代入して,
r'' = h^2/r^3 - k^2/r ln(r/r0)
これが解ければいいわけですが,解析的には無理っぽいですね。
両辺にr'をかけて積分すると
r' = ±√[ h^2(1/r0^2 - 1/r^2) - {k ln(r/r0)}^2 ]
ここから???です。
h=r0=v0=k=1とした数値積分結果を参考として添付します。
![「直線電流周辺での荷電粒子の運動」の回答画像1](http://oshiete.xgoo.jp/_/bucket/oshietegoo/images/media/0/1289033_5497f004421bb/M.jpg)
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 荷電粒子の運動の様子と力のつりあいの式です。qの正負によって、粒子の速度Vがy軸対象に変化しますが、 1 2022/04/17 22:23
- 物理学 磁場と磁束密度についての質問です。 4 2022/10/30 01:09
- 物理学 風車が回転する理由 5 2023/05/08 08:03
- 工学 三相誘導電動機の比例推移に関する問題です。 定格出力200kW、4極、50Hzの巻線形三相誘導電動機 1 2023/05/28 12:29
- 大学・短大 電気電子の問題です。この問題が分からないので教えて頂きたいです 幅2mm、長さ5mm、厚さ0.2mm 1 2023/06/21 22:24
- 物理学 shock drift accelerationの加速方向は衝撃波に対して垂直ですか平行ですか? 1 2022/09/19 15:11
- 工学 直流モーターにかける電圧を上げると、モーターの回転数が上がるってどんな理屈? 7 2023/04/03 06:20
- 物理学 電磁気学の問題がわかりません。 3 2023/07/20 22:13
- 物理学 荷電粒子の運動 1 2023/01/14 15:29
- 物理学 一辺がLの正方形のコイルです。 上の方の図でコイルに直流電源Eを接続したところ、コイルには電流iが流 2 2023/06/03 23:58
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
トンネルの中の質点
-
質量m 半径aの一様な円環の慣性...
-
クリストッフェル記号というか...
-
流体力学 円筒座標系
-
スパン方向とはどの方向ですか?
-
半円周上を移動した時の力がす...
-
角運動量についてです 慣性モー...
-
表式ってなんですか?数学用語?
-
英語で位置は何というんでしょう?
-
慣性モーメントを求める問題が...
-
変位座標から角速度・角度の求め方
-
バネでつながれた二つの球
-
直交座標系で表す熱伝導方程式...
-
特殊相対性理論におけるトンネ...
-
SPEC 時間を止めて、銃で放たれ...
-
d軸インダクタンス・q軸インダ...
-
極座標を用いて数値計算する問...
-
同期電動機の電圧方程式について
-
相対性理論(同時性の破れ)に...
-
段差を乗り越えるのに必要なト...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
スパン方向とはどの方向ですか?
-
高校物理基礎で、変位と位置の...
-
質量m 半径aの一様な円環の慣性...
-
表式ってなんですか?数学用語?
-
測量図の座標の出し方
-
径方向?放射方向?
-
原点中心とする半径10cmの演習...
-
2つのバネに挟まれた物体の振動...
-
段差を乗り越えるのに必要なト...
-
「陽に含まない」について
-
さらに・・4次元距離って?
-
鉛直面内での、円運動を考える...
-
ばね定数の異なる3つのバネと...
-
2物体の運動を重心系で考えると...
-
電磁気学、TEMモードでマクスウ...
-
SPEC 時間を止めて、銃で放たれ...
-
力学の問題です。(1)しかわかり...
-
円筒座標系でのナブラ、ラプラ...
-
大学物理の問題が助けてください><
-
有限要素法の形状関数とは?
おすすめ情報