Access the full text.
Sign up today, get DeepDyve free for 14 days.
ecueto@unizar.es Aragón Institute of Engineering In this paper a non-exhaustive review is made on the existing literature for real-time Research, Universidad de Zaragoza, simulation in the field of computational surgery. Many methods have been proposed Zaragoza, Spain Full list of author information is so far to deal with the very astringent assumption of real-time response in the field, available at the end of the article specially for simulators equipped with haptic peripherals. A special emphasis is made on techniques that respond to the so-called second generation of surgery simulators, that able to adequately model the mechanics of the problem. Techniques employing supercomputing facilities, notably those base upon parallel implementations on GPUs will be covered, while special attention will be paid to techniques based upon model order reduction, a promising technique in the field. Finally, some review is made on techniques able to give some insight in the so-called third generation of surgery simulators, i.e., that able to include physiological details into the simulation. Keywords: Numerical simulation; Computational surgery; Real time; Model order reduction Introduction Computational surgery has been defined recently [1] as “the application of mathematics and algorithm design, enabling imaging, robotics, informatics, and simulation technolo- gies, incorporating biological and physical principles, to improve surgery”. It has been since the early times of development of computers that their promising use in the field of medicine has been investigated [2,3]. With the irruption of endoscopical and minimally invasive procedures, on one side, and robot-operated surgery, on the other, training of surgeons has becoming a task in which computer simulation has acquired a preeminent role in recent years [4]. In essence, surgeons have begun to access organs during surgery in an indirect way, through a screen, and this needs for a period of intensive training to avoid costly errors. R. Satava [5] proposed some fifteen years ago a taxonomy of virtual anatomy that he divided into five different generations. The first generation was composed by systems representing accurately the geometry of the organs at a macroscopic level. The second generation would include an accurate description of the physical dynamics of the body. While it is still hard, more than a decade after, to find a real-time surgical simulator that incorporates accurate, state-of-the-art models for soft tissues at a continuum level, this taxonomy included three more generations. From the third to the fifth one, these virtual © 2014 Cueto and Chinesta; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 2 of 18 http://www.amses-journal.com/content/1/1/11 descriptions of the patient should include, respectively, accurate descriptions of physiol- ogy, microscopic anatomy (at a neurovascular level, for instance) and, finally, biochemical systems. Many computational methods can be found in the biomechanics literature that deal with different levels of description, see for instance [4,6-11] among many others. The last three, however, have not been yet fully incorporated into virtual reality simulators due to the impressive computational cost that they involve, since they incorporate descriptions at the molecular level. To understand the challenge of including state-of-the-art constitutive equations in the simulator, let us consider the structure of a laparoscopic surgery simulator for training purposes [4]. In essence, a virtual surgery simulator is composed by, see Figure 1, a haptic peripheral that communicates with a simulation engine (see also Figure 2). This simula- tion core is responsible of detecting contact between the virtual tool and the organ and, by integrating the stresses produced by such a contact, give back to the user hand a realistic sensation of touch. This is precisely one of the bottlenecks of such systems: to provide the user with a realistic sensation of touch some 500 Hz-1 kHz of feedback response (which coincides with the free hand gesture frequency) should be achieved [4]. For a detailed explanation on haptics the reader can consult [12,13]. It can be easily noticed how this feedback rate becomes a true challenge when we deal, as in the case of living soft tissues, with large deformation settings, frequently in a hyperelastic framework [14,15]. This type of simulations take hours in today’s computers, while we are now requiring to be solved from 500 to one thousand times per second. Sometimes real-time constraints are less astringent, and require lower feedback rates but, on the contrary, are expected to provide long-time responses such as, for instance, a simulated outcome of surgery. This is typical in surgery planning simulators [16,17], that should provide a forecast of the results of a given decision taken in the operating room. In augmented reality applications, for instance, only visual continuity is required, and then feedback at some 25 Hz are enough to perceive this sensation of continuous movement [18]. In this review paper a chronological order is pursued, that coincides roughly with the generations of surgical simulators defined in [5]. It can be said that a true example of second-generation simulator, that which includes state-of-the-art constitutive modeling for soft living tissues, has never been reported in the literature. However, some attempts Figure 1 Schematic architecture of a real-time simulator for surgery. 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 3 of 18 http://www.amses-journal.com/content/1/1/11 Figure 2 Example of a surgery simulator equipped with a haptic peripheral. have been made to develop simulators including accurate descriptions of large deforma- tions. These will be reviewed in Section ‘Refining the appearance: the second generation of surgery simulation’. In Section ‘Surgery simulation by model order reduction’ particu- lar emphasis is made on methods based upon some form of model reduction, which has revealed as one promising alternative for this challenging problem. In Section ‘Adding physiological details: third to fifth generations of surgery simulators’ a review is made of those works aiming at incorporating physiological, microscopic or biochemical details into the simulation. Review Early times of surgery simulation: the first generation The origins of the so-called “deformable models” in the world of computer imaging date back to the mid eighties [19,20]. Basically, they deal with the problem of deforming a solid following the laws of continuum mechanics. However, restrictions posed by real- time requirements did not allow, at that time, to perform these free-form deformations in real-time yet. In 2005, a review by Monserrat et al. [21] divided existing deformable models for surgery simulation into heuristic, hybrid and continuum mechanical approaches, nearly all belonging by that time to the first generation of models. Under the heuristic category, the most relevant approaches are that of spring-masses. Since real-time finite element descriptions of continuum mechanics were not achievable by that time, most approaches tried to simplify (sometimes over-simplify) the equations of linear elasticity [22] by sub- stituting the continuum by a series of springs and masses obeying Newton’s second law. In general, these approaches did not obtain much success, since the deformations they produced were far from realistic. Much more interesting are the approaches based upon continuum mechanics [23]. By that time, most approaches were still based upon linear elasticity. Particularly noteworthy was the attempt to solve linear elasticity by employing boundary element methods, thus avoiding to mesh the interior of the domains [24]. In [25], the authors, on the contrary, computed the equations of motion for a linearly elastic material using FEM and then by using coarse volumetric meshes to do simulations at interactive rates. Simulation of surgical cutting deserves a special comment. From the early times of cut- ting spring-mass systems [26], continuum-based simulations have therefore preferred the use of X-FEM techniques [27,28] or GPU implementations [29]. 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 4 of 18 http://www.amses-journal.com/content/1/1/11 In general, this first generation of surgery simulation lacks today of interest, since some prototypes are now present that incorporate, at least, non-linear strain measures into the simulation. These are reviewed in the next section. Refining the appearance: the second generation of surgery simulation During the fast developments in the first generation of surgery simulators, it was com- monly accepted that existing simulators did not provide with a realistic sensation [4]. This had an origin obviously on the poor constitutive models that real-time constraints allowed to solve by that time. At a visual level it was soon realized that linear elasticity gave raise to a poor perception of deformation, due to the apparent gain in volume that it produces when employed to solve large deformations. In Figure 3 it can be noticed how when solv- ing a beam bending problem undergoing large strains, if linear elasticity is employed, a large gain of volume is perceived. Thus, it is of utmost importance to consider, at least, Saint Venant-Kirchhoff (SVK) models (linear elasticity under large deformation settings) [30,31]. Although SVK models are known to be unstable under compression, for instance [32], they have been considered in a number of recent works as a means to introduce at least a good visual description of large deformations, see [33]. Methods other than finite elements, for instance meshfree methods have also been used for this purpose [34]. Under this same rationale, De and coworkers applied the point collocation-based method of finite spheres (PCMFS) tech- nique to simulate tissue deformations that are geometrically nonlinear [13,35]. They also developed a technique coined as Point-Associated Finite Field (PAFF) [36] that includes a local approximation to the large deformation field near the surgical tool. Even a system based upon the use of neural networks has been presented in [37]. In it, the system is trained with a large set of possible load states in order to achieve real-time performance in the execution loop. A recent survey on the field could also provide the reader with some more insight on the field [38]. An even more surprising approach to the problem of including non-linear strain measures can be found in the work of Delingette [39]. In it, a non-linear (quadratic) spring-mass model is developed so as to be able to reproduce SVK models. The result- ing formulation, in addition, eliminates the typical instabilities in compression of this constitutive model. However, it is known that soft tissues are often composed of collagen fibers and include considerable amount of water [14], and can therefore be modeled as quasi- incompressible. This complex composition cannot be well simulated using only geo- metrically non-linear models and other constitutive laws should be used for a better approximation of the behavior of these tissues. Among them the so-called Holzapfel model [14] and neo-Hookean strain energy functions [32] are very common, see also [15]. Also in [40] a comparison of six different constitutive models was made for the human liver, concluding that an Ogden viscohyperelastic model provided the best result for the registered deformation of a liver on a CT scanner. Furthermore, in [41] a viscoelastic model under large deformations was considered. It can be considered as a particu- lar implementation of a series of Lagrangian codes coined as Total Lagrangian Explicit Dynamics (TLED), developed for that purpose [42]. Basically, it consists of an implemen- tation of a traditional explicit, lumped-mass, finite element code for structural dynamics, in which the constitutive law is evaluated also explicitly. 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 5 of 18 http://www.amses-journal.com/content/1/1/11 Briefly speaking, TLED approaches consider the equations of solid dynamics (including non-linear strain measures) that, after finite element semi-discretisation give rise to a system of equations ¨ ˙ MU + DU + K (U)U = F where M is a lumped mass matrix and D = αM is a dumping matrix, also diagonal. This allows for a computation of the equilibrium equations at an element level, without the Uz -20 -40 -60 -80 -100 -120 -140 -160 -180 -200 -220 -240 -260 Uz -20 -40 -60 -80 -100 -120 -140 -160 -180 -200 -220 -240 -260 Figure 3 Cantilever beam undergoing large deflections modeled using linear elasticity laws (top) and geometrically non-linear strain measures (bottom). Note that the scale in the displacement legend is the same. The apparent gain in volume when linear strain measures are employed is noteworthy. 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 6 of 18 http://www.amses-journal.com/content/1/1/11 need for an assembly procedure. This fact is specially important in a GPU implementa- tion, since each elemental contribution to the equation could be computed separately and independently. In general, these methods do not work well for large integration times, as is well known for classical explicit finite elements. With these techniques, medium-sized models (some 16000 tetrahedra) can be solved at real-time (500 Hz) rates. Of particular interest are also methods based on the use of implicit time integration. Particularly, we highlight works based upon asynchronous preconditioning on CPUs and GPUs [43]. Surgery simulation by model order reduction In view of the above, it seems that the only possibility to date to consider constitutive laws other than SVK was to implement them in Lagrangian finite element codes in which the chosen constitutive law is evaluated explicitly. The main limitation of this approach comes from the fact that the dynamics of the organ is evaluated explicitly during large time intervals (in the order of hours of real user time), therefore lacking of stability in energy conservation, for instance. The possibility of employing implicit energy and momentum conserving integration schemes for the dynamics, along with a consistent evaluation of tangent stiffness matrices in a Newton-Raphson scheme, for instance, seems to be out of reach for the moment. One appealing possibility to solve these limitations comes from model order reduction (MOR) techniques. Roughly speaking, MOR provides the simulation with an “optimal” (in some sense) set of global basis functions, specific to the considered problem, so as to perform simulations with a minimum of degrees of freedom. MOR techniques have been discovered and re-discovered in almost all scientific disciplines in the last century, giving rise to a plethora of different names for essentially similar techniques [44-48]. Thefirst to applyatechnique similartomodel orderreduction in thefield of visu- alization of solid deformations seems to be James [33]. In this work, SVK models were considered, with an heuristic choice of the global (Ritz) basis functions to be employed. In general, in real-time simulation, [8], “... We do notcareabout thetimetaken forone- time pre-calculation such as setting up equations, inverting matrices, etc.”. This means that the possibility exists of employing as much time as needed in off-line computa- tions whose results are to be used on-line. This paves the way for a consistent off-line determination of the optimal set of basis functions for a given problem. Proper orthogonal decomposition In Proper Orthogonal Decomposition (POD) techniques, it is assumed that the evolu- tion of the unknown field (here assumed scalar, for simplicity of the exposition), u(x, t) is known. In practical applications (assume that some numerical simulations have been performed off-line), this field is expressed in a discrete form which is known at the nodes m m of a spatial mesh and for some instants of time t . Thus, we consider that u (x , t ) = m m m m u (x ) ≡ u (t = m × t) are known. We can also write u for the vector contain- ing the nodal degrees of freedom at time t . The main idea of the Karhunen-Loève (KL) decomposition is to obtain the most typical or characteristic structure φ(x) among these u (x), ∀m. This is equivalent to obtain a function that maximizes the functional α: 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 7 of 18 http://www.amses-journal.com/content/1/1/11 m=M i=N φ (x ) u (x ) i i m=1 i=1 α = ,(1) i=N (φ (x )) i=1 where N represents the number of nodes of the complete model and M the number of computed time steps. The maximization leads to: ⎡ ⎛ ⎞ ⎤ j=N m=M i=N i=N m m ⎣ ˜ ⎝ ⎠ ⎦ ˜ ˜ φ (x ) u (x ) φ x u x = α φ (x ) φ (x ) ; ∀φ, i i j j i i m=1 i=1 j=1 i=1 which can be rewritten in the form ⎧ ⎫ j=N i=N m=M i=N ⎨ ⎬ m m ˜ ˜ ˜ u (x ) u x φ x φ (x ) = α φ (x ) φ (x ) ; ∀φ.(2) i j j i i i ⎩ ⎭ i=1 j=1 m=1 i=1 Defining the vector φ such that its i-th component is φ(x ), Eq. (2) takes the following matrix form T T ˜ ˜ ˜ φ c φ = αφ φ; ∀φ ⇒ c φ = αφ,(3) where the two-point correlation matrix is given by m=M m=M m m m m c = u (x ) u x ⇔ c = u u , ij i j m=1 m=1 which is symmetric and positive definite. If we define the matrix Q containing the discrete field history: ⎛ ⎞ 1 2 M u u ··· u 1 1 1 ⎜ ⎟ 1 2 M u u ··· u ⎜ ⎟ 2 2 2 ⎜ ⎟ Q = ,(4) . . . ⎜ . ⎟ . . . . ⎝ . ⎠ . . . 1 2 M u u ··· u N N N then it is easy to verify that the matrix c in Eq. (3) results in c = QQ . Assume now that some direct simulations have been carried out (off-line, this process may take as much time as needed). From them, we can determine u , ∀i ∈ [1, ··· , N] and ∀m ∈ [1, ··· , M], and from these solutions the n eigenvectors related to the n-highest eigenvalues that are expected to contain the most important information (at least, sta- tistically speaking) about the problem solution. For this purpose the eigenvalue problem defined by Eq. (3) is solved and all the eigenvalues φ belonging to some pre-defined inter- val are stored. In practice n is much lower than N, and this constitutes the main advantage of the technique. POD techniques use therefore these n eigenfunctions φ as Ritz-like (thus, globally supported) basis functions for approximating the solution of a problem slightly different to the one that has served to define u . For this purpose we need to define the matrix A = [φ ··· φ ] 1 n ⎛ ⎞ φ (x ) φ (x ) ··· φ (x ) 1 1 2 1 n 1 ⎜ ⎟ φ (x ) φ (x ) ··· φ (x ) ⎜ ⎟ 1 2 2 2 n 2 ⎜ ⎟ A = . . . . ⎜ . ⎟ . . . . ⎝ . ⎠ . . . φ (x ) φ (x ) ··· φ (x ) 1 N 2 N n N 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 8 of 18 http://www.amses-journal.com/content/1/1/11 Consider now the linear system of equations coming from the discretization of a generic problem, in the form: m m−1 Ku = F,(5) where the superscript refers to the time step. By assuming that the unknown vector contains the nodal degrees of freedom, it can be expressed as: i=n m m m u = ζ φ = A ζ , i=1 from which Eq. (5) results m m−1 m m−1 Ku = F ⇒ KA ζ = H , and by multiplying both terms by A T m T m−1 A KA ζ = A F , is obtained, which proves that the final system of equations is of low order, i.e. the dimension of A GA is n × n,with n N. It is well known that reduced-order models allow for larger time steps in structural dynamics simulations [49]. This is the approach followed in [50,51], where POD reduced models were employed trying to enlarge admissible time steps in explicit dynamics solved by finite element methods. In [52] a local/global approach to the problem of model reduc- tion, where a reduced model is employed far from the “conflict” zones, i.e., those with crack sharp gradients, and a more classical finite element implementation is done where the crack influence is bigger. A similar approach (a coupling between POD and FEM) was previously accomplished in [53]. The reader interested in dynamical aspects of model order reduction could also consult an excellent paper by Matthies [54]. In sum, from the results of complete-order models solved previously and off-line, it is possible to obtain a basis which is optimal (in the sense that it minimizes the functional in Eq. (1)) to be used in problems slightly different to the original ones. This is essentially the approach followed in [55-57]. As can be noticed from Figure 4, where an example of palpation of a human cornea is shown, there are clear differences between both results, even if the reduced model, that employs only six degrees of freedom, could be refined. The origin of these discrepancies comes from the fact that a non-linear, hyperelastic, constitutive equation was used for the corneal tissue [15], while the strong requirements Figure 4 Vertical displacement field for a particular load position in the palpation of a cornea. Complete model (left) vs. reduced model (six d.o.f., right). 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 9 of 18 http://www.amses-journal.com/content/1/1/11 of real-time feedback rates prevent from computing a consistently linearized tangent stiff- ness matrix. Thus, the simulation is actually linear, and the tangent stiffness matrix at the origin is actually being employed at any time in the simulation. Other source of problems, which is classical in POD-based model order reduction, is the choice of the complete models from which to extract the snapshots matrix, Eq. (4). More precisely, it is of utmost importance how to interpolate between reduced models when they are obtained after complete models reproducing conditions (load positions, for instance), which are not exactly reproduced in the on-line phase of the method. In [58] it was pointed out that the set of empirical eigenfunctions given by Eq. (3) for a given model forms the so-called Grassman manifold G(n, N ). Therefore, in order to interpolate the set of basis functions, that consequently do not form a vector space, we must move to the tangent plane at a point of the manifold, which is a “flat” space, interpolate there, and project back to the manifold. A less rigorous, but much simpler, method to interpolate among previously computed reduced models was established in [59] and coined as POD with interpolation (PODI). Although in the standard PODI technique the POD procedure is applied to the complete set of snapshots (for different load positions, say) of the system to obtain an orthonormal basis A = [φ ··· φ ], in [55] POD was just applied to each complete model (i.e., to each 1 n load position). Thus, an orthonormal basis for each system’s parameter value is obtained. Basis were then interpolated for intermediate positions of the load. Although it is clear that the interpolation of orthogonal sets of functions does not yield, in general, new orthogonal basis functions, the technique works well if the “distance” of reduced models in the Grassman manifold is not too large (an aspect always very difficult to determine and covered very elegantly in [60]). Further developments in the real-time simulation based upon POD techniques include the employ of asymptotic numerical methods (ANM) [61-64] in combination with POD [65-67]. By means of an asymptotic expansion of the variables of interest, POD-ANM techniques allow to solve non-linear problems without the burden associated to the update and inversion of tangent stiffness matrices. Instead, a series of linear problems, all with the same tangent operator, are solved, one for each level of expansion. These meth- ods will be detailed in Section ‘Methods based on proper generalized decomposition’ below. Another way to deal with non-linear problems in a reduced order context is to consider the Empirical Interpolation Method (EIM) [68] or its discrete counterpart, the DEIM [69]. Some problems, in general, presents special difficulties for their reduction. These are problems in which advancing fronts are present, in general, that lead to hardly separable models that need for many modes (or degrees of freedom). A very interesting prob- lem that poses difficulties in this sense for simulation in general, but very specially for real-time simulation, is that of electrophysiology (see [70,71] and references therein, for instance). The main cause of these difficulties comes from the fact that electrophysiol- ogy is governed by anisotropic reaction-diffusion equations with a very rapidly varying reaction term. Accurate schemes need for very fine meshes and extremely small time step size. The fact that the solution of the problem is roughly composed by an advancing front (very similar to a soliton wave) makes this problem specially difficult for model reduc- tion methods, that do not capture properly the solution. In this framework, the work by Gerbeau and coworkers [72-74] is noteworthy. They use mainly POD techniques to 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 10 of 18 http://www.amses-journal.com/content/1/1/11 that end, although very recently they employed the concept of Lax pairs to avoid the mentioned problems, characteristic of electrophysiology [74]. Methods based on proper generalized decomposition Model reduction methods based on POD or related techniques are known as “a poste- riori” methods, since they need for a number of complete model simulations in order to obtain the snapshots. It is possible, however, to develop “a priori” methods in which there is no need for a complete-model solution. Proper Generalized Decomposition (PGD) methods belong to this class of methods. Although initially designed for a space-time decomposition of the solution within the LATIN method [75] or, independently, for high dimensional problems in statistical mechanics [76,77], PGD methods provide a very appealing alternative to solve in real time complex computational problems. In fact, PGD approaches to the problem at hand constitute a change of paradigm, since they are based upon the off-line computation of a general (high dimensional) solution to the problem, that is after evaluated (rather than simulated) in real time in the on-line phase [78,79]. So to speak, it generates a sort of computational vademecum [80] that pro- vides the solution for any possible situation during the on-line phase of the simulation. Of course, obtaining such a general solution and storing it efficiently is not an easy task. PGD methods tackle these difficulties by means of the following ingredients: 1. The vademecum solution is in fact a high dimensional solution to the problem, depending on many different parameters such as initial or boundary conditions, material parameters, etc. To deal with the so-called curse of dimensionality (the number of degrees of freedom of the solution grows exponentially with the number of dimensions in traditional mesh-based methods), PGD proposes a separated representation. Consider, for the sake of simplicity in the exposition, that the solution of an organ palpation depends solely on the physical coordinates and the load position: n k k u (x, s) = X (x) · Y (s),(6) j j j k=1 where u represents the j-th component of the displacement field, depending on the physical position x and the mentioned position of the load exerted by the scalpel, s. Dependence of the solution on further parameters can be treated under this same rationale without difficulties. Problems defined in spaces up to 100 dimensions have been solved successfully in the literature, see [78] and references therein. 2. Superscript n in Eq. (6) indicates that PGD proceeds actually in an iterative way. Improvements in the approximation are obtained by applying a greedy algorithm, i.e., n+1 u (x, s) = u (x, s) + R (x) · S (s), j j where R(x) and S(s) are the sought functions that improve the approximation. To obtain these functions, consider first the admissible variation of the displacement, given by ∗ ∗ ∗ u (x, s) = R (x) · S (s) + R (x) · S (s).(7) j j j j j 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 11 of 18 http://www.amses-journal.com/content/1/1/11 Several options are at hand so as to determine the new pair of functions R and S. The most frequently used, due to both its ease of implementation and good convergence properties, in general, is a fixed-point algorithm in which functions R and S are sought iteratively so as to minimize the residual. See the extensive bibliography on PGD for more details. 3. PGD methods are in fact a model order reduction technique, since the order of the approximation, n, can be fixed so as to have a minimum of degrees of freedom in the model. To this end, error indicators are indispensable [81,82]. Functions X and Y in Eq. 6 can actually be seen as the PGD counterpart of the POD modes φ . However, as opposed to POD, PGD does not always provide with optimal modes. Analytical proofs of optimality exist for elliptic problems with symmetric operators [83]. Another difficult problem, common to virtually all model reduction methods, is that of the non-linearity of the problem. Since most living tissues are inherently non-linear, these problems must be solved within the PGD framework. In [84,85] an explicit approach and one based upon ANM are presented, respectively. Both approaches can lead to similar levels of accuracy, although obtained with very different modes, see Figure 5. Y Y Z Z X X Uz Uz -0.01 0.03 -0.02 0.025 -0.03 0.02 -0.04 0.015 -0.05 0.01 -0.06 0.005 -0.07 0 -0.08 -0.005 -0.09 -0.01 -0.1 -0.015 -0.11 -0.02 -0.12 -0.025 -0.03 -0.035 -0.04 (a) (b) Y Y Z X Z X Uy Uy 0.015 0.015 0.01 0.01 0.005 0.005 0 0 -0.005 -0.005 -0.01 -0.01 -0.015 -0.015 -0.02 -0.02 -0.025 -0.025 -0.03 -0.03 -0.035 -0.035 -0.04 -0.04 -0.045 -0.045 -0.05 -0.05 -0.055 -0.055 -0.06 -0.06 -0.065 -0.065 (c) (d) 1,2 Figure 5 Two first modes X (x) in the simulation of liver palpation. Top: explicit version in [84]. (a) first mode, (b) second mode. Bottom: PGD-ANM approach in [85]. (c) first mode, (d) second mode. 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 12 of 18 http://www.amses-journal.com/content/1/1/11 Since, in principle, it is only necessary to store these modes in the form of vectors (and, in thecaseoffunctions Y (s) they are only defined on the boundary of the organ = ∂), the storage in memory of these results is very efficient. As an example, this method allows for an interative implementation on deployed, handheld devices such as smartphones and tables (see Figure 6) or web-based demonstrators (see Figure 7). Under the model reduction rationale, aspects such as cutting are much less investigated, but in [67], for instance, an approach based upon a combination of X-FEM and POD is developed. Regarding PGD applications, a similar approach could be envisaged, based on a PGD-XFEM coupling strategy [86]. Adding physiological details: third to fifth generations of surgery simulators Understanding and quantitatively simulating physiology of the human has been a chal- lenge since some fifteen years ago. The Physiome project [87-89] is a clear paradigm of this, opening the era of computational physiology. But an effective, real-time simulation of physiology is nowadays even further than a successful, purely mechanical, real-time simulation at a continuum level of living tissues, as discussed in Section ‘Refining the appearance: the second generation of surgery simulation’ before. While introducing physiology in the simulation is not strictly necessary for many appli- cations (most minimally invasive surgery training procedures, for instance), simulation based medical or surgical planning [90] very often needs it. For instance, in [91] a novel haptic simulation system is developed that includes not only force but thermal feedback, useful for the location of tumors, for instance. In other problems, including physiology in the simulation is simply unavoidable. For instance, predicting the outcome of vein graft surgery is one of these problems that com- bines the need for simulation at a macroscopic level, but also at a gene regulatory network level [1]. It is hypothesized that blood shearing forces modulates a specific gene regulatory network determining the adaptive response of the vein wall. But simulation at a gene Figure 6 Implementation of the PGD approach to the liver palpation problem on a Motorola Xoom tablet. It provides continuous visual feedback without problems. 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 13 of 18 http://www.amses-journal.com/content/1/1/11 Figure 7 Implementation of the PGD approach to the liver palpation problem on a web page using javascript. regulatory network level is a challenge in itself. This is so since at this level a stochastic description of the reactions occurring within the cell must be done, since it is nonsense to speak of concentration at this level, where only few specimens of each component are present [9,92,93]. It is commonly accepted that under some weak hypothesis (well stirred mixture, fixed volume and temperature), the system can be considered as Markovian, and can be consequently modeled by the so-called Chemical Master Equation (CME), [94], which is in fact no more than an ordinary differential equations stating the conservation of the probability density function P in time: ∂P(z, t|z , t ) 0 0 = a z − v P z − v , t|z , t − a (z)P (z, t|z , t ),(8) j j j 0 0 j 0 0 ∂t where P (z, t|z , t ) represents the probability of being at a state in which there are a num- 0 0 ber of molecules of each species stored in the vector z at time t when we started from a state z at time t . a represents the propensity (i.e., the probability) of reaction j to occur, 0 0 j while v represents the change in the number of molecules of each species if reaction j takes place. This change is given, of course, by the stoichiometry of the reaction at hand. This set of equations are defined in a state space which possess as many dimensions as the number of different species involved in the regulatory network. Under this challenging framework, if we consider N different species, present at a number n of copies, the num- berofdifferent possiblestatesofthe system is n .Thisnumbercan take theastronomical 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 14 of 18 http://www.amses-journal.com/content/1/1/11 value of 10 if we consider some types of proteins, for instance [94].This phenomenon is known as the curse of dimensionality and has been already mentioned in this survey, see Section ‘Refining the appearance: the second generation of surgery simulation’. For instance, Nobel prize winner R. B. Laughlin said, when talking about this problem [95] that “No computer existing, or that will ever exist, can break this barrier because it is a catastropheofdimension”. Monte Carlo-like algorithms (the so-called stochastic simulation algorithm, SSA [94,96,97]) are often the preferred approach to these problems. But Monte Carlo tech- niques need for as many as possible individual realizations of the problem that compro- mise its simple application in inverse identification, leading to excessive time consuming simulations, together with great variance in the results. In [98], however, a PGD technique is used to avoid precisely the curse of dimension- ality. Again, an off-line/on-line strategy is developed that allows to obtain vademecum solutions easy to evaluate on-line under very astringent real-time restrictions, for exam- ple when considering the propensities as extra-coordinates. As commented before, PGD approximates the variable of interest, the probability P,asafinite sumofseparable functions: j j j j P (z, t) = F (z ) · F (z ) · ... · F (z ) · F (t), 1 2 D 1 2 D j=1 where, as mentioned before, the variables z represent the number of molecules of species i present at a given time instant. Note the presence of a function depending solely on time, F (t). This means that the algorithm is not incremental. Instead, it solves for the whole time history of the chemical species at each iteration of the method, so as to be available to be evaluated at any point in the multiscale simulation pipeline. With this technique, in [98] cascades of up to 20 terms were successfully simulated giving results in good accor- dance with the best ones available in the literature. This approach allows to compute the steady state of a system, impossible when considering stochastic approaches. Generally speaking, adding physiological details to real-time simulations for surgery planning is still out of reach in many cases. Additional efforts of research are needed to successfully incorporate many aspects such as bleeding, microscopical details such as glandular anatomy or, finally, biochemical aspects such as those just commented. Conclusions Real-time simulation in the field of computational surgery has evolved a lot in the last fifteen years. Although not exhaustive, in this review paper we have presented a sur- vey of some of the most prominent techniques nowadays. Particular attention has been paid to techniques based on model order reduction, a framework that has somewhat revolutionize the field and that seems the only way to incorporate to the field detailed physiological descriptions, with microscopic details and micro-macro approaches to the always challenging work of simulating the human body. Other aspects, such as patient-specific modeling, deserve additional attention. Due to their inherent complexity and the infancy of their state of development, they have not been covered in this review. Competing interests The authors declare that they have no competing interests. 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 15 of 18 http://www.amses-journal.com/content/1/1/11 Authors’ contributions Both authors, EC and FC, participated equally in the manuscript writing. Both read and approved the final manuscript. Acknowledgments This work has been partially supported by the Spanish Ministry of Economy and Competitiveness, through grant number CICYT-DPI2011-27778-C02-01. Author details 1 2 Aragón Institute of Engineering Research, Universidad de Zaragoza, Zaragoza, Spain. Ecole Centrale de Nantes, Nantes, France. Institut Universitaire de France, Paris, France. Received: 26 August 2013 Accepted: 4 April 2014 Published: 29 April 2014 References 1. Garbey M, Bass B, Berceli S (2012) Multiscale mechanobiology modeling for surgery assessment. Acta Mechanica Sinica 28:1186–1202. http://dx.doi.org/10.1007/s10409-012-0133-4, doi:10.1007/s10409-012-0133-4 2. Champion HR, Gallagher AG (2003) Surgical simulation: a good idea whose time has come. Br J Surg 90(7):767–768. doi:10.1002/bjs.4187, http://dx.doi.org/10.1002/bjs.4187 3. Winslow RL, Trayanova N, Geman D, Miller MI (2012) Computational medicine: translating models to clinical care. Sci Transl Med 4(158). doi:10.1126/scitranslmed.3003528 4. Delingette H, Ayache N (2004) Soft tissue modeling for surgery simulation. In: Ayache N (ed) Computational models for the human body, handbook of numerical analysis (Ph. Ciarlet, Ed.) Elsevier, Idem, pp 453–550 5. Satava RM (1996) Medical virtual reality: the current status of the future. In: Sieburg H, Weghorst S, Morgan K (eds) Healthcare in the information age. IOS Press 6. Delingette H, Ayache N (2005) Hepatic surgery simulation. Comm ACM 48:31–36. doi:http://doi.acm.org/10.1145/1042091.1042116, http://doi.acm.org/10.1145/1042091.1042116 7. Cotin S, Delingette H, Ayache N (1999) Real-time elastic deformations of soft tissues for surgery simulation. In: Hagen H (ed) IEEE transactions on visualization and computer graphics, vol 5 (1). IEEE Computer Society, pp 62–73. citeseer. ist.psu.edu/cotin98realtime.html 8. Bro-Nielsen M, Cotin S (1996) Real-time volumetric deformable models for surgery simulation using finite elements and condensation. Comput Graph Forum 15(3):57–66 9. Hasty J, McMillen D, Isaacs F, Collins JJ (2001) Computational studies of gene regulatory networks: in numero molecular biology. Nat Rev Genet 2:268–279 10. Hegland M, Burden C, Santoso L, MacNamara S, Boothm H (2007) A solver for the stochastic master equation applied to gene regulatory networks. J Comput Appl Math 205:708–724 11. Sasai M, Wolynes PG (2003) Stochastic gene expression as a many-body problem. Proc Natl Acad Sci 100(5): 2374–2379 12. Kim J, De S, Srinivasan MA (2002) Proceedings 10th symposium on haptic interfaces for virtual environment and teleoperator systems. haptics 2002 In: Haptic Interfaces for Virtual Environment and Teleoperator Systems, 2002. HAPTICS 2002. Proceedings. 10th Symposium on, pp 51–57. doi:10.1109/HAPTIC.2002.998933 13. Basdogan C, De S, Kim J, Muniyandi M, Kim H, Srinivasan M, 2004 Haptics in minimally invasive surgical simulation and training. IEEE Comput Graph Appl Mag 24(2):56–64. doi:10.1109/MCG.2004.1274062 14. Holzapfel GA, Gasser TC (2000) A new constitutive framework for arterial wall mechanics and a comparative study of material models. J Elasticity 61:1–48 15. Alastrue V, Calvo B, Pena E, Doblare M (2006) Biomechanical modeling of refractive corneal surgery. J Biomecha Eng-Trans ASME 128:150–160 16. Berceli SA, Tran-Son-Tay R, Garbey M, Jiang Z (2009) Hemodynamically driven vein graft remodeling: A systems biology approach. Vascular 17(suppl 1):S2–S9. doi:10.2310/6670.2008.00083, http://vas.sagepub.com/content/17/ suppl_1/S2.abstract, http://vas.sagepub.com/content/17/suppl_1/S2.full.pdf+html 17. Vilanova G, Colominas I, Gomez H (2013) Capillary networks in tumor angiogenesis: From discrete endothelial cells to phase-field averaged descriptions via isogeometric analysis. Int J Numer Methods Biomed Eng. doi:10.1002/cnm.2552, http://dx.doi.org/10.1002/cnm.2552 18. Azuma R (1997) A survey of augmented reality. Presence-Teleoperators Virtual Environ 6(4):355–385 19. Sederberg TW, Parry SR (1986) Free-form deformation of solid geometric models In: Proceedings of the 13th annual conference on Computer graphics and interactive techniques, SIGGRAPH ‘86. New York, pp 151–160. doi:10.1145/15922.15903, http://doi.acm.org/10.1145/15922.15903 20. Terzopoulos D, Fleischer K (1988) Deformable models. Vis Comput 4(6):306–331. doi:10.1007/BF01908877, http://dx. doi.org/10.1007/BF01908877 21. Meier U, Lopez O, Monserrat C, Juan M, Alcaniz M (2005) Real-time deformable models for surgery simulation: a survey. Comput Meth Programs Biomed 77(3):183–197. doi:10.1016/j.cmpb.2004.11.002, http://www.sciencedirect. com/science/article/pii/S0169260704002093 22. Lee Y, Terzopoulos D, Waters K (1995) Realistic modeling for facial animation In: Proceedings of the 22nd annual conference on computer graphics and interactive techniques, SIGGRAPH ‘95. ACM, New York, pp 55–62. doi:10.1145/218380.218407, http://doi.acm.org/10.1145/218380.218407 23. Famaey N, Stolen J (2008) Soft tissue modeling for applications in virtual surgery and surgical robotics. Comput Methods Biomech Biomed Engin 11(4):351–366 24. Monserrat C, Meier U, Alcaniz M, Chinesta F, Juan M (2001) A new approach for the real-time simulation of tissue deformations in surgery simulation. Comput Methods Programs Biomed 64(2):77–85. doi:10.1016/S0169-2607(00)00093-6, http://www.sciencedirect.com/science/article/pii/S0169260700000936 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 16 of 18 http://www.amses-journal.com/content/1/1/11 25. Capell S, Green S, Curless B, Duchamp T, Popovic Z (2002) Interactive skeleton-driven dynamic deformations In: ACM SIGGRAPH. ACM, New York 26. Cotin S, Delingette H, Ayache N (2000) A hybrid elastic model for real-time cutting, deformations, and force feedback for surgery training and simulation. Vis Comput 16(8):437–452. doi:10.1007/PL00007215, http://dx.doi.org/ 10.1007/PL00007215 27. Vigneron L, Verly J, Warfield S On extended finite element method (xfem) for modelling of organ deformations associated with surgical cuts. In: Cotin S, Metaxas D (eds) Medical Simulation, Lecture Notes in Computer Science, vol 3078. Springer Berlin Heidelberg, pp 134–143. doi:10.1007/978-3-540-25968-8_15, http://dx.doi.org/10.1007/ 978-3-540-25968-8_15 28. Jeˇrábková L, Kuhlen T (2009) Stable cutting of deformable objects in virtual environments using xfem. IEEE Comput Graph Appl 29(2):61–71. doi:10.1109/MCG.2009.32, http://dx.doi.org/10.1109/MCG.2009.32 29. Courtecuisse H, Jung H, Allard J, Duriez C, Lee DY, Cotin S (2010) GPU-based real-time soft tissue deformation with cutting and haptic feedback. Prog Biophys Mol Biol 103(2–3):159–168. doi: http://dx.doi.org/10.1016/j.pbiomolbio. 2010.09.016, http://www.sciencedirect.com/science/article/pii/S0079610710000842, special Issue on Biomechanical Modelling of Soft Tissue Motion 30. Zhong H, Peters T (2007) A real time hyperelastic tissue model. Comput Methods Biomech Biomed Engin 10(3): 185–193. doi:10.1080/10255840701292732, http://www.tandfonline.com/doi/abs/10.1080/10255840701292732, pMID: 17558647, http://www.tandfonline.com/doi/pdf/10.1080/10255840701292732 31. Zhong H, Wachowiak MP, Peters TM (2005) A real time finite element based tissue simulation method incorporating nonlinear elastic behavior. Comput Methods Biomech Biomed Engin 8(3):177–189. doi:10.1080/10255840500295852, http://www.tandfonline.com/doi/abs/10.1080/10255840500295852, pMID: 16214712, http://www.tandfonline.com/doi/pdf/10.1080/10255840500295852 32. Bonet J, Wood RD (2008) Nonlinear continuum mechanics for finite element analysis. Cambridge University Press, Cambridge 33. Barbicˇ J, James DL (2005) Real-time subspace integration for St. Venant-Kirchhoff deformable models. ACM Trans Graph 24(3):982–990 34. Horton A, Wittek A, Joldes GR, Miller K (2010) A meshless total lagrangian explicit dynamics algorithm for surgical simulation. Int J Numer Methods Biomed Eng 26(8):977–998. doi:10.1002/cnm.1374, http://dx.doi.org/10.1002/cnm. 35. Lim YJ, De S (2007) Real time simulation of nonlinear tissue response in virtual surgery using the point collocation-based method of finite spheres. Comput Methods Appl Mech Eng 196:3011–3024 36. De S, Lim YJ, Manivannan M, Srinivasan MA (2006) Physically realistic virtual surgery using the point-associated finite field (paff) approach. Presence: Teleoperators Virtual Environ 15(3):294–308. doi:10.1162/pres.15.3.294, http://dx.doi. org/10.1162/pres.15.3.294 37. Deo D, De S (2009) Phyness: a physics-driven neural networks-based surgery simulation system with force feedback In: Third joint eurohaptics conference and symposium on haptic interfaces for virtual environment and teleoperator systems, Salt Lake City, UT, USA 38. Misra S, Ramesh KT, Okamura AM (2008) Modeling of tool-tissue interactions for computer-based surgical simulation: a literature review. Presence: Teleoper Virtual Environ 17(5):463–491. doi:10.1162/pres.17.5.463, http://dx. doi.org/10.1162/pres.17.5.463 39. Delingette H (2008) Biquadratic and quadratic springs for modeling st venant kirchhoff materials In: Proceedings of the 4th international symposium on biomedical simulation, ISBMS ‘08. Springer-Verlag, Berlin, Heidelberg, pp 40–48. doi:10.1007/978-3-540-70521-5_5, http://dx.doi.org/10.1007/978-3-540-70521-5_5 40. Martinez-Martinez F, Ruperez M, Martin-Guerrero J, Monserrat C, Lago M, Pareja E, Brugger S, Lopez-Andujar R (2013) Estimation of the elastic parameters of human liver biomechanical models by means of medical images and evolutionary computation. Comput Methods Programs Biomed 111(3):537–549. doi:http://dx.doi.org/10.1016/j.cmpb.2013.05.005, http://www.sciencedirect.com/science/article/pii/ S0169260713001491 41. Comas O, Taylor ZA, Allard J, Ourselin S, Cotin S, Passenger J (2008) Efficient nonlinear FEM for soft tissue modelling and its GPU implementation within the open source framework SOFA In: Proceedings of the 4th international symposium on biomedical simulation, ISBMS ‘08. Springer-Verlag, Berlin, Heidelberg, pp 28–39. doi:10.1007/978-3-540-70521-5_4, http://dx.doi.org/10.1007/978-3-540-70521-5_4 42. Miller K, Joldes G, Lance D, Wittek A (2007) Total lagrangian explicit dynamics finite element algorithm for computing soft tissue deformation. Comm Numer Meth Eng 23(2):121–134. doi:10.1002/cnm.887, http://dx.doi.org/ 10.1002/cnm.887 43. Courtecuisse H, Allard J, Kerfriden P, Bordas SP, Cotin S, Duriez C (2014) Real-time simulation of contact and cutting of heterogeneous soft-tissues. Med Image Anal 18(2):394–410. doi:http://dx.doi.org/10.1016/j.media.2013.11.001, http://www.sciencedirect.com/science/article/pii/S1361841513001692 44. Karhunen KK (1947) Über lineare Methoden in der Wahrscheinlichkeitsrechnung. Ann Acad Sci Fennicae Ser A I Math-Phys 37:1–79 45. Loève MM (1963) Probability theory. The University Series in Higher Mathematics. 3rd ed. Van Nostrand, Princeton 46. Lorenz EN (1956). MIT, Departement of Meteorology, Scientific Report Number 1, Statistical Forecasting Project 47. Park HM, Cho DH (1996) The use of the Karhunen-Loève decomposition for the modeling of distributed parameter systems. Chem Eng Sci 51(1):81–98 48. Sirovich L (1987) Turbulence and the dynamics of coherent structures part I: coherent structures. Q Appl Math XLV: 561–571 49. Krysl P, Lall S, Marsden J (2001) Dimensional model reduction in non-linear finite element dynamics of solids and structures. Int J Numer Meth in Engng 51:479–504 50. Taylor Z, Crozier S, Ourselin S (2011) A reduced order explicit dynamic finite element algorithm for surgical simulation. IEEE Trans Med Imaging 30(9):1713–1721. doi:10.1109/TMI.2011.2143723 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 17 of 18 http://www.amses-journal.com/content/1/1/11 51. Taylor Z, Ourselin S, Crozier S (2010) A reduced order finite element algorithm for surgical simulation In: Engineering in Medicine and Biology Society (EMBC). 2010 Annual International Conference of the IEEE, pp 239–242. doi:10.1109/IEMBS.2010.5627720 52. Kerfriden P, Passieux JC, Bordas SPA (2012) Local/global model order reduction strategy for the simulation of quasi-brittle fracture. Int J Numer Methods Eng 89(2):154–179. doi:10.1002/nme.3234, http://dx.doi.org/10.1002/ nme.3234 53. Ammar A, Pruliere E, Ferec J, Chinesta F, Cueto E (2009) Coupling finite elements and reduced approximation bases. Eur J Comput Mech 18(5–6):445–463 54. Meyer M, Matthies HG (2003) Efficient model reduction in non-linear dynamics using the Karhunen-Loève expansion and dual-weighted-residual methods. Comput Mech 31(1–2):179–191. doi:10.1007/s00466-002-0404-1, http://dx.doi.org/10.1007/s00466-002-0404-1 55. Niroomandi S, Alfaro I, Cueto E, Chinesta F (2008) Real-time deformable models of non-linear tissues by model reduction techniques. Comput Methods Programs Biomed 91(3):223–231. doi:10.1016/j.cmpb. 2008.04.008, http://www.sciencedirect.com/science/article/B6T5J-4SNPPVY-2/2/ 8a417e7f1371768b4c928d1f12fc7a0f 56. Dogan F, Celebi MS (2011) Real-time deformation simulation of non-linear viscoelastic soft tissues. Simulation 87(3): 179–187 57. Radermacher A, Reese S (2013) Proper orthogonal decomposition-based model reduction for nonlinear biomechanical analysis. Int J Mater Eng Innovat 4(4):149–165. doi:10.1504/IJMATEI.2013.054393 58. Amsallem D, Farhat C (2008) An interpolation method for adapting reduced-order models and application to aeroelasticity. AIAA J 46:1803–1813 59. Ly HV, Tran HT (2005) Modeling and control of physical processes using proper orthogonal decomposition. Math Comput Model 33:223–236 60. Son NT (2013) A real time procedure for affinely dependent parametric model order reduction using interpolation on grassmann manifolds. Int J Numer Meth Eng 93(8):818–833. doi:10.1002/nme.4408, http://dx.doi.org/10.1002/ nme.4408 61. Cochelin B, Damil N, Potier-Ferry M (1994) Asymptotic-numerical methods and Padéapproximants for non-linear elastic structures. Int J Numer Methods Eng 37:1187–1213 62. Abichou H, Zahrouni H, Potier-Ferry M (2002) Asymptotic numerical method for problems coupling several nonlinearities. Comput Meth Appl Mech Eng 191(51–52):5795–5810 63. Cao HL, Potier-Ferry M (1999) An improved iterative method for large strain viscoplastic problems. Int J Numer Methods Eng 44:155–176 64. Cochelin B, Damil N, Potier-Ferry M (1994) The asymptotic numerical method: an efficient perturbation technique for nonlinear structural mechanics. Revue Europeenne des Elements Finis 3:281–297 65. Niroomandi S, Alfaro I, Cueto E, Chinesta F (2010) Model order reduction for hyperelastic materials. J Numer Meth Eng 81(9):1180–1206. doi:10.1002/nme.2733, http://dx.doi.org/10.1002/nme.2733 66. Niroomandi S, Alfaro I, Cueto E, Chinesta F (2012) Accounting for large deformations in real-time simulations of soft tissues based on reduced-order models. Comput Methods Programs Biomed 105(1):1–12. doi:DOI: 10.1016/j.cmpb.2010.06.012, http://www.sciencedirect.com/science/article/B6T5J-50VGHDD-1/2/ 1201566766c0d280af9195bf07bfaf91 67. Niroomandi S, Alfaro I, Gonzalez D, Cueto E, Chinesta F (2012) Real-time simulation of surgery by reduced-order modeling and x-fem techniques. Int J Numer Methods Biomed Eng 28(5):574–588. doi:10.1002/cnm.1491, http://dx. doi.org/10.1002/cnm.1491 68. Barrault M, Maday Y, Nguyen N, Patera A (2004) An ‘empirical interpolation’ method: application to efficient reduced-basis discretization of partial differential equations. Comptes Rendus Mathematique 339(9):667–672. doi:10.1016/j.crma.2004.08.006 69. Chaturantabut S, Sorensen DC (2010) Nonlinear model reduction via discrete empirical interpolation. SIAM J Sci Comput 32:2737–2764. doi:10.1137/090766498, http://dx.doi.org/10.1137/090766498 70. Boulakia M, Cazeau S, Fernandez M, Gerbeau JF, Zemzemi N (2010) Mathematical modeling of electrocardiograms: A numerical study. Ann Biomed Eng 38(3):1071–1097. doi:10.1007/s10439-009-9873-0, http://dx.doi.org/10.1007/ s10439-009-9873-0 71. Heidenreich EA, Ferrero JM, Doblare M, Rodriguez JF (2010) Adaptive macro finite elements for the numerical solution of monodomain equations in cardiac electrophysiology. Ann Biomed Eng 38(7):2331–2345. doi:10.1007/s10439-010-9997-2, http://dx.doi.org/10.1007/s10439-010-9997-2 72. Boulakia M, Schenone E, Gerbeau JF (2012) Reduced-order modeling for cardiac electrophysiology. application to parameter identification. Int J Numer Methods Biomed Eng 28(6–7):727–744. doi:10.1002/cnm.2465, http://dx.doi. org/10.1002/cnm.2465 73. McLeod K, Caiazzo A, Fernandez MA, Mansi T, Vignon-Clementel IE, Sermesant M, Pennec X, Boudjemline Y, Gerbeau JF (2010) Atlas-based reduced models of blood flows for fast patient-specific simulations. In: Camara O, Pop M, Rhode K, Sermesant M, Smith N, Young A (eds) Statistical Atlases and Computational Models of the Heart, Lecture Notes in Computer Science, vol 6364. Springer Berlin Heidelberg, pp 95–104. doi:10.1007/978-3-642-15835-3_10, http://dx.doi.org/10.1007/978-3-642-15835-3_10 74. Gerbeau JF, Lombardi D (2012) Reduced-order modeling based on approximated lax pairs. J Comput Phys 265(2014):246–269 75. Ladeveze P (1999) Nonlinear computational structural mechanics. Springer, New York 76. Ammar A, Mokdad B, Chinesta F, Keunings R (2006) A new family of solvers for some classes of multidimensional partial differential equations encountered in kinetic theory modeling of complex fluids. J Non-Newtonian Fluid Mech 139:153–176 2014, 1:11 Cueto and Chinesta Advanced Modeling and Simulation in Engineering Sciences Page 18 of 18 http://www.amses-journal.com/content/1/1/11 77. Ammar A, Mokdad B, Chinesta F, Keunings R (2007) A new family of solvers for some classes of multidimensional partial differential equations encountered in kinetic theory modeling of complex fluids. part ii: transient simulation using space-time separated representations. J Non-Newtonian Fluid Mech 144:98–121 78. Chinesta F, Ammar A, Cueto E (2010) Recent advances in the use of the proper generalized decomposition for solving multidimensional models. Arch Comput Meth Eng 17(4):327–350 79. Chinesta F, Ladeveze P, Cueto E (2011) A short review on model order reduction based on proper generalized decomposition. Arch Comput Meth Eng 18:395–404 80. Chinesta F, Leygue A, Bordeu F, Aguado J, Cueto E, Gonzalez D, Alfaro I, Ammar A, Huerta A (2013) Pgd-based computational vademecum for efficient design, optimization and control. Arch Comput Meth Eng 20(1):31–59. doi:10.1007/s11831-013-9080-x, http://dx.doi.org/10.1007/s11831-013-9080-x 81. Ammar A, Chinesta F, Diez P, Huerta A (2010) An error estimator for separated representations of highly multidimensional models. Comput Meth Appl Mech Eng 199(25–28):1872–1880. doi:http://dx.doi.org/10.1016/j.cma.2010.02.012, http://www.sciencedirect.com/science/article/pii/ S0045782510000708 82. Ladeveze P, Chamoin L (2011) On the verification of model reduction methods based on the proper generalized decomposition. Comput Methods Appl Mech Eng 200(23–24):2032–2047. doi:http://dx.doi.org/10.1016/j.cma.2011.02.019, http://www.sciencedirect.com/science/article/pii/ S0045782511001058 83. Nouy A (2010) A priori model reduction through proper generalized decomposition for solving time-dependent partial differential equations. Comput Meth Appl Mech Eng 199(23–24):1603–1626. doi:http://dx.doi.org/10.1016/j.cma.2010.01.009, http://www.sciencedirect.com/science/article/pii/ S0045782510000186 84. Niroomandi S, Gonzalez D, Alfaro I, Bordeu F, Leygue A, Cueto E, Chinesta F (2013b) Real-time simulation of biological soft tissues: a pgd approach. Int J Numer Methods Biomed Eng 29(5):586–600. doi:10.1002/cnm.2544, http://dx.doi.org/10.1002/cnm.2544 85. Niroomandi S, Alfaro I, Gonzalez D, Cueto E, Chinesta F (2013) Model order reduction in hyperelasticity: a proper generalized decomposition approach. Int J Numer Meth Eng. doi:10.1002/nme.4531, http://dx.doi.org/10.1002/nme. 86. Ammar A, Chinesta F, Cueto E (2011) Coupling finite elements and proper generalized decompositions. Int J Multiscale Comput Eng 9(1):17–33 87. Hunter P, Nielsen P (2005) A strategy for integrative computational physiology. Physiology 20(5): 316–325. doi:10.1152/physiol.00022.2005, http://physiologyonline.physiology.org/content/20/5/316.abstract, http:// physiologyonline.physiology.org/content/20/5/316.full.pdf+html 88. Hunter P (2006) Modeling human physiology: The iups/embs physiome project. Proc IEEE 94(4):678–691. doi:10.1109/JPROC.2006.871767 89. Bassingthwaighte J, Hunter P, Noble D (2009) The cardiac physiome: perspectives for the future. Exp Physiol 94(5): 597–605. doi:10.1113/expphysiol.2008.044099, http://ep.physoc.org/content/94/5/597.abstract, http://ep.physoc. org/content/94/5/597.full.pdf+html 90. Wilson N, Wang K, Dutton RW, Taylor C (2001) A software framework for creating patient specific geometric models from medical imaging data for simulation based medical planning of vascular surgery. In: Niessen W, Viergever M (eds) Medical Image Computing and Computer-Assisted Intervention, MICCAI 2001, Lecture Notes in Computer Science, vol 2208. Springer Berlin Heidelberg, pp 449–456. doi:10.1007/3-540-45468-3_54, http://dx.doi.org/10. 1007/3-540-45468-3_54 91. Guiatni M, Riboulet V, Duriez C, Kheddar A, Cotin S (2013) A combined force and thermal feedback interface for minimally invasive procedures simulation. IEEE ASME Trans Mechatron 18(3):1170–1181. doi:10.1109/TMECH.2012.2197862 92. Turner TE, Schnell S, Burrage K (2004) Stochastic approaches for modelling in vivo reactions. Comput Biol Chem 28: 165–178 93. Sreenath SN, Cho KH, Wellstead P (2008) Modelling the dynamics of signalling pathways. Essays Biochem 45:1–28 94. Munsky B, Khammash M (2006) The finite state projection algorithm for the solution of the chemical master equation. J Chem Phys 124(4):044104. doi:10.1063/1.2145882, http://link.aip.org/link/?JCP/124/044104/1 95. Laughlin RB, Pines D (2000) The theory of everything. Proc Natl Acad Sci 97(1): 28–31. doi:10.1073/pnas.97.1.28, http://www.pnas.org/content/97/1/28.abstract, http://www.pnas.org/content/97/1/28.full.pdf+html 96. Gillespie DT (1977) Exact stochastic simulation of coupled chemical reactions. J Phys Chem 81(25):2340–2361. doi:10.1021/j100540a008, http://pubs.acs.org/doi/abs/10.1021/j100540a008, http://pubs.acs.org/doi/pdf/10.1021/ j100540a008 97. Gillespie DT (2001) Approximate accelerated stochastic simulation of chemically reacting systems. J Chem Phys 115: 1716–1733 98. Ammar A, Cueto E, Chinesta F (2012) Reduction of the chemical master equation for gene regulatory networks using proper generalized decompositions. Int J Numer Methods Biomed Eng 28(9):960–973 doi:10.1186/2213-7467-1-11 Cite this article as: Cueto and Chinesta: Real time simulation for computational surgery: a review. Advanced Modeling and Simulation in Engineering Sciences 2014 1:11.
"Advanced Modeling and Simulation in Engineering Sciences" – Springer Journals
Published: Dec 1, 2014
You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!
Read and print from thousands of top scholarly journals.
Already have an account? Log in
Bookmark this article. You can see your Bookmarks on your DeepDyve Library.
To save an article, log in first, or sign up for a DeepDyve account if you don’t already have one.
Copy and paste the desired citation format or use the link below to download a file formatted for EndNote
Access the full text.
Sign up today, get DeepDyve free for 14 days.
All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.