• 締切済み

ベッセル関数を用いた解析について

円柱状の筒に分布した電界強度を、ベッセル関数で解析しようとしています。円筒の底をXY平面としたXYZ直交座標を設定し、あるZでのXY平面を解析面としています。 しかしそもそも、具体的にベッセル関数で解析するとはどういう事なのでしょうか?? また電界強度分布が径方向により異なる場合、ベッセル関数で解析することが出来るのか、出来るのであればどう解析を進めればよいのでしょうか??? どのような事柄でも構いません。皆さんのご意見をお聞かせ下さい。

みんなの回答

回答No.1

ベッセル関数は円柱の三角関数ではないかと思います。フーリエ展開みたいなのが円柱でできるのでは。

mashu75
質問者

お礼

回答を頂きありがとうございます。 ベッセル関数への理解を深めるため、フーリエ展開について再度学習してみたいと思います。

関連するQ&A

  • 解析学(ベッセル関数について)おねがいします!

    大学の解析学の授業でベッセル関数について課題が出され困っているところです。その内容ですがまずどんな関数か定義を書くのですがこれは本で調べたのでなんとなくわかったのですが次の直交関数系の証明、完全系かどうかの証明がまったくわかりません。解き方あるいは参考となるようなアドバイスを頂ければありがたいのですが・・・どなたかお願い致します!

  • 平面波exp(-jx)をベッセル関数を用いてあらわすと・・・

    円筒波動関数について勉強しています。 今、平面波exp(-jx)を円筒座標系で表そうとしているのですが、その変換式が Σa*Jn(ρ)*exp(jnΦ) (Σはn=-∞~n=∞まで)と表されています。 aは定数、Jn(ρ)はn次の第一種ベッセル関数、ρは円筒座標系の原点から外に伸びていく変数、Φは円筒座標系のxy面上の角度 この式について、わからないことがあります。 なぜこの式がx方向に進む平面波を表すのでしょうか?定性的なことが理解できません。 しかも第一種ベッセル関数は進行波でないのに進行波をあらわしている。 このこともさらに混乱を深めています。 どのように理解すればよいのでしょうか? あまりベッセル関数に関する知識がないのでできれば優しくおねがしします。 勉強している本はR.F.Harringtonのtime-harmonic electromagnetic fields です。

  • 複素数を引数とする(?)ベッセル関数

    皆様,お忙しい所回答お願いいたします。 一度他の方が似たような質問があったのですが,それだけでは解決できなかったので再度質問させていただきます。 質問はタイトルの通り,複素数を引数zとするベッセル関数Jn(z)がどうなるのか分かりません。 さまざまな参考書では実数もしくは純虚数をを引数とするベッセル関数の理論やプログラムのサブルーチンはあるのですが,複素数に関しては見つかりません。 大きな数学辞典を見ても,載っている数表は引数が実数のものばかりです。 どうしてこのような関数を必要としているのかというと,電磁界の円柱散乱問題の所で円柱媒質が導電率σを持つ損失性媒質の場合,波数kが複素数となり円柱内部電磁界の解析解に含まれるベッセル関数の引数が複素数となってしまうからです。(Jn(kr)という風に) 複素数を引数zとするベッセル関数Jn(z)の理論について,ご教授の方何卒お願いいたします。 また,参考文献等ありましたらそちらもご教授お願いいたします。

  • ANSYSでの静電場解析

    つい先日ANSYSで2D解析をすることになりました。 電極間に導電体(金属)が挿入されたときの電界強度分布が 見たいのですが、未熟なもので境界条件などの設定がわかりません。 導電体の比誘電率などは設定しなくてもいいのでしょうか?? 2Dの静電場解析ですが何から手を付けていいのかよくわかりません。 どなたか、お願いします。

    • 締切済み
    • CAE
  • 有限要素法による2次元解析について

    有限要素法による2次元解析について 「RF Interstitial加温方式における温度分布計算-2次元モデルと3次元モデルとの比較-」 という論文の要旨を読んでいたのですが、p10の2.1の項目で 電界強度分布Eを有限要素法を用いて数値的に解き・・ とあるのですが、 実際どのように解いているのでしょうか。

  • 光の回折の問題です

    光の回折の問題です たとえばHechtの光学の教科書を読むと,x方向に伸びる1次元のスリットのフラウンホーファー回折像はsinc(x)関数で表現できるとあります.一方,xy平面に空く2次元の円形開口の場合は円筒座標で考えると,半径r方向に対してJ(r)で表現できるベッセル関数になるとあります.確かに式を追えばベッセル関数になるのは数学的には理解できるのですが,直観的に納得ができません.1次元スリットは2次元の円形開口のy=0の中心部分を切り取った場合に相当するはずなので,2次元の場合でもr方向に対してはsinc関数になってもよい気がします.なぜ1次元の場合はsinc関数なのに2次元の場合は円筒方向に対してsinc関数にならずベッセル関数になるのか,直観的な説明はできるのでしょうか.確かにsincもbessel関数も形状は似ていますが.

  • ラプラス方程式の解析解

    電磁気学を勉強しているのですが,分からないことがあるので質問させてもらいます. 静電場内にある電荷が作る電位分布を示す方程式としてラプラス方程式(∇^2*V=0)があると思います. ラプラス方程式とポアソン方程式の違いまでは理解できていると思います. 2次元のラプラス方程式は以下の式を変数分離法を用いて解くことで,直交座標系や球面座標系として考えることで,解析解が得られると理解しています. (ここまではたどり着くことが出来ました) (∂^2/∂^2x)V(x,y)+(∂^2/∂^2y)V(x,y)=0 分からないのは,ここから実際の電位分布を求める方法です. 具体的には,xy平面上の原点にポテンシャルV0がある場合,このV0による電位分布を求めることが出来ません. 直交座標系で考えると一般解は,A,B,C,D,kを定数として,次のようになると思います. V(x,y)=(A*exp(kx)+B*exp(-kx))*(Csinky+Dcosky) 境界条件から未知定数を求めたいのですが,うまくいきません・・・. 原点にポテンシャルがあるので,x→∞でV→0,y→∞でV→0,x=0,y=0でV=V0が境界条件になると思ったのですが,y→∞で(Csinky+Dcosky)は0に収束しません. 境界条件の設定が間違っているのでしょうか? 数値解では原点にポテンシャルを設定している解説は見つけられたのですが,解析解では資料がなく,どうすればいいか困っています. すみませんが,教えてください.

  • 電磁気学

    xy平面に平行で、z=±αに存在する2枚の無限大平面上に均一な面電荷密度で電荷が分布している。 各部の電束密度と電場を求めよ。 と 内半径a、外半径bの無限長中空円筒に均一な電荷密度で電荷が分布している。各部の電場を求めよ。 がいくら考えてもわかりません(*_*) 教えてください(___)

  • 電磁気学に関する問題についての回答お願いします。

    電磁気学の問題になります 真空中に図のような半径a[m]の円柱導体とそれを取り囲む半径c[m]の円柱導体よりなる無限長同軸導体がある。円柱導体の周囲は中心より半径b[m]の範囲まで正の一様な電荷密度ρ[C/m^3]で満たされている。円筒導体は接地されており、厚さは考えなくて良い。中心軸からの距離をr[m]として、次の問いに答えよ。ただし、真空の誘電率をε0とする。 (1)中心軸を垂直に横切る断面において、r>b範囲での電気力線の様子を描け。 (2)中心軸を円筒座標系のz軸にとり、rにおける電界E(r)[V/m]を求めよ。 (3)電界の強さE(r)を縦軸、rを横軸にとり、電界のrに対する変化の様子をグラフに描け。 (4)中心から距離aおよび、bの点における電位Va[V]および、Vb[V]を求めよ。 (1)については、円柱の上面、底面には電界はなく、放射状に単位ベクトルer方向の電界が出来る感じでいいのでしょうか? (2)は、ガウスの法則の積分系と微分系の2つのやり方でやってみたのですが、2つとも考え方が正しいのかご指摘お願いします。細かいところまで指摘してほしいため、出来るだけ詳細に考え方を書きました。 ・積分系によるやり方 r<aの範囲について r<aの範囲には電荷分布がないので、E=0 a<r<bの範囲について 導体なので、内部には電荷はなく表面に電荷が分布して、E=0 ちなみに、電荷は外側の表面だけに集まるのでしょうか?この場合、円筒導体は接地されているので、静電気的な力で。もし接地されていなければ、外面と内面に半分ずつに電荷が分布するのでしょうか?どのような割合で分布するのでしょうか? b<r<cの範囲について ガウスの積分法則∫E・dA=(Q/ε0)を用いて、Aは円柱の表面積なので、2πrz。 Qは、ρが表面に集まるという考え方が正しいとすると、Q=(2πbz)ρ。 ただ、このやり方だと、b*zということで、次元は、[M^2]で、ρは、[C/M^3]の次元で、Q=[C/M]の次元となり、間違っているとも思うのですが。。。 よって、E(2πrz)=(1/ε0)(2πbzρ)より、E=(b*ρ)/(ε0*r) r>cの範囲について 円筒導体が接地されていることにより、円柱導体の周りの電荷と等量の負の符号の電荷が分布すると考えたので、円柱状のガウス面を考えてやると、相殺して、外側の電界は、E=0となりました。 ・微分系によるやり方 b<r<cの範囲について divD=ρという微分系のガウスの法則を用いる。φ方向や、z方向は電界は一定で、rだけに依存すると考えられるので、円柱座標のdivDの式をつかって、(1/r)d(rDr)/dr=0より、rDr=C(ただし、Cは、積分定数)。 r=bでは、導体と真空という誘電体の境界なので、Dr=ρとなって、b*ρ=C よって、rDr=b*ρとなり、Dr=(b*ρ)/rとなりました。 よって、E=(b*ρ)/(ε0*r) 勘違いしているところや、おかしな考え方をしているところはどこでしょうか?ご指摘と、なぜ間違っているか理由を教えてください。 (3)については、円柱と円筒の間では、1次関数的に電界の強さは、下がっていき、その他の場所ではE(r)=0というグラフでよろしいでしょうか? (4)は、E=-gradVから得られEを積分して求めればよいとは思うのですが、積分範囲をどうすればいいのかよく分かりません。 回答よろしくお願いいたします。

  • ベクトル解析の面積ベクトルを学習しているのですが

    ベクトル解析の面積ベクトルの正射影の面積について XYZ空間内に平面πを定めてこのπ上に平曲線cで囲まれる図形をDとおきその面積をsとおく。  このとき、平面πに垂直で大きさ1の正の向きのベクトルを単位法線ベクトルと呼び、これをnと表すことにする。すると面積ベクトルS=snとなる。このときn=[cosα,cosβ,cosγ](0≦α≦π,0≦β≦π,0≦γ≦π)さらに、基本ベクトルi=[1,0,0],j=[0,1,0],k=[0,0,1]とするとDのxy平面への正射影の面積は|i・S|=s|cosα|となる。 (jS,kSは省略) ここで、平面πの定め方について疑問があります。まずxy平面と平行な平面πを考えます。 このとき単位法線ベクトルnはz軸と平行です。 そしてここからが問題ですが、平面πを生成するベクトルを考えます。このベクトルの中のひとつをaベクトルとしてaベクトルとx軸との角度はαとします。そして、aベクトルを回転軸に平面πを回転させます。こうすると、この平面πはαβγだけで表すことができるのでしょうか? また、正射影を考えたときにその面積は、|i・S|=s|cosα|にはならないと思うのですが勘違いしているかもしれないので、どなたか詳しく教えて頂けないでしょうか?