固有値法を用いた異方性繊維強化シリコンにおける結合熱・光弾性場の時空間解析

サマリー

本研究は、異方性繊維強化シリコンにおける結合された熱・光弾性場を、ノーマルモード法と固有価値法を用いて解析します。結果は空間減衰と時間依存の場の進化を示し、異方性に対する強い感度を示しています。ヒートマップはフィールドの分布と局在を示しています。

要約

本研究は異方性繊維強化シリコン半導体媒体における結合熱光弾性系を調査し、熱場、キャリア場、機械場間の相互作用を捉えることを目的としています。このような異方性繊維強化材料は、マイクロエレクトロニクスや光電子デバイス、レーザー技術、センサー、高度な複合構造など、方向性や機械性能の向上が求められる現代の工学的応用において重要な役割を果たしています。特に熱的および光学的負荷がかかる半導体部品の設計において重要であり、結合電界挙動の正確な予測が信頼性と性能最適化に不可欠です。支配方程式は結合物理モデルに基づいて定式化され、解析を簡素化し、関与するパラメータの相対的な影響を強調するために無次元形式に変換されます。この問題は正規モード手法を用いて解かれ、一次ベクトル行列微分系に簡約され、その後固有値アプローチで半無限領域内で課された境界条件を満たす解析解を得る。数値解析を用いて時間変動がすべての物理場に与える影響を調べ、強い空間減衰と異方性およびファイバー強化による結合挙動を明らかにします。時空的ヒートマップ表現は、フィールドの進化と局所化を可視化するために用いられ、多物理的相互作用の物理的洞察を提供し、分析手法の有効性を示します。

概要

熱光弾性は、半導体材料における光熱励起下での熱、機械、光学効果の相互作用を説明する重要な学際的分野として浮上しています。これらの場の結合は、レーザー加熱や光励起に関わる現代の応用で特に顕著になります。例えば、Saeed¹ は双曲型二温度モデルを用いて半導体における熱・光弾性相互作用を調査し、熱緩和効果がシステム挙動の正確な予測に重要であることを示しました。繊維強化複合材料の熱機械的応答性も、その向上した機械的強度と異方性の特性により大きな注目を集めています。LiとLambros² は、このような複合材料の動的熱機械的挙動を解析し、高度な工学的応用への適性を強調しました。同様に、Kalkalらは磁場下で回転する機能段階の繊維強化媒体における二次元(2D)変^形を研究 し、異方性や外部の影響がシステム応答に強い影響を与えることを示しました。さらに、ピタレシらは熱弾性挙動における巨視的異質性^{の役割を検証} し、複合材料の正確なモデリングの必要性を強調しました。近年の発展により、これらの研究は結合多物理効果を持つ半導体媒体にも拡大されました。Mondalら5は、記憶応答と磁場を考慮した強化半導体における波動伝播^を調査し 、複雑な結合相互作用を明らかにしました。Akaiらの実験研究6^では、熱弾性温度変化を用いた繊維強化複合材料の疲労損傷を評価し、その実用的重要性が確認されています。さらに、Luら7によって示されたように、効果的な熱弾性特性を推定するためにミクロ力学的手法が用い^{られています。}外部荷重や環境条件の影響も広く研究されています。BarakとDhankhar⁸ は機能的に段階的に調整された繊維強化媒体における傾斜荷重を解析し、KunduとKalkal⁹ は重力下および移動熱荷重下での光熱相互作用を調査しました。Chaudharyら¹⁰ は二重相遅延モデルを用いて温度依存性を調査し、Panditら¹¹ は非局所変形挙動を捉えるために分数順序ひずみモデルを適用しました。これらの研究は、熱弾性解析において現実的な荷重や材料条件を考慮する重要性を強調しています。

動的および振動に関連する現象も半導体システムで探求されています。Songら¹²は半導体構造中の光熱振動を研究し、Mondalと^Sur13は記憶効果を伴う直交対立媒体における波動伝播を解析しました。粘弾性および微細構造効果はAbouelregalら¹⁴によって検討され、スランプ型加熱効果はHobinyら¹⁵によって調査されています。さらに、予測精度向上のためにZenkour¹⁶によって、三相遅延理論のような一般化された熱弾性モデルが開発されました。電磁効果や微小極性効果を取り入れた高度な多物理モデルは、熱光弾性挙動の理解をさらに深めています。Al-Hazaemhら¹⁷は回転半導体媒体中の光電磁熱弾性励起を研究し、Nazirと^Kumar18は微小極熱弾性相互作用を解析しました。Songら¹⁹は非散逸性熱弾性相互作用も調査し、NasrとAbouelregal20はキャビティを持つ半導体における光吸収過程を調査しました。非局所的および分数順序モデルは、近年の発展において重要な役割を果たしています。Guptaら²¹は非局所的多孔質媒体中の光熱励起を研究し、HobinyとAbbas²²は半導体における分数順序波の伝播を解析しました。HafedとZenkour²³は傾斜荷重効果を調査し、Oliinykら²⁴は非定常熱光弾性効果を調査しました。また、Awwadら²⁵はレーザー励起下での機能的に段階化された半導体挙動を探求し、Guptaら²⁶はメモリ依存モデルを用いた熱圧電光電結合を調査しました。GreeneとPatterson(²⁷)、BaroneとPatterson(²⁸)によって熱的および光学的手法を組み合わせた古典的実験的アプローチが確立され、応力解析の信頼できる手法を提供しています。さらに、KaurとSingh⁽²⁹)は半導体共振器の非局所的メモリ依存モデルを開発し、Abbasらは可変熱伝導率を持つ光熱^{相互作用を解析}しました。強化複合膜の実験的および数値的調査はLuらによって行われました。³¹、PurkaitとKanori³²は回転する繊維強化媒体における記憶応答を研究しました。Abo-Dahabら^{はさらに、}繊維強化熱弾性媒体における荷重条件下での波の反射を調査しました。近年では、複雑な材料を記述するために高度な分数モデルや非局所的熱弾性モデルが提案されています。Abouelregalら³⁴は分数モデルを用いて生体組織の熱応答を調査し、Selvamaniらはナノビームにおける非局所波伝播と振動挙^動を研究しました。さらに、Abouelregalらによって二相遅延および粘弾性モデルが微細構造に適用^{され、複雑な}結合現象を捉えるためにメモリカーネルを用いた分数熱弾性定式化^{も開発されています}。

これらの広範な発展にもかかわらず、既存の多くの研究は主に解析的または数値的解に焦点を当てており、物理場の詳細な時空間的可視化を提供していません。多くの実用的応用、特に異方性繊維強化半導体媒体では、システムの応答は空間的および時間的な変動に強く依存するため、正確な解釈のために可視化が不可欠です。さらに、数値的および実験的アプローチは貴重な洞察を提供しますが、正規モードや固有値技術に基づく解析手法は、半無限領域や結合多物理系に関わる問題において大きな利点を提供します。これらの手法は閉形式解を可能にし、波の伝搬、減衰、結合メカニズムの物理的理解を深め、数値的および実験的結果の検証に信頼できるベンチマークとして機能します。

本研究では、異方性繊維強化半導体媒体における結合熱光弾性挙動の包括的な調査を提示します。本研究の新規性は、固有値に基づく分析手法と時空的ヒートマップ可視化を統合し、物理場の進化と局在化についてより深い洞察を提供している点にあります。本研究の目的は、異方性とファイバー強化が熱場、機械場、キャリア場間の相互作用に与える影響を分析し、半導体デバイス、センサー、レーザー技術など現代の工学応用に関連する複雑な多物理現象の解釈における提案手法の適用可能性を示すことです。

プロトコル

この研究は完全に理論的モデリングと数値シミュレーションに基づいており、人間参加者、動物被験者、生物学的標本は関与していません。したがって、倫理的な承認やインフォームド・コンセントは必要ありませんでした。

繊維強化異方性媒体における光熱弾性の数学的定式化
本研究では、2次元ファイバー強化異方性半導体半空間を表面光励起にかけることを検討しました。媒質は x ≥0 領域を占めており、x = 0 の境界は露出した表面を表します。座標系は、x軸が媒質に伸び、y軸が表面に沿って配置され、面内挙動を記述するように定義されました。材料は均質だが異方的であると仮定されていたのは、整列した補強繊維の存在により弾性および結合特性に方向性依存性が生じたためである。表面での光学吸収は局所的な加熱および過剰荷荷キャリアを生み出し、熱場、機械場、キャリア場間の完全結合相互作用をもたらしました。したがって、系の状態は温度θ(x, y, t)(K)、キャリア密度N(x, y, t)(^m-3)、変位成分u(x, y, t)およびv(x, y, t)(m)によって記述され、小さな変形を仮定しました。物理領域、座標系、ファイバーの向き、適用された光励起の回路図は図1に示されています。すべての記号計算および数値計算はWolfram Mathematica(バージョン12.0)を使用して行われました。

図1。境界 x = 0 で光励起を受けた半無限ファイバー強化半導体媒体の回路図表現。 座標系(x, y)を示し、ファイバーの向きは x方向(a = (1, 0))に沿って整列しており、媒質の幾何学的配置と方向依存の異方性を示しています。 この図の拡大版はこちらをクリックしてご覧ください。

繊維強化異方性熱弾性半導体媒体における応力テンソルの本源関係は、式1 ^1,5を用いて一般形で表されました。この定式化では、θは基準温度T₀に対する温度増分を表し、Tは該当する場合の絶対温度を表します。

.  (1)

ここで、C_ijkl は弾性剛性係数、e_kl はひずみテンソル、β_ij と η_ij はそれぞれ熱弾性テンソルとキャリア結合テンソルを表します。繊維補強が存在する場合、材料応答は方向依存となり、ファイバーの向きベクトル a = (a_i) によって支配され、弾性項と結合項の両方に異方的寄与が導入されました。したがって、構成関係はファイバー強化の効果を明示的に^2,3として組み込むように拡張されました。

.(2)

ここで、λとμ_τはラメ定数であり、_{μ L}は繊維方向に沿った縦方向のせん断率です。パラメータαはファイバー補強効果を表し、熱膨張係数を示す α_ij とは異なります。単位ベクトルはファイバーの配向を定義し、応力-ひずみ応答に方向依存性を導入しました。現在の2次元定式では、ファイバーはx軸に沿って整列していると仮定されていました。したがって、向きベクトルは明示的にa = (1, 0)として取られました。この仕様はファイバー方向の明確なパラメータ化を提供し、異方性寄与が支配方程式に一貫して組み込まれることを保証し、ファイバー強化によって誘導される方向挙動に直接対処しました。現在の2次元構成では、支配応力成分は以下に還元されます:

, (3)

、( 4)

.(5)

これらの方程式は異方性、ファイバー強化、多物理的結合効果の複合的な影響を示しています。係数β_ijおよびη_ijは、材料パラメータに基づいて次のように定義されました。

,

,

,

.

ここで、係数 _{A ij} は繊維強化異方性媒体の有効弾性定数を表し、次のように定義されました。

.(6)

ここで、λ、μ_L、μ_Tは異方性繊維強化媒体の弾性定数であり、α_ijとξ_ijはそれぞれ熱膨張係数とキャリア膨張係数を表します。熱光弾性半導体媒体における弾性波の伝播は、動的熱弾性解析の基礎となる線運動量保存の原理によって制御されていました。物体力が存在しない場合、変形可能な連続体の一般運動方程式は^1,15に基づいて次のように表されます。

.(7)

ここで、ρは質量密度、σ_ijは応力テンソルです。本研究では、定式化はx-y平面上の2次元配置に限定され、変位場はu(x, y, t)とv(x, y, t)で表されました。熱光弾性媒体における標準的な定式に従い、二次元における運動の支配方程式は次のように表記されました。

, (8)

.(9)

異方性繊維強化本構関係を上記の方程式に代入すると、次のように結合された偏微分方程式系(DE)が得られます。

, (10)

.(11)

ここで添字は空間変数および時間変数に対する部分微分を示します。これらの方程式は、異方性、ファイバー強化、温度勾配、キャリア拡散が媒質の動的応答に与える結合的影響を強調しています。光励起が存在する場合、半導体内部の熱場はキャリア密度や機械的変形との相互作用によって強く影響を受け、完全結合されたエネルギー輸送過程が形成されました。古典的な熱伝導とは異なり、このような媒体の温度変化はキャリアの再結合や熱弾性効果から生じる追加の源項によって制御され、熱伝播特性を大きく変化させました。一般化熱弾性の枠組みにおける熱伝導方程式は次のように表^{されました。}

.(12)

ここで、C_Eは一定ひずみでの比熱であり、材料の熱容量を表し、T₀は平衡状態の媒質の基準絶対温度を表します。現在の2次元構成では、この式は^16,20に簡約されます。

.(13)

この方程式は、温度場が方向性熱伝導率だけでなく、キャリアの 項による再結合や熱弾性結合項による時間依存変形にも影響を受けていることを示しています。この定式化は異方性繊維強化半導体における熱伝達を支配する本質的な多物理的相互作用を捉え、キャリアの動力学と機械的応答が系の熱挙動を変化させる役割を強調しました。半導体媒体に光励起を加えると、入射放射線の吸収により多くの電荷キャリアが生成されました。これらのキャリアは空間拡散、再結合、熱駆動生成などの輸送過程を経ており、これらはすべて物質内の温度場と本質的に結びついていました。その結果、キャリア密度は結合した熱・光弾性応答を支配する主要な変数の一つとなりました。

本定式では、キャリア濃度の進化は拡散機構、崩壊効果、熱活性化過程のバランスを通じて記述され、次の支配関係^1,5"

.(14)

ここで、D_Eはキャリア拡散係数を表し、はx-y平面における2次元ラプラシアン演算子を表します。項は緩和時間τに伴う組換え効果を説明し、kはとして定義される熱キャリア結合係数であり、平衡キャリア濃度N ₀の温度変化に対する感度を示します。この関係は、キャリア生成の駆動機構としての温度の役割を強調し、異方性ファイバー強化半導体媒体における熱場と電子場の直接結合を確立します。

結合光熱弾性キャリア系の支配方程式と数学的定式化が確立されています。シリコン(Si)媒体に対応する物理的および材料パラメータは、それらの数値、単位、対応する参照資料とともに表 1にまとめられています。これらのパラメータはその後、数値計算や無次元化過程で用いられます。

シンボル	価値	部隊	参考文献
λ	3.64 × 10¹⁰	N/m²	12
μT	5.46 × 10¹⁰	N/m²	12
μL	3.20 × 10¹⁰	N/m²	12
ρ	2330	kg/m³	13
CE	695	J/(kg·K)	30
幼稚園11	0.0921 × 10³	W/(m·K)	30
K22	0.0963 × 10³	W/(m·K)	30
DE	2.5 × 10⁻³	m²/s	22
τ	5 × 10⁻⁵	s	15
T₀	300	K	15
Eg	1.11 × 10⁻¹⁹	J	12
α11	3.1 × 10⁻⁶	K⁻¹	30
α22	3.5 × 10⁻⁶	K⁻¹	30
ξ11	−7 × 10⁻³¹	m³	21
ξ22	−9 × 10⁻³¹	m³	21
κ	2.16 × 10²¹	m⁻³·s⁻¹·K⁻¹	21
α	−1.28 × 10¹⁰	N/m²	28
β	220.90 × 10¹⁰	N/m²	28
ω	2.95 + 1i	S⁻¹	12
a	1	—(無次元)	13
y	0.6	m	13
θ₀	1	—(無次元)	15
N₀	1	—(無次元)	15

表1。異方性繊維強化半導体媒体の数値解析に用いられる材料特性とパラメータ。特に指定 がない限り、すべての量はSI単位で表されます。無次元パラメータはそれに応じて示されています。これらの値は、引用文献から得られたシリコンベースの材料特性およびモデルパラメータに対応しています。熱キャリア結合係数κは、熱・光弾性半導体モデルの標準的な定式化に従い、κ = (∂N₀/∂T)(1/τ)として定義されます。

結合異方光・熱弾性モデルの無次元定式化
支配方程式を簡略化し、結合された熱光弾性系の一貫した無次元表現を得るために、空間座標x、y、時間t、変位成分u、v、温度T、キャリア密度N、応力σに適した特性スケールが導入されました.これらのスケーリングパラメータは、媒体の本質的な物理特性および熱場、機械場、キャリア場間の結合機構に基づき、文献^16,21で確立された定式化に従って一貫して選定されました。したがって、無次元変数は次のように定義されました。

、 、 、 、 、 、 、 、 、 。

この変換により独立した材料パラメータの数が減り、結合系の正規化された表現が得られました。上記の無次元変数を、以前導出された支配方程式に代入することで、システムは無次元形式に書き換えられました。簡便化のため、無次元変数に関連する素数記法は後に省略されました。この手順により、無次元部分微分方程式のコンパクトな集合が得られ、以下の形で表すことができます。

, (15)

, (16)

, (17)

.(18)

無次元変換を適用した後、系の応力成分は以下の正規化形式で表記されました。

, (19)

, (20)

.(21)

無次元パラメータ a_i 結合異方光熱弾性挙動を支配する物理的および材料的性質のコンパクトな組み合わせを表すために導入されました。各係数はシステム内の特定の相互作用機構を反映し、基礎となる物理過程の相対的な影響を洞察しました。 は法線結合剛性と主弾性剛性との比率を表し、2つの変位成分間の異方性相互作用の度合いを反映しています。 は横方向変形が垂直応力成分にどれほど寄与するかを特徴付けます。 熱弾性結合の方向変化を測定し、熱膨張効果における異方性を示します。 キャリア密度が誘導弾性変形に与える異方的な影響を記述します。 は正規化されたせん断剛性を表し、せん断変形の法線変形に対する寄与を定量化します。 制御変位方程式における法線変形とせん断変形の結合を考慮します。 横方向剛性とせん断剛性の比率を表し、異方性変形挙動を強調します。 は正規化された慣性パラメータを表し、波の伝播効果とせん断剛性を関連付けます。 異なる空間方向における変位勾配の結合を特徴付けます。 横方向における変位場に対する熱効果の相対的な寄与を定量化します。 キャリア誘起変形がせん断剛性に対してどのように影響するかを測定します。 は異なる空間方向に沿った熱伝導率の異方性を表します。 キャリアの再結合が媒体内の熱発生に与える影響を特徴付けます。 熱的影響と時間依存弾性変形の結合を表します。 は、両方向における異方熱膨張の複合的な影響を説明しています。 はキャリア輸送速度を制御する正規化拡散パラメータを表します。 キャリア組換え効果の相対的な強さを特徴付ける。 熱変動とキャリア生成過程の結合を記述します。

正規モード技術を用いた解析解
結合異方熱光弾性系の解析解を得るために、支配的な部分微分方程をより扱いやすい通常の微分方程式系に還元する効果があったため、ノーマルモード技術が用いられました。この手法は、分散や減衰を含む波の伝播現象の解析に広く用いられています。したがって、場変数の時間的および横方向の空間方向における調和変動は^1、12、23と仮定されました。したがって、変位成分、温度、キャリア密度、応力は指数関数的に次のように表されます。

.(22)

ここで、ω は場の時間的挙動を支配する複素周波数を表し、 a は y方向に沿った空間変動に関連する波数を表します。これらのパラメータは安定性要件を満たし、半無限領域内で物理的に許容される有界解を保証するために選ばれました。上記の仮定形式を、以前に導出された無次元支配方程式に代入し、得られる式を簡略化することで、元の部分微分方程式の結合系は空間座標に関する通常の微分方程式系に還元され、次のように書けます。

, (23)

, (24)

, (25)

.(26)

さらに、変換領域内の対応する応力成分は次のように表記されました。

, (27)

, (28)

.(29)

ここで 、D は微分作用素 を表します。これらの方程式は正規モード領域における支配系の還元形を表し、特性方程式の導出や後の一般的な解析解の構成の基礎を提供します。係数は次のように定義されました: 、 、 、 、 、 、 、 、 、 。

マトリックス微方方程式の定式化と固有値解析
正規モード変換の適用後、式23–26で示された支配系は空間座標に関する2階の通常の微分方程式の集合に還元されました。体系的な解を容易にするために、このシステムは体量の一次微分に対応する補助変数を導入することで、同等の一階系に変換されました。具体的には、以下の変数が定義されました。

、 。(30)

これらの定義を用いて、式23–26は次の8階微分方程式の系に書き換えられました。

, (31)

, (32)

, (33)

, (34)

.(35)

上記の系は、コンパクト行列 A 形式で次のように表現されました。

.(36)

状態ベクトルは次のように与えられました。

.(37)

そしてシステム行列は明示的な形をとります:

.(38)

この定式化は元のシステムを固有^{値問題1,15}に変換しました。特性方程式は次のように得られた。

.(39)

これにより固有値を支配する8次多項式が得られます。約約形式では、特性多項式は次のように書けます。

.(40)

ここでZ_iは系数の関数であり、以下で明示的に定義されます。得られる固有値は解の空間的挙動、減衰や伝搬特性を決定します。Re(m) > 0を満たす固有値のみが保持され、x → ∞に伴う指数関数的に減衰する物理的に許容可能な解を保証します。

.(41)

特性多項式の根は固有値 mを定義し、解の空間的振る舞いを支配します。これらの固有値は、 特性多項式関数 を用いて特徴多項式を構築し、 NSolveを用いて代数方程式を解くことでMathematicaを用いて数値的に計算されました。この問題は半無限定義域(x ≥ 0)で定式化されているため、 x → ∞に有界な物理的に許容可能な解のみが考慮されます。したがって、Re(m)>0を満たす固有値のみが保持され、x→ ∞としてexp(−mx)の形を指数関数的に減衰する解が保証された。残りの根は、モデルの物理的要件と整合しない非崩壊または有界なしの解に対応するため、破棄されました。

保持される固有値 mに対して、対応する固有ベクトルは対応する代数系から得られた

, (42)

および次の形で表現されました:

.(43)

上記の行列方程式を展開すると、以下の線形方程式系が得られました。

, (44)

, (45)

, (46)

, (47)

.(48)

固有値問題の均質性により、固有ベクトルは任意の乗法定数まで定義されました。一意かつ一貫した表現を得るために、固有ベクトルの1成分を固定することで正規化条件が課されました。本研究では、最初の成分がq1 = 1となり、残りの成分は上記の方程式系から順次決定されました。計算的観点からは、この正規化は1つの成分に単位値を割り当て、その結果得られた線形方程式系を解いて残りの成分を評価することで実装されました。この手法により、各許容される固有値に関連する固有ベクトルを体系的かつ再現可能に計算する方法が提供されました。

.(49)

そして残りの構成要素はシステム関係からそれに従って導かれます。これらの固有ベクトルは、各モード内での温度、キャリア密度、変位場の相対的な寄与を表します。したがって、問題の一般解は許容される固有モードの線形結合として構成され、それぞれが固有値と対応する固有ベクトルに対応し、半空間媒体における結合異方光熱弾性挙動の完全な解析的記述を提供した。したがって、この体系の一般的な解は次のように表されます。

.(50)

ここで、 C_i は境界条件から定められた定数です。上記のベクトル式を展開することで、場の変数は次のように得られます。

, (51)

, (52)

, (53)

.(54)

この表現は、解が指数関数的なモードの重ね合わせで構成され、各固有値-固有ベクトルのペアが全体の物理応答に独立して寄与していることを示しています。許容される固有値は、実部が正になるように選択され、 x → ∞ の時に有界かつ物理的に意味のある解が保証されます。

境界条件と物理的制約
一般解を x = 0 に代入することで、定数 C_i を用いた線形代数方程式系が得られた。具体的には、各境界条件(温度、キャリア密度、変位制約)は固有モード展開で表現され、係数C_iを関連付ける方程式の集合が生成されました。この手順により、BC = Dと行列形式で書ける線形系が生まれました。ここでBは境界で評価された固有ベクトルの成分から構成される係数行列、C = (C₁,_{C 2}, C₃, C₄)^は未知定数のベクトルであり、θ₀などの境界値から決定されます。 N₀、そして変位制約。得られた線形系は、係数行列と右辺ベクトルを明示的に組み立て、未知の定数はLinearSolveルーチンで求めるMathematicaを用いて計算的に解かれました。これらの定数は一般解に代入され、物理場の完全な式を構築し、その後数値評価や結果の図示に用いられました。

課された境界条件は以下の通りです。

温度制約:

.(55)

この条件は、周期的な光学加熱によって誘導される調和的に変化する表面温度を表しています。これは媒質内の結合された熱弾性およびキャリア輸送過程を駆動する一次熱励起として機能します。振幅θ0は加えられた熱荷重の強度を示します。

キャリア密度制約:

.(56)

この境界条件は、光照明によって生じる光生成キャリア密度を記述します。これは光子吸収による電子励起と、入射光場と整合した調和変調を反映します。

変位制約:

.(57)

この条件は、境界が横方向に機械的に制約されていることを示します。したがって、表面の v方向に変位は起こりません。

せん断応力制約:

.(58)

この条件はせん断応力に関してトラクションフリー境界に対応します。これは接線方向における機械的に自由な境界と整合的に表面に接力が作用しないことを保証します。 x = 0 での境界条件に加え、無限遠での物理的要件は次のように課されました。 半無限領域内で有界物理解を保証すること。数値結果を提示する前に、本研究で採用された全体の計算手順を 図2にまとめています。計算に用いられる励起パラメータθ₀、N₀、複素周波数ω、波数aの数値値は 表1に示されています。 表1 に記載されているパラメータには、正規化された定式化で用いられる次元材料定数と無次元パラメータの両方が含まれます。数値評価のために、空間領域は 、横座標は y = 0.6に固定され、時間領域は の範囲内で考慮されました。これらの範囲はすべての数値計算およびグラフィカル表現に使用されました。

図2。提案手法の計算ワークフロー。 図は定式化から数値結果までのステップの順序を示しています:支配方程式、無次元化、正規モード法の適用、一階系への変換、行列の定式化、固有値および固有ベクトル解析、境界条件の適用、定数の決定、数値プロットの生成。 この図の拡大版はこちらをクリックしてご覧ください。

結果

数値結果
このセクションでは、異方性繊維強化半導体媒体における結合熱光弾性キャリア系の挙動を数値計算で解析しました。検討された材料はシリコン(Si)であり、その物理的および材料パラメータは 表1に記載されています。 これらの物質定数は直接支配方程式に代入され、場変数を評価する数値計算に実装されました。すべてのパラメータは定義された無次元スケーリングを用いて一貫して組み込まれました。数値計算は、 Mathematicaのデフォルトの精度と収束設定を使って行われました。空間的および時間的方向の両方でステップサイズ0.01の均一計算グリッドが実装されました。

まず、支配方程式は無次元形式で定式化され、一次場変数の観点から表現されました。正規モードと固有変数のアプローチで得られた解析式は、記号的および数値的評価のためにMathematicaで実装されました。物質定数は 表1 に従って割り当てられ、導...

ディスカッション

得られた結果は、異方性繊維強化半導体媒体における結合された熱光弾性挙動について明確な物理的洞察を提供します。本研究は、熱荷重、キャリア生成、弾性変形の相互作用を調査するための固有値ベースの解析フレームワークを提案します。観測される反応は、これらの物理過程間の強い結合によって根本的に支配されます。境界での光学エネルギーの吸収は局所的な加熱とキャリア励起を引き起こし、これが媒体内で熱弾性応力や変位を誘発します。この結合機構は、温度、キャリア密度、力場が同時に進化し、それぞれの場が独立して振る舞うのではなく、互いに積極的に影響を与える現象を説明しています。

すべての物理量の空間的減衰は、半無限媒質内のエネルギー散逸と波の減衰の直接的な結果です。崩壊挙動は課された境界条件と物理的に整合しており、擾乱が励起領域付近に限定されるシステムの安定性を反映しています。これに対し、時間的進化は熱エネルギーとキャリアダイナミクスの連続的な相互作用により、場の振幅を増強し、これらは内部源としてシステムの応答を時間とともに強化します。これら...

材料

この記事で使用された材料の一覧

名前	会社	カタログ番号	コメント
計算ソフトウェア(記号解析および数値解析)	ウルフラム研究	Wolfram Mathematica(バージョン12.0)が使用されました	固有値計算、解析解の実装、数値評価に使用
データ可視化ツール(輪郭線およびヒートマップ生成)	ウルフラム研究	Wolfram Mathematica(バージョン12.0)は2D等高線プロットや時空ヒートマップの作成に使用されました	2D等高線プロットや時空間ヒートマップの作成に使用
材料パラメータデータセット(シリコン半導体特性)	様々な文献資料	該当なし	計算で使用される物理定数(弾性定数、熱定数、搬送波関連定数)(表1)
パーソナルコンピュータ/ワークステーション	HP	該当なし	計算は、数値シミュレーションに十分なメモリを持つWindows OSを搭載した標準的なパーソナルコンピュータで行われました
式エディタ	Microsoft WordとMathType	該当なし	原稿内の数学式のフォーマット作成や提示に使用
参考文献管理ソフトウェア	エルゼビア	該当なし	参考文献管理および引用のフォーマット作成(バンクーバースタイル)