An integrated 1D breathing lung simulation with relative hysteresis of airway structure and regional pressure for healthy and asthmatic human lungs

This study aims to develop a one-dimensional (1D) computational fluid dynamics (CFD) model with dynamic airway geometry that considers airway wall compliance and acinar dynamics. The proposed 1D model evaluates the pressure distribution and the hysteresis between the pressure and tidal volume (Vtidal) in the central and terminal airways for healthy and asthmatic subjects. Four-dimensional CT images were captured at 11-14 time points during the breathing cycle. The airway diameter and length were reconstructed using a volume-filling method and a stochastic model at respective time points. The obtained values for the airway diameter and length were interpolated via the Akima spline to avoid unboundedness. A 1D energy balance equation considering the effects of wall compliance and parenchymal inertance was solved using the efficient aggregationbased algebraic multigrid solver, a sparse matrix solver, reducing the computational costs by around 90% when compared with the generalized minimal residual solver. In the Vtidal versus displacement in the basal direction (z-coordinate), the inspiration curve was lower than the expiration curve, leading to relative hysteresis. The dynamic deformation model was the major factor influencing the difference in the workload in the central and terminal airways. In contrast, wall compliance and parenchymal inertance appeared only marginally to affect the pressure and workload. The integrated 1D model mimicked dynamic deformation by predicting airway diameter and length at each time point, describing the effects of wall compliance and parenchymal inertance. This computationally efficient model could be utilized to assess breathing mechanism as an alternative to pulmonary function tests.