스무스 입자 유체역학, 엄마가 알려줬다면 좋았을 것들

초록

본 논문은 스무스 입자 유체역학(SPH)의 핵심 특성을 정리하고, 전통적인 격자 기반 방법과의 차이점을 강조한다. 입자 기반 밀도 계산, 해밀토니안 유도 방정식, 정확한 대류와 무내재 점성, 전량 보존, 최소 에너지 상태 등 장점을 설명하고, 인공 점성·전도성 도입 필요성, 접촉 불연속면 처리, 텐서 불안정성 등 실용적 한계와 해결책을 제시한다.

상세 분석

SPH는 입자들의 질량을 고정하고, 커널 가중합을 통해 밀도를 추정한다는 가장 단순한 아이디어에서 출발한다. 이때 사용되는 스무딩 커널 W는 거리 |r–r_j| 에만 의존하므로 절대 위치·회전에 대해 불변이며, 따라서 밀도는 전역적인 좌표계와 무관하게 정의된다. 질량 보존은 입자 자체가 질량을 운반하므로 자동적으로 만족된다.

핵심은 밀도 합을 라그랑지안 L =∑ m_j(½v_j²–u_j) 에 삽입하고, 오일러–라그랑주 방정식을 적용해 운동 방정식을 도출한다는 점이다. 이 과정에서 얻어지는 가속도식은
d v_i/dt = –∑ m_j (P_i/ρ_i² + P_j/ρ_j²) ∇W_ij,
즉 연속 방정식의 정확한 해와 동일한 형태를 가진다. 여기서 파생되는 주요 특성은 다음과 같다.

1. 정확한 대류(Advection) – 입자 자체가 유체를 따라 움직이므로, 물질·에너지·자기장 등 모든 보존량이 수치적 확산 없이 정확히 이동한다. 그림 1의 자기 전류 루프가 10 000번을 넘는 도메인 횡단 후에도 초기 형태를 유지하는 것이 전형적인 예다.

2. 무내재 점성(Zero Intrinsic Dissipation) – 해밀토니안 시스템은 에너지와 엔트로피를 자체적으로 보존한다. 따라서 인공 점성 없이도 파동의 진폭이 감소하지 않는다(그림 1 오른쪽 Alfvén 파). 그러나 실제 물리적 충격파는 엔트로피가 증가해야 하므로, 인공 점성·전도성을 명시적으로 추가해야 한다.

3. 인공 점성의 역할과 스위치 – 전통적인 Monaghan 점성은 α c_s h 형태의 전단·체적 점성을 제공한다. 해상도에 따라 점성 계수가 선형적으로 변하기 때문에, 과도한 점성으로 인한 인공적인 확산을 방지하려면 Cullen & Dehnen(2010) 스위치와 같은 동적 제어가 필요하다.

4. 접촉 불연속면(Contact Discontinuity) 처리 – 압력은 연속이지만 밀도·온도가 불연속인 경우, 표준 SPH는 압력 ‘블립’ 현상으로 혼합을 억제한다. 인공 전도성 항을 도입하면 이 블립을 제거하고 Kelvin‑Helmholtz와 같은 전단 불안정을 정상적으로 재현한다(그림 2).

5. 전량 보존 – 라그랑지안이 변위·회전·시간에 대해 대칭성을 갖기 때문에, 질량·선형·각운동량·에너지·엔트로피가 모두 정확히 보존된다. 이는 복잡한 궤도 역학이나 비대칭 구조를 다룰 때 큰 장점이다. 다만, 보존성이 지나치게 강하면 수치적 잡음이 축적돼 입자 배치가 ‘유리상태’로 재정렬되는 과정에서 불안정성이 나타날 수 있다(그림 4).

6. 최소 에너지 상태와 텐서 불안정 – 해밀토니안은 입자들이 서로를 반발시켜 ‘유리‑같은’ 최소 에너지 배치를 형성하게 만든다. 이는 격자 재구성 없이도 입자 간 간격을 자동으로 유지하게 하지만, 압력이 음수인 경우(예: MHD의 텐서 힘)에는 ‘텐서 불안정’이 발생해 잡음이 급증한다. 비보존적이지만 안정적인 포뮬레이션을 선택하거나, 인공적인 안정화 기법을 적용해야 한다.

이와 같이 SPH는 물리적 보존과 대류 정확도에서 뛰어나지만, 불연속면 처리와 인공 점성·전도성 설계, 텐서 불안정 억제 등 실용적인 구현상의 고민이 반드시 동반된다.