화학생물공학부: 공정유체역학

안녕하세요. 화학생물공학부에 재학 중인 17기 유재하입니다. 화학공학은 반응기, 분리탑, 배관, 펌프, 열교환기 등 거의 모든 공정 장치 안에서 유체가 끊임없이 흐르는 학문입니다. 그렇기에 유체의 흐름을 정량적으로 다루는 능력은 화학공학도에게 필수적입니다. 오늘 소개해드릴 공정유체역학(Process Fluid Mechanics)은 이러한 유체의 흐름을 몇 개의 기본 보존 법칙에서 출발해 "유도하고 푸는" 법을 배우는 과목입니다. 이 과목에서 무엇을 배우는지, 어떻게 공부하면 좋을지, 그리고 진로와 어떻게 연결되는지 차례대로 소개해드리겠습니다.

---

1. 과목에서 배울 수 있는 내용

1.1 과목의 전반적인 개요

공정유체역학은 한마디로 유체역학의 지배방정식을 직접 유도하고, 그것을 풀어 흐름을 예측하는 과목입니다. 많은 경우 Middleman의 An Introduction to Fluid Dynamics 교재를 따라가며, 학부 유체역학을 한 단계 더 엄밀하게 다룹니다.

수업의 큰 흐름은 다음과 같습니다. 먼저 유체정역학과 표면장력에서 출발해, 검사체적(control volume) 단위로 보존량을 따지는 거시적(적분) 수지와 Bernoulli 방정식을 배웁니다. 이어서 점성과 마찰계수 개념을 통해 층류 관유동을 다룬 뒤, 미소 부피에 대해 수지를 세우는 미분 수지로 넘어가 유명한 Navier–Stokes 방정식을 유도합니다. 마지막으로 경계조건과 평행평판·윤활근사 같은 대표 문제, 그리고 시동(start-up)이나 제트 분열 같은 비정상 유동과 안정성 해석으로 마무리됩니다.

이 과목의 핵심 서사는 두 가지 관점의 연결입니다. 전반부는 시스템 전체를 "점 하나"로 보고 총량만 따지는 (lumped parameter) 관점이고, 후반부는 공간의 각 지점에서 속도장·압력장을 다루는 미소 요소(differential) 관점입니다. 같은 보존 법칙이 어떻게 적분형과 미분형으로 표현되는지를 이해하는 것이 이 수업의 묘미입니다. 또한 여기서 익히는 수학적 구조는 열전달·물질전달과 완전히 동일하기 때문에, 이 과목을 탄탄히 다지면 전달현상(transport phenomena) 과목 전체가 한결 수월해집니다.

1.2 키워드 별 개념 설명

수식적인 부분을 모두 다루기에는 너무 길기에, 핵심 개념과 흐름 위주로 설명하겠습니다. (당연히 실제 수업에서는 각 유도 과정과 수식이 시험에 직접적으로 중요합니다.)

(1) 유체의 정의와 점성, Newtonian 구성방정식

유체는 아무리 작은 전단응력에도 연속적으로 변형하는 물질로 정의됩니다. 고체와 달리 정적 평형 상태에서 전단응력을 버티지 못하고, 오직 수직응력(압력)만 지탱할 수 있습니다. 이때 변형에 저항하는 정도를 나타내는 물성이 점성(viscosity)으로, 물리적으로는 운동량이 확산되는 척도(momentum diffusion)입니다. 전단응력 $\tau$가 변형률(전단속도)에 정비례하는 가장 단순한 유체를 Newtonian fluid라 하며, 그 구성방정식(constitutive equation)은 다음과 같습니다.

$$\tau = \mu \frac{du}{dr}$$

여기서 비례상수 $\mu$가 바로 점도입니다. 점도는 온도에 의존하는데, 액체는 온도가 오를수록 점도가 낮아지고($\ln\mu \approx a + b/T$ 꼴), 기체는 반대로 점도가 높아지는($\mu \propto T^{n}$) 경향을 보입니다. 한 가지 더 기억할 개념은 벽에서 유체의 속도가 벽의 속도와 같다는 점착 조건(no-slip condition)입니다. 19세기에 실험으로 확립된 이 조건은 이후 모든 유동 문제를 푸는 출발점이 됩니다.

(2) 거시적 수지와 Bernoulli 방정식

복잡한 흐름도 검사체적(control volume)을 잡고 그 안팎으로 드나드는 보존량의 수지(balance)를 세우면 정량적으로 분석할 수 있습니다. 임의의 보존량 $X$에 대해 일반적인 수지식은 다음과 같이 쓸 수 있습니다.

$$\frac{dX}{dt} = \dot{X}_{in} - \dot{X}_{out} + \dot{X}_{gen}$$

이 틀을 질량에 적용하면 연속(질량보존), 운동량에 적용하면 Newton 제2법칙의 유체판(운동량보존), 에너지에 적용하면 열역학 제1법칙(에너지보존)이 됩니다. 특히 비압축성 유체에서는 밀도가 일정하므로 질량보존이 곧 부피보존이 된다는 점이 중요합니다. 역학적 에너지 수지를 정리하면 화학공학에서 가장 많이 쓰이는 도구 중 하나인 (일반화) Bernoulli 방정식이 나옵니다.

$$\Delta\!\left(\frac{u^2}{2}\right) + g\,\Delta z + \frac{\Delta P}{\rho} + w_s + f = 0$$

여기서 $w_s$는 펌프 등이 한 일, $f$는 마찰로 인한 손실입니다. 마찰과 일이 없는 이상적인 경우($w_s = f = 0$)가 바로 우리가 잘 아는 고전적 Bernoulli 방정식입니다. 이렇게 시스템을 점 하나로 보고 총량만 다루는 접근을 점값 모델(lumped parameter model)이라 하며, 정상상태에서는 대수방정식, 비정상상태에서는 상미분방정식(ODE)으로 귀결됩니다.

(3) 층류 관유동과 마찰계수 (Hagen–Poiseuille 법칙)

완전발달된 층류 관유동에서는 흐름 방향으로 더 이상 속도가 변하지 않으므로, 힘 평형으로부터 전단응력이 압력구배와 정확히 균형을 이룬다는 결론이 나옵니다. 여기에 Newtonian 구성식과 점착 조건을 결합해 풀면, 그 유명한 포물선 속도분포를 얻습니다. 관 반지름을 $a$, 양단 압력차를 $\Delta P = P_1 - P_2$(입구 − 출구)라 하면

$$v_z(r) = \frac{\Delta P}{4\mu L}\,(a^2 - r^2)$$

이고, 이를 단면에 대해 적분하면 유량과 압력차의 관계인 Hagen–Poiseuille 법칙이 됩니다.

$$Q = \frac{\pi a^4}{8\mu L}\,\Delta P$$

유량이 관 반지름의 무려 4제곱에 비례한다는 점이 핵심 직관입니다. 관이 조금만 굵어져도(혹은 가늘어져도) 유량이 급격히 변하므로, 배관 설계에서 결정적인 결과입니다. 이 마찰손실을 무차원화하면 Fanning 마찰계수 $f_F$가 나오며, 층류에서는 다음의 깔끔한 관계가 성립합니다.

$$f_F = \frac{16}{Re}, \qquad Re = \frac{\rho u D}{\mu}$$

레이놀즈 수 $Re$는 관성력과 점성력의 비로, 층류와 난류를 가르는 핵심 무차원수입니다.

(4) 미분 수지와 Navier–Stokes 방정식

이제 전체 시스템이 아니라 미소 부피 $dV = dx\,dy\,dz$에 대해 보존량의 수지를 세우면, 흐름을 공간의 각 지점에서 기술하는 편미분방정식(장 방정식)을 얻습니다. 질량보존으로부터 연속방정식이,

$$\frac{\partial \rho}{\partial t} + \nabla\cdot(\rho\mathbf{v}) = 0 \quad\xrightarrow{\text{비압축성}}\quad \nabla\cdot\mathbf{v} = 0$$

운동량보존으로부터 Cauchy 운동방정식이 나옵니다.

$$\rho\frac{D\mathbf{v}}{Dt} = \nabla\cdot\mathbf{T} + \rho\mathbf{g}$$

여기서 $\mathbf{T}$는 9개 성분을 가진 응력텐서입니다. 미지수(속도 3개 + 응력 9개 = 12개)가 방정식(운동량 3개 + 질량 1개 = 4개)보다 많아 아직 풀 수 없는데, 이 간극을 메우는 것이 바로 Newtonian 구성방정식입니다.

$$\boldsymbol{\tau} = \mu\!\left(\nabla\mathbf{v} + (\nabla\mathbf{v})^{T}\right)$$

이를 Cauchy 방정식에 대입해 응력을 속도구배로 표현하면, 마침내 유체역학의 가장 유명한 방정식인 Navier–Stokes 방정식이 완성됩니다.

$$\rho\frac{D\mathbf{v}}{Dt} = -\nabla P + \mu\nabla^2\mathbf{v} + \rho\mathbf{g}$$

좌변의 물질미분(material derivative) $\dfrac{D}{Dt} = \dfrac{\partial}{\partial t} + \mathbf{v}\cdot\nabla$ 는 "유체 입자를 따라가며 측정한 변화율"을 뜻합니다. 이 유도 과정에서 Reynolds 수송정리, 발산정리 등 벡터·텐서 해석 도구를 배우게 되는데, 처음엔 낯설지만 후반부 전체를 떠받치는 언어이므로 꼭 익숙해져야 합니다.

(5) 경계조건과 대표 문제 (평행평판·윤활근사)

Navier–Stokes 방정식은 타원형(elliptic) 편미분방정식이라 모든 경계에서 경계조건이 필요합니다. 대표적으로 고체 벽에서의 점착 조건, 유체–유체 계면에서의 응력 연속 및 압력 점프(Young–Laplace 식), 출구에서의 완전발달 조건 등이 있습니다. 이 경계조건을 어떻게 합리적으로 세우느냐가 문제 풀이의 절반을 차지합니다.

가장 기본적인 예제인 평행평판 사이 유동에서는 적절한 가정(정상·비압축·완전발달 등)으로 Navier–Stokes를 대폭 단순화해 단일 ODE로 환원합니다. 위판이 움직이면 선형 속도분포(Couette 유동)가, 압력구배가 걸리면 포물선 속도분포(Poiseuille 유동)가 나옵니다.

$$v_x(y) = \frac{\Delta P}{2\mu L}\,(d^2 - y^2), \qquad \Delta P = P_1 - P_2$$

한 걸음 더 나아가, 채널이 거의 평행하고 얇을 때 차수 분석(order-of-magnitude analysis)으로 방정식을 단순화하는 윤활근사(lubrication approximation)를 배웁니다. 이는 베어링의 유막, 코팅 공정, 와이퍼 등 실제 산업 문제에 직접 응용되는 기법입니다. "어떤 가정을 세우면 어떤 항이 사라지는가"를 판단하는 훈련이 이 과목 후반부의 핵심 역량입니다.

(6) 비정상 유동과 선형 안정성 해석

마지막으로 시동(start-up)처럼 시간에 따라 변하는 과도 상태($\partial/\partial t \neq 0$)를 다룹니다. 보통 해를 정상해와 섭동(perturbation)으로 분리하고, 무차원화($\pi$-정리) 후 Sturm–Liouville 방정식을 풀어 Bessel 함수 해와 지수적 감쇠를 얻습니다. 이로부터 점성 유동의 특성시간이 $t \sim \rho R^2/\mu$ 임을 알게 됩니다.

깊게 다루지는 않지만 층류 제트의 안정성(Rayleigh–Plateau 불안정성)으로 강의가 마무리됩니다. 수도꼭지에서 나온 물줄기가 일정 거리 뒤에서 방울로 끊어지는 현상인데, 교란을 정규모드(normal mode)로 가정해 풀면, 교란의 시간적 성장률 $\alpha$와 공간적 파수 $k$ 사이의 분산관계(dispersion relation)가 나옵니다.

$$\alpha^2 = \frac{\sigma k}{\rho R^2}\,(1 - k^2 R^2)\,\frac{I_1(kR)}{I_0(kR)}$$

파장이 $\Lambda > 2\pi R$ 일 때 제트가 불안정해져 방울로 분열하며, 가장 빠르게 자라는 모드(most dangerous mode)가 방울 크기를 결정합니다($\Lambda_{max} \approx 9.1R \approx 4.55D$). 강의의 중요한 내용들을 모두 이해하면 이 예제를 마지막으로 학부에서 중요한 유체 역학 강의를 마무리할 수 있습니다.

2. 선배의 조언

공정유체역학은 특별한 선이수 과목 없이도 들을 수 있지만, 공학수학(미분방정식·벡터미적분)에 익숙하면 후반부가 훨씬 편합니다. 편미분방정식, 벡터·텐서 표기, 발산정리가 자주 등장하므로 미리 겁먹지 말고 부딪쳐 보시길 권합니다.

이 과목의 본질은 암기가 아니라 유도입니다. 보존 법칙에서 출발해 가정을 세우고 불필요한 항을 지워가며 방정식을 푸는 흐름을 이해하고 관유동, 평행평판, 윤활, 제트 같은 대표 문제는 "어떤 가정으로 어떤 항이 0이 되는지"를 반복하는 것이 가장 효과적입니다.

관련 있는 이동현상경시대회도 좋은 경험이 될 수 있습니다. 열 및 물질전달 과목과 공정유체역학에서 다루는 내용과 관련된 전국경시대회 시험으로, 경시대회를 먼저 준비하고 과목을 수강해도 큰 도움이 될 수 있고, 과목을 수강한 후 경시대회를 준비하셔도 좋을 것 같습니다.

3. 진로 선택에 도움되는 점

화학생물공학부에서 공정(분리·반응·열유체) 분야로 진학한다면 유체역학은 가장 직접적인 기초가 됩니다. 그뿐만 아니라 전산유체역학(CFD), 미세유체(microfluidics), 레올로지·복잡유체, 고분자 가공, 코팅·박막 공정 등 다양한 연구 분야의 출발점이기도 합니다.

이 과목에서 익힌 "보존 법칙 → 구성방정식 → 지배방정식"이라는 사고의 틀은 공학과 관련된 다른 연구들에서도 그대로 통용됩니다. 어떤 물리를 가정하고, 어떤 항을 무시하며, 결과를 실험과 비교하는 과정 자체가 모든 정량적 모델링의 공통 언어이기 때문입니다.

관련 분야로 진학을 고민한다면 미세유체, 다상유동, 레올로지, 전산열유체 등을 다루는 연구실을 찾아보시길 권합니다. (구체적인 교수님이나 연구실은 본인의 관심사와 학과 상황에 맞게 직접 확인해 보시기 바랍니다. 공정 과목의 강의를 맡으시는 거의 모든 교수님들의 연구들에서 공정유체역학 과목은 매우 중요합니다.) 산업 측면에서도 정유·석유화학, 반도체 공정, 바이오 공정 등 유체를 다루지 않는 분야가 거의 없으므로, 유체역학에 대한 탄탄한 이해는 어디서든 든든한 자산이 됩니다.

4. 맺음말

유체역학은 화학공학의 가장 기본적이면서도 가장 강력한 언어입니다. 처음 마주하면 수식이 많아 부담스러울 수 있지만, 결국 모든 것이 몇 개의 보존 법칙에서 출발한다는 사실을 기억하면 길을 잃지 않을 수 있습니다. 한 번 제대로 익혀 두면 대학원 연구와 실무 전반에서 두고두고 꺼내 쓰게 될 든든한 기초가 되어 줄 것입니다. 화학생물공학부에 진학하셔 공정유체역학 강의를 필수로 수강하셔야 한다면 흐름을 정량적으로 바라보는 눈을 꼭 길러 보시기를 강력히 추천드립니다!