스톤헨지에서의 전파 산란 현상
1. 무선주파수 산란 개요
무선주파수 파동은 그 크기와 물질적 특성이 전자기파의 파장과 유사한 물체들과 상호작용한다. 파동이 특정 구조물에 부딪히면 일부 에너지는 반사되고 일부는 물질 내부로 투과하며 또 다른 일부는 산란 및 회절 현상을 통해 방향이 바뀐다. 이러한 현상들은 안테나와 검출기에 수신되는 신호를 결정하며 무선파 전파, 레이더 및 무선 통신 분야에서 핵심적인 역할을 한다.
본 튜토리얼에서는 영국 남부에 위치한 선사시대 석조 기념물인 스톤헨지를 기하학적으로 복잡한 산란체의 대표적인 예시로 활용합니다. 이곳의 세워진 돌기둥들, 문턱석들, 쓰러진 돌덩이들 및 빈 공간들은 단일 구체나 직사각형 블록에 비해 훨씬 복잡한 무선주파수 환경을 형성합니다. 본 예시는 실제 현존하는 기념물에 대한 단순화된 수치 모델일 뿐, 고고학적 측량에 기반한 정밀 재현 모델은 아닙니다.
이 시뮬레이션에서는 OghmaNano의 3차원 유한차분 시간영역법(FDTD) 해석기를 활용하여 전기장 E와 자기장 H를 계산한다. 일반적으로 FDTD법은 광학 파장 영역에서 소개되지만 맥스웰 방정식은 무선주파수 영역에서도 적용 가능하다. 단, 물리적 치수, 여기 조건 및 재료 데이터가 해당 무선주파수 문제에 적합해야 한다. 배경 지식을 원한다면 FDTD법이란? 및 FDTD 소스를 참조하기 바란다.
여기서는 두 가지 상호보완적인 결과를 살펴보겠습니다. 즉, 시뮬레이션 영역 내에서 전자기장이 어떻게 형성되는지와 공기만으로 구성된 기준 계산 결과와 비교했을 때 스톤헨지가 검출기에서 측정되는 주파수 의존적 전력 흐름에 어떤 영향을 미치는지입니다. 이후 스톤헨지를 단순한 주거용 거리로 대체하여 또 다른 무선주파수 산란 환경을 예시로 보여줄 것입니다.
2. 스톤헨지 시뮬레이션 생성하기
OghmaNano를 실행한 후 ‘새 시뮬레이션’을 선택합니다. 그러면 ??에 표시된 템플릿 라이브러리가 나타납니다. ‘RF 예제’ 메뉴를 열어 ??에 나와 있는 Stonehenge 시뮬레이션을 선택합니다. 파일 기록이 가능한 디렉터리에 새 시뮬레이션을 저장합니다. 본 튜토리얼에서는 미리 준비된 예제를 활용하므로 Stonehenge의 기하학적 구조를 직접 생성할 필요가 없습니다.
💡 팁: 시뮬레이션 결과는 네트워크 폴더나 동기화 폴더가 아닌 로컬 드라이브에 저장하시기 바랍니다. FDTD 계산 과정에서는 다수의 중간 출력 파일이 생성될 수 있습니다. 자세한 내용은 시뮬레이션 가속화 방법을 참조하십시오.
??에 있는 예제 선택 화면에서는 RF 연구를 위한 기본 예제들이 표시된다.
3. 기하학적 구조를 살펴보고 시뮬레이션 실행하기
시뮬레이션이 생성된 후 메인 창에는 계산 영역 내의 스톤헨지가 표시된다(??). 붉은색 격자는 돌들을 나타내며, 주변의 격자 구조와 색상이 입혀진 평면들은 공간 격자 및 무선주파수 시뮬레이션 요소들을 의미한다. 시점을 회전시키거나 확대하여 서 있는 돌기둥, 삼각형 형태의 석조 구조물 및 쓰러진 돌덩이들을 살펴보도록 한다. 주요 석조 원형 구조의 직경은 약 30미터이며 그 너머로는 보다 멀리 떨어진 힐 스톤도 위치해 있다.
이 예제에서는 직사각형 블록 대신 다각형 수가 적은 표면 모델을 사용합니다. 불규칙한 형태의 돌 모양 덕분에 OghmaNano가 일반적인 3차원 삼각형 격자와도 호환된다는 점을 쉽게 확인할 수 있습니다. 돌의 형상은 해당 격자에 의해 정의되며, 재료의 특성은 유전체 재료 데이터를 통해 별도로 지정됩니다.
장면을 살펴본 후 ‘시뮬레이션 실행’ 버튼(파란색 재생 아이콘)을 클릭하거나 F9 키를 누르세요. OghmaNano가 계산 진행 상황을 ?? 및 ??를 통해 알려줍니다. 계산이 완료되면 ‘출력’ 탭(??)을 열어 결과를 확인하시면 됩니다. 제공된 예제를 실행하기 위해 메쉬나 재료 파라미터를 변경할 필요는 없습니다.
4. 전자기장의 전파 현상을 관찰하기
출력 브라우저에서 시뮬레이션의 스냅샷 결과를 열어보세요. 스냅샷 뷰어를 이용하면 저장된 전자기장 상태들 간을 이동하며 도메인 전반에 걸쳐 전파가 어떻게 전파되는지 확인할 수 있습니다. 그림 ??, ??, ??은 이 과정의 세 단계를 나타냅니다. 전자기장이 석조 구조물에 부딪혀 통과함에 따라 그 분포는 점점 더 복잡해집니다.
이러한 스냅샷들은 여기서 발생하는 여기 신호가 의도한 영역으로 전파되는지, 그리고 시뮬레이션이 산란체에 어떻게 반응하는지를 확인하는 데 유용하다. 교차되는 전기장 패턴은 공간적 간섭 현상 및 파동의 위상을 반영하지만, 단일 스냅샷만으로는 돌들에 의해 반사되는 전력량을 분리해낼 수는 없다.
5. 무선주파수 대역에서의 재료 특성
스톤헨지는 완벽한 금속 반사체가 아닌 유전체 구조물로 모델링된다. 단순한 건조석 구조물을 예시로 설명할 때는 상대 유전율 값이 약 5인 경우가 설명에 유용하나, 실제 석재의 특성은 그 구성 성분, 수분 함량 및 주파수에 따라 달라진다. 무선주파수 영역에서는 상대 유전율의 실수부와 허수부 모두가 전파 전파에 영향을 미칠 수 있다.
??에 제시된 재료 데이터베이스 영역 및 ??의 무선주파수 재료 뷰를 열어 주파수에 따른 전자기 특성 데이터를 확인하십시오. 파장에 따라 표시되는 광학 데이터셋과 달리 무선주파수 데이터는 일반적으로 Hz, MHz 또는 GHz 단위의 주파수에 따라 표시됩니다. OghmaNano는 상대 유전율의 실수부와 허수부인 εr = ε'r + iε''r 값을 저장하며 이는 해석기의 부호 규칙에 따릅니다. 유전율의 실수부는 물질 내에서의 위상 속도와 파장에 주로 영향을 미치는 반면, 허수부는 해당 부호 규칙에 따라 유전 손실을 나타냅니다. 전도성 물질의 경우 저주파 영역에서 매우 큰 유효 허수 상대 유전율 값을 보일 수 있습니다.
예를 들어, 어떤 재료가 전기 전도도 σ로 표현될 경우 전도도와 관련된 허수 상대 유전율의 크기는 \[|\varepsilon_r''(f)|=\frac{\sigma}{2\pi f\varepsilon_0}.\]입니다. 따라서 주파수가 낮아질수록 이 값은 증가하게 됩니다. 금속의 사례는 ??에 제시되어 있습니다. 다만 스톤헨지에 대해서는 금속 곡선을 사용해서는 안 되며, 돌 재료의 경우 별도의 유전체 재료 데이터가 필요합니다.
참고로 30 MHz에서 진공 상태에서의 파장은 약 10 m입니다. εr값이 약 5인 손실이 적은 유전체 내에서는 파장이 약 4.5 m로 줄어듭니다. 따라서 미터 크기의 석조물들은 전자기파의 전파 패턴을 변화시키기에 충분한 크기를 갖게 됩니다. 기본적인 전자기 방정식에 대해서는 FDTD theory를 참조하시기 바랍니다.
6. 시뮬레이션 영역 및 메쉬
수치 해석 영역에는 기념물, 입사파 및 검출기 위치가 포함되어야 하며, 전자기장이 충분히 전파될 수 있도록 넉넉한 공간이 필요합니다. 시뮬레이션 크기 설정 창은 ??에 나타나 있습니다. 메쉬 상태는 ??에 있는 메쉬 제어 기능을 통해 확인할 수 있습니다. 처음에는 예제에서 제공된 기본 설정값을 사용하시고, 정밀도 검증이나 경계 조건의 영향을 조사하고자 할 때에만 영역 크기나 메쉬 밀도를 증가시키기 바랍니다.
격자 간격은 유전체 내부의 전자기파 파장과 가져온 석조 구조물의 주요 형상 특징들을 모두 충분히 표현할 수 있어야 한다. 격자가 미세할수록 기하학적 형상 및 전자기장 해상도가 향상되지만 3차원 FDTD 계산에 필요한 메모리 용량과 연산 시간도 증가한다. 단지 하나의 격자만 사용한다고 해서 수렴성이 보장되는 것은 아니다. 정량적 연구를 위해서는 더 미세한 격자를 사용해 시뮬레이션을 반복 수행한 뒤 검출기에서 얻어지는 스펙트럼의 변화가 미미한지 확인해야 한다.
경계 조건은 개방 공간에서의 산란 계산 시 특히 중요하다. 흡수 경계는 메쉬 가장자리에서 발생하는 인위적인 반사파를 최소화해야 한다. 관련 출력 결과는 ??에 나타나 있다. 이 결과의 정확한 해석은 표시 대상으로 선택된 물리량에 따라 달라진다.
💡 팁: 가장 유용한 메쉬 검증 방법은 색상 플롯의 외관만으로 수렴성을 판단하는 대신, 점차 세밀해지는 격자들에서 동일한 검출기 관측값을 비교하는 것입니다.
7. 주파수에 따른 산란 전력 측정
전자기장의 시각화는 유용하지만, 무선주파수 시뮬레이션에서는 정량적 결과도 필요하다. 본 예제에서는 전기장 및 자기장의 x, y, z 성분을 모두 기록하는 검출기를 사용한다. 이 정보를 통해 단위 면적당 전력 흐름을 나타내는 전자기 포인팅 벡터를 산출할 수 있다. 순간적인 전자기장의 경우 \[\mathbf{S}(t)=\mathbf{E}(t)\times\mathbf{H}(t).\] 주파수 영역에서는 복소 위상자에 기반한 시간 평균 전력 흐름 밀도가 \[\langle\mathbf{S}(f)\rangle=\tfrac{1}{2}\operatorname{Re}\!\left[\mathbf{E}(f)\times\mathbf{H}^{*}(f)\right].\] 된다. 검출기와 관련된 물리량은 바로 이 벡터의 검출기 평면에 수직인 성분이며 그 부호는 검출기 방향 설정에 따라 결정된다. 검출기 출력값에 관한 추가 정보는 FDTD detectors를 참조하기 바란다.
기념물이 미치는 영향을 확인하려면 그 외 조건이 동일한 두 가지 계산 결과를 비교해야 한다. 하나는 스톤헨지가 포함된 경우이고 다른 하나는 공기만 존재하는 경우이다. 이때 검출기에서 측정되는 전력 차이는 \[\Delta P(f)=P_{\mathrm{Stonehenge}}(f)-P_{\mathrm{air}}(f).\]이며 각각의 \(P(f)\) 값은 전기장 및 자기장 데이터를 바탕으로 적절한 포인팅 벡터 성분(필요 시 검출기 표면에 대한 적분도 포함)을 이용해 산출된다. 음수 값을 나타내는 차이는 기념물이 존재할 때 검출기가 설정한 양의 방향으로 흐르는 순 전력량이 감소했음을 의미한다. 그러나 이 결과만으로는 소실된 에너지가 돌에 의해 흡수되었다고 단정할 수 없다. 반사, 방향 전환 및 간섭 현상 모두 측정값에 영향을 줄 수 있기 때문이다.
산출된 주파수 영역 스펙트럼이 ??에 나타나 있다. 뚜렷하게 나타난 스펙트럼 특성은 기념물의 영향이 주파수에 따라 달라짐을 보여준다. 이에 대한 해석은 입력 스펙트럼, 검출기의 위치 및 방향, 그리고 선택된 정규화 방식에 따라 달라지며, 이를 단순히 구조적 공명 현상으로 간주해서는 안 된다.
💡 해석 방법: 전력은 방향성을 가집니다. 검출기의 법선 방향이 반전될 경우, 부호가 붙은 포인팅 벡터 성분도 함께 반전됩니다. 결과를 비교할 때는 검출기의 위치, 면적, 법선 방향 및 기준 시뮬레이션 조건을 일관되게 유지해야 합니다.
8. 스톤헨지에서 주거용 거리까지
동일한 접근법은 다양한 무선주파수 기하구조에도 적용 가능하다. 두 번째 예시로, 스톤헨지 대신 단순화된 영국 주택가 거리를 사용해보자. 이 장면은 여러 채의 주택들과 경사진 지붕들, 그리고 작은 기하학적 요소들로 구성되어 있다(??). 스톤헨지와 마찬가지로 이 역시 해석기의 성능을 시연하기 위한 저다각형 기하구조이며, 실제 건축 모델의 세부 정보까지 반영한 것은 아니다. 해당 거리는 대략 동일한 크기의 시뮬레이션 영역을 차지하도록 설계되어 있어 유사한 여기 신호 및 검출기 배치를 그대로 재사용할 수 있다.
주택가 거리의 기하학적 구조를 검토하거나 변경하려면 3차원 뷰(??)에서 해당 거리 객체를 마우스 오른쪽 버튼으로 클릭한 뒤 기하학적 구조 및 격자 제어 기능(??)을 열어주십시오. 이러한 편집 도구들은 초기 시뮬레이션 실행 전 스톤헨지 구조를 수정해야 했던 상황과 달리, 현재는 다른 외부 모델을 다루고 있기 때문에 소개되는 것입니다. 관련 기하학적 조작 방법에 대해서는 step editor 설명서를 참조하시기 바랍니다.
이 입문용 비교 실험에서는 건물들에 단일 유전체 재료를 할당할 수 있습니다. 따라서 이 건물들은 내부가 공기로 채워져 있고 유리창과 여러 층으로 된 지붕 및 벽돌벽이 있는 실제 주택이 아니라 단순화된 고체 유전체 구조물로 이해하는 것이 적절합니다. 그럼에도 불구하고 이 모델들을 통해 실제 건축 환경에서 전자기장의 전파, 반사 및 공간적 그림자 현상을 연구할 수 있습니다.
거주지 도로 모델 시뮬레이션을 실행한 후 ??에 나타난 정상상태 전력밀도 결과를 열어보면, 건물들 주변 및 사이에 전력밀도가 높거나 낮은 영역들이 존재함을 알 수 있습니다. 이러한 공간적 분포는 장애물들이 배열된 경우 무장애 공간과 비교하여 전자기장이 어떻게 변화하는지를 보여줍니다.
거리 시뮬레이션과 스톤헨지 시뮬레이션은 각각 다른 물리적 목적을 갖고 있지만, 둘 다 동일한 유한차분 시간영역법(FDTD) 방정식을 기반으로 합니다. 일반적인 삼각형 격자 기하학 구조를 활용할 때의 장점 중 하나는 격자 및 유전체 재료의 특성이 정의되면 전자기 해석기를 다양한 무선주파수 문제에 적용할 수 있다는 점입니다.
9. 결과 해석 및 한계점
이러한 시뮬레이션들은 무선주파수 파동이 복잡한 3차원 유전체 구조물들과 어떻게 상호작용하는지를 보여준다. 스톤헨지 예시에서는 불규칙하게 배치된 개별 석조물들에 의한 파동 산란 현상을 강조하며, 도심가 예시에서는 보다 큰 건물들이 반복적으로 배치된 환경에서의 파동 전파 현상을 강조한다. 두 경우 모두 시각화 결과를 통해 전력 밀도의 분포 변화를 확인할 수 있으며, 검출기 스펙트럼을 통해 주파수에 따른 시스템 응답의 변화 양상도 측정할 수 있다.
몇 가지 중요한 제약사항이 존재한다. 실제 건축 재료들은 불균질하며 스톤헨지의 기하학적 모델 역시 근사치에 불과하고, 공간 격자는 메쉬 간격보다 충분히 큰 구조만을 해석할 수 있다. 따라서 본 시뮬레이션 결과는 실제 스톤헨지나 특정 거리에서 측정된 무선주파수 특성을 검증한 예측값이 아니라 단지 시뮬레이션 데모로 간주되어야 한다. 정량적인 공학적 분석을 위해서는 유전체 파라미터를 검증하고 메쉬 수렴성을 확인하며 경계면 반사 현상을 점검한 후 기존의 벤치마크 데이터나 실측값과 비교해야 한다.
축하합니다! 사용자는 3차원 무선주파수 유한차분 시간영역법(FDTD) 시뮬레이션을 수행하고, 전자기장의 전파 현상을 관찰했으며, 주파수 영역에서 방향성 전력 흐름의 변화를 분석했습니다. 또한 기하학적 구조를 변경함으로써 서로 다른 전자기 환경이 형성되는 과정도 확인하였습니다.
10. 유용한 시뮬레이션 출력 결과
출력 탭에서는 스냅샷, 검출기 측정값 및 기타 해석기 결과에 접근할 수 있다. 실제 파일 목록은 시뮬레이션 설정에 따라 달라지나, 무선주파수 산란 현상을 연구할 때 유용한 항목들은 다음과 같다. 출력 브라우저의 모습은 ??에 제시되어 있다.
| 출력 결과 | 목적 |
|---|---|
| 전자기장의 순간 이미지 | 시간에 따른 전자기장의 전파 및 산란 현상을 표시합니다. |
| 전기장 성분 | 전기장 계산 및 전력 계산을 위해 Ex, Ey, Ez 값을 제공해야 한다. |
| 자기장 성분 | 포인팅 벡터를 계산하기 위해 Hx, Hy, Hz 값을 제공해야 한다. |
| 검출기 스펙트럼 | 지정된 검출기 평면을 통과하는 주파수 의존적 전력 흐름을 측정한다. |
| 공기만을 사용한 기준 조건 | 차분 수신 전력 계산의 기준값을 제공한다. |
| 정상상태에서의 전기장 분포 또는 전력 분포도 | 물체 주변에서 전자기 에너지의 공간적 재분포 현상을 보여주십시오. |