복합 영양액 시스템을 위한 ODE 네트워크 모델

본 논문은 폐쇄형 수경재배 시스템에서 영양액의 이온 농도 변화를 정확히 예측하고, 필요한 영양소 투입량을 역산할 수 있는 연속 미분방정식(ODE) 기반 네트워크 모델을 제안한다. 기존의 전도도(EC)·pH만을 측정하고 소량의 농축 영양액을 반복 투입하는 전통적 방법은 이온 흡수 속도 차이로 인해 영양액 내 이온 비율이 서서히 붕괴되고, Na⁺·Cl⁻·SO₄²⁻ 등 흡수율이 낮은 이온이 과다 축적되는 문제를 야기한다. 최근에는 이온 선택 전극(ISE)을 이용해 개별 이온 농도를 측정하고 부족한 이온만 보충하는 방식이 제안되었지만, 다중 이온 간 간섭 현상으로 인해 측정 정확도가 떨어지고, ‘희석‑보충’ 과정을 반복해야 하는 비효율성이 존재한다. 이에 저자들은 영양액을 ‘화학 종(node)’과 ‘반응(edge)’으로 구성된 양방향 네트워크로 모델링하였다. 네트워크 토폴로지는 Yamazaki식 영양액을 기반으로 N, P, K와 미량원소(Fe, Mn, Zn 등)를 포함한 30여 종을 선정하고, 물·수소·수산화 이온, 질산·인산·칼륨 등 주요 이온을 추가하였다. 반응은 크게 ‘증강(촉진)’, ‘억제(제어)’, ‘비상호작용(무관)’ 세 가지 유형으로 구분했으며, Cytoscape를 이용해 시각화하였다. 반응 속도 상수는 실험값이 아닌 문헌에 보고된 평형 상수(K_eq, K_sp, K_a)를 이용해 계산하였다. 평형 상수는 전방·후방 반응 속도 비(k_f/k_b)와 직접 연결되며, 이를 통해 반응이 수렴해야 할 최종 상태를 보장한다. 실제 반응 속도 차이를 반영하기 위해 1~1000 사이의 보정 계수(γ)를 도입했으며, 예를 들어 질산의 용해·이온화는 매우 빠르므로 γ=1000을, Fe‑EDTA의 용해는 느리므로 γ=1을 적용하였다. ODE는 각 종의 농도