(57)【要約】  本発明は、少なくとも１個の波と、各物体の表面が少なくとも２個の媒質間の界面を画定している、少なくとも１個の物体との間の系内の相互作用をモデル化する方法に関し、本方法は以下のステップを含むことを特徴とする。  −対象とする適用分野に対応すべく基本特性関数の集合を選択する（Ｅ１）  −系を構成すると考えられる各媒質の物理特性を定義する（Ｅ２）  −系の各物体をメッシュとしてモデル化されることにより生成し、各メッシュ要素のいずれかの側に少なくとも１個の基本点源を関連付ける  −各界面について境界条件の種類を決定する  −境界条件の種類、媒質の特性、および系の構成に依存して各種の物体間の相互作用を表わす大域行列を構築する  −大域行列の逆行列を計算する  −逆行列に、励起境界条件の値を含んでいる列行列を乗算する  −全ての基本点源の値を含む列行列を求める  −系内の全ての点において、相互作用を表わす物理量を計算し、系内の相互作用の解析的なモデルを得る。【選択図】図１
2009521031.tif 000002
【特許請求の範囲】
【請求項１】
少なくとも１個の波と、各物体の表面が少なくとも２個の媒質間の界面を画定している、少なくとも１個の物体との間の系内の相互作用をモデル化する方法であって、
−対象とする適用分野に対応すべく基本特性関数の集合を選択するステップ（Ｅ１）、
−系を構成すると考えられる各媒質の物理特性を定義するステップ（Ｅ２）、
−系の各物体の幾何学的構造がメッシュとしてモデル化することにより画定され、各々のメッシュ要素の表面に少なくとも１個の試験点を配置するステップ、
−各試験点において、境界条件に対して連続式を決定すべく考慮される各媒質について少なくとも１個の試験量を定義するステップ、
−各メッシュ要素のいずれかの側に少なくとも１個の基本点源を関連付けるステップ、
−物体を空間内で互いに相対的に配置するステップ、
−物体により境界が画定されたボリュームに媒質を関連付けるステップ、
−各界面について境界条件の種類を決定するステップ（Ｅ７）、
−各種の物体間の相互作用を表わす大域行列を構築するステップ（Ｅ８）であって、前記行列が一対にされた物体間の相互作用を特徴付ける少なくとも１個の行列ブロックから形成され、これらの相互作用は選択された特性基本関数を介して対象とする媒質中の波の伝播に関連付けられ、大域行列は多くても、一対にされた全ての物体間の可能な組合せの数のブロックを含み、各ブロックの内容は試験点において指定された境界条件の種類、対象とする２個の物体に共通な媒質の特性、およびこれらの物体の幾何学的形状に依存するステップ、
−大域行列の逆行列を計算するステップ、
−前記逆行列に、ユーザが課した励起境界条件の値および必要に応じて固有境界条件に対応するゼロを含んでいる列行列を乗算するステップ、
−全ての基本点源の値を含む列行列を求めるステップ（Ｅ９）、
−系内の全ての観測点において、系内の相互作用を表わす物理量を、点源について考慮する影響の範囲に応じて計算するステップ、
−系内の相互作用の解析的なモデルを得るステップ
の各ステップを含むことを特徴とする方法。
【請求項２】
前記試験点が、伝播の少なくとも一つの特定方向に偏ることを避けるように、あるメッシュ要素から次のメッシュ要素まで無作為に分布していることを特徴とする、請求項１に記載の方法。
【請求項３】
前記系の少なくとも一部において巨視的量が計算されることを特徴とする、請求項１または２に記載の方法。
【請求項４】
前記基本点源の全てにより生成された物理量が表示される（Ｅ１３）ことを特徴とする、請求項１〜３のいずれか一項に記載の方法。
【請求項５】
表面が常に閉じていて前記物体外の媒質と前記物体内の媒質との間に境界を形成する閉ボリューム物体が画定されていると共に、半無限媒質との界面を表わすことに起因してその表面が開いている開ボリューム物体が画定される（Ｅ６）ことを特徴とする、請求項１〜４のいずれか一項に記載の方法。
【請求項６】
前記界面が生成される際に、内部での相互作用が調べられているボリュームの境界の集合により横方向に境界が画定されることを特徴とする、請求項５に記載の方法。
【請求項７】
２種類の境界条件、すなわち固有境界条件およびユーザが指定した境界条件が用いられ、第１の境界条件が、決定されたまたは決定可能な特性を有する２個の媒質間の界面を越える際のスカラーおよび／またはベクトル量の連続性を表わし、第２の境界条件が、界面のスカラーおよび／またはベクトル量についてユーザが有する先験的知識を表わすことを特徴とする、請求項１〜６のいずれか一項に記載の方法。
【請求項８】
ユーザが指定した境界条件が、界面一方の側に位置するソースの単一の組により満たされることを特徴とする、請求項７に記載の方法。
【請求項９】
前記特性関数が異方性媒質内で計算されることを特徴とする、請求項１〜８のいずれか一項に記載の方法。
【請求項１０】
前記基本特性関数のうち少なくとも１個がグリーン関数であることを特徴とする、請求項１〜９のいずれか一項に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、少なくとも１個の波と、各物体の表面が少なくとも２個の媒質間の界面を画定している、少なくとも１個の物体との間の相互作用をモデル化する汎用的方法に関する。
【背景技術】
【０００２】
仏国特許第２８４７０５１号明細書によれば、波と障害物との間における調査対象である相互作用の物理量で表現可能にする方法が知られている。この方法の実装は、非破壊試験等の用途に極めて有用である。しかし、これは１個の媒質を別の媒質から分離する単一界面の場合にしか適用できない簡便法である。実際、この方法は特に、界面を越える際の入射超音波により生じる反射および透過が極めて多いため、複数の媒質を含む物体との相互作用に一般化することができない。
【発明の開示】
【発明が解決しようとする課題】
【０００３】
本発明の目的は、少なくとも１個の波と、各物体の表面が物理特性の異なる少なくとも２個の媒質間の界面を画定している、少なくとも１個の物体との間の相互作用をモデル化する方法である。本方法は、実装が簡単であって研究対象である系の構成如何にかかわらず、様々な物理学の分野における連続、正弦、またはパルス時間領域に対して適用可能である。
【０００４】
以下の説明全体を通じて、本発明における意味において「系」という概念を、モデル化したい相互作用が生じる空間領域として用いる。
【課題を解決するための手段】
【０００５】
本発明による方法は、少なくとも１個の波と、各物体の表面が少なくとも２個の媒質間の界面を画定している、少なくとも１個の物体との間の系内の相互作用をモデル化する方法であって、以下のステップを含むことを特徴とする。
−対象とする適用分野に対応すべく基本特性関数の集合を選択する（Ｅ１）
−系を構成すると考えられる各媒質の物理特性を定義する
−系の各物体の幾何学的構造がメッシュとしてモデル化することにより画定され、各々のメッシュ要素の表面に少なくとも１個の試験点を配置する
−各試験点において、境界条件に対して連続式を決定すべく考慮される各媒質について少なくとも１個の試験量を定義する
−各メッシュ要素のいずれかの側に少なくとも１個の基本点源を関連付ける
−物体を空間内で互いに相対的に配置する
−物体により境界が画定されたボリュームに媒質を関連付ける
−各界面について境界条件の種類を決定する
−各種の物体間の相互作用を表わす大域行列を構築する、すなわち当該行列が一対にされた物体間の相互作用を特徴付ける少なくとも１個の行列ブロックから形成され、これらの相互作用は選択された特性基本関数を介して対象とする媒質中の波の伝播に関連付けられ、大域行列は多くても、一対にされた全ての物体間の可能な組合せの数のブロックを含み、各ブロックの内容は試験点において指定された境界条件の種類、対象とする２個の物体に共通な媒質の特性、およびこれらの物体の幾何学的形状に依存するように構築する
−大域行列の逆行列を計算する
−当該逆行列に、ユーザが課した励起境界条件の値および必要に応じて固有境界条件に対応するゼロを含んでいる列行列を乗算する
−全ての基本点源の値を含む列行列を求める
−系内の全ての観測点において、系内の相互作用を表わす物理量を、点源について考慮する影響の範囲に応じて計算する
−系内の相互作用の解析的なモデルを得る。
【０００６】
本発明の別の特徴によれば、試験点は、伝播の少なくとも１つの特定方向に偏ることを避けるように、あるメッシュ要素から次のメッシュ要素まで無作為に分布している。
【０００７】
本発明の別の特徴によれば、系の少なくとも一部において巨視的量が計算される。本発明の別の特徴によれば、基本点源の全てにより生成された物理量が表示される。
【０００８】
本発明の別の特徴によれば、表面が常に閉じていて物体外の媒質と物体内の媒質との間に境界を形成する閉ボリューム物体が画定されていると共に、半無限媒質との界面を表わすことに起因してその表面が開いている開ボリューム物体が画定されている。この界面が生成される際に、内部での相互作用が調べられている作用空間の境界の集合により横方向に境界が画定される。
【０００９】
本発明の別の特徴によれば、２種類の境界条件、すなわち固有境界条件およびユーザが指定した境界条件が用いられ、第１の境界条件が、決定されたまたは決定可能な特性を有する２個の媒質間の界面を越える際のスカラーおよび／またはベクトル量の連続性を表わし、第２の境界条件が、界面のスカラーおよび／またはベクトル量についてユーザが有する先験的知識を表わす。
【００１０】
本発明は、非限定的な例として添付の図面に示す実施形態の詳細な記述を精査することにより更に理解が深まるであろう。
【００１１】
入射波と均一媒質との間の相互作用をモデル化する方法の原理は、上述の仏国特許第２８４７０５１号明細書で詳述されているため、本明細書で再び記述することは行なわず、単にその略称ＤＰＳＭ（分散点源方法）で言及するに留める。本方法は、非均一媒質への応用を可能にすべく当該公知方法の改良として提示する。
【発明を実施するための最良の形態】
【００１２】
以下に、本発明の方法のいくつかの実施形態について、超音波およびフーコー電流を用いる基準非破壊試験と共に記述するが、本発明がこれらの応用だけに限定されないこと、および電磁気、静磁気、音響、光学、地球物理、フーコー電流、感熱性現象等、極めて広範な場において少なくとも１個の界面を含む媒質における波伝播現象を実装するあらゆる用途に実装できることを理解されたい。
【００１３】
上述したように、本発明の中核を形成する重要な特徴は、考慮している媒質内の量の統合を担う点源の層を用いて、連続的に境界条件を伝播させることにあり、これらのソースの層は各々の界面のいずれかの側に配置されている。
【００１４】
波または波の群が伝播することにより上述の相互作用を引き起こす系の波の性質および物体と媒質の特性に応じて、これらの相互作用のモデリングは、本発明によれば、個別に点状励起を起こされた各々の媒質の応答を定式化することにあり、結果的に考慮する媒質毎に１個の基本特性関数が得られる。当然ながら、使用するソースの物理的な性質は合成しようとする波の性質に対応している。
【００１５】
励起への応答の計算から一般に、スカラー、ベクトル、またはテンソルのポテンシャルおよび波の特性量が得られる。系の全応答は、対象とする媒質内で放射している全ての点源の寄与を合計することにより、媒質毎に計算される。
【００１６】
大域的な解は、対象とする界面に法線方向におけるポテンシャルおよびその導関数の連続性を課すことにより、各々の界面を越える際に境界条件が定数内で満たされることを示唆する。この条件の集合により、上述の大域行列が得られる。
【００１７】
各々の試験点で得られる方程式の数は、界面のいずれかの側に配置される点源の数の選択を条件付ける。
【００１８】
静電気関連における本発明の一例示的実施形態によれば、電位はスカラー・ポテンシャルＶに対応し、および対象とする界面の表面に法線方向におけるその導関数は、電場
【数１】
2009521031.tif 000003
の法線成分に対応しており、連続性を保障すべき量は、電気変位ＤＮ＝εＥＮの界面の表面の法線成分である。１個のソースが各メッシュ要素（すなわちＶおよびＤＮ）のいずれかの側に配置されている場合、このように各試験点で２個の未知数を有する２個の方程式が得られる。これらのソースは静電荷と同等である。
【００１９】
電磁気学において、対象とする媒質が導電性である場合、対象とする界面の表面の法線方向成分と同様に、ポテンシャルがベクトルであることに注意されたい。このように、各試験点に対して、各試験点においてメッシュ要素のいずれかの側に三重項が存在するという事実により、６個の未知数を有する６個の方程式（すなわち、
【数２】
2009521031.tif 000004
の３個の成分と
【数３】
2009521031.tif 000005
の３個の成分）が得られる。これらの三重項を構成するソースは、各々がこれらの三重項の直交座標軸に沿う電流要素と同等である。
【００２０】
系を方程式の形式で記述することは、このように各種の媒質間の全ての相互作用を考慮した単一の大域行列を用いて大域的に定式化される。次いで、各ソースから放射されるエネルギーの値が前述のように逆大域行列およびユーザが指定した境界条件を用いて計算される。
【００２１】
好都合なことに、上述の基本特性関数はグリーン関数である。
【００２２】
以下の３個の説明図において、直線の区間が界面の一部を表わし、そのいずれかの側に、または一方の側に、少なくとも１個の点源を含む球面を示し、これらの表面は界面に接している。それらを含むこれらの点源および球面が、トランスデューサまたは界面のいずれかを、上述の仏国特許第２８４７０５１号明細書の意味においてモデル化するメッシュから得られる。以下の図において、明灰色の円は固有境界条件に対応し、暗灰色の円はユーザが指定する境界条件に対応する。各メッシュ要素は、少なくとも１個の点源を含む球面に関連付けられている。
【００２３】
１つの特定の場合において、これらの球面と界面間の接点が試験点であるが、各々のメッシュ要素内に複数の試験点が分布していてもよいことが明らかに理解されよう。
【００２４】
ソースの未知パラメータの総数は、利用できる方程式の総数に等しくなければならず、方程式の数は対応して選択された試験点の数により調整される。
【００２５】
固有境界条件（ＩＢＣ）
【数４】
2009521031.tif 000006
界面を越える際に、ＩＢＣにより「自然な」連続性が保証される。媒質ｎ＋１に示すソースは媒質ｎに波を伝播する役割だけを果たし、その逆も同様である。
【００２６】
ユーザ境界条件（ＵＢＣ）
【数５】
2009521031.tif 000007
ＵＢＣは特に、系内の励起の存在を表わす。例えば、以下の図は媒質ｎ＋１へ向かって有効なＵＢＣを示す。この例では、媒質ｎ内で量が計算されないため媒質ｎ＋１のソースは存在しないが、これにより、ソースが媒質ｎ内へ放射するように媒質ｎ＋１内にソースを配置する可能性を排除しないことが明確に理解されよう。例えば、静電学または電気力学によれば、コンデンサの１個の電極上のボルトで表わされる静電ポテンシャルは一定であって通常はユーザにより指定される。
【００２７】
混合境界条件（ＩＢＣ／ＵＢＣ）
【数６】
2009521031.tif 000008
【００２８】
これらの境界条件は、上述の画定を採用している。界面をメッシュモデル化する他の点源は、媒質ｎと媒質ｎ＋１の間にＩＢＣを確立する役割を果たすだけである。
【００２９】
境界条件の如何にかかわらず、全ての点源は個別に設定することができ、および固有の放射パターンを有していてよい。
【００３０】
調査対象の相互作用が内部で生じるＮ個の物体の系を形成する各々の媒質の物理特性が定義され、各種の界面がメッシュモデル化されたならば、調査対象の相互作用を表わす大域行列が次式のように構成される。
【数７】
2009521031.tif 000009
【００３１】
この系はＮ個の物体を含んでいるため、行列はＮ×Ｎ個のブロックを含んでいる。これらのブロックすなわち基本行列は、各種の物体間を通過するための条件を表わす。座標（ｉ，ｊ）を有するブロックは、物体ｊの影響に起因して計算された物体ｉの界面上での境界条件の割合を表わす（物体ｉとｊは共通の媒質内で相互作用する）。
【００３２】
本発明の一特徴によれば、ソース成分の総数を試験点の総数に合わせる必要がある。
・物体ｉ（行列の第ｉ行）上にＵＢＣがある場合、ブロック・ベクトルＣＯＮＤＩＴＩＯＮＳｎ○ｉは対象とするＵＢＣの値を含んでいる。
・ユーザ自身がＵＢＣに関連付けられたソースの値を指定する場合、対応するブロックＣＯＮＤＩＴＩＯＮＳはこれらのソースの値に等しい。対応するＵＢＣ行列は単位行列である。
・ＩＢＣだけが物体に適用可能である場合、ベクトルＣＯＮＤＩＴＩＯＮＳの対応するブロックはゼロである。
【００３３】
本発明の一特徴によれば、基本行列［ＩＢＣ］は、以下に述べる例でＭおよびＱと表記された部分行列から形成されている。これらの部分行列は、各々界面における量、例えばベクトル量の法線および接線成分、さもなければ、スカラー量および法線ベクトル成分の連続性を保証するために、適用分野（超音波、電磁気、音響等）に応じて計算される。
【００３４】
同じことがＵＢＣ行列に当てはまるが、量の１つはユーザにより指定される。
【００３５】
本発明の方法の基本原理を提示した後で、次いで図１のフロー図を参照しつつ、異なる特性を有する媒質の境界を画定する物体間の相互作用のモデリングを実際に解くために本発明のプログラムにより実装されている各種のステップについて記述する。このプログラムは対話的であるため、各々のステップにおいてユーザ主導、または実装されているコンピュータ主導のいずれにより実行されているかを記述し、後者は原則的に方程式の計算と表示ステップを実行する。
【００３６】
本発明が電磁気、静磁気、音響、フーコー電流、または熱的現象等の極めて多様な物理場に適用できることから、第１ステップ（ステップＥ１）では対象とする適用分野に対応する基本特性関数の全てをユーザが選択することを求める。
【００３７】
ユーザは次いで、モデル化したい相互作用が内部で生じる媒質の特性を指定する（ステップＥ２）。これらの特性は、媒質の全ての関連パラメータを含んでいる。媒質の特性の種類は、調査対象である問題に依存する。例えば、静電気の場合、比誘電率は、不可欠なパラメータである。このステップでは、各々の物体に隣接する媒質を示す一致指数を計算することができる。物体には、定式化に応じて２個の媒質が関連付けられている。すなわちｉ（α、β）、ここにｉは物体、αとβはｉに隣接する媒質を示す。各々の物体は一致指数を有し、これによりいくつかの簡素化された例に関して以下に示すように相互作用において考慮するソースを決定することができる。
【００３８】
ステップＥ２から出発して、二通りの可能性がある。すなわちステップＥ３、Ｅ４およびＥ５を実行するか、または直接ステップＥ５へ進むかである。ステップＥ３およびＥ４は、波伝播現象を考慮に入れる場合にのみ実行される。ステップＥ３において、ユーザはトランスデューサ・エミッタから発せられる波の周波数を指定する。この周波数が決定されたならば、プログラムは、相互作用の特徴に応じて、点源間の最大距離（換言すれば、ソースの規則的に構成されていればそれらのピッチ）を計算する（ステップＥ４）。このピッチは好適には、これらの相互作用の原因である波の半波長より短い。
【００３９】
ステップＥ５において、ユーザは閉ボリューム物体（以下ＣＶＯと表記）を指定する。このステップは本質的に、そのような物体を幾何学的に指定して、本発明の方法の実装が可能なようにメッシュのピッチを決定するものであり、メッシュは前記仏国ＤＰＳＭ特許の記載に従い生成される。
【００４０】
本発明によれば、閉ボリューム物体（ＣＶＯ）は、表面が必ず閉じている（例えば球体、角錐または立方体）。この表面は、外部媒質と物体の内部にある媒質との間の異なる特性を有する境界である。ユーザが、この媒質内の量の計算を何ら行ないたくない場合、境界条件（ユーザ自身により指定され、以下にＵＢＣと表記する）を当該表面に適用しなければならない。この物体のボリュームはゼロに向かう傾向がある（例えば、コンデンサの平坦な電極の場合）。
【００４１】
再び本発明によれば、開ボリューム物体（ここではＯＶＯと表記）は、この物体の表面が開いている（理論的には、そのボリュームは無限遠で閉じている）。これは、半無限大媒質間の境界を表わす。この構成を表わすために界面が生成されている場合、一組の境界により必ず横方向の境界が画定されている。これらの境界もまた調査対象となって、これらの境界により境界が画定されている空間の外側の量を観測することは一切できなくなる。この作用空間は従って、自動的に境界が画定されている。本発明が考慮する物体および媒質は必ず作用空間に属する。
【００４２】
解決すべき問題が閉ボリューム物体だけを含む場合、作用空間の境界を画定する必要はない。例示的に、２個の帯電された球体により空間内に生成された電場を調べる場合に作用空間を必要としない。同様に、図６〜８を参照しつつ以下に記述する例は空間内に境界が存在しない２個の系を示す。
【００４３】
系全体が物体の内部ＵＢＣを有するＣＶＯ内に含まれている場合、（図４、５を参照しつつ上述した例のように）その表面が作用空間の境界を画定する。
【００４４】
実際に、媒質はその均一な物理特性により画定される。媒質は、それがＣＶＯ内に含まれるため、あるいは１個以上の開ボリューム物体（ここではＯＶＯと表記）および作用空間を囲む境界により境界を画定されるため、必ず閉じている。媒質を画定するには、その局在性および物理特性が既知でなければならない。以下の説明において、媒質を番号で参照する。
【００４５】
物体の表面は、２個の媒質間の界面を画定している。物体を画定するには以下を知る必要がある。
１．次元および局在性
２．表面に適用可能な境界条件
３．物体が接触する媒質
【００４６】
本明細書の記述において、物体を文字により参照する。物体の生成は、これらの物体の表面の、系内で伝播する波を合成する役割を担う点源の集合が関連付けられている点Ｐの集合への離散化に対応している。これらの物体は、境界条件（ＢＣ）が適用される界面を画定する。
【００４７】
ステップＥ５から出発して、ユーザは、ステップＥ７に到達するためにＯＶＯ（そのような物体が存在する場合）を画定すべくステップＥ６経由で進むか、さもなければ直接ステップＥ７へ進む。
【００４８】
系の全体または一部の画定および／またはメッシュ・モデリングは、本発明対象外の装置またはソフトウェア・アプリケーションにより提供することができる。
【００４９】
ステップＥ７は、境界条件の決定からなり、これは本発明の方法を実装する基本条件の１つをなす。これらの境界条件は、上述のように、異なる特性を有する２個の媒質間を分離する表面である界面に配置されている。本発明の別の特徴は、点源の層を物体（ＯＶＯまたはＣＶＯのいずれであろうと）の表面だけに関連付けることである。本発明は、以下の２種類の境界条件を考慮に入れる。
１．既知の特性を有する２個の媒質間の界面を越える際にスカラーおよび／またはベクトル量の連続性を表わす固有境界条件（ＩＢＣ）。本発明の方法の基本原理によれば、これらのＩＢＣは、表面のいずれかの側に位置しているソースの集合により満たされる。大多数の場合において、ユーザはこれらＩＢＣを気にする必要がない。
２．界面上のスカラーおよび／またはベクトル量についてユーザが有する先見的知識を表わす、ユーザが指定した境界条件（ＵＢＣ）。本発明の方法の基本原理によれば、これらのＵＢＣは、界面の一方の側に位置しているソースの単一の集合により満たされる。これらの境界条件が系内にエネルギー源が存在することを表わすことに注意されたい。例えばコンデンサ内にあるそのようなＵＢＣは、電極間のスカラー電位差であり、ユーザにより指定されるため既知である場合が多い。
【００５０】
全ての界面について境界条件が決定されたならば、プログラムは、その実行に必要な全てのデータを受信した後で、ステップＥ８へ進んで大域行列を記述する。
【００５１】
ステップＥ８は、最初に、マスキングとして知られる動作を必要とする界面間の結合行列を決定することから始まる。マスキングは、図２〜８を参照しつつ以下に記述するように、系内のいくつかの物体が互いに遮蔽し合う（射影）する場合があるという事実を表わす。物体を互いにマスキングするこの動作は試験機能の適用を必要とする。後者は、行列を充填する間に、あるいはマスキング行列（０または１で埋められた）を構成するために別個に適用することができる。元の結合行列が試験なしで計算された状態でこのマスクを一要素毎に乗算することにより、アルゴリズムで効果的に使用される行列を得ることが可能になる。
【００５２】
どの結合行列を計算すべきかを決定するために、以下に述べる試験手順が確立されている。全ての結合行列が計算されたならば、大域行列をブロック編成する。媒質画定ステップにおいて、主情報、すなわち媒質と物体の相対位置と共にそれらの特性を含む行列が、一致指数を用いて生成されている。この行列から出発して、物体を結合可能にする試験を用い、および決定すべきＢＣの種類に応じて、大域行列は、全ての相互作用を系統的に考慮するような仕方でブロックにより記述される。ソース群を相互作用の順にソートすることにより、ブロックは原理的に対角線上およびその周辺に配置される。
【００５３】
先のステップで生成された物体対媒質の一致指数により、２個以上の物体が共通の媒質を有するか否かを容易に調べることができる。次いで、同一媒質内で相互作用する物体の結合行列が記述される。例えば、共通の相互作用媒質αを有する２個の物体ｉ、ｊに対して、得られた行列が定式化：
【数８】
2009521031.tif 000010
に応じて記述される。これは、ｊのソースが媒質β内にある状態で共通の媒質α内における物体ｉの物体ｊとの相互作用として読まれる。同様に、媒質β内にあって媒質α内へ放射している物体ｊに関連付けられたソースをＡｊα＿βと表記する。
【００５４】
同じ共通の媒質を有する全ての結合行列ＭおよびＱは、大域行列内で同符号でなければならない。
【００５５】
ＩＢＣを有する各々の界面に対して、スカラー量（行列Ｑ）と、ベクトル量（行列Ｍ）の法線成分との連続性が、少なくともこれらの例において、保障されている。例えば、静電気の場合、これらの量は場の電位および法線成分
【数９】
2009521031.tif 000011
であり、超音波の場合は圧力および速度の法線成分である。
【００５６】
ＵＢＣを有する各々の界面に対して、ユーザにより表面上の境界条件（例えば、所与の電圧が供給される電極）が課されているため、単一の行列方程式（行列ＭまたはＱ）が適用される。
【００５７】
大域行列が記述されたならば、プログラムは行列方程式を解く（ステップＥ９）。このステップは本質的に大域行列の逆行列を求めることからなる。
【００５８】
ある種の量（例えばコンデンサの容量）を計算するためにはＵＢＣの数値が既知である必要がないという事実により、従ってこれらは逆行列計算の後で直接得ることができ、プログラムがこれを直接実行する（ステップＥ１０）。そのような量が存在しない場合、あるいはそのような量に以外に、他の物理量の計算が必要である場合、ユーザが介入してプログラムに欠落データを供給する（ステップＥ１１）。
【００５９】
相互作用のモデリングに必要な物理量を計算するための欠落データは、ステップＥ７で指定された境界条件（例えば金属電極上の定常ポテンシャルの値）の種類に応じてユーザから要求される。
【００６０】
この欠落データをユーザにより供給されたならば、次いで、表示すべき観測場を画定するステップ（ステップＥ１２）へ進む。計算および観測に対するユーザの選択の自由度は、ユーザがステップＥ６を含めるか否かに依存する。このステップを行なう間作用空間の境界が画定されていない（相互作用を受け入れている系はＣＶＯしか含んでいない）場合、表示の領域に何らの制約がない。逆の場合には、ユーザは自動的に境界が画定される作用空間を画定しなければならない。
【００６１】
その後、プログラムはユーザが与えた空間をメッシュモデル化すること（ステップＥ１３）から始まる。次いで、結合行列の計算に用いたのと同じ手順に従い、今回は表示空間内の位置に適用して、調査対象の相互作用に関わる物理量を計算する。最後に、プログラムは系内の相互作用現象のグラフ表現をトレースする。図１に示す横方向の矢印は、ユーザが構造に戻ってこれを修正し、得られた新たな結果を観測できるようにすべくプログラムがループバック可能であることを示す。
【００６２】
可能な変更は、ステップＥ１０およびＥ１３から出て、ステップＥ２、Ｅ３、Ｅ５、Ｅ６、Ｅ７、Ｅ１１およびＥ１２へ向かう矢印により図１に示す定義（物体、媒質、境界条件等）に関する。これらの変更は、（ＯＶＯおよび／またはＯＶＦの）変位、または研究対象である物体または媒質の特性の変化による影響を許す。
【００６３】
図２〜８を参照しつつ以下に記述する３通りの例は、本発明による大域行列の記述技術を例示する。
【００６４】
図２、３において平面多層系１の場合を示す。図２は極めて簡素化された斜視図であり、図３は図２でＩＩＩ〜ＩＩＩに沿った断面図である。
【００６５】
図２、３に示す系１は、六面が作用空間の境界である直方体の一般的な形状を有している。これは、各々が媒質１〜媒質５として参照される５個の媒質、および各々がａ〜ｆで参照される６個の物体を含むマルチドメイン構造を有している。例えば、媒質１と媒質２との間の界面をｅで、媒質２と媒質３との間の界面をｂで示す。ここではコンデンサ電極である２個のＣＶＯａおよびｄは、作用空間の下面および上面と接触している。媒質２は、界面ｆと接触している第１の包含物ｂおよび第２の包含物ｃを含んでいる。媒質０は、未知およびユーザの関心対象外である作用空間の外側を画定する。
【００６６】
６個の異なる界面が以下のように参照される。
ａ：媒質０と１の間
ｂ：媒質２と３の間
ｃ：媒質２と４の間
ｄ：媒質５と０の間
ｅ：媒質１と２の間
ｆ：媒質２と５の間
【００６７】
電極ａに関するＵＢＣをＰａ（この電極により界面ａで課されるポテンシャル）と記述し、電極ｄにより課されるＵＢＣをＰｄ（界面ｄにおいて電極により印加されるポテンシャル）と記述する。
【００６８】
界面ａ〜ｆの各々で生成されるソース・ベクトルは以下のように記述される。
界面ａ：Ａａ１＿０１（電極ａ側）
界面ｂ：Ａｂ３＿２（媒質２内）およびＡｂ２＿３（媒質３内）
界面ｃ：Ａｃ４＿２（媒質２内）およびＡｃ２＿４（媒質４内）
界面ｄ：Ａｄ５＿０２（電極ｄ側）
界面ｅ：Ａｅ２＿１（媒質１内）およびＡｅ１＿２（媒質２内）
界面ｆ：Ａｆ５＿２（媒質２内）およびＡｆ２＿５（媒質５内）
【００６９】
上示の事例によれば、例えば界面ｃに関する一致指数をｃ（４，２）と表記する。以下の例において、一致指数は、２個の縦長の括弧の間に現れて、決定すべき非ゼロ結合行列の存在を許す。
【数１０】
2009521031.tif 000012
【００７０】
従って物体２の境界条件は以下の通りである。
【数１１】
2009521031.tif 000013
および対応する行列方程式は以下の通りである。
【数１２】
2009521031.tif 000014
【００７１】
図４（簡素化された斜視図）および５（図４のＶ〜Ｖに沿った断面図）において閉鎖系２の例を示す。
【００７２】
系２は、断面が矩形であって、本明細書で物体ａと表記されていて媒質１および媒質２（物体ｂで境界が画定される）を含んでいる中空管３からなる。電極４、５が物体ａの２個の対向面に配置されている。図５の横断面に、点源の各種の層を示す。すなわち、物体ａの外面上のＡａ１＿０、媒質１および媒質２の各々における媒質１と２の間の界面のいずれかの側のＡｂ２−１およびＡｂ１＿２である。電極４、５の配置位置において、点源は能動面の存在を示す一次ソースであるため、他と区別して示していることに注意されたい。境界条件ベクトルＰａは従って、ＵＢＣおよびＩＢＣを含んでいる。
【００７３】
系２の結合行列および一致指数は以下のように記述される。
【数１３】
2009521031.tif 000015
【００７４】
従って境界条件は以下のように記述される。
【数１４】
2009521031.tif 000016
【００７５】
最後の２例を図６〜８に示す。これらは容量が可変または調整可能なマイクロ・コンデンサの２個の変形実施形態であり、各々が、例えば「ＲＦ用途のＶ字微細構造調整可能コンデンサ」（Ａ．Ｃｒｕａｕ他、ＤＴＩＰ２００４、モントルー、および「ＲＦ用途の可動誘電体を備えた高品質の調整可能微細構造コンデンサ」（Ｊ−Ｂ．Ｙｏｏｎ他、国際ＥＤＭ２０００、サンフランシスコ）に記述されている。本発明の方法の応用に関係する特徴のみを本明細書に示す。
【００７６】
図６に示す第１の変形は、各々が平面に配置された３個の電極の４組を含んでいる。これらの４個の平面は全て共通の平面（図示せず）に対して垂直であれ、この平面との交線は、一方がもう一方より僅かに上方に配置されている２個の「Ｖ字」を形成する。電極の各組は、３個の同一の電極を、長軸が互いに平行および前記共通平面に平行な矩形のリボンの形式で含んでいる。２組の内側の「Ｖ字」のうち１組の電極は、この「Ｖ次」のもう一方の組に対向しており、同じことが外側の「Ｖ字」の電極の組にあてはまる。内側「Ｖ字」の電極の平面は、外側「Ｖ字」の平面の電極に対して関して可動であるが、各々の「Ｖ字」内では、２組の電極は互いに対して固定されている。
【００７７】
電極の組の１個（図面に示す左側の組）を物体ａ、内側「Ｖ字」の電極の２個の組を物体ｂ、外側「Ｖ字」の電極のもう一方の組を物体ｃと表記している。当然ながら、各々の電極上の境界条件は、これらの電極に印加される各々のポテンシャルにより指定される。物体ａ、ｂ、およびｃに各々印加されるポテンシャルはＰａ、ＰｂおよびＰｃと表記されている。
【００７８】
図７に示す変形において、可変コンデンサは、各々が平面に配置された２組の電極および１組の誘電体棒を含んでいて、これらの３個の平面は互いに平行である。内側平面の誘電体棒の組は、互いに対して固定された他の２組に対して可動である。外側平面（図面に示す低い方の平面）のうち１個の電極の組を物体ａ、中央平面の電極の組を物体ｂ、もう一方の外側平面の電極の組を物体ｃと各々表記している。これらの要素の全てが媒質１により囲まれていて、誘電体は媒質２を含んでいる。図６の変形の場合と同様に、図７の装置の電極上の境界条件は、電極に印加されるポテンシャルにより固定される。物体ａおよびｃに印加されるポテンシャルはまた各々Ｐａ、Ｐｃと表記される。
【００７９】
各種物体ａ、ｂおよびｃ上の各種点源を図８に示す。
【００８０】
第１の変形（図６）に対する境界条件は以下のように記述される。
【数１５】
2009521031.tif 000017
これは、次の結合行列を与え、
【数１６】
2009521031.tif 000018
一方、第２の変形（図７）に対する一致指数は以下のように記述される。
【数１７】
2009521031.tif 000019
【００８１】
境界条件は以下のように記述される。
【数１８】
2009521031.tif 000020
【００８２】
これらのＢＣにおいて、０１および０２は電極の内部を表わし、作用空間の外側にあることに注意されたい。
【００８３】
この第２の変形に対して結合行列は以下の通りである。
【数１９】
2009521031.tif 000021
【００８４】
物体ｂが物体ａとｃの間で部分的に遮蔽物として機能するため、最初の場合は行列Ｑａｃが、この最後の場合は行列Ｍａｃ（Ｑａｃ）およびＭｃａ（Ｑｃａ）が部分的にマスキングされることに注意されたい。
【００８５】
本発明によれば、複数の物体を含む系の場合、この系の大域行列を確定することができ、従ってソース・ベクトル
【数２０】
2009521031.tif 000022
の全ての成分を見つけることができる。系内の相互作用を表わす物理量は、系の観測のあらゆる時点で計算することができる。電磁気に対する応用の場合、メッシュ要素毎のソースの未知パラメータの数は、利用可能な方程式の数（ここでは６個）に等しいことに注意されたい。
【００８６】
例えば、図２の場合、記述されている系は、物体ｄにより生成された励起場の前に配置された層状の媒質１，２および５からなる導電物体であってよく、物体ｂまたはｃは亀裂または包含物を表わすことができる。
【００８７】
ゲージに関する注記：マクスウェル方程式が積分微分方程式であるため、
【数２１】
2009521031.tif 000023
から場
【数２２】
2009521031.tif 000024
がゲージとして知られる定数内で決定される。カールの発散が常にゼロであるため、
【数２３】
2009521031.tif 000025
が得られる。関数Θの勾配のカールが常にゼロであるため、
【数２４】
2009521031.tif 000026
勾配項は追加してもよい。従って場
【数２５】
2009521031.tif 000027
は以下のようになる。
【数２６】
2009521031.tif 000028
【数２７】
2009521031.tif 000029
の分散は任意に選択できることが示されている。次いでゲージとして知られる条件により
【数２８】
2009521031.tif 000030
を完全に決定することができる。クーロンゲージまたはローレンツ・ゲージとして知られるゲージが広く用いられている。ポテンシャルに対する他の条件も共存する。
【００８８】
ＤＰＳＭ法の場合、ユーザが特定のゲージ状態の検証を望んだ場合、追加的な点源を追加することにより系を完全に特徴付けて、ポテンシャルのあらゆる不確実性も除去することができる。
【００８９】
本発明による非均一媒質のモデリング方法の基本原理について、図９〜１９を参照しつつ以下に記述する。
【００９０】
非均一なアセンブリのいずれかの側に配置されている２個のトランスデューサＴおよびＳを含んでいる、換言すれば複数の物体および異なる媒質を含んでいる系を図９に示す。最初に、これらの媒質が異なる特性および／または材料（隣接層間で異なるが、各々の層は構成および特性に関して均一であると仮定する）を有する平面層であると仮定しており、従って系６は多層系として分類されている。上で述べたように、これらのトランスデューサはこの場合、超音波プローブである。この系６はＣ１〜Ｃｎで示すｎ平面の異なる媒質の層で構成されていて、Ｉ１〜Ｉｎ−１で示す（ｎ−１）個の界面により分離されていると仮定する。説明を簡明にすべく、トランスデューサを、能動面（系６の反対側の）が平面であって界面と平行であるとして示す。しかし、これらの能動面が界面と平行でない場合でも本発明が実装可能できる点は明確に理解されよう。一次ソースと呼ばれるソース層がトランスデューサをモデル化しており、界面は第２のソースと呼ばれるソースの２層によりモデル化されている。
【００９１】
図９の構成から出発して、ＤＰＳＭの拡張を本発明による多層系に実装することにより、２個のトランスデューサＳおよびＴをモデル化する２層の一次ソースＡＳ、ＡＴを、（ｎ−１）個の界面に関係する二次ソースの２（ｎ−１）個の層と共に図１０に示す。これらのソースは、対応する界面に接する小さい円の形で表わされている。各々の界面のいずれかの側に、二次ソースの層が当該界面の全面にわたり分布している。このように、例えば第１の界面Ｉ１の場合、ソースＡ＊１がトランスデューサＳに対向する界面の表面上に分布し、もう一方の面にソースＡ１が分布していて、同様にランクｎの界面まで、ランク（ｎ−１）の界面と反対側の面上にＡ＊ｎ−１、およびトランスデューサＴと反対側の面上にソースＡｎ−１が分布している。図を簡明にすべく、これら一次および二次ソースの一部だけを示す。ＡＳおよびＡＴは、以下でソース・ベクトルと称する点源の値を含むベクトルを表わす。ランク１〜ｎの各々の界面の各々の表面に、二次ソース・ベクトルＡｍ（図からわかるように各々の界面の上面に対して）およびＡ＊ｍ（各々の界面の下面の場合）のアセンブリを示す。
【００９２】
系６の各平面媒質内の超音波場全体が、対象とする媒質の直上および直下にあるソースの２個のアセンブリにより生成されたベクトルの重ね合わせにより得られる。従って、この全ベクトルは以下の方法で各種の媒質について計算される。
媒質Ｃ１：ＡＳおよびＡ１により生成された場の合計
媒質Ｃ２：Ａ＊１およびＡ２により生成された場の合計
媒質Ｃ３：Ａ＊２およびＡ３により生成された場の合計
・・・
媒質Ｃ（ｎ−１）：Ａ＊ｎ−２およびＡｎ−１により生成された場の合計
媒質Ｃｎ：Ａ＊ｎ−１およびＡＴにより生成された場の合計
【００９３】
ソース・ベクトルの値を決定するために満たされなければならない、界面における制限的な境界条件について以下で調べる。トランスデューサＳおよびＴの能動面上で、ユーザは各点において、ＶＳＯおよびＶＴＯと呼ばれる速度ベクトルに各々グループ化される（およびＵＢＣである）速度値を指定する。（ｎ−１）個の界面を越える際に、圧力Ｐおよび界面の平面と局所的に直交する方向への速度Ｖの成分（以下単に速度Ｖと称する）は、超音波物理（ＩＢＣ）から生じる理由により連続でなければならない。同じことが他の分野の物理量の連続性の条件にもあてはまり、その場合で圧力および速度ではない他のパラメータを量ＰおよびＶで表わす。
【００９４】
項Ｖ、Ｐは、関係式Ｖ＝Ｍ．ＡおよびＰ＝Ｑ．Ａを介して基本ソース・ベクトルに紐付けられている。ここに、ＭおよびＱは同じ媒質内で相互作用している物体の結合行列である。行列Ｑは、スカラー量（超音波波の場合は音響圧、静電気の場合は電気スカラー・ポテンシャル、あるいは静磁気の場合は磁気スカラー・ポテンシャル）を特徴付けるのに対し、行列Ｍはベクトル量（本例では超音波波の伝播速度、静電気の電場、または静磁気の磁場）を特徴付ける。境界条件および連続性の条件から次式が得られる。
【数２９】
2009521031.tif 000031
これらの式は、添え字（ｎ−ｊ）内の値をｎｊ、ｊ＝１、２、３の形式で以下のように書くことにより行列形式で表わすことができる。
【数３０】
2009521031.tif 000032
この行列は、記号的な形式で書くこともできる。
【数３１】
2009521031.tif 000033
従って、一次および二次ソースの全てをグループ化して、大域ソースと呼ばれるベクトル｛Ａ｝が得られる。
【００９５】
トランスデューサの能動面が、対向する物体の端部界面上に配置される場合を以下で調べる。異なる媒質の（ｎ＋２）層を分離する（ｎ＋１）個の界面（Ｉ０、Ｉ１、・・・Ｉｎ）を考える。これらの層は（０、１、２、・・・、ｎ、ｎ＋１）と番号付けられたランクを有する。トランスデューサＳおよびＴの能動面は、図１１に示すように、物体１の下部（Ｉ０）および上部（Ｉｎ）界面上に各々配置される。これら２個の界面は各々、下半分の空間および上半分の空間の境界を画定する。ランク０の媒質が界面Ｉ０と接触しており、一方、ランクｎ＋１の媒質は界面Ｉｎと接触してある。本発明の方法に関して上で述べたように、ソースの層はｎ＋１個の界面Ｉ０、Ｉ１．．．Ｉｎのいずれかの側に分布している。ソース・ベクトルのこれら２層は、Ａｍ（界面上方のソースの場合）およびＡ＊ｍ（界面下方のソースの場合）と表記されており、ここに０＜ｍ＜ｎである。しかし、最初（Ｉ０）と最後（Ｉｎ）界面は、３個の異なるソースの層と隣接している。実際、追加的なソースＡＳ、ＡＴの２個のアセンブリは、トランスデューサの能動面がこれらの最初および最後の界面に対して適用されると仮定されていることに起因している。界面Ｉ０はトランスデューサＳの能動面に一致しているが、本発明の方法の図示の便宜上、ここではＩ０とＳを２個の異なる界面として考える。Ｉ０は受動領域（エネルギー源無し）を表わし、Ｓは表能動領域（１個のエネルギー源が存在する）を表わす。同様に、ＩｎとＴを２個の異なる界面と考える。このように、ここに示す場合において、いずれかの側にソースの１個のアセンブリが配置されているｎ＋１個の受動界面Ｉ０、Ｉ１、．．．Ｉｎ、および各々が追加的なソースＡＳ、ＡＴの２個のアセンブリに対応する２個の能動面ＳおよびＴだけが存在する。
【００９６】
２個の端部媒質を除いて、全ての媒質における超音波場は、以下に注記されるように、ソースの２個のアセンブリにより場が生成されると仮定することにより得られる。
媒質２：Ａ＊１およびＡ２により生成された場の合計
媒質３：Ａ＊２およびＡ３により生成された場の合計
．．．．．．．．．．
．．．．．．．．．．
媒質ｎ−１：Ａ＊ｎ−２およびＡｎ−１により生成された場の合計
【００９７】
しかし、界面Ｉｎ上方および界面Ｉ０下方の各半空間内で、超音波場はソースの単一層によってのみ生成される。
媒質０：場はソースＡ０により生成される。
媒質ｎ＋１：場はソースＡ＊ｎにより生成される。
【００９８】
上側および下側の媒質（ランク１およびｎの媒質）において、場は３個のソース層により生成される。
媒質１：Ａ＊０、Ａ１およびＡＳにより生成された場の合計
媒質ｎ：Ａ＊ｎ−１、ＡｎおよびＡＴにより生成された場の合計
【００９９】
界面における制限的な境界条件から出発してソースの値を決定するために、以下の界面における条件が最初に満たされなければならない。ｎ−１個の中間的な受動界面（Ｉ１、Ｉ２、．．．Ｉｎ）を越える際に、界面の平面の局所的に法線方向にある超音波圧力Ｐおよび速度Ｖは連続でなければならない。
【０１００】
能動面Ｓ、Ｔが除外されている場合、界面Ｉ０、Ｉｎは受動領域とのみ隣接している。これらの受動界面を越える際に、他のｎ−１個の受動界面を越える場合と全く同様に、圧力Ｐおよび速度Ｖは連続でなければならない。
【０１０１】
トランスデューサの表面Ｓ、Ｔ上の速度ベクトルＶは各々ＶＳ０およびＶＴ０と表記されている。上で述べたように、Ｖ＝Ｍ．ＡおよびＰ＝Ｑ．Ａとすれば、連続性の条件および境界条件（ＵＢＣおよびＩＢＣ）から次式が得られる。
【数３２】
2009521031.tif 000034
【０１０２】
これらの方程式は、以下のように再構成できる。
【数３３】
2009521031.tif 000035
【０１０３】
行列表現で、（ｎ−ｊ）を形式ｎｊ、ｊ＝１、２、３・・・で表わすことにより、次式が得られる。
【数３４】
2009521031.tif 000036
【０１０４】
方程式の数および未知数の数は、以下のように数えられる。各々の界面Ｉ１、Ｉ２、．．．Ｉｎ−１に対して、界面の上方にＮ個のソースがあり、界面の下方にＮ個のソースがあって、下位Ｉ０個および上位Ｉｎ＋１個の界面において、トランスデューサＳおよびＴの２個の能動面上に各々Ｍ１およびＭ２個のソースがあり、これらの２個の界面の受動領域上に各々２（Ｎ−Ｍ１）個および２（Ｎ−Ｍ２）個のソースであると仮定する。従って、ｎ＋１個の界面（トランスデューサの能動面を含む）に対するソースの総数は、２（ｎ＋１）Ｎ−Ｍ１−Ｍ２に等しい。方程式の総数も２（ｎ＋１）Ｎ−Ｍ１−Ｍ２に等しい。その結果、この方程式系は一意な解を有する。
【０１０５】
行列方程式（６）は、以下のように記号的な形式で記述することができる。
【数３５】
2009521031.tif 000037
行列Ｍの次元が、
｛２（ｎ＋１）Ｎ−Ｍ１−Ｍ２．＊２（ｎ＋１）Ｎ−Ｍ１−Ｍ２．｝
【０１０６】
トランスデューサの能動面が界面の一部を形成する場合について以下で調べる。ここで相互作用のモデリングは上で述べたものとは僅かに異なる方法で行なわれる。先の場合とは対照的に、トランスデューサの一次ソースは、図１２に示すように界面二次ソースのいくつかで置換される。ソースの層が、ｎ＋１個の界面Ｉ０、Ｉ１、．．．Ｉｎのいずれかの側に配置されている。界面Ｉｍの一方の面のソース・ベクトルのうち、界面上方のものをＡｍと表記し、界面下方のものをＡ＊ｍと表記する。図１２において、トランスデューサＳ、Ｔを、ｎ個の異なる媒質を含む物体の界面Ｉ０およびＩｎと接触している状態で示す。これらの２個の界面は、媒質０およびｎ＋１の側に位置する半空間の境界を画定する。図１１と１２の違いは、界面Ｉ０とＩｎの近傍における点源の分布にある。
【０１０７】
各媒質内の超音波場の総計は、ソースの２個のアセンブリにより生成された場の重ね合わせにより以下のように得られる。
媒質１：Ａ＊０とＡ１により生成された場の合計
媒質２：Ａ＊１とＡ２により生成された場の合計
媒質３：Ａ＊２とＡ３により生成された場の合計
媒質ｎ：Ａ＊ｎ−１とＡｎにより生成された場の合計
しかし、上述の各半空間において、場はソースの単一の層により生成される。
媒質０：場はソースＡ０により生成される。
媒質ｎ＋１：場はソースＡ＊ｎにより生成される。
【０１０８】
境界条件から出発してソース・ベクトルの値を決定するために、界面上において以下の条件が満たされなければならない。ｎ−１個の受動界面（Ｉ１、Ｉ２、．．．Ｉｎ−１）を越える際に、対象とする界面に対して局所的に法線方向にある圧力Ｐおよび速度Ｖは、一切の不連続性を示してはならない。界面Ｉ０、Ｉｎは、能動領域および受動領域に隣接している。トランスデューサＳ、Ｔの能動面は能動領域に配置されている。これらの領域を各々
【数３６】
2009521031.tif 000038
と表記する。トランスデューサの能動面を含んでいない受動領域をＩ０、Ｉｎ（指数無しの）と表記する。
【０１０９】
トランスデューサの能動面
【数３７】
2009521031.tif 000039
上における速度ベクトルを各々ＶＳ０、ＶＴ０と表記し、受動界面Ｉ０およびＩｎと交差する際に、対象とする界面に対して局所的に法線方向にある圧力Ｐおよび速度Ｖは、他のｎ−１個の受動界面と全く同様に、一切の不連続性を示してはならない。
【０１１０】
Ｖ＝ＭＡおよびＰ＝ＱＡであることから、界面における制限的な境界条件から以下の方程式が得られる。
【数３８】
2009521031.tif 000040
【０１１１】
これらの方程式は、（ｎ−ｊ）をｎｊ、ｊ＝１、２、３．．．の形式に書くことにより、行列形式で記述することができる。
【数３９】
2009521031.tif 000041
ランク０およびｎの界面に対して中カッコ
【数４０】
2009521031.tif 000042
で囲まれた項がゼロで置換されるのに対し、能動領域
【数４１】
2009521031.tif 000043
およびランク０およびｎの界面に対して山形カッコ
【数４２】
2009521031.tif 000044
で囲まれた項がゼロで置換される。
上述の行列（９）は以下の記号的な形式で記述することができる。
【数４３】
2009521031.tif 000045
従って、大域ソース・ベクトル｛Ａ｝を得ることができる。
【０１１２】
単一界面および単一トランスデューサの特定の場合について図１３を参照しつつ調べる。この場合は、図１０に示す場合においてトランスデューサＡＴおよびＩ１を除く全ての界面が後者から除外されている特定の場合と見なすことができる。この場合、反射係数ＲをＭとＱの値に関して計算することができ、それ故に、行列ＭおよびＱの物理的重要性に関する概念を有している。この場合方程式（２）は以下のように簡約される。
【数４４】
2009521031.tif 000046
（１１）の２番目と３番目の方程式から以下が得られる。
【数４５】
2009521031.tif 000047
（１２）の２個の方程式を減算することにより次式が得られる。
【数４６】
2009521031.tif 000048
但し、
【数４７】
2009521031.tif 000049
ここから、（１１）に従い、
【数４８】
2009521031.tif 000050
以下を記述することができる。
【数４９】
2009521031.tif 000051
媒質２内の界面における速度ベクトルは、以下のように透過係数Ｔの行列形式で表わすことができる。
ＴＭ１ＳＡＳ
また、媒質１における反射係数Ｒの行列形式で以下のように表わされる。
Ｍ１ＳＡＳ＋ＲＭ１ＳＡＳ
これにより以下が得られる。
【数５０】
2009521031.tif 000052
次いで以下が得られる。
【数５１】
2009521031.tif 000053
当該問題が反射および透過係数行列に関して定式化されている場合、前記行列が考慮されているならば速度の連続性が圧力の連続性を保証するため、圧力の連続性の条件は必ずしも個別に満足される訳ではないことに注意されたい。
【０１１３】
方程式（１３）および（１７）から出発して、以下を記述することができる。
【数５２】
2009521031.tif 000054
【０１１４】
ＤＰＳＭの教示から出発して、トランスデューサの表面に分布しているＮ個の点源の値を例えば以下のように計算する。ここで、トランスデューサは流体と接触状態にあると仮定しているが、トランスデューサが固体と接触している場合にも本記述を容易に適用できることが明確に理解されよう。Ａｍをｍ番目の点源の値とする。この点源から距離ｒｍにおける圧力ｐｍは次式で示される。
【数５３】
2009521031.tif 000055
上式でｉ２＝−１、ｋｆ＝ω／ｃｆであり、ωは超音波のパルス周波数、ｃｆは波の伝播速度である。ｍ番目の点源から距離ｒｍにある点のｘにおける超音波圧力ｐ（ｘ）は次式で示される。
【数５４】
2009521031.tif 000056
波の超音波圧力と伝播速度との既存の関係から出発すれば、点ｘにおける波の速度は次式により計算することができる。
【数５５】
2009521031.tif 000057
また、時間の関数としてｅ−ｉωＴの形式で表わされる速度の場合、その導関数はｖｎに正または負の項を単に乗算することで得られる。この項が負の場合には次式が得られる。
【数５６】
2009521031.tif 000058
従って、ｍ番目の点源から距離ｒにおける超音波波の動径速度ｖｎの値は次式で得られる。
【数５７】
2009521031.tif 000059
方向ｘ３における速度成分は次式で得られる。
【数５８】
2009521031.tif 000060
Ｎ個の点源の寄与が加算される場合、方向ｘ３における速度の総和は次式で得られる。
【数５９】
2009521031.tif 000061
ｖ０が方向ｘ３におけるトランスデューサの能動面上の波の速度である場合、この能動面上のｘの全ての値について、方向ｘ３における速度はｖ０に等しく、従って以下が得られる。
【数６０】
2009521031.tif 000062
トランスデューサの能動面にＮ個の点を取ることにより、Ｎ個の未知数（Ａ１、Ａ２、Ａ３、．．．ＡＮ）を有するＮ個の一次方程式系が得られる。しかし、点源が対象点ｘに一致する場合、ｒｍがゼロになって方程式（１８−Ｈ）のｖ３ｍが無制限になるため、問題が生じる。当然ながら、点源と対象点ｘがトランスデューサの能動面に配置されている場合、これら２個の点だけが一致してｒｍがゼロになり得る。この難点を除去すべく、点源はトランスデューサの能動面の僅かに後ろ側に配置される。このような場合、ｒｍの最小値をｒｓと表記する。
【０１１５】
点ｘがトランスデューサの能動面上に配置されている場合、その速度成分ｘ３と当該能動面上の波の速度ｖ０との相関が計算される。次いで、他の２個の速度成分ｖ１およびｖ２の相殺が望ましい場合、３個の方程式、すなわち上記の方程式（１８−Ｈ）および以下の２個の方程式（１８−Ｉ）が同時に満たされなければならない。
【数６１】
2009521031.tif 000063
従って、前記仏国ＤＰＳＭ特許に示すように、球面上のＮ個の試験点に対して３Ｎ個の方程式が得られる。しかし、未知数より多くの方程式が得られる。この問題を解決するために、当該仏国特許は、各々の点源を点源の三重項で置換することにより、３Ｎ個の未知数を得るために未知数の数を増やすことを推奨している。これら３個のソースは異なる値を有し、トランスデューサの能動面に平行な同一平面上に、能動面から距離ｘ３＝−ｒｓの位置に配置されている。３個の点源は、波が伝播する材料の等方性を維持して、放射の選択配向を避けるために、ランダムな向きの正三角形の角に配置されている。従って、３Ｎ個の一次方程式（三重項の場合）系またはＮ個の一次方程式（単一点源の場合）系を解くことにより、全ての点源に関連付けられたＡｍの値が得られる。Ａｍの値が得られたならば、方程式（１８Ｂ）を用いて、トランスデューサの能動面上または向こう側に拘らず全ての点での圧力ｐ（ｘ）が計算される。非粘性完全流体の場合、流体と固体の界面（トランスデューサの能動面上の）の法線方向にある速度（ｖ３）成分だけが連続でなければならないことに注意されたい。能動面と流体との間に滑りが生じる可能性があるため、この能動面と平行な速度成分は必ずしも連続である訳ではない。しかし、粘性流体の場合にはそのよう滑りが生じることはなく、前記界面で速度の３成分は連続でなければならない。
【０１１６】
流体内における超音波波の伝播に関して、行列解析を用いたソースのパラメータ計算を以下に示すが、この計算が本明細書で述べたような他の応用分野に容易に適用できることが明確に理解されよう。この計算は、トランスデューサの能動面上の速度の３個の直交座標成分が、粘性流体の場合と同様にトランスデューサとの界面における流体のそれらに一致する場合の三重項ソースについて示すものである。逆の場合には、計算における速度成分ｖ１およびｖ２を除外することにより、三重項の代わりに単一点源が用いられる。次いで行列の次元およびベクトルの次元が３ＮからＮまで減らされる。
【０１１７】
方程式（１８−Ｈ）と（１８−Ｉ）を組み合わせて次式が得られる。
【数６２】
2009521031.tif 000064
この方程式において、ＶｓはＮ個の点ｘの速度成分の（３Ｎ×１）次元ベクトルであり、Ａｓは３Ｎ個の点源に関する物理量を含む（３Ｎ×１）次元ベクトルである。Ｍｓｓは、２個のベクトルＶｓおよびＡｓを結び付けている（３Ｎ×３Ｎ）次元の行列である。方程式（１８−Ｈ）および（１８−Ｉ）から次式も得られる。
【数６３】
2009521031.tif 000065
列ベクトルＶｓの転置は（１×３Ｎ）次元の行ベクトルであることに注意されたい。このベクトルの要素を
【数６４】
2009521031.tif 000066
表記する。ここに、添え字ｊの値は１、２または３であって、波の伝播速度の直交座標成分の方向を示す。指数ｎは１〜Ｎの間の任意の値を取り得て、速度成分が定義されるトランスデューサの能動面上の点に対応している。
【０１１８】
大多数の超音波トランスデューサの場合、ｊ＝１または２ならば（換言すれば界面に平行な速度成分について）
【数６５】
2009521031.tif 000067
は０、およびｊ＝３ならば（速度の法線成分）
【数６６】
2009521031.tif 000068
である。この場合、方程式（１８−Ｋ）は次式となる。
【数６７】
2009521031.tif 000069
および、ベクトルＡｓは次式で示される。
【数６８】
2009521031.tif 000070
方程式（１８−Ｈ）、（１８−Ｉ）から出発して正方行列ＭＳＳが得られる。
【数６９】
2009521031.tif 000071
但し
【数７０】
2009521031.tif 000072
方程式（１８−Ｏ）において、ｘの第１の添え字ｊは、値１、２または３を取ることができ、ｘがどの方向（ｘ１、ｘ２またはｘ３）に関係するかを示す。ｘおよびｒの添え字ｍは考慮している点源に応じて１〜３Ｎの値を取ることができるのに対し、添え字ｎは、速度成分を計算しているトランスデューサの能動面の対象点に応じて、１〜Ｎの間の任意の値をとることができる。
【０１１９】
ｘｎを点ｘの（トランスデューサの能動面の背後に配置されている）１つであるとする。ｍ番目の点源をこの点に接続しているベクトルは
【数７１】
2009521031.tif 000073
と表記されていて、方向ｘ１、ｘ２およびｘ３のその３個の成分は方程式（１８−Ｎ）、（１８−Ｏ）内で値
【数７２】
2009521031.tif 000074
、ｊ＝１，２，３を取る。
【０１２０】
方程式（１８−Ｊ）を用いて、行列ＭＳＳの逆行列を求めることにより、点源に関係する量が得られる。
【数７３】
2009521031.tif 000075
【０１２１】
ソース・ベクトルＡＳがこのように計算されたならば、方程式（１８−Ｂ）圧力ｐ（ｘ）が、あるいは方程式（１８−Ｈ）および（１８−Ｉ）から全ての点（トランスデューサの能動面上、またはその向こう側）における速度ベクトルＶ（ｘ）が得られる。以下に、速度および圧力ベクトルが計算される任意の点または流体（あるいはより一般的に、波と媒質の材料との間に相互作用が生じる媒質上の任意の点）を「観測点」と称する。圧力および速度の成分は以下の関係から得られる。
【数７４】
2009521031.tif 000076
上式において、ＰＴは（Ｍ×１）次元のベクトルであってＭ個の観測点における圧力の値を含んでおり、ＶＴは（３Ｍ×１）次元ベクトルであって任意の観測点における３成分を含んでいる。ＶＴの式は、式（１８−Ｋ）から取られたＶＳの式に類似している。唯一の違いは、（３Ｎ×１）次元ではなく（３Ｍ×１）次元である事実にある。観測点が方程式（１８−Ｐ）内のソース・ベクトルＡＳを得るように速度成分が選択されているトランスデューサの能動面上の点と同一である場合、行列ＭＴＳは方程式（１８−Ｎ）におけるＭＳＳと同一である。しかし、さまざまな点で速度ベクトルを計算するために、ＭＴＳを与える式は、方程式（１８−Ｎ）から取られたＭＳＳの式とは僅かに異なる。従って、その次元は以下の方程式（１８−Ｒ）で示すように（３Ｍ×３Ｎ）である。
【数７５】
2009521031.tif 000077
上式において、項
【数７６】
2009521031.tif 000078
は方程式（１８−Ｏ）に現れるものと同一である。添え字ｊおよびｍの定義は方程式（１８−Ｏ）のものと同一である。変数ｘおよびｒの指数ｎは、対象点に応じて１〜Ｍの間の任意の値を取ることができる。ＭがＮと異なる場合、ＭＴＳが正方行列でないことに注意されたい。
【０１２２】
行列ＱＴＳは、３Ｎ個の点源およびＭ個の点を考慮している場合に方程式（１８−Ｂ）から得られる。
【数７７】
2009521031.tif 000079
この方程式において、
【数７８】
2009521031.tif 000080
の定義は方程式（１８−Ｒ）および（１８−Ｑ）におけるものと同一であって、ｍ番目の点源とｎ番目の対象点との間の距離である。
【０１２３】
方程式（１８）に戻り、方程式（１８−Ｏ）および（１８−Ｓ）から出発して、この方程式の行列Ｑ（音響圧力とソース・ベクトルの関係を確立する行列）の要素と、その行列Ｍの要素（速度およびソース・ベクトルの法線成分の関係を確立する）と共に得られ、次式を与える。
【数７９】
2009521031.tif 000081
【数８０】
2009521031.tif 000082
【０１２４】
ここで２個の行列
【数８１】
2009521031.tif 000083
を方程式（１８）の行列Ｑ＊およびＭ＊と同様に定義するが、唯一の違いは、行列Ｑ＊およびＭ＊がソース分布値Ａ１＊として媒質２内の圧力および速度を与えるため媒質２の特性に関係するのに対し、行列
【数８２】
2009521031.tif 000084
は媒質１の特性を用いることである。その結果、行列
【数８３】
2009521031.tif 000085
は、媒質２が媒質１により置換された場合に、点Ａ１＊に応じて分布している等価なソースの圧力および速度の値を与える。方程式１９は以下のように記述することができる。
【数８４】
2009521031.tif 000086
同様に以下が得られる。
【数８５】
2009521031.tif 000087
ここに、行列
【数８６】
2009521031.tif 000088
が方程式（１８）の行列Ｑ１１およびＭ１１と同様であって、唯一の違いは、
【数８７】
2009521031.tif 000089
がソース分布構成Ａ１＊の界面における位置の圧力および速度の値を与えるのに対し、行列Ｑ１１およびＭ１１はソース分布構成Ａ１の界面の位置の圧力および速度の値を与えることを注記する。これらの２個の行列が媒質１の材料特性を用いることに注意されたい。分布構成Ａ１が界面の直上に存在することから、分布構成Ａ１＊は界面の直下に存在するのに対し、界面のいずれかの側に存在する点に対する法線速度ベクトルは同じ大きさながら逆向きである。これら点源の２個の層により生成された圧力は同じ大きさである。従って以下が得られる。
【数８８】
2009521031.tif 000090
方程式（２２）を方程式（１８）に代入することにより以下が得られる。
【数８９】
2009521031.tif 000091
【０１２５】
界面の位置Ａ１＊における
【数９０】
2009521031.tif 000092
ソース分布構成が、対象点におけるソースＡＳのそれに等しい場合、ソースの２個のアセンブリの速度の式を等しくすることにより以下が得られる。
【数９１】
2009521031.tif 000093
同様に、圧力についても以下が得られる。
【数９２】
2009521031.tif 000094
【数９３】
2009521031.tif 000095
は、同一位置に存在するにもかかわらず異なることに注意されたい。
【数９４】
2009521031.tif 000096
は、２個の媒質１および２が同一である場合の点源の等しい分布を表わす。換言すれば、媒質２を媒質１で置換したならば、
【数９５】
2009521031.tif 000097
の場合には真の界面は存在しない。しかし、Ａ１＊の場合には媒質１と２の間に真の界面が存在する。
【０１２６】
方程式（２４）と（２５）を比較すれば以下が得られる。
【数９６】
2009521031.tif 000098
方程式（２６）を方程式（２３）に代入すれば以下が得られる。
【数９７】
2009521031.tif 000099
方程式（２７）値でα１、α２、β１、β２の値を方程式（２０）および（２１）で所与の値で代入すれば以下が得られる。
【数９８】
2009521031.tif 000100
波を高周波（ｋｍｒ＞＞１）と見なせば以下が得られる。
【数９９】
2009521031.tif 000101
対照的に、低周波（ｋｍｒ＜＜１）の場合は以下が得られる。
【数１００】
2009521031.tif 000102
【０１２７】
高周波での反射係数が法線入射における平面波の反射係数と同様であることに注意されたい。一方、低周波数では僅かに異なる。その理由は、ここに示す例が点源を有する球面波に関係するためである。平面波は、より多くの点源を重ね合わせることによりモデル化されている。
【０１２８】
本発明の方法の特徴について、多層系およびトランスデューサの能動面上の相互作用との関係において示す。ここで考慮する幾何形状的問題は、図１０に関する問題と同一である。今回の場合と、説明の始めに示したように境界ＳおよびＴにおける波の速度が指定されていて各々ＶＳ０およびＶＴ０と表記した場合における境界条件の違いは、今回示す場合には速度が指定されていないことである。全ての他の界面、およびトランスデューサが存在しない場合、表面ＳおよびＴ上の当該速度は各々ＶＳ０およびＶＴ０である。別の界面および第２のトランスデューサが導入されたならば直ちに、表面Ｓ上の波の速度はＶＳ０でなくなり、また表面Ｔ上の速度はＶＴ０でなくなる。トランスデューサと界面の相互作用に起因するトランスデューサ面上の速度変化について以下で考慮する。
【０１２９】
この問題は２段階で解決される。第１段階において、能動トランスデューサのソース・ベクトルＡＳ、ＡＴの値は、表面速度条件から出発して得られる。
【数１０１】
2009521031.tif 000103
第２段階において、図１４に示すようにＳおよびＴを含んでいる点源の２個の層が各々の界面で導入される。これらの新たなソースの導入により、界面Ｓの下方のソースの合計値はＡＳ＋Ａ＊０に等しく、界面Ｔの上方でＡＴ＋Ａｎに等しい。図１４において、界面ＴをＩ０で表わし、界面ＴをＩｎで表わす。その他の界面は各々Ｉ１〜Ｉｎ−１で示す。
【０１３０】
大多数の層における超音波場は、点源の２個の層により生成された場の重ね合わせにより、以下のように得られる。
媒質２：Ａ＊１およびＡ２により生成された場の合計
媒質３：Ａ＊２およびＡ３により生成された場の合計
・・・・・・・・・
・・・・・・・・・
媒質ｎ−１：Ａ＊ｎ−２およびＡｎ−１により生成された場の合計
対照的に、下側半空間（表面Ｓ上方の媒質１）および上側半空間（表面Ｔ下方の媒質ｎ）において、全場は３個の点源により生成される。
媒質１：Ａ＊０、Ａ１、およびＡＳにより生成された場の合計
媒質ｎ：Ａ＊ｎ−１、Ａｎ、およびＡＴにより生成された場の合計
トランスデューサにより（表面Ｓの下側、および表面Ｔの上側に）形成された媒質内で、場は点源の単一層により生成される。
トランスデューサＳ内部：場はＡ０ソースだけにより生成される
トランスデューサＴ内部：場はＡ＊ｎソースだけにより生成される
【０１３１】
ソースの値は、界面における境界条件から以下のように決定される。
【０１３２】
界面において以下の境界条件（ＩＢＣ）が満たされなければならない。ｎ−１個の中央受動界面（Ｉ１、Ｉ２、．．．Ｉｎ−１）を越える際に、圧力（Ｐ）および法線速度（Ｖ）、すなわち、対象とする界面に対し局所的に法線方向にある速度の成分は連続でなければならない。同様に、界面Ｓ（Ｉ０）およびＴ（Ｉｎ）を越える際に、他の全てのｎ−１個の界面の場合と同様に、対象とする界面に対し局所的に法線方向にある圧力（Ｐ）および速度の成分（Ｖ）は連続でなければならない。
【０１３３】
速度をＶ＝ＭＡおよび圧力をＰ＝ＱＡと表記すれば、連続性の条件から以下の方程式が得られる。
【数１０２】
2009521031.tif 000104
行列形式で、（ｎ−ｊ）をｎｊ、ｊ＝１、２、３．．．の形式で記述すれば、これらの方程式は以下のように記述することができる。
【数１０３】
2009521031.tif 000105
【０１３４】
方程式と未知数の数を数えるために、界面Ｉ１、Ｉ２、．．．Ｉｎ−１の各々において、対象とする界面の上方にＮ個の点源、また当該界面の下方にＮ個の点源であると仮定する。下側（Ｉ０）および上側（Ｉｎ＋１）界面に、表面Ｓの上方および表面Ｔの下方のソースの２層を含めて各々２Ｍ１個および２Ｍ２個の点源がある。従って、ｎ＋１個の界面（表面ＳおよびＴを含む）における未知ソースの総数は、２［（ｎ−１）Ｎ＋Ｍ１＋Ｍ２］に等しく、またこれらの未知ソースに加え、方程式（３１）を用いて得られたＭ１＋Ｍ２個の既知ソースＡＳおよびＡＴがある。方程式の総数もまた、対象とする界面の各々の共通点において、２個の条件を満足しなければならず、および全ての界面について考慮した共通点の総数が［（ｎ−１）Ｎ＋Ｍ１＋Ｍ２］に等しいため、２［（ｎ−１）Ｎ＋Ｍ１＋Ｍ２］に等しい。従って、この方程式系は一意な方法で解かれる。
【０１３５】
行列方程式（３３）は以下の記号的な形式で記述することができる。
【数１０４】
2009521031.tif 000106
【０１３６】
上で述べたように、行列Ｍの次元は２［（ｎ−１）Ｎ＋Ｍ１＋Ｍ２］×２［（ｎ−１）Ｎ＋Ｍ１＋Ｍ２］に等しい。未知ソース・ベクトルの数もまた、２［（ｎ−１）Ｎ＋Ｍ１＋Ｍ２］に等しい。この結果、大域ソース・ベクトルが一意な方法で得られることになる。
【０１３７】
ここで、２個のトランスデューサ間で生じる干渉現象、より具体的には２個の連続する界面間における複数反射の現象に関して図１５を参照しつつ調べる。
【０１３８】
最初に図１５に示す極めて簡単な例、すなわち対向して配置されていて均一な媒質ＭＬにより分離されている２個のトランスデューサＳおよびＴについて調べる。各トランスデューサの能動面は、能動領域と受動領域を含んでいてよい。能動領域において、界面の影響が一切ない状態での波の表面速度Ｖ０Ｓ（またはＶ０Ｔ）が指定されている。受動領域において、他のいかなるソースにも影響されない場合、表面速度はゼロである。
【０１３９】
ＡＳ０およびＡＴ０は、トランスデューサＳとＴの間の相互作用の影響が無視される場合に、トランスデューサＳおよびＴ全体にわたり分布している点源である。換言すれば、ＡＳ０およびＡＴ０は、他に一切のトランスデューサも散乱体も存在しない状態で、均一な媒質内にあるトランスデューサにより生成された場をモデル化するために用いる点源の分布を表わす。
【０１４０】
場ＡＳ０はトランスデューサＴにより反射される。この場は、トランスデューサＴ上に分布するソースＡＴ１の新しい層を導入することによりモデル化される。ソースＡＴ１により生成される表面Ｓ上の、および場ＡＳ０の反射に起因する速度ベクトルの方程式を記述することにより以下が得られる。
【数１０５】
2009521031.tif 000107
【０１４１】
ＡＴ１により生成された場は、続いてトランスデューサＳにより反射される。この場は、トランスデューサＳの表面上に分布するソースＡＳ２の層を導入することによりモデル化される。ソースＡＳ２により生成される表面Ｓ上の、および場ＡＴ１の反射に起因する速度ベクトルの方程式を記述することにより以下が得られる。
【数１０６】
2009521031.tif 000108
これらの動作を他の全ての反射について繰り返すことにより以下が得られる。
【数１０７】
2009521031.tif 000109
同様に、ソースＡＴ０から出発して、ＡＳ１から来る表面Ｓ上の場と、ＡＴ０により反射されたものを結合することにより以下が得られる。
【数１０８】
2009521031.tif 000110
および
【数１０９】
2009521031.tif 000111
この結果、トランスデューサＳおよびＴのソース・ベクトルを次式のように表わすことができる。
【数１１０】
2009521031.tif 000112
次いで、速度ベクトルおよび圧力の分布を以下の方程式から得ることができる。
【数１１１】
2009521031.tif 000113
【０１４２】
トランスデューサから発せられる波の速度が不変の場合、ソース・ベクトルは方程式（４１）から得られる。
【数１１２】
2009521031.tif 000114
【０１４３】
能動領域の表面積が、トランスデューサの総表面積と異なり、および能動領域の表面上で速度ベクトルだけが指定されている場合、次式が得られる。
【数１１３】
2009521031.tif 000115
【数１１４】
2009521031.tif 000116
は単位行列でなければならないことに注意されたい。方程式（４４）はまた、以下のように記述することもできる。
【数１１５】
2009521031.tif 000117
但し、
【数１１６】
2009521031.tif 000118
従って、
【数１１７】
2009521031.tif 000119
【０１４４】
ＡＳ０およびＡＴ０が受動領域においてゼロであることから、この方程式の第２項はゼロである。ＷＭ０ＤＷの式におけるゼロ値の列を除去して、これを行列Ｅと定義することにより、以下の簡素化された方程式が得られる。
【数１１８】
2009521031.tif 000120
【０１４５】
本発明の方法はまた、散乱の場合にも適用可能であって、特に、散乱体が任意の所与の形状を有し、図１６に示すように均一であると仮定される場合に適用可能である。そのような散乱体は、均一な媒質（媒質１と呼ぶ）、または材料（媒質２と呼ぶ）の層内の包含物であってよい。そのような場合、点源の２個の層が実装される。すなわち、散乱体内にあってその外面と接触している第１層（Ａ１）と、媒質１内で当該散乱体の外面上にある第２層（Ａ１＊）である。ソースのこれら２層に加えて、トランスデューサＳ、Ｔをモデル化する点源ＡＳ、ＡＴの２層を考慮する。媒質１内の場の合計値は、ソースＡＳ、ＡＴおよびＡ１の３層を重ね合わせることにより得られるのに対し、媒質２内の場はソースＡ１＊の単一の層に起因する。２個のトランスデューサの表面上の境界条件および対象とする２個の媒質間の界面を越える際の連続性条件を適用することにより、以下が得られる。
【数１１９】
2009521031.tif 000121
この方程式は、以下の行列形式に記述することができる。
【数１２０】
2009521031.tif 000122
この行列方程式がソース・ベクトルＡＳ、ＡＴ、Ａ１およびＡ１＊について解かれる。
【０１４６】
上に示した方法により原理的には、散乱体が存在する状態で相互作用をモデル化する問題を解決することができる。しかし、そのような包含物に起因するマスキングの問題は、解決されないままである。従って、媒質１内のどのような所与の点（以下、「対象点」と称する）も、層Ａ１の全ての点源を「見る」訳ではなく、および２個のトランスデューサの全ての点源を見ることが可能にはならないだろう。先の例（図９〜１５に示す）において、平面界面だけを含む物体の場合、この問題は、対象点が遮蔽された領域には決してなかったという事実のために存在していなかった。この困難さを回避すべく、点源は一方向、すなわち１個の半空間だけにエネルギーを放射するものと仮定されているのに対し、他の半空間はマスキングされたままであり、トランスデューサの点源により発される光線の経路が当該対象点に到達するのを包含物が妨げる場合、当該点源の寄与は無視される。この手順を、図１６と同じ参照要素を再び用いて図１７に示している。
【０１４７】
図１７において、Ｐと表記されて媒質１内に位置している対象点を、点Ｐで計算された場の生成に寄与する点源の３個の層（ＡＳ、Ａ１およびＡＴ）と共に示す。これら点源の３個の層が前の半空間だけに放射して、後の半空間がマスキングされたままである場合、これら点源の３個の層に関係する放射パターンは、図１７の小さい三日月形により表わされるようになる。点Ｐは次いで、点ＡとＢの間にある点源ＡＳの全てからの、また点ＣとＤ（媒質２上の）の間にある点源Ａ１の一部からの、および点ＦとＧの間にある点源ＡＴの全てからの放射を受ける。しかし、実際には、点Ｐと、点ＦとＧの間にある点源ＡＴとの一部との間に散乱体が存在するため、点Ｐは、点ＥとＧの間である点源ＡＴによってのみ励起される（ＰとＧを結ぶ光線は媒質２に接している）。従って、そのような場合に本発明の方法を実装するために、点Ｐにおいてどの点源が場を励起できるかを判定しなければならず、続いて、この場に対する他の全ての点源の寄与がゼロに設定される。
【０１４８】
散乱体が包含物であるこの同一ケースに対して、計算プログラムを利用するためにより簡便な代替方法について図１８を参照しつつ述べる。
【０１４９】
図１８において、トランスデューサＳとＴの間に、図１７での２個の媒質の代わりに、３個の媒質を示す。媒質２は、図１７の場合と同じ包含物を含んでいて、平面層により拡張されている。媒質２は、包含物を迂回してトランスデューサＳおよびＴの能動面に平行な平坦面に沿って伸びていて、包含物の最大厚さ（包含物の最大厚さとは、例えばトランスデューサの能動面の垂線に沿って決定される）より距離Ｄだけ互いに離れている２個の界面Ｉ１、Ｉ２により境界が画定されている。界面Ｉ１、Ｉ２のこれらの平坦面の中央の面は、ほぼ媒質２の中心Ｃを通る（図ではＣを通る線ＴＲにより辿られる）。媒質１は、トランスデューサＳと界面Ｉ１の間に置かれ、媒質３はトランスデューサＴと界面Ｉ２の間に位置している。
【０１５０】
媒質１、２および３内の場の値は、以下のように点源の２個の層の寄与を加えることにより得られる。
媒質１：ＡＳおよびＡ１により生成された場の合計
媒質２：Ａ＊１およびＡ２により生成された場の合計
媒質３：Ａ＊２およびＡＴにより生成された場の合計
【０１５１】
２個のトランスデューサＳ、Ｔの能動面上の境界条件、および各々が一方では媒質１と２を、他方では媒質２と３を分離している２個の界面を越える際の連続性条件を遵守することにより、以下が得られる。
【数１２１】
2009521031.tif 000123
【０１５２】
図１８に示す例において、遮蔽された領域では小さい領域の僅かな点だけを配置することができ、これを考慮しなくても無視できる誤差しか生じないことに注意されたい。
【０１５３】
媒質３の特性が、媒質１の特性と同じで、および距離Ｄがゼロである場合、図１８の例は図１９に示す例のようになる。図１９において、図１８のものと同じ要素に同じ参照符号が割り当てられていて、仮定により媒質１と３が同じ特性を有し、距離Ｄはゼロである事実により、媒質２は媒質１に完全に囲まれていて、２個の界面Ｉ１およびＩ２が、媒質２の中心Ｃを通っていてＴとＳの能動面に平行な１個の仮想界面Ｉ０に簡約される。図１６に示す場合との違いは点源の決定にある。従って、図１６において、ソースの４個の層はＡＳ、ＡＴ、Ａ１およびＡ＊１であるのに対し、図１９では層ＡＳとＡＴは図１６のものと同一であり、一方、図１６の層Ａ１は図１９では２個の層Ａ１およびＡ＊２に分かれていて、図１６の層Ａ＊１は図１９では２個の層Ａ２およびＡ＊１に分かれている。
【０１５４】
方程式（５１）は、図１９の場合に適合すべく以下のように変形される。
【数１２２】
2009521031.tif 000124
これにより以下が得られる。
【数１２３】
2009521031.tif 000125
方程式（５３）は、以下の行列形式で記述することができる。
【数１２４】
2009521031.tif 000126
【０１５５】
方程式（５０）の解とは対照的に、方程式（５４）の解により、遮蔽された領域に付随する問題を解決することができる。
【０１５６】
系が複数の領域を有する場合に、電磁気場、より具体的にはフーコー電流場における本発明の方法の一実施形態について以下に述べる。
【０１５７】
グリーンの理論は、非均一な波［ネイ］に関する以下のような方程式を解く趣旨で電磁気場において利用できる多くの計算技術の１つを構成する。
【数１２５】
2009521031.tif 000127
この方程式において、
【数１２６】
2009521031.tif 000128
は磁気ベクトルポテンシャル、
【数１２７】
2009521031.tif 000129
は電流密度を表わすと共に、
【数１２８】
2009521031.tif 000130
である。上式において、
ε＝ε０・εｒ（但し、ε０＝１／３６π・１０−９である）であり、
μ＝μ０・εｒ（但し、μ０＝４π・１０−７である）であり、
ここに、導電媒質の場合、ωσμ＞＞ω２μεであり、これによりσ＞＞ωεであり、これにより
【数１２９】
2009521031.tif 000131
が得られ、これにより表皮効果の深さを
【数１３０】
2009521031.tif 000132
と推論することができる。次いで散乱方程式が得られる。
逆の場合、媒質が導電性でなく、
【数１３１】
2009521031.tif 000133
である場合（ｃは波の速度）、これらは伝播方程式であって、波長はλ＝２π／ｋ＝２πｃ／ωと定義される。
【０１５８】
グリーンの理論により、問題を大域的に定式化することができ、これは本発明のＤＰＳＭ法に完全に合致し、上述のように、本発明のＤＰＳＭ法に従い対象とする媒質内の放射している全ての点源からの寄与が合計されている。
【０１５９】
【数１３２】
2009521031.tif 000134
を、ヘルムホルツ方程式のディラック・パルスに対する媒質のパルス応答とすれば、その境界条件は以下のように記述される。
【数１３３】
2009521031.tif 000135
【０１６０】
従って特別解は以下のようになる。
【数１３４】
2009521031.tif 000136
但し、ψは未知であり、ｆが励起ならばψは次式で示される。
【数１３５】
2009521031.tif 000137
【数１３６】
2009521031.tif 000138
はグリーンの関数または方程式（Ｆ−４）に付随する基本解であり、ｆはボリュームＶ内のソースである。
【０１６１】
今度は均一な自由空間の場合について調べる。この場合、ボリュームｖ内に電流密度
【数１３７】
2009521031.tif 000139
が存在する。次いで磁気ベクトルポテンシャルが方程式（Ｆ−１）の解である。方程式（５７）から出発して、グリーン関数を球面座標において得ることができる［ネイ］。
【数１３８】
2009521031.tif 000140
ここに、
【数１３９】
2009521031.tif 000141
は観測点と励起点との間の距離である。
【０１６２】
グリーン関数、方程式（６０）の解は、以下の通りである。
【数１４０】
2009521031.tif 000142
この解は球面波と見なすことができ、１／Ｒずつ減少して点Ｒ→０で特異点を有する。
【０１６３】
励起の
【数１４１】
2009521031.tif 000143
特別解は、重畳積分（５９）を適用することにより得ることができる。この解は方程式（５５）の右辺およびグリーン関数（６１）を用いて以下のように記述できる。
【数１４２】
2009521031.tif 000144
【０１６４】
磁場
【数１４３】
2009521031.tif 000145
および電場
【数１４４】
2009521031.tif 000146
は方程式（６２）から推論することができる。また、この方程式は、観測点
【数１４５】
2009521031.tif 000147
におけるボリュームｖ’内のソース全体の寄与を表わす。グリーンの定式化とＤＰＳＭによる定式化の間には、ＤＰＳＭソースが点源であるという事実を除いて、明らかな類似点がある。方程式（６２）において、
【数１４６】
2009521031.tif 000148
はＡ／ｍ２単位で表わされる電流密度であることに注意されたい。
【０１６５】
ＤＰＳＭ法の定式化によれば、項
【数１４７】
2009521031.tif 000149
はＮｓ個の点源の有限な合計を表わす。
【数１４８】
2009521031.tif 000150
方程式（Ｆ−９）を方程式（Ｆ−８）に代入すれば以下が得られる。
【数１４９】
2009521031.tif 000151
【０１６６】
方程式（６４）において、方程式（６２）から考慮して、
【数１５０】
2009521031.tif 000152
は基本ボリュームｄｖに組み込まれていて基本ソース・ベクトル（三重項ソース、すなわち現在の単位での３個の基本直交座標成分にｍ．ｋ．ｓ単位系のＡ．ｍ単位の長さを乗算したもの）を表わす電流密度
【数１５１】
2009521031.tif 000153
に等しいことに注意されたい。
【０１６７】
直交座標系において、座標
【数１５２】
2009521031.tif 000154
を有する点で計算されたベクトルポテンシャルは、座標（Ｃｘｉ，Ｃｙｉ，Ｃｚｉ）を有する点源および座標（ｘ，ｙ，ｚ）を有する観測点Ｐを用いれば、次式で表わされる。
【数１５３】
2009521031.tif 000155
但し、
【数１５４】
2009521031.tif 000156
および
【数１５５】
2009521031.tif 000157
【０１６８】
本発明のＤＰＳＭ法は、観測点Ｐの所与の配列における全ての点源の寄与の合算を実行する。例えばＮｐ個の観測点が存在する場合、ｒｊの各々の値に対して、磁気ベクトルポテンシャル
【数１５６】
2009521031.tif 000158
は３個のベクトル
【数１５７】
2009521031.tif 000159
を含んでいる。例えば、
【数１５８】
2009521031.tif 000160
は次式で示される。
【数１５９】
2009521031.tif 000161
この式において、励起点ｉと観測点ｊとの間の距離はＲｊｉである。
【０１６９】
最後に、Ａの全ての座標ｘ、ｙおよびｚを含むベクトルＡに対して以下が得られる。
【数１６０】
2009521031.tif 000162
【０１７０】
基本三重項の成分Ｊｘｉ，Ｊｙｉ，Ｊｚｉは正確に同じ点
【数１６１】
2009521031.tif 000163
に置かれておらず、および行列Ｗｘｘ，Ｗｙｙ，Ｗｚｚが僅かに異なり得ることも理解されたい。
【０１７１】
従って、所与の点
【数１６２】
2009521031.tif 000164
における磁場
【数１６３】
2009521031.tif 000165
は以下の方程式により得られる。
【数１６４】
2009521031.tif 000166
【０１７２】
次いで、
【数１６５】
2009521031.tif 000167
但し、
【数１６６】
2009521031.tif 000168
である。
【数１６７】
2009521031.tif 000169
を以下のような曲線座標であるとする。
【数１６８】
2009521031.tif 000170
但し、
【数１６９】
2009521031.tif 000171
単一のソースｉの場合、偏導関数から以下が得られる。
【数１７０】
2009521031.tif 000172
または、
【数１７１】
2009521031.tif 000173
但し、
【数１７２】
2009521031.tif 000174
従って、座標ｒｏｔＡｘについて以下が得られる。
【数１７３】
2009521031.tif 000175
同じことが座標ｒｏｔＡｙおよびｒｏｔＡｚについて成り立つ。
【数１７４】
2009521031.tif 000176
最後に以下が得られる。
【数１７５】
2009521031.tif 000177
但しこの場合も
【数１７６】
2009521031.tif 000178
この方程式はまた、Ｎｐ個の観測点の配列に対して、方程式（Ｆ−１２ｂ）の場合と同様に、行列形式で記述することができる。
【数１７７】
2009521031.tif 000179
この方程式において、Ｂαβは、添え字αおよびβが各々微分軸および基本点源の成分方向の軸を示す行列ブロックである。
【０１７３】
静的な場合（ω＝０）において、ビオ・サバール方程式と正確に同じ結果が得られる。
【数１７８】
2009521031.tif 000180
【０１７４】
最後に、静電気成分が無い場合、静電場
【数１７９】
2009521031.tif 000181
は以下の関係により得られる。
【数１８０】
2009521031.tif 000182
【０１７５】
異なるゲージ、例えばローレンツ・ゲージ（上における、図９の説明直前のゲージに関する段落を参照）を用いて以下が得られる。
【数１８１】
2009521031.tif 000183
【０１７６】
方程式（６８）および（６９）は、本発明のＤＰＳＭ法による電磁場を表わす。点源の概念を実装するためにグリーンの理論をＤＰＳＭ法に適用した。以下、ＤＰＳＭ／グリーン定式化において表面を越える際にポテンシャルの連続性を要求する趣旨で境界条件が記述されている。ＤＰＳＭ定式化の重要な特徴の１つが、当該界面により分離された媒質内の場を総合するために、界面のいずれかの側に仮想点源の配列（電磁気の場合は電流源の三重項）を配置することに注意されたい。大域行列を用いる定式化のために、本発明の方法は、単純な幾何学的形状に限定されず、あらゆる種類の表面、および複数の界面を含む物体に適用することができる。
【０１７７】
境界条件
本発明のＤＰＳＭ法によれば、超音波場、静磁気場、または静電場のいずれにあろうが、ポテンシャルおよび分離面に対する法線
【数１８２】
2009521031.tif 000184
に沿った一次導関数の連続性を要求することにより境界条件が満たされる。電磁気への応用において、この定式化により、磁気ベクトルポテンシャル
【数１８３】
2009521031.tif 000185
および界面の各点における法線
【数１８４】
2009521031.tif 000186
に沿った一次導関数
【数１８５】
2009521031.tif 000187
の連続性が要求される。
【０１７８】
例えば、図２０、２１に示す場合において、磁気ポテンシャル
【数１８６】
2009521031.tif 000188
は、ソース
【数１８７】
2009521031.tif 000189
の効果を重ね合わせることにより媒質１で計算されるのに対し、媒質２における磁気ポテンシャル
【数１８８】
2009521031.tif 000190
はソース
【数１８９】
2009521031.tif 000191
だけから生じる。従って、以下を記述することができる。
【数１９０】
2009521031.tif 000192
【０１７９】
界面に対する法線
【数１９１】
2009521031.tif 000193
に沿ったベクトルポテンシャル
【数１９２】
2009521031.tif 000194
の変型を、一般的な場合に以下のように表わすことができる。
【数１９３】
2009521031.tif 000195
【０１８０】
図２０に示す特定の場合において、ベクトル
【数１９４】
2009521031.tif 000196
はｚ軸に対応し、これは先の方程式が以下のようになることを意味する。
【数１９５】
2009521031.tif 000197
【数１９６】
2009521031.tif 000198
とする。
【数１９７】
2009521031.tif 000199
であるこの特定の場合について、Ｎｐ個の観測点の所与の配列に対して、方程式（６６ｂ）および（６８ｂ）と同一の導関数および表記を用いて以下のように記述することができる。
【数１９８】
2009521031.tif 000200
【０１８１】
図２１において、ソース
【数１９９】
2009521031.tif 000201
（媒質１内に位置している）と、ＤＰＳＭソース
【数２００】
2009521031.tif 000202
（各々、媒質１内へ放射すべく外向き、および媒質２へ放射すべく内向きである）の関連付けられた層を有する媒質間の界面が存在する構成を示す。従って、本発明の方法による境界条件は以下のようになる。
−ポテンシャル
【数２０１】
2009521031.tif 000203
および界面の表面「１」で誘導され、ソース
【数２０２】
2009521031.tif 000204
（インダクタ）により生成されたその一次導関数
【数２０３】
2009521031.tif 000205
【数２０４】
2009521031.tif 000206
−ポテンシャル
【数２０５】
2009521031.tif 000207
および界面の表面「１」で誘導され、ソース
【数２０６】
2009521031.tif 000208
（外部ＤＰＳＭ層）により生成された一次導関数
【数２０７】
2009521031.tif 000209
【数２０８】
2009521031.tif 000210
−ポテンシャル
【数２０９】
2009521031.tif 000211
および界面の表面「１」に誘導され、ソース
【数２１０】
2009521031.tif 000212
（内部ＤＰＳＭ層）により生成された一次導関数
【数２１１】
2009521031.tif 000213
界面に沿う
【数２１２】
2009521031.tif 000214
【数２１３】
2009521031.tif 000215
の値は以下の方程式を満足しなければならない。
【数２１４】
2009521031.tif 000216
【０１８２】
ＪＳが公知であると仮定すれば、以下はＪＳ＝Ｉ・ＪＳと記述することができる。
Ｉを３Ｎｓ×３Ｎｓの単位行列であるとすれば、以下が得られる。
【数２１５】
2009521031.tif 000217
【０１８３】
上により詳細に述べたように、未知数
【数２１６】
2009521031.tif 000218
は方程式（７４ｂ）の反転により得られ、ここに、Ｗ１Ｓが３Ｎｐ×３Ｎｓ次元の行列、Ｗ１１が３Ｎｐ×３Ｎｐ次元の行列、Ｗ２２が３Ｎｐ×３Ｎｐ次元の行列である事実により、大域行列は常に正方行列であり、これは大域行列が正方行列であって３（Ｎｓ＋Ｎｐ＋Ｎｐ）×３（Ｎｓ＋Ｎｐ＋Ｎｐ）次元であることを意味する。
【０１８４】
点源の全ての値が既知であるとの前提で、本発明のモデリングの結果を観察するために空間の所与の離散化された領域の場が計算される。
【０１８５】
図２０の実験装置で得られた結果をここに記述する。しかし、この明らかに簡単な問題を解析的に解くのは誘導ループが界面の平坦面に平行な場合でさえ極めて困難である［ドッド６８］ことに注意されたい。更に、他のシナリオにおいて、例えば誘導ループが界面に対して傾斜している場合、解析的計算が更に複雑になるのに対し、ＤＰＳＭ解は簡単なままである。
【０１８６】
図２０に、半無限大の厚さを有すると考えられる導電材料ＭＥ１のブロックの上面から距離Ｄに、電圧Ｖ（ｔ）の電力が供給される導電ループＢ１が配置されているようなブロックを示す。当該ブロックＭＥ１は、直方体の形状を有し、その上面がループＢ１に向き合い、その寸法は電圧Ｖ（ｔ）の波長より大きいことから、半無限大であると考えられる。ループＢ１は、ブロックＭＥ１内で磁場ＣＭを誘導する。
【０１８７】
図２２、２３が関係する一例示的な実施形態によれば、ループＢ１は直径５０ｍｍであって、ＮＳ個の電流要素ＪＳ（仮想励起ソース）によりメッシュモデル化されていて、後者からＤ＝２０ｍｍの距離で第１にブロックＭＥ１の上面と平行であり、ブロックＭＥ１の上面は３００×３００ｍｍの正方形であって、空気／ブロック界面はＮＰ＝２２×２２＝４８４個の基本ソース
【数２１７】
2009521031.tif 000219
の三重項により当該界面のいずれかの側でメッシュモデル化されている。
【０１８８】
図２１において、本発明の方法による仮想点源の分布を模式的に示す。ループをモデル化すべく、先に定まったように、ループＢ１の平面内の励起一次点源ＪＳの層、および空気／ブロック界面ＭＥ１のいずれかの側の２次点源
【数２１８】
2009521031.tif 000220
の層を示し、上で述べたように、ソース
【数２１９】
2009521031.tif 000221
は外部空気の側にある。ループＢ１が単一ターンであるために、ソースＪＳのピッチが他のソースのピッチより大きいことに注意されたい。
【０１８９】
図２２、２３に、空気／ブロック界面ＭＥ１に平行な平面内における各々ソース
【数２２０】
2009521031.tif 000222
の実部の分布の上面図を示し、ブロックＢ１はアルミニウム製であって電圧Ｖ（ｔ）の周波数が１０００Ｈｚである。同様に、図２４および２５に、ループＢ１が１０°傾斜していて、ループの中心とブロックの間の距離が５０ｍｍ、電圧Ｖ（ｔ）の周波数が１００Ｈｚである場合におけるこれらのソースの分布を示す。
【０１９０】
最後に、図２６、２７において、図２４および２５に対応する場合について、磁場Ｈの実部および虚部を３種の異なるレベル、すなわち、ブロックＭＥ１の上面の平面、ブロック内で当該上面から５ｍｍに位置している平面、同様にブロック内部で当該上面から１０ｍｍに位置している表面において示す。図２２〜２７におけるこれら全ての座標はミリメーター尺度である。
【０１９１】
当然ながら、本発明は異方性媒質に実装可能であり、その場合の基本特性関数も解析的に得られる。
【０１９２】
ＤＰＳＭソースは大多数の場合、モデル化したい系を構成する物体の真の影響の代わりに、所与の空間領域内の物理量を総合することを唯一の目的として配置（特に界面のいずれかの側に）されているため、「架空の」ソースである。
【０１９３】
しかし、これらのＤＰＳＭソースがこの程度には架空でない、いくつかの場合が存在する。例えば、静電気において、表面の近傍におけるＤＰＳＭ電荷分布は時折、物体の表面に蓄積された静電荷の真の分布を相当忠実に表わすことができ、コンデンサ電極の場合がこの特性をよく例示している。これらの電荷は次いで、真の静電荷と同じ法則に従う。
【０１９４】
この興味深い特性は、電気的に絶縁されているため全電荷がゼロのままであることを観察することにより、静電場において「浮動的ポテンシャルで」あるといわれる物体が獲得したポテンシャルを計算する例を利用して示すことができる。
【０１９５】
以下に、２個の導電性物体ＡおよびＢが均一な媒質内に配置されている簡素化された例を示す。ポテンシャルｖ１が物体Ａ（０基準電位が無限遠に置かれている）に印加されている。Ａに出現している電荷が、導体内部においてゼロ電場を保証すべく導体Ｂの自由電荷をその表面へ移動させる静電場を生成する。導体Ｂは、従って、均一なポテンシャルｖ２を有し、全電荷をゼロに維持する。
【０１９６】
理論的には、ポテンシャルｖ２は系の静電容量係数により決定される値を有する（「０．５〜２０ＧＨｚの間の周波数を処理するＲＦ  ＭＥＭＳ可変静電容量の研究、設計、製造および試験」Ａｕｒｅｌｉｅ  Ｃｒｕａｕ論文、２００５年１月、第４章、ｐ１２０およびｐ１２５参照）。
【数２２１】
2009521031.tif 000223
【０１９７】
ＤＰＳＭ法の場合、各導体は点源Ａｉのアセンブリ（外側へ放射している内側ソースを明灰色の正方形で示す）および、境界条件が適用される（黒丸）表面に配置された試験点の対応するアセンブリで表わされる。
【数２２２】
2009521031.tif 000224
【０１９８】
この例では、物体ＡおよびＢは、各々Ｎ１およびＮ２個の試験点によりメッシュモデル化されていて、各々がこの場合は静電荷であるＮ１およびＮ２個の点源（層Ａ１、Ａ２）を含んでいる。Ｐ１を物体Ａ上の全ての試験点におけるポテンシャルの値を含むベクトル（Ｐ１はユーザが指定）とし、Ｐ２を、物体Ｂ上の全ての試験点におけるポテンシャルの値を含むベクトルとする。ここで、物体Ｂが導電性であって、その結果Ｂの表面におけるポテンシャルが均一であると仮定するとＰ２は以下のように記述することができる。
【数２２３】
2009521031.tif 000225
（ユーザが物体Ａ上の全ての試験点に同じポテンシャルを印加することに決めたならばＰ１についても同じことが当てはまる）。
【０１９９】
この問題を解くために、ユーザが指定した物体Ａ上での境界条件を考慮にいれて、ソースＡ１およびＡ２の値を物体Ｂ上のポテンシャルＰ２と共に計算しければならない。従って、Ｎ１＋Ｎ２＋１個の未知数があるため、Ｎ１＋Ｎ２＋１個の方程式が必要である。
【０２００】
ＤＰＳＭ法は、上述の特性を考慮に入れながら、この問題を以下の方程式系により表わす。
【数２２４】
2009521031.tif 000226
ここに、Ａｉは導体ｉの内部電荷ベクトル、ＱｉｊはＤＰＳＭ結合行列、およびＰｉは導体ｉのポテンシャルベクトルである。
【０２０１】
未知の値はＡ１、Ａ２およびｖ２である。
【０２０２】
以下の行列が得られ、その逆行列を計算するだけで上記結果が得られる。
【数２２５】
2009521031.tif 000227
【０２０３】
引用文献
［ＤＯＤＤ  ６８］：Ｃ．Ｖ．Ｄｏｄｄ  ａｎｄ  Ｗ．Ｅ．Ｄｅｅｄ，‘Ａｎａｌｙｔｉｃａｌ  ｓｏｌｕｔｉｏｎｓ  ｔｏ  ｅｄｄｙ  ｃｕｒｒｅｎｔ  ｐｒｏｂｅ−ｃｏｉｌ  ｐｒｏｂｌｅｍｓ’，Ｊｏｕｒｎａｌ  ｏｆ  Ａｐｐｌｉｅｄ  Ｐｈｙｓｉｃｓ，３９（６），ｐａｇｅｓ  ２８２９  ａ  ２８３８，１９６８．
［Ｎｅｙ］：  Ｍｉｃｈｅｌ  Ｎｅｙ，‘Ｂａｓｅｓ  ｄｅ  ｌ’ｅｌｅｃｔｒｏｍａｇｎｅｔｉｓｍｅ’，Ｔｅｃｈｎｉｑｕｅｓ  ｄｅ  ｌ’Ｉｎｇｅｎｉｅｕｒ，ｒｅｆ．Ｅ１０２０．
【図面の簡単な説明】
【０２０４】
【図１】本発明の方法の一実施形態の簡素化されたフロー図である。
【図２−８】本発明の方法が実装されている系の簡素化された模式図である。
【図９−１９】本発明の方法が実装されている非均一媒質の極めて簡素化された模式図である。
【図２０−２７】フーコー電流を用いる非破壊試験への本発明の方法の応用に関する模式図である。
【図１】
2009521031.tif 000228
【図２】
2009521031.tif 000229
【図３】
2009521031.tif 000230
【図４】
2009521031.tif 000231
【図５】
2009521031.tif 000232
【図６】
2009521031.tif 000233
【図７】
2009521031.tif 000234
【図８】
2009521031.tif 000235
【図９】
2009521031.tif 000236
【図１０】
2009521031.tif 000237
【図１１】
2009521031.tif 000238
【図１２】
2009521031.tif 000239
【図１３】
2009521031.tif 000240
【図１４】
2009521031.tif 000241
【図１５】
2009521031.tif 000242
【図１６】
2009521031.tif 000243
【図１７】
2009521031.tif 000244
【図１８】
2009521031.tif 000245
【図１９】
2009521031.tif 000246
【図２０】
2009521031.tif 000247
【図２１】
2009521031.tif 000248
【図２２】
2009521031.tif 000249
【図２３】
2009521031.tif 000250
【図２４】
2009521031.tif 000251
【図２５】
2009521031.tif 000252
【図２６】
2009521031.tif 000253
【図２７】
2009521031.tif 000254
【国際調査報告】
2009521031.tif 000255
2009521031.tif 000256
2009521031.tif 000257
2009521031.tif 000258
2009521031.tif 000259
2009521031.tif 000260