유한차분 시간영역법(FDTD) 전자기 검출기
1. 서론
유한차분 시간영역법(FDTD) 시뮬레이션에서 검출기는 계산 영역 내에 배치되는 가상 측정 프로브 역할을 합니다. 이들은 광학 및 무선주파수(RF) 시뮬레이션 모두에 활용되는데, 파장과 응용 분야는 매우 다를 수 있음에도 전기장 및 자기장에 대한 측정 원리는 동일합니다. 이에 대한 적절한 비유로는 ??에 제시된 전압계가 있습니다. 전압계가 전기적 전위를 측정하는 것과 달리, FDTD 검출기는 시뮬레이션 내 특정 영역을 통과하는 전자기파를 통해 전자기장을 샘플링합니다.
검출기들은 FDTD 주파수 메쉬에 의해 결정된 지점들에서 전기장 값을 샘플링한다. 검출기는 시간에 따른 신호 변화를 추적하거나, 3차원 시뮬레이션에서는 표면 전체에 걸친 전자기 에너지 분포를 측정하는 데 활용될 수 있다. OghmaNano에서는 ??에 나타난 바와 같이 검출기 표면들이 3차원 뷰에서 격자 형태의 평면으로 표시된다. 분석 목적에 따라 측정된 전기장 값을 이용해 투과율, 반사율 또는 물체에 의해 산란된 전력 등을 계산할 수 있다.
FDTD 시뮬레이션이 실행되는 동안 검출기는 시간에 따라 변화하는 전자기장을 측정할 수 있습니다. OghmaNano는 복소 주파수 영역장을 직접 누적할 수도 있어, 전체 시간 데이터를 저장한 뒤 FFT를 적용할 필요 없이 원하는 주파수에서 스펙트럼 측정이 가능합니다. 이러한 특성 덕분에 동일한 검출기 프레임워크가 광학 소자, 안테나 및 RF 산란 문제 분석에도 적용될 수 있습니다.
검출기의 위치 지정 및 설정에 관한 일반적인 정보는 광학 검출기 항목을 참조하시기 바랍니다. 본 페이지에서는 FDTD 측정 방법과 그러한 측정값이 어떻게 주파수 영역 데이터 및 전력 값으로 변환되는지에 대해 중점적으로 설명합니다.
2. 검출기 설정
광학 리본 메뉴에서 ??에 나타난 검출기 편집기를 열거나, 3D 뷰에서 검출기 평면을 마우스 오른쪽 버튼으로 클릭한 뒤 ??에 나와 있듯이 ‘객체 편집’을 선택하십시오. 하나의 시뮬레이션에서는 여러 개의 검출기를 사용할 수 있는데, 예를 들어 소자 앞뒤에 각각 검출기를 배치하거나 RF 산란체를 향해 검출기를 설치하는 경우가 있습니다.
검출기 편집기의 모습은 ??에 나와 있습니다. 여기서 구성 탭을 통해 검출기의 기하학적 형태, 위치 및 정규화된 스펙트럼과 같은 옵션들을 설정할 수 있습니다. 참조값이나 입력 측정값을 기준으로 출력값을 정규화하는 작업은 광원 스펙트럼이 균일하지 않을 수 있는 투과 특성 연구에 특히 유용합니다. 검출기 표면을 통과하는 전력 측정 시에는 해당 표면의 방향 역시 중요합니다. 측정값의 부호는 검출기 법선 방향을 기준으로 결정됩니다.
?? 내의 ‘주파수 메쉬’ 탭에서는 복소장 정보가 누적되는 주파수 값들을 지정할 수 있습니다. 광학 응용 분야에서는 이러한 값들이 대개 파장 단위로 표현되며, RF 시뮬레이션에서는 MHz 또는 GHz 단위로 지정되는 경우가 많습니다. 이러한 주파수 샘플링 설정은 검출기의 파라미터이므로 기하학적 구조 및 전자기장을 해석하는 데 사용되는 공간적 FDTD 메쉬와 혼동해서는 안 됩니다.
3. FDTD 주파수 영역 누적기
FDTD 방식은 연속적인 시간 단계를 통해 전자기장을 계산하지만, 많은 유용한 소자 특성들은 주파수 함수 형태로 표현하는 것이 가장 적절합니다. 모든 시간 단계의 데이터를 저장한 뒤 처리하는 대신, 검출기는 시뮬레이션 진행 중에 사용자가 지정한 주파수에 대한 응답값을 산출할 수 있습니다. 이러한 작업은 복소 푸리에 누적기를 통해 구현됩니다.
시간 \(t_n=n\Delta t\)에 샘플링된 전기장 성분에 대해, 각주파수 \(\omega=2\pi f\)에서 누적된 값은
\[\widetilde{E}(\omega)=\sum_{n=0}^{N-1} E(t_n)\,e^{-i\omega t_n}\,\Delta t\]
각 시간 단계마다 지정된 주파수에 대한 복소값이 누적됩니다. 해당 주파수로 진동하는 성분은 일관되게 누적되는 반면, 다른 주파수의 성분들은 충분히 긴 관측 시간 동안 상쇄되는 경향을 보입니다. 물결표 기호는 푸리에 적분값을 나타내며, 아직 정규화된 정현파 형태의 전기장 진폭은 아닙니다. 동일한 계산을 전기장의 세 성분 \((E_x,E_y,E_z)\) 및 자기장 \((H_x,H_y,H_z)\)에 대해서도 수행할 수 있습니다.
해당 누적기는 전기장의 진폭과 위상 정보를 모두 보존한다. 이러한 차이점은 매우 중요하다. 즉, 투과 스펙트럼은 전기장 세기에 기반하여 계산될 수 있으나, 방향성 있는 전자기 전력의 계산에는 전기장과 자기장 간의 위상 관계가 필수적이기 때문이다. 충분히 긴 시간 \(T\) 동안 누적된 정상 상태의 연속파 응답에 대해 복소 피크 전기장 진폭은 대략 \(2\widetilde{E}/T\)가 되며, 이와 유사한 식이 \(\widetilde{H}\)에도 적용된다. 그러나 유한한 관측 시간 및 과도 현상이 이러한 근사값에 영향을 줄 수 있으므로 여기서의 누적 간격은 입력 신호의 특성에 맞게 적절히 설정되어야 한다.
요청된 각 주파수에 대해 별도의 FDTD 누적기가 유지되므로, 방대한 시간 영역 데이터를 저장하지 않고도 스펙트럼 정보를 얻을 수 있다. 주파수 점들은 주파수 메쉬 탭(??)에서 설정된다. 해당 수치적 방법에 관해서는 FDTD 이론을 참조하기 바란다.
4. 검출기 출력값
시뮬레이션이 완료된 후에는 OghmaNano의 출력 탭에서 검출기 결과를 직접 열어볼 수 있습니다. 사용 가능한 파일 목록은 활성화된 검출기 측정 항목에 따라 달라지며, 시간 영역 측정값, 주파수 영역 스펙트럼, 그리고 교정된 산란 시뮬레이션의 경우 산란 전력 데이터도 포함됩니다.
시뮬레이션 실행 후 검출기의 출력 결과는 출력 탭에 작은 CCD 형태의 아이콘으로 표시됩니다(??). 이 아이콘을 더블클릭하면 ??에 나타난 출력 파일들을 확인할 수 있습니다.
파일 lam_E.csv 및 lam_E_norm.csv는 주파수 영역 FDTD 누적기를 통해 생성됩니다. 첫 번째 파일은 원본 전기장 관련 스펙트럼을 기록하며, 두 번째 파일은 입력 측정값을 기준으로 정규화된 스펙트럼을 기록합니다. 이들 파일은 투과 특성이나 파장 선택성이 있는 소자들을 연구할 때 유용합니다. 관련 예시들은 ?? 및 ??에 제시되어 있습니다.
파일 power.csv에는 시간 영역에서의 검출기 신호가 포함되어 있으며 이는 ??에 도시되어 있다. 시간 영역 및 주파수 영역 결과는 동일한 전자기 시뮬레이션을 나타내지만, 주파수 영역 스펙트럼은 여기에 표시된 전력 곡선에 FFT를 적용하여 생성되는 것이 아니다. 대신 시뮬레이션 과정에서 전기장 및 자기장으로부터 직접 누적된 값들이며 이에 대해서는 3절에서 설명한다.
RF 산란 시뮬레이션에서는 scattered_power.csv 출력값도 생성될 수 있습니다. 이 값은 물체가 있을 때와 없을 때 측정된 총 전력의 차이값에 불과한 것이 아닙니다. 이 값은 주 시뮬레이션 결과와 기준 교정 시뮬레이션 결과를 비교하여 얻어진 복소 산란 전기장 및 자기장을 이용해 계산됩니다. 해당 측정 방법에 대해 아래에서 설명합니다.
5. 산란 전자기 전력 측정
전자기파가 골프공, 금속 물체 또는 기타 3차원 구조물에 부딪히면, 생성되는 전자기장은 원래 입사파와 물체에 의해 산란된 파의 합이 됩니다. 산란 현상을 측정하기 위해 OghmaNano는 물체가 없는 상태에서 교정 시뮬레이션을 수행하는데, 이때 광원, 계산 영역, 검출기 위치, 주파수 샘플링 및 누적 시간은 동일하게 유지됩니다. 그 후 물체가 존재하는 상태에서 또 다른 시뮬레이션을 실시합니다. 이러한 비교 과정을 통해 검출기에서 기록된 데이터로부터 입사장과 산란장을 분리해낼 수 있습니다.
각 검출기 위치 및 주파수에서, 교정된 전기장 성분들은 물체가 존재하는 경우에 기록된 전기장 성분들로부터 차감된다.
\[\widetilde{\mathbf E}_{\mathrm{sc}}=\widetilde{\mathbf E}_{\mathrm{object}}-\widetilde{\mathbf E}_{\mathrm{cal}}\qquad \widetilde{\mathbf H}_{\mathrm{sc}}=\widetilde{\mathbf H}_{\mathrm{object}}-\widetilde{\mathbf H}_{\mathrm{cal}}\]
여기서 \(\widetilde{\mathbf E}\) 및 \(\widetilde{\mathbf H}\)는 3절에서 설명한 복소 푸리에 적분값입니다. 필드의 크기가 아닌 복소 필드 자체를 빼줌으로써 필드의 위상을 보존할 수 있으며 이는 산란파를 분리해내기 위해 필수적입니다. 계산 과정에서는 6개의 직교 좌표계 필드 성분 모두가 사용됩니다.
전자기 에너지 흐름은 포이팅 벡터로 기술된다. 복소 정현파 형태의 장 진폭에 대해 그 시간 평균값은
\[\langle\mathbf S_{\mathrm{sc}}\rangle=\frac{1}{2}\operatorname{Re}\!\left(\mathbf E_{\mathrm{sc}}\times\mathbf H_{\mathrm{sc}}^{*}\right)\]
벡터 외적은 에너지 흐름의 방향을 나타내며, 복소 켤레 \((*)\)는 전기장과 자기장의 상대적 위상을 반영한다. 국부적인 포이팅 벡터의 단위는 W/m²이다. 안정된 정현파 형태의 장의 경우, 푸리에 적분값을 직접 활용하여 검출기 표면에 대해 적분하면
\[\boxed{P_{\mathrm{sc}}(f)\simeq\frac{2}{T^2}\sum_{k} \operatorname{Re}\!\left(\widetilde{\mathbf E}_{\mathrm{sc},k}\times \widetilde{\mathbf H}_{\mathrm{sc},k}^{*}\right) \cdot\hat{\mathbf n}\,\Delta A_k}\]
이 식에서 \(T\)은 필드 누적 시간을 나타내며, \(k\)는 검출기의 샘플링 지점을 지정합니다. \(\Delta A_k\)은 해당 지점이 나타내는 면적이고 \(\hat{\mathbf n}\)는 검출기 표면의 법선 벡터입니다. \(2/T^2\)는 누적된 전기장과 자기장을 각각 정현파 진폭으로 변환한 뒤 시간 평균 곱을 구함으로써 얻어지는 계수입니다. 검출기 표면 전체에 걸쳐 포이팅 벡터의 법선 성분을 합산하면 단위는 와트인 산란 전력 스펙트럼이 산출되며, 이 값은 scattered_power.csv에 저장됩니다.
💡 부호 해석 방법: 산란 전력은 검출기를 통과하는 유도된 플럭스 값이며, 필드 강도의 비음수 척도가 아닙니다. 음수 값은 순 산란 전력의 방향이 지정된 검출기 법선 방향과 반대임을 의미하며, 물체를 향해 배치된 검출기의 경우 이는 광원 쪽으로 되돌아가는 역산란 현상을 나타낼 수 있습니다. 또한 보고된 값은 전체 전력에서 공기만 존재할 때의 전력을 단순히 빼낸 결과와도 다릅니다. 이는 두 전력 값 모두에 입사-산란 간 간섭 항목이 포함되어 있으나, 필드 값을 먼저 빼낸 경우에는 그러한 항목이 존재하지 않기 때문입니다.
이 교정 방법은 주변 시뮬레이션 환경을 그대로 유지한 상태에서 다양한 형태의 물체로부터 발생하는 산란 현상의 주파수 의존성을 비교하는 데 특히 유용하다. 구체적인 예시로는 10 GHz에서 골프공에 의한 전파 산란 및 스톤헨지에서의 전파 산란을 참조할 수 있다. 신뢰성 있는 정량적 결과를 얻으려면 충분한 공간적 및 시간적 해상도, 일관된 교정 설정, 그리고 적절한 검출기 배치를 활용해야 한다.