1
|
Obermeier L, Vellguth K, Schlief A, Tautz L, Bruening J, Knosalla C, Kuehne T, Solowjowa N, Goubergrits L. CT-Based Simulation of Left Ventricular Hemodynamics: A Pilot Study in Mitral Regurgitation and Left Ventricle Aneurysm Patients. Front Cardiovasc Med 2022; 9:828556. [PMID: 35391837 PMCID: PMC8980692 DOI: 10.3389/fcvm.2022.828556] [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/03/2021] [Accepted: 02/03/2022] [Indexed: 12/30/2022] Open
Abstract
Background Cardiac CT (CCT) is well suited for a detailed analysis of heart structures due to its high spatial resolution, but in contrast to MRI and echocardiography, CCT does not allow an assessment of intracardiac flow. Computational fluid dynamics (CFD) can complement this shortcoming. It enables the computation of hemodynamics at a high spatio-temporal resolution based on medical images. The aim of this proposed study is to establish a CCT-based CFD methodology for the analysis of left ventricle (LV) hemodynamics and to assess the usability of the computational framework for clinical practice. Materials and Methods The methodology is demonstrated by means of four cases selected from a cohort of 125 multiphase CCT examinations of heart failure patients. These cases represent subcohorts of patients with and without LV aneurysm and with severe and no mitral regurgitation (MR). All selected LVs are dilated and characterized by a reduced ejection fraction (EF). End-diastolic and end-systolic image data was used to reconstruct LV geometries with 2D valves as well as the ventricular movement. The intraventricular hemodynamics were computed with a prescribed-motion CFD approach and evaluated in terms of large-scale flow patterns, energetic behavior, and intraventricular washout. Results In the MR patients, a disrupted E-wave jet, a fragmentary diastolic vortex formation and an increased specific energy dissipation in systole are observed. In all cases, regions with an impaired washout are visible. The results furthermore indicate that considering several cycles might provide a more detailed view of the washout process. The pre-processing times and computational expenses are in reach of clinical feasibility. Conclusion The proposed CCT-based CFD method allows to compute patient-specific intraventricular hemodynamics and thus complements the informative value of CCT. The method can be applied to any CCT data of common quality and represents a fair balance between model accuracy and overall expenses. With further model enhancements, the computational framework has the potential to be embedded in clinical routine workflows, to support clinical decision making and treatment planning.
Collapse
Affiliation(s)
- Lukas Obermeier
- Institute of Computer-Assisted Cardiovascular Medicine, Charité - Universitätsmedizin Berlin, Berlin, Germany
| | - Katharina Vellguth
- Institute of Computer-Assisted Cardiovascular Medicine, Charité - Universitätsmedizin Berlin, Berlin, Germany
| | - Adriano Schlief
- Institute of Computer-Assisted Cardiovascular Medicine, Charité - Universitätsmedizin Berlin, Berlin, Germany
| | - Lennart Tautz
- Institute of Computer-Assisted Cardiovascular Medicine, Charité - Universitätsmedizin Berlin, Berlin, Germany
- Fraunhofer Institute for Digital Medicine MEVIS, Bremen, Germany
| | - Jan Bruening
- Institute of Computer-Assisted Cardiovascular Medicine, Charité - Universitätsmedizin Berlin, Berlin, Germany
| | - Christoph Knosalla
- Department of Cardiothoracic and Vascular Surgery, German Heart Center Berlin, Berlin, Germany
- DZHK (German Centre for Cardiovascular Research), Partner Site Berlin, Berlin, Germany
- Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Humboldt - Universität zu Berlin and Berlin Institute of Health, Berlin, Germany
| | - Titus Kuehne
- Institute of Computer-Assisted Cardiovascular Medicine, Charité - Universitätsmedizin Berlin, Berlin, Germany
- DZHK (German Centre for Cardiovascular Research), Partner Site Berlin, Berlin, Germany
- Department of Congenital Heart Disease, German Heart Center Berlin, Berlin, Germany
| | - Natalia Solowjowa
- Department of Cardiothoracic and Vascular Surgery, German Heart Center Berlin, Berlin, Germany
| | - Leonid Goubergrits
- Institute of Computer-Assisted Cardiovascular Medicine, Charité - Universitätsmedizin Berlin, Berlin, Germany
- Einstein Center Digital Future, Berlin, Germany
| |
Collapse
|
2
|
Abe H, Caracciolo G, Kheradvar A, Pedrizzetti G, Khandheria BK, Narula J, Sengupta PP. Contrast echocardiography for assessing left ventricular vortex strength in heart failure: a prospective cohort study. ACTA ACUST UNITED AC 2013; 14:1049-60. [DOI: 10.1093/ehjci/jet049] [Citation(s) in RCA: 89] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
|
3
|
Zhang H, Liu L, Chen L, Ma N, Zhou L, Liu Y, Li Z, Liu C, Hou R, Zhu S. The evolution of intraventricular vortex during ejection studied by using vector flow mapping. Echocardiography 2012; 30:27-36. [PMID: 22985123 DOI: 10.1111/j.1540-8175.2012.01806.x] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022] Open
Abstract
AIMS The purpose of this study was to assess the evolution of intraventricular vortex during left ventricular (LV) ejection. METHODS Vector flow mapping was performed in 51 patients with coronary artery disease and LV ejection fraction (EF) >50%, 70 patients with EF <50% (13 with coronary artery disease and 57 with dilated cardiomyopathy), and 62 healthy volunteers. RESULTS In normals and patients with EF >50%, the intraventricular vortex dissipated quickly during early ejection. In patients with EF <50%, the vortex stayed mainly at apex and persisted for a significantly longer time. The evolution of vortex during ejection was significantly correlated with QRS width, EF, fractional shortening, LV outflow velocity time integral, wall motion score index (WMSI), LV dimensions, left atrial diameter, and diastolic mitral annular velocities. LV end-diastolic short diameter and WMSI were the independent determinants of the duration of vortex (R(2) = 0.482, P < 0.001). End-systolic short diameter and apical WMSI were the independent determinants of duration of vortex corrected for ejection time (R(2) = 0.565, P < 0.001). End-systolic short diameter was the independent determinant of percentage change in vortex area during early ejection (R(2) = 0.355, P < 0.001). End-systolic short diameter and ejection time were the independent determinants of percentage change in vortex flow volume (R(2) = 0.415, P < 0.001). CONCLUSIONS In patients with LV systolic dysfunction, the vortex persists during ejection and stays mainly at apex. The vortex evolution during ejection is closely associated with LV dimensions and functions.
Collapse
Affiliation(s)
- Haibin Zhang
- Department of Ultrasound, PLA 210th Hospital, Dalian, China
| | | | | | | | | | | | | | | | | | | |
Collapse
|
4
|
Eriksson J, Dyverfeldt P, Engvall J, Bolger AF, Ebbers T, Carlhäll CJ. Quantification of presystolic blood flow organization and energetics in the human left ventricle. Am J Physiol Heart Circ Physiol 2011; 300:H2135-41. [PMID: 21421820 DOI: 10.1152/ajpheart.00993.2010] [Citation(s) in RCA: 97] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Intracardiac blood flow patterns are potentially important to cardiac pumping efficiency. However, these complex flow patterns remain incompletely characterized both in health and disease. We hypothesized that normal left ventricular (LV) blood flow patterns would preferentially optimize a portion of the end-diastolic volume (LVEDV) for effective and rapid systolic ejection by virtue of location near and motion towards the LV outflow tract (LVOT). Three-dimensional cine velocity and morphological data were acquired in 12 healthy persons and 1 patient with dilated cardiomyopathy using MRI. A previously validated method was used for analysis in which the LVEDV was separated into four functional flow components based on the blood's locations at the beginning and end of the cardiac cycle. Each component's volume, kinetic energy (KE), site, direction, and linear momentum relative to the LVOT were calculated. Of the four components, the LV inflow that passes directly to outflow in a single cardiac cycle (Direct Flow) had the largest volume. At the time of isovolumic contraction, Direct Flow had the greatest amount of KE and the most favorable combination of distance, angle, and linear momentum relative to the LVOT. Atrial contraction boosted the late diastolic KE of the ejected components. We conclude that normal diastolic LV flow creates favorable conditions for ensuing ejection, defined by proximity and energetics, for the Direct Flow, and that atrial contraction augments the end-diastolic KE of the ejection volume. The correlation of Direct Flow characteristics with ejection efficiency might be a relevant investigative target in cardiac failure.
Collapse
Affiliation(s)
- Jonatan Eriksson
- Dept. of Medical and Health Sciences, Division of Cardiovascular Medicine, Linköping Univ., SE-581 85 Linköping .
| | | | | | | | | | | |
Collapse
|
5
|
Eriksson J, Carlhäll CJ, Dyverfeldt P, Engvall J, Bolger AF, Ebbers T. Semi-automatic quantification of 4D left ventricular blood flow. J Cardiovasc Magn Reson 2010; 12:9. [PMID: 20152026 PMCID: PMC2831022 DOI: 10.1186/1532-429x-12-9] [Citation(s) in RCA: 138] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2009] [Accepted: 02/12/2010] [Indexed: 12/02/2022] Open
Abstract
BACKGROUND The beating heart is the generator of blood flow through the cardiovascular system. Within the heart's own chambers, normal complex blood flow patterns can be disturbed by diseases. Methods for the quantification of intra-cardiac blood flow, with its 4D (3D+time) nature, are lacking. We sought to develop and validate a novel semi-automatic analysis approach that integrates flow and morphological data. METHOD In six healthy subjects and three patients with dilated cardiomyopathy, three-directional, three-dimensional cine phase-contrast cardiovascular magnetic resonance (CMR) velocity data and balanced steady-state free-precession long- and short-axis images were acquired. The LV endocardium was segmented from the short-axis images at the times of isovolumetric contraction (IVC) and isovolumetric relaxation (IVR). At the time of IVC, pathlines were emitted from the IVC LV blood volume and traced forwards and backwards in time until IVR, thus including the entire cardiac cycle. The IVR volume was used to determine if and where the pathlines left the LV. This information was used to automatically separate the pathlines into four different components of flow: Direct Flow, Retained Inflow, Delayed Ejection Flow and Residual Volume. Blood volumes were calculated for every component by multiplying the number of pathlines with the blood volume represented by each pathline. The accuracy and inter- and intra-observer reproducibility of the approach were evaluated by analyzing volumes of LV inflow and outflow, the four flow components, and the end-diastolic volume. RESULTS The volume and distribution of the LV flow components were determined in all subjects. The calculated LV outflow volumes [ml] (67 +/- 13) appeared to fall in between those obtained by through-plane phase-contrast CMR (77 +/- 16) and Doppler ultrasound (58 +/- 10), respectively. Calculated volumes of LV inflow (68 +/- 11) and outflow (67 +/- 13) were well matched (NS). Low inter- and intra-observer variability for the assessment of the volumes of the flow components was obtained. CONCLUSIONS This semi-automatic analysis approach for the quantification of 4D blood flow resulted in accurate LV inflow and outflow volumes and a high reproducibility for the assessment of LV flow components.
Collapse
Affiliation(s)
- Jonatan Eriksson
- Division of Cardiovascular Medicine, Department of Medical and Health Sciences, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
| | - Carl Johan Carlhäll
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
- Department of Clinical Physiology, Linköping University Hospital, Linköping, Sweden
| | - Petter Dyverfeldt
- Division of Cardiovascular Medicine, Department of Medical and Health Sciences, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
- Division of Applied Thermodynamics and Fluid Mechanics, Department of Management and Engineering, Linköping University, Linköping, Sweden
| | - Jan Engvall
- Division of Cardiovascular Medicine, Department of Medical and Health Sciences, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
- Department of Clinical Physiology, Linköping University Hospital, Linköping, Sweden
| | - Ann F Bolger
- Division of Cardiovascular Medicine, Department of Medical and Health Sciences, Linköping University, Linköping, Sweden
- Department of Medicine, University of California, San Francisco, California, USA
| | - Tino Ebbers
- Division of Cardiovascular Medicine, Department of Medical and Health Sciences, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
- Division of Applied Thermodynamics and Fluid Mechanics, Department of Management and Engineering, Linköping University, Linköping, Sweden
| |
Collapse
|
6
|
Tanaka M, Sakamoto T, Sugawara S, Nakajima H, Katahira Y, Ohtsuki S, Kanai H. Blood flow structure and dynamics, and ejection mechanism in the left ventricle: Analysis using echo-dynamography. J Cardiol 2008; 52:86-101. [DOI: 10.1016/j.jjcc.2008.05.005] [Citation(s) in RCA: 52] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/28/2008] [Revised: 05/15/2008] [Accepted: 05/16/2008] [Indexed: 11/16/2022]
|
7
|
|
8
|
Fyrenius A, Wigström L, Ebbers T, Karlsson M, Engvall J, Bolger AF. Three dimensional flow in the human left atrium. Heart 2001; 86:448-55. [PMID: 11559688 PMCID: PMC1729920 DOI: 10.1136/heart.86.4.448] [Citation(s) in RCA: 130] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/03/2022] Open
Abstract
BACKGROUND Abnormal flow patterns in the left atrium in atrial fibrillation or mitral stenosis are associated with an increased risk of thrombosis and systemic embolisation; the characteristics of normal atrial flow that avoid stasis have not been well defined. OBJECTIVES To present a three dimensional particle trace visualisation of normal left atrial flow in vivo, constructed from flow velocities in three dimensional space. METHODS Particle trace visualisation of time resolved three dimensional magnetic resonance imaging velocity measurements was used to provide a display of intracardiac flow without the limitations of angle sensitivity or restriction to imaging planes. Global flow patterns of the left atrium were studied in 11 healthy volunteers. RESULTS In all subjects vortical flow was observed in the atrium during systole and diastolic diastasis (mean (SD) duration of systolic vortex, 280 (77) ms; and of diastolic vortex, 256 (118) ms). The volume incorporated and recirculated within the vortices originated predominantly from the left pulmonary veins. Inflow from the right veins passed along the vortex periphery, constrained between the vortex and the atrial wall. CONCLUSIONS Global left atrial flow in the normal human heart comprises consistent patterns specific to the phase of the cardiac cycle. Separate paths of left and right pulmonary venous inflow and vortex formation may have beneficial effects in avoiding left atrial stasis in the normal subject in sinus rhythm.
Collapse
Affiliation(s)
- A Fyrenius
- Department of Medicine and Care, Clinical Physiology, Linköping University, Linköping, Sweden
| | | | | | | | | | | |
Collapse
|
9
|
Firstenberg MS, Vandervoort PM, Greenberg NL, Smedira NG, McCarthy PM, Garcia MJ, Thomas JD. Noninvasive estimation of transmitral pressure drop across the normal mitral valve in humans: importance of convective and inertial forces during left ventricular filling. J Am Coll Cardiol 2000; 36:1942-9. [PMID: 11092668 DOI: 10.1016/s0735-1097(00)00963-3] [Citation(s) in RCA: 50] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Abstract
OBJECTIVES We hypothesized that color M-mode (CMM) images could be used to solve the Euler equation, yielding regional pressure gradients along the scanline, which could then be integrated to yield the unsteady Bernoulli equation and estimate noninvasively both the convective and inertial components of the transmitral pressure difference. BACKGROUND Pulsed and continuous wave Doppler velocity measurements are routinely used clinically to assess severity of stenotic and regurgitant valves. However, only the convective component of the pressure gradient is measured, thereby neglecting the contribution of inertial forces, which may be significant, particularly for nonstenotic valves. Color M-mode provides a spatiotemporal representation of flow across the mitral valve. METHODS In eight patients undergoing coronary artery bypass grafting, high-fidelity left atrial and ventricular pressure measurements were obtained synchronously with transmitral CMM digital recordings. The instantaneous diastolic transmitral pressure difference was computed from the M-mode spatiotemporal velocity distribution using the unsteady flow form of the Bernoulli equation and was compared to the catheter measurements. RESULTS From 56 beats in 16 hemodynamic stages, inclusion of the inertial term ([deltapI]max = 1.78+/-1.30 mm Hg) in the noninvasive pressure difference calculation significantly increased the temporal correlation with catheter-based measurement (r = 0.35+/-0.24 vs. 0.81+/-0.15, p< 0.0001). It also allowed an accurate approximation of the peak pressure difference ([deltapc+I]max = 0.95 [delta(p)cathh]max + 0.24, r = 0.96, p<0.001, error = 0.08+/-0.54 mm Hg). CONCLUSIONS Inertial forces are significant components of the maximal pressure drop across the normal mitral valve. These can be accurately estimated noninvasively using CMM recordings of transmitral flow, which should improve the understanding of diastolic filling and function of the heart.
Collapse
Affiliation(s)
- M S Firstenberg
- Cardiovascular Imaging Center, Department of Cardiology, Cleveland Clinic Foundation, Ohio 44195, USA
| | | | | | | | | | | | | |
Collapse
|
10
|
Lemmon JD, Yoganathan AP. Three-dimensional computational model of left heart diastolic function with fluid-structure interaction. J Biomech Eng 2000; 122:109-17. [PMID: 10834150 DOI: 10.1115/1.429648] [Citation(s) in RCA: 57] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Aided by advancements in computer speed and modeling techniques, computational modeling of cardiac function has continued to develop over the past twenty years. The goal of the current study was to develop a computational model that provides blood-tissue interaction under physiologic flow conditions, and apply it to a thin-walled model of the left heart. To accomplish this goal, the Immersed Boundary Method was used to study the interaction of the tissue and blood in response to fluid forces and changes in tissue pathophysiology. The fluid mass and momentum conservation equations were solved using Patankar's Semi-Implicit Method for Pressure Linked Equations (SIMPLE). A left heart model was developed to examine diastolic function, and consisted of the left ventricle, left atrium, and pulmonary flow. The input functions for the model included the pulmonary driving pressure and time-dependent relationship for changes in chamber tissue properties during the simulation. The results obtained from the left heart model were compared to clinically observed diastolic flow conditions for validation. The inflow velocities through the mitral valve corresponded with clinical values (E-wave = 74.4 cm/s, A-wave = 43 cm/s, and E/A = 1.73). The pressure traces for the atrium and ventricle, and the appearance of the ventricular flow fields throughout filling, agreed with those observed in the heart. In addition, the atrial flow fields could be observed in this model and showed the conduit and pump functions that current theory suggests. The ability to examine atrial function in the present model is something not described previously in computational simulations of cardiac function.
Collapse
Affiliation(s)
- J D Lemmon
- Schools of Mechanical and Biomedical Engineering, Georgia Institute of Technology, Atlanta 30332, USA.
| | | |
Collapse
|
11
|
Houlind K, Schroeder AP, Egeblad H, Pedersen EM. Age-dependent changes in spatial and temporal blood velocity distribution of early left ventricular filling. Magn Reson Imaging 1999; 17:859-68. [PMID: 10402593 DOI: 10.1016/s0730-725x(99)00018-1] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
This study describes early diastolic inflow dynamics based on three-directional magnetic resonance velocity data and investigates age-dependent changes in early diastolic inflow characteristics. We examined 26 young healthy volunteers age 25 (3) years (mean, SD), and 23 healthy older volunteers age 63 (8) years. Three-directional magnetic resonance velocity mapping was performed in a long axis plane through the heart. Transverse velocity profiles were read in five different positions in the early diastolic inflow stream of the left ventricle. The size and timing of the maximum velocities at each level were recorded and the repeatability of the method was tested. Compared with the younger group, the older group was characterized by: 1) lower maximum velocity in all positions, 2) increased deceleration of blood downstream from the mitral leaflet tips, and 3) delayed velocity propagation. The described method was repeatable and enabled detection of the age-dependent differences between groups of normal subjects. In conclusion, the early diastolic inflow pattern changes with age, probably reflecting changes in the diastolic function of the myocardium.
Collapse
Affiliation(s)
- K Houlind
- Department of Cardiothoracic and Vascular Surgery, Skejby Sygehus, Aarhus University Hospital, Denmark.
| | | | | | | |
Collapse
|
12
|
Wigström L, Ebbers T, Fyrenius A, Karlsson M, Engvall J, Wranne B, Bolger AF. Particle trace visualization of intracardiac flow using time-resolved 3D phase contrast MRI. Magn Reson Med 1999; 41:793-9. [PMID: 10332856 DOI: 10.1002/(sici)1522-2594(199904)41:4<793::aid-mrm19>3.0.co;2-2] [Citation(s) in RCA: 146] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The flow patterns in the human heart are complex and difficult to visualize using conventional two-dimensional (2D) modalities, whether they depict a single velocity component (Doppler echocardiography) or all three components in a few slices (2D phase contrast MRI). To avoid these shortcomings, a temporally resolved 3D phase contrast technique was used to derive data describing the intracardiac velocity fields in normal volunteers. The MRI data were corrected for phase shifts caused by eddy currents and concomitant gradient fields, with improvement in the accuracy of subsequent flow visualizations. Pathlines describing the blood pathways through the heart were generated from the temporally resolved velocity data, starting from user-specified locations and time frames. Flow trajectories were displayed as 3D particle traces, with simultaneous demonstration of morphologic 2D slices. This type of visualization is intuitive and interactive and may extend our understanding of dynamic and previously unrecognized patterns of intracardiac flow.
Collapse
Affiliation(s)
- L Wigström
- Department of Clinical Physiology, Linköping University, Sweden.
| | | | | | | | | | | | | |
Collapse
|
13
|
Wigstr�m L, Ebbers T, Fyrenius A, Karlsson M, Engvall J, Wranne B, Bolger AF. Particle trace visualization of intracardiac flow using time-resolved 3D phase contrast MRI. Magn Reson Med 1999. [DOI: 10.1002/(sici)1522-2594(199904)41:4%3c793::aid-mrm19%3e3.0.co;2-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/08/2023]
|