A Coupled Fluid Naghdi Shell Windkessel Model for Pulmonary Artery Banding: Partitioned Numerical Scheme and Parameter Estimation

기술 요약: 폐동맥 결찰술을 위한 결합된 유체-Naghdi 쉘-Windkessel 모델

문제 정의 및 동기
본 연구는 선천성 심장 질환, 특히 거대 심실 중격 결손이 있는 신생아에게 과도한 폐혈류를 제한하기 위해 사용되는 완화적 수술 기법인 폐동맥 결찰술(PAB)과 관련된 혈역학적 및 기계적 과제를 다룹니다. 현재의 임상 관행은 경험적인 크기 결정과 수술 중 측정에 의존하고 있으며, 이는 최적의 결찰 구성이나 원위부 폐 순환의 특성을 고유하게 결정할 수 없습니다. 기존의 계산 모델들은 종종 강체 벽 가정을 사용하거나 규정된 출구 조건을 사용하여, 동맥벽의 변형과 혈류 사이의 피드백을 포착하지 못하거나 국소적인 결찰 효과와 하류 혈관 임피던스를 구분하지 못하는 한계가 있습니다.

저자들은 결찰의 국소적인 기계적 효과를 해결되지 않은 하류 혈관의 영향으로부터 분리하기 위한 수학적 프레임워크를 제안합니다. 다루고자 하는 구체적인 과제는 다음과 같습니다:

1. 국소적 수축(불연속적인 표면 기울기)이 있는 조각별로 매끄러운 기하학적 구조를 처리할 수 있는 얇고 변형 가능한 구조로서 동맥벽을 모델링하는 것.
2. 3D 유체 영역을 2D 쉘 모델 및 0D 하류 순환 모델(Windkessel)과 결합하되, 첨가 질량 불안정성(added-mass instability)을 겪지 않도록 하는 것.
3. 제한된 입구 압력 관측값으로부터 미지의 하류 혈관 파라미터(근위 저항, 원위 저항, 컴플라이언스)를 추정하는 것.

방법론
본 논문은 0D Windkessel 모델과 통합된 결합된 3차원(3D) 유체-구조 상호작용(FSI) 모델을 개발합니다.

• 유체 영역: 혈류는 움직이는 벽에 의한 유체 영역의 변형을 처리하기 위해 Arbitrary Lagrangian-Eulerian (ALE) 프레임워크 내에서 비압축성 나비에-스토크스(Navier-Stokes) 방정식을 사용하여 모델링됩니다.
• 구조 영역 (Shell): 동맥벽은 2D 참조 영역에 정의된 3D Naghdi 쉘 모델로 표현됩니다. 고전적인 Koiter 쉘과 달리, Naghdi 정식화는 중간면 변위와 쉘 디렉터 회전을 독립적인 미지수로 취급합니다. 이를 통해 모델은 막(membrane), 굽힘(bending), 그리고 횡방향 전단(transverse-shear) 변형 메커니즘을 유지할 수 있습니다. 결정적으로, 이 정식화는 참조 중간면의 $W^{1,\infty}$ 정규성만을 요구하므로, 결찰의 전이 구역(geometric corners)과 같은 조각별로 매끄러운 기하학적 구조를 모델링할 수 있습니다.
• 하류 순환: 원위부 폐 순환은 근위 저항($R_p$), 원위 저항($R_d$), 그리고 컴플라이언스($C$)를 통해 유출 유량과 압력을 연결하는 3-요소 Windkessel 모델을 사용하여 모델링됩니다.
• 결합 방식: 시스템은 비반복적(non-iterative)으로 운동학적으로 결합된 분할 기법(partitioned scheme)을 사용하여 해결됩니다. 첨가 질량 불안정성을 완화하기 위해, 이 기법은 Robin 유형의 인터페이스 조건을 통해 쉘의 병진 관성을 유체 부문제에 포함합니다. 알고리즘은 각 타임 스텝마다 세 단계로 진행됩니다: (1) 이전 유체 트랙션의 일부를 사용하여 구조 부문제를 풀고, (2) ALE 메쉬를 업데이트하며, (3) 업데이트된 도메인에서 유체와 Windkessel 부문제를 함께 풉니다.
• 파라미터 추정: 하류 파라미터의 불확실성을 해결하기 위해 축소된 상태의 무향 칼만 필터(reduced-state unscented Kalman filter, UKF)가 채택되었습니다. 필터 상태는 양수성을 보장하기 위해 세 가지 Windkessel 파라미터($R_p, R_d, C$)의 로그값으로 구성됩니다. 필터는 접선 또는 어드조인트 방정식의 필요 없이, 입구 압력의 관측값과 시뮬레이션 값 사이의 차이를 바탕으로 이러한 추정치를 업데이트합니다.

주요 기여

1. 저정규성 Naghdi 쉘 정식화: 저자들은 불연속적인 기울기를 가진 조각별로 매끄러운 중간면을 수용할 수 있는 Naghdi 쉘의 약형식(weak formulation)을 유도했습니다. 이는 결찰의 기하학적 전이를 정확하게 모델링하는 데 필수적입니다.
2. 안정적인 분할 기법: Robin 인터페이스 조건을 통해 쉘의 관성을 포함하는 비반복적 분할 기법이 개발되었습니다. 선형화된 고정 도메인 문제에 대한 이산 에너지 추정치가 도출되었으며, 이는 분할 파라미터 $\beta \in [0, 1]$에 대해 무조건적인 안정성을 입증합니다.
3. 축소된 상태의 데이터 동화: 축소된 상태의 UKF 통합을 통해 입구 압력 데이터만으로 하류 혈관 특성을 추정할 수 있으며, 이를 통해 국소적 결찰 효과를 원위부 순환의 불확실성과 분리합니다.
4. 비교 분석: 표준 Naghdi 모델을 수정된, 주로 반경 방향(radial-dominant)인 Koiter형 쉘 모델과 비교하여 횡방향 전단 및 축 방향 변형이 압력파 전파 속도에 미치는 영향을 정량화했습니다.

수치 결과
프레임워크는 세 가지 수치 실험을 통해 검증되었습니다:

• 직선 튜브 벤치마크: 구현을 검증하기 위해 3D 확장된 압력 펄스 벤치마크를 수행했습니다. 공간 및 시간 수렴 연구 모두 일관된 압력파 전파와 벽 변형을 보여주었습니다.
• 모델 비교: 표준 Naghdi 쉘과 펜널티(penalized)가 적용된 반경 방향 우세 모델을 비교한 결과, 두 모델 모두 예상된 압력파 및 벽 팽창 거동을 재현했으나, 수정된(펜널티 적용) 모델이 표준 Naghaki 모델보다 더 빠른 압력 전면(pressure front)을 생성함을 확인했습니다. 이론적 분석에 따르면, 3D 모델은 동일한 반경 방향 강성 가정 하에서 2D 모델보다 약 $1/\sqrt{2}$ 배 느린 파동 속도를 보입니다.
• 폐동맥 결찰 및 추정: 이상적인 결찰 폐동맥을 Windkessel 출구와 결합했습니다. 축소된 상태의 UKF는 일정한 입구 압력 관측으로부터 근위 저항, 원위 저항, 컴플라이언스를 성공적으로 추정했습니다. 결과는 관측된 압력 강하를 맞추기 위해 추정된 근위 저항이 크게 조정된 것(약 22배 증가)을 보여준 반면, 원위 저항과 컴플라이언스는 상대적으로 안정적이었습니다. 시뮬레이션된 입구 압력은 목표 관측값의 2.69% 이내로 수렴했습니다.

의의 및 주장
저자들은 본 연구를 환자 맞춤형 폐동맥 결찰 최적화를 위한 기초적인 단계로 위치시킵니다. 그들은 제안된 프레임워크가 원위부 혈관 특성을 명시적으로 고려하면서 국소적인 혈관 수축의 혈역학적 및 기계적 효과를 조사할 수 있는 토대를 제공한다고 주장합니다.

저자들은 수치 실험을 완성된 환자 맞춤형 최적화 연구라기보다는 결합 솔버 및 파라미터 추정 프레임워크에 대한 "개념 증명(proof of concept)"으로 겸허하게 규정합니다. 본 연구의 의의는 국소적 결찰 역학과 원위부 임피던스를 구분할 수 있는 능력에 있으며, 저자들은 이것이 기존의 강체 벽 또는 규정된 출구 모델에는 결여된 기능이라고 주장합니다. 저자들이 구상하는 향후 응용 분야에는 맥동하는 환자 특정 데이터와 영상 기반 기하학적 구조를 사용하여 결찰 구성 및 동맥벽 특성을 체계적으로 조사하는 것이 포함됩니다.