스트레스 구동 이방성 확산과 활성 변형 매체에서의 역할

초록

본 논문은 기계적 응력이 확산 텐서를 직접 조절한다는 가정 하에, 활성 변형 매체, 특히 심장 조직에서의 전기‑기계 피드백을 모델링한다. 스트레스‑보조 확산 텐서를 도입한 반응‑확산‑역학 방정식을 제시하고, 혼합‑프라이멀 유한요소법으로 수치해석을 수행한다. 실험적 증거와 시뮬레이션을 통해 응력에 의해 초기의 등방성·균일한 확산 텐서가 이방성·비균일하게 변하고, 이는 파동 전도 속도와 전도 방향에 중대한 영향을 미침을 확인한다.

상세 분석

이 연구는 기존의 전기‑기계 결합 모델이 확산(전도) 텐서를 변형에 무관하게 고정된 상수로 취급한 점을 비판하고, 물리적·수학적 일관성을 확보하기 위해 스트레스‑보조 확산 텐서 dᵢⱼ = D₀(δᵢⱼ + D₁σᵢⱼ + D₂σᵢₖσₖⱼ) 를 제안한다. 여기서 σᵢⱼ는 Cauchy 응력 텐서이며, D₁, D₂는 실험적으로 조정 가능한 계수이다. 이 식은 Landau‑Lifshitz의 비선형 전기역학 이론을 기반으로 하여, 무한소 변형에서도 전기적 물성(전도도)이 변한다는 고전적 결과를 일반화한다.

수학적으로는 확산 텐서가 양정정(positive‑definite)이어야 하는 타당성(ellipticity) 조건을 도출한다. 2차 다항식 y = 1 + D₁σ + D₂σ²가 모든 σ에 대해 양수이어야 하며, 이는 D₁과 D₂ 사이의 불평등 ‑D₁ ± √(D₁²‑4D₂) ≥ 0 등으로 표현된다. 이러한 조건을 만족하지 않으면 확산 연산자가 비타원형이 되어 수치 해석이 발산하거나 물리적으로 비현실적인 전도 현상이 발생한다.

모델 구현에서는 비선형 Mooney‑Rivlin 재료를 사용해 비압축성( J=1) 유한 변형을 기술하고, 활성 응력 Tₐ를 Heaviside 함수 Γ(V) 로 전압에 의존적으로 전환한다. 반응‑확산 부분은 FitzHugh‑Nagumo 형태의 두 변수(RD) 시스템을 채택했으며, 전압 V와 회복 변수 r의 동역학을 기존 모델과 동일하게 유지한다.

수치 해석은 혼합‑프라이멀 유한요소법을 적용해 응력·변형 텐서를 직접 계산한다. Raviart‑Thomas 4 차원 요소와 BDM(1) 요소를 각각 응력·변위에 사용함으로써, 후처리 없이도 정확한 텐서 값을 얻는다. 시간 적분은 1차 후진 오일러와 고정점 반복을 결합해 전기‑역학 연산을 분리하고, 비선형 역학 방정식은 Newton 방법으로 해결한다.

시뮬레이션 결과는 두 가지 주요 현상을 보여준다. 첫째, 스트레칭에 의해 확산 텐서가 이방성으로 변하면서 파동 전도 속도가 스트레칭 방향에 비해 크게 증가하거나 감소한다. D₁, D₂가 양수일 때 전도 속도는 x‑방향(스트레칭 방향)에서 최대 5배까지 빨라지는 반면, 음수 조합에서는 오히려 감소한다. 둘째, 특정 파라미터 영역에서는 타원성 조건이 위배되어 전도 파동이 멈추거나 비정상적인 전도 패턴(예: 파동 회전, 전도 불연속)으로 전환된다. 이는 스트레스‑보조 확산이 전도 불안정성(예: 심방세동)과 직접 연관될 수 있음을 시사한다.

또한, 실험적 증거로 스폰지를 한쪽으로 늘렸을 때 염료가 원형이 아닌 타원형으로 퍼지는 현상을 제시한다. 비록 생물학적 조직을 완벽히 재현한 것은 아니지만, 기계적 변형이 확산 경로를 왜곡한다는 직관적 근거를 제공한다.

전체적으로 이 논문은 스트레스‑보조 확산 텐서를 도입함으로써, 기존 전기‑기계 모델이 놓쳤던 중요한 피드백 메커니즘을 수학적으로 정형화하고, 물리적 실험과 수치 시뮬레이션을 통해 그 타당성을 검증한다. 이는 심장 전기역학 모델링뿐 아니라, 변형이 중요한 다른 연성 매체(연골, 조직 공학 스캐폴드 등)에도 적용 가능성을 열어준다.