부력 여과식 정수 시스템의 동적 수리 저항 및 여과 효율에 대한 물리수학적 모델링 연구

저자: 이론 물리 수리 모델링 연구팀

요약: 본 논문은 부력 여과(Buoyancy Filtration) 시스템의 핵심 메커니즘인 여재(Filter Media) 및 오염물질의 부력 효과를 통합한 새로운 동적 물리수학 방정식을 설계하고 분석한다. 기존의 Kozeny-Carman 방정식을 확장하고, 오염물질 입자의 부력-중력 평형 속도와 유속의 비를 나타내는 **부력 영향 매개변수 ($\Pi_B$)**를 도입하여, 여과층 내 수리 저항($R_H$)의 시간 의존적 변화를 예측하는 정량적 모델을 제시한다. 이 모델은 부력이 여과 효율에 미치는 동적 안정화 효과를 물리적으로 엄밀하게 설명한다.

---

1. 서론 및 이론적 배경

전통적인 심층 여과(Deep Bed Filtration) 방식은 중력 방향으로의 흐름을 이용하며, 여재의 침강된 무게와 압축성에 의해 여과층이 안정화된다. 그러나 부력 여과식 정수 시스템, 특히 팽창성 여과층을 갖는 상향류(Upflow) 시스템에서는 여재의 유효 밀도가 유체 밀도보다 낮거나 오염물질 입자의 밀도가 유체 밀도와 유사하여 부력 효과가 수리 역학적 안정성에 결정적인 영향을 미친다. 이러한 시스템의 동적 거동을 예측하고 최적화하기 위해서는 단순한 Darcy 법칙을 넘어선, **부력력($F_B$)**, **항력($F_D$)**, 그리고 **여과층 내 퇴적물(Deposit)**의 상호작용을 포괄하는 전문적인 물리수학적 방정식이 필수적이다.

본 연구는 여과층 내 오염물질 퇴적(Straining and Adhesion)으로 인한 수두 손실 증가를 예측하는 기존 모델의 한계를 극복하고, 부력이 퇴적물의 안정화/재분산에 미치는 영향을 **부력 영향 매개변수 ($\Pi_B$)**를 통해 정량화하는 것을 목표로 한다.

---

2. 부력-안정화 동적 수리 저항 설계 방정식

부력 여과 시스템에서 수리 저항($R_H$)은 단순히 여과층의 초기 특성($R_{H0}$)과 퇴적된 오염물질의 양($\sigma(t)$)뿐만 아니라, 오염물질 입자의 부력 성분과 유체 흐름 간의 동적 평형에 의해 조절된다. 우리는 여과층의 단위 길이당 수리 저항($R_H(t)$)이 시간에 따라 변화하는 동적 프로세스를 다음과 같은 **비선형 상미분 방정식(Nonlinear ODE)**으로 설계한다.

2.1. 동적 수리 저항 방정식 ($R_H(t)$)

$$ R_H (t) = R_{H0} \cdot \left[ 1 + \Phi_{KC} \cdot \left(\frac{\sigma(t)}{\epsilon_0 - \epsilon_{\min}}\right)^{\beta} \cdot \exp\left( - \Pi_B \cdot \frac{v_{\text{crit}} t}{L} \right) \right] $$

이 방정식은 다음과 같은 주요 구성 요소로 분해된다:

• $R_{H0}$: 초기 수리 저항. Kozeny-Carman 관계에서 유도된 초기 여과층 저항: $$ R_{H0} = \frac{150 \mu L (1-\epsilon_0)^2}{\rho_f g d_{m}^2 \epsilon_0^3} $$
• $\sigma(t)$: 시간 $t$에서의 여과층 단위 부피당 오염물질의 비(比) 퇴적량 ($\text{mass}/\text{volume}$).
• $\epsilon_0, \epsilon_{\min}$: 초기 및 최소 달성 가능 공극률 (Porosity).
• $\Phi_{KC}$: 퇴적물 형태와 여재 표면 특성에 의한 Kozeny-Carman 수정 계수 ($\Phi_{KC} > 0$).
• $\beta$: 퇴적물이 수리 저항 증가에 미치는 영향의 비선형 지수 ($\beta \ge 1$).
• $\exp\left( - \Pi_B \cdot \frac{v_{\text{crit}} t}{L} \right)$: 부력에 의한 수리 저항 안정화/감쇄 항.

2.2. 부력 영향 매개변수 ($\Pi_B$) 설계

**부력 영향 매개변수 ($\Pi_B$)**는 여과층 내에서 오염물질 입자가 흐름에 대항하여 가지는 상대적인 부력 안정화 능력을 나타내는 **무차원 수(Dimensionless Number)**로 정의된다. 이 매개변수는 입자 크기와 유체 특성을 통해 부력 여과의 핵심 물리량을 통합한다:

$$ \Pi_B = \frac{|v_{\text{term}, c}|}{v_s} \cdot \frac{L}{d_m} \cdot \left(\frac{\rho_f}{\rho_p}\right)^{\alpha} $$

여기서:

• $|v_{\text{term}, c}|$: 오염물질 입자의 부력-중력 평형 **종단 속력** (Terminal Velocity). 입자가 유체 밀도와 다를 때, Stokes 법칙 또는 Newton 법칙에 의해 결정됨. $$ v_{\text{term}, c} = \frac{g d_c^2 |\rho_c - \rho_f|}{18 \mu} \quad (\text{Stokes Regime}) $$
• $v_s$: 여과층의 표면 유속 (Superficial Velocity).
• $L$: 여과층의 깊이(높이), $d_m$: 여재의 평균 입자 지름.
• $\rho_f$: 유체 밀도, $\rho_p$: 여재 입자의 밀도, $\rho_c$: 오염물질 입자의 밀도.
• $\alpha$: 밀도 비에 대한 경험적 민감도 지수 ($\alpha > 0$).

$\Pi_B$는 퇴적물의 부력으로 인한 재분산이나 여재의 부력에 의한 팽창 안정화 효과가 전체 수리 저항 증가에 미치는 영향을 정량화한다. $\Pi_B$ 값이 클수록 부력에 의한 안정화 효과가 지배적이어서 $R_H$의 증가율이 시간에 따라 크게 감소한다.

---

3. 동적 질량 보존 및 여과 효율 모델

부력 여과 시스템 내에서의 오염물질 농도 변화는 연속 방정식(Continuity Equation)과 침전/부유 및 여과 메커니즘을 포함하는 동적 질량 보존 법칙을 따라야 한다. 여과층 깊이 $z$와 시간 $t$에 따른 유체 내 오염물질 농도 $C(z, t)$에 대한 기본 방정식을 다음과 같이 정의한다.

3.1. 오염물질 동적 질량 보존 방정식

$$ v_s \frac{\partial C}{\partial z} + \epsilon \frac{\partial C}{\partial t} + \frac{\partial \sigma}{\partial t} = 0 $$

여기서 $v_s \frac{\partial C}{\partial z}$는 이류(Advection) 항, $\epsilon \frac{\partial C}{\partial t}$는 공극 부피 내 농도 축적 항, $\frac{\partial \sigma}{\partial t}$는 여재 표면으로의 침착(Deposition) 속도이다. 이 침착 속도는 부력 효과를 통합한 새로운 동적 수집 효율($\eta_B$)을 통해 정의된다.

3.2. 부력-안정화 침착 속도 ($\partial \sigma / \partial t$)

침착 속도는 유체 흐름에 수직한 단위 면적당 여과 효율($\lambda$)에 비례한다. $\lambda$는 부력에 의해 수정된 동적 여과 효율 계수 $\lambda_B$를 포함한다. 특히 상향류 부력 여과 시스템에서 입자의 침강(혹은 부유) 속도 $v_{\text{term}, c}$가 흐름에 역행하는 경우, 수집 효율은 중력 침강 효율 $\eta_{\text{grav}}$의 반대 방향으로 작용하여 재분산 메커니즘을 유발할 수 있다.

$$ \frac{\partial \sigma}{\partial t} = \frac{3 (1-\epsilon)}{2 d_m} \cdot \eta_B \cdot v_s \cdot C $$

여기서 **부력 기반 총 수집 효율 ($\eta_B$)**은 다음과 같이 모델링된다. $\eta_B$는 확산($\eta_{\text{diff}}$), 관성($\eta_{\text{int}}$), 그리고 부력(수정된 $\eta_{\text{buoy}}$)의 기여도를 포함한다:

$$ \eta_B = \eta_{\text{diff}} + \eta_{\text{int}} + \eta_{0} \cdot \left[ 1 - \Gamma \cdot \left(\frac{v_{\text{term}, c}}{v_s / \epsilon}\right)^{\gamma} \right] \cdot f(\sigma) $$

새롭게 설계된 항 $\left[ 1 - \Gamma \cdot \left(\frac{v_{\text{term}, c}}{v_s / \epsilon}\right)^{\gamma} \right]$은 오염물질 입자의 종단 속도와 공극 속도($v_s / \epsilon$)의 비가 증가함에 따라 침착 효율이 감소하는 **동적 재분산(Dynamic Redispersion)** 효과를 모델링한다. 즉, 입자 부력이 흐름의 전단력보다 강할 때 여과층에서 입자가 이탈하는 현상을 물리적으로 표현한다.

• $\eta_{0}$: 초기 단일 수집 효율 (Initial Single Collector Efficiency).
• $\Gamma$: 부력-유속 비에 대한 재분산 민감도 계수 ($\Gamma > 0$).
• $\gamma$: 비선형 재분산 지수 ($\gamma > 0$).
• $f(\sigma)$: 퇴적된 양에 따른 효율 변화 함수 ($f(\sigma) = 1 + c_1 \sigma - c_2 \sigma^2$, 여기서 $c_1, c_2$는 계수).

---

4. 결론 및 향후 연구

본 연구는 부력 여과 시스템의 동적 거동을 예측하는 데 있어, 기존의 정적인 모델을 뛰어넘는 **부력 영향 매개변수 ($\Pi_B$)**와 **동적 재분산 효율 항**을 도입한 전문적인 물리수학 방정식을 성공적으로 설계하였다. 특히, 수리 저항($R_H$) 모델은 퇴적물의 축적 효과와 부력에 의한 안정화 효과를 지수 함수 항을 통해 비선형적으로 결합함으로써 시스템의 안정적 장기 운영 가능성을 수리적으로 분석할 수 있는 기반을 마련하였다.

향후 연구는 제시된 $\Pi_B$ 모델의 예측 능력을 실제 부력 여과 시스템의 실험 데이터와 비교하여 $\Phi_{KC}$, $\beta$, $\Gamma$, $\gamma$와 같은 매개변수들의 물리적 의미를 정량화하고 시스템 최적 설계에 활용하는 데 초점을 맞출 것이다. 또한, 여재의 미세한 움직임(Micro-movement)과 부력에 의한 국소적인 공극률 변화($\epsilon(z, t)$)를 포함하는 텐서 기반의 다차원 모델로 확장하여 모델의 엄밀성을 높이는 것이 목표이다.

---

5. 기호 목록

• $R_H$: 동적 수리 저항 ($\text{unitless}$ or $\text{pressure drop per velocity}$)
• $R_{H0}$: 초기 수리 저항
• $\Phi_{KC}$: Kozeny-Carman 수정 계수 ($\text{unitless}$)
• $\sigma(t)$: 시간 $t$에서의 비 퇴적량 ($\text{kg}/\text{m}^3$)
• $\epsilon_0$: 초기 공극률 ($\text{unitless}$)
• $\Pi_B$: 부력 영향 매개변수 ($\text{unitless}$)
• $v_{\text{crit}}$: 임계 속도 ($\text{m}/\text{s}$)
• $L$: 여과층 높이 ($\text{m}$)
• $|v_{\text{term}, c}|$: 오염물질 종단 속력 ($\text{m}/\text{s}$)
• $v_s$: 표면 유속 ($\text{m}/\text{s}$)
• $\rho_f$: 유체 밀도 ($\text{kg}/\text{m}^3$)
• $\rho_p$: 여재 입자 밀도 ($\text{kg}/\text{m}^3$)
• $\mu$: 동적 점성 계수 ($\text{Pa}\cdot\text{s}$)
• $\eta_B$: 부력 기반 총 수집 효율 ($\text{unitless}$)
• $\Gamma$: 재분산 민감도 계수 ($\text{unitless}$)

부력 역중력 여과 정수기 시스템에 대한 동역학적 물리수학 분석

물리학적 모델링 및 제어 방정식 설계에 관한 연구

저자: 김민준, 이서연

소속: 첨단 정수 시스템 연구소

초록

본 논문은 부력을 이용하여 역중력 방향으로 물을 여과하는 혁신적인 정수기 시스템의 동역학적 거동을 물리수학적으로 분석한다. 전통적인 중력 기반 여과 방식과 달리, 이 시스템은 부력과 유체역학적 저항, 그리고 다공성 매체 내 미립자 상호작용을 통합적으로 고려해야 하는 복잡한 비선형 시스템이다. 본 연구는 나비에-스토크스 방정식, 다르시의 법칙, 그리고 부력을 통합한 새로운 운동 방정식을 제안하고, 이를 통해 시스템의 여과 효율, 압력 변화, 그리고 부유체의 안정성을 예측하는 모델을 개발한다. 제안된 모델은 시스템의 최적 설계 및 제어 전략 수립에 중요한 이론적 기반을 제공한다.

키워드

부력, 역중력 여과, 나비에-스토크스, 다르시의 법칙, 다공성 매체, 비선형 동역학, 정수 시스템

1. 서론

수질 오염 문제 해결을 위한 정수 기술은 인류의 중요한 과제이다. 기존의 정수 시스템은 주로 중력에 의한 하향식 여과 방식에 의존하지만, 이는 여과 속도와 효율성 면에서 한계가 있다. 본 논문에서 제시하는 부력 역중력 여과 시스템은 물의 부력을 활용하여 여과 과정을 상향식으로 진행함으로써 기존 방식의 단점을 극복하고자 한다. 이 시스템의 물리적 거동을 정확하게 이해하고 예측하기 위해서는 단순히 부력만 고려하는 것이 아니라, 유체의 운동, 다공성 매체 내의 압력 손실, 그리고 미세 입자의 포획 및 배출 과정까지 종합적으로 다루는 엄밀한 물리수학적 모델이 필수적이다.

2. 물리적 모델링

부력 역중력 여과 시스템은 크게 세 가지 주요 물리적 원리에 의해 지배된다: 부력($\vec{F}_B$), 유체역학적 항력($\vec{F}_D$), 그리고 필터 매체 내의 압력 손실($\Delta P$). 물이 시스템 하부에서 상부로 이동하면서 부유하는 필터 매체($\rho_p < \rho_w$)를 통과한다. 이 과정은 다음과 같은 동역학적 방정식을 통해 기술될 수 있다.

3. 시스템의 동역학 방정식 설계

부유하는 필터 매체($\mathbb{P}$)의 부피($V_p$)와 밀도($\rho_p$)를 고려하고, 물($\mathbb{W}$)의 밀도($\rho_w$)를 이용하여 부력 $\vec{F}_B$를 정의한다. 미립자($\mathbb{M}$)의 부피 농도($\phi_m$)와 여과 매체 내의 공극률($\epsilon$)을 고려하여 유체의 실제 속도($\vec{u}$)를 계산한다. 전체 시스템의 운동 방정식은 나비에-스토크스 방정식과 다르시의 법칙을 결합하여 다음과 같이 표현된다. $$ \rho_p V_p \frac{d^2\vec{x}}{dt^2} = \vec{F}_B + \vec{F}_D + \vec{F}_G + \vec{F}_{pf} $$

여기서 각 항은 다음과 같이 정의된다:

$$ \vec{F}_B = -V_p \nabla P_w = -\int_S P_w d\vec{S} $$

이는 부력의 일반적인 정의이며, 중력 가속도($\vec{g}$)를 고려하면 다음과 같이 쓸 수 있다:

$$ \vec{F}_B = - \rho_w g V_p \hat{k} $$

항력 $\vec{F}_D$는 필터 매체의 형태와 속도에 따라 변하며, 비선형적인 항으로 나타난다. 일반적으로는 다음과 같은 형태로 모델링할 수 있다:

$$ \vec{F}_D = - \frac{1}{2} C_D A \rho_w |\vec{v}| \vec{v} $$

여기서 $C_D$는 항력 계수, $A$는 필터 매체의 단면적, $\vec{v}$는 필터 매체의 속도이다. $\vec{F}_G$는 중력으로 $\vec{F}_G = \rho_p V_p \vec{g}$로 표현된다.

가장 중요한 항은 다공성 매체 내부의 유체-고체 상호작용력 $\vec{F}_{pf}$이다. 이는 폰-카르만-다르시 방정식을 통해 모델링할 수 있다:

$$ \vec{F}_{pf} = - V_p \left( \frac{\mu}{K_D} + \frac{\rho_w C_F}{\sqrt{K_I}} |\vec{u}| \right) \vec{u} $$

여기서 $\mu$는 물의 동점성 계수, $K_D$와 $K_I$는 각각 다르시와 폰-카르만 투과율 계수, $C_F$는 비선형 항력 계수이다. $\vec{u}$는 다공성 매체 내부의 유체 유효 속도이며, 매크로 유체 속도 $\vec{U}$와 공극률 $\epsilon$ 사이의 관계는 $\vec{u} = \vec{U} / \epsilon$이다.

4. 여과 매체 내 미립자 포획 모델

미립자 포획 과정은 필터 매체의 공극률 $\epsilon$과 투과율 $K_D$의 시간에 따른 변화로 모델링할 수 있다. 미립자($\mathbb{M}$)의 부피 농도 $\phi_m$는 다음과 같은 보존 방정식을 따른다:

$$ \frac{\partial(\epsilon \rho_w)}{\partial t} + \nabla \cdot (\epsilon \rho_w \vec{U}) = 0 $$

그리고 미립자의 질량 보존 방정식은 다음과 같다:

$$ \frac{\partial(\epsilon \phi_m)}{\partial t} + \nabla \cdot (\epsilon \vec{U} \phi_m) = 0 $$

미립자 포획으로 인해 공극률이 감소하는 효과는 Kozeny-Carman 방정식을 기반으로 한 수정된 모델로 표현할 수 있다. 필터 매체 내부에 포획된 미립자 질량 $m_f(t)$는 다음과 같이 정의된다:

$$ m_f(t) = \int_{0}^{t} \dot{m}_{in}(\tau) d\tau - \int_{0}^{t} \dot{m}_{out}(\tau) d\tau $$

여기서 $\dot{m}_{in}$과 $\dot{m}_{out}$은 각각 시스템으로 유입 및 유출되는 미립자 질량 유량이다. 공극률 $\epsilon$의 시간에 따른 변화는 다음과 같은 비선형 미분방정식으로 나타낼 수 있다:

$$ \frac{d\epsilon}{dt} = -k_{dep} (\vec{U}, \phi_m) \cdot (1-\epsilon)^n $$

여기서 $k_{dep}$는 침착 속도 계수이며, 유체의 속도와 미립자 농도에 의존하는 복잡한 함수이다. 이 방정식은 여과 과정에서 매체의 투과성이 점진적으로 감소하는 현상을 모델링한다.

5. 통합 시스템의 제어 방정식

시스템의 안정적이고 효율적인 운전을 위해서는 필터 매체($\mathbb{P}$)의 위치($h$)와 유량($Q$)을 제어하는 것이 중요하다. 이는 시스템에 가해지는 압력($P_{ext}$)을 조절하여 달성할 수 있다. 시스템의 동역학 방정식과 미립자 포획 모델을 결합하면 다음과 같은 통합 제어 방정식을 얻을 수 있다:

$$ \rho_p V_p \frac{d^2h}{dt^2} = (\rho_w - \rho_p) g V_p - \left( \frac{\mu}{K_D(\epsilon(t))} + \frac{\rho_w C_F}{\sqrt{K_I(\epsilon(t))}} \left| \frac{Q}{A_s \epsilon(t)} \right| \right) \frac{Q}{A_s} V_p - \vec{F}_{friction} $$

여기서 $Q$는 시스템의 유량, $A_s$는 시스템의 단면적, $\vec{F}_{friction}$은 시스템 벽면과의 마찰력을 나타낸다. 공극률 $K_D(\epsilon(t))$와 $K_I(\epsilon(t))$는 시간에 따라 변하며, 이는 미립자 포획 과정에 의해 결정된다.

이 방정식은 시스템의 유량 $Q$를 독립 변수로, 필터 매체의 높이 $h$를 종속 변수로 갖는 비선형 2차 미분방정식이다. 이 방정식을 통해 우리는 주어진 유량에 대해 필터 매체의 안정적인 평형 높이를 예측하고, 시스템이 막히기 시작하는 임계점을 계산할 수 있다. 또한, 실시간으로 압력 $P_{ext}$를 조절하여 $Q$를 제어함으로써 시스템의 여과 효율을 최적화하는 동적 제어 전략을 수립할 수 있다.

6. 결론 및 향후 연구

본 논문은 부력 역중력 여과 정수기 시스템의 복잡한 물리적 거동을 설명하는 엄밀한 물리수학적 모델을 제시하였다. 나비에-스토크스 방정식과 다르시의 법칙을 결합한 새로운 운동 방정식을 통해 부력, 항력, 그리고 필터 매체 내의 유체-고체 상호작용을 통합적으로 분석하였다. 특히, 미립자 포획에 따른 필터 매체의 특성 변화를 모델링함으로써 시스템의 장기적인 성능 저하를 예측할 수 있는 기반을 마련했다. 향후 연구는 이 모델을 기반으로 한 수치 시뮬레이션과 실제 실험 데이터를 비교하여 모델의 정확성을 검증하고, 최적 제어 이론을 적용하여 시스템의 자동화 및 효율 극대화 방안을 모색할 것이다. 또한, 다양한 오염 물질에 대한 이 모델의 확장성을 탐구하여 부력 역중력 여과 기술의 응용 범위를 넓힐 계획이다.