A mathematical model of blood,cerebrospinal fluid and brain dynamics

Using first principles of fluid and solid mechanics a comprehensive model of human intracranial dynamics is proposed. Blood, cerebrospinal fluid (CSF) and brain parenchyma as well as the spinal canal are included. The compartmental model predicts intracranial pressure gradients, blood and CSF flows and displacements in normal and pathological conditions like communicating hydrocephalus. The system of differential equations of first principles conservation balances is discretized and solved numerically. Fluid–solid interactions of the brain parenchyma with cerebral blood and CSF are calculated. The model provides the transitions from normal dynamics to the diseased state during the onset of communicating hydrocephalus. Predicted results were compared with physiological data from Cine phase-contrast magnetic resonance imaging to verify the dynamic model. Bolus injections into the CSF are simulated in the model and found to agree with clinical measurements.

---

     

Computing the ankle-brachial index with parallel computational fluid dynamics

The ankle-brachial index (ABI), a ratio of arterial blood pressure in the ankles and upper arms, is used to diagnose and monitor circulatory conditions such as coarctation of the aorta and peripheral artery disease. Computational simulations of the ABI can potentially determine the parameters that produce an ABI indicative of ischemia or other abnormalities in blood flow. However, 0- and 1-D computational methods are limited in describing a 3-D patient-derived geometry. Thus, we present a massively parallel framework for computational fluid dynamics (CFD) simulations in the full arterial system. Using the lattice Boltzmann method to solve the Navier–Stokes equations, we employ highly parallelized and scalable methods to generate the simulation domain and efficiently distribute the computational load among processors. For the first time, we compute an ABI with 3-D CFD. In this proof-of-concept study, we investigate the dependence of ABI on the presence of stenoses, or narrowed regions of the arteries, by directly modifying the arterial geometry. As a result, our framework enables the computation a hemodynamic factor characterizing flow at the scale of the full arterial system, in a manner that is extensible to patient-specific imaging data and holds potential for treatment planning.     

Airflow behavior changes in upper airway caused by different head and neck positions: Comparison by computational fluid dynamics

The feasibility of computational fluid dynamics (CFD) to evaluate airflow characteristics in different head and neck positions has not been established. This study compared the changes in volume and airflow behavior of the upper airway by CFD simulation to predict the influence of anatomical and physiological airway changes due to different head–neck positions on mechanical ventilation. One awake volunteer with no risk of difficult airway underwent computed tomography in neutral position, extension position (both head and neck extended), and sniffing position (head extended and neck flexed). Three-dimensional airway models of the upper airway were reconstructed. The total volume (V) and narrowest area (A_min) of the airway models were measured. CFD simulation with an Spalart–Allmaras model was performed to characterize airflow behavior in neutral, extension, and sniffing positions of closed-mouth and open-mouth ventilation. The comparison result for V was neutral <extension≈sniffing, and for A_min was neutral<extension<sniffing. A_min in sniffing position was nearly 3.0 times that in neutral position and 1.7 times that in extension position. The pressure drop and velocity increasing were more obvious in neutral than sniffing or extension position at the same airflow rate. In sniffing position, pressure differences decreased and velocity remained almost constant. Recirculation airflow was generated near the subglottic region in neutral and extension positions. Sniffing position improves airway patency by increasing airway volume and decreasing airway resistance, suggesting that sniffing position may be the optimal choice for mask ventilation.     

Assessment of human left ventricle flow using statistical shape modelling and computational fluid dynamics

Blood flow patterns in the human left ventricle (LV) have shown relation to cardiac health. However, most studies in the literature are limited to a few patients and results are hard to generalize. This study aims to provide a new framework to generate more generalized insights into LV blood flow as a function of changes in anatomy and wall motion. In this framework, we studied the four-dimensional blood flow in LV via computational fluid dynamics (CFD) in conjunction with a statistical shape model (SSM), built from segmented LV shapes of 150 subjects. We validated results in an in-vitro dynamic phantom via time-resolved optical particle image velocimetry (PIV) measurements. This combination of CFD and the SSM may be useful for systematically assessing blood flow patterns in the LV as a function of varying anatomy and has the potential to provide valuable data for diagnosis of LV functionality.     

Estimation of valvular resistance of segmented aortic valves using computational fluid dynamics

Aortic valve stenosis is associated with an elevated left ventricular pressure and transaortic pressure drop. Clinicians routinely use Doppler ultrasound to quantify aortic valve stenosis severity by estimating this pressure drop from blood velocity. However, this method approximates the peak pressure drop, and is unable to quantify the partial pressure recovery distal to the valve. As pressure drops are flow dependent, it remains difficult to assess the true significance of a stenosis for low-flow low-gradient patients. Recent advances in segmentation techniques enable patient-specific Computational Fluid Dynamics (CFD) simulations of flow through the aortic valve. In this work a simulation framework is presented and used to analyze data of 18 patients. The ventricle and valve are reconstructed from 4D Computed Tomography imaging data. Ventricular motion is extracted from the medical images and used to model ventricular contraction and corresponding blood flow through the valve. Simplifications of the framework are assessed by introducing two simplified CFD models: a truncated time-dependent and a steady-state model. Model simplifications are justified for cases where the simulated pressure drop is above 10 mmHg. Furthermore, we propose a valve resistance index to quantify stenosis severity from simulation results. This index is compared to established metrics for clinical decision making, i.e. blood velocity and valve area. It is found that velocity measurements alone do not adequately reflect stenosis severity. This work demonstrates that combining 4D imaging data and CFD has the potential to provide a physiologically relevant diagnostic metric to quantify aortic valve stenosis severity.     

Functional imaging using computational fluid dynamics to predict treatment success of mandibular advancement devices in sleep-disordered breathing

Mandibular advancement devices (MADs) have emerged as a popular alternative for the treatment of sleep-disordered breathing. These devices bring the mandibula forward in order to increase upper airway (UA) volume and prevent total UA collapse during sleep. However, the precise mechanism of action appears to be quite complex and is not yet completely understood; this might explain interindividual variation in treatment success. We examined whether an UA model, that combines imaging techniques and computational fluid dynamics (CFD), allows for a prediction of the treatment outcome with MADs. Ten patients that were treated with a custom-made mandibular advancement device (MAD), underwent split-night polysomnography. The morning after the sleep study, a low radiation dose CT scan was scheduled with and without the MAD. The CT examinations allowed for a comparison between the change in UA volume and the anatomical characteristics through the conversion to three-dimensional computer models. Furthermore, the change in UA resistance could be calculated through flow simulations with CFD. Boundary conditions for the model such as mass flow rate and pressure distributions were obtained during the split-night polysomnography. Therefore, the flow modeling was based on a patient specific geometry and patient specific boundary conditions. The results indicated that a decrease in UA resistance and an increase in UA volume correlate with both a clinical and an objective improvement. The results of this pilot study suggest that the outcome of MAD treatment can be predicted using the described UA model.     

Analysis of thoracic aorta hemodynamics using 3D particle tracking velocimetry and computational fluid dynamics

Parallel to the massive use of image-based computational hemodynamics to study the complex flow establishing in the human aorta, the need for suitable experimental techniques and ad hoc cases for the validation and benchmarking of numerical codes has grown more and more.     

Computational fluid dynamics study of common stent models inside idealised curved coronary arteries

The haemodynamic behaviour of blood inside a coronary artery after stenting is greatly affected by individual stent features as well as complex geometrical properties of the artery including tortuosity and curvature. Regions at higher risk of restenosis, as measured by low wall shear stress (WSS < 0.5 Pa), have not yet been studied in detail in curved stented arteries. In this study, three-dimensional computational modelling and computational fluid dynamics methodologies were used to analyse the haemodynamic characteristics in curved stented arteries using several common stent models. Results in this study showed that stent strut thickness was one major factor influencing the distribution of WSS in curved arteries. Regions of low WSS were found behind struts, particularly those oriented at a large angle relative to the streamwise flow direction. These findings were similar to those obtained in studies of straight arteries. An uneven distribution of WSS at the inner and outer bends of curved arteries was observed where the WSS was lower at the inner bend. In this study, it was also shown that stents with a helical configuration generated an extra swirling component of the flow based on the helical direction; however, this extra swirl in the flow field did not cause significant changes on the distribution of WSS under the current setup.     

Integration of microbial kinetics and fluid dynamics toward model‐driven scale‐up of industrial bioprocesses

Scale‐up of bioprocesses is hampered by open questions, mostly related to poor mixing and mass transfer limitations. Concentration gradients of substrate, carbon dioxide, and oxygen in time and space, especially in large‐scale high‐cell density fed‐batch processes, are likely induced as the mixing time of the fermentor is usually longer than the relevant cellular reaction time. Cells in the fermentor are therefore repeatedly exposed to dynamic environments or perturbations. As a consequence, the heterogeneity in industrial practices often decreases either yield, titer, or productivity, or combinations thereof and increases by‐product formation as compared to well‐mixed small‐scale bioreactors, which is summarized as scale‐up effects. Identification of response mechanisms of the microorganism to various external perturbations is of great importance for pinpointing metabolic bottlenecks and targets for metabolic engineering. In this review, pulse response experimentation is proposed as an ideal way of obtaining kinetic information in combination with scale‐down approaches for in‐depth understanding of dynamic response mechanisms. As an emerging tool, computational fluid dynamics is able to draw a holistic picture of the fluid flow and concentration fields in the fermentor and finds its use in the optimization of fermentor design and process strategy. In the future, directed strain improvement and fermentor redesign are expected to largely depend on models, in which both microbial kinetics and fluid dynamics are thoroughly integrated.     

Shoulder joint kinetics and dynamics during underwater forward arm elevation

Aquatic exercises are widely implemented into rehabilitation programs. However, both evaluating their mechanical demands on the musculoskeletal system and designing protocols to provide progressive loading are difficult tasks. This study reports for the first time shoulder joint kinetics and dynamics during underwater forward arm elevation performed at speeds ranging from 22.5 to 90°/s. Net joint moments projected onto anatomical axes of rotation, joint power, and joint work were calculated in 18 participants through a novel approach coupling numerical fluid flow simulations and inverse dynamics. Joint dynamics was revealed from the 3D angle between the joint moment and angular velocity vectors, identifying three main functions—propulsion, stabilization, and resistance. Speeds <30°/s necessitated little to no power at all, whereas peaks about 0.20 W⋅kg⁻¹ were seen at 90°/s. As speed increased, peak moments were up to 61 × higher at 90 than at 22.5°/s, (1.82 ± 0.12%BW⋅AL vs 0.03 ± 0.01%BW⋅AL, P < 0.038). This was done at the expense of a substantial decrease in the joint moment contribution to joint stability though, which goes against the intuition that greater stabilization is required to protect the shoulder from increasing loads. Slow arm elevations (<30°/s) are advantageous for joint mobility gain at low mechanical solicitation, whereas the intensity at 90°/s is high enough to stimulate muscular endurance improvements. Simple predictive equations of shoulder mechanical loading are provided. They allow for easy design of progressive protocols, either for the postoperative shoulder or the conditioning of athlete targeting very specific intensity regions.     

Flow in the well: computational fluid dynamics is essential in flow chamber construction

A perfusion system was developed to generate well defined flow conditions within a well of a standard multidish. Human vein endothelial cells were cultured under flow conditions and cell response was analyzed by microscopy. Endothelial cells became elongated and spindle shaped. As demonstrated by computational fluid dynamics (CFD), cells were cultured under well defined but time varying shear stress conditions. A damper system was introduced which reduced pulsatile flow when using volumetric pumps. The flow and the wall shear stress distribution were analyzed by CFD for the steady and unsteady flow field. Usage of the volumetric pump caused variations of the wall shear stresses despite the controlled fluid environment and introduction of a damper system. Therefore the use of CFD analysis and experimental validation is critical in developing flow chambers and studying cell response to shear stress. The system presented gives an effortless flow chamber setup within a 6-well standard multidish.     

A multiscale 0-D/3-D approach to patient-specific adaptation of a cerebral autoregulation model for computational fluid dynamics studies of cardiopulmonary bypass

Neurological complication often occurs during cardiopulmonary bypass (CPB). One of the main causes is hypoperfusion of the cerebral tissue affected by the position of the cannula tip and diminished cerebral autoregulation (CA). Recently, a lumped parameter approach could describe the baroreflex, one of the main mechanisms of cerebral autoregulation, in a computational fluid dynamics (CFD) study of CPB. However, the cerebral blood flow (CBF) was overestimated and the physiological meaning of the variables and their impact on the model was unknown. In this study, we use a 0-D control circuit representation of the Baroreflex mechanism, to assess the parameters with respect to their physiological meaning and their influence on CBF. Afterwards the parameters are transferred to 3D-CFD and the static and dynamic behavior of cerebral autoregulation is investigated.     

Torsional motion of the left ventricle does not affect ventricular fluid dynamics of both foetal and adult hearts

Left ventricular torsion is caused by shortening and relaxation of the helical fibres in the myocardium, and is thought to be an optimal configuration for minimizing myocardial tissue strains. Characteristics of torsional motion has also been proposed to be markers for cardiac dysfunction. However, its effects on fluid and energy dynamics in the left ventricle have not been comprehensively investigated. To investigate this, we performed image-based flow simulations on five healthy adult porcine and two healthy human foetal left ventricles (representing two different length scales) at different degrees of torsional motions. In the adult porcine ventricles, cardiac features such as papillary muscles and mitral valves, and cardiac conditions such as myocardial infarctions, were also included to investigate the effect of twist. The results showed that, for all conditions investigated, ventricular torsional motion caused minimal changes to flow patterns, and consistently accounted for less than 2% of the energy losses, wall shear stresses, and ejection momentum energy. In contrast, physiological characteristics such as chamber size, stroke volume and heart rate had a much greater influence on flow patterns and energy dynamics. The results thus suggested that it might not be necessary to model the torsional motion to study the flow and energy dynamics in left ventricles.     

Characterizing the fluid dynamics in the flow fields of cylindrical orbitally shaken bioreactors with different geometry sizes

Orbitally shaken bioreactors (OSRs) are commonly used for the cultivation of mammalian cells in suspension. To aid the geometry designing and optimizing of OSRs, we conducted a three‐dimensional computational fluid dynamics (CFD) simulation to characterize the flow fields in a 10 L cylindrical OSR with different vessel diameters. The liquid wave shape captured by a camera experimentally validated the CFD models established for the cylindrical OSR. The geometry size effect on volumetric mass transfer coefficient (k_La) and hydromechanical stress was analyzed by varying the ratio of vessel diameter (d) to liquid height at static (h_L), d/h_L. The highest value of k_La about 30 h^?1 was observed in the cylindrical vessel with the d/h_L of 6.35. Moreover, the magnitudes of shear stress and energy dissipation rate in all the vessels tested were below their minimum values causing cells damage separately, which indicated that the hydromechanical‐stress environment in OSRs is suitable for cells cultivation in suspension. Finally, the CFD results suggested that the d/h_L higher than 8.80 should not be adopted for the 10 L cylindrical OSR at the shaking speed of 180 rpm because the “out of phase” state probably will happen there.     

Wind modelling for wind erosion research by open source computational fluid dynamics

The open source computational fluid dynamics (CFD) wind model (CFD-WEM) for wind erosion research in the Xilingele grassland in Inner Mongolia (autonomous region, China) is compared with two open source CFD models Gerris and OpenFOAM. The evaluation of these models was made according to software technology, implemented methods, handling, accuracy and calculation speed. All models were applied to the same wind tunnel data set. Results show that the simplest CFD-WEM has the highest calculation speed with acceptable accuracy, and the most powerful OpenFOAM produces the simulation with highest accuracy and the lowest calculation speed. Gerris is between CFD-WEM and OpenFOAM. It calculates faster than OpenFOAM, and it is capable to solve different CFD problems. CFD-WEM is the optimal model to be further developed for wind erosion research in Inner Mongolia grassland considering its efficiency and the uncertainties of other input data. However, for other applications using CFD technology, Gerris and OpenFOAM can be good choices. This paper shows the powerful capability of open source CFD software in wind erosion study, and advocates more involvement of open source technology in wind erosion and related ecological researches.     

Survivorship of juvenile barnacles and mussels: spatial dependence and the origin of alternative communities

Mussels, barnacles, and rockweeds often form a distinct mosaic of patches on rocky intertidal shores, and it has been suggested that these communities may represent alternative community states. One way that alternative community states can arise is if early successional events are scale-dependent, but it is not known if juvenile survivorships of mussels and barnacles are, in fact, scale-dependent. Scale-dependence of barnacles (Semibalanus balanoides (L.)) and mussels (Mytilus edulis L.) was tested in the Gulf of Maine, USA. In winter 1997, clearings of 1, 2, 4 and 8 m in diameter and uncleared controls were made in stands of the rockweed Ascophyllum nodosum (L.) Le Jolis at 12 sites spread evenly across four bays on Swan's Island, ME. Summer and fall-winter survivorship of barnacles, which recruited in spring 1997, were estimated by tracking the 1997 cohort until late winter 1998. Survivorship of mussels was estimated from following the fate of transplanted juveniles over 4 days in late August 1997. Both barnacles and mussels showed better survival in 4 and 8 m clearings than in small clearings and controls. There was also significant variation in survivorship among sites. Densities of gastropods in the clearings did not reflect survivorship patterns of barnacles and mussels. Barnacle survivorship increased in fall and winter, and in large clearings was comparable to survivorship in barnacle-dominated habitats. Mussel survivorship was low in all clearing sizes suggesting that mussel beds develop slowly.     

Computational fluid dynamics modeling of the paddle dissolution apparatus: Agitation rate,mixing patterns,and fluid velocities

The purpose of this research was to further investigate the hydrodynamics of the United States Pharmacopeia (USP) paddle dissolution apparatus using a previously generated computational fluid dynamics (CFD) model. The influence of paddle rotational speed on the hydrodynamics in the dissolution vessel was simulated. The maximum velocity magnitude for axial and tangential velocities at different locations in the vessel was found to increase linearly with the paddle rotational speed. Path-lines of fluid mixing, which were examined from a central region at the base of the vessel, did not reveal a region of poor mixing between the upper cylin-drical and lower hemispherical volumes, as previously speculated. Considerable differences in the resulting flow patterns were observed for paddle rotational speeds between 25 and 150 rpm. The approximate time required to achieve complete mixing varied between 2 to 5 seconds at 150 rpm and 40 to 60 seconds at 25 rpm, although complete mixing was achievable for each speed examined. An analysis of CFD-generated velocities above the top surface of a cylindrical compact positioned at the base of the vessel, below the center of the rotating paddle, revealed that the fluid in this region was undergoing solid body rotation. An examination of the velocity boundary layers adjacent to the curved surface of the compact revealed large peaks in the shear rates for a region within∼3 mm from the base of the compact, consistent with a ‘grooving’ effect, which had been previously seen on the surface of compacts following dissolution, associated with a higher dissolution rate in this region.     

An industrial perspective on bioreactor scale-down: what we can learn from combined large-scale bioprocess and model fluid studies

For industrial bioreactor design, operation, control and optimization, the scale-down approach is often advocated to efficiently generate data on a small scale, and effectively apply suggested improvements to the industrial scale. In all cases it is important to ensure that the scale-down conditions are representative of the real large-scale bioprocess. Progress is hampered by limited detailed and local information from large-scale bioprocesses. Complementary to real fermentation studies, physical aspects of model fluids such as air-water in large bioreactors provide useful information with limited effort and cost. Still, in industrial practice, investments of time, capital and resources often prohibit systematic work, although, in the end, savings obtained in this way are trivial compared to the expenses that result from real process disturbances, batch failures, and non-flyers with loss of business opportunity. Here we try to highlight what can be learned from real large-scale bioprocess in combination with model fluid studies, and to provide suitable computation tools to overcome data restrictions. Focus is on a specific well-documented case for a 30-m(3) bioreactor. Areas for further research from an industrial perspective are also indicated.     

The complex barnacle perfume: identification of waterborne pheromone homologues in Balanus improvisus and their differential expression during settlement

A key question in barnacle biology is the nature of cues that induce gregarious settlement. One of the characterised cues is the waterborne settlement pheromone (WSP). This study aimed to identify WSP homologues in Balanus improvisus and to investigate their expression during settlement. Six WSP homologues were identified, all containing an N-terminal signal peptide, a conserved core region, and a variable C-terminus comprising several -GR- and -HDDH- motifs. The B. improvisus WSP homologues were expressed in all settlement stages but showed different expression patterns. The homologue most similar to the B. amphitrite WSP was the most abundant and was constantly expressed during settlement. In contrast, several of the other WSP homologues showed the greatest expression in the juvenile stage. The presence of several WSP homologues suggests the existence of a pheromone mix, where con-specificity might be determined by a combination of sequence characteristics and the concentration of the individual components.