Next Article in Journal
Biochemical Analysis of Urine Samples from Diabetic and Hypertensive Patients without Renal Dysfunction Using Spectrophotometry and Raman Spectroscopy Techniques Aiming Classification and Diagnosis
Next Article in Special Issue
Freezing of Solute-Laden Aqueous Solutions: Kinetics of Crystallization and Heat- and Mass-Transfer-Limited Model
Previous Article in Journal
Diagnostic and Therapeutic Roles of the “Omics” in Hypoxic–Ischemic Encephalopathy in Neonates
Previous Article in Special Issue
Synergetic Thermal Therapy for Cancer: State-of-the-Art and the Future
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Real-Time Temperature Rise Estimation during Irreversible Electroporation Treatment through State-Space Modeling

by
Sabrina N. Campelo
*,
Edward J. Jacobs IV
,
Kenneth N. Aycock
and
Rafael V. Davalos
Virginia Tech-Wake Forest School of Biomedical Engineering and Sciences, Virginia Tech Department of Biomedical Engineering and Mechanics, Virginia Tech, Blacksburg, VA 24061, USA
*
Author to whom correspondence should be addressed.
Bioengineering 2022, 9(10), 499; https://doi.org/10.3390/bioengineering9100499
Submission received: 1 September 2022 / Revised: 14 September 2022 / Accepted: 18 September 2022 / Published: 23 September 2022
(This article belongs to the Special Issue Multiscale Thermal Engineering for Biomedical Applications)

Abstract

:
To evaluate the feasibility of real-time temperature monitoring during an electroporation-based therapy procedure, a data-driven state-space model was developed. Agar phantoms mimicking low conductivity (LC) and high conductivity (HC) tissues were tested under the influences of high (HV) and low (LV) applied voltages. Real-time changes in impedance, measured by Fourier Analysis SpecTroscopy (FAST) along with the known tissue conductivity and applied voltages, were used to train the model. A theoretical finite element model was used for external validation of the model, producing model fits of 95.8, 88.4, 90.7, and 93.7% at 4 mm and 93.2, 58.9, 90.0, and 90.1% at 10 mm for the HV-HC, LV-LC, HV-LC, and LV-HC groups, respectively. The proposed model suggests that real-time temperature monitoring may be achieved with good accuracy through the use of real-time impedance monitoring.

Graphical Abstract

1. Introduction

Irreversible electroporation (IRE) and second-generation high-frequency irreversible electroporation (H-FIRE) are techniques currently being developed to treat malignant tumors when other treatment methods, such as surgical resection or thermal ablation, are not appropriate [1]. It is also being developed for the treatment of arrhythmogenic cardiac tissues [2,3]. This technique uses short (1–100 μs) high-magnitude electric pulses (1–3 kV) to produce an electric field resulting in an increase in transmembrane potential (TMP), leading to cell death within the target zone [4,5,6]. When this TMP limit is reached (∼1 V), nanoscale defects, or “pores”, form in the cellular membrane. When numerous pulses are administered, the formation of large, long-lived pores disrupts cellular homeostasis, resulting in cell death via different mechanisms [7,8]. Because cell death with IRE is primarily dependent on the generated TMP, IRE is classified as a non-thermal ablation method and may thus be used in a variety of settings. An electric field is produced between two or more needle electrodes placed directly into the tumor-containing area prior to pulse administration. In some cases, a single-insertion bipolar probe (two electrodes integrated into the same cylindrical shaft) is employed.
It has been demonstrated with IRE that ablation and thermal effects can be mutually exclusive phenomena [9,10,11]. Innumerable follow-up studies have found appropriate parameters for ablation in various contexts. However, unfavorable side effects caused by temperature increases continue to be the most important issue limiting the size of the IRE ablations [12]. Clinically, discrepancies in pulse delivery methods and clinical expertise might result in variable degrees of thermal tissue damage, which is most likely to blame for the wide range of complication rates and oncological outcomes documented in the literature [13,14,15,16]. Even at low voltages, the steep potential gradient at the electrode borders can generate quite large electric fields (and consequently temperatures), causing tissue coagulation, extracellular protein denaturation, and bleeding [17,18]. This may compromise the safety of therapy if the electrodes are placed in close proximity to or in direct contact with specific anatomical structures sensitive to thermal damage — for instance, in the liver, thermal injury to the bile ducts can lead to life-threatening consequences [17,19].
Despite these reservations, present clinical pulsing methods do not actively take temperature increases into account during IRE treatments. In the absence of real-time feedback during pulse delivery, the risk of thermal injury varies greatly across patients depending on the tissue being treated, electrode design, pulsing paradigm, and patient-specific tissue parameters. A recent study, for example, suggested that around 30% of the typical IRE ablation volume suffers from moderate hyperthermia (40–50 ° C), with 5% exposed to temperatures over 50 ° C [20]. Maintaining the nonthermal elements of IRE will become more challenging as clinicians and researchers strive to increase ablation volumes in order to treat bigger tumors. Other factors to consider include that IRE performed in proximity to a metal stent has demonstrated higher temperatures surrounding the electrodes and viable tissue remnants [21]. Furthermore, current research shows that, even when thermal damage is improbable, lowering mild-to-moderate thermal impacts may encourage greater immune activation, which may improve long-term therapeutic effectiveness and patient outcomes [22,23,24]. As a result, strategies to reduce and/or manage thermal effects in the proximity of the electrodes are critical to improving IRE and other pulsed electric field (PEF)-based ablation modalities for wider clinical usage.
Several thermal mitigation (TM) solutions have been proposed to decrease the temperature increase during PEF therapy. Innovative probe designs that include heat-dissipating technology such as phase-change materials or active cooling [25,26], as well as pulse paradigm changes that enable tissue perfusion to distribute heat between successive pulses [27,28], have demonstrated an ability to reduce thermal damage. However, the advantages of these techniques must be weighed against the increased time, cost, and effort required by modified probes [20,29,30]. While previous studies have utilized computational modeling for prospective and retrospective estimations of temperature rises during treatment, no studies have been conducted to implement real-time temperature monitoring without requiring the insertion of additional sensor devices.
In this study, we develop a mathematical model to estimate real-time tissue temperature changes at distances 4 mm ( T 4 mm ) and 10 mm ( T 10 mm ) from the electrode surface through the integration of Fourier Analysis SpecTroscopy (FAST) [31], a real-time method for making bioimpedance measurements during electroporation-based treatments. FAST performs electrical impedance spectroscopy (EIS) across a large frequency range in real time (i.e., between H-FIRE bursts) using custom, low-voltage, high-bandwidth rectangular waveforms. One advantage to using FAST over conventional EIS is that the same electronics used to deliver the H-FIRE pulses can also be used for data acquisition for FAST, eliminating the need for a separate set of electronics. In addition to the expenses required with additional equipment, the impedance spectrum acquisition time frame of commercial EIS equipment (∼10 s) is significantly longer than that between therapeutic pulses (∼1 s) employed by electroporation-based therapies. Diagnostic FAST waveforms are interlaced between high voltage therapeutic H-FIRE bursts using a pulse generator. This diagnostic waveform consists of a high-frequency 1-50-1-50 μ s (energized time 164 μ s) sequence appended to a low-frequency 250-10-250-10 μ s (energized time 1 ms) sequence. This concatenated waveform allows for an impedance spectra sweep of 1.8 kHz–4.93 MHz which is critical in the realm of tissue electroporation-based applications. Low-frequency pulses are primarily confined to the extracellular regions. Once pores begin to form in the cellular membranes, more current pathways are formed, lowering the effective impedance until the tissue is fully electroporated. At higher frequencies, currents short the membrane reactance and can penetrate the cell more easily; subsequently, impedance changes at these high frequencies are minimally affected by pore formation. However, in the case of this specific study, no living cells or tissue are used, therefore, considerations of electroporation effects are not necessitated.
The model is based on the development of a state-space model (SSM) that considers tissue as a system with inputs (applied voltage, tissue conductivity, and impedance changes) and outputs (tissue temperature at 4 mm and 10 mm). The state-space model would be easy to implement for estimating temperature from real-time impedance measurements during treatments while keeping all of the crucial aspects of more complicated models. Hence, it may be used for real-time assessment of tissue response or as a theoretical foundation for more sophisticated temperature controls. In this work, the data-driven model was built from experimental results in an agar tissue phantom and then validated using established theoretical models based on numerical modeling.

2. Materials and Methods

2.1. AGAR Model

The experimental design was based on a tissue-mimic agar phantom (Figure 1). Four groups were evaluated to create a comprehensive model: high-voltage high-conductivity (HV-HC), low-voltage low-conductivity (LV-LC), high-voltage low-conductivity (HV-LC), and low-voltage high-conductivity (LV-HC). Data acquisition of n = 5 for each group was conducted (N = 20). A separate agar phantom was created for each sample. Deionized water combined with 1% agar (w/v) and 0.034% or 0.24% NaCl (w/v), to achieve 0.1220 or 0.526 S/m conductivities, respectively, were mixed and heated to 90 ° C. Mixed agar solution was poured into a plastic cup mold with a top diameter of 71 mm, a bottom diameter of 46 mm, and a height of 66 mm and allowed to cool to room temperature. Conductivities of the agar were measured at 37 ° C to ensure accurate conductive properties of tissues at physiological temperatures. The low conductivity model was set to mimic the conductivity of brain tissue while the high conductivity model was set to mimic initial conductivity of pancreatic tissue at a characteristic frequency of 1.8 kHz [32,33].

2.2. Treatment Parameters

Biphasic pulsed electric fields were administered using a custom pulse generator (VoltMed Inc., Blacksburg, Virginia). Treatment voltages were applied at either 1250 V or 2500 V (achieving voltage to distance ratios of 834 and 1666 V/cm, respectively), and delivered at a rate of 1 burst per second for 300 s. A 5-5-5-5 μ s waveform with 100 μ s of energized time was delivered through a single insertion bipolar probe (diameter = 1.65 mm; electrode exposure = 7 mm; insulation = 8 mm) with voltage and current waveforms monitored with a WaveSurfer 3024 z oscilloscope (Teledyne LeCroy, Chestnut Ridge, NY, USA) equipped with a 1000 × high voltage probe (BTX Enhancer 3000, Holliston, MA, USA) and a 10 × current probe (3972, Pearson Electronics, Palo Alto, CA, USA).

2.3. Real-Time Impedance Acquisition

FAST impedance measurements were recorded using a 3024z Oscilloscope (Teledyne LeCroy, Chesnut Ridge, NY, USA) and a 1 × current probe (2877, Pearson Electronics, Palo Alto, CA, USA). The recorded voltage and current waveforms from the custom FAST waveform (Figure 2A) were converted into the frequency domain (Figure 2B), and the complex tissue impedance was calculated as a ratio of the two following Ohm’s law. For this study, impedance measurements were calculated from peak values at the 1.8 kHz frequency as the agar phantom is unaffected by tissue electroporation effects. The low-frequency impedance decreases as a function of temperature increase. After a considerable number of pulses, the low-frequency impedance mimics the stable high-frequency impedance. Further details on this technique can be found at [31].

2.4. Temperature Acquisition

To achieve reproducible conditions between electroporation treatments, a custom-made electrode and temperature sensor holder was developed to precisely position the small fiber-optic temperature sensors (Luxtron m3300; LumaSense, Santa Clara, CA, USA) and single bipolar electrode. The holder was printed on a Form3 printer (FormLabs, USA) with a 25 μ m print resolution using FormLabs Clear Resin. The holder consists of a 2.1 mm canal in the center for the electrode, two 1.1 mm canals for the fiber-optic temperature sensors with a center to electrode surface distance of 4 mm and 10 mm, and a depth of 30 mm, placing the temperature probes at a depth parallel with the active top electrode. Prior to treatment, a slice parallel to and at a 10 mm distance from the electrode surface was made to create a flat surface for generation of an immediate post-treatment thermal profile, recorded on a FLIR A325SC thermal camera (FLIR, Wilsonville, OR, USA). Shortly after treatment (within ∼5 s), the agar was sliced along the electrode surface to obtain a thermal profile of the highest temperature areas.

2.5. Data Processing

Voltage and current processing and impedance extractions were conducted in MATLAB v.R2022a (MathWorks Inc., Natick, MA, USA). Impedance measurements were transmitted to GraphPad Prism (GraphPad Software, San Diego, CA, USA) for fitting to a one-phase decay line to remove noise, and then sent back to MATLAB for difference calculations. Similarly, collected temperature and impedance measurements were loaded into GraphPad Prism for quadratic line fitting to remove noise and then imported back into MATLAB for calculation of temperature changes. All further data processing and statistical analysis were conducted in GraphPad Prism.

2.6. Development of the Mathematical Model

A state-space model was chosen because of its capacity to decompose a complex higher-order differential equation into a sequence of first-order equations that can be easily solved, hence reducing computational burden. State-space control utilizes differential equations describing the time domain of the system using state variables in vector form. This makes it easy to evaluate the system using simple matrix algebra, which also enables the evaluation of complex multi-input multi-output (MIMO) systems. Other control system approaches require complicated Laplace transforms and Fourier transforms to transfer the system’s time domain representation—supplied as a complex set of differential equations—into the frequency domain—given as a collection of algebraic equations.
The agar tissue mimic was seen as a dynamic system with the applied voltage (V), tissue identification by electrical conductivity (S) at the applied characteristic frequency, and intra-treatment impedance change measurements ( Δ Z ) as the input signals, and the temperatures measured at depths of 4 mm ( T 4 ) and 10 mm ( T 10 ) as the output signals (Figure 3). State-space modeling only requires unique state equations for values influencing the system that cannot be deduced from other model states already included; thus, a separate state for applied waveform is not necessary as its influence can be deduced from the tissue conductivity.
We consider a system governed by the following continuous-time identified state-space model:
x ˙ = A x ( t ) + B u ( t ) + K e ( t )
x ( 0 ) = 0
y = C x ( t ) + D u ( t ) + e ( t )
where x is the state vector, u is the input vector, y is the output vector, and K is the disturbance component. A, B, C and D are the state-space matrices which describe the system dynamics. By default, D is set to 0 for dynamic systems, indicating that the system has no feedthrough.
The state-space model was estimated by utilizing the numerical algorithms for subspace state-space identification (N4SID). Implementation of this method in MATLAB’s System Identification toolbox uses input–output data to construct the model’s weighting matrices. This method can be viewed as a type of multi-step, constrained prediction error method utilizing linear regression to solve the system matrix. The algorithm was configured to compute a fourth order SSM. Further details on the methodology for the N4SID method are outlined in Appendix A.

2.7. Validation of the Mathematical Model against FEM

After developing an SSM based on results from experimentally collected data, a finite element model was developed to create new test data to assess the validity of the state-space model. A three-dimensional Comsol Multiphysics 6.0 (Comsol, Stockholm, Sweden) model was designed to replicate the experimental conditions of the agar phantom.
The distribution of electrical fields and thermal effects were computed using standard methods [34,35]. The electric potential distribution at the end of a given pulse was calculated with a modified Laplace equation under the electroquasistatic approximation (Equation (4)), and the resulting electric field distribution was calculated with Equation (5).
· ( σ ( T ) Φ ) = 0
E = Φ
where σ is the electrical conductivity, E is the local field magnitude, T is the temperature, and Φ is the local electric potential.
Following the computation of the electric field distribution, tissue temperature was computed using the general heat conduction equation with the addition of a Joule heating term ( Q J ),
ρ c p T t = · ( k T ) + Q J
Q J = σ | ϕ | 2 · p τ
where ρ is the density of the medium; c p is the specific heat; k is the thermal conductivity; p is the burst on-time ( 100 × 10 6 s ) and τ is the period of the burst delivery (1 s). The initial conditions were set to match those of the experimental conditions with the outer boundaries of the geometry being assigned a convective heat flux due to air convective cooling with an external temperature of T e x t = 18.8 ° C. Due to the discrepancies between agar conductivity being measured at physiological temperatures (37 ° C), set conductivities ( σ 0 ) were scaled by the temperature of the agar at the time of pulsing ( T a g a r ) to match the conductivity of the agar at the time of pulsing ( σ f ) by using the following:
σ f ( T ) = σ 0 ( 1 + α ( T a g a r 37 ° C )
producing initial conductivities of 0.331 S/m, 0.085 S/m, 0.093 S/m, and 0.331 S/m for the HV-HC, LV-LC, HV-LC and LV-HC groups, respectively. Material properties of the agar phantom were set to those of water. All other thermal and electrical conditions may be found in Table 1.

3. Results

3.1. Experimental Data Results

The recorded temperature and impedance change data points were used as training data along with known parameters such as voltage and tissue conductivity (Figure 4).

3.2. Model Validation with Training and Validation Data Sets

Selection of the final model was chosen through analysis of the accuracy of the training data being tested through the model, as well as a validation set in which specific data sets from each group were omitted from training the model.
The training and validation data set accuracies are shown in the heat map tiles (Figure 5) for T 4 and T 10 . Specification of the grouping where a data set was omitted for use as the validation data set is specified by the column titles. The fit accuracy was calculated by:
f i t = 100 ( 1 | y p r e d y | | y m e a n ( y ) | )
where y is the validation data point value and y p r e d is the predicted model from the output.

3.3. Model Implementation with Test Data

We validated our data-driven state-space model against a COMSOL finite element model, replicating our experimental setup. The calculated outputs of the COMSOL model were impedance and temperature at 4mm and 10mm for each voltage and conductivity condition. By importing the extracted impedance into MATLAB, we were then able to obtain predicted temperatures with our state-space model, trained on experimental data. These predicted temperatures were then used to calculate fit accuracy and maximum absolute error (Table 2).
We found that the mathematical model provided excellent estimations of temperature rise, with above 90% fit accuracy for all conditions, except for the LV-LC condition. The discrepancy for this point can be attributed to the low absolute change in temperature. At the 10 mm distance, the temperature increase is negligible, and the model was not able to fit such a small change. Though the LV-LC had the lowest fit accuracy, the absolute error in temperature prediction was under 0.5 ° C. Further, all the test conditions had a maximum absolute error of under 0.5 ° C, except for the HV-HC, which has a maximum absolute error of 1.45 and 0.52 ° C for the 4 mm and 10 mm positions, respectively. The state-space model temperature predictions and finite element temperature simulations at both distances are given in Figure 6 for each voltage and conductivity condition.
The overall model quality was found to be 9.116e-10 as calculated by Akaike’s Final Prediction Error (FPE):
F P E = V ( 1 + 2 d N )
where d is the number of estimated parameters and the loss function, V, is represented by:
V = d e t ( 1 N 1 N ϵ ( t , θ N ) ( ϵ ( t , θ N ) ) T )
where N represents the number of values in the estimation data set, ϵ ( t ) is the respective prediction error, and ( θ ) N is the estimated parameters.

4. Discussion

The goal of this study was to evaluate the feasibility of real-time temperature monitoring during an IRE or H-FIRE procedure without external temperature probes or devices. This feasibility study demonstrates a proof-of-concept construction of a black-box model to forecast tissue temperature rise during the ablation process based on observations of real-time impedance changes. In addition to not requiring additional sensors, one of the greatest potential advantages to implementing a state-space model is its ability to provide spatial data at several locations which would otherwise be challenging, if not impossible, in vivo. The model yielded a reasonably good estimate of temperature rise in most cases; however, in certain instances, considerable inaccuracies were detected (e.g., those in which the total treatment temperature rise was on average < 1.5   ° C). Despite these discrepancies at low temperatures, temperature rises of this magnitude are not relevant in the context of clinical temperature monitoring. From Figure 7, it is evident that the conductivity of the agar plays the greatest role in temperature fluctuations and thus suggests that thermal monitoring may only be necessary in the treatment of tissues with high conductivities (e.g., pancreas, prostate, etc.).
Due to the fact that the electrical and thermal properties of the gel may differ from those of normal tissue—for instance, the gel was cooler and was not perfused—the absence of our model accounting for flow-rate perfusion, which in the clinical setting should permit a considerable degree of variability, means that these temperature curves cannot be directly applied to the clinical environment. Even though we have focused on the development of wholly data-driven prediction models, future computational modeling could be used to incorporate more biophysical information in the model’s input parameters as a means for expanding the utility of the model without making it entirely data-driven. This strategy could improve the accuracy of predictions, as we believe the proposed technique could be enhanced by combining new data from a range of electrical and thermal tissue properties. Moreover, state-space models built on computationally driven data could make it simpler to create multi-output models with more locations for temperature predictions. Future studies may also consider developing a model which considers a multi-electrode setup. However, such geometries would introduce additional variables to be fit including electrode spacing and exposure.
In addition, the process could be more complicated in the tissue domain due to the complex impedance spectrum in living tissue [38,39,40,41,42]. Because low-frequency impedance measurements are primarily confined to the extracellular region prior to tissue electroporation, membrane permeabilization during treatment has a substantial effect on impedance alterations. However, high-frequency impedance measurements, which correspond to currents that short the membrane reactance and penetrate the cell membrane, are less vulnerable to membrane pore formation. Thus, in the case of a tissue model, high-frequency impedance measurements may be used as a benchmark to identify impedance changes caused primarily by thermal effects.

5. Conclusions

State-space models can be created to predict tissue temperature increases during H-FIRE treatment by utilizing impedance changes calculated from rapid impedance spectroscopy, applied voltage, and the conductivity of the tissue as described at the characteristic frequency of the applied waveform. The best model produces a reasonably accurate prediction of tissue temperature with an overall model output to test data fit of 95.8% and 93.22% for T 4 and T 10 , respectively. The largest model discrepancies are observed in the scenario where the least amount of temperature rise is observed overall (temperature increases < 1.5   ° C). Test data validation of the model suggests that the model is successful at reasonably predicting temperature increases in cases where temperature increases are of clinical relevance.

Author Contributions

Conceptualization, S.N.C. and R.V.D.; Data curation, S.N.C., E.J.J.IV and K.N.A.; Formal analysis, S.N.C. and K.N.A.; Funding acquisition, R.V.D.; Investigation, S.N.C., E.J.J.IV, K.N.A. and R.V.D.; Methodology, S.N.C., E.J.J.IV, K.N.A. and R.V.D.; Project administration, S.N.C. and R.V.D.; Resources, S.N.C. and R.V.D.; Software, S.N.C., E.J.J.IV and K.N.A.; Supervision, S.N.C. and R.V.D.; Validation, S.N.C., E.J.J.IV and K.N.A.; Visualization, S.N.C., E.J.J.IV, K.N.A. and R.V.D.; Writing—original draft, S.N.C., E.J.J.IV, K.N.A. and R.V.D.; Writing—review and editing, S.N.C., E.J.J.IV, K.N.A. and R.V.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was in part funded by the National Institutes of Health/National Cancer Institute R01 CA240476. Sabrina Campelo was further funded by the ICTAS Doctoral Fellowship at Virginia Tech and the P.E.O Scholar Award.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Acknowledgments

The authors would like to acknowledge support from the Institute of Critical Technology and Sciences (ICTAS) and its Center for Engineered Health (CEH) at Virginia Tech. Additionally the authors thank Melvin Lorenzo for training on FAST technologies. The authors also would like to thank Zaid Salameh and Kailee David for assistance with reviewing and editing the manuscript.

Conflicts of Interest

Kenneth N. Aycock and Rafael V. Davalos have issued patents and/or patents pending in the area of irreversible electroporation and FAST technology and may receive royalties.

Abbreviations

The following abbreviations are used in this manuscript:
FASTFourier Analysis SpecTroscopy
FEMFinite Element Modeling
H-FIREHigh-Frequency Irreversible Electroporation
IREIrreversible Electroporation
N4SIDNumerical Algorithms for Subspace State-Space Identification
PEFPulsed Electric Field
SSMState-Space Model

Appendix A

N4sid Oblique Projection by LQ Decomposition

The N4SID subspace method uses LQ decomposition, singular value decomposition, and applies an oblique projection approach for input and output Hankel blocks to produce subspace weighting matrices [43,44]. Past and future block Hankel matrices for inputs U p , U f and outputs Y p , Y f are constructed from experimental input and output data,
U p = u 0 u 1 u N 1 u 1 u 2 u N u 2 k 1 u 2 k u 2 k + N 2 U f = u k u k + 1 u k + N 1 u k + 1 u k + 2 u k + N u k 1 u k u k + N 2
Y p = y 0 y 1 y N 1 y 1 y 2 y N y 2 k 1 y 2 k y 2 k + N 2 Y f = y k y k + 1 y k + N 1 y k + 1 y k + 2 y k + N y k 1 y k y k + N 2
in which N gives us the dimensions of the matrix size and k gives us the number of elements in each matrix.
Our output state-space equation can be defined as
y ( t ) = C x ( t ) + D u ( t )
The N4SID subspace method determines that our output matrix, Y k ( t ) , can thus be found through the following relationship:
y ( t ) y ( t + 1 ) y ( t + k 1 ) Y k ( t ) = C C A C A k 1 O k x ( t ) X 0 + D 0 0 C B D 0 C A k 2 B C B D Ψ k u ( t ) u ( t + 1 ) u ( t + k 1 ) U k ( t )
with O k representing the observability matrix, X 0 representing the state initial state, Ψ k representing the Toeplitz matrix designating the shift-invariant properties of the system, and U k is the input matrix. Thus, the output matrix can be represented by
Y f = O k X f + Ψ k U f .
In order to solve for the solution to our subspace weighting factors [A B C D], we must perform LQ decomposition of the data matrix ( W p ) into the product of a lower triangular matrix (L) and a unitary matrix (Q):
W p = U p Y p
U f W p Y f = R 11 0 0 R 21 R 22 0 R 31 R 32 0 Q 1 T Q 2 T Q 3 T
The resulting equations are as follows:
U f = R 11 Q 1 T
W p = R 21 Q 1 T + R 22 Q 2 T
Y f = R 31 Q 1 T + R 32 Q 2 T
Through substitution of terms, we end up with the following relationship:
Y f = ( R 31 R 32 R 22 # R 21 ) R 11 1 U f + R 32 R 22 # W p
suggesting that the two terms U f and W p exist in two different sub-spaces with no overlapping bases. Thus from this realization, we can conclude that from Equation (A2),
O k X f = R 32 R 22 # W p
and use singular value decomposition of R 32 R 22 # W p to solve for O k and X f .
R 32 R 22 # W p = U 1 U 2 Σ 1 0 0 0 V 1 T V 2 T = U 1 Σ 1 V 1 T
The output of Equation (A10) can therefore be split between O k and X f :
O k = U 1 Σ 1 1 / 2 T
X f = T 1 Σ 1 1 / 2 V 1 T
By further decomposing X f , we can define the series of states:
x ( k ) , x ( k + 1 ) , x ( k + N 1 )
By combining the states and input–output data, we end up with the following matrices:
X ¯ k = ( x ( k ) x ( k + N 2 ) )
X ¯ k + 1 = ( x ( k + 1 ) x ( k + N 1 ) )
U ¯ k | k = ( u ( k ) u ( k + N 2 ) )
Y ¯ k | k = ( y ( k ) y ( k + N 2 ) )
and presented in the following state and measurement equations:
X ¯ k + 1 Y ¯ k | k = A B C D X ¯ k U ¯ k | k
through which our weighting matrices [A B C D] can be solved through a least squares regression:
A ^ B ^ C ^ D ^ = X ¯ k + 1 Y ¯ k | k X ¯ k U ¯ k | k T X ¯ k U ¯ k | k X ¯ k U ¯ k | k T 1
For our specific system, the following weighting matrices were produced:
A= x 1 x 2 x 3 x 4 x 1 0.002906 0.0004279 0.001115 0.001509 x 2 70.003136 0.003182 0.008007 0.0006536 x 3 0.002613 0.005233 0.006247 0.0004292 x 4 0.002367 0.00598 0.001739 0.00521
B= d Z V S x 1 2.941 e 07 4.348 e 08 8.484 e 06 x 2 1.581 e 07 3.313 e 08 5.155 e 05 x 3 5.258 e 07 3.26 e 08 3.914 e 05 x 4 2.326 e 06 1.432 e 07 8.647 e 06
C= x 1 x 2 x 3 x 4 T 4 m m 238.5 74.14 0.6752 0.2305 T 10 m m 98.88 6.118 0.1381 0.1084
D= d Z V S T 4 m m 0 0 0 T 10 m m 0 0 0
K= T 4 m m T 10 m m x 1 6.139 e 05 0.002894 x 2 0.003739 0.009211 x 3 0.1299 0.2822 x 4 0.1849 0.7076

References

  1. Geboers, B.; Scheffer, H.J.; Graybill, P.M.; Ruarus, A.H.; Nieuwenhuizen, S.; Puijk, R.S.; van den Tol, P.M.; Davalos, R.V.; Rubinsky, B.; de Gruijl, T.D.; et al. High-voltage electrical pulses in oncology: Irreversible electroporation, electrochemotherapy, gene electrotransfer, electrofusion, and electroimmunotherapy. Radiology 2020, 295, 254–272. [Google Scholar] [CrossRef]
  2. Maor, E.; Sugrue, A.; Witt, C.; Vaidya, V.R.; DeSimone, C.V.; Asirvatham, S.J.; Kapa, S. Pulsed electric fields for cardiac ablation and beyond: A state-of-the-art review. Heart Rhythm 2019, 16, 1112–1120. [Google Scholar] [CrossRef] [PubMed]
  3. Wittkampf, F.H.; van Es, R.; Neven, K. Electroporation and its relevance for cardiac catheter ablation. JACC Clin. Electrophysiol. 2018, 4, 977–986. [Google Scholar] [CrossRef] [PubMed]
  4. Weaver, J.C.; Smith, K.C.; Esser, A.T.; Son, R.S.; Gowrishankar, T. A brief overview of electroporation pulse strength–duration space: A region where additional intracellular effects are expected. Bioelectrochemistry 2012, 87, 236–243. [Google Scholar] [CrossRef]
  5. Yarmush, M.L.; Golberg, A.; Serša, G.; Kotnik, T.; Miklavčič, D. Electroporation-based technologies for medicine: Principles, applications, and challenges. Annu. Rev. Biomed. Eng. 2014, 16, 295–320. [Google Scholar] [CrossRef]
  6. Weaver, J.C. Electroporation: A general phenomenon for manipulating cells and tissues. J. Cell. Biochem. 1993, 51, 426–435. [Google Scholar] [CrossRef]
  7. Mercadal, B.; Beitel-White, N.; Aycock, K.N.; Castellví, Q.; Davalos, R.V.; Ivorra, A. Dynamics of cell death after conventional IRE and H-FIRE treatments. Ann. Biomed. Eng. 2020, 48, 1451–1462. [Google Scholar] [CrossRef]
  8. Napotnik, T.B.; Polajžer, T.; Miklavčič, D. Cell death due to electroporation–a review. Bioelectrochemistry 2021, 141, 107871. [Google Scholar] [CrossRef] [PubMed]
  9. Davalos, R.V.; Mir, L.; Rubinsky, B. Tissue ablation with irreversible electroporation. Ann. Biomed. Eng. 2005, 33, 223–231. [Google Scholar] [CrossRef] [PubMed]
  10. Faroja, M.; Ahmed, M.; Appelbaum, L.; Ben-David, E.; Moussa, M.; Sosna, J.; Nissenbaum, I.; Goldberg, S.N. Irreversible electroporation ablation: Is all the damage nonthermal? Radiology 2013, 266, 462–470. [Google Scholar] [CrossRef]
  11. Van Den Bos, W.; Scheffer, H.J.; Vogel, J.A.; Wagstaff, P.G.; De Bruin, D.M.; De Jong, M.C.; Van Gemert, M.J.; De La Rosette, J.J.; Meijerink, M.R.; Klaessens, J.H.; et al. Thermal energy during irreversible electroporation and the influence of different ablation parameters. J. Vasc. Interv. Radiol. 2016, 27, 433–443. [Google Scholar] [CrossRef] [PubMed]
  12. Yang, Y.; Moser, M.; Zhang, E.; Zhang, W.; Zhang, B. Optimization of electrode configuration and pulse strength in irreversible electroporation for large ablation volumes without thermal damage. J. Eng. Sci. Med. Diagn. Ther. 2018, 1, 021002. [Google Scholar] [CrossRef]
  13. Martin, R.C.; Kwon, D.; Chalikonda, S.; Sellers, M.; Kotz, E.; Scoggins, C.; McMasters, K.M.; Watkins, K. Treatment of 200 locally advanced (stage III) pancreatic adenocarcinoma patients with irreversible electroporation: Safety and efficacy. Ann. Surg. 2015, 262, 486–494. [Google Scholar] [CrossRef]
  14. Philips, P.; Hays, D.; Martin, R.C. Irreversible electroporation ablation (IRE) of unresectable soft tissue tumors: Learning curve evaluation in the first 150 patients treated. PLoS ONE 2013, 8, e76260. [Google Scholar] [CrossRef] [PubMed]
  15. Ansari, D.; Kristoffersson, S.; Andersson, R.; Bergenfeldt, M. The role of irreversible electroporation (IRE) for locally advanced pancreatic cancer: A systematic review of safety and efficacy. Scand. J. Gastroenterol. 2017, 52, 1165–1171. [Google Scholar] [CrossRef] [PubMed]
  16. Silk, M.T.; Wimmer, T.; Lee, K.S.; Srimathveeravalli, G.; Brown, K.T.; Kingham, P.T.; Fong, Y.; Durack, J.C.; Sofocleous, C.T.; Solomon, S.B. Percutaneous ablation of peribiliary tumors with irreversible electroporation. J. Vasc. Interv. Radiol. 2014, 25, 112–118. [Google Scholar] [CrossRef]
  17. Dunki-Jacobs, E.; Philips, P.; Martin Ii, R. Evaluation of thermal injury to liver, pancreas and kidney during irreversible electroporation in an in vivo experimental model. J. Br. Surg. 2014, 101, 1113–1121. [Google Scholar] [CrossRef]
  18. Kurata, K.; Nomura, S.; Takamatsu, H. Three-dimensional analysis of irreversible electroporation: Estimation of thermal and non-thermal damage. Int. J. Heat Mass Transf. 2014, 72, 66–74. [Google Scholar] [CrossRef]
  19. Wagstaff, P.G.; de Bruin, D.M.; van den Bos, W.; Ingels, A.; van Gemert, M.J.; Zondervan, P.J.; Verdaasdonk, R.M.; van Lienden, K.P.; van Leeuwen, T.G.; de la Rosette, J.J.; et al. Irreversible electroporation of the porcine kidney: Temperature development and distribution. In Urologic Oncology: Seminars and Original Investigations; Elsevier: Amsterdam, The Netherlands, 2015; Volume 33, pp. 168.e1–168.e7. [Google Scholar]
  20. Agnass, P.; van Veldhuisen, E.; van Gemert, M.J.; van der Geld, C.W.; van Lienden, K.P.; van Gulik, T.M.; Meijerink, M.R.; Besselink, M.G.; Kok, H.P.; Crezee, J. Mathematical modeling of the thermal effects of irreversible electroporation for in vitro, in vivo, and clinical use: A systematic review. Int. J. Hyperth. 2020, 37, 486–505. [Google Scholar] [CrossRef] [PubMed]
  21. Scheffer, H.J.; Vogel, J.A.; Van Den Bos, W.; Neal, R.E.; Van Lienden, K.P.; Besselink, M.G.; Van Gemert, M.J.; Van Der Geld, C.W.; Meijerink, M.R.; Klaessens, J.H.; et al. The influence of a metal stent on the distribution of thermal energy during irreversible electroporation. PLoS ONE 2016, 11, e0148457. [Google Scholar] [CrossRef] [PubMed]
  22. Brock, R.M.; Beitel-White, N.; Davalos, R.V.; Allen, I.C. Starting a fire without flame: The induction of cell death and inflammation in electroporation-based tumor ablation strategies. Front. Oncol. 2020, 10, 1235. [Google Scholar] [CrossRef]
  23. Zhao, J.; Wen, X.; Tian, L.; Li, T.; Xu, C.; Wen, X.; Melancon, M.P.; Gupta, S.; Shen, B.; Peng, W.; et al. Irreversible electroporation reverses resistance to immune checkpoint blockade in pancreatic cancer. Nat. Commun. 2019, 10, 899. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, N.; Li, Z.; Han, X.; Zhu, Z.; Zhao, Y.; Liu, Z.; Lv, Y. Irreversible Electroporation: An Emerging Immunomodulatory Therapy on Solid Tumor. Front. Immunol. 2022, 12, 811726. [Google Scholar] [CrossRef] [PubMed]
  25. Arena, C.B.; Mahajan, R.L.; Nichole Rylander, M.; Davalos, R.V. An experimental and numerical investigation of phase change electrodes for therapeutic irreversible electroporation. J. Biomech. Eng. 2013, 135, 111009. [Google Scholar] [CrossRef] [PubMed]
  26. O’Brien, T.J.; Bonakdar, M.; Bhonsle, S.; Neal, R.E.; Aardema Jr, C.H.; Robertson, J.L.; Goldberg, S.N.; Davalos, R.V. Effects of internal electrode cooling on irreversible electroporation using a perfused organ model. Int. J. Hyperth. 2018, 35, 44–55. [Google Scholar] [CrossRef] [PubMed]
  27. O’Brien, T.J.; Lorenzo, M.F.; Zhao, Y.; Neal II, R.E.; Robertson, J.L.; Goldberg, S.N.; Davalos, R.V. Cycled pulsing to mitigate thermal damage for multi-electrode irreversible electroporation therapy. Int. J. Hyperth. 2019, 36, 952–962. [Google Scholar] [CrossRef] [PubMed]
  28. Kaufman, J.D.; Fesmire, C.C.; Petrella, R.A.; Fogle, C.A.; Xing, L.; Gerber, D.; Sano, M.B. High-frequency irreversible electroporation using 5000-V waveforms to create reproducible 2-and 4-cm ablation zones—A laboratory investigation using mechanically perfused liver. J. Vasc. Interv. Radiol. 2020, 31, 162–168. [Google Scholar] [CrossRef]
  29. Aycock, K.N.; Campelo, S.N.; Davalos, R.V. A Comparative Modeling Study of Thermal Mitigation Strategies in Irreversible Electroporation Treatments. J. Heat Transf. 2022, 144, 031206. [Google Scholar] [CrossRef] [PubMed]
  30. Cindrič, H.; Mariappan, P.; Beyer, L.; Wiggermann, P.; Moche, M.; Miklavčič, D.; Kos, B. Retrospective study for validation and improvement of numerical treatment planning of irreversible electroporation ablation for treatment of liver tumors. IEEE Trans. Biomed. Eng. 2021, 68, 3513–3524. [Google Scholar] [CrossRef] [PubMed]
  31. Lorenzo, M.F.; Bhonsle, S.P.; Arena, C.B.; Davalos, R.V. Rapid impedance spectroscopy for monitoring tissue impedance, temperature, and treatment outcome during electroporation-based therapies. IEEE Trans. Biomed. Eng. 2020, 68, 1536–1546. [Google Scholar] [CrossRef] [PubMed]
  32. Gabriel, S.; Lau, R.; Gabriel, C. The dielectric properties of biological tissues: II. Measurements in the frequency range 10 Hz to 20 GHz. Phys. Med. Biol. 1996, 41, 2251. [Google Scholar] [CrossRef] [PubMed]
  33. IT’IS Database for thermal and electromagnetic parameters of biological tissues, Version 4.0. Available online: https://itis.swiss/virtual-population/tissue-properties/downloads/database-v4-1/ (accessed on 15 August 2022).
  34. Garcia, P.A.; Davalos, R.V.; Miklavcic, D. A numerical investigation of the electric and thermal cell kill distributions in electroporation-based therapies in tissue. PLoS ONE 2014, 9, e103083. [Google Scholar] [CrossRef]
  35. Mohammadi, A.; Bianchi, L.; Asadi, S.; Saccomandi, P. Measurement of ex vivo liver, brain and pancreas thermal properties as function of temperature. Sensors 2021, 21, 4236. [Google Scholar] [CrossRef] [PubMed]
  36. Popiel, C.; Wojtkowiak, J. Simple formulas for thermophysical properties of liquid water for heat transfer calculations (from 0 C to 150 C). Heat Transf. Eng. 1998, 19, 87–101. [Google Scholar] [CrossRef]
  37. Ngo, I.L.; Jeon, S.; Byon, C. Thermal conductivity of transparent and flexible polymers containing fillers: A literature review. Int. J. Heat Mass Transf. 2016, 98, 219–226. [Google Scholar] [CrossRef]
  38. Castellví, Q.; Mercadal, B.; Ivorra, A. Assessment of electroporation by electrical impedance methods. Handb. Electroporation 2016, 67, 1–20. [Google Scholar]
  39. Ivorra, A.; Rubinsky, B. In vivo electrical impedance measurements during and after electroporation of rat liver. Bioelectrochemistry 2007, 70, 287–295. [Google Scholar] [CrossRef] [PubMed]
  40. Dean, D.; Ramanathan, T.; Machado, D.; Sundararajan, R. Electrical impedance spectroscopy study of biological tissues. J. Electrost. 2008, 66, 165–177. [Google Scholar] [CrossRef] [PubMed]
  41. Schwan, H.P. Electrical properties of tissue and cell suspensions. In Advances in Biological and Medical Physics; Elsevier: Amsterdam, The Netherlands, 1957; Volume 5, pp. 147–209. [Google Scholar]
  42. Voyer, D.; Silve, A.; Mir, L.M.; Scorretti, R.; Poignard, C. Dynamical modeling of tissue electroporation. Bioelectrochemistry 2018, 119, 98–110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Jamaludin, I.; Wahab, N.; Khalid, N.; Sahlan, S.; Ibrahim, Z.; Rahmat, M.F. N4SID and MOESP subspace identification methods. In Proceedings of the 2013 IEEE 9th International Colloquium on Signal Processing and its Applications, Kuala Lumpur, Malaysia, 8–10 March 2013; pp. 140–145. [Google Scholar]
  44. Asada, H. Subspace Methods for System Identification: Realization; Department of Mechanical Engineering MIT: Cambridge, MA, USA, 2020. [Google Scholar]
Figure 1. Schematic of the physical experimental setup for recording temperature and impedance measurements during treatment delivered from a high-voltage (HV) generator. A custom 3D printed electrode and fiber-optic temperature holder was fit on top of the agar tissue phantom. Two fiber-optic temperature probes were fed into the canals at a distance of 4 and 10 mm away from the electrode surface. Impedance measurements were calculated from voltage and current measurements collected on the oscilloscope. A thermal camera was placed at a distance 25 cm from the flat edge surface of the tissue phantom.
Figure 1. Schematic of the physical experimental setup for recording temperature and impedance measurements during treatment delivered from a high-voltage (HV) generator. A custom 3D printed electrode and fiber-optic temperature holder was fit on top of the agar tissue phantom. Two fiber-optic temperature probes were fed into the canals at a distance of 4 and 10 mm away from the electrode surface. Impedance measurements were calculated from voltage and current measurements collected on the oscilloscope. A thermal camera was placed at a distance 25 cm from the flat edge surface of the tissue phantom.
Bioengineering 09 00499 g001
Figure 2. (A) The custom diagnostic FAST acquisition waveform delivered between therapeutic voltage bursts creates (B) a normalized frequency spectrum of the absolute FFT voltage magnitude, which was then used to collect real-time impedance measurements from our agar phantom. Delivered waveforms for the respective low-voltage (C) and high-voltage (D) groups are also shown.
Figure 2. (A) The custom diagnostic FAST acquisition waveform delivered between therapeutic voltage bursts creates (B) a normalized frequency spectrum of the absolute FFT voltage magnitude, which was then used to collect real-time impedance measurements from our agar phantom. Delivered waveforms for the respective low-voltage (C) and high-voltage (D) groups are also shown.
Bioengineering 09 00499 g002
Figure 3. Schematic proposing the model inputs of change in impedance ( Δ Z), applied voltage (V), and initial tissue conductivity under the influence of a specific frequency waveform ( σ T i s s u e , f ) to be fed into a fit state-space model to give real-time estimation of the temperature rise at distances 4 mm and 10 mm away from the electrode surface.
Figure 3. Schematic proposing the model inputs of change in impedance ( Δ Z), applied voltage (V), and initial tissue conductivity under the influence of a specific frequency waveform ( σ T i s s u e , f ) to be fed into a fit state-space model to give real-time estimation of the temperature rise at distances 4 mm and 10 mm away from the electrode surface.
Bioengineering 09 00499 g003
Figure 4. Results of the experimental temperature changes for the (A) High Voltage-High Conductivity (HV-HC), (B) Low Voltage-Low Conductivity (LV-LC), (C) High Voltage-Low Conductivity (HV-LC), and (D) Low Voltage-High Conductivity (LV-HC) groups as well as the (E) recorded changes in impedance measurements from the agar tissue phantom. Shaded regions relate to standard deviations within each group.
Figure 4. Results of the experimental temperature changes for the (A) High Voltage-High Conductivity (HV-HC), (B) Low Voltage-Low Conductivity (LV-LC), (C) High Voltage-Low Conductivity (HV-LC), and (D) Low Voltage-High Conductivity (LV-HC) groups as well as the (E) recorded changes in impedance measurements from the agar tissue phantom. Shaded regions relate to standard deviations within each group.
Bioengineering 09 00499 g004
Figure 5. A heat map distribution of the percent accuracies for the respective High Voltage-High Conductivity (HV-HC), Low Voltage-High Conductivity (LV-HC), High Voltage-Low Conductivity (HV-LC), and Low Voltage-Low Conductivity (LV-LC) training and validation sets during the training phase of the model formation.
Figure 5. A heat map distribution of the percent accuracies for the respective High Voltage-High Conductivity (HV-HC), Low Voltage-High Conductivity (LV-HC), High Voltage-Low Conductivity (HV-LC), and Low Voltage-Low Conductivity (LV-LC) training and validation sets during the training phase of the model formation.
Bioengineering 09 00499 g005
Figure 6. Impedance change measurements extracted from the FEM were used as test data fed through the resulting state-space model. Temperature rise from the FEM and the resulting outputs from the state-space model for the High Voltage-High Conductivity (HV-HC), Low Voltage-Low Conductivity (LV-LC), High Voltage-Low Conductivity(HV-LC), and Low Voltage-Low Conductivity (LV-LC) groups are plotted. Note: y-axis is scaled to view the largest absolute errors in model performance.
Figure 6. Impedance change measurements extracted from the FEM were used as test data fed through the resulting state-space model. Temperature rise from the FEM and the resulting outputs from the state-space model for the High Voltage-High Conductivity (HV-HC), Low Voltage-Low Conductivity (LV-LC), High Voltage-Low Conductivity(HV-LC), and Low Voltage-Low Conductivity (LV-LC) groups are plotted. Note: y-axis is scaled to view the largest absolute errors in model performance.
Bioengineering 09 00499 g006
Figure 7. Snapshots taken from the FLIR thermal camera immediately after pulsing at the phantom surface (a distance 10 mm away from the electrode), and shortly after pulsing (within ∼5 s) at the center of the phantom in the plane of the electrode surface.
Figure 7. Snapshots taken from the FLIR thermal camera immediately after pulsing at the phantom surface (a distance 10 mm away from the electrode), and shortly after pulsing (within ∼5 s) at the center of the phantom in the plane of the electrode surface.
Bioengineering 09 00499 g007
Table 1. Electrical and Thermal Properties for Computational Modeling.
Table 1. Electrical and Thermal Properties for Computational Modeling.
MaterialParameterValueUnitsRef
AgarDensity, ρ
Specific heat, cp
Thermal conductivity, k
Temperature coefficient, α
998
4182
0.598
2
kg/m3
J/(kg·K)
W/(m·K)
%/ ° C
[36]
[36]
[36]
[36]
InsulationDensity, ρ
Specific heat, cp
Thermal conductivity, k
Electrical conductivity, σ
2329
700
0.2
1 × 10−12
kg/m3
J/(kg·K)
W/(m·K)
S/m
[26]
[26]
[37]
[26]
Stainless SteelDensity, ρ
Specific heat, cp
Thermal conductivity, k
Electrical conductivity, σ
7900
500
15
2.22 × 106
kg/m3
J/(kg·K)
W/(m·K)
S/m
[26]
[26]
[26]
[26]
Table 2. Quantitative metrics on the validity of the model following feedthrough of FEM test data.
Table 2. Quantitative metrics on the validity of the model following feedthrough of FEM test data.
Temperature DistanceHV-HCLV-LCHV-LCLV-HC
Fit Accuracy (%)T 4 mm 95.888.490.793.7
T 10 mm 93.258.990.090.9
Maximum Absolute Error ( ° C)T 4 mm 1.450.060.300.19
T 10 mm 0.520.350.310.17
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Campelo, S.N.; Jacobs, E.J., IV; Aycock, K.N.; Davalos, R.V. Real-Time Temperature Rise Estimation during Irreversible Electroporation Treatment through State-Space Modeling. Bioengineering 2022, 9, 499. https://doi.org/10.3390/bioengineering9100499

AMA Style

Campelo SN, Jacobs EJ IV, Aycock KN, Davalos RV. Real-Time Temperature Rise Estimation during Irreversible Electroporation Treatment through State-Space Modeling. Bioengineering. 2022; 9(10):499. https://doi.org/10.3390/bioengineering9100499

Chicago/Turabian Style

Campelo, Sabrina N., Edward J. Jacobs, IV, Kenneth N. Aycock, and Rafael V. Davalos. 2022. "Real-Time Temperature Rise Estimation during Irreversible Electroporation Treatment through State-Space Modeling" Bioengineering 9, no. 10: 499. https://doi.org/10.3390/bioengineering9100499

APA Style

Campelo, S. N., Jacobs, E. J., IV, Aycock, K. N., & Davalos, R. V. (2022). Real-Time Temperature Rise Estimation during Irreversible Electroporation Treatment through State-Space Modeling. Bioengineering, 9(10), 499. https://doi.org/10.3390/bioengineering9100499

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop