교차빔 에너지 전달과 자극성 브릴루앙 산란을 위한 파동 기반 시뮬레이션 모델

초록

WEBS(웨브 기반 시뮬레이션) 코드는 슈뢰딩거형 포락선 방정식과 Du Fort‑Frankel 스킴을 이용해 교차빔 에너지 전달(CBET)과 자극성 브릴루앙 산란(SBS)을 동시에 풀어내는 효율적인 파동 기반 수치 모델이다. 단일 스킴으로 에너지 보존성을 유지하면서도 거친 격자에서도 정확한 결과를 얻으며, PIC과 기존 유체 이론을 검증한다. 고강도 레이저 조건에서 CBET 이득이 SBS 반사율 증가에 의해 억제되고, 두 빔 사이에 비대칭적인 SBS가 나타나는 등 상호 결합 현상을 규명한다.

상세 분석

WEBS는 레이저와 이온음파를 동일한 슈뢰딩거형 포락선 형태로 기술함으로써, 기존의 전자기 파동 방정식과 유체식 이온음파 방정식을 하나의 연립 시스템으로 통합한다. 이 접근법은 파동 진폭 a₀, a₁(펌프와 시드 레이저)와 밀도 변동 ˜nₑ를 복소수 포락선 변수로 정의하고, 시간 2차 미분을 무시한 비선형 1차 방정식(식 3‑5)으로 전개한다. 핵심은 Du Fort‑Frankel 스킴을 적용해 시간‑공간 격자에서 무조건적인 안정성을 확보하고, O(Δt²/Δx²) 정확도를 유지하면서도 에너지 누수를 최소화한다는 점이다. 특히, 스킴은 현재 시점 n, 이전 시점 n‑1, 그리고 반시점 n‑½을 동시에 사용해 업데이트를 수행하므로, 인접 격자값만 필요해 고도의 병렬화가 가능하고, 대규모 2‑D 혹은 3‑D 시뮬레이션에 적합하다.

코드 구현에서는 ξ₀, ξ₁ 등 복소수 계수를 통해 레이저 파동의 전파와 비선형 결합항을 정확히 반영한다. 이온음파 방정식(식 10‑11)에서는 플라즈마 흐름 Vₓ, V_y, Landau 감쇠 ν, 비선형 주파수 이동 δ ω 등을 모두 포함해 실제 실험 플라즈마 조건을 재현한다. 특히, 비선형 주파수 이동은 전자·이온 기여(αₑ)와 밀도 변동 비율 |˜nₑ/n₀|¹ᐟ²를 통해 동적 변조를 구현한다.

시뮬레이션 검증에서는 동일 파라미터를 갖는 PIC 코드(EPOCH)와 비교했을 때, 10 cells/λ₀ 격자에서도 펌프와 시드 레이저의 진폭 감소·증가 패턴이 거의 일치한다. 격자 해상도를 15, 20 cells/λ₀로 늘려도 결과 차이는 미미해, WEBS가 거친 격자에서도 물리적 정확성을 유지함을 보여준다. 또한, 플라즈마 흐름 속도 V_y를 0.5 Cₛ, 0.9 Cₛ, Cₛ 로 변화시켰을 때 CBET 전력 전환이 이론적 공명 조건(V_y≈Cₛ)에서 가장 크게 나타나는 것을 확인했다.

고강도 레이저(I≈3×10¹⁴ W/cm²) 조건에서 CBET와 SBS가 동시에 활성화되면, 초기 CBET에 의해 펌프 레이저 에너지가 시드 레이저로 전달되고, 이 과정에서 이온음파가 증폭돼 SBS 반사율이 급격히 상승한다. 결과적으로, SBS에 의해 손실된 에너지가 CBET 이득을 억제하고, 두 빔 사이에 비대칭적인 SBS 반사(한쪽 빔은 강하게, 다른 빔은 약하게)가 나타난다. 이는 기존 유체 기반 CBET 이론이 예측한 이득보다 30 % 이상 낮은 값을 보이며, 플라즈마 흐름과 Landau 감쇠가 약한 경우에 특히 두드러진다. 이러한 상호 결합 현상은 ICF 캡슐 대칭성 유지와 레이저 에너지 효율에 직접적인 영향을 미치므로, 실험 설계 단계에서 CBET와 SBS를 동시에 고려해야 함을 시사한다.

WEBS의 장점은 (1) 단일 스킴으로 에너지 보존성을 확보해 장시간 시뮬레이션이 가능, (2) 격자 해상도에 대한 민감도가 낮아 계산 비용을 크게 절감, (3) CBET와 SBS를 선택적으로 켜고 끌 수 있어 상호 작용 메커니즘을 명확히 분리해 분석 가능, (4) PIC 및 기존 유체 이론과의 정량적 일치를 통해 신뢰성을 검증받음이다. 향후에는 SRS와 TPD 같은 다른 레이저 플라즈마 불안정성을 추가 구현하고, 3‑D 비정상 플라즈마 흐름 및 비균일 밀도 구배를 포함한 보다 복잡한 ICF 환경을 모사할 계획이다.