Computational fluid dynamics benchmark dataset of airflow in tracheas.

Computational Fluid Dynamics (CFD) is fast becoming a useful tool to aid clinicians in pre-surgical planning through the ability to provide information that could otherwise be extremely difficult if not impossible to obtain. However, in order to provide clinically relevant metrics, the accuracy of the computational method must be sufficiently high. There are many alternative methods employed in the process of performing CFD simulations within the airways, including different segmentation and meshing strategies, as well as alternative approaches to solving the Navier-Stokes equations. However, as in vivo validation of the simulated flow patterns within the airways is not possible, little exists in the way of validation of the various simulation techniques. The data presented here consists of very highly resolved flow data. The degree of resolution is compared to the highest necessary resolutions of the Kolmogorov length and time scales. Therefore this data is ideally suited to act as a benchmark case to which cheaper computational methods may be compared. A dataset and solution setup for one such more efficient method, large eddy simulation (LES), is also presented.

The effects of curvature and constriction on airflow and energy loss in pathological tracheas.

This paper considers factors that play a significant role in determining inspiratory pressure and energy losses in the human trachea. Previous characterisations of pathological geometry changes have focussed on relating airway constriction and subsequent pressure loss, however many pathologies that affect the trachea cause deviation, increased curvature, constriction or a combination of these. This study investigates the effects of these measures on tracheal flow mechanics, using the compressive goitre (a thyroid gland enlargement) as an example. Computational fluid dynamics simulations were performed in airways affected by goitres (with differing geometric consequences) and a normal geometry for comparison. Realistic airways, derived from medical images, were used because idealised geometries often oversimplify the complex anatomy of the larynx and its effects on the flow. Two mechanisms, distinct from stenosis, were found to strongly affect airflow energy dissipation in the pathological tracheas. The jet emanating from the glottis displayed different impingement and breakdown patterns in pathological geometries and increased loss was associated with curvature.

Power loss mechanisms in pathological tracheas.

The effort required to inhale a breath of air is a critically important measure in assessing airway function. Although the contribution of the trachea to the total flow resistance of the airways is generally modest, pathological alterations in tracheal geometry can have a significant negative effect. This study investigates the mechanisms of flow energy loss in a healthy trachea and in four geometries affected by retrosternal goitre which can cause significant distortions of tracheal geometry including constriction and deviation with abnormal curvature. By separating out the component of energy loss related to the wall shear (frictional loss), striking differences are found between the patterns of energy dissipation in the normal and pathological tracheas. Furthermore the ratio of frictional to total loss is dramatically reduced in the pathological geometries.

Theory of mind in the early course of schizophrenia: stability, symptom and neurocognitive correlates, and relationship with functioning.

BACKGROUND: Numerous studies have reported links between theory of mind (ToM) deficits, neurocognition and negative symptoms with functional outcome in chronic schizophrenia patients. Although the ToM deficit has been observed in first-episode patients, fewer studies have addressed ToM as a possible trait marker, neurocognitive and symptom correlations longitudinally, and associations with later functioning. METHOD: Recent-onset schizophrenia patients (n = 77) were assessed at baseline after reaching medication stabilization, and again at 6 months (n = 48). Healthy controls (n = 21) were screened, and demographically comparable with the patients. ToM was assessed with a Social Animations Task (SAT), in which the participants' descriptions of scenes depicting abstract visual stimuli 'interacting' in three conditions (ToM, goal directed and random) were rated for degree of intentionality attributed to the figures and for appropriateness. Neurocognition, symptoms and role functioning were also assessed. RESULTS: On the SAT, patients had lower scores than controls for both intentionality (p < 0.01) and appropriateness (p < 0.01) during the ToM condition, at baseline and 6 months. The ToM deficit was stable and present even in remitted patients. Analyses at baseline and 6 months indicated that for patients, ToM intentionality and appropriateness were significantly correlated with neurocognition, negative symptoms and role functioning. The relationship between ToM and role functioning was mediated by negative symptoms. CONCLUSIONS: The ToM deficit was found in recent-onset schizophrenia patients and appears to be moderately trait-like. ToM is also moderately correlated with neurocognition, negative and positive symptoms, and role functioning. ToM appears to influence negative symptoms which in turn makes an impact on role functioning.

Dynamics of airflow in a short inhalation.

During a rapid inhalation, such as a sniff, the flow in the airways accelerates and decays quickly. The consequences for flow development and convective transport of an inhaled gas were investigated in a subject geometry extending from the nose to the bronchi. The progress of flow transition and the advance of an inhaled non-absorbed gas were determined using highly resolved simulations of a sniff 0.5 s long, 1 l s⁻¹ peak flow, 364 ml inhaled volume. In the nose, the distribution of airflow evolved through three phases: (i) an initial transient of about 50 ms, roughly the filling time for a nasal volume, (ii) quasi-equilibrium over the majority of the inhalation, and (iii) a terminating phase. Flow transition commenced in the supraglottic region within 20 ms, resulting in large-amplitude fluctuations persisting throughout the inhalation; in the nose, fluctuations that arose nearer peak flow were of much reduced intensity and diminished in the flow decay phase. Measures of gas concentration showed non-uniform build-up and wash-out of the inhaled gas in the nose. At the carina, the form of the temporal concentration profile reflected both shear dispersion and airway filling defects owing to recirculation regions.

Cognitive functioning in first-episode schizophrenia: MATRICS Consensus Cognitive Battery (MCCB) Profile of Impairment.

BACKGROUND: Although many studies have assessed cognitive functioning in first-episode schizophrenia (FESz), the pattern and severity of impairment across cognitive domains remain unclear. Moreover, few studies have directly compared the pattern of cognitive performance between FESz and chronic schizophrenia (CSz). In this study we examined the cognitive impairment profile in FESz using a standardized neurocognitive battery (MATRICS Consensus Cognitive Battery; MCCB). METHODS: MCCB data were compared from 105 FESz patients, 176 CSz patients and 300 non-psychiatric (NP) participants. Mixed model analysis evaluated group differences in MCCB profiles and relative strengths and weaknesses in the MCCB profiles of patients. Clinical implications of MCCB performance were also examined; we compared the proportion of participants from each group who exhibited clinically-significant global cognitive impairment based on the MCCB Overall Composite score. RESULTS: FESz and CSz showed impaired performance across all MCCB domains relative to NP. With the exception of relative preservation of working memory and social cognition in FESz, the MCCB domain scores were similar in FESz and CSz. The distribution of impairment on the Overall Composite score did not significantly differ between FESz and CSz; compared to NP, both patient groups were overrepresented in moderate and severe impairment categories. CONCLUSION: The pattern, magnitude, and distribution of severity of impairment in FESz were similar to that observed in CSz. However, early in the illness, there may be relative sparing of working memory and social cognition.

Lamination and mixing in three fundamental flow sequences driven by electromagnetic body forces.

This article pursues the idea that the degree of striations, called lamination, could be engineered to complement stretching and to design new sequential mixers. It explores lamination and mixing in three new mixing sequences experimentally driven by electromagnetic body forces. To generate these three mixing sequences, Lorentz body forces are dynamically controlled to vary the flow geometry produced by a pair of local jets. The first two sequences are inspired from the "tendril and whorl" and "blinking vortex" flows. The third novel sequence is called the "cat's eyes flip." These three mixing sequences exponentially stretch and laminate material lines representing the interface between two domains to be mixed. Moreover, the mixing coefficient (defined as 1-σ(2)/σ(0)(2) where σ(2)/σ(0)(2) is the rescaled variance) and its rate grow exponentially before saturation. This saturation of the mixing process is related to the departure of the mixing rate from an exponential growth when the striations' thicknesses reach the diffusive length scale of the measurements or species and dyes. Incidentally, in our experiments, for the same energy or forcing input, the cat's eyes flip sequence has higher lamination, stretching, and mixing rates than the tendril and whorl and the blinking vortex sequences. These features show that bakerlike in situ mixers can be conceived by dynamically controlling a pair of local jets and by integrating lamination during stirring stages with persistent geometries. Combined with novel insights provided by the quantification of the lamination, this paper should offer perspectives for the development of new sequential mixers, possibly on all scales.

Decomposition and description of the nasal cavity form.

Patient-specific studies of physiological flows rely on anatomically realistic or idealized models. Objective comparison of datasets or the relation of specific to idealized geometries has largely been performed in an ad hoc manner. Here, two rational procedures (based respectively on Fourier descriptors and medial axis (MA) transforms) are presented; each provides a compact representation of a complex anatomical region, specifically the nasal airways. The techniques are extended to furnish average geometries. These retain a sensible anatomical form, facilitating the identification of a specific anatomy as a set of weighted perturbations about the average. Both representations enable a rapid translation of the surface description into a virtual model for computation of airflow, enabling future work to comprehensively investigate the relation between anatomic form and flow-associated function, for the airways or for other complex biological conduits. The methodology based on MA transforms is shown to allow flexible geometric modeling, as illustrated by a local alteration in airway patency. Computational simulations of steady inspiratory flow are used to explore the relation between the flow in individual vs. averaged anatomical geometries. Results show characteristic flow measures of the averaged geometries to be within the range obtained from the original three subjects, irrespective of averaging procedure. However the effective regularization of anatomic form resulting from the shape averaging was found to significantly reduce trans-nasal pressure loss and the mean shear stress in the cavity. It is suggested that this may have implications in attempts to relate model geometries and flow patterns that are broadly representative.

Inflow boundary profile prescription for numerical simulation of nasal airflow.

Knowledge of how air flows through the nasal passages relies heavily on model studies, as the complexity and relative inaccessibility of the anatomy prevents detailed in vivo measurement. Almost all models to date fail to incorporate the geometry of the external nose, instead employing a truncated inflow. Typically, flow is specified to enter the model domain either directly at the nares (nostrils), or via an artificial pipe inflow tract attached to the nares. This study investigates the effect of the inflow geometry on flow predictions during steady nasal inspiration. Models that fully replicate the internal and external nasal airways of two anatomically distinct subjects are used as a reference to compare the effects of common inflow treatments on physiologically relevant quantities including regional wall shear stress and particle residence time distributions. Inflow geometry truncation is found to affect flow predictions significantly, though slightly less so for the subject displaying more pronounced passage area contraction up to the internal nasal valve. For both subject geometries, a tapered pipe inflow provides a better approximation to the natural inflow than a blunt velocity profile applied to the nares. Computational modelling issues are also briefly outlined, by comparing quantities predicted using different surface tessellations, and by evaluation of domain-splitting techniques.

Computational modeling of flow and gas exchange in models of the human maxillary sinus.

The present study uses numerical modeling to increase the understanding of sinus gas exchange, which is thought to be a factor in sinus disease. Order-of-magnitude estimates and computational fluid dynamics simulations were used to investigate convective and diffusive transport between the nose and the sinus in a range of simplified geometries. The interaction between mucociliary transport and gas exchange was modeled and found to be negligible. Diffusion was the dominant transport mechanism for small ostia and large concentration differences between the sinus and the nose, whereas convection was important for larger ostia or smaller concentration differences. The presence of one or more accessory ostia can increase the sinus ventilation rate by several orders of magnitude, because it allows a net flow through the sinus. Estimates of nitric oxide (NO) transport through the ostium based on measured sinus and nasal NO concentrations suggest that the sinuses cannot supply all the NO in nasally exhaled air.

Mixing through stirring of steady flow in small amplitude helical tubes.

In this paper we numerically simulate flow in a helical tube for physiological conditions using a co-ordinate mapping of the Navier-Stokes equations. Helical geometries have been proposed for use as bypass grafts, arterial stents and as an idealized model for the out-of-plane curvature of arteries. Small amplitude helical tubes are also currently being investigated for possible application as A-V shunts, where preliminary in vivo tests suggest a possibly lower risk of thrombotic occlusion. In-plane mixing induced by the geometry is hypothesized to be an important mechanism. In this work, we focus mainly on a Reynolds number of 250 and investigate both the flow structure and the in-plane mixing in helical geometries with fixed pitch of 6 tube diameters (D), and centerline helical radius ranging from 0.1 D to 0.5 D. High-order particle tracking, and an information entropy measure is used to analyze the in-plane mixing. A combination of translational and rotational reference frames are shown to explain the apparent discrepancy between flow field and particle trajectories, whereby particle paths display a pattern characteristic of a double vortex, though the flow field reveals only a single dominant vortex. A radius of 0.25 D is found to provide the best trade-off between mixing and pressure loss, with little increase in mixing above R=0.25 D, whereas pressure continues to increase linearly.

Mechanics of airflow in the human nasal airways.

The mechanics of airflow in the human nasal airways is reviewed, drawing on the findings of experimental and computational model studies. Modelling inevitably requires simplifications and assumptions, particularly given the complexity of the nasal airways. The processes entailed in modelling the nasal airways (from defining the model, to its production and, finally, validating the results) is critically examined, both for physical models and for computational simulations. Uncertainty still surrounds the appropriateness of the various assumptions made in modelling, particularly with regard to the nature of flow. New results are presented in which high-speed particle image velocimetry (PIV) and direct numerical simulation are applied to investigate the development of flow instability in the nasal cavity. These illustrate some of the improved capabilities afforded by technological developments for future model studies. The need for further improvements in characterising airway geometry and flow together with promising new methods are briefly discussed.

Experimental investigation of nasal airflow.

The airway geometry of the nasal cavity is manifestly complex, and the manner in which it controls the airflow to accomplish its various physiological functions is not fully understood. Since the complex morphology and inaccessibility of the nasal passageways precludes detailed in-vivo measurements, either computational simulation or in-vitro experiments are needed to determine how anatomical form and function are related. The fabrication of a replica model of the nasal cavity, of a high optical clarity and derived from in-vivo scan data is described here, together with characteristics of the flow field investigated using particle image velocimetry (PIV) and flow visualization. Flow visualization is shown to be a capable and convenient technique for identifying key phenomena. Specifically the emergence of the jet from the internal nasal valve into the main cavity, how it impacts on the middle turbinate, and the large enhancement of dispersion that accompanies the initial appearance of flow instability are revealed as particularly significant features. The findings from the visualization experiments are complemented by PIV imaging, which provides quantitative detail on the variations in velocity in different regions of the nasal cavity. These results demonstrate the effectiveness of the cavity geometry in partitioning the flow into high shear zones, which facilitate rapid heat transfer and humidification from the nasal mucosa, and slower zones affording greater residence times to facilitate olfactory sensing. The experimental results not only provide a basis for comparison with other computational modelling but also demonstrate an alternative and flexible means to investigate complex flows, relevant to studies in different parts of the respiratory or cardiovascular systems.

Nasal architecture: form and flow.

Current approaches to model nasal airflow are reviewed in this study, and new findings presented. These new results make use of improvements to computational and experimental techniques and resources, which now allow key dynamical features to be investigated, and offer rational procedures to relate variations in anatomical form. Specifically, both replica and simplified airways of a single subject were investigated and compared with the replica airways of two other individuals with overtly differing geometries. Procedures to characterize and compare complex nasal airway geometry are first outlined. It is then shown that coupled computational and experimental studies, capable of obtaining highly resolved data, reveal internal flow structures in both intrinsically steady and unsteady situations. The results presented demonstrate that the intimate relation between nasal form and flow can be explored in greater detail than hitherto possible. By outlining means to compare complex airway geometries and demonstrating the effects of rational geometric simplification on the flow structure, this work offers a fresh approach to studies of how natural conduits guide and control flow. The concepts and tools address issues that are thus generic to flow studies in other physiological systems.

Numerical simulations of phase contrast velocity mapping of complex flows in an anatomically realistic bypass graft geometry.

Combined in vitro experiments and numerical simulations were performed to study flow artifacts in phase contrast (PC) velocity mapping of steady flow through an anatomically realistic aortocoronary bypass graft model. The geometry was obtained through imaging and computational reconstruction of a left anterior descending (LAD) coronary artery of a porcine heart. Simulated images of through-plane velocity were obtained at selected slices of the geometry. These were then compared and contrasted with velocity images of corresponding sites that were obtained from in vitro experiments. The shift and distortion of the measured velocity profile was well predicted by the simulation, while trajectories obtained from particle tracking were shown to be useful in understanding the origins of the flow artifacts that were observed.

Local and global geometric influence on steady flow in distal anastomoses of peripheral bypass grafts.

We consider the effect of geometrical configuration on the steady flow field of representative geometries from an in vivo anatomical data set of end-to-side distal anastomoses constructed as part of a peripheral bypass graft. Using a geometrical classification technique, we select the anastomoses of three representative patients according to the angle between the graft and proximal host vessels (GPA) and the planarity of the anastomotic configuration. The geometries considered include two surgically tunneled grafts with shallow GPAs which are relatively planar but have different lumen characteristics, one case exhibiting a local restriction at the perianastomotic graft and proximal host whilst the other case has a relatively uniform cross section. The third case is nonplanar and characterized by a wide GPA resulting from the graft being constructed superficially from an in situ vein. In all three models the same peripheral resistance was imposed at the computational outflows of the distal and proximal host vessels and this condition, combined with the effect of the anastomotic geometry, has been observed to reasonably reproduce the in vivo flow split. By analyzing the flow fields we demonstrate how the local and global geometric characteristics influences the distribution of wall shear stress and the steady transport of fluid particles. Specifically, in vessels that have a global geometric characteristic we observe that the wall shear stress depends on large scale geometrical factors, e.g., the curvature and planarity of blood vessels. In contrast, the wall shear stress distribution and local mixing is significantly influenced by morphology and location of restrictions, particular when there is a shallow GPA. A combination of local and global effects are also possible as demonstrated in our third study of an anastomosis with a larger GPA. These relatively simple observations highlight the need to distinguish between local and global geometric influences for a given reconstruction. We further present the geometrical evolution of the anastomoses over a series of follow-up studies and observe how the lumen progresses towards the faster bulk flow of the velocity in the original geometry. This mechanism is consistent with the luminal changes in recirculation regions that experience low wall shear stress. In the shallow GPA anastomoses the proximal part of the native host vessel occludes or stenoses earlier than in the case with wide GPA. A potential contribution to this behavior is suggested by the stronger mixing that characterizes anastomoses with large GPA.

Automated classification of peripheral distal by-pass geometries reconstructed from medical data.

Abnormal haemodynamic conditions are implicated in the development of anastomotic myointimal hyperplasia (MIH). However, these conditions are difficult to determine in vivo, prompting research using ex vivo idealised models. To relate the understanding gained in idealised geometries to anatomically correct conditions we have investigated a reproducible approach to classify in vivo distal graft anastomoses and their inter-patient variability. In vivo distal anastomotic geometries were acquired by magnetic resonance (MR) angiography from 13 patients who had undergone infrageniculate autologous venous by-pass surgery. On average, the images were acquired 2 weeks post-operatively. Five patients also underwent repeat examinations 2 to 7 weeks later. For each geometry, the surface of the arterial lumen is represented by the zero level set of an implicit function constructed from radial basis functions that minimise curvature. The three-dimensional binary image created from the interpolated surface is processed using a skeletonisation algorithm to obtain the centreline of each branch in the geometry. This allows for the measurement of the branching angles between straight line approximations of the centrelines of each vessel, averaging them over a characteristic length of each anastomosis. The main finding in the application of the proposed classification methodology to this set of patients is that the spectrum of anastomoses can be reduced to a small subset of cases characterised by two angles: the angle between the graft and the plane of the host artery and the angle between the graft and the proximal branch of the artery.

The influence of out-of-plane geometry on pulsatile flow within a distal end-to-side anastomosis.

We present an experimental and computational investigation of time-varying flow in an idealized fully occluded 45 degrees distal end-to-side anastomosis. Two geometric configurations are assessed, one where the centerlines of host and bypass vessels lie within a plane, and one where the bypass vessel is deformed out of the plane of symmetry, respectively, termed planar and non-planar. Flow experiments were conducted by magnetic resonance imaging in rigid wall models and computations were performed using a high order spectral/hp algorithm. Results indicate a significant change in the spatial distribution of wall shear stress and a reduction of the time-averaged peak wall shear stress magnitude by 10% in the non-planar model as compared to the planar configuration. In the planar geometry the stagnation point follows a straight-line path along the host artery bed with a path length of 0.8 diameters. By contrast in the non-planar case the stagnation point oscillates about a center that is located off the symmetry plane intersection with the host artery bed wall, and follows a parabolic path with a 0.7 diameter longitudinal and 0.5 diameter transverse excursion. A definition of the oscillatory shear index (OSI) is introduced that varies between 0 and 0.5 and that accounts for a continuous range of wall shear stress vector angles. In both models, regions of elevated oscillatory shear were spatially associated with regions of separated or oscillating stagnation point flow. The mean oscillatory shear magnitude (considering sites where OSI>0.1) in the non-planar geometry was reduced by 22% as compared to the planar configuration. These changes in the dynamic behavior of the stagnation point and the oscillatory shear distribution introduced by out-of-plane graft curvature may influence the localization of vessel wall sites exposed to physiologically unfavorable flow conditions.

Unsteady near wall residence times and shear exposure in model distal arterial bypass grafts.

Building on previous studies of unsteady flow within model distal bypass grafts we analyse the near wall residence times and shear exposure in a 45 degrees anastomosis under symmetrical and symmetry breaking geometric configurations. We define residence time as the minimum time for a particle to exit a spherical region and shear exposure as a temporal integral of the Huber-Henky-von-Mises criterion along a particle path over a fixed time interval. Decomposing the pulsatile cycle into four equal intervals we find that the interval of peak residence time in the host vessel is from mid-deceleration to peak diastole and peak diastole to mid-acceleration. The asymmetric model is shown to have a significantly lower residence time during these intervals. Considering the shear exposure prior to the residence time evaluation we determine that a higher average shear exposure exists in the asymmetric model associated with the upstream geometry modification. Analysis of the regions of high residence time and shear exposure suggests that the "toe" region and the interface between the "heel" and bulk flow are more significant than the bed and heel region. Although the asymmetric model considered in this study reduces residence times in the host artery, the product of the measure of shear exposure and residence time is not found to be preferable. If shear exposure were to be considered as an important factor in particle activation, the findings imply that for junction optimisation, greater consideration needs to be given both to the local junction asymmetry and upstream influence on the shear history.