산소 공급과 신경 전기활동의 결합 모델링 및 수치 해석

초록

본 논문은 저산소증(허혈) 상황에서 산소 공급이 신경세포의 흥분성에 미치는 영향을 다중스케일 수학 모델로 구현한다. 혈류 감소와 세포 부종을 산소 농도와 연계한 부피 변동 법칙을 도입하고, 이를 모노도메인 전기 모델과 이온·대사 모델에 결합한다. 고차 p‑적응 불연속 Galerkin(DG) 방법과 Crank‑Nicolson 시간 적분을 이용해 2차원 정사각형 및 실제 뇌 형상에서 급속 파동과 복잡한 기하를 효율적으로 시뮬레이션한다. 결과는 경미한 허혈부터 심각한 허혈까지 다양한 병리 시나리오에서 전기 활동과 대사 농도의 변화를 정량적으로 제시한다.

상세 분석

이 연구는 뇌 허혈이라는 복합 병리를 수학적으로 정량화하기 위해 세 가지 핵심 요소를 통합한다. 첫째, 혈류 감소와 산소 농도에 의존하는 부피 변동 법칙을 통해 세포 부종을 기계적으로 모델링한다. 이 법칙은 현지 산소 농도와 혈류 감소 비율을 입력으로 받아 세포 내·외부 부피 비율을 동적으로 조정함으로써, 부종이 전기 전도도와 확산 토리오시티에 미치는 영향을 직접적으로 반영한다. 둘째, 전기‑대사 결합 모델은 기존 Barreto‑Cressman 이온 전도 모델에 astrocyte의 K⁺ 흡수와 ATP 의존성 Na⁺/K⁺ 펌프 조절을 추가한다. 대사 모델은 혈액, 세포외공간(ECS), 뉴런, astrocyte 네 개의 구획으로 구성된 다중 구획 ODE 체계이며, 산소 전달을 Fick‑법칙과 Hill식으로 기술한다. 특히 ATP와 ADP, NADH/NAD⁺ 비율을 상태 변수로 두어, 전기 활동이 에너지 소모를 통해 산소 고갈을 촉진하고, 반대로 산소 공급이 ATP 재생을 통해 이온 펌프 활성을 유지하도록 상호 피드백을 구현한다. 셋째, 수치 해석 측면에서는 고차 p‑적응 불연속 Galerkin(DG) 방법을 polygonal/polyhedral 격자에 적용함으로써 복잡한 뇌 해부학적 경계와 급격한 전위 파동을 효율적으로 포착한다. 공간 적응은 전위 급등 구역에만 차수를 높여 연산 비용을 최소화하고, 시간 적응은 Crank‑Nicolson 스키마와 결합해 ms 수준의 전기 시간 스케일과 초~분 수준의 대사 시간 스케일을 동시에 안정적으로 통합한다. 수치 실험에서는 이상적인 정사각형 도메인과 실제 MRI 기반 뇌 형상 두 가지를 사용했으며, 혈류 감소 비율을 10 %에서 70 %까지 변화시켜 경미한 저산소증부터 심각한 허혈까지의 전기·대사 반응을 비교하였다. 결과는 혈류 감소가 클수록 세포 부피가 팽창해 ECS가 축소되고, K⁺ 확산 저항이 증가해 전위 파동 전도가 비정상적으로 빨라지며, 동시에 ATP 고갈로 인해 Na⁺/K⁺ 펌프 효율이 급격히 감소해 지속적인 탈분극 상태가 유발된다는 것을 보여준다. 이러한 다중 스케일·다중 물리 결합 모델은 허혈성 발작, 확산 억제, 그리고 치료적 재관류 전략의 효과를 정량적으로 예측하는 데 유용한 도구가 될 것으로 기대된다.