Bridging Three Orders of Magnitude: Multiple Scattered Waves Sense Fractal Microscopic Structures via Dispersion.

Wave scattering provides profound insight into the structure of matter. Typically, the ability to sense microstructure is determined by the ratio of scatterer size to probing wavelength. Here, we address the question of whether macroscopic waves can report back the presence and distribution of microscopic scatterers despite several orders of magnitude difference in scale between wavelength and scatterer size. In our analysis, monosized hard scatterers 5 µm in radius are immersed in lossless gelatin phantoms to investigate the effect of multiple reflections on the propagation of shear waves with millimeter wavelength. Steady-state monochromatic waves are imaged in situ via magnetic resonance imaging, enabling quantification of the phase velocity at a voxel size big enough to contain thousands of individual scatterers, but small enough to resolve the wavelength. We show in theory, experiments, and simulations that the resulting coherent superposition of multiple reflections gives rise to power-law dispersion at the macroscopic scale if the scatterer distribution exhibits apparent fractality over an effective length scale that is comparable to the probing wavelength. Since apparent fractality is naturally present in any random medium, microstructure can thereby leave its fingerprint on the macroscopically quantifiable power-law exponent. Our results are generic to wave phenomena and carry great potential for sensing microstructure that exhibits intrinsic fractality, such as, for instance, vasculature.

Multi-scale parameterisation of a myocardial perfusion model using whole-organ arterial networks.

A method to extract myocardial coronary permeabilities appropriate to parameterise a continuum porous perfusion model using the underlying anatomical vascular network is developed. Canine and porcine whole-heart discrete arterial models were extracted from high-resolution cryomicrotome vessel image stacks. Five parameterisation methods were considered that are primarily distinguished by the level of anatomical data used in the definition of the permeability and pressure-coupling fields. Continuum multi-compartment porous perfusion model pressure results derived using these parameterisation methods were compared quantitatively via a root-mean-square metric to the Poiseuille pressure solved on the discrete arterial vasculature. The use of anatomical detail to parameterise the porous medium significantly improved the continuum pressure results. The majority of this improvement was attributed to the use of anatomically-derived pressure-coupling fields. It was found that the best results were most reliably obtained by using porosity-scaled isotropic permeabilities and anatomically-derived pressure-coupling fields. This paper presents the first continuum perfusion model where all parameters were derived from the underlying anatomical vascular network.

Myocardial perfusion distribution and coronary arterial pressure and flow signals: clinical relevance in relation to multiscale modeling, a review.

Coronary artery disease, CAD, is associated with both narrowing of the epicardial coronary arteries and microvascular disease, thereby limiting coronary flow and myocardial perfusion. CAD accounts for almost 2 million deaths within the European Union on an annual basis. In this paper, we review the physiological and pathophysiological processes underlying clinical decision making in coronary disease as well as the models for interpretation of the underlying physiological mechanisms. Presently, clinical decision making is based on non-invasive magnetic resonance imaging, MRI, of myocardial perfusion and invasive coronary hemodynamic measurements of coronary pressure and Doppler flow velocity signals obtained during catheterization. Within the euHeart project, several innovations have been developed and applied to improve diagnosis-based understanding of the underlying biophysical processes. Specifically, MRI perfusion data interpretation has been advanced by the gradientogram, a novel graphical representation of the spatiotemporal myocardial perfusion gradient. For hemodynamic data, functional indices of coronary stenosis severity that do not depend on maximal vasodilation are proposed and the Valsalva maneuver for indicating the extravascular resistance component of the coronary circulation has been introduced. Complementary to these advances, model innovation has been directed to the porous elastic model coupled to a one-dimensional model of the epicardial arteries. The importance of model development is related to the integration of information from different modalities, which in isolation often result in conflicting treatment recommendations.

Parameterisation of multi-scale continuum perfusion models from discrete vascular networks.

Experimental data and advanced imaging techniques are increasingly enabling the extraction of detailed vascular anatomy from biological tissues. Incorporation of anatomical data within perfusion models is non-trivial, due to heterogeneous vessel density and disparate radii scales. Furthermore, previous idealised networks have assumed a spatially repeating motif or periodic canonical cell, thereby allowing for a flow solution via homogenisation. However, such periodicity is not observed throughout anatomical networks. In this study, we apply various spatial averaging methods to discrete vascular geometries in order to parameterise a continuum model of perfusion. Specifically, a multi-compartment Darcy model was used to provide vascular scale separation for the fluid flow. Permeability tensor fields were derived from both synthetic and anatomically realistic networks using (1) porosity-scaled isotropic, (2) Huyghe and Van Campen, and (3) projected-PCA methods. The Darcy pressure fields were compared via a root-mean-square error metric to an averaged Poiseuille pressure solution over the same domain. The method of Huyghe and Van Campen performed better than the other two methods in all simulations, even for relatively coarse networks. Furthermore, inter-compartment volumetric flux fields, determined using the spatially averaged discrete flux per unit pressure difference, were shown to be accurate across a range of pressure boundary conditions. This work justifies the application of continuum flow models to characterise perfusion resulting from flow in an underlying vascular network.

A finite-element approach to the direct computation of relative cardiovascular pressure from time-resolved MR velocity data.

The evaluation of cardiovascular velocities, their changes through the cardiac cycle and the consequent pressure gradients has the capacity to improve understanding of subject-specific blood flow in relation to adjacent soft tissue movements. Magnetic resonance time-resolved 3D phase contrast velocity acquisitions (4D flow) represent an emerging technology capable of measuring the cyclic changes of large scale, multi-directional, subject-specific blood flow. A subsequent evaluation of pressure differences in enclosed vascular compartments is a further step which is currently not directly available from such data. The focus of this work is to address this deficiency through the development of a novel simulation workflow for the direct computation of relative cardiovascular pressure fields. Input information is provided by enhanced 4D flow data and derived MR domain masking. The underlying methodology shows numerical advantages in terms of robustness, global domain composition, the isolation of local fluid compartments and a treatment of boundary conditions. This approach is demonstrated across a range of validation examples which are compared with analytic solutions. Four subject-specific test cases are subsequently run, showing good agreement with previously published calculations of intra-vascular pressure differences. The computational engine presented in this work contributes to non-invasive access to relative pressure fields, incorporates the effects of both blood flow acceleration and viscous dissipation, and enables enhanced evaluation of cardiovascular blood flow.