[论文解读] Benchmark problems for transcranial ultrasound simulation: Intercomparison of compressional wave models

Jean‐François Aubry, Oscar Bates+20 more|arXiv (Cornell University)|Feb 9, 2022

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衍生材料参数的模拟。