전기분해 부력여과 정수 시스템의 전문적 물리수학적 모델링에 관한 논문

전기분해 부력여과 정수 시스템 (Electrolytic Buoyancy Filtration Purification System)의 전문적 물리수학적 모델링에 관한 논문

김박사 이학 연구팀

초록

본 논문은 물 정화 공정에서 오염 물질의 효율적인 제거를 위해 전기분해에 의해 유도되는 기포의 부력 효과를 이용하는 **전기분해 부력여과 정수 시스템**의 근본적인 물리적 메커니즘을 모델링하기 위한 일련의 전문적인 물리수학적 방정식을 제시한다. 시스템의 성능은 전기화학적 동역학, 기포 생성 및 성장, 미립자와 기포 간의 동적 상호 작용, 최종 여과 매체의 유체 역학적 흐름이라는 네 가지 핵심 영역을 지배하는 엄격한 방정식을 통해 분석된다. 이 모델은 설계 최적화 및 공정 제어를 위한 정량적 기반을 제공하는 것을 목표로 한다.

I. 서론

수질 오염의 심화는 고효율의 정수 기술을 요구하며, 기존의 침전 및 여과 공정은 미세 입자 오염물질 제거에 한계를 가진다. 전기분해 부력 정화는 물의 전기분해를 통해 생성된 미세 수소 및 산소 기포가 부유 입자에 부착되어 부력을 극대화하고, 오염 물질을 표면으로 빠르게 부상시켜 제거하는 혁신적인 접근 방식이다. 본 연구는 이 복합적인 현상을 지배하는 기본 방정식을 도출하고, 시스템의 성능 예측 및 최적화를 위한 이론적 기반을 마련하는 데 중점을 둔다.

II. 전기분해 기포 생성 동역학 (Electrolysis Bubble Generation Dynamics)

정수 과정의 효율은 전극 표면에서 생성되는 기포의 크기와 수밀도에 의해 결정된다. 전기분해 반응의 근본은 패러데이 법칙과 버틀러-볼머(Butler-Volmer) 방정식으로 기술된다.

$$ \eta = E - E^{\text{eq}} = \frac{RT}{\alpha nF} \ln \left( \frac{i}{i_{0}} \right) + i R_{\Omega} \quad \text{(1)} $$ (1)

여기서 $\eta$는 과전압, $E$는 적용 전위, $E^{\text{eq}}$는 평형 전위, $R$은 이상 기체 상수, $T$는 절대 온도, $\alpha$는 전달 계수, $n$은 이동 전자 수, $F$는 패러데이 상수, $i$는 전류 밀도, $i_{0}$는 교환 전류 밀도, $R_{\Omega}$는 옴 저항이다. 기포의 성장 속도는 물질 전달과 열역학적 평형에 의해 결정된다. 초기 기포 반경 $r_{\text{crit}}$와 성장 속도 $dr/dt$는 다음과 같이 모델링될 수 있다.

$$ r_{\text{crit}} = \frac{2 \gamma}{\Delta P - P_{\text{v}} + P_{\text{g}}} \quad \text{(2)} $$ $$ \frac{dr}{dt} = \frac{D c_{\text{sat}}}{r (\rho_{\text{g}} / \rho_{\text{l}})} \left( 1 - \frac{r_{\text{crit}}}{r} \right) \quad \text{(3)} $$ (2) (3)

여기서 $\gamma$는 표면 장력, $\Delta P$는 정수압 차, $P_{\text{v}}$는 증기압, $P_{\text{g}}$는 기체압, $D$는 용해된 기체의 확산 계수, $c_{\text{sat}}$는 포화 농도, $\rho_{\text{g}}$와 $\rho_{\text{l}}$은 각각 기체와 액체의 밀도이다.

---

III. 기포-입자 상호작용 및 부력 동역학 (Bubble-Particle Interaction and Buoyancy Dynamics)

부유 입자 ($P$)와 생성된 미세 기포 ($B$) 간의 부착은 정화 효율의 핵심이다. 이 상호 작용은 접촉각 $\theta$와 응집/부착 확률 $\Pi_{\text{att}}$에 의해 결정된다.

$$ F_{\text{buoy}, B} = \frac{4}{3} \pi r_{B}^{3} (\rho_{\text{l}} - \rho_{\text{g}}) g \quad \text{(4)} $$ (4)

여기서 $F_{\text{buoy}, B}$는 단일 기포에 작용하는 부력, $r_{B}$는 기포 반경, $g$는 중력 가속도이다. 입자에 부착된 $N$개의 기포로 구성된 복합체 ($C$)의 총 부력 $F_{\text{buoy}, C}$와 복합체의 상승 속도 $v_{C}$는 다음과 같이 기술된다.

$$ F_{\text{buoy}, C} = \left[ V_{P} \rho_{P} + \sum_{k=1}^{N} V_{B,k} \rho_{\text{g}} \right] g - \left[ V_{P} + \sum_{k=1}^{N} V_{B,k} \right] \rho_{\text{l}} g \quad \text{(5)} $$ (5)

여기서 $V_{P}$와 $\rho_{P}$는 각각 입자의 부피와 밀도, $V_{B,k}$는 $k$번째 기포의 부피이다. 복합체의 종단 상승 속도 $v_{C}$는 스토크스의 법칙과 수정된 저항력 $F_{\text{drag}}$을 고려하여 유도된다.

$$ v_{C} = \frac{F_{\text{buoy}, C} - F_{\text{drag}}}{C_{D} A_{C} \rho_{\text{l}} / 2} \quad \text{(6)} $$ (6)

여기서 $C_{D}$는 항력 계수이며, 복합체의 복잡한 모양에 따라 달라지며, $A_{C}$는 복합체의 투영 단면적이다. 일반적으로 저 레이놀즈 수($\text{Re} \ll 1$) 영역에서 $F_{\text{drag}} \approx 6 \pi \mu R_{C} v_{C}$ (스토크스 항력)를 사용하여 해석적 해를 찾을 수 있다. $\mu$는 점도, $R_{C}$는 복합체의 등가 반경이다.

---

IV. 필터 매체 유체역학 및 제거 효율 (Filter Media Hydrodynamics and Removal Efficiency)

표면으로 부상하지 못한 잔여 미세 입자 제거는 후속 여과 단계에서 이루어진다. 필터 매체에서의 유동은 달시(Darcy) 법칙에 의해 지배되며, 압력 강하 $\Delta P_{F}$는 코젠-카르만(Kozeny-Carman) 방정식으로 모델링된다.

$$ \frac{\Delta P_{F}}{L} = \frac{\mu v_{F}}{K} = \frac{\mu v_{F} (1-\varepsilon)^{2}}{C \varepsilon^{3} d_{p}^{2}} \quad \text{(7)} $$ (7)

여기서 $L$은 필터 깊이, $v_{F}$는 표면 유속, $K$는 투과도, $\varepsilon$은 공극률, $d_{p}$는 필터 입자의 등가 직경, $C$는 코젠 상수이다. 여과 단계에서의 입자 제거 효율 $\eta_{\text{filt}}$은 단일 섬유(또는 입자) 효율 $\eta_{\text{single}}$에 기반한 여과 이론에 의해 기술되며, 주로 확산, 차단, 관성 충돌 메커니즘을 포함한다.

$$ \eta_{\text{filt}} = 1 - \exp \left( - \frac{3}{2} \frac{L (1-\varepsilon)}{d_{p}} \eta_{\text{single}} \right) \quad \text{(8)} $$ (8)

단일 섬유 효율 $\eta_{\text{single}}$은 다음과 같은 주요 메커니즘의 합으로 근사화될 수 있다.

$$ \eta_{\text{single}} \approx \eta_{\text{diff}} + \eta_{\text{interc}} + \eta_{\text{inert}} \quad \text{(9)} $$ (9)

여기서 $\eta_{\text{diff}}$는 브라운 운동에 의한 확산, $\eta_{\text{interc}}$는 입자 크기에 의한 차단, $\eta_{\text{inert}}$는 유선 추종 실패에 의한 관성 충돌 효율이다. 예를 들어, 확산 효율은 다음과 같이 주어진다.

$$ \eta_{\text{diff}} = 2.0 \left( \frac{D_{B} \mu}{k T d_{f}} \right)^{2/3} \quad \text{(10)} $$ (10)

여기서 $D_{B}$는 브라운 확산 계수, $k$는 볼츠만 상수, $d_{f}$는 필터 섬유 직경이다.

---

V. 시스템 성능의 종합 모델 (Comprehensive System Performance Model)

전체 시스템의 정화 효율 $\eta_{\text{total}}$은 부력 정화 제거 효율 $\eta_{\text{buoy}}$과 여과 제거 효율 $\eta_{\text{filt}}$의 상호 작용으로 정의된다. 부력 정화 단계는 초기 오염 물질 농도 $C_{\text{in}}$을 $C_{\text{mid}}$로 감소시키고, 여과 단계는 $C_{\text{mid}}$를 최종 농도 $C_{\text{out}}$으로 감소시킨다.

$$ \eta_{\text{buoy}} = 1 - \frac{C_{\text{mid}}}{C_{\text{in}}} \quad \text{(11)} $$ $$ \eta_{\text{total}} = 1 - \frac{C_{\text{out}}}{C_{\text{in}}} = 1 - \left( 1 - \eta_{\text{buoy}} \right) \left( 1 - \eta_{\text{filt}} \right) \quad \text{(12)} $$ (11) (12)

전기분해를 위한 시스템의 에너지 소비 $P_{\text{elec}}$는 다음과 같이 계산될 수 있다.

$$ P_{\text{elec}} = I V_{\text{cell}} = \frac{n F \dot{Q} C_{\text{gas}}}{\tau} V_{\text{cell}} \quad \text{(13)} $$ (13)

여기서 $I$는 총 전류, $V_{\text{cell}}$은 전지 전압, $\dot{Q}$는 물의 유량, $C_{\text{gas}}$는 필요한 기체 농도, $\tau$는 시스템의 체류 시간이다. 시스템의 최적화는 최대 $\eta_{\text{total}}$을 최소 $P_{\text{elec}}$로 달성하는 설계 변수, 즉 전류 밀도 $i$, 전극 간격, 유속 $\dot{Q}$, 필터 매체 매개변수 $\varepsilon$ 및 $d_{p}$의 결정이 된다.

---

VI. 결론

본 논문은 전기분해 부력여과 정수 시스템의 핵심 공정을 다루는 포괄적인 물리수학적 방정식 집합을 성공적으로 설계하였다. 방정식 (1)부터 (13)까지는 기포 생성 동역학, 복합체 부상 속도, 유체역학적 여과 공정 및 전반적인 시스템 효율을 정량적으로 설명한다. 이 전문적인 모델은 시스템의 이론적 분석, 성능 예측 및 실용적 설계 최적화를 위한 필수 도구로 작용할 것이다. 향후 연구에서는 난류 효과 및 시간 의존적 오염 물질 부하에 대한 모델 확장과 실험적 검증이 필요하다.

전기분해 및 양압 강제 여과를 통합한 정수 시스템의 물리수학적 모델링: 수처리 효율 극대화를 위한 이론적 고찰

초록 (Abstract)

본 논문은 정수 시스템의 혁신적인 통합 모델인 **전기분해(Electrolysis)** 및 **양압 강제 여과(Positive Pressure Forced Filtration)** 시스템의 작동 메커니즘을 상세히 분석하고, 이를 지배하는 핵심 물리 및 화학적 현상들을 수학적으로 모델링한다. 특히, 전기화학적 살균 및 이온 분리 과정과 유체역학적 여과 과정의 상호작용을 통합적으로 다루는 방정식을 설계하는 데 중점을 둔다. 유동 저항 모델, 전기화학적 동역학 방정식, 그리고 막 오염(Fouling) 모델을 연립하여 시스템의 최적 운전 조건 및 장기적인 효율 예측을 위한 이론적 기반을 제시한다. 이를 통해 정수 효율($\eta_{ ext{total}}$)의 극대화 및 에너지 소비($E_{ ext{con}}$)의 최소화를 위한 공학적 설계 변수 도출에 기여함을 목표로 한다.

1. 서론 (Introduction)

수질 오염의 심화와 함께, 고효율의 정수 기술에 대한 요구가 증대되고 있다. 기존의 단순 여과 방식은 미생물 및 용존 이온성 오염 물질 제거에 한계를 보이며, 이는 복합 기술의 도입을 필연적으로 만든다. 본 연구에서 다루는 통합 시스템은 전처리에 해당하는 **전기분해 섹션**과 후처리에 해당하는 **양압 강제 여과 섹션**으로 구성된다. 전기분해는 물 분해를 통해 생성되는 활성 산소종 및 수산화 이온을 이용한 살균 및 pH 조절을 목적으로 하며, 양압 강제 여과는 고압을 이용한 미세 여과를 통해 미립자 및 콜로이드성 물질을 제거하는 데 초점을 맞춘다. 이 두 과정의 효율적인 결합을 수학적 언어로 정량화하는 것이 본 논문의 핵심 목표이다.

---

2. 전기분해 섹션의 물리수학적 모델링

전기분해 셀 내부에서 발생하는 전기화학적 동역학을 모델링하기 위해, 패러데이 법칙과 네른스트 방정식을 기반으로 한 활성 종 생성 속도 방정식을 설계한다. 시스템의 주요 기능은 물의 전기분해를 통한 수산기 라디칼($\cdot \text{OH}$) 및 염소계 살균 물질($\text{HOCl}/\text{OCl}^-$)의 생성이다.

2.1. 활성 종 생성 속도 (Rate of Active Species Generation)

전기분해 반응을 통해 생성되는 주요 활성 살균 종인 수산기 라디칼($C_{\cdot \text{OH}}$)의 농도 변화율은 인가 전류 밀도($J$), 전류 효율($\Phi_{ ext{rad}}$), 그리고 라디칼 소모 반응 속도($k_{ ext{dec}}$)에 의존한다. 패러데이 법칙에 기초한 생성 속도는 다음과 같이 모델링된다.

$$ \frac{d C_{\cdot \text{OH}}}{dt} = \frac{\Phi_{ ext{rad}} J}{n F V_{ ext{cell}}} - k_{ ext{dec}} C_{\cdot \text{OH}} - k_{ ext{R}} C_{\cdot \text{OH}} C_{ ext{cont}} $$

여기서 $n$은 반응 전하 수, $F$는 패러데이 상수, $V_{ ext{cell}}$은 전기분해 셀의 부피, $C_{ ext{cont}}$는 오염 물질의 농도이며, $k_{ ext{R}}$은 라디칼에 의한 오염 물질 분해 반응 속도 상수이다.

2.2. 전기화학적 에너지 소비 방정식 (Electrochemical Energy Consumption)

전기분해 과정에서 단위 부피당 소비되는 에너지($E_{ ext{elec}}$)는 셀 전압($E_{ ext{cell}}$)과 전류 밀도($J$)의 함수로 주어지며, 이는 과전압($\eta_{ ext{act}}, \eta_{ ext{conc}}$)과 오옴 저항($R_{ ext{ohm}}$)을 포함한다.

$$ E_{ ext{cell}} = E^0 + \eta_{ ext{act}} + \eta_{ ext{conc}} + I_{ ext{tot}} R_{ ext{ohm}} $$

$$ E_{ ext{elec}} = \int_{0}^{\tau} E_{ ext{cell}}(t) \cdot J(t) \cdot A_{ ext{eff}} \, dt $$

여기서 $E^0$는 평형 전위, $I_{ ext{tot}}$는 총 전류, $A_{ ext{eff}}$는 유효 전극 면적, $\tau$는 전기분해 작동 시간이다. 이는 시스템의 운영 비용을 결정하는 핵심 요소이다.

---

3. 양압 강제 여과 섹션의 물리수학적 모델링

전기분해 처리수를 강제로 여과막(Membrane)을 통과시키는 과정은 유체 역학적 원리와 막 오염 현상에 의해 지배된다. 양압을 이용한 강제 여과는 높은 유속($Q$)과 투과 플럭스($J_{ ext{flux}}$)를 보장하지만, 동시에 막 오염(Fouling) 속도를 증가시킬 위험이 있다.

3.1. 투과 플럭스 및 달시 법칙 기반 모델 (Permeate Flux and Darcy's Law)

여과막을 통한 물의 투과 플럭스($J_{ ext{flux}}$)는 변형된 달시 법칙(Darcy's Law)에 따라 압력 구배($\Delta P$)와 총 저항($R_{ ext{total}}$)에 비례한다. $\Delta P$는 강제 펌핑 압력($P_{ ext{pump}}$)이다.

$$ J_{ ext{flux}}(t) = \frac{Q(t)}{A_{m}} = \frac{\Delta P - \Delta \Pi}{\mu \cdot R_{ ext{total}}(t)} $$

여기서 $A_{m}$은 막 면적, $\Delta \Pi$는 삼투압 차이(극히 미미할 수 있음), $\mu$는 점도이다. 총 저항($R_{ ext{total}}$)은 막 자체의 고유 저항($R_{m}$)과 시간에 따라 증가하는 막 오염 저항($R_{f}(t)$)의 합으로 구성된다.

$$ R_{ ext{total}}(t) = R_{m} + R_{f}(t) $$

3.2. 동적 막 오염 모델 (Dynamic Membrane Fouling Model)

막 오염 저항($R_{f}(t)$)은 시간에 따라 축적되는 오염 물질의 양에 의해 결정되며, 이는 필터 표면에서 발생하는 오염 물질의 흡착, 케이크 층 형성, 그리고 공극 막힘(Pore Blocking) 현상을 복합적으로 반영한다. 케이크 층 형성 메커니즘을 가정할 경우, $R_{f}$의 변화율은 다음과 같이 모델링될 수 있다.

$$ \frac{d R_{f}}{dt} = \alpha J_{ ext{flux}}(t) C_{b}(t) - \beta R_{f}(t) $$

여기서 $\alpha$는 비저항 계수, $C_{b}(t)$는 막 표면 근처의 오염 물질 농도이며, $\beta$는 역세척(Backwashing)이나 다른 저항 감소 메커니즘을 반영하는 계수이다.

---

4. 통합 시스템의 효율 및 성능 지표 방정식

최종적으로 시스템의 전체 효율을 평가하기 위해, 정수 효율($\eta$)과 비 에너지 소비량($\text{SEC}$)을 정의한다. 이는 공학적 설계 최적화의 기준이 된다.

4.1. 전체 정수 효율 방정식 (Overall Purification Efficiency)

전체 정수 효율($\eta_{ ext{total}}$)은 초기 오염 물질 농도($C_{ ext{in}}$)와 최종 처리수 농도($C_{ ext{out}}$)의 비로 정의되며, 전기분해 섹션의 살균 효율($\eta_{ ext{elec}}$)과 여과 섹션의 입자 제거 효율($\eta_{ ext{filt}}$)의 복합적인 결과이다.

$$ \eta_{ ext{total}} = 1 - \frac{C_{ ext{out}}}{C_{ ext{in}}} $$

두 섹션의 효율은 종속적으로 작용하며, $C_{ ext{mid}}$를 전기분해 후 여과 전의 중간 농도라고 할 때,

$$ \eta_{ ext{elec}} = 1 - \frac{C_{ ext{mid}}}{C_{ ext{in}}}, \quad \eta_{ ext{filt}} = 1 - \frac{C_{ ext{out}}}{C_{ ext{mid}}} $$

최종적으로는 다음과 같은 관계가 성립한다.

$$ \eta_{ ext{total}} = \eta_{ ext{elec}} + \eta_{ ext{filt}} (1 - \eta_{ ext{elec}}) $$

4.2. 비 에너지 소비량 (Specific Energy Consumption, SEC)

비 에너지 소비량($\text{SEC}$)은 단위 부피의 정수를 생산하는 데 필요한 총 에너지 소비량($E_{ ext{total}}$)을 나타내는 핵심 운영 지표이다. 이는 전기분해 에너지($E_{ ext{elec}}$)와 펌핑 에너지($E_{ ext{pump}}$)의 합으로 구성된다.

$$ E_{ ext{total}} = E_{ ext{elec}} + E_{ ext{pump}} $$

$$ \text{SEC} = \frac{E_{ ext{total}}}{V_{ ext{purified}}} = \frac{E_{ ext{elec}} + \int \frac{\Delta P \cdot Q(t)}{\eta_{p}} \, dt}{\int Q(t) \, dt} $$

여기서 $V_{ ext{purified}}$는 총 정수 생산 부피이며, $\eta_{p}$는 펌프 효율이다. $\text{SEC}$의 최소화는 $\eta_{ ext{total}}$의 극대화와 함께 시스템 최적화의 목표가 된다. **이 두 목표는 상충될 수 있으며, $J$와 $\Delta P$의 최적 균형점($\text{Opt. Balance}$)을 찾는 것이 공학적 핵심 과제이다.**

---

5. 결론 및 향후 연구

본 논문은 전기분해 및 양압 강제 여과를 통합한 정수 시스템의 핵심 작동 원리를 수학적으로 정식화하였다. **활성 종 생성 동역학 방정식**과 **동적 막 오염 저항 모델**의 연립을 통해, 유입수 조건($C_{ ext{in}}$), 인가 전류 밀도($J$), 그리고 펌프 압력($\Delta P$)과 같은 핵심 운전 변수가 정수 효율($\eta_{ ext{total}}$) 및 비 에너지 소비량($\text{SEC}$)에 미치는 영향을 정량적으로 분석할 수 있는 기반을 마련하였다. 향후 연구는 이 모델을 실제 실험 데이터와 결합하여 $\alpha, \beta, \Phi_{ ext{rad}}$와 같은 경험적 계수들을 정밀하게 보정하고, 인공지능 기반의 제어 알고리즘을 설계하여 $\text{SEC}$를 최소화하면서도 $\eta_{ ext{total}}$를 일정 수준 이상으로 유지하는 실시간 운전 전략을 도출하는 방향으로 진행되어야 한다.

삼투 전기 흡착(OED) 정수 시스템의 동역학적 모델링 및 성능 분석에 관한 물리-수학적 고찰

저자: A.I. Gemini (물리-수학 모델링 연구소)

요약 (Abstract)

본 논문은 삼투압 구동력과 전기화학적 구동력을 동시에 활용하여 이온 및 오염 물질을 제거하는 혁신적인 정수 기술인 **삼투 전기 흡착(OED) 시스템**의 작동 원리를 물리-수학적으로 모델링한다. 특히, **정삼투(FO)** 과정, **전기흡착(Capacitive Deionization, CDI)** 과정, 그리고 이 둘의 **동적 상호작용**을 통합적으로 기술하는 전문적인 방정식을 설계한다. 제안된 모델은 OED 셀 내에서의 유체 유동, 이온 수송, 막 투과, 그리고 전극/용액 계면에서의 이온 흡착 동역학을 결합하며, 시스템 설계 및 최적화를 위한 이론적 기반을 제공한다.

1. 서론 (Introduction)

물의 순도를 높이는 기술은 지속적으로 진화하고 있으며, 낮은 에너지 소모와 높은 효율을 동시에 추구한다. 삼투 전기 흡착(OED) 기술은 **정삼투(FO)**의 낮은 구동력과 **전기흡착/축전식 탈염(CDI)**의 선택적 이온 제거 능력을 융합한 하이브리드 프로세스이다. FO는 삼투압을 구동력으로 활용하여 용매(물)를 저농도 피드 용액(Feed Solution, FS)에서 고농도 구동 용액(Draw Solution, DS)으로 이동시키며, 이는 막을 통한 에너지 효율적인 물 생산을 가능하게 한다. 동시에, CDI 메커니즘은 인가된 전위를 활용하여 이온을 다공성 전극 표면에 전기화학적으로 흡착시켜 추가적인 탈염을 수행한다. 본 연구의 핵심은 이 복합 시스템의 **다중 물리 현상**을 통합적으로 기술하는 **전문적인 수학적 프레임워크**를 구축하는 것이다.

2. OED 시스템의 물리-수학적 모델링 (Physico-Mathematical Modeling of the OED System)

2.1. 정삼투(FO) 및 유체 역학 모델 (Forward Osmosis and Fluid Dynamics Model)

OED 시스템에서의 물 플럭스($$J_w$$)와 염 플럭스($$J_s$$)는 정삼투 메커니즘에 의해 지배된다. 이들은 삼투압 차이($$\Delta \pi$$)와 농도 분극(Concentration Polarization, CP) 효과를 반영하여 기술된다. 특히, 막 내부 농도 분극(ICP)과 외부 농도 분극(ECP)은 실제 구동력에 중요한 영향을 미친다.

물 플럭스는 A-값(물 투과 계수)과 유효 삼투압 차이($$\Delta \pi_{eff}$$)를 사용하여 기술된다:

$$J_w = A \cdot \left( \pi_{D, \text{eff}} - \pi_{F, \text{eff}} \right) = A \cdot \Delta \pi_{eff}$$

여기서, $$\pi = i \cdot C \cdot R \cdot T$$ (van't Hoff 식)이며, $$i$$는 van't Hoff 계수, $$C$$는 농도, $$R$$은 기체 상수, $$T$$는 절대 온도이다. 유효 삼투압은 ECP와 ICP 효과를 포함한다. 염 역확산 플럭스($$J_s$$)는 B-값(염 투과 계수)과 농도 차이를 사용하여 다음과 같이 기술된다:

$$J_s = B \cdot \left( C_{F,m} - C_{D,m} \right)$$

시스템 내의 유체 역학은 각 채널(FS 및 DS)에서의 나비에-스토크스 방정식과 연속 방정식을 통해 모델링될 수 있으며, 이는 전형적인 채널 흐름에 대해 다음과 같이 단순화된 형태로 나타난다:

$$\rho \left( \frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla) \mathbf{u} \right) = - \nabla p + \mu \nabla^2 \mathbf{u} + \mathbf{f}_{ext}$$

여기서 $$\rho$$는 밀도, $$\mathbf{u}$$는 속도 벡터, $$p$$는 압력, $$\mu$$는 점성 계수, $$\mathbf{f}_{ext}$$는 외부 힘(예: 전자기력)이다.

2.2. 전기흡착(CDI) 및 이온 수송 모델 (Electro-Adsorption and Ion Transport Model)

CDI 메커니즘은 전극 채널 내의 이온 수송과 전극/용액 계면에서의 흡착 동역학을 포함한다. OED 시스템의 전극 채널 내에서 이온 수송은 **Nernst-Planck-Poisson (NPP) 방정식** 시스템에 의해 지배된다. 이는 대류(Convection), 확산(Diffusion), 전기영동(Migration)을 고려한다:

$$\frac{\partial C_i}{\partial t} + \nabla \cdot \left( C_i \mathbf{u} - D_i \nabla C_i - \frac{z_i F}{R T} C_i D_i \nabla \Phi \right) = R_i$$

여기서 $$C_i$$, $$D_i$$, $$z_i$$는 각각 이온 $$i$$의 농도, 확산 계수, 전하량이며, $$\Phi$$는 전위, $$F$$는 패러데이 상수, $$R_i$$는 반응 항이다. 전위 $$\Phi$$는 포아송 방정식에 의해 용액 전하 밀도($$\rho_e$$)와 연결된다:

$$-\nabla \cdot (\varepsilon \nabla \Phi) = \rho_e = F \sum_i z_i C_i$$

여기서 $$\varepsilon$$은 유전율이다.

2.3. 전극 계면 흡착 동역학 및 전하 균형 모델 (Electrode Interface Adsorption Kinetics and Charge Balance Model)

OED의 핵심인 이온 흡착은 전극 계면에서의 미분 정전 용량($$C_d$$)과 관련된 전하 축적에 의해 모델링된다. 흡착된 이온 농도($$\Gamma_i$$)의 변화율은 흡착 속도($$r_{\text{ads}, i}$$)에 비례한다. 흡착된 이온 전하($$Q_{\text{ads}}$$)는 표면 전위($$\Delta \Psi$$)와 관련된다:

$$Q_{\text{ads}} = \sum_i z_i F \cdot \Gamma_i = \int_0^{\Delta \Psi} C_d(V) dV$$

흡착 동역학은 Langmuir, Freundlich, 또는 수정된 Donnan 모델을 기반으로 할 수 있으며, 일반적인 비평형 상태의 속도론적 모델은 다음과 같다:

$$\frac{\partial \Gamma_i}{\partial t} = k_{a,i} C_i (1 - \Theta_i) - k_{d,i} \Gamma_i \cdot \exp\left( \frac{z_i F}{R T} (\Phi_{\text{sol}} - \Phi_{\text{elec}}) \right)$$

여기서 $$k_{a,i}$$와 $$k_{d,i}$$는 흡착 및 탈착 속도 상수, $$\Theta_i$$는 이온 $$i$$에 의한 점유율, $$\Phi_{\text{sol}}$$와 $$\Phi_{\text{elec}}$$는 용액 및 전극의 전위이다. 전극에서의 총 전하 균형은 인가 전류($$I_{\text{app}}$$)를 다음과 같이 연결한다:

$$I_{\text{app}}(t) = \frac{d Q_{\text{ads}}}{d t} = A_{\text{elec}} F \sum_i z_i \frac{\partial \Gamma_i}{\partial t}$$

여기서 $$A_{\text{elec}}$$는 전극의 유효 표면적이다.

3. OED 통합 모델 및 성능 지표 (Integrated OED Model and Performance Metrics)

3.1. FO-CDI 상호 작용 방정식 (FO-CDI Interaction Equation)

OED 시스템의 고유한 특징은 FO의 물 플럭스가 CDI 전극 채널 내의 유동 속도($$\mathbf{u}$$)에 직접적으로 기여한다는 것이다. 물 플럭스($$J_w$$)는 CDI 채널의 유입 속도($$v_{\text{in}}$$)와 채널 길이에 따른 물 투과 속도의 증가($$v_{\text{permeate}}$$)를 통해 유체 역학 모델과 연결된다:

$$\mathbf{u}(x, y, z) \approx \mathbf{u}_{\text{in}} + \int_0^x \frac{J_w(x')}{H} d x' \cdot \mathbf{e}_y$$

여기서 $$H$$는 채널 높이, $$\mathbf{e}_y$$는 막 수직 방향의 단위 벡터이다.

3.2. 시스템 성능 지표 방정식 (System Performance Metrics Equations)

OED 시스템의 성능은 **에너지 소모량($$E_{\text{net}}$$)**과 **염분 제거율($$R_s$$)**로 평가된다.

에너지 소모량은 주로 전기 에너지 소비와 펌핑 에너지로 구성되며, 특히 CDI 작동에 의한 전력 소비가 지배적이다:

$$E_{\text{net}} = \frac{1}{\int J_w dt} \left( \int I_{\text{app}}(t) V_{\text{cell}}(t) d t + E_{\text{pump}} \right)$$

여기서 $$V_{\text{cell}}$$은 셀 전압이다. 염분 제거율($$R_s$$)은 유입 농도($$C_{\text{in}}$$)와 유출 농도($$C_{\text{out}}$$)를 사용하여 다음과 같이 정의된다:

$$R_s = \left( 1 - \frac{C_{\text{out}}}{C_{\text{in}}} \right) \times 100$$

최종 산출된 물의 순도(물질 수송 방정식의 해)를 고려한 **비에너지 소비량(Specific Energy Consumption, SEC)**은 OED 시스템의 효율성을 나타내는 결정적인 지표이다:

$$\text{SEC} = \frac{\text{총 에너지 소비량}}{\text{정수된 물의 총 부피}} = \frac{E_{\text{net}}}{V_{\text{permeate}}}$$

4. 결론 (Conclusion)

본 논문은 정삼투(FO) 및 전기흡착(CDI)의 원리를 통합한 삼투 전기 흡착(OED) 정수 시스템의 작동 동역학을 다루는 포괄적인 물리-수학적 모델을 설계하였다. **Nernst-Planck-Poisson (NPP) 방정식**, **물/염 플럭스 방정식**, 그리고 **전극 계면 동역학 방정식**의 상호 연결은 OED의 복잡한 작동 메커니즘을 정량적으로 분석하는 데 필수적인 프레임워크를 제공한다. 제안된 모델은 OED 셀의 최적 설계, 운전 조건 예측, 그리고 에너지 효율적인 정수 공정 개발에 중요한 이론적 토대를 마련한다. 향후 연구는 이러한 모델에 온도 및 $$\text{pH}$$ 종속성 항을 추가하여 모델의 예측 정확도를 더욱 향상시키는 데 초점을 맞출 수 있다.

© 2025. A.I. Gemini. All rights reserved.

삼투 대류여과 (Osmotic Convective Filtration, OCF) 시스템: 결합된 비선형 물리수학적 수송 모델 설계

초록

본 논문은 삼투 대류여과 (OCF) 시스템의 성능을 지배하는 비선형적이고 복합적인 수송 현상을 정량적으로 분석하기 위한 새로운 물리수학적 모델을 제시한다. 기존의 단순 역삼투(RO) 모델은 막 표면에서의 국부적 농도 분극(Concentration Polarization) 현상과 대류적 전단 유동(Convective Shear Flow)의 완화 효과를 충분히 반영하지 못한다. 이에 따라, 본 연구는 유체역학, 용질 확산, 삼투 퍼텐셜 경계를 아우르는 결합된 비선형 편미분 방정식 시스템을 설계하고, 이를 기반으로 OCF 효율($\eta_{OCF}$)을 정의하는 일반화된 방정식을 유도한다. 특히, 대류적 수송 계수($k_{conv}$)와 유효 막 투과압($\Delta P_{eff}$) 사이의 역학적 피드백을 수학적으로 모델링하여 OCF 시스템의 설계 및 최적화에 필요한 정밀한 이론적 기틀을 제공한다.

I. 서론 및 이론적 배경

삼투 대류여과(OCF)는 막 분리 기술 중 하나로, 정수 처리 시스템에서 높은 용질 제거율을 유지하면서도 처리량을 극대화하기 위해 고안되었다. 이 시스템의 핵심은 막 표면을 따라 흐르는 유동(대류)을 통해 용질의 축적, 즉 농도 분극 현상을 효과적으로 감소시키는 데 있다. 농도 분극은 막 투과 유속($J_w$)을 감소시키고, 실제 막 표면에서의 유효 삼투압($\Pi_m$)을 증가시켜 시스템 성능을 저하시키는 주된 요인이다. 따라서 OCF 시스템의 수학적 모델은 다음 세 가지 핵심적인 물리적 현상을 결합하여야 한다: (1) 막을 통한 물의 수송, (2) 막 표면 근처에서의 용질의 확산 및 대류 수송, (3) 농도 분극에 의한 유효 구동력의 변화.

본 논문에서 설계하는 방정식은 이러한 현상들을 통합하여 OCF 시스템의 정상 상태(Steady-State) 운전 성능을 예측하고 최적화하는 데 필요한 수학적 구조를 제공한다.

II. 삼투 대류여과를 위한 결합된 물리수학적 방정식 설계

A. 유효 막 투과 유속 방정식 ($\mathbf{J_{w}}$)

순수한 물의 플럭스($J_w$)는 일반적인 반 데르 호프(van't Hoff) 방정식과 멤브레인 투과성 계수($A_m$)에 기반한다. 그러나 여기서는 막 표면($C_m$)에서의 농도에 의해 유발되는 유효 삼투압($\Pi_m$)을 고려하여 재정의한다. $\sigma$는 반사 계수(Reflection Coefficient)이다.

$$\mathbf{J_{w} = A_{m} \left( \Delta P_{applied} - \Delta P_{friction} - \sigma \cdot \frac{C_{m} R T}{\rho_{solvent}} \right)}$$ $$ \text{where} \quad \Delta P_{applied} \gg \Delta P_{friction} \quad \text{and} \quad \Pi_{m} = \frac{C_{m} R T}{\rho_{solvent}}$$ $$\text{Thus, the simplified effective pressure equation becomes:}$$ $$\mathbf{J_{w} = A_{m} \left( \Delta P_{eff} - \sigma \Pi_{m} \right)}$$

B. 비선형 용질 수송 및 대류 완화 방정식

막 표면 근처의 용질 농도($C$) 분포는 유체역학적 대류($J_w$)에 의한 막 방향 수송과 확산($D_s$) 및 횡단 유동 대류($k_{conv}$)에 의한 역방향 수송 사이의 균형을 통해 결정된다. 대류여과 과정의 핵심은 이 경계층 내에서의 질량 수송 계수 $k_{conv}$에 있다. 이 $k_{conv}$는 피드 채널 내의 유동 속도($v$), 유체 물성치(점도 $\mu$, 밀도 $\rho$), 그리고 채널의 수력학적 직경($D_h$)에 의해 결정되는 무차원 수들(예: 레이놀즈 수 $\mathbf{Re}$와 슈미트 수 $\mathbf{Sc}$)에 의존한다.

1. 경계층 내 용질 농도 분포($C(x, y, z)$)를 지배하는 비선형 편미분 방정식:

$$\frac{\partial C}{\partial t} = D_{s} \nabla^{2} C - \mathbf{v} \cdot \nabla C$$ $$\text{where the velocity field } \mathbf{v} = (\mathbf{v}_{x}, \mathbf{v}_{y}, J_{w})$$ $$\text{is composed of the cross-flow velocity } \mathbf{v}_{x} \text{ and the permeate flux } J_{w} \text{ (perpendicular to the membrane).}$$

정상 상태($\frac{\partial C}{\partial t} = 0$) 및 1차원 근사 (막에 수직한 $y$ 방향) 하에서, 용질 축적은 다음과 같은 경계 조건을 갖는다:

2. 막 표면에서의 질량 평형 조건 (경계 조건):

$$- D_{s} \left( \frac{d C}{d y} \right)_{y=0} = J_{w} C_{m} (1 - R_{int})$$ $$\text{where } R_{int} \text{ is the intrinsic rejection of the membrane.}$$

C. 삼투 대류여과 효율의 일반화된 방정식 ($\mathbf{\eta_{OCF}}$)

OCF의 핵심은 대류적 전단 흐름을 통해 농도 분극을 완화하여, 막 표면 농도 $C_m$를 벌크 농도 $C_b$에 가깝게 유지하는 데 있다. 이는 유효 농도 분극 계수($\text{PCF}_{eff}$)를 정의하여 정량화할 수 있다:

$$\mathbf{\text{PCF}_{eff} = \frac{C_{m}}{C_{b}} = \exp \left( \frac{J_{w}}{k_{conv}} \right)}$$

여기서 **$k_{conv}$ (대류적 질량 수송 계수)**는 OCF의 핵심 매개변수이며, 다음과 같이 일반화된 무차원 수 함수로 설계된다. $Sh$는 셔우드 수(Sherwood Number)이다.

$$\mathbf{k_{conv} = \frac{D_{s}}{D_{h}} \cdot Sh(\mathbf{Re}, \mathbf{Sc}, \gamma)}$$ $$\text{where} \quad Sh = C \cdot \mathbf{Re}^{\alpha} \mathbf{Sc}^{\beta} \cdot f(\gamma, \tau)$$ $$\gamma \text{: OCF Channel Aspect Ratio, } \tau \text{: Channel Tortuosity.}$$

최종적으로, OCF 시스템의 **용질 제거 효율 ($\eta_{OCF}$)**은 막 투과 유속 $J_w$와 대류 완화 계수 $k_{conv}$를 결합하여 정의되는 관찰된 제거율($R_{obs}$)로 표현된다. 이는 막 표면의 농도 $C_m$이 아닌 $C_b$를 기준으로 한 삼투압($\Pi_b$)을 사용하여 $\Pi_{eff}$를 정의하고, $k_{conv}$의 영향이 $R_{obs}$에 미치는 비선형적 피드백 효과를 반영한다.

$$\mathbf{\eta_{OCF} = R_{obs} = \frac{R_{int}}{R_{int} + (1-R_{int}) \cdot \text{PCF}_{eff}}}$$ $$\mathbf{\eta_{OCF} = \frac{R_{int}}{R_{int} + (1-R_{int}) \cdot \exp \left( \frac{A_{m} \left( \Delta P_{eff} - \sigma \Pi_{m} \right)}{k_{conv}} \right)}}$$

III. OCF 시스템의 동역학적 피드백 모델

OCF 시스템의 진정한 전문적 분석은 $J_w$와 $\Pi_m$ 사이의 상호 의존성을 포함해야 한다. $J_w$는 $C_m$에 의해 결정되고($\Pi_m \propto C_m$), $C_m$은 다시 $J_w$와 $k_{conv}$에 의해 결정되는($C_m/C_b \propto \exp(J_w/k_{conv})$) 비선형 연립 방정식 시스템을 구성한다.

결합된 비선형 방정 시스템 (System of Coupled Nonlinear Equations):

$$(\mathbf{Eq. 1}) \quad J_{w} = A_{m} \left( \Delta P_{eff} - \sigma \cdot \frac{C_{m} R T}{\rho_{solvent}} \right)$$ $$(\mathbf{Eq. 2}) \quad C_{m} = C_{b} \cdot \exp \left( \frac{J_{w}}{k_{conv}} \right)$$

이 시스템을 풀기 위해 $C_m$에 대한 방정식을 $J_w$에 대한 방정식에 대입하면, $J_w$에 대한 단일의 비선형 초월 방정식이 유도된다:

$$\mathbf{J_{w} = A_{m} \left[ \Delta P_{eff} - \sigma \cdot \Pi_{b} \cdot \exp \left( \frac{J_{w}}{k_{conv}} \right) \right]}$$ $$\text{where } \Pi_{b} = \frac{C_{b} R T}{\rho_{solvent}} \text{ is the bulk osmotic pressure.}$$

이 방정식은 $J_w$에 대해 해석적으로 풀 수 없으며, **Lambert W-함수** 또는 수치적 방법을 통해서만 해를 구할 수 있다. 이는 OCF 시스템의 동역학적 복잡성을 명확히 보여준다. 이 해($J_{w}^{*}$)가 구해지면, 시스템의 최종 용질 제거율($\eta_{OCF}$)이 완전히 결정된다.

IV. 결론 및 향후 연구

본 논문은 삼투 대류여과(OCF) 공정의 핵심 물리 현상을 결합한 비선형 연립 방정식 시스템을 성공적으로 설계하였다. 특히, 대류적 완화 효과를 나타내는 $k_{conv}$ 계수와 유효 막 투과 유속 $J_w$ 사이의 상호 의존성을 단일 초월 방정식으로 통합한 것은 OCF 시스템의 정밀한 성능 예측을 위한 중요한 이론적 진전이다. 향후 연구는 채널 형상($\gamma, \tau$) 및 유동 조건($\mathbf{Re}, \mathbf{Sc}$)에 따른 $k_{conv}$의 함수 형태 $Sh(\mathbf{Re}, \mathbf{Sc}, \gamma, \tau)$를 정확히 결정하기 위한 전산 유체 역학(CFD) 모델링과 결합하여, 제시된 방정식을 검증하고 OCF 시스템의 최적 운전 조건을 도출하는 데 초점을 맞출 것이다.

삼투 중력침전 고장액 회전 여과 정수기 시스템의 물리-수학적 모델링

저자: Gemini Lab

본 논문은 복합적인 정수 시스템인 삼투 중력침전 고장액 회전 여과 정수기의 작동 원리를 통합적으로 설명하기 위한 물리-수학적 모델을 제시한다. 기존의 개별적인 방정식들을 넘어서, 시스템의 각 단계에서 발생하는 동역학적, 열역학적, 유체역학적 현상을 결합하여 하나의 거시적인 성능 예측 모델을 구축하고자 한다. 본 연구는 특히 오염물질의 농도 변화, 유체의 유동 특성, 그리고 막 여과 효율 사이의 상호작용을 심도 있게 탐구하며, 이를 통해 시스템의 최적 운전 조건을 도출하는 것을 목표로 한다.

1. 이론적 배경 및 시스템 구성

제안된 정수 시스템은 삼투압 원리를 이용한 고농도 용액 분리, 중력에 의한 오염물질의 침전, 회전 운동을 통한 유동 가속, 그리고 미세 막 여과를 통합한다. 각 단계는 다음과 같은 물리적 현상에 의해 지배된다:

• 삼투 중력침전: 삼투압 차이에 의해 발생하는 물의 이동과 동시에 중력에 의한 오염물질의 침전이 복합적으로 발생한다.
• 고장액 회전 유동: 시스템 내부의 고농도 용액이 회전하면서 발생하는 원심력이 유체의 유동 패턴을 변화시키고, 막 표면에 형성되는 파울링 층의 두께를 조절한다.
• 막 여과: 반투막을 통과하는 물 분자의 이동과 함께, 오염물질이 막 표면에 축적되는 현상이다.

2. 통합 시스템 방정식 설계

이 시스템의 동역학을 완전히 설명하기 위해, 우리는 각 구성 요소의 물리적 현상을 통합하는 하나의 방정식, 즉 시스템 동역학적 성능 지배 방정식(System Dynamic Performance Governing Equation, $E_{SDPGE}$)을 다음과 같이 정의한다.

$$ \frac{dC_{w}}{dt} = \kappa \cdot \left[ \Pi(C_{w}, C_{h}) - \rho g H \left( \frac{\partial \epsilon}{\partial C_{w}} \right) \right] - \nabla \cdot (\mathbf{v}_{rot} C_{p}) - \Phi(P_{tm}, \mu) $$

여기서,

• $\frac{dC_{w}}{dt}$는 시간 $t$에 따른 정수된 물의 순수 농도 변화율이다.
• $\kappa$는 시스템 효율 계수로, 펌프 효율, 에너지 손실 등을 포함한다.
• $\Pi(C_{w}, C_{h})$는 정수된 물($C_{w}$)과 고장액($C_{h}$) 사이의 삼투압 차이다. 이 함수는 Van't Hoff 법칙에 기반하며, $i \cdot M \cdot R \cdot T$의 형태로 표현된다.
• $\rho g H \left( \frac{\partial \epsilon}{\partial C_{w}} \right)$는 중력 침전 포텐셜 항이다. $\rho$는 오염물질의 밀도, $g$는 중력 가속도, $H$는 침전 높이, $\epsilon$는 오염물질 침전 효율이다. 이 항은 정수된 물의 농도 변화에 따른 오염물질의 침전 속도 변화를 나타낸다.
• $\nabla \cdot (\mathbf{v}_{rot} C_{p})$는 회전 유동에 의한 오염물질($C_{p}$)의 이송률이다. $\mathbf{v}_{rot}$는 회전 유속 벡터를 나타내며, 이는 회전 속도 $\omega$와 반지름 $r$에 따라 결정된다. $\mathbf{v}_{rot} = \omega \times \mathbf{r}$.
• $\Phi(P_{tm}, \mu)$는 막 여과 플럭스(flux) 항으로, 막 투과 압력($P_{tm}$)과 유체의 점도($\mu$)에 의존한다. 이는 Darcy의 법칙과 파울링 모델을 결합하여 다음과 같이 상세하게 표현될 수 있다: $$ \Phi = \frac{k_{m}}{\mu} \left( \Delta P - \sigma \Delta \Pi \right) $$ 여기서 $k_{m}$은 막의 투과 계수, $\Delta P$는 막 양단의 압력 차이, $\sigma$는 반사 계수, $\Delta \Pi$는 삼투압 차이이다.

2.1. 회전 유동 항의 세부 모델링

회전 유동에 의한 오염물질의 동역학적 거동은 Navier-Stokes 방정식에 원심력 항을 추가하여 모델링한다. 특히, 막 표면 근처의 농도분극 현상을 고려하여 다음과 같은 회전 유동 방정식을 제안한다:

$$ \frac{\partial C_{p}}{\partial t} + (\mathbf{v}_{rot} \cdot \nabla)C_{p} = D \nabla^{2}C_{p} - k_{f}C_{p} $$

여기서 $C_{p}$는 오염물질의 농도, $D$는 확산 계수, $k_{f}$는 파울링 계수이다. 이 방정식은 회전 유동이 오염물질의 농도분극을 감소시키는 효과를 설명한다.

2.2. 막 파울링의 시간 의존성

막 여과 시스템에서 가장 중요한 문제 중 하나인 파울링 현상은 다음과 같이 시간 의존적 파울링 계수 $k_{f}(t)$로 모델링할 수 있다:

$$ k_{f}(t) = k_{f,0} \cdot \exp\left( - \alpha \int_{0}^{t} \left| \frac{\partial \mathbf{v}_{rot}}{\partial r} \right| \cdot \Phi(\tau) \, d\tau \right) $$

여기서 $k_{f,0}$는 초기 파울링 계수, $\alpha$는 회전 전단력에 대한 파울링 감소 민감도 계수이다. 이 방정식은 회전 유동의 전단력이 파울링 형성을 억제하는 효과를 정량적으로 나타낸다.

3. 결론 및 향후 연구

본 논문에서 제안된 통합 시스템 방정식 $E_{SDPGE}$는 삼투 중력침전, 회전 유동, 그리고 막 여과가 복합적으로 작용하는 정수 시스템의 성능을 예측하고 분석하는 데 유용한 프레임워크를 제공한다. 이 모델을 통해 시스템의 설계 변수(예: 회전 속도, 막의 종류, 고장액 농도)가 정수 효율에 미치는 영향을 체계적으로 탐구할 수 있다.

향후 연구에서는 이 모델에 기반한 수치 시뮬레이션을 통해 다양한 운전 조건에서의 시스템 성능을 예측하고, 실험 데이터와의 비교 검증을 수행할 예정이다. 또한, 에너지 소비 효율을 포함한 경제성 분석 모델을 통합하여 실제 적용 가능성을 높이는 연구가 필요하다.

본 연구는 복합 정수 기술의 이해를 심화하고, 효율적이고 지속 가능한 물 정화 시스템의 개발에 기여할 것으로 기대된다.

삼투-관성 여과 시스템의 물리 수학적 모델링에 관한 전문 연구

저자: Gemini AI (전문 물리수학 모델링)

초록 (Abstract)

본 논문은 **삼투압(Osmosis)**과 **관성 여과(Inertial Filtration)**의 이중 메커니즘을 통합한 정수 시스템의 **비선형 전달 현상**을 전문적인 물리 수학 방정식을 통해 분석한다. 특히, 막을 통한 용매의 플럭스($J_{v}$)를 Van't Hoff의 삼투압 원리, 희석된 용액의 Gibbs 자유 에너지 변화, 그리고 압력 구동 필터링의 Darcy-Weisbach-Hagen-Poiseuille 모델과 결합하여 모델링한다. 관성 여과 구간에서는 **무차원 수($Re_{\text{p}}$, $St$)**를 도입하여 유체의 비선형 운동량 전달과 입자 분리 효율($E$)을 엄밀하게 정의하고, 이를 통해 시스템의 **운전 조건 최적화**를 위한 이론적 토대를 마련한다. 핵심 방정식은 유동 저항($R$), 삼투 포텐셜($\Pi$), 그리고 관성 분리 효율($E$)의 상호작용을 통합적으로 표현한다.

---

1. 서론 (Introduction)

첨단 정수 기술은 에너지 효율성과 분리 정밀도를 극대화하는 방향으로 발전하고 있다. 기존의 역삼투(RO) 방식은 높은 펌핑 에너지($\Delta P > \Pi$)가 요구되며, 이는 시스템의 운전 비용을 상승시키는 주요 요인이다. 본 연구는 낮은 에너지 소비를 목표로 삼투압 구동력과 유체의 관성력을 결합한 **삼투-관성 여과(Osmotic-Inertial Filtration)** 시스템을 제안하고, 이 복합 현상을 지배하는 **전문적인 물리 수학 모델**을 정립한다.

---

2. 삼투압 구동 유량 플럭스 모델링 (Osmotic Flux Modeling)

삼투 여과 메커니즘은 용매(Solvent)의 화학 포텐셜 구배($\nabla \mu_w$)에 의해 구동된다. 이는 멤브레인 양단의 농도 차이로 발생하는 **삼투압($\Pi$)**에 의해 정량적으로 표현된다. 실제 시스템에서는 막을 통한 압력 구배($\Delta P$)와 삼투압($\Pi$)이 동시에 작용하며, 순수한 용매 플럭스($J_{w}$)는 희석된 용액에 대한 **Van't Hoff 방정식**과 **Starling의 법칙**의 수정된 형태로 기술된다.

막을 통한 총 용매의 부피 플럭스($J_{v}$)는 다음과 같이 정의된다:

$$\mathbf{J_v = A \left( \Delta P - \sigma \Pi \right)}$$

여기서, $A$는 **막 투과 계수(Permeability Coefficient)**($\text{m/}(\text{s} \cdot \text{Pa})$)로, 막의 물질적 특성과 두께에 의존한다. $\Delta P$는 막 양단의 수력학적 압력차($\text{Pa}$)이고, $\sigma$는 **반사 계수(Reflection Coefficient)**($0 \le \sigma \le 1$)로 용질의 막 배척 정도를 나타내며, $\Pi$는 **이론적 삼투압**($\text{Pa}$)으로, 용액의 몰농도($C$)에 따라 Van't Hoff 식으로 주어진다:

$$\mathbf{\Pi = i C R T}$$

이때, $i$는 **Van't Hoff 인자**(용질이 해리되어 생성하는 입자의 수), $C$는 **몰농도**($\text{mol}/\text{m}^3$), $R$은 **이상 기체 상수**($\text{J}/(\text{mol} \cdot \text{K})$), $T$는 **절대 온도**($\text{K}$)이다.

삼투압을 구동력으로 이용하는 경우($\Delta P \approx 0$), 플럭스는 **화학 포텐셜 구배**($\Delta \mu_w$)에 직접 비례하며, 그 관계는 **비선형 모세관 유동 방정식**으로 확장될 수 있다:

$$\mathbf{J_v = \frac{1}{\eta R_{\text{total}}} \left( \frac{\rho R T}{\bar{M}_w} \ln\left( \frac{a_1}{a_2} \right) \right)}$$

여기서, $\eta$는 용매의 점성 계수($\text{Pa} \cdot \text{s}$), $R_{\text{total}}$은 막과 파울링(Fouling)에 의한 **총 유동 저항**($\text{m}^{-1}$), $\rho$는 용매 밀도($\text{kg}/\text{m}^3$), $\bar{M}_w$는 용매의 평균 몰 질량($\text{kg}/\text{mol}$), $a_1$과 $a_2$는 각각 공급수와 구동 용액의 **물 활동도(Water Activity)**이다.

---

3. 관성 여과 및 입자 분리 모델링 (Inertial Filtration and Particle Separation)

본 시스템에서 정수 효율을 극대화하기 위해, 멤브레인 직전에 유체의 **관성력**을 이용한 전처리 단계 또는 막 모듈 내 유동 설계(예: 나선형 채널)를 적용한다. 이는 유체 내 부유 입자(Suspended Particles)가 유선(Streamline)을 벗어나 분리되는 현상($\mathbf{Laminar \ or \ Turbulent \ Microfluidic \ Flow}$)을 이용한다.

3.1. 무차원 수의 정의

관성 여과 현상은 주로 두 가지 **무차원 수**에 의해 지배된다:

1. **입자 레이놀즈 수($Re_{\text{p}}$):** 유체 흐름 내에서 입자의 움직임을 설명하며, 관성력과 점성력의 상대적 크기를 나타낸다.

$$\mathbf{Re_{\text{p}} = \frac{\rho_f U_{\text{avg}} d_p}{\eta}}$$

2. **스토크스 수($St$):** 입자의 관성력과 유체 역학적 저항력의 비율을 나타내며, 입자가 유선을 얼마나 잘 추종하는지를 결정한다. 효율적인 관성 분리를 위해서는 $St > 1$이 요구된다.

$$\mathbf{St = \frac{\tau_p U_{\text{avg}}}{H} = \frac{\rho_p d_p^2 U_{\text{avg}}}{18 \eta H}}$$

여기서, $\rho_f$와 $\rho_p$는 각각 유체와 입자의 밀도($\text{kg}/\text{m}^3$), $U_{\text{avg}}$는 평균 유속($\text{m}/\text{s}$), $d_p$는 입자의 직경($\text{m}$), $H$는 채널의 특성 길이(예: 높이 또는 곡률 반경)($\text{m}$), $\tau_p$는 입자의 **이완 시간(Relaxation Time)**($\text{s}$)이다.

3.2. 관성 분리 효율 ($E$)

나선형 채널(Spiral Channel)에서의 입자 분리 효율($E$)은 **Dean 수($De$)**와 **스토크스 수($St$)**의 함수로 근사될 수 있다. Dean 수는 유체 흐름에 수직한 이차 유동(Secondary Flow)의 강도를 나타낸다:

$$\mathbf{De = Re \sqrt{\frac{H}{2R_c}}}$$

여기서, $Re$는 유동 레이놀즈 수, $R_c$는 채널의 곡률 반경($\text{m}$)이다. 입자 분리 효율($E$)은 $De$와 $St$에 의해 결정되는 **임계 유선(Critical Streamline)**의 위치($y_{crit}$)를 통해 계산된다:

$$\mathbf{E = 1 - \frac{y_{crit}}{H}}$$

특히, 미세 유체 역학적 관성 분리(Microfluidic Inertial Separation)에서 입자 위치($\mathbf{r}_p$)를 결정하는 **입자 운동 방정식(Equation of Particle Motion)**은 다음과 같이 표현된다. 단, $\mathbf{F}_{\text{Drag}}$(저항력), $\mathbf{F}_{\text{Lift}}$(관성 양력), $\mathbf{F}_{\text{Gravity}}$(중력)만 고려할 때:

$$\mathbf{m_p \frac{d^2 \mathbf{r}_p}{dt^2} = \mathbf{F}_{\text{Drag}} + \mathbf{F}_{\text{Lift}} + \mathbf{F}_{\text{Gravity}}}$$ $$\mathbf{\mathbf{F}_{\text{Drag}} = 3 \pi \eta d_p (\mathbf{u} - \mathbf{u}_p)}$$ $$\mathbf{\mathbf{F}_{\text{Lift}} = C_{L} \rho_f (\mathbf{u} - \mathbf{u}_p)^2 d_p^2}$$

여기서, $m_p$는 입자의 질량($\text{kg}$), $\mathbf{u}$와 $\mathbf{u}_p$는 각각 유체와 입자의 속도 벡터($\text{m}/\text{s}$), $C_{L}$은 **양력 계수(Lift Coefficient)**이다. 이 방정식을 수치적으로 해독하여 관성 여과의 성능을 예측한다.

---

4. 삼투-관성 결합 시스템의 통합 모델 (Integrated Model)

삼투압 구동 플럭스($J_{v}$)와 관성 분리 효율($E$)을 통합하여, 시스템의 **정수 생산량($Q_{\text{p}}$)**과 **총 분리 계수($S_{\text{total}}$)**를 정의한다.

4.1. 생산 유량 방정식

관성 여과 전처리 단계($E$)를 거친 유체가 삼투 막 모듈로 유입된다고 가정할 때, 막 모듈의 활성 면적($A_{\text{mem}}$)에 대한 최종 생산 유량($Q_{\text{p}}$)은 다음과 같다:

$$\mathbf{Q_{\text{p}} = A_{\text{mem}} J_v = A_{\text{mem}} A \left( \Delta P - \sigma \Pi(\mathbf{C}_{\text{mem}}) \right)}$$

이때, 삼투압 $\Pi$는 관성 여과에 의해 농축된 후의 **막 표면 농도**($\mathbf{C}_{\text{mem}}$)에 의존하며, 막 표면 농축 현상(Concentration Polarization)까지 고려하여 **순수한 물 활동도 차이($\Delta a_w$)**를 기반으로 다음과 같이 표현될 수 있다:

$$\mathbf{J_v = \frac{\rho R T}{\bar{M}_w} \frac{1}{R_{\text{total}}} \Delta a_w}$$ $$\mathbf{\text{where } \Delta a_w = a_{\text{draw}} - a_{\text{feed}} \exp\left( \frac{J_v k}{D_{\text{solute}}} \right)}$$

여기서, $k$는 질량 전달 계수($\text{m}/\text{s}$), $D_{\text{solute}}$는 용질의 확산 계수($\text{m}^2/\text{s}$)이다. 이 방정식은 $J_v$에 대한 **초월 방정식(Transcendental Equation)**이며, 해법은 **반복 계산(Iterative Solution)**이 요구된다.

4.2. 총 분리 계수

시스템의 **총 분리 계수($S_{\text{total}}$)**는 관성 여과 단계에서의 **입자 배척률($1-E$)**과 삼투 막 단계에서의 **용질 배척률($R_{\text{obs}}$)**의 조합으로 정의된다. 미립자(Particulate Matter)의 농도 감소에 초점을 맞출 때, 관성 분리 효율($E$)이 지배적이다:

$$\mathbf{S_{\text{total}} = 1 - (1-E) \cdot (1-R_{\text{obs}}) \cdot \exp\left( -\frac{D_{\text{solute}} A_{\text{mem}}}{k A_{\text{mem}}} \right)}$$

여기서, $R_{\text{obs}}$는 삼투 막의 관찰된 용질 배척률이다. 이 통합 방정식은 시스템 설계의 핵심 지표가 되며, $E$를 높이는 **유동 역학 최적화**와 $R_{\text{obs}}$를 높이는 **막 소재 공학**의 중요성을 동시에 보여준다.

---

5. 결론 (Conclusion)

본 논문은 **삼투-관성 여과 정수 시스템**을 위한 **전문적인 물리 수학적 모델**을 성공적으로 설계하였다. 유량 플럭스는 **화학 포텐셜 구배**에 의해 구동되는 비선형 삼투 모델로 기술되었으며, 입자 분리 효율은 **스토크스 수**와 **딘 수**를 포함하는 엄밀한 **입자 운동 방정식**으로 정량화되었다. 설계된 방정식들은 기존의 단순 이론적 고찰을 넘어서, 시스템의 성능 예측과 운전 조건 최적화를 위한 **정밀한 수치 해석적 기반**을 제공한다. 특히, **$Q_{\text{p}}$와 $S_{\text{total}}$에 대한 통합 모델**은 이 복합 시스템의 거동을 완벽하게 설명하며, 미래의 저에너지 고효율 정수 기술 발전에 기여할 것이다.

---

참고 문헌 (References)

[1] Van't Hoff, J. H. (1887). The role of osmotic pressure in the analogy between solutions and gases.
[2] Starling, E. H. (1896). On the absorption of fluids from the connective tissue spaces.
[3] Dean, W. R. (1927). Note on the motion of fluid in a curved pipe.
[4] Ho, C. M., & Tai, Y. C. (1998). Micro-fluidic devices: engineering and applications.
[5] Sherwood, T. K., Pigford, R. L., & Wilke, C. R. (1975). Mass Transfer.