Patient-Specific 3-Dimensional Model for High-Intensity Focused Ultrasound Treatment Through the Rib Cage: A Preliminary Study

A comprehensive numerical procedure for high-intensity focused ultrasound ablation of breast tumour on an anatomically realistic breast phantom

High-Intensity Focused Ultrasound (HIFU) as a promising and impactful modality for breast tumor ablation, entails the precise focalization of high-intensity ultrasonic waves onto the tumor site, culminating in the generation of extreme heat, thus ablation of malignant tissues. In this paper, a comprehensive three-dimensional (3D) Finite Element Method (FEM)-based numerical procedure is introduced, which provides exceptional capacity for simulating the intricate multiphysics phenomena associated with HIFU. Furthermore, the application of numerical procedures to an anatomically realistic breast phantom (ARBP) has not been explored before. The integrity of the present numerical procedure has been established through rigorous validation, incorporating comparative assessments with previous two-dimensional (2D) simulations and empirical data. For ARBP ablation, the administration of a 0.1 MPa pressure input pulse at a frequency of 1.5 MHz, sustained at the focal point for 10 seconds, manifests an ensuing temperature elevation to 80°C. It is noteworthy that, in contrast, the prior 2D simulation using a 2D phantom geometry reached just 72°C temperature under the identical treatment regimen, underscoring the insufficiency of 2D models, ascribed to their inherent limitations in spatially representing acoustic energy, which compromises their overall effectiveness. To underscore the versatility of this numerical platform, a simulation of a more clinically relevant HIFU therapy procedure has been conducted. This scenario involves the repositioning of the ultrasound focal point to three separate lesions, each spaced at 3 mm intervals, with ultrasound exposure durations of 6 seconds each and a 5-second interval for movement between focal points. This approach resulted in a more uniform high-temperature distribution at different areas of the tumour, leading to the ablation of almost all parts of the tumour, including its verges. In the end, the effects of different abnormal tissue shapes are investigated briefly as well. For solid mass tumors, 67.67% was successfully ablated with one lesion, while rim-enhancing tumors showed only 34.48% ablation and non-mass enhancement tumors exhibited 20.32% ablation, underscoring the need for multiple lesions and tailored treatment plans for more complex cases.

Anatomic thermochromic tissue-mimicking phantom of the lumbar spine for pre-clinical evaluation of MR-guided focused ultrasound (MRgFUS) ablation of the facet joint

OBJECTIVE

To develop a thermochromic tissue-mimicking phantom (TTMP) with an embedded 3D-printed bone mimic of the lumbar spine to evaluate MRgFUS ablation of the facet joint and medial branch nerve.

MATERIALS AND METHODS

Multiple 3D-printed materials were selected and characterized by measurements of speed of sound and linear acoustic attenuation coefficient using a through-transmission technique. A 3D model of the lumbar spine was segmented from a de-identified CT scan, and 3D printed. The 3D-printed spine was embedded within a TTMP with thermochromic ink color change setpoint at 60 °C. Multiple high energy sonications were targeted to the facet joints and medial branch nerve anatomical location using an ExAblate MRgFUS system connected to a 3T MR scanner. The phantom was dissected to assess sonication targets and the surrounding structures for color change as compared to the expected region of ablation on MR-thermometry.

RESULTS

The measured sound attenuation coefficient and speed of sound of gypsum was 240 Np/m-MHz and 2471 m/s, which is the closest to published values for cortical bone. Following sonication, dissection of the TTMP revealed good concordance between the regions of color change within the phantom and expected areas of ablation on MR-thermometry. No heat deposition was observed in critical areas, including the spinal canal and nerve roots from either color change or MRI.

CONCLUSION

Ablated regions in the TTMP correlated well with expected ablations based on MR-thermometry. These findings demonstrate the utility of an anatomic spine phantom in evaluating MRgFUS sonication for facet joint and medial branch nerve ablations.

In silico feasibility assessment of extracorporeal delivery of low-intensity pulsed ultrasound to intervertebral discs within the lumbar spine

Low intensity pulsed ultrasound (LIPUS) may have utility for non-invasive treatment of discogenic lower back pain through stimulating, remodeling and accelerating healing of injured or degenerated intervertebral disc (IVD) tissues. This study investigates the feasibility of delivering LIPUS to lumbar IVDs between L2 and S1 spine vertebra using a planar extracorporeal phased array (8 × 8 cm, 1024 elements, 500 kHz). Three 3D anatomical models with heterogenous tissues were generated from patient CT image sets and used in the simulation-based analysis. Time-reversal acoustic modeling techniques were applied to optimize posterior-lateral placement of the array with respect to the body to facilitate energy deposition in discrete target regions spanning the annulus fibrosus and central nucleus of each IVD. Forward acoustic and biothermal simulations were performed with time-reversal optimized array placements and driving amplitude/phase settings to predict LIPUS intensity distributions at target sites and to investigate off-target energy deposition and heating potential. Simulation results demonstrate focal intensity gain of 5-168 across all IVD targets and anatomical models, with greater average intensity gain (>50) and energy localization in posterior, posterolateral, and lateral target sites of IVDs. Localized LIPUS delivery was enhanced in thinner patient anatomies and in the high lumbar levels (L2-L3 and L3-L4). Multiple amplitude/phasing illumination patterns could be sequenced at a fixed array position for larger regional energy coverage in the IVD. Biothermal simulations demonstrated that LIPUS-appropriate exposures of 100 mW cm^-2 I_SPTA to the target disc region would result in <1 °C global peak temperature elevation for all cases. Hence, simulations suggest that spatially-precise extracorporeal delivery of therapeutically relevant LIPUS doses to discrete regions of lumbar IVDs is feasible and may be useful in clinical management of discogenic back pain.