흐름 관류 조직공학 스캐폴드 전단응력 근사 정확도 비교: 경계조건의 영향

초록

본 연구는 3차원 다공성 스캐폴드에서 흐름 관류에 의해 발생하는 전단응력을 전체 스캐폴드 시뮬레이션과 비교했을 때, 대표 부피 요소(RVE)에 적용되는 두 가지 경계조건(벽 경계조건 WBC, 주기적 경계조건 PBC)의 정확도를 평가한다. Lattice‑Boltzmann Method를 이용해 μCT 기반 스캐폴드 모델을 계산했으며, RVE 크기가 기공 크기의 6배 이상이어도 평균 전단응력 오차가 20~80%에 달한다는 것을 발견했다. 오차는 스캐폴드 기공률이 높을수록 커졌으며, PBC가 WBC보다 일관되게 정확하고 연산 효율이 높아 권장된다.

상세 분석

이 논문은 조직공학에서 흔히 사용되는 흐름 관류 배양 시스템에서 세포에 가해지는 전단응력을 정량화하려는 시도이다. 전체 스캐폴드에 대한 전산유체역학(CFD) 시뮬레이션은 기하학적 복잡성 때문에 계산 비용이 크게 증가한다. 이를 회피하기 위해 연구자들은 전체 구조에서 작은 대표 부피 요소(RVE)를 추출하고, 경계조건을 단순화한다. 여기서는 두 가지 전형적인 경계조건, 즉 물리적 벽을 가정하는 WBC와 인접 RVE와의 연속성을 강제하는 PBC를 비교한다.

Lattice‑Boltzmann Method(LBM)를 선택한 이유는 복잡한 다공성 매체에 대한 흐름 해석에 높은 정확도와 병렬 처리 효율을 제공하기 때문이다. 저자들은 고해상도 μCT 스캔으로부터 3D 스캐폴드 모델을 재구성하고, 기공률을 70%에서 90%까지 변화시키는 여러 샘플을 제작하였다. 각 샘플에 대해 전체 스캐폴드와 RVE(크기 ≥ 6× 평균 기공 직경) 시뮬레이션을 수행하고, 표면 전단응력 분포를 비교하였다.

결과는 두 경계조건 모두 전단응력을 크게 낮게 예측한다는 점에서 일관되었다. 평균 오차는 WBC에서 3080%, PBC에서 2060%로, 특히 기공률이 높을수록 오차가 증가하였다. 이는 높은 기공률일수록 작은 기공이 전체 흐름장을 지배하게 되어, 인위적인 경계가 실제 흐름 패턴을 왜곡시키는 효과가 커짐을 의미한다. 반면, 기공률이 낮은 경우(즉, 구조가 더 촘촘한 경우)에는 작은 기공이 흐름을 제한하므로 경계조건의 부정확성이 상대적으로 감소한다.

연산 효율 측면에서 PBC는 경계면에서의 속도와 압력 연속성을 자동으로 보장함으로써, 추가적인 경계 레이어나 인위적 구속조건을 설정할 필요가 없어 메모리 사용량과 계산 시간을 크게 절감한다. WBC는 고정된 벽면에서 유속이 급격히 감소하는 현상을 강제하기 때문에, 동일한 격자 해상도에서도 더 많은 반복과 수렴 과정을 요구한다.

이러한 결과는 기존 문헌에서 제시된 “RVE 크기는 평균 기공 직경의 6배 이상이어야 한다”는 경험적 가이드라인이 전단응력 정확도 보장을 충분히 설명하지 못함을 시사한다. 실제로, RVE 크기가 충분히 크더라도 경계조건 자체가 흐름장에 비현실적인 제약을 가해 전단응력을 과소평가한다. 따라서 연구자는 PBC를 기본 선택으로 채택하고, 필요 시 경계조건에 따른 보정 계수를 도입하거나, 전체 스캐폴드 시뮬레이션을 병행하는 전략을 권고한다.

이 논문의 의의는 조직공학 분야에서 RVE 기반 전단응력 해석에 대한 정량적 정확도 기준을 최초로 제시하고, 경계조건 선택이 실험 설계와 결과 해석에 미치는 영향을 명확히 밝힌 점이다. 향후 연구에서는 다중 스케일 모델링, 비정상 흐름(펄스 흐름) 및 세포-기질 상호작용을 포함한 복합 모델을 구축함에 있어, PBC 기반 RVE가 제공하는 효율성과 정확성을 기반으로 보다 정교한 생물학적 예측이 가능할 것이다.