윤활 이론

유체역학 전산유체역학 마지막 수정: 2026-07-24 04:49:41

1. 개요[편집]

윤활 이론 Lubrication Theory
정립	Osborne Reynolds (1886)
핵심 가정	막 두께 ≪ 유동 방향 길이
지배 방정식	레이놀즈 윤활 방정식
모태	나비에-스토크스 방정식
대표 무차원수	조머펠트 수, 감소 레이놀즈수
대표 응용	저널 베어링, 씰, 기어, 관절

윤활 이론(lubrication theory)은 두 표면 사이에 낀 아주 얇은 유체 막의 유동을 다루는 근사 이론으로, 막 두께 $h$ 가 유동 방향 길이 스케일 $L$ 보다 압도적으로 작다는 사실( $h/L \ll 1$ )을 이용해 나비에-스토크스 방정식을 압력 하나에 대한 스칼라 편미분방정식으로 축약한 것이다. 그 결과물이 1886년 오즈번 레이놀즈가 유도한 레이놀즈 윤활 방정식이며, 100년 넘게 베어링·씰·기어·관절 설계의 뼈대 노릇을 하고 있다.

핵심은 “3차원 벡터 문제를 2차원 스칼라 문제로 접었다”는 것이다. 막 두께 방향의 속도 분포는 손으로 풀 수 있는 쿠에트 유동+푸아죄유 유동의 중첩으로 못박아 두고, 남은 미지수는 막 평면 위의 압력장뿐이다. 미지수가 하나로 줄어드니 3D 전산유체역학 대비 계산량이 수백~수천 배로 줄어든다. 요즘 워크스테이션이면 저널 베어링 하나쯤 3D로 풀 수도 있지만, 수만 번 돌려야 하는 설계 최적화 루프에서는 여전히 레이놀즈 방정식이 국룰이다.¹

레이놀즈수 60의 2D 비압축성 유동 — 점성 확산이 관성을 억눌러 후류가 얌전히 붙어 있는, 점성 지배 영역의 분위기를 보여준다. 다만 정직하게 밝히면 이건 일반적인 저Re 원기둥 유동 데모이지 얇은 막의 레이놀즈 윤활 방정식을 푸는 게 아니다. 윤활 막은 여기 보이는 것보다 훨씬 극단적인 종횡비(h/L ~ 1/1000)를 가지며, 그 극단이 바로 이 문서가 말하는 근사의 출발점이다.

2. 스케일 해석 — 왜 접히는가[편집]

막 두께 방향을 $y$ , 유동 방향을 $x$ 라 하고 $\varepsilon \equiv h/L \ll 1$ 을 작은 매개변수로 잡자. 연속 방정식에서 두께 방향 속도는 $v \sim \varepsilon u$ 수준으로 억눌린다. 이 스케일을 운동량 방정식에 넣으면 관성항과 점성항의 비가 레이놀즈수 자체가 아니라 감소 레이놀즈수(reduced Reynolds number)로 나온다.

Re^* = \frac{\rho U h}{\mu}\cdot\frac{h}{L} = Re\,\varepsilon

즉 관성이 무시되는 조건은 $Re \ll 1$ 이 아니라 $Re\,\varepsilon \ll 1$ 이다. 이게 실무적으로 어마어마한 차이를 만든다. 고속 저널 베어링의 막 레이놀즈수는 수백에 이르지만, $\varepsilon$ 이 $10^{-3}$ 수준이라 $Re^*$ 는 여전히 1보다 훨씬 작다. 관성을 버려도 되는 것이다.

두께 방향 운동량 방정식은 같은 스케일 논리로 $\partial p/\partial y \approx 0$ 을 준다. 압력은 막 두께 방향으로 일정하다는 이 결론이 윤활 이론의 진짜 심장이다. 남는 것은 앙상한 한 줄이다.

\frac{\partial p}{\partial x} = \mu \frac{\partial^2 u}{\partial y^2}

$p$ 가 $y$ 에 무관하므로 이 식은 $y$ 에 대해 두 번 그냥 적분되고, 벽면 점착 조건을 먹이면 속도 분포가 확정된다 — 움직이는 벽이 끌고 오는 선형 쿠에트 성분과 압력구배가 밀어내는 포물선 푸아죄유 성분의 합이다.

3. 레이놀즈 방정식[편집]

앞의 속도 분포를 두께 방향으로 적분해 유량을 만들고, 그 유량에 질량 보존을 강제하면 압력에 대한 방정식 하나가 떨어진다. 이것이 레이놀즈 윤활 방정식이다.

\frac{\partial}{\partial x}\!\left(\frac{h^3}{12\mu}\frac{\partial p}{\partial x}\right) + \frac{\partial}{\partial z}\!\left(\frac{h^3}{12\mu}\frac{\partial p}{\partial z}\right) = \frac{U}{2}\frac{\partial h}{\partial x} + \frac{\partial h}{\partial t}

좌변은 압력이 유체를 밀어내는 확산성 항이고, 우변은 압력을 만들어내는 두 개의 원천이다. 읽어낼 것이 세 가지다.

$h^3$ 의존성. 계수가 두께의 세제곱이라 막이 절반으로 얇아지면 같은 유량을 흘리는 데 여덟 배의 압력이 필요하다. 베어링이 얇은 막으로 수백 기압을 만들어 하중을 떠받치는 비결이 여기 있다. 동시에 해석 정확도가 막 두께 측정 오차에 잔인하게 민감하다는 뜻이기도 하다.
쐐기(wedge) 항 $\frac{U}{2}\partial h/\partial x$ . 벽이 유체를 끌고 들어가는데 통로가 좁아지면 압력이 솟는다. 뒤집어 말해 완벽히 평행한 두 평판은 아무리 빨리 미끄러져도 하중을 못 버틴다. 미끄럼 베어링에 굳이 경사를 주거나(틸팅 패드), 계단을 파거나(레일리 스텝), 저널을 편심시키는 이유가 전부 이 항을 살리기 위해서다.
스퀴즈(squeeze) 항 $\partial h/\partial t$ . 두 면이 서로 눌러 붙을 때 유체가 빠져나가며 버티는 효과. 물속에서 젖은 유리판 두 장을 떼기 힘든 이유이자, 충격 하중을 받는 베어링이 순간적으로 살아남는 이유다.

방정식은 압력에 대한 가변계수 타원형 방정식이라, 유한차분법·유한요소법·유한체적법 어느 쪽으로 이산화해도 희소행렬 선형계 하나가 나온다. 다중격자법과 궁합이 아주 좋기로 유명하다.

4. 스토크스 유동과 어떻게 다른가[편집]

두 이론 모두 관성을 버리지만, 버리는 근거와 남는 구조가 다르다. 스토크스 유동은 $Re \ll 1$ 이라는 조건 하나로 비선형 관성항을 지워 방정식을 선형화한다. 기하학에 대한 요구는 없다 — 구든 박테리아든 형상이 무엇이든 성립한다. 대신 남는 것은 여전히 벡터 방정식이라 3차원 속도장을 다 풀어야 한다.

윤활 이론은 그 위에 긴 파장(long-wavelength) 가정을 하나 더 얹는다. 형상이 한 방향으로 극단적으로 얇아야 한다는 제약을 받아들이는 대가로, 두께 방향 구조를 해석적으로 소거하고 미지수를 압력 하나로 줄인다. 요약하면 이렇다.

항목	스토크스 유동	윤활 이론
관성 무시 조건	$Re \ll 1$	$Re\,\varepsilon \ll 1$
형상 제약	없음	$h/L \ll 1$ 필수
미지수	속도 벡터 + 압력	압력 스칼라
대표 결과	스토크스 항력	레이놀즈 방정식

따라서 윤활 이론은 스토크스 유동의 부분집합이 아니라, 교집합이 크지만 조건이 서로 어긋나기도 하는 이웃이다. 고속 베어링은 $Re \gg 1$ 이라 스토크스 영역이 아니지만 윤활 이론은 멀쩡히 성립한다. 반대로 점성이 지배하는 미생물 주위 유동은 스토크스지만 얇은 막이 아니라 윤활 이론을 못 쓴다.

5. 확장과 응용[편집]

저널 베어링. 회전축이 베어링 중심에서 편심되면 자동으로 쐐기가 생겨 유막 압력이 축을 띄운다. 이 세계의 주인공은 무차원수 중 조머펠트 수이며, 유막이 만드는 강성·감쇠 8계수는 회전체 동역학 해석의 필수 입력이 된다. 오일 휩(oil whip) 같은 불안정도 여기서 나온다.
탄성유체윤활(EHL). 기어 이빨이나 구름 베어링 접촉부에서는 압력이 1 GPa를 넘어 두 가지가 동시에 깨진다 — 점도가 압력에 지수적으로 치솟고(Barus 관계), 표면이 탄성적으로 눌린다. 그래서 레이놀즈 방정식과 헤르츠 접촉 탄성 방정식을 연성해 풀어야 하며, 결과는 그 유명한 말발굽 압력 스파이크와 출구 목(constriction)이다. 이 유막이 구름 접촉 피로 수명을 좌우한다.²
캐비테이션. 발산 구간에서 레이놀즈 방정식은 태연히 음압을 뱉는데, 실제 오일은 그 전에 찢어진다. 그래서 압력을 0 이하로 못 내려가게 하는 상보성 조건(Half-Sommerfeld, Reynolds, JFO 모델)을 얹는데, 이 순간 문제가 선형에서 선형 상보성 문제로 바뀐다. 캐비테이션 처리를 어떻게 하느냐가 하중 지지력 예측을 통째로 흔든다.
Hele-Shaw와 성형 해석. 같은 얇은 막 논리가 사출성형 해석의 Hele-Shaw 근사, 코팅·박막 도포, 미세유체 채널, 심지어 빙하 흐름과 용암류에도 그대로 쓰인다. 생체 쪽에서는 관절 활액막과 각막 눈물막이 단골 소재다.

6. 언제 무너지는가[편집]

윤활 이론은 얌전히 쓰면 강력하지만, 가정을 벗어나면 조용히 틀린 답을 준다. 대표적인 이탈 조건은 다음과 같다.

막 두께가 표면 거칠기와 비슷해질 때. 매끄러운 연속체 가정이 깨지고 혼합윤활·경계윤활 영역으로 넘어간다. 유막 두께 대 거칠기 비( $\lambda$ 비)가 3 이하면 이미 위험 신호다.
급격한 형상 변화·모서리·큰 편심으로 $h/L \ll 1$ 이 국소적으로 깨질 때. 이 경우 국소 영역만 3D CFD로 풀고 나머지는 레이놀즈 방정식으로 푸는 하이브리드가 실용적이다.
막 레이놀즈수와 종횡비의 곱이 커져 관성·난류가 살아날 때. 대형 터빈 베어링은 실제로 막 안에서 테일러 소용돌이나 난류가 생겨 유효 점도 보정(난류 모델링 계열의 경험 모델)을 얹는다.
열이 무시 못 할 때. 전단이 유체를 데우고 데워진 유체는 점도가 떨어지고 그래서 유막이 얇아진다. 에너지 방정식을 함께 푸는 열유체윤활(THL) 해석이 필요해진다.³

정리하면 윤활 이론은 “가정을 지키는 한 거의 공짜에 가까운 정확도”를 주는, 근사 이론의 모범 사례다. 검증 및 확인 관점에서 중요한 건 결과 숫자보다 내가 지금 그 가정 안에 있는가를 매번 확인하는 습관이다.⁴

7. 관련 문서[편집]

3D CFD로 저널 베어링 한 케이스 푸는 데 수십 분이 걸린다면, 같은 조건의 레이놀즈 방정식은 노트북에서 수십 밀리초에 끝난다. 설계 공간 탐색이 수만 케이스인 순간 이 차이는 “가능/불가능”의 차이가 된다. ↩
EHL 최소 유막 두께 공식(Hamrock–Dowson 등)은 하중 지수가 $-0.073$ 정도로 놀랍도록 둔감하다. 하중을 두 배로 올려도 유막은 겨우 5% 얇아진다는 뜻이라, 처음 보면 계산 실수인 줄 안다. 유막을 진짜로 죽이는 건 하중이 아니라 속도와 점도다. ↩
이쯤 되면 압력·온도·점도·형상이 서로를 물고 도는 완전한 다물리 연성해석이 된다. “얇은 막 하나 푸는 건데요”라고 시작한 프로젝트가 어느새 연성 반복 수렴을 기도하는 현생으로 바뀌는 전형적인 경로. ↩
레이놀즈 방정식은 가정이 깨져도 발산하지 않고 그럴듯한 압력장을 내놓는다는 점에서 특히 위험하다. 발산해서 멈추는 솔버보다, 태연하게 틀린 답을 주는 솔버가 늘 더 무섭다. ↩