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.

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.

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.

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.

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.

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.

Comparison of the antioxidant status in tracheal and bronchoalveolar epithelial lining fluids in recurrent airway obstruction.

REASONS FOR PERFORMING STUDY: Following a period of airway inflammation the clearance of inflammatory cells along the mucociliary escalator may impose a considerable oxidant load on the trachea. OBJECTIVES: To determine the degree of oxidative stress in tracheal epithelial lining fluid (ELF) in comparison to that present in peripheral airways after an acute exposure to organic dust. METHODS: Tracheal wash fluid and bronchoalveolar lavage fluid (BALF) were collected for cytology and antioxidant analyses from 6 recurrent airway obstruction (RAO)-affected horses and 6 healthy control horses before and after stabling on straw bedding for 24 h. RESULTS: In RAO-affected horses, organic dust exposure resulted in a significant decrease in ascorbic acid concentration in tracheal ELF (P<0.0001), which was greater than the decrease in bronchoalveolar ELF (P = 0.0003). The percentage decrease in tracheal ELF ascorbic acid correlated with the percentage decrease in bronchoalveolar ELF ascorbic acid (r = 0.76; P = 0.004) following exposure. CONCLUSIONS: Acute organic dust exposure results in significant antioxidant depletion in the trachea, which may reflect inflammation and oxidative processes in peripheral airways. POTENTIAL RELEVANCE: Further work is required to evaluate the role of ascorbic acid depletion in the pathogenesis of RAO.

A model of non-uniform lung parenchyma distortion.

A finite element model of mammalian lung parenchyma is used to study the effect of large non-uniform distortions on lung elastic behaviour. The non-uniform distortion is a uni-axial stretch from an initial state of uniform pressure expansion. For small distortions, the parenchymal properties are linearly isotropic and described by two elastic moduli. However, for large distortions, the parenchyma has anisotropic non-linear elastic properties described by five independent elastic moduli dependent on the degree of distortion; they are computed for a range of distortions and initial pressures. Ez, the Young's modulus in the direction of stretch, increases significantly with distortion (epsilon(z)) while Ex, the Young's modulus in the plane perpendicular to the stretch, is approximately constant. The greater the initial pressure, the bigger the difference between the two moduli at larger distortion strains. The shear modulus G(xz) is approximately independent of degree of distortion except at the highest initial pressure. The Poisson's ratio, nu(xz) is approximately constant with distortion strain for lower initial pressures, but increases significantly with epsilon(z) at higher pressures. Model predictions of the relation between G(xz) and initial uniform inflation pressure show a good correlation with reported experimental data for small distortion strains in a range of species. The model also exhibits similar behaviour to the experimentally measured uni-axial large deformations of a tri-axially pre-loaded block of parenchyma (Hoppin et al., 1975, Journal of Applied Physiology 39, 742-751).

Changes in selected physiological and laboratory measurements in elite horses competing in a 160 km endurance ride.

REASON FOR PERFORMING STUDY: Limited information exists about the physiological changes and clinical problems that occur in elite horses competing in high-speed 160 km endurance races. OBJECTIVES: To provide initial data describing changes in physiological and laboratory measurements in horses competing in a high-speed, 160 km endurance race under temperate conditions and to compare data between horses that successfully completed the race and those that failed to finish. METHODS: Body mass (BM) was measured, blood samples were collected, and veterinary examinations performed on horses before, during, and at the finish of a CEI*** 160 km endurance race. RESULTS: Of 36 horses participating in the study, 22 (61%) completed the race. Twelve horses were eliminated for lameness and 2 for persistent heart rate elevation. Mean speed of finishers was 15.2 km/h. Mean +/- s.d. BM loss of finishers at the end of the race (5.7 +/- 2.6%) was not different (P = 0.58) from BM loss of nonfinishers at elimination (6.7 +/- 34%). Similarly, there were no significant differences in heart rate or veterinary assessment of hydration at the race end for finishers as compared to the elimination point for nonfinishers. PCV increased while sodium, chloride and potassium concentrations decreased with exercise but differences between finishers and nonfinishers were not detected. In contrast, both total and ionised calcium concentrations decreased in successful horses but remained unchanged in nonfinishers. CONCLUSIONS: Elite endurance horses are more likely to be eliminated from competition for lameness than metabolic problems; however, it remains unclear whether these conditions are entirely distinct. The magnitude of the decrease in sodium concentration in both finishers and nonfinishers was greater than in previous reports of 160 km rides. POTENTIAL RELEVANCE: These data should be of use for both organisers and participants in elite 160 km endurance races. The tendency toward hyponatraemia as well as the difference in calcium concentrations between finishers and nonfinishers warrant further study.

The effect of inspired gas density on pulmonary artery transmural pressure and exercise induced pulmonary haemorrhage.

REASONS FOR PERFORMING STUDY: Pulmonary capillary stress failure, largely as a result of high pulmonary vascular pressures, has been implicated in the aetiology of EIPH. However, the role of the respiratory system in determining the magnitude of EIPH has received little attention. HYPOTHESIS: Horses breathing a gas of greater density than air will exhibit greater transmural pulmonary arterial pressures (TPAP) and more severe EIPH, and horses breathing a gas of lower density than air will exhibit lower TPAP and less severe EIPH, both compared with horses breathing air. METHODS: Following a warm-up, 8 Thoroughbred horses were exercised for 1 min at 10, 11 and 12 m/sec (5 degrees incline) breathing air or 21% oxygen/79% helium or 21% oxygen/79% argon in a randomised order. Heart rate, respiratory rate, pulmonary arterial pressure and oesophageal pressure were measured during exercise. Bronchoalveolar lavage fluid (BALF) was collected from the dorsocaudal regions of the left and right lungs 40 min post exercise and red blood cell (RBC) counts were performed. RESULTS: The exercise tests induced mild EIPH. Maximum changes in oesophageal pressure were lower on helium-oxygen compared to argon-oxygen (P<0.001). TPAP and median RBC counts did not differ between gas mixtures. BALF RBC counts from the left lung correlated with counts from the right lung (P<0.0001). However BALF RBC counts from the left lung were higher than those from the right lung (P = 0.004). CONCLUSION: As alterations in pulmonary arterial and oesophageal pressure caused by changes in inspired gas density were of similar magnitude, TPAP remained unchanged and there was no significant effect on EIPH severity. POTENTIAL RELEVANCE: Manipulations that decrease swings in intrapleural pressure may only decrease the degree of EIPH in horses severely affected by the condition.

Antioxidant and inflammatory responses of healthy horses and horses affected by recurrent airway obstruction to inhaled ozone.

REASONS FOR PERFORMING STUDY: Inhaled ozone can induce oxidative injury and airway inflammation. Horses affected by recurrent airway obstruction (RAO) have a decreased pulmonary antioxidant capacity, which may render them more susceptible to oxidative challenge. It is currently unknown whether RAO-affected horses are more susceptible to oxidative stress than those unaffected by RAO. OBJECTIVES: To determine whether ozone exposure induces greater oxidative stress and airway inflammation in RAO-affected horses in remission than in healthy horses. METHODS: Seven healthy control horses and 7 RAO-affected horses were exposed to 0.8 ppm ozone for 2 h at rest. RESULTS: At baseline, bronchoalveolar lavage fluid (BALF) ascorbic acid concentrations were lower in RAO-affected horses than healthy controls. Ozone appeared to preferentially oxidise glutathione rather than ascorbic acid 6 h after exposure. Individual healthy and RAO-affected horses demonstrated oxidation of BALF glutathione after ozone exposure. Overall, RAO-affected horses did not demonstrate increased oxidative stress following ozone exposure, compared with healthy horses. Ozone did not induce significant airway inflammation in either group. CONCLUSIONS: RAO-affected horses in remission are not more sensitive to ozone despite a decreased pulmonary antioxidant capacity. Sensitivity to ozone appears to be independent of initial pulmonary antioxidant status. POTENTIAL RELEVANCE: Horses with high susceptibility to oxidative stress may benefit from antioxidant supplementation.

Recording of ECG signals on a portable MiniDisc recorder for time and frequency domain heart rate variability analysis.

Analysis of heart rate variability (HRV) is a non-invasive technique useful for investigating autonomic function in both humans and animals. It has been used for research into both behaviour and physiology. Commercial systems for human HRV analysis are expensive and may not have sufficient flexibility for appropriate analysis in animals. Some heart rate monitors have the facility to provide inter-beat interval (IBI), but verification following collection is not possible as only IBIs are recorded, and not the raw electrocardiogram (ECG) signal. Computer-based data acquisition and analysis systems such as Po-Ne-Mah and Biopac offer greater flexibility and control but have limited portability. Many laboratories and veterinary surgeons have access to ECG machines but do not have equipment to record ECG signals for further analysis. The aim of the present study was to determine whether suitable HRV data could be obtained from ECG signals recorded onto a MiniDisc (MD) and subsequently digitised and analysed using a commercial data acquisition and analysis package. ECG signals were obtained from six Thoroughbred horses by telemetry. A split BNC connecter was used to allow simultaneous digitisation of analogue output from the ECG receiver unit by a computerised data acquisition system (Po-Ne-Mah) and MiniDisc player (MZ-N710, Sony). Following recording, data were played back from the MiniDisc into the same input channel of the data acquisition system as previously used to record the direct ECG. All data were digitised at a sampling rate of 500 Hz. IBI data were analysed in both time and frequency domains and comparisons between direct recorded and MiniDisc data were made using Bland-Altman analysis. Despite some changes in ECG morphology due to loss of low frequency content (primarily below 5 Hz) following MiniDisc recording, there was minimal difference in IBI or time or frequency domain analysis between the two recording methods. The MiniDisc offers a cost-effective approach to intermediate recording of ECG signals for subsequent HRV analysis and also provides greater flexibility than use of human Holter systems.

Antioxidant supplementation and pulmonary function at rest and exercise.

Antioxidants have been implicated in the reduction and prevention of oxidative stress during exercise. We hypothesised that a dietary supplement containing a mixture of natural antioxidants together with vitamins E, C and selenium, given for 4 weeks, would increase the systemic and pulmonary antioxidant capacity leading to a reduction in markers of oxidative damage and an improvement in pulmonary function during exercise. In 6 healthy horses studied, the antioxidant supplement significantly increased plasma concentrations of ascorbic acid (from mean +/- s.d. 16 +/- 7 to 23 +/- 4 micromol/l; P = 0.007) and alpha-tocopherol (from 10 +/- 3 to 14 +/- 3 micromol/l; P = 0.02) and increased the bronchoalveolar lavage pulmonary epithelial lining fluid (ELF) concentration of ascorbic acid compared to a placebo, but not significantly (2.0 +/- 0.9 mmol/l and 1.2 +/- 0.9 mmol/l, respectively; P>0.05). Alpha-tocopherol was not detected in ELF either before or after supplementation or exercise. The mean concentration of malondialdehyde (MDA) in ELF was lower following antioxidant supplementation compared to placebo and control periods, but not significantly. An intermittent exercise test consisting of 2 min at 70, 80 and 90% of the horses' individual maximum oxygen uptake, failed to induce significant systemic or pulmonary oxidative stress (based on the glutathione redox ratio (GRR) and the ascorbic acid redox ratio (ARR)) and lipid peroxidation (based on the concentration of thiobarbituric acid reactive substances in plasma and MDA in ELF) either for placebo or antioxidant treatments. There was a strong correlation between GRR and ARR in the pulmonary epithelial lining fluid (r = 0.89; P<0.0001). In healthy horses on a diet containing adequate levels of antioxidants, additional antioxidant supplementation has no apparent beneficial or detrimental effect on pulmonary function during moderate intensity exercise. The importance of antioxidant supplementation may only become apparent if the diet is deficient in antioxidants, if exercise intensity is higher or more prolonged, or if disease or additional stresses are present.

Modelling the oxygen cost of transport in competitions over ground of variable slope.

This study provides an objective method for estimating the oxygen consumption of horses while running on variable slopes so that realistic comparisons may be made of the locomotory transport cost involved in 3-day events, particularly the Speed and Endurance Test, at sites of differing terrain. A knowledge of the work profile over a particular course would enable competitors to plan speed and interval times appropriately along its length. We have developed a semi-empirical, but mechanistically based, model to calculate the oxygen cost of transport [COTpath in ml O2/kg/m path] for running on the flat, up or down a slope of given gradient (from -0.3 to +0.3). The model is then used to calculate the overall effort of running on a number of 3-day event courses of differing standard; the model does not assess the energetic cost of jumping. The cost of transport over the range of gradient of -0.3 to +0.3 was modelled using the following equations: On the flat or uphill: COTpath = 0.123 + 1.561(gradient); Downhill: COTpath = 0.123 + 1.591(gradient) + 9.762(gradient)2 + 14.0(gradient)3.

Quantification of the response of equine apocrine sweat glands to beta2-adrenergic stimulation.

The aim of the present study was to characterise the quantitative sweating response of the horse to beta2-adrenergic stimulation. The sweating responses of 6 horses to the randomised infusion of 8 different adrenaline concentrations (0.025, 0.05, 0.075, 0.1, 0.2, 0.4, 1.0 or 2.0 microg/kg bwt/min), was investigated. Sweating rate (SR) and skin temperature (TSK) on the neck (N) and gluteal region (G), and plasma adrenaline and noradrenaline concentrations were measured. Peak SR was approximately 15 (N) and approximately 9 g/m2/min (G) during infusion of both 1.0 and 2.0 microg/kg bwt/min adrenaline. Sweat produced per nmol/l plasma adrenaline peaked during the infusion of 0.075 microg/kg bwt/min adrenaline. Higher adrenaline infusion concentrations resulted in a progressive decrease in the amount of sweat produced per nmol/l plasma adrenaline and a plateau of 6 g/m2/(nmol/l) plasma adrenaline was reached for infusions between 1.0 and 2.0 microg/kg bwt/min. Peak SR were far lower than we have previously reported during exercise. There was no evidence of sweat gland fatigue or vasoconstriction during infusion, suggesting saturation of sweat gland beta2 receptors. We conclude that sweating in the horse is under dual control from a combination of hormonal and neural mechanisms.

Recovery from transport and acclimatisation of competition horses in a hot humid environment.

The aims of the present field-based study were to investigate changes in fit horses undergoing acclimatisation to a hot humid environment and to provide data on which to base recommendations for safe transport and acclimatisation. Six horses (age 7-12 years) were flown from Europe to Atlanta and underwent a 16 day period of acclimatisation. Exercise conditions during acclimatisation (wet bulb globe temperature index 27.6+/-0.0 [mean +/- s.e.]) were more thermally stressful compared with the European climate from which the horses had come (22.0+/-1.8, P<0.001). Following the flight, weight loss was 4.1+/-0.8% bodyweight and took around 7 days to recover. Water intake during the day was significantly increased (P<0.05) compared with night during acclimatisation. Daily mean exercise duration was 72+/-12 min and the majority of work was performed with a heart rate below 120 beats/min. Respiratory rate (fR) was increased (P<0.05) throughout acclimatisation compared with in Europe, but resting morning (AM) and evening (PM) rectal temperature (TREC), heart rate (fC) and plasma volume were unchanged. White blood cell (WBC) count was significantly increased at AM compared with in Europe on Days 4 and 10 of acclimatisation (P<0.01), but was not different by Day 16. In conclusion, horses exposed to hot humid environmental conditions without prior acclimatisation are able to accommodate these stresses and, with appropriate management, remain fit and clinically healthy, without significant risk of heat illness or heat-related disorders, provided they are allowed sufficient time to recover from transport, acclimatisation is undertaken gradually and they are monitored appropriately.

Analysis of mechanical stresses within the alveolar septa leading to pulmonary edema.

Mechanical ventilation has been associated with pulmonary edema in the clinical setting, but the pathophysiological mechanisms of this process have not been clearly defined. Experimental studies have shown that high transpulmonary pressures resulting from ventilation may damage the capillary walls, thereby leading to edema. Knowledge of the stress distribution within the alveolar septa would be an important step in understanding this phenomenon. A newly developed saline-filled alveolar sac model was utilized for analysis of septal stresses in young and aging healthy lungs, in order to examine their vulnerability to pulmonary edema during ventilation. Significant stress concentrations were shown to develop near highly curved regions (small local radii of less than 4 mum in a lung inflated to 80% could be as high as 25 times that of average septal stresses. The combination of elevated stress sites that are formed in the stiffer parenchyma of the aging lung, together with the cyclic loading of ventilation, may explain the gaps and breaks previously observed in pulmonary edema.