[論文レビュー] Radiative transfer in plane-parallel media and Cauchy integral equations III. The finite case

B. Rutily, L. Chevallier+1 more|arXiv (Cornell University)|Jan 16, 2006

Radiative Heat Transfer Studies参考文献 13被引用数 6

ひとこと要約

本稿では、等方的散乱を伴う有限で平面平行な媒質における放射線伝搬の新しい解析的解法を提示する。この手法では、一般化された因数分解関係を満たす二つの補助関数 $\zeta_{+}$ と $\zeta_{-}$ を導入し、これらは非結合なフレッドホルム積分方程式を介して計算される。この方法により、複素平面全域における $X$-関数および $Y$-関数の高精度な計算が可能となり、半無限媒質の解を有限な平板に拡張する際の数値的安定性が向上する。

ABSTRACT

We come back to the Cauchy integral equations occurring in radiative transfer problems posed in finite, plane-parallel media with light scattering taken as monochromatic and isotropic. Their solution is calculated following the classical scheme where a Cauchy integral equation is reduced to a couple of Fredholm integral equations. It is expressed in terms of two auxiliary functions $ζ_+$ and $ζ_-$ we introduce in this paper. These functions show remarkable analytical properties in the complex plane. They satisfy a simple algebraic relation which generalizes the factorization relation of semi-infinite media. They are regular in the domain of the Fredholm integral equations they satisfy, and thus can be computed accurately. As an illustration, the X- and Y-functions are calculated in the whole complex plane, together with the extension in this plane of the so-called Sobouti's functions.

研究の動機と目的

• 有限で平面平行な平板における等方的散乱を伴う放射線伝搬に生じる特異なコーシー積分方程式を解くこと。
• 半無限媒質の因数分解関係を、新しい補助関数を用いて有限幾何に一般化すること。
• 有限平板における複素平面全域における $X$-, $Y$-, $\xi_X$-, および $\xi_Y$-関数の高精度な計算を可能にすること。
• $\zeta_{\pm}$ 関数を介して伝搬方程式をフレッドホルム方程式に還元することにより、有限媒質における解法を体系的かつ包括的に確立すること。

提案手法

• 標準的な積分方程式技法を用いて、元のコーシー積分方程式を $[-1, +1]$ 上から $[0,1]$ 上の二つの結合した方程式に変換する。
• 半無限の場合を一般化するため、一般化された代数的因数分解関係を満たす二つの補助関数 $\zeta_{+}$ と $\zeta_{-}$ を導入する。
• $[0,1]$ 上で二つの非結合な第二種フレッドホルム積分方程式を用いて $\zeta_{+}$ と $\zeta_{-}$ を解き、正則性と数値的精度を保証する。
• $\zeta_{+}$ と $\zeta_{-}$ の解を用いて、明示的な代数的式（式 50–53）により $X$-関数および $Y$-関数を再構築する。
• $c_0 = 1$ の場合に基本的な $X$ および $Y$ 関数を導出し、その後式 (57) の重ね合わせ公式を用いて一般の $c_0$ へ拡張する。
• $\zeta_{+}$ および $\zeta_{-}$ が $\mathbb{C}^*$ 上で正則であり、$[0,1]$ 上で正則であることを示し、関連するフレッドホルム方程式の数値的評価の安定性を保証する。

実験結果

リサーチクエスチョン

• RQ1 有限平板における放射線伝搬の特異なコーシー積分方程式は、補助関数を用いてどのように解ける形に還元できるか？
• RQ2 複素平面における新しい関数 $\zeta_{+}$ および $\zeta_{-}$ が有する解析的性質は何か？ また、それらは半無限媒質の因数分解関係をどのように一般化するか？
• RQ3 $\zeta_{+}$ および $\zeta_{-}$ を用いて、有限平板における $X$-および $Y$-関数を高精度に計算できるか？ また、境界 $u \to 0^+$ および $u \to 1^-$ における挙動はいかなるものか？
• RQ4 一般の源 $c_0$ に対する解と $c_0 = 1$ の基本解との間の関係は何か？
• RQ5 $H_+(a,b,z)$ および $H_-(a,b,z)$ が $H(a,z)\zeta_{\pm}(a,b,z)$ として定義されるとき、$b \to \infty$ の極限で古典的 $H$-関数にどのように回復するか？

主な発見

• $\zeta_{+}$ および $\zeta_{-}$ は、半無限媒質の因数分解 $T(a,z)H(a,z)H(a,-z)=1$ を有限平板に拡張する一般化された因数分解関係（式 40）を満たす。
• $X$-および $Y$-関数は、明示的な代数的式（式 50–53）を用いて $\zeta_{\pm}$ 関数から高精度に計算可能であり、$u=0$ や $u=1$ の近傍における数値的不安定性を回避できる。
• 補助関数 $\zeta_{+}$ および $\zeta_{-}$ は $[0,1]$ 上で正則であり、$\mathbb{C}^*$ 上で正則であるため、関連するフレッドホルム方程式の数値的解法が信頼性を保証される。
• 一般の $c_0$ に対する解は、式 (57) を用いて $\zeta_{\pm}$ と源関数 $c_0$ を用いた重ね合わせとして表現され、$X_0$ の完全な再構築が可能になる。
• $H_+(a,b,z) = H(a,z)\zeta_+(a,b,z)$ および $H_-(a,b,z) = H(a,z)\zeta_-(a,b,z)$ は一般化された因数分解（式 98）を満たし、$b \to \infty$ の極限で古典的 $H$-関数に還元される。これは既知の極限と整合していることを確認する。
• 本手法は $X$ および $Y$ 関数を中間ステップとして直接用いないため、従来のアプローチよりもより安定的かつ洗練された解法への道筋を提供する。