The effect of airway motion and breathing phase during imaging on CFD simulations of respiratory airflow

RATIONALE

Computational fluid dynamics (CFD) simulations of respiratory airflow can quantify clinically useful information that cannot be obtained directly, such as the work of breathing (WOB), resistance to airflow, and pressure loss. However, patient-specific CFD simulations are often based on medical imaging that does not capture airway motion and thus may not represent true physiology, directly affecting those measurements.

OBJECTIVES

To quantify the variation of respiratory airflow metrics obtained from static models of airway anatomy at several respiratory phases, temporally averaged airway anatomies, and dynamic models that incorporate physiological motion.

METHODS

Neonatal airway images were acquired during free-breathing using 3D high-resolution MRI and reconstructed at several respiratory phases in two healthy subjects and two with airway disease (tracheomalacia). For each subject, five static (end expiration, peak inspiration, end inspiration, peak expiration, averaged) and one dynamic CFD simulations were performed. WOB, airway resistance, and pressure loss across the trachea were obtained for each static simulation and compared with the dynamic simulation results.

RESULTS

Large differences were found in the airflow variables between the static simulations at various respiratory phases and the dynamic simulation. Depending on the static airway model used, WOB, resistance, and pressure loss varied up to 237%, 200%, and 94% compared to the dynamic simulation respectively.

CONCLUSIONS

Changes in tracheal size and shape throughout the breathing cycle directly affect respiratory airflow dynamics and breathing effort. Simulations incorporating realistic airway wall dynamics most closely represent airway physiology; if limited to static simulations, the airway geometry must be obtained during the respiratory phase of interest for a given pathology.

Assessing Changes in Airflow and Energy Loss in a Progressive Tracheal Compression Before and After Surgical Correction

The energy needed to drive airflow through the trachea normally constitutes a minor component of the work of breathing. However, with progressive tracheal compression, patient subjective symptoms can include severe breathing difficulties. Many patients suffer multiple respiratory co-morbidities and so it is important to assess compression effects when evaluating the need for surgery. This work describes the use of computational prediction to determine airflow resistance in compressed tracheal geometries reconstructed from a series of CT scans. Using energy flux analysis, the regions that contribute the most to airway resistance during inhalation are identified. The principal such region is where flow emerging from the zone of maximum constriction undergoes breakup and turbulent mixing. Secondary regions are also found below the tongue base and around the glottis, with overall airway resistance scaling nearly quadratically with flow rate. Since the anatomical extent of the imaged airway varied between scans-as commonly occurs with clinical data and when assessing reported differences between research studies-the effect of sub-glottic inflow truncation is considered. Analysis shows truncation alters the location of jet breakup and weakly influences the pattern of pressure recovery. Tests also show that placing a simple artificial glottis in the inflow to a truncated model can replicate patterns of energy loss in more extensive models, suggesting a means to assess sensitivity to domain truncation in tracheal airflow simulations.