{"id":{"repo_id":"cambridge","oai_identifier":"oai:www.repository.cam.ac.uk:1810/380520"},"canonical_url":"https://search.dev.ndltd.org/etd/cambridge/oai:www.repository.cam.ac.uk:1810/380520","repository":{"repo_id":"cambridge","name":"Cambridge University","base_url":"https://api.repository.cam.ac.uk/server/oai/request"},"display":{"title":"Fast finite element methods for focused ultrasound applications","abstract":"This thesis presents the development of a three-dimensional, time-domain, fast finite element solver for focused ultrasound (FUS) applications. FUS is a rapidly evolving medical technology used for tissue heating, ablation, and neuromodulation. In recent years, interest in FUS has grown significantly, especially for the treatment of cancers and brain diseases. Effective patient-specific treatment planning can be supported by accurate simulations. However, this poses a challenge from a computational standpoint since the ratio of the domain size to the wavelength is typically large, therefore computationally demanding. This is made even more challenging for nonlinear models such as the Westervelt model, commonly used in high-intensity focused ultrasound (HIFU) treatments, as higher harmonics need to be resolved. To address this challenge, a three-dimensional, time-domain acoustic solver for FUS applications is developed. The solver is based on the high-order finite element method and the explicit Runge-Kutta method. It is well known that low-order finite element methods suffer from pollution errors, especially for high-frequency problems. However, it is demonstrated that with high-order finite elements, the pollution error is significantly reduced. Moreover, the availability of modern open-source finite element software, such as FEniCSx, eases the implementation of high-order finite element methods. A mass-lumped finite element scheme is implemented, which, combined with the explicit Runge-Kutta method avoids the need to solve a linear system of equations. To further accelerate the solver, a fast finite element integration algorithm is employed. In particular, an algorithm with the lowest known computational cost complexity is implemented, namely the sum-factorisation algorithm. Furthermore, the algorithm is carefully implemented to ensure optimal performance. The solver is implemented for both CPU and GPU architectures. On both architectures, it is shown that the finite element integration implementations achieve a good fraction of peak hardware performance. Additionally, it is shown that the solver exhibits excellent parallel scalability. Contrary to common perceptions, it is demonstrated that a high-order finite element method is well-suited for realistic acoustic simulation for focused ultrasound applications. The novel contributions in this thesis include the development of a high-order finite element method with a fast finite element integration algorithm. The solver is validated against other acoustic wave solvers for several benchmark problems. Overall, it is shown that the developed solver is fast, accurate, scalable, and is suitable for large-scale linear and nonlinear time-domain acoustic simulations.","abstract_html":"This thesis presents the development of a three-dimensional, time-domain, fast finite element solver for focused ultrasound (FUS) applications. FUS is a rapidly evolving medical technology used for tissue heating, ablation, and neuromodulation. In recent years, interest in FUS has grown significantly, especially for the treatment of cancers and brain diseases. Effective patient-specific treatment planning can be supported by accurate simulations. However, this poses a challenge from a computational standpoint since the ratio of the domain size to the wavelength is typically large, therefore computationally demanding. This is made even more challenging for nonlinear models such as the Westervelt model, commonly used in high-intensity focused ultrasound (HIFU) treatments, as higher harmonics need to be resolved. To address this challenge, a three-dimensional, time-domain acoustic solver for FUS applications is developed. The solver is based on the high-order finite element method and the explicit Runge-Kutta method. It is well known that low-order finite element methods suffer from pollution errors, especially for high-frequency problems. However, it is demonstrated that with high-order finite elements, the pollution error is significantly reduced. Moreover, the availability of modern open-source finite element software, such as FEniCSx, eases the implementation of high-order finite element methods. A mass-lumped finite element scheme is implemented, which, combined with the explicit Runge-Kutta method avoids the need to solve a linear system of equations. To further accelerate the solver, a fast finite element integration algorithm is employed. In particular, an algorithm with the lowest known computational cost complexity is implemented, namely the sum-factorisation algorithm. Furthermore, the algorithm is carefully implemented to ensure optimal performance. The solver is implemented for both CPU and GPU architectures. On both architectures, it is shown that the finite element integration implementations achieve a good fraction of peak hardware performance. Additionally, it is shown that the solver exhibits excellent parallel scalability. Contrary to common perceptions, it is demonstrated that a high-order finite element method is well-suited for realistic acoustic simulation for focused ultrasound applications. The novel contributions in this thesis include the development of a high-order finite element method with a fast finite element integration algorithm. The solver is validated against other acoustic wave solvers for several benchmark problems. Overall, it is shown that the developed solver is fast, accurate, scalable, and is suitable for large-scale linear and nonlinear time-domain acoustic simulations.","abstract_has_math":false,"creators":["Mohd Arif Kor, Mohd Adeeb Bin"],"institution":"University of Cambridge","degree_name":null,"degree_level":"Doctoral","degree_discipline":null,"degree_department":null,"school":null,"contributors":[],"advisors":["Wells, Garth"],"committee_chairs":[],"committee_members":[],"year":2024,"date_issued":"2024-08-09","date_published":"2024-08-09","updated_at":"2026-07-22T22:23:56Z","subjects":["Finite element method","Focused ultrasound","Computational acoustic","High-performance computing","Numerical algorithm"],"languages":["eng"],"rights":[],"rights_urls":["https://www.repository.cam.ac.uk/bitstreams/1af91afa-a50c-4c03-ae56-d351a0bd6bfa/download","https://creativecommons.org/licenses/by/4.0/"],"identifier_entries":[]},"links":{"outbound_url":"https://doi.org/10.17863/CAM.116100","outbound_label":"DOI","outbound_source":"dc:identifier.doi"},"metadata_groups":[{"id":"people","label":"People","entries":[{"key":"dc:contributor.advisor","label":"Advisor","values":["Wells, Garth"]},{"key":"dc:contributor.sponsor","label":"Sponsor","values":["Ministry of Higher Education Malaysia Scholarship, University of Cambridge Physics Department"]},{"key":"dc:creator","label":"Author","values":["Mohd Arif Kor, Mohd Adeeb Bin"]}]},{"id":"academic_context","label":"Academic Context","entries":[{"key":"dc:date.issued","label":"Date","values":["2024-08-09"]},{"key":"dc:publisher.institution","label":"Dc Publisher Institution","values":["University of Cambridge"]},{"key":"dc:relation.isreferencedby.uri","label":"Dc Relation Isreferencedby URI","values":["https://www.repository.cam.ac.uk/handle/1810/380520"]},{"key":"dc:type","label":"Dc Type","values":["Thesis"]},{"key":"dc:type.qualificationlevel","label":"Dc Type Qualificationlevel","values":["Doctoral"]}]},{"id":"subjects_keywords","label":"Subjects and Keywords","entries":[{"key":"dc:subject","label":"Dc Subject","values":["Finite element method","Focused ultrasound","Computational acoustic","High-performance computing","Numerical algorithm"]}]},{"id":"language_rights","label":"Language and Rights","entries":[{"key":"dc:language","label":"Dc Language","values":["eng"]},{"key":"dc:rights","label":"Dc Rights","values":["https://www.repository.cam.ac.uk/bitstreams/1af91afa-a50c-4c03-ae56-d351a0bd6bfa/download","https://creativecommons.org/licenses/by/4.0/"]}]},{"id":"identifiers","label":"Identifiers","entries":[{"key":"dc:identifier.doi","label":"DOI","values":["https://doi.org/10.17863/CAM.116100"]},{"key":"dc:identifier.uri","label":"Identifier URI","values":["https://www.repository.cam.ac.uk/bitstreams/6f16b7c4-cdfb-4920-b734-50737b45924a/download"]}]},{"id":"additional","label":"Additional Metadata","entries":[{"key":"dc:description.abstract","label":"Abstract","values":["This thesis presents the development of a three-dimensional, time-domain, fast finite element solver for focused ultrasound (FUS) applications. FUS is a rapidly evolving medical technology used for tissue heating, ablation, and neuromodulation. In recent years, interest in FUS has grown significantly, especially for the treatment of cancers and brain diseases. Effective patient-specific treatment planning can be supported by accurate simulations. However, this poses a challenge from a computational standpoint since the ratio of the domain size to the wavelength is typically large, therefore computationally demanding. This is made even more challenging for nonlinear models such as the Westervelt model, commonly used in high-intensity focused ultrasound (HIFU) treatments, as higher harmonics need to be resolved. To address this challenge, a three-dimensional, time-domain acoustic solver for FUS applications is developed. The solver is based on the high-order finite element method and the explicit Runge-Kutta method. It is well known that low-order finite element methods suffer from pollution errors, especially for high-frequency problems. However, it is demonstrated that with high-order finite elements, the pollution error is significantly reduced. Moreover, the availability of modern open-source finite element software, such as FEniCSx, eases the implementation of high-order finite element methods. A mass-lumped finite element scheme is implemented, which, combined with the explicit Runge-Kutta method avoids the need to solve a linear system of equations. To further accelerate the solver, a fast finite element integration algorithm is employed. In particular, an algorithm with the lowest known computational cost complexity is implemented, namely the sum-factorisation algorithm. Furthermore, the algorithm is carefully implemented to ensure optimal performance. The solver is implemented for both CPU and GPU architectures. On both architectures, it is shown that the finite element integration implementations achieve a good fraction of peak hardware performance. Additionally, it is shown that the solver exhibits excellent parallel scalability. Contrary to common perceptions, it is demonstrated that a high-order finite element method is well-suited for realistic acoustic simulation for focused ultrasound applications. The novel contributions in this thesis include the development of a high-order finite element method with a fast finite element integration algorithm. The solver is validated against other acoustic wave solvers for several benchmark problems. Overall, it is shown that the developed solver is fast, accurate, scalable, and is suitable for large-scale linear and nonlinear time-domain acoustic simulations."]},{"key":"dc:format.checksum.md5","label":"Dc Format Checksum Md5","values":["59f8285504f22dac30da023e55a97a55","87eda9de84448d1f82354d60eee3eb5f"]},{"key":"dc:title","label":"Title","values":["Fast finite element methods for focused ultrasound applications"]}]}],"canonical_facts":{"dc:contributor.advisor":["Wells, Garth"],"dc:contributor.sponsor":["Ministry of Higher Education Malaysia Scholarship, University of Cambridge Physics Department"],"dc:creator":["Mohd Arif Kor, Mohd Adeeb Bin"],"dc:date.issued":["2024-08-09"],"dc:description.abstract":["This thesis presents the development of a three-dimensional, time-domain, fast finite element solver for focused ultrasound (FUS) applications. FUS is a rapidly evolving medical technology used for tissue heating, ablation, and neuromodulation. In recent years, interest in FUS has grown significantly, especially for the treatment of cancers and brain diseases. Effective patient-specific treatment planning can be supported by accurate simulations. However, this poses a challenge from a computational standpoint since the ratio of the domain size to the wavelength is typically large, therefore computationally demanding. This is made even more challenging for nonlinear models such as the Westervelt model, commonly used in high-intensity focused ultrasound (HIFU) treatments, as higher harmonics need to be resolved. To address this challenge, a three-dimensional, time-domain acoustic solver for FUS applications is developed. The solver is based on the high-order finite element method and the explicit Runge-Kutta method. It is well known that low-order finite element methods suffer from pollution errors, especially for high-frequency problems. However, it is demonstrated that with high-order finite elements, the pollution error is significantly reduced. Moreover, the availability of modern open-source finite element software, such as FEniCSx, eases the implementation of high-order finite element methods. A mass-lumped finite element scheme is implemented, which, combined with the explicit Runge-Kutta method avoids the need to solve a linear system of equations. To further accelerate the solver, a fast finite element integration algorithm is employed. In particular, an algorithm with the lowest known computational cost complexity is implemented, namely the sum-factorisation algorithm. Furthermore, the algorithm is carefully implemented to ensure optimal performance. The solver is implemented for both CPU and GPU architectures. On both architectures, it is shown that the finite element integration implementations achieve a good fraction of peak hardware performance. Additionally, it is shown that the solver exhibits excellent parallel scalability. Contrary to common perceptions, it is demonstrated that a high-order finite element method is well-suited for realistic acoustic simulation for focused ultrasound applications. The novel contributions in this thesis include the development of a high-order finite element method with a fast finite element integration algorithm. The solver is validated against other acoustic wave solvers for several benchmark problems. Overall, it is shown that the developed solver is fast, accurate, scalable, and is suitable for large-scale linear and nonlinear time-domain acoustic simulations."],"dc:format.checksum.md5":["59f8285504f22dac30da023e55a97a55","87eda9de84448d1f82354d60eee3eb5f"],"dc:identifier.doi":["https://doi.org/10.17863/CAM.116100"],"dc:identifier.uri":["https://www.repository.cam.ac.uk/bitstreams/6f16b7c4-cdfb-4920-b734-50737b45924a/download"],"dc:language":["eng"],"dc:publisher.institution":["University of Cambridge"],"dc:relation.isreferencedby.uri":["https://www.repository.cam.ac.uk/handle/1810/380520"],"dc:rights":["https://www.repository.cam.ac.uk/bitstreams/1af91afa-a50c-4c03-ae56-d351a0bd6bfa/download","https://creativecommons.org/licenses/by/4.0/"],"dc:subject":["Finite element method","Focused ultrasound","Computational acoustic","High-performance computing","Numerical algorithm"],"dc:title":["Fast finite element methods for focused ultrasound applications"],"dc:type":["Thesis"],"dc:type.qualificationlevel":["Doctoral"]},"updated_at":"2026-07-22T22:23:56Z"}