고유값 방법을 이용한 이방성 섬유 강화 실리콘에서 결합된 열-광탄성장의 시공간 분석

요약

본 연구는 이방성 섬유 강화 실리콘에서 결합된 열-광탄성장을 정상 모드 및 고유값 방법을 사용하여 분석합니다. 결과는 공간적 붕괴와 시간 의존적인 장의 진화, 그리고 이방성에 대한 강한 민감도를 보인다. 히트맵은 필드 분포와 위치 파악을 보여줍니다.

초록

본 연구는 이방성 섬유 강화 실리콘 반도체 매질에서 결합된 열-광탄성 시스템을 조사하여 열, 캐리어, 기계적 장 간의 상호작용을 포착하는 것을 목표로 합니다. 이러한 이방성 섬유 강화 재료는 마이크로전자 및 광전자 장치, 레이저 기반 기술, 센서, 첨단 복합 구조 등 현대 공학 응용 분야에서 방향성 특성과 향상된 기계적 성능이 요구되는 분야에서 중요한 역할을 합니다. 특히 열 및 광학 부하를 받는 반도체 부품 설계에서 중요하며, 결합장 거동의 정확한 예측이 신뢰성과 성능 최적화에 필수적입니다. 지배 방정식은 결합된 물리 모델을 기반으로 공식화되며, 분석을 단순화하고 관련 매개변수의 상대적 영향을 강조하기 위해 무차원 형태로 변환됩니다. 이 문제는 정규 모드 기법을 사용하여 1차 벡터-행렬 미분계로 환원한 뒤, 고유값 접근법으로 반무한 영역 내에서 부과된 경계 조건을 만족하는 해석적 해를 얻습니다. 수치 분석을 통해 시간 변화가 모든 물리적 장에 미치는 영향을 조사하며, 강한 공간 감쇠와 이방성 및 섬유 강화에 의해 지배되는 결합 거동이 드러난다. 시공간적 히트맵 표현은 장의 진화와 국소화를 시각화하는 데 사용되며, 다물리적 상호작용에 대한 물리적 통찰을 제공하고 분석 접근법의 효과를 입증합니다.

서론

열광탄성은 반도체 재료에서 광열 여기 하에서 열, 기계, 광 효과 간의 상호작용을 설명하는 중요한 다학제 분야로 부상했습니다. 이 필드 간의 결합은 레이저 가열과 광 여기가 관련된 현대 응용에서 특히 두드러집니다. 예를 들어, Saeed¹ 은 쌍곡선 2열 모델 사용으로 반도체의 열-광탄성 상호작용을 조사하여, 열 이완 효과가 시스템 행동을 정확히 예측하는 데 중요함을 입증했습니다. 섬유 강화 복합재의 열기계적 반응성도 향상된 기계적 강도와 이방성 특성 때문에 상당한 주목을 받고 있습니다. Li와 Lambros² 는 이러한 복합재의 동적 열기계적 거동을 분석하여 첨단 공학 응용에 적합함을 강조했습니다. 마찬가지로 Kalkal 등^은 자기장 하에서 회전하는 기능적 등급 섬유 강화 매질에서 2차원(2D) 변형을 연구하여, 이방성과 외부 효과가 시스템 반응에 강한 영향을 미친다는 것을 보여주었습니다. 더 나아가 Pitarresi 등은 열탄성 거동에서 거시적 이질성이 어떤 역할을 ^{하는지 살펴보며} 복합재료의 정확한 모델링의 필요성을 강조했습니다. 최근 발전은 이러한 연구를 결합된 다물리 효과를 가진 반도체 매체로 확장했습니다. Mondal 등은 메모리 반응과 자기장을 고려하여 강화 반도체의 파동 전파^{를 조사하여} 복잡한 결합 상호작용을 밝혀냈습니다. Akai 등⁶의 실험 연구들은 열탄성 온도 변화를 이용해 섬유 강화 복합재의 피로 손상을 평가하여 그 실용적 중요성을 확인했습니다. 또한, Lu 등7에 의해 입증된 바와 같이 미세역학적 접근법이 효과적인 열탄성 특성을 추정하는 데 사용^{되었습니다.} 외부 하중과 환경 조건의 영향도 널리 연구되어 왔습니다. Barak과 Dhankhar⁸ 은 기능적으로 등급화된 섬유 강화 매체에서의 경사 하중을 분석했고, Kundu와 Kalkal⁹ 는 중력 하와 이동하는 열 하중 하에서의 광열 상호작용을 조사했다. Chaudhary 등은 이중 상 지연 모델을 사용해 온도 의존적 ^{특성을 조사했고} , Pandit 등^은 비국소적 변형 거동을 포착하기 위해 분수 차수 변형 모델을 적용했습니다. 이 연구들은 열탄성 분석에서 현실적인 하중과 재료 조건을 고려하는 것의 중요성을 강조합니다.

역동적이고 진동 관련 현상도 반도체 시스템에서 탐구되어 왔습니다. Song 등^은 반도체 구조의 광열 진동을 연구했고, Mondal과^{Sur 13}은 메모리 효과를 가진 직교 매질에서의 파동 전파를 분석했습니다. 점탄성 및 미세구조 효과는 Abouelregal 등¹⁴에 의해 고려되었고, 경사로형 가열 효과는 Hobiny 등¹⁵에 의해 연구되었습니다. 더욱이, Zenkour¹⁶은 예측 정확도를 높이기 위해 삼상 지연 이론과 같은 일반화된 열탄성 모델을 개발했습니다. 전자기 및 미세 극성 효과를 포함한 고급 다물리 모델은 열광탄성 거동에 대한 이해를 더욱 크게 향상시켰습니다. Al-Hazaemh 등은 회전하는 반도체 매질에서의 광전자기-열탄성 여기^{를 연구했고}, Nazir와 Kumar¹⁸은 미세 극성 열탄성 상호작용을 분석했습니다. Song 등^은 비소산성 열탄성 상호작용도 연구했으며, Nasr와 Abouelregal²⁰은 공동이 있는 반도체에서 빛 흡수 과정을 조사했다. 비국소 및 분수 차수 모델은 최근 발전에서 중요한 역할을 해왔습니다. 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))을 따라 정렬되어 있어 매질의 기하학적 구성과 방향에 따른 이방성(anisotropy)을 보여줍니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

섬유 강화 이방성 열탄성 반도체 매질에서의 응력 텐서의 구성 관계는 방정식 1 ^1,5를 사용하여 일반 형태로 표현되었습니다. 이 공식에서 θ는 기준 온도 T₀에 대한 온도 증가량을 나타내며, T는 해당 시 절대 온도를 나타냅니다.

.   (1)

여기서 C_ijkl 은 탄성 강성 계수, e_kl 은 변형 텐서, β_ij와 η_ij 는 각각 열탄성 텐서와 캐리어 결합 텐서를 나타냅니다. 섬유 강화가 존재할 때, 재료 반응은 방향에 따라 달라지며, 섬유 방향 벡터 a = (a_i)에 의해 지배되었고, 이는 탄성 및 결합 항 모두에 이방성 기여를 도입했다. 이에 따라 구성 관계는 광섬유 강화의 효과를 명시적으로 포함하도록 확장^되었다:

. (2)

여기서 λ와 μ_τ는 라메 상수이고, μ_L은 섬유 방향을 따라 종방향 전단 계수입니다. 매개변수 α는 섬유 강화 효과를 나타내며, 열팽창 계수를 나타내는 α_ij와는 구별됩니다. 단위 벡터는 섬유 방향을 정의하고 응력-변형 반응에 방향 의존성을 도입했다. 현재의 2D 공식에서는 섬유가 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)로 표현되었다. 열-광탄성 매체에서의 표준 공식에 따라, 2차원에서의 운동 지배 방정식은 다음과 같이 작성되었습니다:

, (8)

. (9)

위의 방정식에 이방성 섬유 강화 구성 관계를 대입하면, 다음과 같이 결합된 편미분방정식(DE)계가 얻어졌습니다:

, (10)

. (11)

여기서 첨자는 공간 및 시간적 변수에 대한 부분 미분을 나타냅니다. 이 방정식들은 이방성, 섬유 강화, 온도 구배, 캐리어 확산이 매질의 동적 반응에 미치는 연합적 영향을 강조합니다. 광학 여기가 존재하는 상황에서 반도체 내부의 열장은 반송반자 밀도와 기계적 변형과의 상호작용에 의해 강하게 영향을 받아 완전 결합된 에너지 전달 과정을 이루게 되었습니다. 고전적 열 전도와 달리, 이러한 매질의 온도 변화는 캐리어 재결합과 열탄성 효과에서 비롯된 추가 원원 항에 의해 조절되어 열 전파 특성을 크게 변화시켰다. 일반화된 열탄성 체계 내에서 열 전도 방정식은 다음과 같이 표현되었다:^16,20:

. (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
K11	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
예g	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 결합된 이방성 광-열탄성 거동을 지배하는 물리적 및 재료 특성의 콤팩트 조합을 나타내기 위해 도입되었습니다. 각 계수는 시스템 내 특정 상호작용 메커니즘을 반영하며, 근본적인 물리적 과정의 상대적 영향력에 대한 통찰을 제공했습니다. 는 두 변위 성분 간의 이방성 상호작용 정도를 반영하는 정상 결합 강성과 주 탄성 강성의 비율을 나타냅니다. 는 횡방향 변형이 수직응력 성분에 대한 상대적 기여도를 특징짓습니다. 열탄성 결합의 방향 변화를 측정하여 열팽창 효과의 이방성을 나타냅니다. 운반체 밀도가 유도된 탄성 변형에 미치는 이방성 영향을 설명합니다. 정규화된 전단 강성을 나타내며, 전단 변형이 정상 변형에 비해 기여하는 바를 정량화합니다. 이는 지배 변위 방정식에서 수직 변형과 전단 변형의 결합을 설명합니다. 횡방향 강성과 전단 강성의 비율을 나타내며, 이방성 변형 거동을 강조합니다. 는 파동 전파 효과와 전단 강성을 연관시키는 정규화된 관성 매개변수를 나타냅니다. 서로 다른 공간 방향의 변위 구배 간의 결합을 특징지어. 횡방향 변위장에 대한 열 효과의 상대적 기여도를 정량화합니다. 캐리어에 의해 유발된 변형이 전단 강성에 미치는 영향을 측정합니다. 는 서로 다른 공간 방향에서의 열전도율 이방성을 나타냅니다. 매질 내 열 발생에 미치는 캐리어 재결합의 영향을 특성화합니다. 열 효과와 시간에 따른 탄성 변형 간의 결합을 나타냅니다. 이는 양방향 공간에서의 이방성 열팽창의 복합적 영향을 설명한다. 는 캐리어 수송 속도를 제어하는 정규화된 확산 매개변수를 나타냅니다. 는 캐리어 재조합 효과의 상대적 강도를 특징지어집니다. 열 변동과 캐리어 생성 과정 간의 결합을 설명합니다.

정상 모드 기법을 이용한 해석적 해법
결합된 이방성 열광탄성 시스템에 대한 해석적 해를 얻기 위해, 지배 부분 DE를 보다 다루기 쉬운 일반 DE 시스템으로 환원하는 효과 때문에 정상 모드 기법이 사용되었습니다. 이 접근법은 분산과 감쇠를 포함한 파동 전파 현상 분석에 널리 사용됩니다. 따라서 장 변수의 시간 및 횡단 공간 방향 모두에서 조화적 변화가 1^{, 12, 23으로 가정}되었다. 따라서 변위 성분, 온도, 운반체 밀도, 응력은 다음과 같이 지수 함수로 표현되었습니다:

. (22)

여기서 ω 는 장의 시간적 거동을 지배하는 복소 주파수를 나타내고, a 는 y-방향을 따라 공간적 변화와 관련된 파수입니다. 이 매개변수들은 안정성 요건을 충족하고 반무한 영역 내에서 물리적으로 허용 가능한 유계 해를 보장하기 위해 선택되었습니다. 위에서 가정한 형식들을 이전에 도출된 무차원 지배 방정식에 대입하고 그 결과의 표현식을 단순화함으로써, 원래의 부분 미정분방정분계 결합계는 공간 좌표 에 관한 일반 미정분방정분들의 계로 환원되었으며, 이는 다음과 같이 쓸 수 있다:

, (23)

, (24)

, (25)

. (26)

또한, 변환된 영역 내 대응하는 응력 성분은 다음과 같이 작성되었습니다:

, (27)

, (28)

. (29)

여기서 D는 미분 연산자 를 나타냅니다. 이 방정식들은 정규 모드 영역에서 지배 시스템의 축소형을 나타내며, 특성 방정식을 도출하고 이후 일반적인 해석 해를 구성하는 기초를 제공합니다. 계수는 다음과 같이 정의되었습니다: , , .

행렬 DE 공식화 및 고유값 분석
정규 모드 변환을 적용한 후, 방정식 23–26에 주어진 지배 체계는 공간 좌표에 대해 2차 일반 미방정분방정 집합으로 축소되었다. 체계적인 해를 용이하게 하기 위해, 이 시스템은 체의 1차 미분에 대응하는 보조 변수를 도입하여 동등한 1차 시스템으로 전환되었다. 구체적으로, 다음과 같은 변수들이 정의되었습니다:

, . (30)

이 정의를 바탕으로 23–26식 방정식은 8개의 1차 미정분방정식(DE)으로 구성된 다음 시스템으로 다시 작성되었습니다:

, (31)

, (32)

, (33)

, (34)

. (35)

위 시스템은 다음과 같이 콤팩트 행렬 A 형태로 표현되었습니다:

. (36)

상태 벡터는 다음과 같이 주어졌습니다:

. (37)

그리고 시스템 행렬은 다음과 같은 명시적 형태를 취했습니다:

. (38)

이 공식화는 원래 시스템을 고유값 문제로 전환시켰다 ^1,15. 특성 방정식은 다음에서 얻어졌습니다.

. (39)

이로 인해 고유값을 지배하는 8차 다항식이 생성된다. 환원형에서 특성 다항식은 다음과 같이 쓸 수 있다

. (40)

여기서 Z_i는 계수들의 함수이며 아래에서 명시적으로 정의된다. 결과적인 고유값은 감쇠 및 전파 특성을 포함한 해의 공간적 거동을 결정합니다. Re(m) > 0을 만족하는 고유값만이 유지되어, x → ∞에 따라 지수적으로 붕괴하는 물리적으로 허용 가능한 해를 보장합니다.

. (41)

특성 다항식의 근은 해의 공간적 거동을 지배하는 고유값 m을 정의한다. 이 고유값들은 Mathematica를 사용해 특성 다항식을 CharacteristicPolynomial 함수 로 구성하고 NSolve를 사용해 대수 방정식을 풀어 수치적으로 계산했습니다. 이 문제는 반무한 영역(x ≥ 0)에서 정식화되므로, x 와 → ∞ 사이에 상계를 유지하는 물리적으로 허용 가능한 해만 고려된다. 따라서 Re(m) > 0을 만족하는 고유값만 유지되어, x → ∞ 함수로 exp(−mx) 형태의 지수적으로 붕괴하는 해가 보장되었다. 나머지 근들은 모델의 물리적 요구사항과 일치하지 않는 붕괴 또는 무한한 해에 해당하므로 폐기되었다.

각 유지된 고유값 m에 대해, 대응하는 고유벡터는 관련 대수계로부터 얻어졌다

, (42)

다음과 같은 형태로 표현되었습니다:

. (43)

위 행렬 방정식을 확장하면 다음과 같은 선형 방정식 체계가 얻어졌습니다:

, (44)

, (45)

, (46)

, (47)

. (48)

고유값 문제의 동차성으로 인해, 고유벡터는 임의의 곱셈 상수까지 정의되었다. 유일하고 일관된 표현을 얻기 위해, 고유벡터의 한 성분을 고정하여 정규화 조건을 가했다. 본 연구에서는 첫 번째 성분이 q₁ = 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₄)입니다. ^T는 미지의 상수 벡터이며, θ₀과 같은 경계값에 의해 결정됩니다. N₀, 그리고 변위 제약 조건. 결과된 선형 시스템은 Mathematica를 사용하여 계산적으로 풀었으며, 계수 행렬과 오른쪽 베크터를 명시적으로 조립하고, 미지의 상수는 LinearSolve 루틴을 사용해 얻었습니다. 이 상수들은 다시 일반 해법에 대입되어 물리장의 완전한 표현식을 구성했으며, 이는 이후 결과의 수치 평가와 그래픽 표현에 사용되었다.

부과된 경계 조건은 다음과 같이 제시되었습니다:

온도 제약:

. (55)

이 조건은 주기적인 광 가열에 의해 고조적으로 변하는 표면 온도를 나타냅니다. 이는 매질 내에서 결합된 열탄성 및 운반체 수송 과정을 구동하는 주요 열 여기 역할을 합니다. 진폭 θ₀은 가해진 열 하중의 강도를 특징지어줍니다.

캐리어 밀도 제약:

. (56)

이 경계 조건은 광학 조명에 의해 생성된 광 생성 운반체 밀도를 설명합니다. 이는 광자 흡수에 의한 전자 여기와 입사 광장과 일치하는 조화 변조를 반영합니다.

변위 제약:

. (57)

이 조건은 경계가 횡방향으로 기계적으로 제약받음을 나타냅니다. 따라서 표면에서 v 방향을 따라 변위는 일어나지 않는다.

전단 응력 제약:

. (58)

이 조건은 전단 응력에 대해 무미끄러운 경계에 해당합니다. 이는 접선 방향에서 기계적으로 자유로운 경계와 일치하여 표면에 접선력이 작용하지 않도록 보장합니다. x = 0에서의 경계 조건 외에도, 무한대에서의 물리적 요구사항은 반무한 영역 내에서 유한한 물리적 해를 보장하는 것으로 부과 되었다. 수치 결과를 발표하기 전에, 본 연구에서 채택된 전체 계산 절차는 그림 2에 요약되어 있습니다. 계산에 사용되는 여기 매개변수 θ₀, N₀, 복소 주파수 ω, 파수 a의 수치값은 표 1에 나열되어 있습니다. 표 1 에 나열된 매개변수에는 정규화된 공식화에 사용되는 차원 재료 상수와 무차원 매개변수가 모두 포함됩니다. 수치 평가를 위해 공간 영역은 로 정의되었고, 횡좌표는 y = 0.6에 고정되었으며, 시간 영역은 로 간주 되었습니다. 이 범위들은 모든 수치 계산과 그래픽 표현에 사용되었습니다.

그림 2. 제안된 방법의 계산 워크플로우. 그림은 공식화에서 수치 결과까지의 단계들을 보여줍니다: 방정식 제어, 무차원화, 정규 모드 기법 적용, 1차 시스템으로의 변환, 행렬 공식화, 고유값 및 고유벡터 분석, 경계 조건 적용, 상수 결정, 수치 플롯 생성. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

결과

수치 결과
이 절에서는 이방성 섬유 강화 반도체 매질에서 결합된 열-광탄성 운반체 시스템의 거동을 분석하기 위한 수치 계산이 수행되었습니다. 고려된 재료는 실리콘(Si)이며, 그 물리적 및 재료 매개변수는 표 1에 나열되어 있습니다. 이 물질 상수는 지배 방정식에 직접 대입되어 수치 계산에서 장 변수를 평가하는 데 구현되었다. 모든 매개변수는 정의된 무차원 스케일링을 사용하여 일관되게 통합되었습니다. 수치 계산은 Mathematica의 기본 정밀도와 수렴 설정을 사용하여 수행되었습니다. 공간적, 시간적 방향 모두에서 0.01의 스텝 크기로 통일된 계산 그리드가 구현되었습니다.

첫째, 지배 방정식을 무차원 형태로 공식화하고 기본 장 변수로 표현했다. 정상 모드와 고유값 접근법으로 얻은 해석식은 Mathematica에서 기호 및 수치 평가를 ...

토론

얻어진 결과는 이방성 섬유 강화 반도체 매질에서 결합된 열-광탄성 거동에 대한 명확한 물리적 통찰을 제공합니다. 본 연구는 이러한 매질에서 열하중, 캐리어 생성, 탄성 변형 간의 상호작용을 조사하기 위한 고유값 기반 분석 틀을 제안합니다. 관찰되는 반응은 이 물리적 과정들 간의 강한 결합에 근본적으로 의해 지배됩니다. 경계에서 광학 에너지 흡수는 국소적인 가열과 캐리어 여기를 유발하며, 이는 매질 내에서 열탄성 응력과 변위를 유도합니다. 이 결합 메커니즘은 온도, 운반체 밀도, 기계적 장이 동시에 진화하는 현상을 설명하며, 각 장이 독립적으로 행동하기보다는 서로에 적극적으로 영향을 미치는 현상을 설명합니다.

모든 물리량의 공간적 감쇠는 반무한 매질 내에서 에너지 소산과 파동 감쇠의 직접적인 결과입니다. 붕괴 거동은 부과된 경계 조건과 물리적으로 일치하며, 교란이 여기 영역 근처에 제한되는 시스템 안정성을 반영합니다. 반...

재료

이 논문에 사용된 재료 목록

이름	회사	카탈로그 번호	댓글
계산 소프트웨어(기호 및 수치 분석)	울프람 연구	Wolfram Mathematica(버전 12.0)가 사용되었습니다	고유값 계산, 해석적 해법 구현, 수치 평가에 사용됩니다
데이터 시각화 도구(등고선 및 히트맵 생성)	울프람 연구	Wolfram Mathematica(버전 12.0)는 2D 등고선 플롯과 시공간 히트맵 생성에 사용되었습니다	2D 등고선 플롯 및 시공간 히트맵 생성에 사용됩니다
재료 파라미터 데이터셋(실리콘 반도체 특성)	다양한 문헌 출처	해당 없음	계산에 사용되는 물리 상수(탄성, 열, 반송파 관련) (표 1)
개인용 컴퓨터/워크스테이션	HP	해당 없음	계산은 수치 시뮬레이션을 위한 충분한 메모리를 갖춘 Windows OS를 실행하는 표준 개인용 컴퓨터에서 수행되었습니다
방정식 편집기	Microsoft  워드와 매스타입	해당 없음	원고에서 수학식 형식 작성과 제시에 사용됩니다
참고문헌 관리 소프트웨어	엘스비어	해당 없음	참고문헌 관리 및 인용 서식 작성에 사용됨 (밴쿠버 스타일)