The application of computational fluid dynamics (CFD) in simulating pulmonary breathing plays a vital role in advancing the understanding of lung physiology and pathology. Beyond providing mechanistic insights, CFD-based models enable the exploration ...
The application of computational fluid dynamics (CFD) in simulating pulmonary breathing plays a vital role in advancing the understanding of lung physiology and pathology. Beyond providing mechanistic insights, CFD-based models enable the exploration of multiple physiological and pathological scenarios at relatively low cost, thereby reducing the dependence on physical prototypes when boundary conditions change. Despite these advantages, the complex branching structure of the airways introduces challenges in mesh generation, often requiring expert intervention to ensure a high-quality mesh suitable for CFD simulations. Furthermore, dense meshes, while necessary to capture the intricate airway geometry, significantly increase computational demands, thus limiting the feasibility of real- time, patient-specific simulations .
To overcome these challenges, this thesis proposes a Deep Operator Network based on Graph Neural Network, namely GraphDeepONet, to develop a surrogate model for the 1D airway breathing simulations. The surrogate model is trained on a benchmark dataset comprising flow rate, pressure, and pleural pressure. Given the hierarchical airway geometry that induces multiscale flow variations, we propose a novel volume-based normalization technique to mitigate these effects during training.
Additionally, three distinct learning strategies are introduced based on GraphDeepONet. Model 1 employs a dual-network GraphDeepONet to predict flow rate and pressure at the same time with physical constraint. Model 2 first predicts flow rate and computes pressure based on the predicted flow rate in the training phase. Similar to Model 2 in computing pressure, but in flow rate prediction, Model 3 first predicts the flow rate at the acinar region, then summing up to obtain the other flow rate. Among them, Model 3 achieved a prediction L2 relative error of 10.44% (5.08%) in flow rate, 9.99% (5.52%) in pressure, and 9.29% (4.10%) in pleural pressure.
Overall, the proposed method significantly speeds up computational time at least 700x compared to conventional solvers. This work presents an efficient surrogate model for 1D airway breathing simulation, offering a potential application for real-time respiratory analysis and the development of personalized treatment strategies.