A coupled discrete/continuum model for describing cancer-therapeutic transport in the lung

A phase-field model for non-small cell lung cancer under the effects of immunotherapy

Formulating mathematical models that estimate tumor growth under therapy is vital for improving patient-specific treatment plans. In this context, we present our recent work on simulating non-small-scale cell lung cancer (NSCLC) in a simple, deterministic setting for two different patients receiving an immunotherapeutic treatment. At its core, our model consists of a Cahn-Hilliard-based phase-field model describing the evolution of proliferative and necrotic tumor cells. These are coupled to a simplified nutrient model that drives the growth of the proliferative cells and their decay into necrotic cells. The applied immunotherapy decreases the proliferative cell concentration. Here, we model the immunotherapeutic agent concentration in the entire lung over time by an ordinary differential equation (ODE). Finally, reaction terms provide a coupling between all these equations. By assuming spherical, symmetric tumor growth and constant nutrient inflow, we simplify this full 3D cancer simulation model to a reduced 1D model. We can then resort to patient data gathered from computed tomography (CT) scans over several years to calibrate our model. Our model covers the case in which the immunotherapy is successful and limits the tumor size, as well as the case predicting a sudden relapse, leading to exponential tumor growth. Finally, we move from the reduced model back to the full 3D cancer simulation in the lung tissue. Thereby, we demonstrate the predictive benefits that a more detailed patient-specific simulation including spatial information as a possible generalization within our framework could yield in the future.

An in silico assessment on the potential of using saline infusion to overcome non-confluent coagulation zone during two-probe, no-touch bipolar radiofrequency ablation of liver cancer

No-touch bipolar radiofrequency ablation (bRFA) is known to produce incomplete tumour ablation with a 'butterfly-shaped' coagulation zone when the interelectrode distance exceeds a certain threshold. Although non-confluent coagulation zone can be avoided by not implementing the no-touch mode, doing so exposes the patient to the risk of tumour track seeding. The present study investigates if prior infusion of saline into the tissue can overcome the issues of non-confluent or butterfly-shaped coagulation. A computational modelling approach based on the finite element method was carried out. A two-compartment model comprising the tumour that is surrounded by healthy liver tissue was developed. Three cases were considered; i) saline infusion into the tumour centre; ii) one-sided saline infusion outside the tumour; and iii) two-sided saline infusion outside the tumour. For each case, three different saline volumes were considered, i.e. 6, 14 and 22 ml. Saline concentration was set to 15% w/v. Numerical results showed that saline infusion into the tumour centre can overcome the butterfly-shaped coagulation only if the infusion volume is sufficient. On the other hand, one-sided infusion outside the tumour did not overcome this. Two-sided infusion outside the tumour produced confluent coagulation zone with the largest volume. Results obtained from the present study suggest that saline infusion, when carried out correctly, can be used to effectively eradicate liver cancer. This presents a practical solution to address non-confluent coagulation zone typical of that during two-probe bRFA treatment.

Multi-compartmental model of glymphatic clearance of solutes in brain tissue

The glymphatic system is the subject of numerous pieces of research in biology. Mathematical modelling plays a considerable role in this field since it can indicate the possible physical effects of this system and validate the biologists' hypotheses. The available mathematical models that describe the system at the scale of the brain (i.e. the macroscopic scale) are often solely based on the diffusion equation and do not consider the fine structures formed by the perivascular spaces. We therefore propose a mathematical model representing the time and space evolution of a mixture flowing through multiple compartments of the brain. We adopt a macroscopic point of view in which the compartments are all present at any point in space. The equations system is composed of two coupled equations for each compartment: One equation for the pressure of a fluid and one for the mass concentration of a solute. The fluid and solute can move from one compartment to another according to certain membrane conditions modelled by transfer functions. We propose to apply this new modelling framework to the clearance of 14C-inulin from the rat brain.

A 1D-0D-3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D-0D-3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D-0D-3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

A computational framework to simulate the thermochemical process during thermochemical ablation of biological tissues

Thermochemical ablation (TCA) is a thermal ablation therapy that utilises heat released from acid-base neutralisation reaction to destroy tumours. This procedure is a promising low-cost solution to existing thermal ablation treatments such as radiofrequency ablation (RFA) and microwave ablation (MWA). Studies have demonstrated that TCA can produce thermal damage that is on par with RFA and MWA when employed properly. Nevertheless, TCA remains a concept that is tested only in a few animal trials due to the risks involved as the result of uncontrolled infusion and incomplete acid-base reaction. In this study, a computational framework that simulates the thermochemical process of TCA is developed. The proposed framework consists of three physics, namely chemical flow, neutralisation reaction and heat transfer. An important parameter in the TCA framework is the neutralisation reaction rate constant, which has values in the order of 10⁸ m³/(mol⋅s). The present study will demonstrate that since the rate constant impacts only the rate and direction of the reaction but has little influence on the extent of reaction, it is possible to replicate the thermochemical process of TCA by employing significantly smaller values of rate constant that are numerically tractable. Comparisons of the numerical results against experimental studies from the literature supports this. The aim of this framework is for researchers to advance and develop TCA to gain an in-depth understanding of the fundamental mechanisms of TCA and to develop a safe treatment protocol of TCA in the hope of advancing TCA into clinical trials.

Structure engineering of alveoli-like ZSM-5 with encapsulated Pt nanoparticles for the enhanced benzene oxidation

Inspired by the alveolar configuration, an alveoli-like ZSM-5 and the corresponding platinum encapsulated nanocomposite (Pt@PZ5) were fabricated via a dual-template method and a controlled selective desilication-recrystallization strategy. The dimensions of the central cavity, interconnected zeolitic vesicles, and mesoporous shell could be tuned by adjusting the synthesis parameters, as verified by scanning electron microscopy, transmission electron microscopy, nitrogen physisorption investigations, X-ray photoelectron spectroscopy, and X-ray diffraction techniques. Thanks to these properties and merits, the alveoli-like Pt@PZ5 showed the highest catalytic performance with excellent stability, obtaining 100% benzene conversion at 180 °C. Adsorption experiments combined with a finite-element simulation study uncovered that the alveolar architecture could expedite the accumulation of reactants and boost mass transfer; the conversion of intermediates in the voids could be further facilitated, giving optimal catalytic performance. Additionally, the alveolar architecture is resistant to metal sintering (5-20 nm) and leaching, even after calcination at 850 °C for 360 min. This work provides an alveolar concept into the rational design of efficient catalysts for fundamental catalytic action.

How does saline backflow affect the treatment of saline-infused radiofrequency ablation?

BACKGROUND AND OBJECTIVE

Saline infusion is applied together with radiofrequency ablation (RFA) to enlarge the ablation zone. However, one of the issues with saline-infused RFA is backflow, which spreads saline along the insertion track. This raises the concern of not only thermally ablating the tissue within the backflow region, but also the loss of saline from the targeted tissue, which may affect the treatment efficacy.

METHODS

In the present study, 2D axisymmetric models were developed to investigate how saline backflow influence saline-infused RFA and whether the aforementioned concerns are warranted. Saline-infused RFA was described using the dual porosity-Joule heating model. The hydrodynamics of backflow was described using Poiseuille law by assuming the flow to be similar to that in a thin annulus. Backflow lengths of 3, 4.5, 6 and 9 cm were considered.

RESULTS

Results showed that there is no concern of thermally ablating the tissue in the backflow region. This is due to the Joule heating being inversely proportional to distance from the electrode to the fourth power. Results also indicated that larger backflow lengths led to larger growth of thermal damage along the backflow region and greater decrease in coagulation volume. Hence, backflow needs to be controlled to ensure an effective treatment of saline-infused RFA.

CONCLUSIONS

There is no risk of ablating tissues around the needle insertion track due to backflow. Instead, the risk of underablation as a result of the loss of saline due to backflow was found to be of greater concern.

Voxelized simulation of cerebral oxygen perfusion elucidates hypoxia in aged mouse cortex

Departures of normal blood flow and metabolite distribution from the cerebral microvasculature into neuronal tissue have been implicated with age-related neurodegeneration. Mathematical models informed by spatially and temporally distributed neuroimage data are becoming instrumental for reconstructing a coherent picture of normal and pathological oxygen delivery throughout the brain. Unfortunately, current mathematical models of cerebral blood flow and oxygen exchange become excessively large in size. They further suffer from boundary effects due to incomplete or physiologically inaccurate computational domains, numerical instabilities due to enormous length scale differences, and convergence problems associated with condition number deterioration at fine mesh resolutions. Our proposed simple finite volume discretization scheme for blood and oxygen microperfusion simulations does not require expensive mesh generation leading to the critical benefit that it drastically reduces matrix size and bandwidth of the coupled oxygen transfer problem. The compact problem formulation yields rapid and stable convergence. Moreover, boundary effects can effectively be suppressed by generating very large replica of the cortical microcirculation in silico using an image-based cerebrovascular network synthesis algorithm, so that boundaries of the perfusion simulations are far removed from the regions of interest. Massive simulations over sizeable portions of the cortex with feature resolution down to the micron scale become tractable with even modest computer resources. The feasibility and accuracy of the novel method is demonstrated and validated with in vivo oxygen perfusion data in cohorts of young and aged mice. Our oxygen exchange simulations quantify steep gradients near penetrating blood vessels and point towards pathological changes that might cause neurodegeneration in aged brains. This research aims to explain mechanistic interactions between anatomical structures and how they might change in diseases or with age. Rigorous quantification of age-related changes is of significant interest because it might aide in the search for imaging biomarkers for dementia and Alzheimer’s disease.

Brain function critically depends on the maintenance of physiological blood supply and metabolism in the cortex. Disturbances to adequate perfusion have been linked to age-related neurodegeneration. However, the precise correlation between age-related hemodynamic changes and the resulting decline in oxygen delivery is not well understood and has not been quantified. Therefore, we introduce a new compact, and therefore highly scalable, computational method for predicting the physiological relationship between hemodynamics and cortical oxygen perfusion for large sections of the cortical microcirculation. We demonstrate the novel mesh generation-free (MGF), multi-scale simulation approach through realistic in vivo case studies of cortical microperfusion in the mouse brain. We further validate mechanistic correlations and a quantitative relationship between blood flow and brain oxygenation using experimental data from cohorts of young, middle aged and old mouse brains. Our computational approach overcomes size and performance limitations of previous unstructured meshing techniques to enable the prediction of oxygen tension with a spatial resolution of least two orders of magnitude higher than previously possible. Our simulation results support the hypothesis that structural changes in the microvasculature induce hypoxic pockets in the aged brain that are absent in the healthy, young mouse.

Role of saline concentration during saline-infused radiofrequency ablation: Observation of secondary Joule heating along the saline-tissue interface

Infusion of saline prior to radiofrequency ablation (RFA) is known to enlarge the thermal coagulation zone. The abundance of ions in saline elevate the electrical conductivity of the saline-saturated region. This promotes greater electric current flow inside the tissue, which increases the amount of RF energy deposition and subsequently enlarges the coagulation zone. In theory, infusion of higher concentration of saline should lead to larger coagulation zone due to the greater number of ions. Nevertheless, existing studies on the effects of concentration on saline-infused RFA have been conflicting, with the exact role of saline concentration yet to be fully elucidated. In this paper, computational models of saline-infused RFA were developed to investigate the role of saline concentration on the outcome of saline-infused RFA. The elevation in tissue electrical conductivity was modelled using the microscopic mixture model, while RFA was modelled using the coupled dual porosity-Joule heating model. Results obtained indicated that the presence of a concentration threshold to which no further elevation in tissue electrical conductivity and enlargement in thermal coagulation can occur. This threshold was determined to be at 15% NaCl. Analysis of the Joule heating distribution revealed the presence of a secondary Joule heating site located along the interface between wet and dry tissue. This secondary Joule heating was responsible for the enlargement in coagulation volume and its rapid growth phase during ablation.

A 3D-1D coupled blood flow and oxygen transport model to generate microvascular networks

In this work, we introduce an algorithmic approach to generate microvascular networks starting from larger vessels that can be reconstructed without noticeable segmentation errors. Contrary to larger vessels, the reconstruction of fine-scale components of microvascular networks shows significant segmentation errors, and an accurate mapping is time and cost intense. Thus there is a need for fast and reliable reconstruction algorithms yielding surrogate networks having similar stochastic properties as the original ones. The microvascular networks are constructed in a marching way by adding vessels to the outlets of the vascular tree from the previous step. To optimise the structure of the vascular trees, we use Murray's law to determine the radii of the vessels and bifurcation angles. In each step, we compute the local gradient of the partial pressure of oxygen and adapt the orientation of the new vessels to this gradient. At the same time, we use the partial pressure of oxygen to check whether the considered tissue block is supplied sufficiently with oxygen. Computing the partial pressure of oxygen, we use a 3D-1D coupled model for blood flow and oxygen transport. To decrease the complexity of a fully coupled 3D model, we reduce the blood vessel network to a 1D graph structure and use a bi-directional coupling with the tissue which is described by a 3D homogeneous porous medium. The resulting surrogate networks are analysed with respect to morphological and physiological aspects.

From tumour perfusion to drug delivery and clinical translation of in silico cancer models

In silico cancer models have demonstrated great potential as a tool to improve drug design, optimise the delivery of drugs to target sites in the host tissue and, hence, improve therapeutic efficacy and patient outcome. However, there are significant barriers to the successful translation of in silico technology from bench to bedside. More precisely, the specification of unknown model parameters, the necessity for models to adequately reflect in vivo conditions, and the limited amount of pertinent validation data to evaluate models' accuracy and assess their reliability, pose major obstacles in the path towards their clinical translation. This review aims to capture the state-of-the-art in in silico cancer modelling of vascularised solid tumour growth, and identify the important advances and barriers to success of these models in clinical oncology. Particular emphasis has been put on continuum-based models of cancer since they - amongst the class of mechanistic spatio-temporal modelling approaches - are well-established in simulating transport phenomena and the biomechanics of tissues, and have demonstrated potential for clinical translation. Three important avenues in in silico modelling are considered in this contribution: first, since systemic therapy is a major cancer treatment approach, we start with an overview of the tumour perfusion and angiogenesis in silico models. Next, we present the state-of-the-art in silico work encompassing the delivery of chemotherapeutic agents to cancer nanomedicines through the bloodstream, and then review continuum-based modelling approaches that demonstrate great promise for successful clinical translation. We conclude with a discussion of what we view to be the key challenges and opportunities for in silico modelling in personalised and precision medicine.

Liver Bioreactor Design Issues of Fluid Flow and Zonation, Fibrosis, and Mechanics: A Computational Perspective

Tissue engineering, with the goal of repairing or replacing damaged tissue and organs, has continued to make dramatic science-based advances since its origins in the late 1980’s and early 1990’s. Such advances are always multi-disciplinary in nature, from basic biology and chemistry through physics and mathematics to various engineering and computer fields. This review will focus its attention on two topics critical for tissue engineering liver development: (a) fluid flow, zonation, and drug screening, and (b) biomechanics, tissue stiffness, and fibrosis, all within the context of 3D structures. First, a general overview of various bioreactor designs developed to investigate fluid transport and tissue biomechanics is given. This includes a mention of computational fluid dynamic methods used to optimize and validate these designs. Thereafter, the perspective provided by computer simulations of flow, reactive transport, and biomechanics responses at the scale of the liver lobule and liver tissue is outlined, in addition to how bioreactor-measured properties can be utilized in these models. Here, the fundamental issues of tortuosity and upscaling are highlighted, as well as the role of disease and fibrosis in these issues. Some idealized simulations of the effects of fibrosis on lobule drug transport and mechanics responses are provided to further illustrate these concepts. This review concludes with an outline of some practical applications of tissue engineering advances and how efficient computational upscaling techniques, such as dual continuum modeling, might be used to quantify the transition of bioreactor results to the full liver scale.

Shape-shifting thermal coagulation zone during saline-infused radiofrequency ablation: A computational study on the effects of different infusion location

BACKGROUND AND OBJECTIVE

The majority of the studies on radiofrequency ablation (RFA) have focused on enlarging the size of the coagulation zone. An aspect that is crucial but often overlooked is the shape of the coagulation zone. The shape is crucial because the majority of tumours are irregularly-shaped. In this paper, the ability to manipulate the shape of the coagulation zone following saline-infused RFA by altering the location of saline infusion is explored.

METHODS

A 3D model of the liver tissue was developed. Saline infusion was described using the dual porosity model, while RFA was described using the electrostatic and bioheat transfer equations. Three infusion locations were investigated, namely at the proximal end, the middle and the distal end of the electrode. Investigations were carried out numerically using the finite element method.

RESULTS

Results indicated that greater thermal coagulation was found in the region of tissue occupied by the saline bolus. Infusion at the middle of the electrode led to the largest coagulation volume followed by infusion at the proximal and distal ends. It was also found that the ability to delay roll-off, as commonly associated with saline-infused RFA, was true only for the case when infusion is carried out at the middle. When infused at the proximal and distal ends, the occurrence of roll-off was advanced. This may be due to the rapid and more intense heating experienced by the tissue when infusion is carried out at the electrode ends where Joule heating is dominant.

CONCLUSION

Altering the location of saline infusion can influence the shape of the coagulation zone following saline-infused RFA. The ability to 'shift' the coagulation zone to a desired location opens up great opportunities for the development of more precise saline-infused RFA treatment that targets specific regions within the tissue.

Correlative 3D Imaging and Microfluidic Modelling of Human Pulmonary Lymphatics using Immunohistochemistry and High-resolution μCT

Lung lymphatics maintain fluid homoeostasis by providing a drainage system that returns fluid, cells and metabolites to the circulatory system. The 3D structure of the human pulmonary lymphatic network is essential to lung function, but it is poorly characterised. Image-based 3D mathematical modelling of pulmonary lymphatic microfluidics has been limited by the lack of accurate and representative image geometries. This is due to the microstructural similarity of the lymphatics to the blood vessel network, the lack of lymphatic-specific biomarkers, the technical limitations associated with image resolution in 3D, and sectioning artefacts present in 2D techniques. We present a method that combines lymphatic specific (D240 antibody) immunohistochemistry (IHC), optimised high-resolution X-ray microfocus computed tomography (μCT) and finite-element mathematical modelling to assess the function of human peripheral lung tissue. The initial results identify lymphatic heterogeneity within and between lung tissue. Lymphatic vessel volume fraction and fractal dimension significantly decreases away from the lung pleural surface (p < 0.001, n = 25 and p < 0.01, n = 20, respectively). Microfluidic modelling successfully shows that in lung tissue the fluid derived from the blood vessels drains through the interstitium into the lymphatic vessel network and this drainage is different in the subpleural space compared to the intralobular space. When comparing lung tissue from health and disease, human pulmonary lymphatics were significantly different across five morphometric measures used in this study (p ≤ 0.0001). This proof of principle study establishes a new engineering technology and workflow for further studies of pulmonary lymphatics and demonstrates for the first time the combination of correlative μCT and IHC to enable 3D mathematical modelling of human lung microfluidics at micrometre resolution.

Angioarchitecture and hemodynamics of microvascular arterio-venous malformations

Introduction

Arteriovenous malformations (AVMs) are characterized by pathological high flow, low resistance connections between arteries and veins. Treatment is critically dependent on correct interpretation of angioarchitectural features. However, some microfistular AVMs do not match the characteristics described in current AVM classification systems. Therefore, we propose a new subgroup of microfistular AVMs, composed of enlarged, fistulous paths on the venous half of capillaries and/or dilated draining venules (hyperdynamic, capillary-venulous malformation [CV-AVM]). CV-AVMs still ensure arterial flow to the periphery and fistulous venous drainage is less pronounced than in classical AVMs such that these lesions are often misinterpreted as venous malformations.

Materials and methods

We developed a computational model to study the effects of microvascular anomalies on local hemodynamics, as well as their impact on angiographic contrast propagation. Flow rates and pressures were computed with a lumped parameter description, while contrast propagation was determined by solving the 1D advection-diffusion equation.

Results and conclusions

For the newly proposed CV-AVM angioarchitecture, the computational model predicts increased arterio-venous contrast agent transit times and highly dispersive transport characteristics, compared to microfistular, interstitial type IV AVMs and high flow type II and III AVMs. We related these findings to time-contrast intensity curves sampled from clinical angiographies and found that there is strong evidence for the existence of CV-AVM.

A multiscale model for heterogeneous tumor spheroid in vitro

In this paper, a novel multiscale method is proposed for the study of heterogeneous tumor spheroid growth in vitro. The entire tumor spheroid is described by an ellipsoid-based model while nutrient and other environmental factors are treated as continua. The ellipsoid-based discrete component is capable of incorporating mechanical effects and deformability, while keeping a minimum set of free variables to describe complex shape variations. Moreover, our purely cell-based description of tumor avoids the complex mutual conversion between a cell-based model and continuum model within a tumor, such as force and mass transformation. This advantage makes it highly suitable for the study of tumor spheroids in vitro whose size are normally less than 800 μm in diameter. In addition, our numerical scheme provides two computational options depending on tumor size. For a small or medium tumor spheroid, a three-dimensional (3D) numerical model can be directly applied. For a large spheroid, we suggest the use of a 3D-adapted 2D cross section configuration, which has not yet been explored in the literature, as an alternative for the theoretical investigation to bridge the gap between the 2D and 3D models. Our model and its implementations have been validated and applied to various studies given in the paper. The simulation results fit corresponding in vitro experimental observations very well.

Multiscale modelling of blood flow in cerebral microcirculation: Details at capillary scale control accuracy at the level of the cortex

Aging or cerebral diseases may induce architectural modifications in human brain microvascular networks, such as capillary rarefaction. Such modifications limit blood and oxygen supply to the cortex, possibly resulting in energy failure and neuronal death. Modelling is key in understanding how these architectural modifications affect blood flow and mass transfers in such complex networks. However, the huge number of vessels in the human brain—tens of billions—prevents any modelling approach with an explicit architectural representation down to the scale of the capillaries. Here, we introduce a hybrid approach to model blood flow at larger scale in the brain microcirculation, based on its multiscale architecture. The capillary bed, which is a space-filling network, is treated as a porous medium and modelled using a homogenized continuum approach. The larger arteriolar and venular trees, which cannot be homogenized because of their fractal-like nature, are treated as a network of interconnected tubes with a detailed representation of their spatial organization. The main contribution of this work is to devise a proper coupling model at the interface between these two components. This model is based on analytical approximations of the pressure field that capture the strong pressure gradients building up in the capillaries connected to arterioles or venules. We evaluate the accuracy of this model for both very simple architectures with one arteriole and/or one venule and for more complex ones, with anatomically realistic tree-like vessels displaying a large number of coupling sites. We show that the hybrid model is very accurate in describing blood flow at large scales and further yields a significant computational gain by comparison with a classical network approach. It is therefore an important step towards large scale simulations of cerebral blood flow and lays the groundwork for introducing additional levels of complexity in the future.

Hemodynamic Characterization of Peripheral Arterio-venous Malformations

Peripheral arterio-venous malformations (pAVMs) are congenital vascular anomalies that require treatment, due to their severe clinical consequences. The complexity of lesions often leads to misdiagnosis and ill-planned treatments. To improve disease management, we developed a computational model to quantify the hemodynamic effects of key angioarchitectural features of pAVMs. Hemodynamic results were used to predict the transport of contrast agent (CA), which allowed us to compare our findings to digital subtraction angiography (DSA) recordings of patients. The model is based on typical pAVM morphologies and a generic vessel network that represents realistic vascular feeding and draining components related to lesions. A lumped-parameter description of the vessel network was employed to compute blood pressure and flow rates. CA-transport was determined by coupling the model to a 1D advection-diffusion equation. Results show that the extent of hemodynamic effects of pAVMs, such as arterial steal and venous hypertension, strongly depends on the lesion type and its vascular architecture. Dimensions of shunting vessels strongly influence hemodynamic parameters. Our results underline the importance of the dynamics of CA-transport in diagnostic DSA images. In this context, we identified a set of temporal CA-transport parameters, which are indicative of the presence and specific morphology of pAVMs.

Hollow Alveolus-Like Nanovesicle Assembly with Metal-Encapsulated Hollow Zeolite Nanocrystals

Inspired by the vesicular structure of alveolus which has a porous nanovesicle structure facilitating the transport of oxygen and carbon dioxide, we designed a hollow nanovesicle assembly with metal-encapsulated hollow zeolite that would enhance diffusion of reactants/products and inhibit sintering and leaching of active metals. This zeolitic nanovesicle has been successfully synthesized by a strategy which involves a one-pot hydrothermal synthesis of hollow assembly of metal-containing solid zeolite crystals without a structural template and a selective desilication-recrystallization accompanied by leaching-hydrolysis to convert the metal-containing solid crystals into metal-encapsulated hollow crystals. We demonstrate the strategy in synthesizing a hollow nanovesicle assembly of Fe2O3-encapsulated hollow crystals of ZSM-5 zeolite. This material possesses a microporous (0.4-0.6 nm) wall of hollow crystals and a mesoporous (5-17 nm) shell of nanovesicle with macropores (about 350 nm) in the core. This hierarchical structure enables excellent Fe2O3 dispersion (3-4 nm) and resistance to sintering even at 800 °C; facilitates the transport of reactant/products; and exhibits superior activity and resistance to leaching in phenol degradation. Hollow nanovesicle assembly of Fe-Pt bimetal-encapsulated hollow ZSM-5 crystals was also prepared.