[논문 리뷰] Benchmark problems for transcranial ultrasound simulation: Intercomparison of compressional wave models

Jean‐François Aubry, Oscar Bates외 20명|arXiv (Cornell University)|2022. 02. 09.

Ultrasound Imaging and Elastography인용 수 4

한 줄 요약

이 논문은 복합성 증가하는 뼈 흉내 모형을 통과하는 종방향 파동 전파를 이용해 11종의 다른 수치 모델을 평가하는 표준화된 벤치마크 세트를 제시한다. 상호 비교 결과, 모델 간에 뛰어난 일치도를 보였으며, 초점 압력의 중앙값 차이는 10% 미만, 초점 위치의 중앙값 차이는 1 mm 미만이었으며, 이는 집중 초음파 치료 및 단계 보정과 같은 임상 적용 분야에서 수치 시뮬레이션의 신뢰성을 검증한다.

ABSTRACT

Computational models of acoustic wave propagation are frequently used in transcranial ultrasound therapy, for example, to calculate the intracranial pressure field or to calculate phase delays to correct for skull distortions. To allow intercomparison between the different modeling tools and techniques used by the community, an international working group was convened to formulate a set of numerical benchmarks. Here, these benchmarks are presented, along with intercomparison results. Nine different benchmarks of increasing geometric complexity are defined. These include a single-layer planar bone immersed in water, a multi-layer bone, and a whole skull. Two transducer configurations are considered (a focused bowl and a plane piston), giving a total of 18 permutations of the benchmarks. Eleven different modeling tools are used to compute the benchmark results. The models span a wide range of numerical techniques, including the finite-difference time-domain method, angular-spectrum method, pseudospectral method, boundary-element method, and spectral-element method. Good agreement is found between the models, particularly for the position, size, and magnitude of the acoustic focus within the skull. When comparing results for each model with every other model in a cross comparison, the median values for each benchmark for the difference in focal pressure and position are less than 10\% and 1 mm, respectively. The benchmark definitions, model results, and intercomparison codes are freely available to facilitate further comparisons.

연구 동기 및 목표

• 경두개 초음파 전파의 수치 모델 평가를 위한 표준화된 벤치마크 프레임워크 수립.
• 경두개 초음파 시뮬레이션에 사용되는 다양한 수치 방법의 일관성과 정확도 평가.
• 다양한 모델링 도구와 기법 간 시뮬레이션 결과의 변동성 원인 규명.
• 집중 초음파 치료 및 뇌내압 예측을 위한 수치 모델링 최선의 실천 방안 지원.
• 재현성과 향후 신규 솔버 검증을 위해 자유롭게 이용 가능한 벤치마크, 데이터, 코드 제공.

제안 방법

• 기하학적 복잡도가 증가하는 아홉 가지 벤치마크 문제 정의: 평면 및 다층 뼈 모형, 전체 뇌골 모형을 포함하며, 모두 물속에 잠겨 있음.
• 두 가지 초음파 전도체 구성—집중형 볼과 평판 피스톤—사용으로 총 18개의 벤치마크 조합 생성.
• 11종의 다른 모델링 도구 활용: 유한차분시간영역, 각도스펙트럼, 가상스펙트럼, 경계요소, 스펙트럼요소 방법 포함.
• 각 모델이 음압장 계산하며, 결과는 초점 위치, 초점 크기, 압력 크기 등의 지표를 사용해 비교.
• 모든 모델 쌍 간의 상호 비교를 수행해 모든 벤치마크에서 모델 간 일치도 정량화.
• 모든 벤치마크 정의, 시뮬레이션 결과, 상호 비교 코드를 공개하여 재현성과 확장성 확보.

실험 결과

연구 질문

• RQ1 다양한 복잡도의 경두개 뼈 구조를 통과하는 초음파를 시뮬레이션할 때, 종방향 파동 전파의 다양한 수치 모델은 어떻게 성능을 발휘하는가?
• RQ2 두개골 통과 후 뇌 내부에서의 초점 압력장 예측에서, 다양한 모델링 도구 간에 어떤 정도의 일치도를 달성할 수 있는가?
• RQ3 수치 방법과 구현 선택의 차이가 초점 위치 및 압력 크기와 같은 주요 시뮬레이션 지표에 어떻게 영향을 미치는가?
• RQ4 모델 간의 오차는 기하학적 복잡도와 초음파 전도체 구성에 따라 어떻게 변화하는가?
• RQ5 표준화된 벤치마크 프레임워크는 경두개 초음파 치료의 계산 예측에 있어 불확실성을 줄이고 신뢰도를 높일 수 있는가?

주요 결과

• 모든 벤치마크에서 어떤 모델도 다른 모델과 비교할 때, 초점 압력의 중앙값 차이가 10% 미만, 초점 위치의 중앙값 차이가 1 mm 미만으로 나타났다.
• 특히 두개골 통과 후 음압 초점의 위치, 크기, 크기에서 양호한 일치도를 관찰했다.
• 전체 필드 지표에서의 더 큰 오차는 주로 전도체와 두개골 사이의 근접 필드 파동 거동의 차이로 인해 발생했으며, 이는 주로 이들 영역에서 발생했다.
• 유한차분, 스펙트럼, 경계요소 방법을 포함한 다양한 수치 기법 간에 강력한 일관성이 입증되었다.
• 벤치마크 정의, 시뮬레이션 결과, 상호 비교 코드가 공개되어 있어 새로운 모델링 도구의 재현 및 향후 검증이 가능하다.
• 본 연구는 향후 상호 비교의 기초가 되며, 탄성파 모델 및 환자 특화 CT 유도 물성 매개변수를 사용한 시뮬레이션 포함 가능.