스테판 방정식은 압착력, 점도, 디스크 크기, 간극 및 접근 속도 사이의 윤활 이론 관계식입니다. 반지름이 (R)인 두 평행 원형 디스크가 서로 접근할 때, 반간극 규약 (H(t))을 사용한 방정식은 다음과 같습니다.

[ F(t)=\frac{3\pi \mu U(t)R^4}{8H(t)^3}, ]

여기서 (F)는 법선력, (\mu)는 동점도, (U=-dH/dt)는 각 디스크의 접근 속도, (H)는 디스크 사이의 반간격입니다. 교육용 파일럿 플랜트에서는 (F), (H), (U)를 측정하여 학생들이 윤활 이론을 검증하고 힘-변위 데이터로부터 점도를 추정할 수 있습니다.

핵심 결과는 간극에 대한 강한 역세제곱 의존성입니다: 반간극을 줄이면 필요한 압착력이 (H^{-3})에 따라 증가하고, 디스크 반지름을 늘리면 힘이 (R^4)에 따라 증가합니다. 이 방정식은 간극이 좁고, 유동이 느리며, 점성 효과가 관성보다 우세할 때 신뢰성이 높습니다.

스테판 방정식이 설명하는 내용

물리적 실험

실험에서는 얇은 유체층으로 분리된 두 개의 평행 원형 디스크를 사용합니다. 디스크가 서로 가까워지면서 액체가 좁아지는 간극을 통해 방사형 바깥쪽으로 흐르게 됩니다.

이 장치는 윤활, 고분자 가공, 코팅, 성형 및 유변학적 측정과 관련된 단순화된 압착 유동 공정을 나타냅니다.

주요 변수

변수는 다음과 같습니다.

(R): 디스크 반지름
(H(t)): 디스크 사이의 반간격
(U(t)): 각 디스크의 접근 속도
(F(t)): 가해진 법선력
(\mu): 동점도
(\rho): 유체 밀도
(r): 디스크 중심으로부터의 방사 거리

전체 유체층 두께는 (2H)입니다. 이 규약은 중요합니다. 다른 참고 문헌에서는 전체 간극을 (h=2H)로 정의하는 경우가 많아 수치 계수가 다르게 보일 수 있기 때문입니다.

지배 가정

유도 과정에서는 다음을 가정합니다.

크리핑 유동, 관성 무시
좁은 간극, (H_0/R \ll 1)
뉴턴 유체 및 비압축성 유체 거동
디스크 표면에서 미끄럼 없음
평행 디스크 및 축대칭 유동
각 디스크 전체에서 일정한 접근 속도
가장자리 효과 및 유체 누출 무시

관련 레이놀즈 수는 다음과 같습니다.

[ Re=\frac{H_0U_0}{\nu}\ll 1, ]

여기서 (\nu=\mu/\rho)는 동점도입니다.

방정식의 유도 방법

방사형 압착 유동에서 시작하기

디스크가 접근하면 유체는 고체 디스크를 통과할 수 없습니다. 따라서 유체는 방사형 바깥쪽으로 이동합니다.

간극이 디스크 반지름에 비해 작기 때문에 간극을 가로지르는 속도 구배가 방사 방향의 속도 구배보다 우세합니다. 이것이 윤활 이론의 핵심 근사입니다.

운동량 방정식 단순화

크리핑 유동 조건에서는 나비에-스토크스 방정식의 관성항을 무시할 수 있습니다. 방사 방향 운동량 수지는 다음과 같이 단순화됩니다.

[ 0=-\frac{\partial P}{\partial r} +\mu\frac{\partial^2 u_r}{\partial z^2}, ]

여기서 (u_r)는 방사 속도이고 (z)는 유체층을 가로지르는 좌표입니다.

압력은 주로 반지름에 따라 변하고, 방사 속도는 좁은 간극을 가로질러 크게 변합니다.

디스크 경계 조건 적용

미끄럼 없음 조건에 따라 두 디스크 표면에서 방사 속도는 0이어야 합니다.

[ u_r=0 \quad \text{at} \quad z=-H \text{ and } z=H. ]

단순화된 운동량 방정식을 간극을 가로질러 두 번 적분하면 포물선형 방사 속도 분포를 얻습니다. 이 속도 분포에 해당하는 체적 유량은 접근하는 디스크에 의해 밀려나는 유체와 일치해야 합니다.

연속 방정식 사용

두 디스크가 대칭적으로 접근할 때 반간극은 속도 (U)로 감소합니다.

[ \frac{dH}{dt}=-U. ]

변하는 간극은 국소 디스크 면적에 비례하는 속도로 유체를 밀어냅니다. 축대칭 연속 방정식은 방사 유동을 간극 감소율과 연관시킵니다.

연속 방정식과 포물선형 속도 분포를 결합하면 방사 방향 압력 구배는 다음과 같습니다.

[ \frac{\partial P}{\partial r}

-\frac{3\mu U r}{2H^3}. ]

음의 부호는 중심에서 바깥쪽 가장자리로 갈수록 압력이 감소함을 나타냅니다.

압력장 적분

디스크 가장자리의 압력을 기준 압력으로 설정하면, (P(R)=0), (r)에서 (R)까지 적분하여 다음을 얻습니다.

[ P(r)=\frac{3\mu U}{4H^3}\left(R^2-r^2\right). ]

따라서 압력은 중심에서 가장 높고 가장자리로 갈수록 이차적으로 감소합니다.

압력을 적분하여 힘 구하기

전체 법선력은 원형 디스크 면적에 대해 압력을 적분한 값입니다.

[ F=2\pi\int_0^R P(r),r,dr. ]

압력 분포를 대입하면 다음과 같습니다.

[ F

2\pi\int_0^R \frac{3\mu U}{4H^3} \left(R^2-r^2\right)r,dr, ]

따라서

[ \boxed{ F=\frac{3\pi\mu U R^4}{8H^3} }. ]

이는 명시된 반간극 및 대칭 접근 규약에 따른 스테판 방정식입니다.

스케일링이 중요한 이유

간극이 닫힐수록 힘이 급격히 증가

방정식에는 다음 인자가 포함됩니다.

[ F\propto \frac{1}{H^3}. ]

(U), (R), (\mu)가 일정하게 유지된다면 반간극을 절반으로 줄일 때 예측되는 힘은 8배 증가합니다.

이러한 강한 의존성 때문에 간극 측정의 작은 오차도 계산된 점도 또는 예측된 힘에 큰 오차를 일으킬 수 있습니다.

디스크 반지름은 더 강한 영향을 미침

힘은 반지름의 4제곱에 비례합니다.

[ F\propto R^4. ]

더 큰 디스크는 더 많은 유체를 밀어내고 압력이 작용하는 방사 거리를 증가시키므로 훨씬 더 큰 저항을 발생시킵니다.

점도와 접근 속도는 선형적으로 영향

힘은 점도와 접근 속도 모두에 정비례합니다.

[ F\propto \mu U. ]

유체가 뉴턴 유체로 유지되고 크리핑 유동 가정이 유효하다면, 점도가 높은 액체나 더 빠른 압축에는 더 큰 법선력이 필요합니다.

파일럿 플랜트에서 방정식 적용

장치 준비 및 점검

원형 디스크를 설치하고 평행 여부를 확인한 다음 디스크 반지름을 측정합니다. 갇힌 공기는 힘 측정을 왜곡할 수 있으므로 공기가 포획되지 않도록 간극을 주의 깊게 채웁니다.

점도는 온도에 따라 크게 변할 수 있으므로 온도를 모니터링해야 합니다.

운동 제어

알려진 접근 속도 프로파일 (U(t))을 설정하거나 제어된 변위 이력 (H(t))을 부여합니다. 제어기는 명령된 위치에만 의존하지 말고 실제 디스크 위치를 기록해야 합니다.

반간극은 측정된 디스크 간격으로부터 다음과 같이 계산합니다.

[ H(t)=\frac{\text{full disk spacing}}{2}. ]

힘과 변위 기록

압착 단계에서 다음 항목을 기록합니다.

법선력 (F(t))
디스크 위치 또는 간극 (H(t))
접근 속도 (U(t))
유체 온도
시간

힘-변위 곡선은 간극이 좁아질수록 저항이 어떻게 증가하는지 보여주므로 특히 유용합니다.

점도 추정

스테판 방정식을 재배열하면 다음과 같습니다.

[ \boxed{ \mu

\frac{8F H^3}{3\pi U R^4} }. ]

학생들은 측정된 각 데이터 지점에서 순간 겉보기 점도를 계산하거나, 전체 힘-시간 또는 힘-간극 데이터셋에 모델을 적합할 수 있습니다.

점도가 일정한 경우 계산된 값은 유효한 운전 범위에서 대략적으로 일관되어야 합니다.

예측된 스케일링 검증

실험에서는 한 번에 하나의 매개변수만 변경하여 모델을 검증할 수 있습니다.

(F)를 (U)에 대해 도시하여 속도에 대한 선형 의존성을 검증합니다.
점도가 서로 다른 유체를 사용하여 (F)를 (\mu)에 대해 도시합니다.
(F)를 (H^{-3})에 대해 도시하여 간극 의존성을 검증합니다.
서로 다른 반지름의 디스크를 비교하여 (R^4) 의존성을 검증합니다.

로그-로그 플롯을 사용하면 예측된 멱법칙 지수를 실험적으로 더 쉽게 식별할 수 있습니다.

절충 관계 이해

좁은 간극 근사가 실패할 수 있음

이 모델은 다음 조건을 요구합니다.

[ \frac{H_0}{R}\ll 1. ]

디스크 반지름에 비해 간극이 충분히 작지 않으면 방사 방향 효과와 가장자리 효과가 중요해져 1차원 윤활 근사의 정확도가 떨어집니다.

고속 시험에서 관성이 영향을 줄 수 있음

크리핑 유동 조건은 다음을 요구합니다.

[ Re=\frac{H_0U_0}{\nu}\ll 1. ]

접근 속도가 높거나 점도가 낮은 유체를 사용하는 경우 관성력이 측정 가능한 수준이 될 수 있습니다. 파일럿 플랜트 실험에서 관성 간섭을 작게 유지하려면 가해지는 압착력은 다음의 실용적 조건을 만족해야 합니다.

[ F\ll \frac{\mu^2}{\rho} \left(\frac{R}{H_0}\right)^4, ]

이는 파일럿 플랜트 실험에서 관성의 영향을 작게 유지하는 데 도움이 됩니다.

가장자리 누출은 측정된 힘을 감소시킴

이상적인 유도에서는 명확한 원형 경계와 디스크 가장자리에서 기준 압력과 동일한 압력을 가정합니다. 유체가 불균일하게 빠져나가면 압력장과 측정된 힘은 예측값보다 낮아집니다.

평행도가 좋지 않아 비균일한 간극이 형성되는 경우에도 같은 유형의 오차가 발생할 수 있습니다.

뉴턴 유체 거동이 적용되지 않을 수 있음

스테판 방정식은 일정한 동점도를 사용합니다. 고분자 용융물, 고농도 현탁액, 겔 및 많은 구조화 유체는 전단 박화, 전단 농화, 점탄성 또는 항복응력 유체일 수 있습니다.

이러한 재료에서는 측정된 힘이 하나의 물질 상수가 아니라 겉보기 점도 또는 속도 의존적 점도를 나타낼 수 있습니다.

간극 측정은 특히 민감함

힘이 (H^{-3})에 따라 변하기 때문에 영점 설정의 작은 오차나 디스크 간격의 불확실성도 큰 점도 오차를 발생시킬 수 있습니다.

데이터를 수집하기 전에 장치를 변위, 힘, 디스크 평행도 및 초기 간극에 대해 교정해야 합니다.

과도 효과가 중요할 수 있음

이 방정식은 준정상 상태 가정에 따른 점성 압착 응답을 설명합니다. 운동 주기의 시작과 끝에서는 구동기의 가속, 기기의 응답 시간, 갇힌 공기 및 유체 관성이 과도한 편차를 일으킬 수 있습니다.

명령된 운동과 측정된 운동이 일치하는 영역에서 데이터를 해석해야 합니다.

프로젝트에 적용하는 방법

가장 신뢰성 높은 교육용 실험은 제어된 운동, 정확한 간극 측정, 온도 제어 및 힘 데이터를 해석하기 전의 무차원 검증을 결합합니다.

주된 목표가 윤활 이론 검증인 경우: 저점도 또는 특성이 잘 알려진 뉴턴 유체, 작은 간극-반지름비 및 여러 접근 속도를 사용하여 (F\propto U) 및 (F\propto H^{-3})를 검증합니다.
주된 목표가 점도 측정인 경우: 방정식을 재배열하여 (\mu=8FH^3/(3\pi U R^4))를 계산한 다음, 독립적인 점도계 측정값과 결과를 비교합니다.
주된 목표가 스케일 효과 연구인 경우: 디스크 반지름 또는 초기 간극을 체계적으로 변경하고 측정된 힘을 예측된 (R^4) 및 (H^{-3}) 의존성과 비교합니다.
주된 목표가 고분자 또는 비뉴턴 가공 연구인 경우: 결과를 겉보기 점도 측정값으로 취급하고 추정된 점도가 접근 속도와 간극에 따라 어떻게 변하는지 조사합니다.
주된 목표가 신뢰성 있는 데이터 확보인 경우: 스테판 모델을 적합하기 전에 (Re\ll1), (H_0/R\ll1), 디스크 평행도, 유체 구속 및 온도 안정성을 확인합니다.

형상과 운전 조건을 가정 범위 내에서 유지하면 스테판 방정식은 측정된 압착력과 좁은 간극 유동의 점도 구동 물리 사이를 직접 연결해 줍니다.

요약 표:

주제	핵심 정보
방정식	F = 3πμUR⁴/(8H³)
주요 변수	F: 힘; μ: 점도; U: 접근 속도; R: 디스크 반지름; H: 반간극
핵심 스케일링	힘은 μ, U 및 R⁴에 따라 증가하며 H⁻³에 따라 급격히 증가
주요 가정	크리핑 유동, 좁은 간극, 뉴턴 유체, 미끄럼 없음, 평행 디스크, 누출 무시
실험 측정 항목	법선력, 디스크 간격, 접근 속도, 온도 및 시간
점도 계산	μ = 8FH³/(3πUR⁴)
검증 시험	속도, 점도, 간극 또는 디스크 반지름을 변경하고 측정된 멱법칙 경향을 비교
일반적인 오차 원인	갇힌 공기, 불량한 평행도, 간극 측정 오차, 가장자리 누출, 관성 및 온도 변화
교육용 적용 분야	윤활 이론, 유변학, 고분자 가공, 코팅, 성형 및 유체역학 교육

LABPARK와 함께 압착 유동 이론 적용

교과서의 스테판 방정식을 제어되고 측정 가능한 실험으로 구현하세요. LABPARK는 화학공학, 바이오공정 및 생명공학, 환경 및 수처리 교육을 위한 교육용 및 직업교육용 단위조작 파일럿 플랜트를 제공합니다.

당사의 파일럿 플랜트는 대학교, 연구기관 및 기업이 제어된 운동, 힘 및 변위 측정, 온도 모니터링 및 반복 가능한 운전 조건을 통해 유체역학을 교육하고 평가할 수 있도록 지원합니다. 학생과 연구자는 윤활 이론을 검증하고, 점도를 추정하며, 뉴턴 및 비뉴턴 거동을 조사하고, 간극, 디스크 반지름 및 접근 속도가 공정 성능에 미치는 영향을 이해할 수 있습니다.

LABPARK에 문의하기를 통해 교육, 연구 또는 산업 교육 목표에 맞춘 파일럿 플랜트 솔루션을 상담해 보세요.