Numerical simulation of blood flow and pressure drop in the pulmonary arterial and venous circulation

A novel multiscale mathematical and computational model of the pulmonary circulation is presented and used to analyse both arterial and venous pressure and flow. This work is a major advance over previous studies by Olufsen et al. (Ann Biomed Eng 28:1281–1299, 2012) which only considered the arterial circulation. For the first three generations of vessels within the pulmonary circulation, geometry is specified from patient-specific measurements obtained using magnetic resonance imaging (MRI). Blood flow and pressure in the larger arteries and veins are predicted using a nonlinear, cross-sectional-area-averaged system of equations for a Newtonian fluid in an elastic tube. Inflow into the main pulmonary artery is obtained from MRI measurements, while pressure entering the left atrium from the main pulmonary vein is kept constant at the normal mean value of 2 mmHg. Each terminal vessel in the network of ‘large’ arteries is connected to its corresponding terminal vein via a network of vessels representing the vascular bed of smaller arteries and veins. We develop and implement an algorithm to calculate the admittance of each vascular bed, using bifurcating structured trees and recursion. The structured-tree models take into account the geometry and material properties of the ‘smaller’ arteries and veins of radii ≥ 50 μ m. We study the effects on flow and pressure associated with three classes of pulmonary hypertension expressed via stiffening of larger and smaller vessels, and vascular rarefaction. The results of simulating these pathological conditions are in agreement with clinical observations, showing that the model has potential for assisting with diagnosis and treatment for circulatory diseases within the lung

Fluids in cosmology

We review the role of fluids in cosmology by first introducing them in General Relativity and then by applying them to a FRW Universe's model. We describe how relativistic and non-relativistic components evolve in the background dynamics. We also introduce scalar fields to show that they are able to yield an inflationary dynamics at very early times (inflation) and late times (quintessence). Then, we proceed to study the thermodynamical properties of the fluids and, lastly, its perturbed kinematics. We make emphasis in the constrictions of parameters by recent cosmological probes.Comment: 34 pages, 4 figures, version accepted as invited review to the book "Computational and Experimental Fluid Mechanics with Applications to Physics, Engineering and the Environment". Version 2: typos corrected and references expande

Assessing animal affect: an automated and self-initiated judgement bias task based on natural investigative behaviour

Scientific methods for assessing animal affect, especially affective valence (positivity or negativity), allow us to evaluate animal welfare and the effectiveness of 3Rs Refinements designed to improve wellbeing. Judgement bias tasks measure valence; however, task-training may be lengthy and/or require significant time from researchers. Here we develop an automated and self-initiated judgement bias task for rats which capitalises on their natural investigative behaviour. Rats insert their noses into a food trough to start trials. They then hear a tone and learn either to stay for 2 s to receive a food reward or to withdraw promptly to avoid an air-puff. Which contingency applies is signalled by two different tones. Judgement bias is measured by responses to intermediate ambiguous tones. In two experiments we show that rats learn the task in fewer sessions than other automated variants, generalise responses across ambiguous tones as expected, self-initiate 4-5 trials/min, and can be tested repeatedly. Affect manipulations generate main effect trends in the predicted directions, although not localised to ambiguous tones, so further construct validation is required. We also find that tone-reinforcer pairings and reinforcement or non-reinforcement of ambiguous trials can affect responses to ambiguity. This translatable task should facilitate more widespread uptake of judgement bias testing

Matrix Development in Self-Assembly of Articular Cartilage

Articular cartilage is a highly functional tissue which covers the ends of long bones and serves to ensure proper joint movement. A tissue engineering approach that recapitulates the developmental characteristics of articular cartilage can be used to examine the maturation and degeneration of cartilage and produce fully functional neotissue replacements for diseased tissue.This study examined the development of articular cartilage neotissue within a self-assembling process in two phases. In the first phase, articular cartilage constructs were examined at 1, 4, 7, 10, 14, 28, 42, and 56 days immunohistochemically, histologically, and through biochemical analysis for total collagen and glycosaminoglycan (GAG) content. Based on statistical changes in GAG and collagen levels, four time points from the first phase (7, 14, 28, and 56 days) were chosen to carry into the second phase, where the constructs were studied in terms of their mechanical characteristics, relative amounts of collagen types II and VI, and specific GAG types (chondroitin 4-sulfate, chondroitin 6-sulfate, dermatan sulfate, and hyaluronan). Collagen type VI was present in initial abundance and then localized to a pericellular distribution at 4 wks. N-cadherin activity also spiked at early stages of neotissue development, suggesting that self-assembly is mediated through a minimization of free energy. The percentage of collagen type II to total collagen significantly increased over time, while the proportion of collagen type VI to total collagen decreased between 1 and 2 wks. The chondroitin 6- to 4- sulfate ratio decreased steadily during construct maturation. In addition, the compressive properties reached a plateau and tensile characteristics peaked at 4 wks.The indices of cartilage formation examined in this study suggest that tissue maturation in self-assembled articular cartilage mirrors known developmental processes for native tissue. In terms of tissue engineering, it is suggested that exogenous stimulation may be necessary after 4 wks to further augment the functionality of developing constructs

Transient integral boundary layer method to calculate the translesional pressure drop and the fractional flow reserve in myocardial bridges

BACKGROUND: The pressure drop – flow relations in myocardial bridges and the assessment of vascular heart disease via fractional flow reserve (FFR) have motivated many researchers the last decades. The aim of this study is to simulate several clinical conditions present in myocardial bridges to determine the flow reserve and consequently the clinical relevance of the disease. From a fluid mechanical point of view the pathophysiological situation in myocardial bridges involves fluid flow in a time dependent flow geometry, caused by contracting cardiac muscles overlying an intramural segment of the coronary artery. These flows mostly involve flow separation and secondary motions, which are difficult to calculate and analyse. METHODS: Because a three dimensional simulation of the haemodynamic conditions in myocardial bridges in a network of coronary arteries is time-consuming, we present a boundary layer model for the calculation of the pressure drop and flow separation. The approach is based on the assumption that the flow can be sufficiently well described by the interaction of an inviscid core and a viscous boundary layer. Under the assumption that the idealised flow through a constriction is given by near-equilibrium velocity profiles of the Falkner-Skan-Cooke (FSC) family, the evolution of the boundary layer is obtained by the simultaneous solution of the Falkner-Skan equation and the transient von-Kármán integral momentum equation. RESULTS: The model was used to investigate the relative importance of several physical parameters present in myocardial bridges. Results have been obtained for steady and unsteady flow through vessels with 0 – 85% diameter stenosis. We compare two clinical relevant cases of a myocardial bridge in the middle segment of the left anterior descending coronary artery (LAD). The pressure derived FFR of fixed and dynamic lesions has shown that the flow is less affected in the dynamic case, because the distal pressure partially recovers during re-opening of the vessel in diastole. We have further calculated the wall shear stress (WSS) distributions in addition to the location and length of the flow reversal zones in dependence on the severity of the disease. CONCLUSION: The described boundary layer method can be used to simulate frictional forces and wall shear stresses in the entrance region of vessels. Earlier models are supplemented by the viscous effects in a quasi three-dimensional vessel geometry with a prescribed wall motion. The results indicate that the translesional pressure drop and the mean FFR compares favourably to clinical findings in the literature. We have further shown that the mean FFR under the assumption of Hagen-Poiseuille flow is overestimated in developing flow conditions