Search

As an emerging desalination technology, flow-electrode capacitive deionisation (FCDI) has the advantages of theoretically infinite adsorption capacity and applicability to high-concentration brine. However, during the operation of FCDI, the flow electrode in the S-shape channel is prone to sedimentation and clogging the channel. This undesirable phenomenon brings low efficiency and security issues. Therefore, a drop-shape channel was designed for FCDI to improve the flow regime of the flow electrode. The flow simulation of the drop-shape channel was performed to select the appropriate geometry to avoid the formation of the vortex and low-velocity region. The simulation results showed that the streamlined design of the drop-shape channel has insignificant velocity gradients. It significantly reduces the low-velocity region and improves the phenomenon of particle sedimentation. The desalination performance with varieties of electrode flow rate, AC content, and voltage was used to investigate the advantage between S-shape and drop-shape channels. It was found that under the conditions of low flow rate, high AC content, and high voltage, the drop-shape channel FCDI system could still obtain better ASRR and CE.

During the production of fiber-reinforced plastics using resin transfer molding (RTM), various characteristic defects and flaws can occur, such as fiber displacement and fiber waviness. Particularly in high-pressure RTM (HP-RTM), fiber misalignments are generated during infiltration by local peaks in the flow rate, leading to a significant reduction in the mechanical properties. To minimize or avoid this effect, the manufacturing process must be well controlled. Simulative approaches allow for a basic design of the mold filling process; however, due to the high number of influencing variables, the real behavior cannot be exactly reproduced. The focus of this work is on flow front monitoring in an HP-RTM mold using phased array ultrasonic testing. By using an established non-destructive testing instrument, the effort required for integration into the manufacturing process can be significantly reduced. For this purpose, investigations were carried out during the production of test specimens composed of glass fiber-reinforced polyurethane resin. Specifically, a phased array ultrasonic probe was used to record individual line scans over the form filling time. Taking into account the specifications of the probe used in these experiments, an area of 48.45 mm was inspected with a spatial resolution of 0.85 mm derived from the pitch. Due to the aperture that had to be applied to improve the signal-to-noise ratio, an averaging of the measured values similar to a moving average over a window of 6.8 mm had to be considered. By varying the orientation of the phased array probe and therefore the orientation of the line scans, it is possible to determine the local flow velocities of the matrix system during mold filling. Furthermore, process simulation studies with locally varying fiber volume contents were carried out. Despite the locally limited measuring range of the monitoring method presented, conclusions about the global flow behavior in a large mold can be drawn by comparing the experimentally determined results with the process simulation studies. The agreement between the measurement and simulation was thus improved by around 70%.

X-ray computed tomography was applied in imaging 3D-printed gyroids used for bioseparation in order to visualize and characterize structures from the entire geometry down to individual nanopores. Methacrylate prints were fabricated with feature sizes of 500 µm, 300 µm, and 200 µm, with the material phase exhibiting a porous substructure in all cases. Two X-ray scanners achieved pixel sizes from 5 µm to 16 nm to produce digital representations of samples across multiple length scales as the basis for geometric analysis and flow simulation. At the gyroid scale, imaged samples were visually compared to the original computed-aided designs to analyze printing fidelity across all feature sizes. An individual 500 µm feature, part of the overall gyroid structure, was compared and overlaid between design and imaged volumes, identifying individual printed layers. Internal subvolumes of all feature sizes were segmented into material and void phases for permeable flow analysis. Small pieces of 3D-printed material were optimized for nanotomographic imaging at a pixel size of 63 nm, with all three gyroid samples exhibiting similar geometric characteristics when measured. An average porosity of 45% was obtained that was within the expected design range, and a tortuosity factor of 2.52 was measured. Applying a voidage network map enabled the size, location, and connectivity of pores to be identified, obtaining an average pore size of 793 nm. Using Avizo XLAB at a bulk diffusivity of 7.00 × 10 -11 m 2 s -1 resulted in a simulated material diffusivity of 2.17 × 10 -11 m 2 s -1  ± 0.16 × 10 -11 m 2 s -1 ., (© 2023. The Author(s).)

The performance of a heat sink is significantly influenced by the type of cooling used: passive or active (forced), the shape of the heat sink, and the material from which it is made. This paper presents a review of the literature on the influence of geometry and surface parameters on effective heat transfer in heat sinks. The results of simulation studies for three different heat sink fin geometries and cooling types are presented. Furthermore, the influence of the surface texture of the heat sink fins on the heat transfer efficiency was determined. It was shown that the best performance in terms of geometries was that of a wave fin heat sink. When the surface texture was analyzed, it was found that an increase in the amplitude values of the texture decreases the heat dissipation efficiency in the case of active cooling, while for passive cooling, an increase in these parameters has a beneficial effect and increases the effective heat transfer to the surroundings. The cooling method was found to be the most important factor affecting heat dissipation efficiency. Forced airflow results in more efficient heat transfer from the heat sink fins to the surroundings.

This paper proposes a pollution diffusion model that accurately assesses changes in instantaneous river pollution in vegetation open channels. The model is established based on cellular automata and lattice Boltzmann method (LBM-CA). Flow influence coefficients are incorporated into cellular automata (CA) to represent the effect of vegetation on pollutant diffusion, while the lattice Boltzmann method (LBM) is utilized to simulate flow in vegetation open channels and obtain the flow influence coefficients for each cellular. The results show that the LBM-CA model has high accuracy and that pollutants tend to accumulate in vegetation areas, thereby extending the residence time of pollutants. The model incorporates pollution limits, allowing the prediction of basin pollution levels at specific times. The LBM-CA model provides a method for simulating pollutant diffusion in natural rivers., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Turbulent filling of molten metal in sand-casting leads to bi-films, porosity and oxide inclusions which results in poor mechanical properties and high scrap rate of sand castings. Hence, it is critical to understand the metal flow in sand-molds, i.e., casting hydrodynamics to eliminate casting defects. While multiple numerical methods have been applied to simulate this phenomenon for decades, harsh foundry environments and expensive x-ray equipment have limited experimental approaches to accurately visualize metal flow in sand molds. In this paper, a novel approach to solve this challenge is proposed using Succinonitrile (SCN) as a more accurate metal analog in place of water. SCN has a long history in solidification research due to its BCC (Body-Centered-Cubic) crystal structure and dendrite-like solidification (melting temperature ~60 °C) like molten aluminum. However, this is the first reported study on applying SCN through novel casting hydrodynamics to accurately visualize melt flow for casting studies. This paper used numerical simulations and experiments using both water and SCN to identify the critical dimensionless numbers to perform accurate metal flow analog testing. Froude's number and wall roughness were identified as critical variables. Experimental results show that SCN flow testing was more accurate in recreating the flow profile of molten aluminum, thus validating its utility as a metal analog for metal flow research. Findings from this study can be used in future metal flow analysis such as: runner, in-gate and integrated filling-feeding-solidification studies.

Complex hydraulic fracture networks are critical for enhancing permeability in unconventional reservoirs and mining industries. However, accurately simulating the fluid flow in realistic fracture networks (compared to the statistical fracture networks) is still challenging due to the fracture complexity and computational burden. This work proposes a simple yet efficient numerical framework for the flow simulation in fractured porous media obtained by 3D high-resolution images, aiming at both computational accuracy and efficiency. The fractured rock with complex fracture geometries is numerically constructed with a cell-based discrete fracture-matrix model (DFM) having implicit fracture apertures. The flow in the complex fractured porous media (including matrix flow, fracture flow, as well as exchange flow) is simulated with a pipe-based cell-centered finite volume method. The performance of this model is validated against analytical/numerical solutions. Then a lab-scale true triaxial hydraulically fractured shale sample is reconstructed, and the fluid flow in this realistic fracture network is simulated. Results suggest that the proposed method achieves a good balance between computational efficiency and accuracy. The complex fracture networks control the fluid flow process, and the opened natural fractures behave as primary fluid pathways. Heterogeneous and anisotropic features of fluid flow are well captured with the present model., Competing Interests: Conflict of interestThe authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© The Author(s) 2023.)

Introduction: In the circulatory system, the vessel branching angle may have hemodynamic consequences. We hypothesized that there is a hemodynamically optimal range for the renal artery's branching angle., Methods: Data on the posttransplant kinetics of estimated glomerular filtration rate (eGFR) were analyzed according to the donor and implant sides (right-to-right and left-to-right position; n = 46). The renal artery branching angle from the aorta of a randomly selected population was measured using an X-ray angiogram (n = 44). Computational fluid dynamics simulation was used to elucidate the hemodynamic effects of angulation., Results and Discussion: Renal transplant patients receiving a right donor kidney to the right side showed faster adaptation and higher eGFR values than those receiving a left donor kidney to the right side (eGFR: 65 ± 7 vs. 56 ± 6 mL/min/1.73 m2; p < 0.01). The average branching angle on the left side was 78° and that on the right side was 66°. Simulation results showed that the pressure, volume flow, and velocity were relatively constant between 58° and 88°, indicating that this range is optimal for the kidneys. The turbulent kinetic energy does not change significantly between 58° and 78°., Conclusion: The results suggest that there is an optimal range for the renal artery's branching angle from the aorta where hemodynamic vulnerability caused by the degree of angulation is the lowest, which should be considered during kidney transplantations., (© 2023 The Author(s). Published by S. Karger AG, Basel.)

The cerebrovasculature and its mural cells must meet brain regional energy demands, but how their spatial relationship with different neuronal cell types varies across the brain remains largely unknown. Here we apply brain-wide mapping methods to comprehensively define the quantitative relationships between the cerebrovasculature, capillary pericytes, and glutamatergic and GABAergic neurons, including neuronal nitric oxide synthase-positive (nNOS + ) neurons and their subtypes in adult mice. Our results show high densities of vasculature with high fluid conductance and capillary pericytes in primary motor sensory cortices compared with association cortices that show significant positive and negative correlations with energy-demanding parvalbumin + and vasomotor nNOS + neurons, respectively. Thalamo-striatal areas that are connected to primary motor sensory cortices also show high densities of vasculature and pericytes, suggesting dense energy support for motor sensory processing areas. Our cellular-resolution resource offers opportunities to examine spatial relationships between the cerebrovascular network and neuronal cell composition in largely understudied subcortical areas., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

Continuous flow injection is a key technology for serial crystallography measurements of protein crystals suspended in the lipidic cubic phase (LCP). To date, there has been little discussion in the literature regarding the impact of the injection process itself on the structure of the lipidic phase. This is despite the fact that the phase of the injection matrix is critical for the flow properties of the stream and potentially for sample stability. Here we report small-angle X-ray scattering measurements of a monoolein:water mixture during continuous delivery using a high viscosity injector. We observe both an alignment and modification of the LCP as a direct result of the injection process. The orientation of the cubic lattice with respect to the beam was estimated based on the anisotropy of the diffraction pattern and does not correspond to a single low order zone axis. The solvent fraction was also observed to impact the stability of the cubic phase during injection. In addition, depending on the distance traveled by the lipid after exiting the needle, the phase is observed to transition from a pure diamond phase (Pn3m) to a mixture containing both gyriod (Ia3d) and lamellar (L α ) phases. Finite element modelling of the observed phase behaviour during injection indicates that the pressure exerted on the lipid stream during extrusion accounts for the variations in the phase composition of the monoolein:water mixture., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Groundwater has become increasingly vulnerable to quality degradation. An elaborate understanding of its flow, draft, recharge and pollutant transport processes needs to be developed to understand its risk to contamination. This paper has discussed different tools and methods that are used to map groundwater vulnerability around the world. To maintain the quality and impact of the study, rigorous search for relevant literature published in high impact scientific journals has been done, and the comprehensive information on groundwater vulnerability assessment methods being used, has been compiled. The GIS based overlay and index-based methods like DRASTIC, GALDIT, GOD, COP and PI takes into consideration various thematic layers, overlays them to calculate weighted index and identifies vulnerability classes. They have been criticised for the lack of numerical basis in their formulation. Therefore, over the years, many of the proposed indices have been modified to provide quantitative estimates of groundwater potential to degrade and deplete. However, where the data and software are not a constraint, the use of numerical based simulation models can be done for more elaborate and numerical based quantification of the vulnerability. These numerical models typically require extensive data and are exceedingly becoming more sophisticated with the introduction of new parameters. This study concludes that integrating the GIS with numerical models offers the advantage of data management and assists to spatially analyse the datasets. The difficulties that are associated with the differences between GIS and numerical model's data structures should be thoroughly understood, prior to coupling, to develop uniform conversion software., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Over recent decades, tremendous advances in the field of scalable numerical tools and mesh immersion techniques have been achieved to improve numerical efficiency while preserving a good quality of the obtained results. In this context, an octree-optimized microstructure generation and domain reconstruction with adaptative meshing is presented and illustrated through a flow simulation example applied to permeability computation of micrometric fibrous materials. Thanks to the octree implementation, the numerous distance calculations in these processes are decreased, thus the computational complexity is reduced. Using the parallel environment of the ICI-tech library as a mesher and a solver, a large scale case study is performed. The study is applied to the computation of the full permeability tensor of a three-dimensional microstructure containing 10,000 fibers. The considered flow is a Stokes flow and it is solved with a stabilized finite element formulation and a monolithic approach.

In order to assess the limits and applicability of Pitot tubes for the measurement of flow velocity in narrow ducts, e.g., biomass burning plants, an optical, dual function device was implemented. This sensor, based on spectroscopic techniques, targets a trace gas, injected inside the stack either in bursts, or continuously, so performing transit time or dilution measurements. A comparison of the two optical techniques with respect to Pitot readings was carried out in different flow conditions (speed, temperature, gas composition). The results of the two optical measurements are in agreement with each other and fit quite well the theoretical simulation of the flow field, while the results of the Pitot measurements show a remarkable dependence on position and inclination of the Pitot tube with respect to the duct axis. The implications for the metrology of small combustors' emissions are outlined.

The LabPET II detection module is the building block of PET scanners for ultra-high-resolution imaging of small to mid-sized animals and the human brain. For optimal performance, it must be operated at a stable temperature. The detection module is composed of four APD-LYSO detector arrays with two flip-chip ASICs mounted on the backside of an interposer generating 550 mW each. Currently, the scanner architecture includes an air cavity around the electronics and smaller cavities close to the detectors. Cooling down the front-end electronics located in these small cavities becomes problematic as the number of modules increases to address the different targeted configurations of the LabPET II scanners from mouse to human brain geometries. A basic knowledge of the heat distribution is necessary to develop an efficient thermal management in all cases. The aim of this work is to build a model of the LabPET II ASIC and associated PCB for enabling heat flow simulations and circumscribe the thermal management requirements. The Flow Simulation module (SolidWorks), was used to build the thermal model. The ASIC and the interconnection with the PCB were reproduced accurately while some adjacent structures were simplified to ease the simulation burden. The model was applied to simulate three different configurations of printed-circuit boards carrying the ASICs and other components where a fan is turned on/off to create a forced airflow. Each simulation was compared to some experimental measurements. A temperature difference of less than 5 degree Celsius between the simulations and experimental measurements is noticed, giving confidence that the thermal model of the ASIC is valid and transferable to different mechanical assemblies.

Chemically driven colloidal motors capable of implementing different movements under a common environment are of great importance for various complex tasks. However, the key parameters underlying different motion behaviors are incompletely understood. Here, we demonstrate that carbonaceous nanoflask (CNF) motors move spontaneously in glucose powered by the cascade reaction of glucose oxidase and catalase, and their directional propulsion can be premeditated by controlling the surface wettability of nanomotors. The hydrophilic CNF motors move from the round-bottom to the opening neck (backward), whereas the hydrophobic CNF motors swim from the opening neck to the round-bottom (forward). We demonstrate that the backward motion of the hydrophilic CNF motors is driven by the local glucose gradient due to self-diffusiophoresis, and the forward movement of the hydrophobic CNF motors is caused by the locally produced glucose acid gradient. The fluid simulation reveals that the hydrophilic and hydrophobic CNF motors correspond to the puller and pusher models, respectively. Our study offers a minimal strategy to manipulate the direction of motion of motors for specific applications and to change the hydrodynamic behaviors of glucose-powered motors.

We analyse a relationship between the bulk microstructure of randomly packed beds, which we quantify through chord length distribution (CLD) analysis of the interparticle void space, and the associated flow heterogeneity, as expressed by the longitudinal dispersion coefficient at a Péclet number of Pe = 10. A random collection of physically reconstructed packings is complemented with a systematic set of computer-generated packings of monosized spheres, for which the packing-generation algorithm has been carefully adjusted to realize a monotonic variation of the bed porosity and microstructural heterogeneity. The most relevant difference in the morphology between these computer-generated and the physically reconstructed packings are structural defects present in the real packings, such as particle oligomers and larger voids as well as contaminations and particle debris. These defects influence the pore space morphology and introduce additional structural heterogeneity. Hydrodynamic dispersion coefficients for all packings are derived by implementing the lattice-Boltzmann method to simulate fluid flow and a random-walk particle tracking technique to record the transport of passive, point-like tracers in the flow fields. We propose a morphological descriptor, σ/μ, based on statistical parameters of a CLD (standard deviation σ and mean chord length μ) that can be used to predict the dispersion coefficient in packed beds, independent from the underlying particle size distribution, packing-generation protocol, bed porosity, and the occurrence of structural defects., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Purpose: Myocardial blood-flow simulation based on laws of fluid mechanics is a valuable tool for understanding tissue behavior. Our aim is to evaluate the ability of a porous-media flow model approach to reflect disturbed washout of contrast media (CM) from the myocardium as observed by cardiovascular MR., Methods: A coupled advection-diffusion model is used to describe the CM flow in the vascular and extravascular space as separate compartments. Their exchange of CM is controlled by the exchange rate ExR , which in turn determines the washout behavior. We fitted simulations to CM concentration measurements, derived from T 1 maps of the midventricular slice. The CM concentration was extracted from 18 patients with myocarditis in the acute phase and during follow-up after 6 months. The results were compared with 18 sex- and age-matched controls. For each subject, the measurements were acquired before and during the first 10 minutes at 5 time points after CM administration, representing CM washout. Image registration was applied to compensate for motion between different time points., Results: Eight matched data sets had to be excluded due to low registration quality. Processing was successful in n = 10 matched data sets of acute and healed myocarditis as well as controls. Significant differences in ExR were observed when comparing patients with acute myocarditis to controls (P < .001), to their follow-up (P < .05), and the follow-up to controls (P < .05)., Conclusion: Our study suggests the feasibility of using the proposed porous-medium flow framework for the simulation of pathologic myocardial tissue., (© 2019 International Society for Magnetic Resonance in Medicine.)

A better understanding of airflow characteristics in the upper airway (UA) is crucial in investigating obstructive sleep apnea (OSA), particle sedimentation, drug delivery, and many biomedical problems. Direct visualization of air flow patterns in in-vitro models with realistic anatomical structures is a big challenge. In this study, we constructed unique half-side transparent physical models of normal UA based on realistic anatomical structures. A smoke-wire method was developed to visualize the air flow in UA models directly. The results revealed that the airflow through the pharynx was laminar but not turbulent under normal inspiration, which suggested that compared with turbulent models, a laminar model should be more suitable in numerical simulations. The flow predicted numerically using the laminar model was consistent with the observations in the physical models. The comparison of the velocity fields predicted numerically using the half-side and complete models confirmed that it was reasonable to investigate the flow behaviors in UA using the half-side model. Using the laminar model, we simulated the flow and evaluated the effects of UA narrowing caused by rostral fluid shift on pharyngeal resistance. The results suggested that fluid shift could play an important role in the formation of hypopnea or OSA during sleep.

One challenge for hepatic flow simulation is to divide the hepatic vasculature into individual Couinaud segments, and to simulate flow at both segmental and organ levels. We propose to integrate a segment simulation algorithm with the flow solver in a Constructive Constraint Optimisation (CCO) algorithm to address this problem. In this way blood flow simulations can be conducted for large segment-specific vasculatures as relevant to surgical procedures. In this short communication we outline the methods and present some preliminary results.

Objective: Theoretical models are essential tools for studying microcirculatory function. Recently, the validity of a well-established phase separation model was questioned and it was claimed that it produces problematically low hematocrit predictions and lack of red cells in small diameter vessels. We conducted a quantitative evaluation of this phase separation model to establish common ground for future research., Methods: Model predictions were validated against a comprehensive database with measurements from 4 mesenteric networks. A Bayesian data analysis framework was used to integrate measurements and network model simulations into a combined analysis and to model uncertainties related to network boundary conditions as well as phase separation model parameters. The model evaluation was conducted within a cross-validation scheme., Results: Unlike the recently reported results, our analysis demonstrated good correspondence in global characteristics between measurements and predictions. In particular, predicted hematocrits for vessels with small diameters were consistent with measurements. Incorporating phase separation model parameter uncertainties further reduced the hematocrit validation error by 17% and led to the absence of red-cell-free segments. Corresponding model parameters are presented as alternatives to standard parameters., Conclusions: Consistent with earlier studies, our quantitative model evaluation supports the continued use of the established phase separation model., (© 2018 John Wiley & Sons Ltd.)

Purpose: The purpose of this study is to develop and evaluate a functionally personalized boundary condition (BC) model for estimating the fractional flow reserve (FFR) from coronary computed tomography angiography (CCTA) using flow simulation (CT-FFR)., Materials and Methods: The CCTA data of 90 subjects with subsequent invasive FFR in 123 lesions within 21 days (range: 0-83) were retrospectively collected. We developed a functionally personalized BC model accounting specifically for the coronary microvascular resistance dependency on the coronary outlets pressure suggested by several physiological studies. We used the proposed model to estimate the hemodynamic significance of coronary lesions with an open-loop physics-based flow simulation. We generated three-dimensional (3D) coronary tree geometries using automatic software and corrected manually where required. We evaluated the improvement in CT-FFR estimates achieved using a functionally personalized BC model over anatomically personalized BC model using k-fold cross-validation., Results: The functionally personalized BC model slightly improved CT-FFR specificity in determining hemodynamic significance of lesions with intermediate diameter stenosis (30%-70%, N = 72), compared to the anatomically personalized model lesions with invasive FFR measurements as the reference (sensitivity/specificity: 0.882/0.79 vs 0.882/0.763). For the entire set of 123 coronary lesions, the functionally personalized BC model improved only the area under the curve (AUC) but not the sensitivity/specificity in determining the hemodynamic significance of lesions, compared to the anatomically personalized model (AUC: 0.884 vs 0.875, sensitivity/specificity: 0.848/0.805)., Conclusion: The functionally personalized BC model has the potential to improve the quality of CT-FFR estimates compared to an anatomically personalized BC model., (© 2018 American Association of Physicists in Medicine.)

Currently there are no available methods for prediction of thrombotic complications in Coronary Artery disease. Additionally, blood coagulation tests are mainly performed in a steady system while coagulation in vivo occurs under flow conditions. In this work, a phenomenological model for coagulation up-to thrombin generation is proposed; the model is mainly based on the results of thrombin generation assays and therefore it can account for the variation of the coagulability that is observed in different individuals. The model is applied on 3 cases of left anterior descending arteries (LAD) with 50% maximum stenosis placed at a different location and have been statistically assessed as of different complication risk. The simulations showed that parameters of thrombin generation assays obtain different values when they refer to thrombin generation under realistic coronary flow conditions. The flow conditions prevailing locally because of the geometric differences among the arterial trees can lead to different initiation times and thrombin production rates and it also alters the spatial distribution of the coagulation products. Similarly, small changes of the coagulation characteristics of blood under identical flow conditions can allow or prevent the initiation of coagulation. The results indicate that combined consideration of geometry and coagulation characteristics of blood can lead to entirely different conclusions compared to independent assessment of each factor., (Copyright © 2017 IPEM. Published by Elsevier Ltd. All rights reserved.)

Blood flow imaging is a novel technology in cardiovascular medicine and surgery. Today, two types of blood flow imaging tools are available: measurement-based flow visualization including 4D flow MRI (or 3D cine phase-contrast magnetic resonance imaging), or echocardiography flow visualization software, and computer flow simulation modeling based on computational fluid dynamics (CFD). MRI and echocardiography flow visualization provide measured blood flow but have limitations in temporal and spatial resolution, whereas CFD flow calculates the flow according to assumptions instead of flow measurement, and it has sufficiently fine resolution up to the computer memory limit, and it enables even virtual surgery when combined with computer graphics. Blood flow imaging provides profound insight into the pathophysiology of cardiovascular diseases, because it quantifies and visualizes mechanical stress on the vessel walls or heart ventricle. Wall shear stress (WSS) is a stress on the endothelial wall caused by the near wall blood flow, and it is thought to be a predictor of atherosclerosis progression in coronary or aortic diseases. Flow energy loss (EL) is the loss of blood flow energy caused by viscous friction of turbulent diseased flow, and it is expected to be a predictor of ventricular workload on various heart diseases including heart valve disease, cardiomyopathy, and congenital heart diseases. Blood flow imaging can provide useful information for developing predictive medicine in cardiovascular diseases, and may lead to breakthroughs in cardiovascular surgery, especially in the decision-making process.

There is higher long-term failure of the saphenous vein graft (SVG) compared with the left internal mammary artery (LIMA) graft, which is affected by the hemodynamic environment. A comprehensive analysis of postoperative structure-function changes is important to study the atherogenesis in the SVG A comparison of morphometric and hemodynamic parameters was carried out between LIMA grafts and SVGs and between different time points postoperatively. Various parameters were obtained from the image reconstruction and flow simulation in patients, who underwent CT exams for ~1 year, 5 and 10 years after revascularization. Morphometric data showed a decrease in lumen size in the entire SVG and anastomosis of different patients in a sequence of ~1 year, 5 and 10 years postoperatively despite negligible changes of LIMA size. Computational results indicated the fourfold increased surface area ratio (SAR) of low time-averaged wall shear stress (TAWSS) in the SVG and anastomosis at postoperative 10 years than that at postoperative ~1 year. The SAR of high TAWSS gradient (TAWSSG) at the distal anastomosis between SVG and coronary arteries was significantly higher (14 ± 9% vs. 6 ± 8%) than that in the LIMA group at postoperative ~1 year. There were strong correlations between morphometric and hemodynamic parameters in the SVG and distal anastomosis at various time points postoperatively, which showed deterioration relevant to persistent diffuse diseases at postoperative ~10 years., (© 2017 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of The Physiological Society and the American Physiological Society.)

Column wall effects are well recognized as major limiting factor in achieving high separation efficiency in HPLC. This is especially important for modern analytical columns packed with small particles, where wall effects dominate the band broadening. Detailed knowledge about the packing microstructure of packed analytical columns has so far not been acquired. Here, we present the first three-dimensional reconstruction protocol for these columns utilizing focused ion-beam scanning electron microscopy (FIB-SEM) on a commercial 2.1mm inner diameter×50mm length narrow-bore analytical column packed with 1.7μm bridged-ethyl hybrid silica particles. Two sections from the packed bed are chosen for reconstruction by FIB-SEM: one from the bulk packing region of the column and one from its critical wall region. This allows quantification of structural differences between the wall region and the center of the bed due to effects induced by the hard, confining column wall. Consequences of these effects on local flow velocity in the column are analyzed with flow simulations utilizing the lattice-Boltzmann method. The reconstructions of the bed structures reveal significant structural differences in the wall region (extending radially over approximately 62 particle diameters) compared to the center of the column. It includes the local reduction of the external porosity by up to 10% and an increase of the mean particle diameter by up to 3%, resulting in a decrease of the local flow velocity by up to 23%. In addition, four (more ordered) layers of particles in the direct vicinity of the column wall induce local velocity fluctuations by up to a factor of three regarding the involved velocity amplitudes. These observations highlight the impact of radial variations in packing microstructure on band migration and column performance. This knowledge on morphological peculiarities of column wall effects helps guiding us towards further optimization of the packing process for analytical HPLC columns., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Objective: In vivo imaging of the microcirculation and network-oriented modeling have emerged as powerful means of studying microvascular function and understanding its physiological significance. Network-oriented modeling may provide the means of summarizing vast amounts of data produced by high-throughput imaging techniques in terms of key, physiological indices. To estimate such indices with sufficient certainty, however, network-oriented analysis must be robust to the inevitable presence of uncertainty due to measurement errors as well as model errors., Methods: We propose the Bayesian probabilistic data analysis framework as a means of integrating experimental measurements and network model simulations into a combined and statistically coherent analysis. The framework naturally handles noisy measurements and provides posterior distributions of model parameters as well as physiological indices associated with uncertainty., Results: We applied the analysis framework to experimental data from three rat mesentery networks and one mouse brain cortex network. We inferred distributions for more than 500 unknown pressure and hematocrit boundary conditions. Model predictions were consistent with previous analyses, and remained robust when measurements were omitted from model calibration., Conclusion: Our Bayesian probabilistic approach may be suitable for optimizing data acquisition and for analyzing and reporting large data sets acquired as part of microvascular imaging studies., (© 2016 John Wiley & Sons Ltd.)

The geometry of the coronary vessel network is believed to play a decisive role in the initiation, progression and outcome of coronary artery disease (CAD) and the occurrence of acute coronary syndromes (ACS). It also determines the flow field in the coronary artery which can be linked to CAD evolution. In this work geometric 3D models of left anterior descending (LAD) coronary arteries associated with either myocardial infarction (MI) or stable (STA) CAD were constructed. Transient numerical simulations of the flow for each model showed that specific flow patterns develop in different extent in the different groups examined. Recirculation zones, present distal the stenosis in all models, had larger extent and duration in MI cases. For mild stenosis (up to 50%) areas with low time averaged wall shear stress TAWSS (<0.15Pa) as well as areas with high TAWSS (>3Pa) appeared only in MI models; in moderate and severe stenosis (>50%) these areas were present in all models but were significantly larger for MI than STA models. These differentiations were expressed via numerical indices based on TAWSS, oscillating shear index (OSI) and relative residence time (RRT). Additionally we introduced the coagulation activation index (CAI), based on the threshold behaviour of coagulation initiation, which exceeded the suggested threshold only for MI models with intermediate stenosis (up to 50%). These results show that numerical simulations of flow can produce arithmetic indices linked with the risk of CAD complications., (Copyright © 2016 IPEM. Published by Elsevier Ltd. All rights reserved.)

The ability of clothing to provide protection against external environments is critical for wearer's safety and thermal comfort. It is a function of several factors, such as external environmental conditions, clothing properties and activity level. These factors determine the characteristics of the different microclimates existing inside the clothing which, ultimately, have a key role in the transport processes occurring across clothing. As an effort to understand the effect of transport phenomena in clothing microclimates on the overall heat transport across clothing structures, a numerical approach was used to study the buoyancy-driven heat transfer across horizontal air layers trapped inside air impermeable clothing. The study included both the internal flow occurring inside the microclimate and the external flow occurring outside the clothing layer, in order to analyze the interdependency of these flows in the way heat is transported to/from the body. Two-dimensional simulations were conducted considering different values of microclimate thickness (8, 25 and 52 mm), external air temperature (10, 20 and 30 °C), external air velocity (0.5, 1 and 3 m s(-1)) and emissivity of the clothing inner surface (0.05 and 0.95), which implied Rayleigh numbers in the microclimate spanning 4 orders of magnitude (9 × 10(2)-3 × 10(5)). The convective heat transfer coefficients obtained along the clothing were found to strongly depend on the transport phenomena in the microclimate, in particular when natural convection is the most important transport mechanism. In such scenario, convective coefficients were found to vary in wavy-like manner, depending on the position of the flow vortices in the microclimate. These observations clearly differ from data in the literature for the case of air flow over flat-heated surfaces with constant temperature (which shows monotonic variations of the convective heat transfer coefficients, along the length of the surface). The flow patterns and temperature fields in the microclimates were found to strongly depend on the characteristics of the external boundary layer forming along the clothing and on the distribution of temperature in the clothing. The local heat transfer rates obtained in the microclimate are in marked contrast with those found in the literature for enclosures with constant-temperature active walls. These results stress the importance of coupling the calculation of the internal and the external flows and of the heat transfer convective and radiative components, when analyzing the way heat is transported to/from the body.

Objectives: Single ventricle with apicocaval juxtaposition (ACJ) is a rare, complex anomaly, in which the optimal position of the conduit for completion of total cavopulmonary connection (TCPC) is still controversial. The purpose of this study was to identify a preoperative method for optimal conduit position using the IVC anatomy and computational fluid dynamics (CFD)., Methods: Twenty-four patients with ACJ (5.3 ± 5.7 years) who underwent TCPC were enrolled. A conduit was placed ipsilateral to the cardiac apex in each of 11 patients, of which 9 were intra-atrial and 2 extracardiac (group A) and, in a further 13 patients, extracardiac on the contralateral side (group B). As control, 10 patients with tricuspid atresia were also enrolled (group C). The location of the IVC in relation to the spine was evaluated from the frontal view of preoperative angiogram, using the following index: IVC-index = IVC width overlapping the vertebra/width of the vertebra × 100%. Energy loss was calculated by CFD simulation., Results: IVC-index of group B was larger than groups A and C (45 ± 26 vs. 20 ± 21 and 28 ± 19%, P = 0.03). Postoperative catheterizations showed that, due to its curvature, conduit length in group B was significantly longer than the others (65 ± 12 vs. 36 ± 14 and 44 ± 10 mm, P < 0.001), although there was no statistical difference in central venous pressure or cardiac output. CFD studies revealed less energy loss in group A conduits compared with group B (1.6 ± 0.3 vs. 3.6 ± 0.6 mW, P = 0.05), although this did not appear to be clinically significant. Moreover, CFD simulation showed significant energy loss within the Fontan circulation when the conduit was either compressed or kinked: 4.9 and 18.2 mW respectively., Conclusions: In patients with ACJ, placement of a straighter and shorter conduit on the ventricular apical side provides better laminar blood flow with less energy loss. However, conduit compression and kinking are far more detrimental to the Fontan circulation. A preoperative IVC-index is pivotal for avoiding these factors and deciding the optimal conduit route.