1
|
Molinari L, Zaltieri M, Massaroni C, Filippi S, Gizzi A, Schena E. Multiscale and Multiphysics Modeling of Anisotropic Cardiac RFCA: Experimental-Based Model Calibration via Multi-Point Temperature Measurements. Front Physiol 2022; 13:845896. [PMID: 35514332 PMCID: PMC9062295 DOI: 10.3389/fphys.2022.845896] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2021] [Accepted: 03/22/2022] [Indexed: 11/13/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Leonardo Molinari
- Department of Mathematics and Computer Science, Emory University, Atlanta, GA, United States
| | - Martina Zaltieri
- Laboratory of Measurement and Biomedical Instrumentation, Department of Engineering, University of Rome Campus Bio-Medico, Rome, Italy
| | - Carlo Massaroni
- Laboratory of Measurement and Biomedical Instrumentation, Department of Engineering, University of Rome Campus Bio-Medico, Rome, Italy
| | - Simonetta Filippi
- Nonlinear Physics and Mathematical Modeling Lab, Department of Engineering, University of Rome Campus Bio-Medico, Rome, Italy
| | - Alessio Gizzi
- Nonlinear Physics and Mathematical Modeling Lab, Department of Engineering, University of Rome Campus Bio-Medico, Rome, Italy
| | - Emiliano Schena
- Laboratory of Measurement and Biomedical Instrumentation, Department of Engineering, University of Rome Campus Bio-Medico, Rome, Italy
| |
Collapse
|
2
|
Gryshkov O, Mutsenko V, Tarusin D, Khayyat D, Naujok O, Riabchenko E, Nemirovska Y, Danilov A, Petrenko AY, Glasmacher B. Coaxial Alginate Hydrogels: From Self-Assembled 3D Cellular Constructs to Long-Term Storage. Int J Mol Sci 2021; 22:3096. [PMID: 33803546 PMCID: PMC8003018 DOI: 10.3390/ijms22063096] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2021] [Accepted: 03/16/2021] [Indexed: 12/22/2022] Open
Abstract
Alginate as a versatile naturally occurring biomaterial has found widespread use in the biomedical field due to its unique features such as biocompatibility and biodegradability. The ability of its semipermeable hydrogels to provide a favourable microenvironment for clinically relevant cells made alginate encapsulation a leading technology for immunoisolation, 3D culture, cryopreservation as well as cell and drug delivery. The aim of this work is the evaluation of structural properties and swelling behaviour of the core-shell capsules for the encapsulation of multipotent stromal cells (MSCs), their 3D culture and cryopreservation using slow freezing. The cells were encapsulated in core-shell capsules using coaxial electrospraying, cultured for 35 days and cryopreserved. Cell viability, metabolic activity and cell-cell interactions were analysed. Cryopreservation of MSCs-laden core-shell capsules was performed according to parameters pre-selected on cell-free capsules. The results suggest that core-shell capsules produced from the low viscosity high-G alginate are superior to high-M ones in terms of stability during in vitro culture, as well as to solid beads in terms of promoting formation of viable self-assembled cellular structures and maintenance of MSCs functionality on a long-term basis. The application of 0.3 M sucrose demonstrated a beneficial effect on the integrity of capsules and viability of formed 3D cell assemblies, as compared to 10% dimethyl sulfoxide (DMSO) alone. The proposed workflow from the preparation of core-shell capsules with self-assembled cellular structures to the cryopreservation appears to be a promising strategy for their off-the-shelf availability.
Collapse
Affiliation(s)
- Oleksandr Gryshkov
- Institute for Multiphase Processes, Leibniz University Hannover, An der Universität 1, Building 8143, 30823 Garbsen, Germany; (V.M.); (D.K.); (B.G.)
- Lower Saxony Centre for Biomedical Engineering, Implant Research and Development, Stadtfelddamm 34, 30625 Hannover, Germany
| | - Vitalii Mutsenko
- Institute for Multiphase Processes, Leibniz University Hannover, An der Universität 1, Building 8143, 30823 Garbsen, Germany; (V.M.); (D.K.); (B.G.)
- Lower Saxony Centre for Biomedical Engineering, Implant Research and Development, Stadtfelddamm 34, 30625 Hannover, Germany
| | - Dmytro Tarusin
- Institute for Problems of Cryobiology and Cryomedicine of the National Academy of Sciences of Ukraine, 23 Pereyaslavsky Street, 61015 Kharkiv, Ukraine; (D.T.); (Y.N.); (A.Y.P.)
| | - Diaa Khayyat
- Institute for Multiphase Processes, Leibniz University Hannover, An der Universität 1, Building 8143, 30823 Garbsen, Germany; (V.M.); (D.K.); (B.G.)
- Lower Saxony Centre for Biomedical Engineering, Implant Research and Development, Stadtfelddamm 34, 30625 Hannover, Germany
| | - Ortwin Naujok
- Institute of Clinical Biochemistry, Hannover Medical School, Carl-Neuberg-Strasse 1, 30625 Hannover, Germany;
| | - Ekaterina Riabchenko
- Institute for Biomedical Systems, National Research University of Electronic Technology, 124498 Moscow, Russia; (E.R.); (A.D.)
| | - Yuliia Nemirovska
- Institute for Problems of Cryobiology and Cryomedicine of the National Academy of Sciences of Ukraine, 23 Pereyaslavsky Street, 61015 Kharkiv, Ukraine; (D.T.); (Y.N.); (A.Y.P.)
| | - Arseny Danilov
- Institute for Biomedical Systems, National Research University of Electronic Technology, 124498 Moscow, Russia; (E.R.); (A.D.)
| | - Alexander Y. Petrenko
- Institute for Problems of Cryobiology and Cryomedicine of the National Academy of Sciences of Ukraine, 23 Pereyaslavsky Street, 61015 Kharkiv, Ukraine; (D.T.); (Y.N.); (A.Y.P.)
| | - Birgit Glasmacher
- Institute for Multiphase Processes, Leibniz University Hannover, An der Universität 1, Building 8143, 30823 Garbsen, Germany; (V.M.); (D.K.); (B.G.)
- Lower Saxony Centre for Biomedical Engineering, Implant Research and Development, Stadtfelddamm 34, 30625 Hannover, Germany
| |
Collapse
|
3
|
Abstract
Motivated by the need to understand better the dynamics of non-Fourier fluids, we examine the linear and weakly nonlinear stabilities of a horizontal layer of fluid obeying the Maxwell–Cattaneo relationship of heat flux and temperature using three different forms of the time derivative of the heat flux. Linear stability mode regime diagrams in the parameter plane have been established and used to summarize the linear instabilities. The energy balance of the system is used to identify the mechanism by which the Maxwell–Cattaneo effect (i) introduces overstability, (ii) leads to preferred stationary modes with the critical Rayleigh and wavelengths either both increasing or both decreasing, (iii) gives rise to instabilities in a layer heated from above, and (iv) enhances heat transfer. A formal weakly nonlinear analysis leads to evolution equations for the amplitudes of linear instability modes. It is shown that the amplitude of the stationary mode obeys an equation of the Landau–Stuart type. The two equally excitable overstable modes obey two equations of the same type coupled by an interaction term. The evolution of the different amplitudes leads to supercritical stability, supercritical instability or subcritical instability depending on the model and parameter values. The results are presented in regime diagrams.
Collapse
|
4
|
Bissell JJ. Thermal convection in a magnetized conducting fluid with the Cattaneo-Christov heat-flow model. Proc Math Phys Eng Sci 2016; 472:20160649. [PMID: 27956886 PMCID: PMC5134317 DOI: 10.1098/rspa.2016.0649] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2016] [Accepted: 10/17/2016] [Indexed: 11/12/2022] Open
Abstract
By substituting the Cattaneo-Christov heat-flow model for the more usual parabolic Fourier law, we consider the impact of hyperbolic heat-flow effects on thermal convection in the classic problem of a magnetized conducting fluid layer heated from below. For stationary convection, the system is equivalent to that studied by Chandrasekhar (Hydrodynamic and Hydromagnetic Stability, 1961), and with free boundary conditions we recover the classical critical Rayleigh number [Formula: see text] which exhibits inhibition of convection by the field according to [Formula: see text] as [Formula: see text], where Q is the Chandrasekhar number. However, for oscillatory convection we find that the critical Rayleigh number [Formula: see text] is given by a more complicated function of the thermal Prandtl number [Formula: see text], magnetic Prandtl number [Formula: see text] and Cattaneo number C. To elucidate features of this dependence, we neglect [Formula: see text] (in which case overstability would be classically forbidden), and thereby obtain an expression for the Rayleigh number that is far less strongly inhibited by the field, with limiting behaviour [Formula: see text], as [Formula: see text]. One consequence of this weaker dependence is that onset of instability occurs as overstability provided C exceeds a threshold value CT(Q); indeed, crucially we show that when Q is large, [Formula: see text], meaning that oscillatory modes are preferred even when C itself is small. Similar behaviour is demonstrated in the case of fixed boundaries by means of a novel numerical solution.
Collapse
Affiliation(s)
- J. J. Bissell
- Department of Physics, University of Bath, BA2 7AY Bath, UK
| |
Collapse
|
5
|
Bissell JJ. On oscillatory convection with the Cattaneo-Christov hyperbolic heat-flow model. Proc Math Phys Eng Sci 2015; 471:20140845. [PMID: 25792960 PMCID: PMC4353040 DOI: 10.1098/rspa.2014.0845] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2014] [Accepted: 01/21/2015] [Indexed: 11/12/2022] Open
Abstract
Adoption of the hyperbolic Cattaneo-Christov heat-flow model in place of the more usual parabolic Fourier law is shown to raise the possibility of oscillatory convection in the classic Bénard problem of a Boussinesq fluid heated from below. By comparing the critical Rayleigh numbers for stationary and oscillatory convection, Rc and RS respectively, oscillatory convection is found to represent the preferred form of instability whenever the Cattaneo number C exceeds a threshold value CT≥8/27π2≈0.03. In the case of free boundaries, analytical approaches permit direct treatment of the role played by the Prandtl number [Formula: see text], which-in contrast to the classical stationary scenario-can impact on oscillatory modes significantly owing to the non-zero frequency of convection. Numerical investigation indicates that the behaviour found analytically for free boundaries applies in a qualitatively similar fashion for fixed boundaries, while the threshold Cattaneo number CT is computed as a function of [Formula: see text] for both boundary regimes.
Collapse
Affiliation(s)
- J. J. Bissell
- Department of Mathematical Sciences, University of Durham, Durham City DH1 3LE, UK
| |
Collapse
|