19.9 나비에-스토크스 방정식의 유도와 의미

1. 개요

나비에-스토크스 방정식은 점성 유체의 운동을 기술하는 운동량 보존 법칙의 미분 형태이며, 고전 유체역학의 지배 방정식이다. 이 방정식은 클로드-루이 나비에가 1822년에 분자 논의를 통해 처음 제시하였고, 조지 가브리엘 스토크스가 1845년에 연속체 논의를 통해 엄밀하게 재유도하였다. 나비에-스토크스 방정식은 뉴턴 유체의 응력-변형률 관계에 기반하며, 점성 항을 통해 유체의 내부 마찰과 운동량 확산을 기술한다. 본 절은 응력 텐서의 구성, 구성 방정식의 유도, 방정식의 일반 형태와 비압축성 형태, 각 항의 물리적 해석, 경계 조건, 에너지 소산, 해의 존재성과 정규성 문제, 로봇 공학적 의의를 학술적으로 정리한다.

2. 코시 운동량 방정식

유체의 운동량 보존 법칙은 연속체 역학의 일반적 결과인 코시 운동량 방정식으로부터 출발한다.

\rho\,\frac{D\mathbf{u}}{Dt} = \nabla\cdot\boldsymbol\sigma + \rho\,\mathbf{b}

여기서 \boldsymbol\sigma는 코시 응력 텐서, \mathbf{b}는 단위 질량당 체적력이다. 이 식은 임의의 연속체에 적용되며, 유체와 고체의 구분은 응력 텐서와 변형률 사이의 구성 관계에 의해 결정된다. 유체의 정의에 따라 응력 텐서는 정지 상태에서 등방성을 가지며, 운동 상태에서는 변형률 속도에 의존하는 점성 응력이 추가된다.

응력 텐서의 분해

유체의 응력 텐서는 등방성 부분과 편차성 부분으로 분해된다.

\sigma_{ij} = -p\,\delta_{ij} + \tau_{ij}

여기서 p는 열역학적 압력, \tau_{ij}는 점성 응력 텐서이다. 음의 부호는 압력이 수축 방향으로 작용함을 나타내며, \delta_{ij}는 크로네커 델타이다. 이 분해는 유체 응력의 물리적 본질을 반영하며, 정지 상태에서는 점성 응력이 영이 되어 압력만 남는다.

3. 뉴턴 유체의 구성 방정식

뉴턴 유체의 가정에 따르면 점성 응력은 변형률 속도 텐서에 선형적으로 의존하며, 등방성 유체의 경우 구성 관계는 다음과 같이 표현된다.

\tau_{ij} = 2\mu\,S_{ij} + \lambda\,(\nabla\cdot\mathbf{u})\,\delta_{ij}

여기서 S_{ij}는 변형률 속도 텐서이며 다음과 같이 정의된다.

S_{ij} = \frac{1}{2}\left(\frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i}\right)

\mu는 동역학 점성 계수(전단 점성), \lambda는 제2 점성 계수(체적 점성)이다. 스토크스 가설은 3\lambda + 2\mu = 0, 즉 \lambda = -2\mu/3의 관계를 제안하며, 이는 대부분의 공학적 유동에서 만족되는 것으로 간주된다. 이 가정 하에서 점성 응력의 대각합은 영이 되며, 열역학적 압력과 평균 수직 응력이 일치한다.

4. 나비에-스토크스 방정식의 일반 형태

응력 텐서의 구성 관계를 코시 운동량 방정식에 대입하면 나비에-스토크스 방정식의 일반 형태가 유도된다.

\rho\,\frac{D\mathbf{u}}{Dt} = -\nabla p + \nabla\cdot(2\mu\,\mathbf{S}) + \nabla(\lambda\,\nabla\cdot\mathbf{u}) + \rho\,\mathbf{b}

점성 계수가 공간적으로 일정한 경우 이 식은 다음과 같이 단순화된다.

\rho\,\frac{D\mathbf{u}}{Dt} = -\nabla p + \mu\,\nabla^2\mathbf{u} + (\lambda + \mu)\,\nabla(\nabla\cdot\mathbf{u}) + \rho\,\mathbf{b}

이 형태는 압축성 뉴턴 유체의 완전한 운동량 방정식이며, 연속 방정식과 에너지 방정식과 결합되어 유동의 지배 방정식 체계를 구성한다.

5. 비압축성 나비에-스토크스 방정식

비압축성 유동의 경우 연속 방정식으로부터 \nabla\cdot\mathbf{u} = 0이 성립하며, 이로 인해 체적 변형 항이 소멸된다. 결과로 얻어지는 비압축성 나비에-스토크스 방정식은 다음과 같이 표현된다.

\rho\,\frac{D\mathbf{u}}{Dt} = -\nabla p + \mu\,\nabla^2\mathbf{u} + \rho\,\mathbf{b}

또는 동역학 점성 \nu = \mu/\rho를 사용하여 다음과 같이 쓸 수 있다.

\frac{\partial\mathbf{u}}{\partial t} + (\mathbf{u}\cdot\nabla)\mathbf{u} = -\frac{1}{\rho}\nabla p + \nu\,\nabla^2\mathbf{u} + \mathbf{b}

이 형태는 액체 유동과 저속 공기 유동의 해석에 광범위하게 적용되며, 로봇 공학의 유체 응용의 대부분에서 표준 지배 방정식으로 사용된다.

6. 각 항의 물리적 해석

나비에-스토크스 방정식의 각 항은 명확한 물리적 의미를 가진다. 국소 가속도 \partial\mathbf{u}/\partial t는 고정된 점에서의 속도의 시간 변화율을 나타낸다. 대류 가속도 (\mathbf{u}\cdot\nabla)\mathbf{u}는 유체 입자가 공간적으로 변화하는 속도장을 따라 이동할 때의 가속을 기술하며, 방정식의 비선형성의 근원이다. 압력 구배 항 -\nabla p/\rho는 압력의 공간적 변화에 의한 힘을 나타낸다. 점성 항 \nu\nabla^2\mathbf{u}는 유체 분자의 무질서 운동에 의한 운동량 확산을 기술하며, 열 확산과 수학적으로 유사한 형태를 가진다. 체적력 항 \mathbf{b}는 중력과 같은 외력을 나타낸다. 이러한 물리적 해석은 유동 현상의 직관적 이해를 제공한다.

7. 점성과 운동량 확산

점성 항은 운동량이 확산 과정을 통해 공간적으로 전달되는 메커니즘을 기술한다. 점성 계수 \mu는 전단 변형에 대한 응력 응답의 척도이며, 온도에 의존한다. 액체의 점성은 온도가 증가함에 따라 감소하는 반면, 기체의 점성은 증가하며, 이는 분자 이동론으로 설명된다. 점성에 의한 운동량 확산은 유체 내부의 속도 경사를 평활화하며, 이는 경계층 현상과 층류의 안정성에 직접 연결된다.

8. 경계 조건

나비에-스토크스 방정식의 해를 결정하기 위해서는 경계 조건이 필요하다. 고체 표면에서의 표준적 조건은 점착 조건(no-slip condition)이며, 유체의 속도가 고체 표면의 속도와 완전히 일치함을 요구한다.

\mathbf{u} = \mathbf{u}_{\text{wall}}

이 조건은 실험적으로 잘 검증된 경험 법칙이며, 이상 유체의 비관통 조건과 구별된다. 점착 조건의 결과로 벽면 근처에서 속도 경사가 급격히 변화하는 경계층이 형성된다. 자유 표면, 입구, 출구에서는 각각 응력 연속 조건, 속도 분포 조건, 압력 조건이 적용된다. 원거리 조건은 무한 영역의 해석에서 속도와 압력의 점근적 거동을 지정한다.

에너지 소산과 비가역성

점성 유동에서는 역학적 에너지가 열로 변환되는 비가역적 에너지 소산이 발생한다. 단위 부피당 점성 소산율은 다음과 같이 표현된다.

\Phi = \tau_{ij}\,\frac{\partial u_i}{\partial x_j} = 2\mu\,S_{ij}S_{ij} + \lambda\,(\nabla\cdot\mathbf{u})^2

이 양은 항상 비음이며, 열역학 제2법칙과 일치한다. 점성 소산은 유동의 운동 에너지를 내부 에너지로 변환하며, 유체의 온도 상승과 관련된다. 이러한 비가역적 효과는 이상 유체의 가역적 모델과 구별되는 점성 유동의 본질적 특성이다.

9. 해의 존재성과 정규성

3차원 비압축성 나비에-스토크스 방정식의 매끄러운 해의 존재성과 유일성은 수학적으로 완전히 해결되지 않은 문제이며, 클레이 수학연구소의 밀레니엄 문제 중 하나로 지정되어 있다. 레르이가 1934년에 약한 해(weak solution)의 존재성을 증명하였으나, 3차원의 경우 고전적 해(classical solution)의 전역적 존재성과 정규성은 미해결이다. 2차원의 경우에는 모든 초기 조건에 대해 매끄러운 해가 전역적으로 존재함이 증명되어 있다. 이러한 수학적 난제에도 불구하고, 나비에-스토크스 방정식은 공학적 유동의 해석에서 극히 정확한 결과를 제공하는 것으로 확립되어 있다.

10. 방정식의 비선형성과 복잡성

나비에-스토크스 방정식의 가장 본질적인 수학적 특성은 대류 항으로부터 오는 비선형성이다. 이 비선형성은 난류의 발생, 다중 해의 존재, 불안정성, 혼돈 거동의 원인이 되며, 유체역학을 수학적으로 풍부하고 동시에 다루기 어려운 분야로 만든다. 선형 점성 항과 비선형 대류 항의 상대적 크기는 레이놀즈 수로 특징지어지며, 이 무차원 수가 유동의 성격을 결정한다.

11. 수치 해법의 개관

나비에-스토크스 방정식의 해석적 해는 극히 제한된 경우에만 얻어지며, 대부분의 실용 문제는 수치 해법으로 해결된다. 유한 체적법, 유한 차분법, 유한 요소법, 스펙트럼 방법, 격자 볼츠만 방법이 대표적 수치 해법이다. 비압축성 해법에서는 속도와 압력의 결합 문제가 특히 중요하며, 투영 방법과 SIMPLE 계열의 반복 알고리즘이 널리 사용된다. 고레이놀즈 수 유동에서는 난류 모델(RANS, LES, DNS)이 계산 비용과 정확도의 균형을 제공한다.

12. 로봇 공학적 응용

나비에-스토크스 방정식은 로봇 공학의 여러 응용에서 기본 지배 방정식으로 사용된다. 드론의 프로펠러와 날개 주위 유동의 정밀 해석, 수중 로봇의 항력과 양력 예측, 소프트 로봇의 유체 구동 시스템의 상세 해석, 생체 모방 수영 로봇의 추진 메커니즘 해석, 미세 로봇의 저레이놀즈 수 유동 해석이 대표적 예이다. 특히 수중 로봇의 설계에서는 나비에-스토크스 방정식 기반의 전산 유체역학 해석이 형상 최적화와 에너지 효율 평가의 표준 도구로 활용된다.

13. 본 절의 의의

본 절은 나비에-스토크스 방정식의 유도와 물리적 의미를 체계적으로 정리한다. 이 방정식은 점성 유체 운동의 가장 일반적인 지배 방정식이며, 이상 유체의 오일러 방정식과 구별되는 점성 효과를 엄밀히 포함한다. 본 절의 내용은 이후 절에서 다룰 단순화 모델, 층류와 난류, 경계층, 항력과 양력의 해석에 필수적인 수학적 기반을 제공한다.

14. 학습 권장사항

독자는 응력 텐서의 구성 관계로부터 나비에-스토크스 방정식을 단계별로 유도해 보고, 각 항의 물리적 의미를 스스로 해석해 볼 것을 권장한다. 또한 단순한 해석적 해(예: Couette 유동, Poiseuille 유동)를 직접 도출하는 실습은 방정식의 구조를 체감하는 유익한 방법이다. 간단한 2차원 비압축성 유동의 수치 해법을 직접 구현하여 속도와 압력의 결합 문제를 관찰하는 것도 권장된다.

15. 참고 문헌

• Navier, C.-L. (1822). Mémoire sur les lois du mouvement des fluides. Mémoires de l’Académie Royale des Sciences de l’Institut de France, 6, 389–440.
• Stokes, G. G. (1845). On the theories of the internal friction of fluids in motion. Transactions of the Cambridge Philosophical Society, 8, 287–319.
• Leray, J. (1934). Sur le mouvement d’un liquide visqueux emplissant l’espace. Acta Mathematica, 63(1), 193–248.
• Batchelor, G. K. (2000). An Introduction to Fluid Dynamics. Cambridge University Press.
• White, F. M. (2011). Fluid Mechanics (7th ed.). McGraw-Hill.
• Kundu, P. K., Cohen, I. M., & Dowling, D. R. (2015). Fluid Mechanics (6th ed.). Academic Press.
• Temam, R. (2001). Navier-Stokes Equations: Theory and Numerical Analysis. AMS Chelsea.
• Panton, R. L. (2013). Incompressible Flow (4th ed.). Wiley.

version: 1.0