[논문 리뷰] A fluid dynamics multidimensional model of biofilm growth: stability, influence of environment and sensitivity

Fabrizio Clarelli, Cristiana Di Russo외 2명|arXiv (Cornell University)|2014. 08. 09.

Particle Dynamics in Fluid FlowsEngineering참고 문헌 37인용 수 22

한 줄 요약

이 논문은 빛과 온도 의존성을 포함한 다차원 유체역학 모델을 제시하며, 생물막 안정성 분석과 모델 매개변수 민감도 분석을 수행한다. 모델은 30일간의 2차원 및 3차원 생물막 발달을 정확하게 시뮬레이션하여 실험적 두께(0.11–0.12 cm)와 빛 감쇠 패턴과 양호한 일치를 보인다.

ABSTRACT

In this article, we study in details the fluid dynamics system proposed in Clarelli et al (2013) to model the formation of cyanobacteria biofilms. After analyzing the linear stability of the unique non trivial equilibrium of the system, we introduce in the model the influence of light and temperature, which are two important factors for the development of cyanobacteria biofilm. Since the values of the coefficients we use for our simulations are estimated through information found in the literature, some sensitivity and robustness analyses on these parameters are performed. All these elements enable us to control and to validate the model we have already derived and to present some numerical simulations in the 2D and the 3D cases.

연구 동기 및 목표

• 물층 아래의 돌 기질 상에서 생물막 형성에 대한 물리적으로 타당한 유체역학 모델을 개발하기 위해.
• 일차원에서 비자명한 평형점 주변의 선형 안정성 분석을 수행하여 이론적 탄탄성을 확보하기 위해.
• 생물막 성장 과정에서 빛과 온도를 핵심 환경 구동력으로 통합하기 위해.
• 문헌 기반 추정치를 활용하여 모델 매개변수의 민감도 및 탄력성 분석을 수행하기 위해.
• 교정된 모델을 2D 및 3D 수치 시뮬레이션을 통해 실험적 생물막 두께와 빛 분포 패턴과 일치시켜 검증하기 위해.

제안 방법

• 질량 및 운동량의 균형 법칙을 기반으로 한 다상 혼합 모델을 구성하며, 세포(세균 B), 사멸 세포(D), EPS(E), 액체(L)의 네 성분을 포함한다.
• 압축성 없고 뉴턴성 유체를 가정한 혼합 이론을 적용하며, 고체(세균 B, 사멸 세포 D, EPS E)와 액체(L) 단위의 별개 속도를 고려한다.
• 질량 균형 방정식에 반응 항(Γ_B, Γ_D, Γ_E, Γ_L)을 도입하며, 총 부피율(B+D+E+L=1)과 총 질량 보존(ΣΓ_φ=0) 조건을 충족시킨다.
• 빛 의존적 성장은 세균 방정식에 광합성 항을 도입하여 모델링하며, 빛 감쇠는 베어-앰베르트 법칙을 사용한다.
• 성장 속도, 밀도, 응력 텐서 계수 등의 문헌 기반 값에 기반한 매개변수 교정을 수행한다.
• 주요 매개변수(예: 성장 속도, 응력 텐서 계수)의 민감도 분석을 수행하고 2D 및 3D 시뮬레이션을 통해 결과를 검증한다.

실험 결과

연구 질문

• RQ1 모델의 비자명한 평형점 주변에서 일차원에서 선형 안정성 분석을 수행했을 때, 유체역학 모델은 어떻게 행동하는가?
• RQ2 빛 강도와 온도는 모델 내 시안오토프일 생물막의 성장 역학에 어떤 영향을 미치는가?
• RQ3 모델 매개변수 중 어떤 것이 변동에 가장 민감한가? 그리고 매개변수 불확실성 하에서 모델은 얼마나 탄력적인가?
• RQ4 모델은 2D 및 3D 기하 구조에서 실험적으로 관측된 생물막 두께와 빛 감쇠 패턴을 재현할 수 있는가?
• RQ5 실제 환경 조건 하에서 30일 동안 생물막 부피와 구조의 정량적 변화는 어떻게 진행되는가?

주요 결과

• 모델은 일차원에서 유일한 비자명한 평형점 주변에서 선형 안정성을 보이며, 이론적 일관성을 확인한다.
• 빛과 온도 의존성을 통합함으로써 모델의 생물학적 관련성과 예측 능력이 크게 향상된다.
• 민감도 분석 결과, 시안오토프일 생물막의 성장 속도와 응력 텐서 계수가 가장 중요한 매개변수로 규명되었다.
• 30일간의 2D 시뮬레이션은 최종 생물막 부피 0.5980 cm², 평균 두께 0.1196 cm를 산출하여 실험 결과와 일치한다.
• 3D 시뮬레이션은 최종 생물막 부피 0.0227 cm³, 최대 두께 2–3 mm, 평균 두께 0.1135 cm를 기록하였으며, 초기 조건이 정렬된 경우에 해당한다.
• 빛 분포도는 생물막 내에서 강한 감쇠가 발생하고 액체층과 생물막층 사이에 명확한 경계가 존재함을 확인하였으며, 실험 데이터와 일치한다.