SDPD 접근법을 통한 Navier‑Stokes 방정식 라그랑지안 탐구

초록

본 논문은 Smoothed Dissipative Particle Dynamics(SDPD) 모델이 라그랑지안을 가질 수 있는지 Helmholtz 조건을 이용해 검증한다. 유한 입자 수에서는 조건이 위배되어 전통적 Euler‑Lagrange 방식으로 유도할 수 없지만, 입자 수가 무한대로 갈 때는 모든 조건이 만족되어 연속 Navier‑Stokes 방정식의 라그랑지안을 원리적으로 도출할 가능성을 제시한다.

상세 분석

논문은 먼저 라그랑지안 존재 여부를 판단하는 수학적 도구인 Helmholtz 조건을 소개한다. 이 조건은 2차 미분항(가속도)과 1차 미분항(속도) 사이의 대칭성, 그리고 시간 미분 연산자를 포함한 일련의 편미분 관계식으로 구성된다. 저자들은 SDPD의 운동 방정식을 H_i 형태로 재작성하고, 각 변수 q_i를 입자 위치의 x, y, z 성분으로 정의한다. 첫 번째 조건(∂H_i/∂¨q_j = ∂H_j/∂¨q_i)은 입자 간 가속도 상호작용이 없으므로 자동으로 만족한다. 그러나 두 번째와 세 번째 조건은 복잡한 상호작용 텐서와 커널 함수 W(r)·F(r) 형태의 항들 때문에 위배된다. 구체적으로, 속도에 대한 편미분을 전개하면 A·Q_xij·…와 같은 비대칭 항이 남으며, 이는 0이 아니다. 이러한 잔차는 입자 수 N이 커짐에 따라 N^(-7/3) 혹은 N^(-2) 등 음의 거듭제곱 형태로 감소한다는 정량적 분석을 제공한다. 저자들은 Lucy 커널을 사용해 입자 밀도와 부피를 추정하고, K≈(4π/3)h^3N/V_0 관계를 통해 h∝N^(-1/3) 스케일을 가정한다. 이 가정 하에 모든 비대칭 항이 N→∞ 한계에서 사라짐을 보이며, 결국 연속적인 Navier‑Stokes 방정식은 Helmholtz 조건을 만족하는 라그랑지안 형태로 표현될 수 있음을 증명한다. 또한, 수치 시뮬레이션을 통해 N을 10^3에서 10^5까지 증가시켰을 때 조건 위반 항의 절댓값이 급격히 감소함을 확인한다. 이 결과는 이론적 분석과 일치하며, SDPD가 실제로는 연속 유체역학의 이산 근사임을 뒷받침한다. 논문은 또한 온도와 점성 등을 일정하게 가정함으로써 복잡성을 줄였지만, 이러한 단순화가 라그랑지안 존재 여부에 근본적인 영향을 주지는 않음을 암시한다. 전체적으로, 저자들은 SDPD가 미시적 입자 모델과 연속 유체 모델 사이의 다리 역할을 하며, 무한 입자 수 한계에서 최소 작용 원리를 적용할 수 있는 수학적 근거를 제공한다는 점에서 의미가 크다.