10 GHz에서 골프공에 의한 전파 산란

1. 10 GHz에서의 전파 산란에 대한 소개

전자기파는 주변 공기와 전기적 특성이 다른 물체에 부딪힐 때 반사되거나 투과되며 산란되기도 한다. 물체의 크기와 파장 간의 상대적 관계도 중요한데, 크기가 작은 물체는 파장이 긴 전자기파에 거의 영향을 주지 않을 수 있으나 파장이 짧은 경우에는 복잡한 산란 패턴을 만들어낼 수 있다. 본 예제에서는 오목한 홈이 있는 골프공을 활용하여 전파가 불규칙한 형상을 가진 3차원 물체와 어떻게 상호작용하는지 분석한다.

우리가 사용하는 여기 신호의 주파수는 약 10 GHz이며, 진공 상태에서 이에 대응하는 파장은 대략 30 mm입니다. 규격화된 골프공의 직경은 약 42.7 mm로, 이는 진공 상태에서의 파장보다 약간 더 큽니다. 따라서 이 골프공은 마이크로파 산란 현상을 연구하기에 적합한 대상이라 할 수 있습니다. 골프공 표면의 딤플들은 기하학적 특징을 부여하지만, 각 딤플이 미치는 영향을 정확히 분석하려면 충분히 세밀한 수치적 메쉬가 필요합니다.

OghmaNano의 유한차분 시간영역법(FDTD) 해석기를 활용하여 전기장 및 자기장의 분포를 기록할 것입니다. 이후 해당 장들의 시간적 변화 양상, 전력밀도 분포, 그리고 검출기에 기록된 전력값을 검토할 예정입니다. 마지막으로 골프공을 찻주전자로 교체함으로써 동일한 RF 시뮬레이션 절차가 다른 삼각형 메쉬 구조를 가진 물체에도 적용 가능함을 보여줄 것입니다. 보다 상세한 이론적 배경은 FDTD 이론 및 스톤헨지 RF 튜토리얼을 참조하시기 바랍니다.

2. 골프공 RF 시뮬레이션 생성하기

OghmaNano를 실행한 후 ‘새 시뮬레이션’을 선택합니다. 사용 가능한 시뮬레이션 카테고리들은 ??에 표시됩니다. ‘RF 예제’ 카테고리(??)를 열고 ‘골프공 예제’를 선택하십시오. 쓰기 권한이 있는 디렉터리에 프로젝트를 저장합니다. 해당 예제에서는 기하학적 구조, 광원 및 검출기가 이미 설정되어 있으므로 별도로 모델을 구성할 필요 없이 바로 시뮬레이션을 실행할 수 있습니다.

3. 골프공 모델을 살펴보고 시뮬레이션을 실행하기

OghmaNano의 메인 창에는 계산 영역 내부에 빨간색으로 표시된 골프공이 나타난다(??). 이 빨간색 표면은 골프공과 그 딤플들을 나타내는 삼각형들로 구성되어 있다. 주변의 와이어프레임은 시뮬레이션 영역 및 측정 평면들을 보여준다. 사용자는 해당 기하학적 구조를 확대·회전시켜 RF 여기 신호에 대한 물체의 위치 관계를 확인할 수 있다.

‘시뮬레이션 실행’ 버튼(파란색 재생 아이콘)을 클릭하거나 F9 키를 누르십시오. 그러면 OghmaNano가 메쉬와 전자기장을 초기화한 뒤 FDTD 계산을 시간 순으로 수행합니다. 계산 진행 상황은 ?? 및 ??에 나타난 바와 같이 출력 콘솔에 표시됩니다. 계산이 완료되면 ‘출력’ 탭(??)을 열어주십시오.

4. 전자기장의 전파 현상을 관찰하기

시뮬레이션의 스냅샷 출력 결과를 열어 스냅샷 뷰어에서 power_density.csv을 확인하십시오. 슬라이더를 이용해 계산 과정을 단계별로 살펴보세요. ?? 및 ??에 수록된 이미지들은 전파가 시뮬레이션 영역 내를 진행함에 따라 전력밀도 패턴이 어떻게 변화하는지를 보여줍니다. 반복적으로 나타나는 밝고 어두운 띠들은 전자기 에너지 흐름의 공간적 구조를 나타내며, 골프공은 입사 복사선을 산란시켜 이러한 구조를 교란시키는 역할을 합니다.

초기 단계에서는 입사 전자기장이 물체에 도달하기 전의 상태를 파악하는 것이 유용합니다. 이후의 스냅샷들을 통해 산란장이 형성되는 과정 및 입사파와 재방향화된 파동 간의 간섭 현상을 확인할 수 있습니다. 색상 스케일은 W/m² 단위의 전력밀도를 나타냅니다. 특정 시점의 스냅샷과 나중에 분석되는 주파수 영역의 검출기 결과를 명확히 구분하는 것이 중요합니다.

5. 검출기의 전력값과 산란 스펙트럼을 확인하기

시각적 스냅샷 외에도 OghmaNano는 검출기에서 전자기장을 기록합니다. 전기장 및 자기장 벡터를 통해 순간적인 포인팅 벡터가 결정되며, 이 벡터는 해당 지점에서의 전자기파 에너지 흐름의 방향과 전력밀도를 나타냅니다.

\(\mathbf{S}(t)=\mathbf{E}(t)\times\mathbf{H}(t)\)

기록된 전력 값은 시간에 따라 또는 주파수 함수로 분석할 수 있다. 시간에 따른 검출기 출력 신호(??)는 연속파 여기가 시작되고 계산이 진행됨에 따라 신호가 어떻게 변화하는지를 보여준다. 산란 스펙트럼(??)은 10 GHz 부근에서의 검출기 반응을 요약한 결과이다.

물체가 없는 상태에서의 기준 시뮬레이션 결과가 존재할 경우, 해당 기준 시뮬레이션에서 얻은 검출기 전력 흐름 값을 골프공 시뮬레이션 결과에서 빼주면 차분 응답값이 산출된다 \(\Delta P(f)=P_{\mathrm{ball}}(f)-P_{\mathrm{air}}(f)\). 이러한 차감 처리는 특정 검출기에서 측정된 전력 흐름의 변화량을 나타내는 것이며, 이를 모든 방향으로 산란된 총 전력량으로 간주해서는 안 된다. 주파수 영역의 복소장의 경우, 시간 평균화된 전력밀도는 \(\tfrac{1}{2}\operatorname{Re}[\mathbf E\times\mathbf H^{*}]\)를 통해 계산된다.

6. 무선주파수 재료와 재료 데이터베이스

물체의 전기적 응답은 메쉬의 형태에 의해 결정되는 것이 아니라 그 물체에 할당된 재료에 의해 결정됩니다. 재료 데이터베이스에서 사용 가능한 RF 재료들을 살펴보십시오(??). 무선주파수 재료 데이터셋은 광학 분야에서 흔히 사용되는 나노미터 단위의 파장 대신 주파수(Hz, MHz 또는 GHz)에 따른 유전율 값을 기술합니다.

RF 데이터셋에는 상대 유전율의 실수부와 허수부가 포함되어 있으며, 일반적으로 채택된 부호 규칙에 따라 각각 \(\varepsilon_r'=\operatorname{Re}(\varepsilon_r)\) 및 \(\varepsilon_r''\)로 표기됩니다. 실수부는 파면의 위상 속도와 경계면에서의 반사 현상에 영향을 미치며, 허수부는 유전체 손실을 나타냅니다. 재료를 선택할 때는 해당 주파수 범위에 10 GHz가 포함되어 있고 허수부가 솔버에서 요구하는 부호 규칙을 따르는지 반드시 확인해야 합니다. 자세한 내용은 재료 데이터베이스 설명서를 참조하십시오.

7. 10 GHz 광원 및 검출기 설정 탐색하기

예제 시뮬레이션을 실행해보았으니 이제 여기서 사용된 여기화 설정들을 살펴보겠습니다. 광학 리본(??)을 열면 광원, 검출기 및 광학 메쉬 등 전자기파 시뮬레이션 도구들을 확인할 수 있습니다. 해당 설정들은 광학 분야 용어를 사용하고 있지만, 실제로는 RF 계산에도 동일한 FDTD 인프라가 활용됩니다.

광원 편집기를 열고 FDTD 탭(??)을 선택하십시오. 예시에서는 연속적인 사인(CW) 파형을 이용한 부드러운(가산적) 주입 방식이 사용되었습니다. 또한 진폭, 위상, 시작 및 종료 시간, 램프/윈도우 설정은 물론 Ex, Ey, Ez, Hx, Hy, Hz에 대한 독립적인 여기 스위치들도 제공됩니다. 스크린샷에서는 Ey만 활성화되어 있고 나머지 성분들은 비활성화된 상태입니다. 이러한 설정들은 전자기장이 어떻게 생성되는지, 특히 그 편광 특성이 어떠한지를 결정합니다.

검출기 편집기(??)를 열어보세요. Frequency Mesh 탭에서는 예시로 제시된 검출기가 9.50~10.50 GHz 범위에서 100개의 주파수 샘플을 사용하도록 설정되어 있습니다. 이를 통해 10 GHz에서 발생하는 여기 신호에 대응하는 주파수 샘플 집합을 얻게 됩니다. 검출기 샘플링 설정과 광원의 파형 설정은 서로 별개의 항목입니다. 따라서 검출기의 주파수 범위를 변경하더라도 시뮬레이션에 입력되는 복사 에너지 자체는 변하지 않습니다.

추가적인 배경 정보는 시간영역유한차분법 광원 및 시간영역유한차분법 검출기를 참조하시기 바랍니다.

8. 광학 메쉬와 시뮬레이션 크기 설정

수치적 메쉬는 OghmaNano가 맥스웰 방정식을 해석하는 데 사용되는 간격을 결정합니다. 광학 리본에서 광학 메쉬 편집기를 열어주십시오. ??에 해당 창이 나타납니다. 이 창에는 X, Y, Z 축 및 파장에 대한 별도의 제어 기능이 존재합니다. 제공된 스크린샷에서 주파수/파장 항목에는 약 10 GHz에 해당하는 단일 지점이 표시되어 있으며(편집기에는 10.0 GHz부터 10.10 GHz까지의 범위가 나타남), 이는 검출기의 9.50–10.50 GHz 범위에 걸친 샘플링 창과는 별개입니다.

공간 설정값들은 각 메쉬 세그먼트의 폭과 그에 할당되는 점의 개수를 지정합니다. 간격을 줄이면 골프공 및 전자기장의 해상도가 향상되지만 3차원 계산에 필요한 메모리와 계산 시간도 증가합니다. 일반적으로 사용되는 기준은 재료 내에서 가장 짧은 관련 파장을 여러 개의 메쉬 셀로 분해하는 것이며, 작은 요철이 있는 경우에는 더욱 세밀한 격자가 필요할 수 있습니다. 메쉬를 세밀화했을 때 결과가 안정적으로 유지되는지 반드시 확인해야 합니다.

또 다른 관련 설정 항목으로는 시뮬레이션 영역의 크기가 있습니다. ??에 나타난 바와 같이 Substrate XY size 버튼을 열면 치수 편집기가 나타납니다. 이곳에서 월드 오프셋 값 및 X, Y, Z축 방향의 총 치수를 변경할 수 있습니다. 여기서 지정한 치수는 광학 메쉬 편집기에서 설정된 전체 공간 범위와 반드시 일치해야 합니다. 그렇지 않을 경우 의도한 영역을 제대로 포괄하지 못하거나 불일치하게 샘플링되는 메쉬가 생성될 위험이 있습니다. 이 스크린샷들에 제시된 예시 값들은 편집기 사용법을 보여주기 위한 것이며, 튜토리얼 실행 시에는 실제 골프볼 프로젝트에 저장된 치수 값을 사용하시기 바랍니다.

9. 골프공을 찻주전자로 교체하기

동일한 RF 시뮬레이션 설정을 전혀 다른 물체에도 재사용할 수 있다. 3D 뷰 화면으로 돌아가서 골프공(??)을 마우스 오른쪽 버튼으로 클릭한다. 제공된 툴을 이용해 골프공의 기하학적 형상을 주어진 찻주전자 형태로 대체한다. 대체된 모델의 예시는 ??에서 확인할 수 있다. 이 찻주전자는 곡선형 몸체와 주둥이, 손잡이, 뚜껑을 갖추고 있어 기존 골프공과는 현저히 다른 형태를 보이며 산란 특성 또한 다르다.

가져온 객체가 시뮬레이션 영역 내부에 위치해 있는지 확인한 뒤 시뮬레이션을 다시 실행하십시오. 새로 적용된 기하학적 형상은 스냅샷 결과(??)에 나타납니다. 입사파가 주전자와 상호작용하며 반사 및 투과된 전자기 에너지가 이전과는 다른 공간적 분포를 보이는 것을 확인할 수 있습니다. 이는 솔버가 특정 내장형 형상에 종속되지 않음을 입증하는 것이며, 동일한 광원, 검출기 및 전자기장 시각화 절차가 일반적인 삼각형 메쉬에도 적용될 수 있음을 의미합니다.

10. 결과 해석 및 한계점

필드 플롯들은 입사하는 전자기파가 물체와 어떻게 상호작용하며 시뮬레이션 영역 전반에서 전력밀도 분포가 어떻게 변화하는지를 보여줍니다. 검출기로부터 얻은 데이터 역시 정량적 측정값을 제공합니다. 동일한 주파수, 광원, 검출기, 재료 설정 및 수치적 해상도 조건 하에서 골프공과 찻주전자를 비교함으로써 기하학적 형상 변화가 미치는 영향을 분리해낼 수 있습니다.

이들 모델은 단지 데모용이므로 그 결과를 실험적으로 검증된 레이더 단면적 값으로 간주해서는 안 됩니다. 예측된 산란 반응은 가져온 고체 물체에 할당된 유전체 또는 전도체 특성, 메쉬 해상도, 경계 조건 및 검출기 배치 위치에 따라 달라집니다. 특히, 도자기로 만들어진 속이 빈 주전자 모델과 달리 고체 유전체로 구성된 주전자 모델의 거동도 다를 수 있습니다. 정량적인 산란 결과를 얻기 위해서는 수렴성 검증을 수행하고 적절한 배경/기준 계산을 활용해야 합니다.

11. 유용한 출력 결과 및 추가 실험

출력 탭(??)에는 계산된 데이터가 포함되어 있습니다. 정확한 파일명은 예제 내용 및 OghmaNano 버전에 따라 달라질 수 있으며, 본 튜토리얼에서 사용된 파일들과 인터페이스에 대해 아래에 요약해 두었습니다.

더 큰 물리적 규모에서의 산란 현상을 살펴보려면 스톤헨지에서의 전파 산란을 참고하십시오. 시간영역유한차분법(FDTD) 기반의 광원 및 검출기 작동 원리에 대해 알아보려면 FDTD 광원과 검출기를 참조하시기 바랍니다.