A computational model of open-irrigated radiofrequency catheter ablation accounting for mechanical properties of the cardiac tissue

Impact of the dispersive patch placement on dissipated power in radiofrequency ablation for pulmonary vein isolation via a virtual patient study

Radiofrequency ablation (RFA) is a minimally invasive technique for treating arrhythmias by interrupting abnormal electrical signals in the heart. Through a catheter tip, it delivers an alternating current that flows through the heart muscle tissue and the blood to a dispersive patch on the patient's skin. This study aims to test the hypothesis that the placement of the dispersive patch affects the efficacy and safety of RFA. By optimizing the patch position, the procedure could be made more effective and less risky for patients. A 3D in-silico model, based on patient imaging data, was developed to examine the effects of dispersive patch (DP) positioning on electric field distribution within cardiac tissue and the torso during RFA. We conducted 80 computer simulations using a CT-segmented torso model, exploring various DP and electrode configurations while applying standard (25 W) and high (90 W) power settings. For each configuration, we assessed the effectiveness of the DP in delivering power to cardiac tissue near the electrode. The main finding indicates that DP efficacy is significantly influenced by the current delivered to cardiac tissue. Notably, using an anterior patch during ablation proved more effective for the posterior left atrium compared to a posterior patch.

Higher-order thermal modeling and computational analysis of laser ablation in anisotropic cardiac tissue

Laser ablation techniques employ fast hyperthermia mechanisms for diseased-tissue removal, characterized by high selectivity, thus preserving the surrounding healthy tissue. The associated modeling approaches are based on classical Fourier-type laws, though a limited predictivity is observed, particularly at fast time scales. Moreover, limited knowledge is available for cardiac tissue compared to radiofrequency approaches. The present work proposes a comprehensive modeling approach for the computational investigation of the key factors involved in laser-based techniques and assessing the outcomes of induced cellular thermal damage in the cardiac context. The study encompasses a comparative finite element study involving various thermal and cellular damage models incorporating optical-thermal couplings, three-state cellular death dynamics, and a second-order heat transfer formulation generalizing the classical Fourier-based heat equation. A parametric investigation of the thermal profiles shows that higher-order models accurately capture temperature dynamics and lesion formation compared with the classical Fourier-based model. The results highlight the critical role of cardiac anisotropy, influencing the shape and extent of thermal damage, while the three-state cell death model effectively describes the transition from reversible to irreversible damage. These findings demonstrate the reliability of higher-order thermal formulations, laying the basis for future investigations of arrhythmia management via in silico approaches.

Computational study on the effect of thermal deformation of myocardium on lesion formation during radiofrequency ablation

Radiofrequency (RF) catheter ablation treats cardiac diseases by inducing thermal lesion of cardiac tissues through radiofrequency energy operating at around 500 kHz. The electromagnetic wavelength is significantly longer than the size of the radiofrequency active electrode, the tissue is heated through resistive heating. During thermal ablation, the coupled thermo-mechanical property of cardiac tissue influencing the contact area between the electrode and tissue plays a crucial role in the formation of thermal lesions, yet the literature often overlooks the effect of thermal deformation. This paper proposes a thermo-hyperelastic constitutive model for myocardium that models thermal contraction and expansion during ablation. Furthermore, a finite element model was established to investigate the effect of the electro-thermo-mechanical coupling property of myocardium on lesion formation under different contact forces. To ensure convergence, we solved the fully coupled electro-thermo-mechanical finite element model using the segregated step method. The computational results demonstrate that thermal deformation, which causes an expansion in the tissue-electrode contact area, increases lesion width and volume, while its influence on lesion depth is negligible. Specifically, after a 30-s ablation under contact forces of 0.1, 0.15, and 0.2 N, the lesion volume increased from 4.53, 7.66, and 10.62 mm3 (without thermo-mechanical coupling) to 5.36, 8.33, and 13.34 mm3 (with thermo-mechanical coupling), respectively. Similarly, the lesion width increased from 2.68, 3.12, and 3.44 mm to 2.78, 3.22, and 3.62 mm. Moreover, both thermal deformation and contact force exert a minimal effect on lesion formation time.

Ablation catheter-induced mechanical deformation in myocardium: computer modeling and ex vivo experiments

Cardiac catheter ablation requires an adequate contact between myocardium and catheter tip. Our aim was to quantify the relationship between the contact force (CF) and the resulting mechanical deformation induced by the catheter tip using an ex vivo model and computational modeling. The catheter tip was inserted perpendicularly into porcine heart samples. CF values ranged from 10 to 80 g. The computer model was built to simulate the same experimental conditions, and it considered a 3-parameter Mooney-Rivlin model based on hyper-elastic material. We found a strong correlation between the CF and insertion depth (ID) (R2 = 0.96, P < 0.001), from 0.7 ± 0.3 mm at 10 g to 6.9 ± 0.1 mm at 80 g. Since the surface deformation was asymmetrical, two transversal diameters (minor and major) were identified. Both diameters were strongly correlated with CF (R2 ≥ 0.95), from 4.0 ± 0.4 mm at 20 g to 10.3 ± 0.0 mm at 80 g (minor), and from 6.4 ± 0.7 mm at 20 g to 16.7 ± 0.1 mm at 80 g (major). An optimal fit between computer and experimental results was achieved, with a prediction error of 0.74 and 0.86 mm for insertion depth and mean surface diameter, respectively.

Effect of dispersive electrode position (anterior vs. posterior) in epicardial radiofrequency ablation of ventricular wall: A computer simulation study

An epicardial approach is often used in radiofrequency (RF) catheter ablation to ablate ventricular tachycardia when an endocardial approach fails. Our objective was to analyze the effect of the position of the dispersive patch (DP) on lesion size using computer modeling during epicardial approach. We compared the posterior position (patient's back), commonly used in clinical practice, to the anterior position (patient's chest). The model considered ventricular wall thicknesses between 4 and 8 mm, and electrode insertion depths between .3 and .7 mm. RF pulses were simulated with 20 W of power for 30 s duration. Statistically significant differences (P < .001) were found between both DP positions in terms of baseline impedance, RF current (at 15 s) and thermal lesion size. The anterior position involved lower impedance (130.8 ± 4.7 vs. 146.2 ± 4.9 Ω) and a higher current (401.5 ± 5.6 vs. 377.5 ± 5.1 mA). The anterior position created lesion sizes larger than the posterior position: 8.9 ± 0.4 vs. 8.4 ± 0.4 mm in maximum width, 8.6 ± 0.4 vs. 8.1 ± 0.4 mm in surface width, and 4.5 ± 0.4 vs. 4.3 ± 0.4 mm in depth. Our results suggest that: (1) the redirection of the RF currents due to repositioning the PD has little impact on lesion size and only affects baseline impedance, and (2) the differences in lesion size are only 0.5 mm wider and 0.2 mm deeper for the anterior position, which does not seem to have a clinical impact in the context of VT ablation.

Computer Simulation of Catheter Cryoablation for Pulmonary Vein Isolation

Cryoablation is a well-established medical procedure for surgically treating atrial fibrillation. Cryothermal catheter therapy induces cellular necrosis by freezing the insides of pulmonary veins, with the goal of disrupting abnormal electrical heart signals. Nevertheless, tissue damage induced by cold temperatures may also lead to other complications after cardiac surgery. In this sense, the simulation of catheter ablation can provide safer environments for training and the performance of cryotherapy interventions. Therefore, in this paper, we propose a novel approach to help better understand how temperature rates can affect this procedure by using computer tools to develop a simulation framework to predict lesion size and determine optimal temperature conditions for reducing the risk of major complications. The results showed that a temperature profile of around -40 °C caused less penetration, reduced necrotic damage, and smaller lesion size in the tissue. Instead, cryotherapy close to -60 °C achieved a greater depth of temperature flow inside the tissue and a larger cross-section area of the lesion. With further development and validation, the framework could represent a cost-effective strategy for providing personalized modeling, better planning of cryocatheter-based treatment, and preventing surgical complications.

Meshless Simulation of Multi-site Radio Frequency Catheter Ablation through the Fragile Points Method

Computational models for radio frequency catheter ablation (RFCA) of cardiac arrhythmia have been developed and tested in conditions where a single ablation site is considered. However, in reality arrhythmic events are generated at multiple sites which are ablated during treatment. Under such conditions, heat accumulation from several ablations is expected and models should take this effect into account. Moreover, such models are solved using the Finite Element Method which requires a good quality mesh to ensure numerical accuracy. Therefore, clinical application is limited since heat accumulation effects are neglected and numerical accuracy depends on mesh quality. In this work, we propose a novel meshless computational model where tissue heat accumulation from previously ablated sites is taken into account. In this way, we aim to overcome the mesh quality restriction of the Finite Element Method and enable realistic multi-site ablation simulation. We consider a two ablation sites protocol where tissue temperature at the end of the first ablation is used as initial condition for the second ablation. The effect of the time interval between the ablation of the two sites is evaluated. The proposed method demonstrates that previous models that do not account for heat accumulation between ablations may underestimate the tissue heat distribution.Clinical relevance- The proposed computational model may be used to build and update a heat map for ablation guidance taking into account the contribution from previously ablated sites. Being a meshless model, it does not require significant input from the user during preprocessing. Therefore, it is suitable for application in a clinical setting.

Comparative Analysis of Temperature Rise between Convective Heat Transfer Method and Computational Fluid Dynamics Method in an Anatomy-Based Left Atrium Model during Pulsed Field Ablation: A Computational Study

The non-thermal effects are considered one of the prominent advantages of pulsed field ablation (PFA). However, at higher PFA doses, the temperature rise in the tissue during PFA may exceed the thermal damage threshold, at which time intracardiac pulsatile blood flow plays a crucial role in suppressing this temperature rise. This study aims to compare the effect of heat dissipation of the different methods in simulating the pulsatile blood flow during PFA. This study first constructed an anatomy-based left atrium (LA) model and then applied the convective heat transfer (CHT) method and the computational fluid dynamics (CFD) method to the model, respectively, and the thermal convective coefficients used in the CHT method are 984 (W/m2*K) (blood-myocardium interface) and 4372 (W/m2*K) (blood-catheter interface), respectively. Then, it compared the effect of the above two methods on the maximum temperature of myocardium and blood, as well as the myocardial ablation volumes caused by irreversible electroporation (IRE) and hyperthermia under different PFA parameters. Compared with the CFD method, the CHT method underestimates the maximum temperature of myocardium and blood; the differences in the maximum temperature of myocardium and blood between the two methods at the end of the last pulse are significant (>1 °C), and the differences in the maximum temperature of blood at the end of the last pulse interval are significant (>1 °C) only at a pulse amplitude greater than 1000 V or pulse number greater than 10. Under the same pulse amplitude and different heat dissipation methods, the IRE ablation volumes are the same. Compared with the CFD method, the CHT method underestimates the hyperthermia ablation volume; the differences in the hyperthermia ablation volume are significant (>1 mm3) only at a pulse amplitude greater than 1000 V, a pulse interval of 250 ms, or a pulse number greater than 10. Additionally, the hyperthermia ablation isosurfaces are completely wrapped by the IRE ablation isosurfaces in the myocardium. Thus, during PFA, compared with the CFD method, the CHT method cannot accurately simulate the maximum myocardial temperature; however, except at the above PFA parameters, the CHT method can accurately simulate the maximum blood temperature and the myocardial ablation volume caused by IRE and hyperthermia. Additionally, within the range of the PFA parameters used in this study, the temperature rise during PFA may not lead to the appearance of additional hyperthermia ablation areas beyond the IRE ablation area in the myocardium.

Influence of pulsating intracardiac blood flow on radiofrequency catheter ablation outcomes in an anatomy-based atrium model

BACKGROUND

Highly consistent cardiac ablation outcomes through radiofrequency catheter ablation (RFCA) under pulsatile and constant flow profiles (PP&CP) of intracardiac blood were previously indicated by computer modeling, with simplified geometry and lossless receipt of inflow for ablation catheters. This study aimed to further investigate the effects of intracardiac blood pulsatility in an anatomy-based atrium model.

METHODS

Four pulmonary veins were blood inflows at 10 mm Hg. The mitral valve was the outflow, with PP based on pulsatile velocity curve from clinical measurements, and CP was obtained by averaging the velocity curve under PP over an ablation time of 30 s. A numerical comparison between ablation results under PP and CP, without experimental validation, was performed.

RESULTS

Temperature fluctuations persisted in mid-myocardium, and most clearly in blood and endocardium under PP. At a constant power of 20 W, marked differences in ablation outcome between PP and CP occurred in the middle of unilateral pulmonary veins and the posterior wall of the left atrium (LA) where the blood velocities were significantly decreased under CP. The mid-myocardial, blood and endocardial temperatures, as well as the effective lesion volume at the former position, were decreased by 4.1%, 15%, 13.6%, and 13.8%, respectively under PP. The extents for the latter position were 11%, 22%, 22.5%, and 55.6%, respectively.

CONCLUSION

Intracardiac flow pulsatility causes a greater reduction in blood and endocardial temperatures at ablation sites away from the main bloodstream, effective cooling of which is more likely to rely on blood velocities approaching peak PP values.

Impact of electrode tip shape on catheter performance in cardiac radiofrequency ablation

Background

The role of catheter tip shape on the safety and efficacy of radiofrequency (RF) ablation has been overlooked, although differences have been observed in clinical and research fields.

Objective

The purpose of this study was to analyze the role of electrode tip shape in RF ablation using a computational model.

Methods

We simulated 108 RF ablations through a realistic 3-dimensional computational model considering 2 clinically used, open-irrigated catheters (spherical and cylindrical tip), varying contact force (CF), blood flow, and irrigation. Lesions are defined by the 50°C isotherm contour and evaluated by means of width, depth, depth at maximum width, and volume. Ablations are deemed as safe, critical (tissue temperature >90°C), and pop (tissue temperature >100°C).

Results

Tissue-electrode contact is less for the spherical tip at low CF but the relationship is inverted at high CF. At low CF, the cylindrical tip generates deeper and wider lesions and a 4-fold larger volume. With increasing CF, the lesions generated by the spherical tip become comparable to those generated by the cylindrical tip. The 2 tips feature different safety profiles: CF and power are the main determinants of pops for the spherical tip; power is the main factor for the cylindrical tip; and CF has a marginal effect. The cylindrical tip is more prone to pop generation at higher powers. Saline irrigation and blood flow effect do not depend on tip shape.

Conclusion

Tip shape determines the performance of ablation catheters and has a major impact on their safety profile. The cylindrical tip shows more predictable behavior in a wide range of CF values.

Computer simulations of consecutive radiofrequency pulses applied at the same point during cardiac catheter ablation: Implications for lesion size and risk of overheating

BACKGROUND AND OBJECTIVES

To study temperature distribution and lesion size during two repeated radiofrequency (RF) pulses applied at the same point in the context of RF cardiac ablation (RFCA).

METHODS

An in-silico RFCA model accounting for reversible and irreversible changes in myocardium electrical properties due to RF-induced heating. Arrhenius damage model to estimate lesion size during the application of two 20 W pulses at intervals (INT) of from 5 to 70 s. We considered two pulse durations: 20 s and 30 s.

RESULTS

INT has a significant effect on lesion size and maximum tissue temperature (TMAX). The shorter the INT the greater the increase in lesion size after the second pulse but also the greater the TMAX. If the second pulse is applied almost immediately (INT=5 s), depth increases 1.4 mm and 1.5 mm for pulses of 20 s and 30 s, respectively. If INT is longer than 30 s it increases 1.1 mm and 1.3 mm for pulses of 20 s and 30 s, respectively. While a single 20 s pulse causes T_MAX=79 ºC, a second pulse produces values of from 92 to 96 ºC (the higher the temperature the shorter the INT). For 30 s pulses, T_MAX=93 ºC for a single pulse, and varied from 98 to 104 ºC for a second pulse.

CONCLUSIONS

Applying a second RF pulse at the same ablation site increases lesion depth by 1 - 1.5 mm more than a single pulse and could lead to higher temperatures (up to 17 ºC). Both lesion depth and maximum tissue temperature increased at shorter inter-pulse intervals, which could cause clinical complications from overheating such as steam pops.

Multiscale and Multiphysics Modeling of Anisotropic Cardiac RFCA: Experimental-Based Model Calibration via Multi-Point Temperature Measurements

Radiofrequency catheter ablation (RFCA) is the mainstream treatment for drug-refractory cardiac fibrillation. Multiple studies demonstrated that incorrect dosage of radiofrequency energy to the myocardium could lead to uncontrolled tissue damage or treatment failure, with the consequent need for unplanned reoperations. Monitoring tissue temperature during thermal therapy and predicting the extent of lesions may improve treatment efficacy. Cardiac computational modeling represents a viable tool for identifying optimal RFCA settings, though predictability issues still limit a widespread usage of such a technology in clinical scenarios. We aim to fill this gap by assessing the influence of the intrinsic myocardial microstructure on the thermo-electric behavior at the tissue level. By performing multi-point temperature measurements on ex-vivo swine cardiac tissue samples, the experimental characterization of myocardial thermal anisotropy allowed us to assemble a fine-tuned thermo-electric material model of the cardiac tissue. We implemented a multiphysics and multiscale computational framework, encompassing thermo-electric anisotropic conduction, phase-lagging for heat transfer, and a three-state dynamical system for cellular death and lesion estimation. Our analysis resulted in a remarkable agreement between ex-vivo measurements and numerical results. Accordingly, we identified myocardium anisotropy as the driving effect on the outcomes of hyperthermic treatments. Furthermore, we characterized the complex nonlinear couplings regulating tissue behavior during RFCA, discussing model calibration, limitations, and perspectives.

Spatial temperature reconstructions in myocardial tissues undergoing radiofrequency ablations by performing high-resolved temperature measurements

BACKGROUND

Radiofrequency (RF) lesion creation is related to the heat propagation induced by RF application on tissues. Thermocouple embedded in the RF antenna are not able to predict deep tissue temperature at various level.

OBJECTIVES

This study aims to investigate the influence of power delivered on radiofrequency catheter ablation (RFCA) effects by means of high resolved 2D temperature maps.

METHODS

Three trials of four ablations (12 applications) were executed on each specimen of healthy excised swine myocardium in different application points at four RF power values (30 W, 40 W, 50 W, and 60 W) for a fixed treatment time. All the data provided by the fiber Bragg gratings (FBGs) were analyzed. Temperature variations (ΔT) in time recorded in the 28 sites of measurements were reported. Also, temperature maps showing the ΔT spatial distribution reached within the tissue at the end of the RFCA were produced and displayed, together with the representation of the lethal isotherm. Moreover, the time of achievement of the lethal isotherm at different tissue depths (from 1 to 8 mm) was evaluated for the four power settings.

RESULTS

Temperature trends reported comparable profiles across the different power settings. ΔT values and ΔT rising times showed dependence on the sensors' proximity to the RF energy source and on the set RF power. Temperature maps confirmed that heat propagation occurs preferentially along the width of the tissue than in the depth. Also, for the adjusted treatment time, no power setting guarantees lesions thicker than 6 mm.

CONCLUSIONS

ΔT maximal values and ΔT rising time strongly depends on the proximity of the tissues to RF energy source, as well as on the RF power setting. A plateau is reached in lesion size, regardless of the power setting. A first correlation between lesion size, power setting, and time to achieve lethal isotherms has been established.

Computer modeling of radiofrequency cardiac ablation: 30 years of bioengineering research

This review begins with a rationale of the importance of theoretical, mathematical and computational models for radiofrequency (RF) catheter ablation (RFCA). We then describe the historical context in which each model was developed, its contribution to the knowledge of the physics of RFCA and its implications for clinical practice. Next, we review the computer modeling studies intended to improve our knowledge of the biophysics of RFCA and those intended to explore new technologies. We describe the most important technical details of the implementation of mathematical models, including governing equations, tissue properties, boundary conditions, etc. We discuss the utility of lumped element models, which despite their simplicity are widely used by clinical researchers to provide a physical explanation of how RF power is absorbed in different tissues. Computer model verification and validation are also discussed in the context of RFCA. The article ends with a section on the current limitations, i.e. aspects not yet included in state-of-the-art RFCA computer modeling and on future work aimed at covering the current gaps.

Systematic Characterization of High-Power Short-Duration Ablation: Insight From an Advanced Virtual Model

Background: High-power short-duration (HPSD) recently emerged as a new approach to radiofrequency (RF) catheter ablation. However, basic and clinical data supporting its effectiveness and safety is still scarce.

Objective: We aim to characterize HPSD with an advanced virtual model, able to assess lesion dimensions and complications in multiple conditions and compare it to standard protocols.

Methods: We evaluate, on both atrium and ventricle, three HPSD protocols (70 W/8 s, 80 W/6 s, and 90 W/4 s) through a realistic 3D computational model of power-controlled RF ablation, varying catheter tip design (spherical/cylindrical), contact force (CF), blood flow, and saline irrigation. Lesions are defined by the 50°C isotherm contour. Ablations are deemed safe or complicated by pop (tissue temperature >97°C) or charring (blood temperature >80°C). We compared HPSD with standards protocols (30–40 W/30 s). We analyzed the effect of a second HPSD application.

Results: We simulated 432 applications. Most (79%) associated a complication, especially in the atrium. The three HPSD protocols performed similarly in the atrium, while 90 W/4 s appeared the safest in the ventricle. Low irrigation rate led frequently to charring (72%). High-power short-duration lesions were 40–60% shallower and smaller in volume compared to standards, although featuring similar width. A second HPSD application increased lesions to a size comparable to standards.

Conclusion: High-power short-duration lesions are smaller in volume and more superficial than standards but comparable in width, which can be advantageous in the atrium. A second application can produce lesions similar to standards in a shorter time. Despite its narrow safety margin, HPSD seems a valuable new clinical approach.

Effect of anisotropy in myocardial electrical conductivity on lesion characteristics during radiofrequency cardiac ablation: a numerical study

BACKGROUND

Traditional computer simulation studies of radiofrequency catheter ablation (RFCA) usually neglect the anisotropy in myocardial electrical conductivity (MEC), which is likely an essential factor in governing the ablation outcome. Here, a numerical study of lesion characteristics during RFCA based on an anatomy-based model incorporating fiber orientation was performed to investigate the anisotropy in MEC.

METHODS

A three-dimensional thorax model including atria, blood, connective tissue, muscle, fat, and skin was constructed. The myocardial fiber was established through a rule-based method (RBM) based on the anatomical structure of the heart. The anisotropic MEC were 0.40 and 0.28 S m^-1 in longitudinal and transverse directions, respectively. The ablation result was compared with the isotropic scenario where the isotropic MEC was the average of the anisotropic conductivities as 0.34 S m^-1.

RESULTS

The complexity of fiber architecture varied with that of the local anatomical structure. At RF power of 20 W for 30 s, the tissue temperature and lesion volume were reduced by 2.8 ± 0.1% and 6.9 ± 0.5%, respectively, under anisotropic MEC around the ostium of the pulmonary vein and left atrial appendage. Those for the posterior wall and roof of the left atrium, and the inside of the superior vena cava were 1.9 ± 0.3% and 5.6 ± 1.2%, respectively.

CONCLUSIONS

Anisotropy in MEC has a greater reduction effect on lesion volume than on tissue temperature during RFCA; this effect tends to be restrained at positions with more uniform fiber distributions and can be enhanced where significant variation in fiber architecture occurred.

Proactive esophageal cooling protects against thermal insults during high-power short-duration radiofrequency cardiac ablation

BACKGROUND

Proactive cooling with a novel cooling device has been shown to reduce endoscopically identified thermal injury during radiofrequency (RF) ablation for the treatment of atrial fibrillation using medium power settings. We aimed to evaluate the effects of proactive cooling during high-power short-duration (HPSD) ablation.

METHODS

A computer model accounting for the left atrium (1.5 mm thickness) and esophagus including the active cooling device was created. We used the Arrhenius equation to estimate the esophageal thermal damage during 50 W/ 10 s and 90 W/ 4 s RF ablations.

RESULTS

With proactive esophageal cooling in place, temperatures in the esophageal tissue were significantly reduced from control conditions without cooling, and the resulting percentage of damage to the esophageal wall was reduced around 50%, restricting damage to the epi-esophageal region and consequently sparing the remainder of the esophageal tissue, including the mucosal surface. Lesions in the atrial wall remained transmural despite cooling, and maximum width barely changed (<0.8 mm).

CONCLUSIONS

Proactive esophageal cooling significantly reduces temperatures and the resulting fraction of damage in the esophagus during HPSD ablation. These findings offer a mechanistic rationale explaining the high degree of safety encountered to date using proactive esophageal cooling, and further underscore the fact that temperature monitoring is inadequate to avoid thermal damage to the esophagus.

A computational model of radiofrequency ablation in the stomach, an emerging therapy for gastric dysrhythmias

Gastric ablation has recently emerged as a promising potential therapy for bioelectrical dysrhythmias that underpin many gastrointestinal disorders. Despite similarities to well-developed cardiac ablation, gastric ablation is in early development and has thus far been limited to temperature-controlled, non-irrigated settings. A computational model of gastric ablation is needed to enable in silico testing and optimization of ablation parameters and techniques. In this study, we developed a computational model of radio-frequency (RF) gastric ablation. Model parameters and boundary conditions were established based on the current in vivo experimental application of serosal gastric ablation with a non-irrigated RF catheter. The Pennes bioheat transfer equation was used to model the thermal component of RF ablation, and Laplace's equation was used to model the Joule heating component. Tissue, blood, and catheter parameters were obtained from literature. The performance of the model was compared to previously established experimental values of temperature measured from various distances from the catheter tip. The model produced temperature estimations that were within 6% of the maximum experimental temperature at 2.5 mm from the catheter, and within 13% of the maximum temperature change at 4.7 mm. This model now provides a computational basis through which to conduct in silico testing of gastric ablation, and can be usefully applied to optimize gastric ablation parameters. In future, the model can be expanded to include irrigation of the catheter tip and power-controlled RF settings.Clinical Relevance- This work presents a computational model of gastric ablation that can now guide the in silico development of effective ablation parameters and therapeutic strategies, expanding the breadth of this promising therapy.

Fluid–Structure Interaction and Non-Fourier Effects in Coupled Electro-Thermo-Mechanical Models for Cardiac Ablation

In this study, a fully coupled electro-thermo-mechanical model of radiofrequency (RF)-assisted cardiac ablation has been developed, incorporating fluid–structure interaction, thermal relaxation time effects and porous media approach. A non-Fourier based bio-heat transfer model has been used for predicting the temperature distribution and ablation zone during the cardiac ablation. The blood has been modeled as a Newtonian fluid and the velocity fields are obtained utilizing the Navier–Stokes equations. The thermal stresses induced due to the heating of the cardiac tissue have also been accounted. Parametric studies have been conducted to investigate the effect of cardiac tissue porosity, thermal relaxation time effects, electrode insertion depths and orientations on the treatment outcomes of the cardiac ablation. The results are presented in terms of predicted temperature distributions and ablation volumes for different cases of interest utilizing a finite element based COMSOL Multiphysics software. It has been found that electrode insertion depth and orientation has a significant effect on the treatment outcomes of cardiac ablation. Further, porosity of cardiac tissue also plays an important role in the prediction of temperature distribution and ablation volume during RF-assisted cardiac ablation. Moreover, thermal relaxation times only affect the treatment outcomes for shorter treatment times of less than 30 s.

Computer simulations of an irrigated radiofrequency cardiac ablation catheter and experimental validation by infrared imaging

PURPOSE

To develop and validate a three-dimensional (3-D) computer model based on accurate geometry of an irrigated cardiac radiofrequency (RF) ablation catheter with microwave radiometry capability, and to test catheter performance.

METHODS

A computer model was developed based on CAD geometry of a RF cardiac ablation catheter prototype to simulate electromagnetic heating, heat transfer, and computational fluid dynamics (blood flow, open irrigation, and natural convection). Parametric studies were performed; blood flow velocity (0-25 cm/s) and irrigation flow (0-40 ml/min) varied, both with perpendicular (PE) and parallel (PA) catheter orientations relative to tissue. Tissue Agar phantom studies were performed under similar conditions, and temperature maps were recorded via infrared camera. Computer model simulations were performed with constant voltage and with voltage adjusted to achieve maximum tissue temperatures of 95-105 °C.

RESULTS

Model predicted thermal lesion width at 5 W power was 5.8-6.4 mm (PE)/6.5-6.6 mm (PA), and lesion depth was 4.0-4.3 mm (PE)/4.0-4.1 mm (PA). Compared to phantom studies, the mean errors of the computer model were as follows: 6.2 °C(PE)/4.3 °C (PA) for maximum gel temperature, 0.7 mm (10.9%) (PE)/0.1 mm (0.8%) (PA) for lesion width, and 0.3 mm (7.7%)(PE)/0.7 mm (19.1%) (PA) for lesion depth. For temperature-controlled ablation, model predicted thermal lesion width was 7-9.2 mm (PE)/8.6-9.2 mm (PA), and lesion depth was 4.3-5.5 mm (PE)/3.4-5.4 mm (PA).

CONCLUSIONS

Computer models were able to reproduce device performance and to enable device evaluation under varying conditions. Temperature controlled ablation of irrigated catheters enables optimal tissue temperatures independent of patient-specific conditions such as blood flow.

Domain Heterogeneity in Radiofrequency Therapies for Pain Relief: A Computational Study with Coupled Models

The objective of the current research work is to study the differences between the predicted ablation volume in homogeneous and heterogeneous models of typical radiofrequency (RF) procedures for pain relief. A three-dimensional computational domain comprising of the realistic anatomy of the target tissue was considered in the present study. A comparative analysis was conducted for three different scenarios: (a) a completely homogeneous domain comprising of only muscle tissue, (b) a heterogeneous domain comprising of nerve and muscle tissues, and (c) a heterogeneous domain comprising of bone, nerve and muscle tissues. Finite-element-based simulations were performed to compute the temperature and electrical field distribution during conventional RF procedures for treating pain, and exemplified here for the continuous case. The predicted results reveal that the consideration of heterogeneity within the computational domain results in distorted electric field distribution and leads to a significant reduction in the attained ablation volume during the continuous RF application for pain relief. The findings of this study could provide first-hand quantitative information to clinical practitioners about the impact of such heterogeneities on the efficacy of RF procedures, thereby assisting them in developing standardized optimal protocols for different cases of interest.