The effects of the no-touch gap on the no-touch bipolar radiofrequency ablation treatment of liver cancer : a numerical study using a two compartment model
- Authors: Yap, Shelley , Cheong, Jason , Foo, Ji , Ooi, Ean Tat , Ooi, Ean Hin
- Date: 2020
- Type: Text , Journal article
- Relation: Applied Mathematical Modelling Vol. 78, no. (2020), p. 134-147
- Full Text: false
- Reviewed:
- Description: The no-touch bipolar radiofrequency ablation (RFA) for cancer treatment is advantageous primarily because of its capability to prevent tumour track seeding (TTS). In this technique, the RF probes are placed at a distance (no-touch gap) away from the tumour boundary. Ideally, the RF probes should be placed sufficiently far from the tumour in order to avoid TTS. However, having a gap that is too large can lead to ineffective ablation. This paper investigates how the selection of the no-touch gap can affect the tissue electrical and thermal responses during the no-touch bipolar RFA treatment. Simulations were carried out on a two compartment model using the finite element method. Results obtained indicated that a gap that is too large may lead to incomplete ablation and failure to achieve significant ablation margin. However, keeping the gap to be too small may not be clinically practical. It was suggested that the incomplete ablation and the insufficient ablation margin observed in some of the cases may require the placement of additional probes around the tumour. The present study stresses on the importance of identifying the optimal no-touch gap that can avoid TTS without compromising the treatment outcome. © 2019 Elsevier Inc.
Comparison between single- and dual-porosity models for fluid transport in predicting lesion volume following saline-infused radiofrequency ablation
- Authors: Ooi, Ean Hin , Chia, Nicholas , Ooi, Ean Tat , Foo, Ji , Liao, Imam , Nair, Shalini , Mohd Ali, Ahmad
- Date: 2018
- Type: Text , Journal article
- Relation: International Journal of Hyperthermia Vol. 34, no. 8 (2018), p. 1142-1156
- Full Text: false
- Reviewed:
- Description: A recent study by Ooi and Ooi (EH Ooi, ET Ooi, Mass transport in biological tissues: Comparisons between single- and dual-porosity models in the context of saline-infused radiofrequency ablation, Applied Mathematical Modelling, 2017, 41, 271-284) has shown that single-porosity (SP) models for describing fluid transport in biological tissues significantly underestimate the fluid penetration depth when compared to dual-porosity (DP) models. This has raised some concerns on whether the SP model, when coupled with models of radiofrequency ablation (RFA) to simulate saline-infused RFA, could lead to an underestimation of the coagulation size. This paper compares the coagulation volumes obtained following saline-infused RFA predicted based on the SP and DP models for fluid transport. Results showed that the SP model predicted coagulation zones that are consistently 0.5 to 0.9 times smaller than that of DP model. This may be explained by the low permeability value of the tissue interstitial space, which causes the majority of the saline to flow through the vasculature. The absence of fluid flow tracking in the vasculature in the SP model meant that any flow of saline into the vasculature is treated as losses and do not contribute to the saline penetration depth of the tissue. Comparisons with experimental results from the literature revealed that the DP models predicted coagulation zone sizes that are closer to the experimental values than the SP models. This supports the hypothesis that the SP model is a poor choice for simulating the outcome of saline-infused RFA.
Crack propagation modelling in concrete using the scaled boundary finite element method with hybrid polygon-quadtree meshes
- Authors: Ooi, Ean Tat , Natarajan, Sundararajan , Song, Chongmin , Ooi, Ean Hin
- Date: 2017
- Type: Text , Journal article
- Relation: International Journal of Fracture Vol. 203, no. 1-2 (2017), p. 135-157
- Full Text: false
- Reviewed:
- Description: This manuscript presents an extension of the recently-developed hybrid polygon-quadtree-based scaled boundary finite element method to model crack propagation in concrete. This hybrid approach combines the use of quadtree cells with arbitrary sided polygons for domain discretization. The scaled boundary finite element formulation does not distinguish between quadtree cells and arbitrary sided polygons in the mesh. A single formulation is applicable to all types of cells and polygons in the mesh. This eliminates the need to develop transitional elements to bridge the cells belonging to different levels in the quadtree hierarchy. Further to this, the use of arbitrary sided polygons facilitate the accurate discretization of curved boundaries that may result during crack propagation. The fracture process zone that is characteristic in concrete fracture is modelled using zero-thickness interface elements that are coupled to the scaled boundary finite element method using a shadow domain procedure. The scaled boundary finite element method can accurately model the asymptotic stress field in the vicinity of the crack tip with cohesive tractions. This leads to the accurate computation of the stress intensity factors, which is used to determine the condition for crack propagation and the resulting direction. Crack growth can be efficiently resolved using an efficient remeshing algorithm that employs a combination of quadtree decomposition functions and simple Booleans operations. The flexibility of the scaled boundary finite element method to be formulated on arbitrary sided polygons also result in a flexible remeshing algorithm for modelling crack propagation. The developed method is validated using three laboratory experiments of notched concrete beams subjected to different loading conditions.
A computational model to investigate the influence of electrode lengths on the single probe bipolar radiofrequency ablation of the liver
- Authors: Cheong, Jason , Yap, Shelley , Ooi, Ean Tat , Ooi, Ean Hin
- Date: 2019
- Type: Text , Journal article
- Relation: Computer Methods and Programs in Biomedicine Vol. 176, no. (2019), p. 17-32
- Full Text: false
- Reviewed:
- Description: Background and objectives: Recently, there have been calls for RFA to be implemented in the bipolar mode for cancer treatment due to the benefits it offers over the monopolar mode. These include the ability to prevent skin burns at the grounding pad and to avoid tumour track seeding. The usage of bipolar RFA in clinical practice remains uncommon however, as not many research studies have been carried out on bipolar RFA. As such, there is still uncertainty in understanding the effects of the different RF probe configurations on the treatment outcome of RFA. This paper demonstrates that the electrode lengths have a strong influence on the mechanics of bipolar RFA. The information obtained here may lead to further optimization of the system for subsequent uses in the hospitals. Methods: A 2D model in the axisymmetric coordinates was developed to simulate the electro-thermophysiological responses of the tissue during a single probe bipolar RFA. Two different probe configurations were considered, namely the configuration where the active electrode is longer than the ground and the configuration where the ground electrode is longer than the active. The mathematical model was first verified with an existing experimental study found in the literature. Results: Results from the simulations showed that heating is confined only to the region around the shorter electrode, regardless of whether the shorter electrode is the active or the ground. Consequently, thermal coagulation also occurs in the region surrounding the shorter electrode. This opened up the possibility for a better customized treatment through the development of RF probes with adjustable electrode lengths. Conclusions: The electrode length was found to play a significant role on the outcome of single probe bipolar RFA. In particular, the length of the shorter electrode becomes the limiting factor that influences the mechanics of single probe bipolar RFA. Results from this study can be used to further develop and optimize bipolar RFA as an effective and reliable cancer treatment technique. (C) 2019 Elsevier B.V. All rights reserved.
A quadtree-polygon-based scaled boundary finite element method for crack propagation modeling in functionally graded materials
- Authors: Chen, Xiaojun , Luo, Tao , Ooi, Ean Tat , Ooi, Ean Hin , Song, Chongmin
- Date: 2018
- Type: Text , Journal article
- Relation: Theoretical and Applied Fracture Mechanics Vol. 94, no. (2018), p. 120-133
- Full Text: false
- Reviewed:
- Description: This paper presents a method to improve the computational efficiency of the scaled boundary finite element formulation for functionally graded materials. Both isotropic and orthotropic functionally graded materials are considered. This is achieved using a combination of quadtree and polygon meshes. This hybrid meshing approach is particularly suitable to be used with the SBFEM for functionally graded materials because of the significant amount of calculations required to compute the stiffness matrices of the polygons/cells in the mesh. When a quadtree structure is adopted, most of the variables required for the numerical simulation can be pre-computed and stored in the memory, retrieved and scaled as required during the computations, leading to an efficient method for crack propagation modeling. The scaled boundary finite element formulation enables accurate computation of the stress intensity factors directly from the stress solutions without any special post-processing techniques or local mesh refinement in the vicinity of the crack tip. Numerical benchmarks demonstrate the efficiency of the proposed method as opposed to using a purely polygon-mesh based approach. © 2018 Elsevier Ltd
The effects of electrical and thermal boundary condition on the simulation of radiofrequency ablation of liver cancer for tumours located near to the liver boundary
- Authors: Ooi, Ean Hin , Lee, Khiy , Yap, Shelley , Khattab, Mahmoud , Liao, Iman , Ooi, Ean Tat , Foo, Ji , Nair, Shalini , Ali, Ahmad
- Date: 2019
- Type: Text , Journal article
- Relation: Computers in Biology and Medicine Vol. 106, no. (2019), p. 12-23
- Full Text: false
- Reviewed:
- Description: Effects of different boundary conditions prescribed across the boundaries of radiofrequency ablation (RFA) models of liver cancer are investigated for the case where the tumour is at the liver boundary. Ground and Robin-type conditions (electrical field) and body temperature and thermal insulation (thermal field) conditions are examined. 3D models of the human liver based on publicly-available CT images of the liver are developed. An artificial tumour is placed inside the liver at the boundary. Simulations are carried out using the finite element method. The numerical results indicated that different electrical and thermal boundary conditions led to different predictions of the electrical potential, temperature and thermal coagulation distributions. Ground and body temperature conditions presented an unnatural physical conditions around the ablation site, which results in more intense Joule heating and excessive heat loss from the tissue. This led to thermal damage volumes that are smaller than the cases when the Robin type or the thermal insulation conditions are prescribed. The present study suggests that RFA simulations in the future must take into consideration the choice of the type of electrical and thermal boundary conditions to be prescribed in the case where the tumour is located near to the liver boundary.
Dynamic fracture simulations using the scaled boundary finite element method on hybrid polygon-quadtree meshes
- Authors: Ooi, Ean Tat , Natarajan, Sundararajan , Song, Chongmin , Ooi, Ean Hin
- Date: 2016
- Type: Text , Journal article
- Relation: International Journal of Impact Engineering Vol. 90, no. (2016), p. 154-164
- Full Text: false
- Reviewed:
- Description: In this paper, we present an efficient computational procedure to model dynamic fracture within the framework of the scaled boundary finite element method (SBFEM). A quadtree data structure is used to discretise the domain, and 2:1 ratio between the cells is maintained. This limits the number of patterns in the quadtree decomposition and allows for efficient computation of the system matrices. The regions close to the boundary are discretised with arbitrary sided polygons so as to facilitate accurate modelling of the curved boundaries. The stiffness and the mass matrix over all the cells are computed by the SBFEM. Moreover, the semi-analytical nature of the SBFEM enables accurate modelling of the asymptotic stress fields in the vicinity of the crack tip. An efficient remeshing algorithm that combines the quadtree decomposition with simple Boolean operations is proposed to model the crack propagation. The remeshing is restricted only to a small region in the vicinity of the crack tip. The efficiency and the convergence properties of the proposed framework are demonstrated with a few benchmark problems. © 2015 Elsevier Ltd. All rights reserved.
Numerical investigation of the meshless radial basis integral equation method for solving 2D anisotropic potential problems
- Authors: Ooi, Ean Hin , Ooi, Ean Tat , Ang, Whye Teong
- Date: 2015
- Type: Text , Journal article
- Relation: Engineering Analysis with Boundary Elements Vol. 53, no. (2015), p. 27-39
- Full Text: false
- Reviewed:
- Description: The radial basis integral equation (RBIE) method is derived for the first time to solve potential problems involving material anisotropy. The coefficients of the anisotropic conductivity require the gradient term to be modified accordingly when deriving the boundary integral equation so that the flux expression can be properly accounted. Analyses of the behavior of the anisotropic fundamental solution and its spatial gradients showed that their variations along the subdomain boundaries may be large and they increase as the diagonal components of the material anisotropy become larger. The accuracy of the anisotropic RBIE was found to depend primarily on the accuracy of the influence coefficients evaluations and this precedes the number of nodes used. Root mean squared errors of less than 10(-4)% can be obtained if evaluations of the influence coefficients are sufficiently accurate. An alternative formulation of the anisotropic RBIE was derived. The levels of accuracy obtained were not significantly different from the standard formulation. (C) 2014 Elsevier Ltd. All rights reserved.
Shape-shifting thermal coagulation zone during saline-infused radiofrequency ablation: A computational study on the effects of different infusion location
- Authors: Kho, Antony , Foo, Ji , Ooi, Ean Tat , Ooi, Ean Hin
- Date: 2020
- Type: Text , Journal article
- Relation: Computer Methods and Programs in Biomedicine Vol. 184, no. (2020), p.
- Full Text: false
- Reviewed:
- Description: Background and objective: The majority of the studies on radiofrequency ablation (RFA) have focused on enlarging the size of the coagulation zone. An aspect that is crucial but often overlooked is the shape of the coagulation zone. The shape is crucial because the majority of tumours are irregularly-shaped. In this paper, the ability to manipulate the shape of the coagulation zone following saline-infused RFA by altering the location of saline infusion is explored. Methods: A 3D model of the liver tissue was developed. Saline infusion was described using the dual porosity model, while RFA was described using the electrostatic and bioheat transfer equations. Three infusion locations were investigated, namely at the proximal end, the middle and the distal end of the electrode. Investigations were carried out numerically using the finite element method. Results: Results indicated that greater thermal coagulation was found in the region of tissue occupied by the saline bolus. Infusion at the middle of the electrode led to the largest coagulation volume followed by infusion at the proximal and distal ends. It was also found that the ability to delay roll-off, as commonly associated with saline-infused RFA, was true only for the case when infusion is carried out at the middle. When infused at the proximal and distal ends, the occurrence of roll-off was advanced. This may be due to the rapid and more intense heating experienced by the tissue when infusion is carried out at the electrode ends where Joule heating is dominant. Conclusion: Altering the location of saline infusion can influence the shape of the coagulation zone following saline-infused RFA. The ability to ‘shift’ the coagulation zone to a desired location opens up great opportunities for the development of more precise saline-infused RFA treatment that targets specific regions within the tissue. © 2019 Elsevier B.V.
A dual scaled boundary finite element formulation over arbitrary faceted star convex polyhedra
- Authors: Ooi, Ean Tat , Saputra, Albert , Natarajan, Sundararajan , Ooi, Ean Hin , Song, Chongmin
- Date: 2020
- Type: Text , Journal article
- Relation: Computational Mechanics Vol. 66, no. 1 (2020), p. 27-47
- Full Text: false
- Reviewed:
- Description: A novel technique to formulate arbritrary faceted polyhedral elements in three-dimensions is presented. The formulation is applicable for arbitrary faceted polyhedra, provided that a scaling requirement is satisfied and the polyhedron facets are planar. A triangulation process can be applied to non-planar facets to generate an admissible geometry. The formulation adopts two separate scaled boundary coordinate systems with respect to: (i) a scaling centre located within a polyhedron and; (ii) a scaling centre on a polyhedron’s facets. The polyhedron geometry is scaled with respect to both the scaling centres. Polygonal shape functions are derived using the scaled boundary finite element method on the polyhedron facets. The stiffness matrix of a polyhedron is obtained semi-analytically. Numerical integration is required only for the line elements that discretise the polyhedron boundaries. The new formulation passes the patch test. Application of the new formulation in computational solid mechanics is demonstrated using a few numerical benchmarks. © 2020, Springer-Verlag GmbH Germany, part of Springer Nature.
Thermal and thermal damage responses during switching bipolar radiofrequency ablation employing bipolar needles : a computational study on the effects of different electrode configuration, input voltage and ablation duration
- Authors: Cheong, Jason , Ooi, Ean Hin , Ooi, Ean Tat
- Date: 2020
- Type: Text , Journal article
- Relation: International Journal for Numerical Methods in Biomedical Engineering Vol. 36, no. 9 (2020), p.
- Full Text: false
- Reviewed:
- Description: Recent studies have demonstrated the effectiveness of switching bipolar radiofrequency ablation (bRFA) in treating liver cancer. Nevertheless, the clinical use of the treatment remains less common than conventional monopolar RFA – likely due to the lack of understanding of how the tissues respond thermally to the switching effect. The problem is exacerbated by the numerous possible switching combinations when bRFA is performed using bipolar needles, thus making theoretical deduction and experimental studies difficult. This article addresses this issue via computational modelling by examining if significant variation in the treatment outcome exists amongst six different electrode configurations defined by the X-, C-, U-, N-, Z- and O-models. Results indicated that the tissue thermal and thermal damage responses varied depending on the electrode configuration and the operating conditions (input voltage and ablation duration). For a spherical tumour, 30 mm in diameter, complete ablation could not be attained in all configurations with 70 V input voltage and 5 minutes ablation duration. Increasing the input voltage to 90 V enlarged the coagulation zone in the X-model only. With the other configurations, extending the ablation duration to 10 minutes was found to be the better at enlarging the coagulation zone. © 2020 John Wiley & Sons, Ltd.
A quadtree-polygon-based scaled boundary finite element method for image-based mesoscale fracture modelling in concrete
- Authors: Guo, H. , Ooi, Ean Tat , Saputra, Albert , Yang, Zhenjun , Natarajan, Sundararajan , Ooi, Ean Hin , Song, Chongmin
- Date: 2019
- Type: Text , Journal article , acceptedVersion
- Relation: Engineering Fracture Mechanics Vol. 211, no. (2019), p. 420-441
- Full Text:
- Reviewed:
- Description: A quadtree-polygon scaled boundary finite element-based approach for image-based modelling of concrete fracture at the mesoscale is developed. Digital images representing the two-phase mesostructure of concrete, which comprises of coarse aggregates and mortar are either generated using a take-and-place algorithm with a user-defined aggregate volume ratio or obtained from X-ray computed tomography as an input. The digital images are automatically discretised for analysis by applying a balanced quadtree decomposition in combination with a smoothing operation. The scaled boundary finite element method is applied to model the constituents in the concrete mesostructure. A quadtree formulation within the framework of the scaled boundary finite element method is advantageous in that the displacement compatibility between the cells are automatically preserved even in the presence of hanging nodes. Moreover, the geometric flexibility of the scaled boundary finite element method facilitates the use of arbitrary sided polygons, allowing better representation of the aggregate boundaries. The computational burden is significantly reduced as there are only finite number of cell types in a balanced quadtree mesh. The cells in the mesh are connected to each other using cohesive interface elements with appropriate softening laws to model the fracture of the mesostructure. Parametric studies are carried out on concrete specimens subjected to uniaxial tension to investigate the effects of various parameters e.g. aggregate size distribution, porosity and aggregate volume ratio on the fracture of concrete at the meso-scale. Mesoscale fracture of concrete specimens obtained from X-ray computed tomography scans are carried out to demonstrate its feasibility.
A computational framework for the multiphysics simulation of microbubble-mediated sonothrombolysis using a forward-viewing intravascular transducer
- Authors: Tan, Zhi , Ooi, Ean Hin , Chiew, Yeong , Foo, Ji , Ng, Eddie , Ooi, Ean Tat
- Date: 2023
- Type: Text , Journal article
- Relation: Ultrasonics Vol. 131, no. (2023), p.
- Full Text: false
- Reviewed:
- Description: Sonothrombolysis is a technique that utilises ultrasound waves to excite microbubbles surrounding a clot. Clot lysis is achieved through mechanical damage induced by acoustic cavitation and through local clot displacement induced by acoustic radiation force (ARF). Despite the potential of microbubble-mediated sonothrombolysis, the selection of the optimal ultrasound and microbubble parameters remains a challenge. Existing experimental studies are not able to provide a complete picture of how ultrasound and microbubble characteristics influence the outcome of sonothrombolysis. Likewise, computational studies have not been applied in detail in the context of sonothrombolysis. Hence, the effect of interaction between the bubble dynamics and acoustic propagation on the acoustic streaming and clot deformation remains unclear. In the present study, we report for the first time the computational framework that couples the bubble dynamic phenomena with the acoustic propagation in a bubbly medium to simulate microbubble-mediated sonothrombolysis using a forward-viewing transducer. The computational framework was used to investigate the effects of ultrasound properties (pressure and frequency) and microbubble characteristics (radius and concentration) on the outcome of sonothrombolysis. Four major findings were obtained from the simulation results: (i) ultrasound pressure plays the most dominant role over all the other parameters in affecting the bubble dynamics, acoustic attenuation, ARF, acoustic streaming, and clot displacement, (ii) smaller microbubbles could contribute to a more violent oscillation and improve the ARF simultaneously when they are stimulated at higher ultrasound pressure, (iii) higher microbubbles concentration increases the ARF, and (iv) the effect of ultrasound frequency on acoustic attenuation is dependent on the ultrasound pressure. These results may provide fundamental insight that is crucial in bringing sonothrombolysis closer to clinical implementation. © 2023 Elsevier B.V.
Fracture analysis of cracked magneto-electro-elastic functionally graded materials using scaled boundary finite element method
- Authors: Nguyen, Duc , Javidan, Fatemeh , Attar, Mohammadmahdi , Natarajan, Sundararajan , Yang, Zhenjun , Ooi, Ean Hin , Song, Chongmin , Ooi, Ean Tat
- Date: 2022
- Type: Text , Journal article
- Relation: Theoretical and Applied Fracture Mechanics Vol. 118, no. (2022), p.
- Full Text: false
- Reviewed:
- Description: This paper develops the scaled boundary finite element method to analyse fracture of functionally graded magneto-electro-elastic materials. Polygon meshes are employed to discretize the domain. No asymptotic solution, local mesh refinement or other special treatments around a crack tip are required to calculate the intensity factors. When the material gradients of the coefficients in the constitutive matrix are expressed as a series of power functions of the scaled boundary coordinates, the stiffness matrices can be integrated analytically. The formulation enables the generalized intensity factors of stress, electric displacement and magnetic induction fields along the radial direction to be represented analytically. This permits the calculation of the generalized intensity factors directly from the scaled boundary finite element solution of the singular stress, electric displacement and magnetic induction fields by following the standard stress recovery procedures in the finite element method. Several numerical benchmarks are presented to validate the proposed technique with the results reported in the literature. © 2022 Elsevier Ltd
Construction of generalized shape functions over arbitrary polytopes based on scaled boundary finite element method's solution of Poisson's equation
- Authors: Xiao, B. , Natarajan, Sundararajan , Birk, Carolin , Ooi, Ean Hin , Song, Chongmin , Ooi, Ean Tat
- Date: 2023
- Type: Text , Journal article
- Relation: International Journal for Numerical Methods in Engineering Vol. 124, no. 17 (2023), p. 3603-3636
- Full Text:
- Reviewed:
- Description: A general technique to develop arbitrary-sided polygonal elements based on the scaled boundary finite element method is presented. Shape functions are derived from the solution of the Poisson's equation in contrast to the well-known Laplace shape functions that are only linearly complete. The application of the Poisson shape functions can be complete up to any specific order. The shape functions retain the advantage of the scaled boundary finite element method allowing direct formulation on polygons with arbitrary number of sides and quadtree meshes. The resulting formulation is similar to the finite element method where each field variable is interpolated by the same set of shape functions in parametric space and differs only in the integration of the stiffness and mass matrices. Well-established finite element procedures can be applied with the developed shape functions, to solve a variety of engineering problems including, for example, coupled field problems, phase field fracture, and addressing volumetric locking in the near-incompressibility limit by adopting a mixed formulation. Application of the formulation is demonstrated in several engineering problems. Optimal convergence rates are observed. © 2023 The Authors. International Journal for Numerical Methods in Engineering published by John Wiley & Sons Ltd.