기본 지배 방정식은 비압축성 연속 방정식과 나비에-스토크스 운동량 방정식입니다. 일정한 밀도(\rho)와 동점도(\mu)를 가진 등온, 비압축성 뉴턴 유체의 경우, 방정식은 다음과 같습니다.

[ \nabla \cdot \mathbf{v}=0 ]

및

[ \rho\left(\frac{\partial \mathbf{v}}{\partial t} +\mathbf{v}\cdot\nabla\mathbf{v}\right)

-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]

이 방정식들은 적절한 장비 형상, 입구, 출구, 벽면 및 초기 조건에 따라 세 가지 속도 성분과 압력장을 결정합니다.

핵심 요약: 파일럿 플랜트 장비에서의 등온, 비압축성 뉴턴 흐름은 주로 질량 보존과 운동량 보존에 의해 모델링됩니다. 에너지, 화학종, 반응 또는 다상 방정식은 단위 조작에 이 두 방정식만으로는 표현할 수 없는 효과가 포함된 경우에만 추가됩니다.

지배 방정식 세트

질량 보존: 연속 방정식

일반적인 질량 보존 방정식은 다음과 같습니다.

[ \frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\mathbf{v})=0. ]

밀도가 일정한 비압축성 유체의 경우,

[ \boxed{\nabla\cdot\mathbf{v}=0} ]

이 조건은 유체 요소가 이동할 때 팽창하거나 수축하지 않음을 나타냅니다. 데카르트 좌표계에서는 다음과 같습니다.

[ \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} =0, ]

여기서 (u), (v), (w)는 각각 (x), (y), (z) 방향의 속도 성분입니다.

운동량 보존: 나비에-스토크스 방정식

선형 운동량 수지는 다음과 같습니다.

[ \boxed{ \rho\frac{D\mathbf{v}}{Dt}

-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g} } ]

여기서 물질 도함수는 다음과 같습니다.

[ \frac{D\mathbf{v}}{Dt}

\frac{\partial\mathbf{v}}{\partial t} + (\mathbf{v}\cdot\nabla)\mathbf{v}. ]

따라서 운동량 방정식은 다음 간의 균형을 나타냅니다.

유체 가속도: (\rho D\mathbf{v}/Dt)
압력력: (-\nabla p)
점성력: (\mu\nabla^2\mathbf{v})
물체력: (\rho\mathbf{g}), 일반적으로 중력

대류항 ((\mathbf{v}\cdot\nabla)\mathbf{v})은 굽은 관, 수축부, 확장부, 혼합기 및 속도가 공간적으로 변하는 기타 영역에서 중요합니다.

데카르트 성분 형식

속도 (\mathbf{v}=(u,v,w))에 대해 운동량 방정식은 다음과 같습니다.

[ \rho\left( \frac{\partial u}{\partial t} +u\frac{\partial u}{\partial x} +v\frac{\partial u}{\partial y} +w\frac{\partial u}{\partial z} \right)

-\frac{\partial p}{\partial x} +\mu\nabla^2u+\rho g_x, ]

[ \rho\left( \frac{\partial v}{\partial t} +u\frac{\partial v}{\partial x} +v\frac{\partial v}{\partial y} +w\frac{\partial v}{\partial z} \right)

-\frac{\partial p}{\partial y} +\mu\nabla^2v+\rho g_y, ]

[ \rho\left( \frac{\partial w}{\partial t} +u\frac{\partial w}{\partial x} +v\frac{\partial w}{\partial y} +w\frac{\partial w}{\partial z} \right)

-\frac{\partial p}{\partial z} +\mu\nabla^2w+\rho g_z. ]

이 세 방정식은 연속 방정식과 결합하여 (u), (v), (w), (p)라는 네 가지 주요 필드에 대한 네 개의 방정식을 제공합니다.

뉴턴 거동이 중요한 이유

뉴턴의 점성 법칙

뉴턴 유체는 변형률에 따라 선형적으로 변하는 점성 응력을 가집니다. 텐서 형식으로 나타내면,

[ \boldsymbol{\tau}

2\mu\mathbf{D}, ]

비압축성 흐름의 경우,

[ \mathbf{D}

\frac{1}{2} \left[ \nabla\mathbf{v} + (\nabla\mathbf{v})^{T} \right] ]

는 변형률 텐서입니다.

일정한 점도와 비압축성 흐름의 경우, 점성 응력의 발산은 다음과 같이 축소됩니다.

[ \nabla\cdot\boldsymbol{\tau}

\mu\nabla^2\mathbf{v}. ]

이 축소 과정이 표준 나비에-스토크스 방정식에서 점성항을 생성하는 원리입니다.

물리적 해석

점성항은 운동량 확산을 설명합니다. 이는 정지된 파이프 벽과 파이프 중앙에서 움직이는 유체 사이의 속도 구배로 인해 발생하는 운동량 전달을 설명합니다.

모델을 닫기 위해 필요한 조건

노슬립(No-slip) 벽면 조건

정지된 고체 벽면에서,

[ \mathbf{v}=0. ]

이 조건은 파이프, 채널, 열교환기 및 충전 장비 내의 벽면 전단 응력, 압력 손실 및 속도 프로파일을 예측하는 데 핵심적입니다.

움직이는 벽의 경우, 벽면에서의 유체 속도는 벽 속도와 같습니다.

입구 조건

입구는 속도 프로파일,

[ \mathbf{v}=\mathbf{v}_{\text{in}}, ]

또는 질량 유량,

[ \dot{m}

\int_A \rho,\mathbf{v}\cdot\mathbf{n},dA. ]

으로 지정할 수 있습니다. 균일한 입구 프로파일은 단순화된 연구에 자주 사용되며, 발달된 프로파일이나 측정된 프로파일은 상세한 파일럿 플랜트 분석에 더 적합합니다.

출구 조건

일반적인 출구 사양에는 알려진 압력,

[ p=p_{\text{out}}, ]

또는 하류 흐름이 충분히 발달했을 때의 적절한 유출 조건이 포함됩니다. 선택된 출구 조건은 예측된 압력장과 재순환에 실질적인 영향을 미칠 수 있습니다.

초기 조건

과도 상태 운전의 경우, 시작 속도장을 정의해야 합니다.

[ \mathbf{v}(\mathbf{x},0)=\mathbf{v}_0(\mathbf{x}). ]

비압축성 압력은 일반적으로 임의의 상수만큼만 결정되므로 압력 기준점도 필요합니다.

파일럿 플랜트 모듈의 일반적인 단순화

정상 흐름

운전 조건이 시간에 따라 크게 변하지 않는 경우,

[ \frac{\partial\mathbf{v}}{\partial t}=0. ]

운동량 방정식은 다음과 같이 됩니다.

[ \rho(\mathbf{v}\cdot\nabla)\mathbf{v}

-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]

이 근사는 정상 상태의 파이프, 채널 및 열교환기 분석에 일반적으로 사용됩니다.

완전히 발달된 내부 흐름

직선 채널이나 파이프에서 완전히 발달된 흐름의 경우, 속도 프로파일은 축 방향으로 변하지 않습니다. 이때 지배 방정식은 압력 구배와 점성 저항 사이의 균형으로 축소될 수 있습니다.

예를 들어, 평행판 사이의 1차원 흐름의 경우,

[ 0=-\frac{dp}{dx}+\mu\frac{d^2u}{dy^2}. ]

압력 구배가 일정할 때 결과적인 속도 프로파일은 포물선 형태가 됩니다.

무시할 수 있는 관성

레이놀즈 수가 충분히 낮으면 관성 효과는 점성 효과에 비해 작습니다. 운동량 방정식은 스토크스 흐름 형태로 축소됩니다.

[ 0=-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]

이 근사는 고점도 유체, 좁은 채널, 낮은 유량의 마이크로 또는 실험실 모듈에 유용합니다.

중력 및 정수압

속도를 무시할 수 있는 경우, 운동량 방정식은 다음과 같이 축소됩니다.

[ \nabla p=\rho\mathbf{g}. ]

이는 정수압 변화를 설명합니다. 흐름 시스템에서 중력은 수직 배관, 충전탑, 탱크 및 자연 대류 관련 구성에서 여전히 중요할 수 있습니다.

방정식이 단위 조작에 적용되는 방식

배관망 및 펌프

이 방정식들은 배관 전체의 속도 분포와 압력 구배를 설명합니다. 통합되거나 축소된 형태는 압력 손실을 추정하고 펌프 요구 사항을 결정하는 데 사용됩니다.

파일럿 플랜트의 경우 전체 필드 방정식을 계산으로 풀 수 있으며, 단순한 파이프 섹션은 이러한 수지에서 유도된 마찰 계수 상관관계를 통해 표현되는 경우가 많습니다.

열교환기

엄격한 등온 유체 흐름 모델링의 경우, 연속 방정식과 운동량 방정식이 속도 분포, 벽면 전단 및 압력 강하를 포함한 유체 역학을 결정합니다.

온도 변화나 열전달 성능이 목표에 포함된다면 에너지 방정식도 추가해야 합니다. 일반적인 스트림 수준의 수지는 다음과 같습니다.

[ Q=\dot{m}C_p\Delta T, ]

교환기 성능은 종종 다음을 통해 관련됩니다.

[ Q=UA\Delta T_{\mathrm{lm}}F. ]

이는 흐름 방정식을 대체하는 것이 아니라 추가적인 열 관계식입니다.

반응기 및 혼합 장비

교반 탱크 및 흐름 반응기에서 나비에-스토크스 방정식은 순환, 재순환 및 혼합 관련 속도장을 예측합니다.

반응, 화학종 전달 및 체류 시간 거동에는 추가적인 보존 방정식이 필요합니다. 유체 흐름 방정식만으로는 전환율이나 농도를 예측할 수 없습니다.

충전층 및 컬럼

충전 장비에서 기공 규모의 나비에-스토크스 방정식은 상세 모델링에 사용될 수 있지만, 공학적 규모의 모델은 일반적으로 이를 다공성 매질 압력 강하 관계로 대체합니다.

기본 원리는 질량 및 운동량 보존을 유지하며, 충전물은 추가적인 저항 항을 통해 표현됩니다.

2유체 또는 다상 모듈

두 개의 섞이지 않는 유체의 경우, 각 상은 계면 조건과 함께 자체 보존 방정식이 필요합니다. 액체-액체 계면에서 일반적인 조건은 다음과 같습니다.

[ \mathbf{v}_1=\mathbf{v}_2 ]

(접선 속도 연속성) 및

[ \boldsymbol{\tau}_1\cdot\mathbf{n}

\boldsymbol{\tau}_2\cdot\mathbf{n} ]

(점성 응력 연속성, 계면 효과가 적절히 처리될 때).

따라서 상 경계, 상간 운동량 전달 또는 상 변화가 중요할 때는 단일 유체 비압축성 방정식만으로는 불충분합니다.

트레이드오프 이해

등온이라고 해서 항상 에너지가 무관한 것은 아님

등온 모델은 온도가 고정되어 있거나 열 효과가 흐름 솔루션에 영향을 미치지 않는다고 가정합니다. 이는 온도 분포, 열전달 계수 또는 열 부하를 제공하지 않습니다.

열교환기, 가열 반응기 또는 점도 변화가 심한 시스템의 경우, 초기 흐름 분석에서 등온 근사를 사용하더라도 에너지 방정식이 필요할 수 있습니다.

비압축성은 정의된 유효 범위를 가짐

액체는 밀도 변화가 작기 때문에 종종 비압축성으로 모델링됩니다. 기체도 속도가 음속보다 충분히 낮고 압력에 의한 밀도 변화가 무시할 수 있는 수준일 때는 비압축성으로 취급될 수 있습니다.

고속 기체 흐름, 큰 압력 변화, 강한 온도 변화 또는 압축성 주도 효과가 관련된 시스템에서는 이 근사를 재고해야 합니다.

일정한 점도는 좋지 않은 근사일 수 있음

표준 방정식은 일정한 (\mu)를 가정합니다. 이는 많은 액체에 대해 적절하지만, 온도 의존성이 강한 유체나 비뉴턴 물질에는 그렇지 않습니다.

점도가 공간적으로 변하는 경우, 점성항은 단순히 (\mu\nabla^2\mathbf{v})로 쓰기보다는 응력 발산 형태로 유지하는 것이 좋습니다.

층류 가정이 나비에-스토크스에 필수적인 것은 아님

나비에-스토크스 방정식은 층류든 난류든 뉴턴 흐름에 적용됩니다. 그러나 난류 척도를 직접 해결하는 것은 계산 비용이 많이 들기 때문에 파일럿 규모 장비에는 난류 모델이나 실험적으로 검증된 상관관계가 필요한 경우가 많습니다.

"점성 소산"은 에너지 수지에 속함

점성 응력은 운동량 방정식에 직접적인 영향을 미칩니다. 점성 소산에 의해 발생하는 열은 에너지 방정식에 나타나며 엄격한 등온 공식에서는 생략됩니다.

이러한 구분은 점성 소산을 운동량 방정식의 추가 항으로 취급하는 일반적인 실수를 방지합니다.

목표에 맞는 올바른 선택

올바른 지배 모델은 파일럿 플랜트 분석에서 무엇을 예측해야 하는지에 따라 다릅니다.

주요 초점이 속도 프로파일과 압력 손실인 경우: 비압축성 연속 방정식과 나비에-스토크스 운동량 방정식을 사용하고, 운전 조건에 적합한 밀도와 점도를 사용하십시오.
주요 초점이 열전달 성능인 경우: 에너지 방정식을 추가하고 흐름 솔루션을 사용하여 대류, 압력 강하 및 열전달 거동을 결정하십시오.
주요 초점이 반응기 전환율이나 화학종 분포인 경우: 유체 역학 방정식에 화학종 또는 반응 수지를 추가하십시오.
주요 초점이 액체-액체 또는 기체-액체 운전인 경우: 상별 보존 방정식을 사용하고 계면 속도, 응력 및 물질 전달 조건을 포함하십시오.
주요 초점이 빠른 장비 규모 설계인 경우: 마찰 계수, 다공성 매질 또는 체류 시간 상관관계와 같은 검증된 축소 모델을 사용하되 그 운전 한계를 인식하십시오.

가정이 명확하게 명시되고 경계 조건이 올바르게 지정되면, 이러한 보존 방정식은 파일럿 플랜트의 유체 흐름 거동을 분석하기 위한 신뢰할 수 있는 기반을 제공합니다.

요약 표:

모델 요소	지배 관계식	주요 목적
질량 보존	∇ · v = 0	비압축성 흐름 연속성 보장
운동량 보존	ρ(∂v/∂t + v · ∇v) = −∇p + μ∇²v + ρg	속도, 압력, 점성 및 물체력 효과 예측
뉴턴 점성	τ = 2μD	점성 응력과 변형률 관계 설정
경계 조건	입구, 출구, 벽면 및 초기 조건	수학적 모델 완성
추가 물리학	에너지, 화학종, 반응 또는 다상 방정식	열전달, 반응 또는 다상을 위한 모델 확장

유체 흐름 교육, 공정 개발 또는 스케일업을 위한 신뢰할 수 있는 파일럿 플랜트 장비가 필요하십니까? LABPARK에 문의하여 요구 사항을 논의하십시오. LABPARK는 화학 공학, 바이오 공정 및 생명 공학, 환경 및 수처리 응용 분야를 위한 교육 및 직업 단위 조작 파일럿 플랜트를 제공합니다. 당사의 모듈식 시스템은 대학, 연구소 및 기업이 실용적이고 응용 중심적인 장비를 통해 압력 강하, 속도 프로파일, 혼합, 열전달, 분리, 반응 및 공정 제어를 연구할 수 있도록 돕습니다. 목표 공정, 운전 조건, 교육 또는 연구 목표, 필요한 규모를 알려주시면 당사 팀이 적합한 구성을 추천하고 파일럿 플랜트 계획을 지원해 드립니다.