S4: Smoothed Solvers for Soft Tissue Simulation

Lead Research Organisation: University of Sheffield
Department Name: Mechanical Engineering

Abstract

In this project a new solution framework and simulation toolkit for modelling soft tissue deformations in time-critical biomedical applications will be established, with the aim of bringing computational biomechanics closer to the clinic. An approach based on the smoothed finite element method (SFEM) will be formulated. The work will cover development of both the numerical solution framework and a high performance computation scheme based on graphics processing units (GPUs). The resulting software library will be released to the biomedical community as open source, to promote dissemination and further development.

Computational biomechanics provides a powerful basis for modelling soft tissues in biomedical applications. In this project, solid biomechanics problems are of interest: analysis and simulation of the motion and mechanical response of deformable solid tissues. Herein, continuum mechanics formalism provides the mathematical basis for analysis, generally formulated as a set of partial differential equations, and the finite element (FE) method is easily the predominant solution approach. With such a framework, one can, in principle, compute the deformations and stresses produced in arbitrarily complicated structures, under the influence of arbitrarily complicated loads. This capability is of central importance in biomechanics. It is also a key enabling technology for initiatives like the Virtual Physiological Human (www.vph-noe.eu) and the Physiome Project (http://physiomeproject.org), which aim, ultimately, to realise the vision of in silico medicine based on personalised computational modelling. These simulation technologies are also essential tools in development of systems for guidance and planning of highly localised and minimally invasive therapies, and in interactive simulators, for example for risk-free surgeon training.

This project is motivated by three key difficulties that inhibit integration of FE-based models of this kind into clinical applications: (i) model construction: FE methods require discretisation of the involved structures into a high quality mesh of "well shaped" elements, which process remains labour intensive and time consuming for complicated biological structures, and which is particularly detrimental when patient-specific models are required; (ii) handling of large deformations: even after a carefully constructed, high quality mesh has been produced, the large deformations that soft tissues may undergo can distort the mesh so much that the solution fails nonetheless; and (iii) computation time: FE methods are computationally intensive, rendering them unsuitable for time-critical applications like surgical guidance and interactive simulators, or limiting the resolution of large scale simulations like bone microstructural models. SFEMs, the focus of developments in this project, are a recent innovation in computational mechanics, which arise from "smoothing" of spatial gradient fields (e.g. strains) over subdomains of the FE mesh. Among other favourable properties, existing formulations are known to reduce mesh sensitivity substantially. In this project, this feature will form the starting point for an algorithm for which meshes are easier to construct in the first place, and which is insensitive to large deformations, subsequently. The approach will also be explicitly formulated to maximise its efficiency in execution on parallel hardware, thus allowing substantial acceleration using cheap and efficient GPUs. By these means, the proposed simulation framework potentially will ameliorate all three of the mentioned difficulties, thus promoting integration of computational biomechanics into clinically-relevant applications.

Planned Impact

S4 will promote national competitiveness in the medical technologies sector. Despite ongoing economic pressure, the UK MedTech sector has defied the recession and experienced continuous growth: according to Department of Business, Innovation and Skills data released in February 2014, sector turnover and employment grew at equivalents (compounded annually) of 2% and 5%, respectively, between 2009 and 2012; in 2012, the sector turned over £17.6bn and employed 76,700 people (the largest sector employer in the UK Life Sciences industry). Globally, the sector is worth £223bn. Given the problems of aging population and concomitant healthcare needs, the market is likely to see even faster growth in future, and the UK's MedTech research base must be supported, if the UK economy is properly to benefit. Simulation already is a central tool in medical device design and its importance will increase as it is further incorporated into clinical workflows (e.g. for surgical planning and guidance) and, in the longer term, as in silico medicine technologies are realised. S4 will support the MedTech industry by overcoming several key challenges of developing patient-specific tissue- and organ-scale simulations. S4's simulation tools will improve soft tissue modelling reliability, significantly accelerate the simulation process itself, and help to automate the process of individualised model generation. Concrete application examples are prosthesis design, for which faster and more accurate assessment of device performance and interactions with tissues will be possible; interactive surgical simulators, in which higher fidelity models of tissue response and motion will be feasible; and surgical planning systems, for which faster and more automated surgical outcomes predictions may be leveraged.

Ultimately, it is the clinical community that drives applications in these areas, and it is they and the patients they serve who will benefit in the medium to longer term. As well as through more immediate improvements in medical devices and systems, improved simulation will benefit clinicians and patients of the future by promoting realisation of the in silico medicine paradigm. The recently published "Roadmap for the Digital Patient" (http://www.digital-patient.net/) sets out a strategy for achievement and widespread adoption of such technologies. These will allow individualisation of clinical decision-making by systematically integrating patient data and personalised physiological models to predict, for example, disease progression and responses to different therapies. From a patient (and, correspondingly, wider societal health) point of view, this is a revolutionary change from the current population-based approach to medicine. At the core of the vision are reliable and efficient patient-specific simulation technologies. Via the methodological advances described above (viz. reliability, acceleration and automation), S4 will support faster translation of simulation-based methods into clinical practice and, in turn, realisation of personalised in silico medicine.

If adopted, S4's outputs will also impact upstream technology developers, such as finite element analysis (FEA) software vendors and GPU manufacturers. The former may achieve competitive advantages by using S4 tools to improve the performance of their software, while the latter would benefit from further penetration of their devices into healthcare and simulation markets. As mentioned previously, the methods developed within S4 certainly have applications in general (non-biomedical) engineering analysis, also, meaning that the impact for the FEA and GPU industries could be significant.

Publications

10 25 50
 
Description Two main outcomes emerged from the project:

1. An efficient and robust method for simulating the mechanical deformation of soft tissues was formulated. The method drew on our group's previous work on explicit dynamic simulation techniques, and extended them with new ideas from so-called smoothed finite element methods. In particular, we produced an original formulation based on bubble enrichment of shape functions. The latter allowed us to circumvent one of the key numerical difficulties when simulating soft tissues: that of 'volumetric locking', which arises from the near incompressibility of the tissues. There exist indeed various approaches to addressing this phenomenon, but all of them either involve significantly higher computational costs - too high for most clinical applications, and certainly for real-time scenarios - , or afford significantly inferior results compared with the new approach. Use of bubble enrichment, moreover, leads to a more elegant and straightforward formulation than was envisaged in the original proposal. Preliminary results using the algorithm were presented at the World Congress on Computational Mechanics, and a full description has been submitted for publication in a high quality international journal (see Publications section).

2. The new simulation method has been incorporated into our open source simulation toolkit NiftySim, as described in the Software & Technical Products (STP) section. NiftySim is a lightweight and compact library of software tools that are aimed at real-time simulation. In particular, it was developed from the beginning to exploit the very high performance of modern Graphics Processing Units (GPUs), as general purpose computational co-processors. The new formulation described above was designed purposely to be highly parallelisable, meaning it can map efficiently to parallel hardware like GPUs, just as do the existing NiftySim algorithms. Optimised parallel execution strategies were developed for this purpose. These tools allow us now to exploit the speed of GPUs, while achieving similar accuracy to the more complex (and slower) methods deployed in commercial finite element packages. As described in the STP section, the new libraries will be released as open source in the next public release of NiftySim, enabling the medical simulation community to leverage the advantages described here.
Exploitation Route The basic work of formulating the method is now complete, and immediate next steps should revolve around exploiting the new tools in clinically-driven applications. We envisage two main application areas: systems and methods for computer-assisted medical treatments (most importantly, model-based techniques for image-guided therapies); and physical modelling for interactive surgical simulators.

We are actively pursuing funding in both streams: we're teaming with colleagues from University College London to develop new systems for optimising and guiding prostate cancer treatments; and working with new clinical partners from Leeds Cancer Research Centre to develop new tools for simulating needle insertion for use in virtual training simulators.

The methods, however, are general, and could form the basis of developments in almost any computer-assisted surgery or surgical simulation application, in which soft tissue deformations play a role (which is most of them). The availability of efficient open source software libraries can facilitate this exploitation.

We also engaged one of the leading European experts on surgical simulation (Prof. Matthias Harders, Uni. Innsbruck) as a scientific advisor on the S4 project, with a view to eventual translation of the new techniques to commercial training systems. This objective is yet premature, but should certainly pursued once further application-focused research, as described, has been undertaken.
Sectors Healthcare

 
Description Sheffield International Mobility Scheme
Amount £10,000 (GBP)
Organisation University of Sheffield 
Sector Academic/University
Country United Kingdom
Start 02/2016 
End 07/2017
 
Description Collaboration with Legato Team, University of Luxembourg 
Organisation University of Luxembourg
Country Luxembourg 
Sector Academic/University 
PI Contribution - Led development of smoothed finite element methods, as described in S4 project documents. - Provided expertise in computational mechanics, and in particular for simulation of soft tissues. - Provided main intellectual input to all aspects of the project. - Undertook development work, including: algorithm design, code implementation, testing, paper writing.
Collaborator Contribution - Expertise in smoothed finite element methods, and other numerical simulation techniques. - Intellectual input to methods development, validation experiments, and paper writing.
Impact Two papers have been the main outputs of the collaboration so far: - Strain-smoothed real-time explicit dynamic (STARTED) algorithm for soft tissue simulation - Strain smoothed explicit dynamic algorithm for soft tissues undergoing large deformations
Start Year 2015
 
Title NiftySim update 
Description New computational methods developed within the S4 project, based on smoothed finite element methods, have been integrated into our existing open source toolkit NiftySim. This enabled us both to leverage the existing infrastructure provided by NiftySim, while substantial enhancing the capabilities of the latter with the new methods. Final testing and refinement of the new codes are nearing completion, and they will be included in the next public NiftySim release. 
Type Of Technology New/Improved Technique/Technology 
Year Produced 2017 
Impact These tools allow fast and robust simulation of nearly incompressible soft tissues using easy-to-generate meshes of low-order tetrahedral elements. The new methods and codes formed an integral part of a new EPSRC proposal submitted in November (awaiting outcome) on simulation of prostate motion for image-guided cancer therapies. 
 
Description CSM-Lux 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact I delivered a talk covering some of the outcomes of the S4 project, and how these are planned to form the basis of a new project on prostate modelling for cancer image-guided therapies. The event was in a workshop format, meaning extensive discussion followed the talk. It is planned to prepare a joint paper with several other workshop participants covering this work and related developments presented by researchers from other fields.
Year(s) Of Engagement Activity 2016
 
Description Schools visit 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach Regional
Primary Audience Schools
Results and Impact My team and I conducted a workshop at Penistone Grammar School aimed at introducing students to the field of computational medicine. In this context, I presented an introductory talk, in which developments from the S4 project played a part. Students then engaged in various practical activities around building computational models of individual patients. The school reported positive feedback, and forwarded an invitation to return the following year.
Year(s) Of Engagement Activity 2016