[論文レビュー] Some integral identities involving products of general solutions of Bessel's equation of integral order

S. K. H. Auluck|arXiv (Cornell University)|Jun 23, 2010

Differential Equations and Boundary Problems参考文献 2被引用数 4

ひとこと要約

本稿は、整数次のベッセル方程式の一般解の積を含む解析的積分恒等式を導出する。特に2関数および3関数の積に注目し、これらは円柱幾何におけるスペクトル分解を可能にし、空間的依存性を排除することで非線形偏微分方程式を時間発展方程式に変換する。主な貢献は、ベッセル関数の三重積積分に対する閉形式の表現であり、数値的に検証され、高次のモード数においても高い正確性を示す。応用分野はプラズマ物理学および流体力学に及ぶ。

ABSTRACT

Spectral decomposition of dynamical equations using curl-eigenfunctions has been extensively used in fluid and plasma dynamics problems using their orthogonality and completeness properties for both linear and non-linear cases. Coefficients of such expansions are integrals over products of Bessel functions in problems involving cylindrical geometry. In this paper, certain identities involving products of two and three general solutions of Bessel's equation have been derived. Some of these identities have been useful in the study of Turner relaxation of annular magnetized plasma [S.K.H. Auluck, Phys. Plasmas, 16, 122504, 2009], where quadratic integral quantities such as helicity and total energy were expressed as algebraic functions of the arbitrary constants of the general solution of Bessel's equation, allowing their determination by a minimization procedure. Identities involving products of three solutions enable expanding a product of two solutions in a Fourier-Bessel series of single Bessel functions facilitating transformation of partial differential equations representing non-linear dynamics problems into time evolution equations by eliminating spatial dependences. These identities are required in an ongoing investigation of fluctuation-driven coherent effects in a nonlinear dynamical system.

研究の動機と目的

• 整数次のベッセル方程式の一般解の積を含む解析的積分恒等式を導出すること。
• 固有関数の積をフーリエ・ベッセル級数として表現することにより、円柱幾何における非線形偏微分方程式のスペクトル分解を可能にすること。
• 高次のモード数領域における数値的困難を克服するため、ベッセル関数の三重積分の解析的表現を提供すること。
• 特に揺動駆動型の秩序的効果および環状磁化プラズマにおけるターナー緩和を含む非線形力学系の研究を支援すること。
• 数値的統合が急激な振動により失敗する高次のモード数において、三重積分係数の近似スケーリング関係を導出すること。

提案手法

• 第一種および第二種ベッセル関数の既知の微分的恒等式を用いて、ベッセル関数の積の積分に関する再帰的関係を導出する。
• 境界条件を伴う環状領域における直交性条件を適用し、ベッセル関数の任意の線形結合を含む定積分恒等式を導出する。
• β → α の極限を用いて、一般解のノルムの二乗を表す式を導出し、任意定数を用いたエネルギーおよびヘリシティの表現を可能にする。
• 大引数におけるベッセル関数の漸近的近似を用いて、三重積分の近似スケーリング関係を導出する。
• フレネル積分および位相一致条件を用いて、三重積分を位相依存の振幅を持つ振動関数の形で表現する。
• 直接積分との比較により、特に振動が激しく収束に困難を伴う高次のモード数において、解析的表現の数値的妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 境界条件を伴う環状幾何において、ベッセル方程式の2つの一般解の積を閉形式の積分としてどのように表現できるか？
• RQ2 非線形偏微分方程式のスペクトル分解を可能にする、ベッセル関数の三重積の解析的表現は何か？
• RQ3 高次のモード数において、三重積分の漸近的近似の正確性はどの程度で、そのスケーリング特性は何か？
• RQ4 導出された恒等式を用いて、エネルギーおよびヘリシティといった二次積分量を一般解の任意定数の代数的関数として表現できるか？
• RQ5 直接積分が急激な振動により失敗する状況において、三重積分係数の近似公式の数値的信頼性はいかがなものか？

主な発見

• 本稿では、一般ベッセル関数解の二乗の積分に対する閉形式の恒等式を導出し、任意定数 a および b を用いた全エネルギーおよびヘリシティの代数的表現を可能にした。
• 隣接する次数のベッセル関数の積の積分を結ぶ再帰的関係が確立され、高次次数の積分を再帰的に計算可能となった。
• n=0 および n=1 の場合に、一般解のノルムを J₀² および J₁² 項の和として明示的に導出し、微分を含む補正項を含めた。
• 大引数におけるベッセル関数の振る舞いとフレネル積分近似を用いて、三重積分係数 Cₘₙₚ¹⁰⁰⁰ の近似スケーリング関係を導出した。これは高次のモード数において有効である。
• 数値的検証では、低次のモードでは 10⁻⁸ から 10⁻²² の範囲で一致が確認され、高次のモード数においてはさらに高い正確性を示したが、振動的被積分関数および近似の限界により若干の乖離が生じた。
• 高次のモード数（例：>150）においては、収束不能となる数値的統合に比べ、近似公式が優れた性能を示し、想定を上回る一致を示したが、残存誤差は依然として存在した。