초음파 충격파 치료를 위한 재료 계면 전산 모델 연구

초록

본 논문은 체외 충격파 치료(ESWT)에서 물과 골 사이의 파동 전파를 고해상도 유한체적법과 적응형 격자 세분화(AMR)를 이용해 3차원 라그랑지안 형태의 등엔트로피 오일러 방정식(유체)과 선형 탄성 방정식(골)으로 시뮬레이션한다. 복잡한 골 형상, 계면에서의 반사·모드 변환, 골 내부 전단응력 발생을 정량적으로 분석하고, 모델의 단순화 가정과 검증 결과를 논의한다.

상세 분석

이 연구는 ESWT의 물 매질에서 발생한 충격파가 골 조직에 전달되는 물리적 메커니즘을 수치적으로 재현하는 데 초점을 맞춘다. 저자들은 유체를 등엔트로피(비열이 일정) 가정하에 라그랑지안 형태의 오일러 방정식으로 기술하고, 골은 선형 탄성체로서 변형률-응력 관계를 선형 탄성 방정식으로 모델링한다. 라그랑지안 형태를 선택함으로써 고체와 유체 사이의 경계에서 변위와 속도 연속성을 자연스럽게 적용할 수 있다.

수치 해법은 고해상도 유한체적법을 기반으로 하며, 충격파와 같은 급격한 불연속을 정확히 포착하기 위해 TVD(총변동감소) 제한자를 사용한다. 또한, 적응형 격자 세분화(AMR)를 도입해 충격 전파 경로와 계면 근처에만 격자를 집중시켜 계산 비용을 크게 절감하면서도 해상도를 유지한다. 특히 3차원 계산을 수행함으로써 실제 골 형상의 복잡성을 반영할 수 있다.

계면에서 발생하는 반사와 전이(모드 변환)는 유체와 고체의 파동 속도 차이와 경계 조건에 의해 결정된다. 저자들은 압축파가 골에 입사할 때 전단파가 생성되는 현상을 관찰했으며, 이는 골 내부에 전단응력이 축적되어 조직 재생을 촉진할 가능성을 시사한다. 전단응력의 분포와 크기는 골의 기하학적 특성(곡률, 두께)과 물-골 경계의 매끄러움에 민감하게 반응한다.

모델 검증을 위해 실험적 충격파 프로파일과 기존 문헌의 전파 속도 데이터를 비교했으며, 전반적인 파형과 압력 피크가 실험과 일치함을 확인했다. 그러나 등엔트로피 가정, 선형 탄성 모델, 그리고 점성·비선형 효과를 무시한 단순화가 고주파·고압 영역에서 오차를 야기할 수 있음을 인정한다. 이러한 제한에도 불구하고, 제안된 전산 프레임워크는 ESWT의 파동-구조 상호작용을 정량적으로 분석하고, 치료 파라미터 최적화 및 개인 맞춤형 시뮬레이션에 활용될 수 있는 기반을 제공한다.