Biomechanical Models of Human Upper and Tracheal Airway Functionality

The respiratory tract, in other words, the airway, is the primary airflow path for several physiological activities such as coughing, breathing, and sneezing. Diseases can impact airway functionality through various means including cancer of the head and neck, Neurological disorders such as Parkinson\u27s disease, and sleep disorders and all of which are considered in this study. In this dissertation, numerical modeling techniques were used to simulate three distinct airway diseases: a weak cough leading to aspiration, upper airway patency in obstructive sleep apnea, and tongue cancer in swallow disorders. The work described in this dissertation, therefore, divided into three biomechanical models, of which fluid and particulate dynamics model of cough is the first. Cough is an airway protective mechanism, which results from a coordinated series of respiratory, laryngeal, and pharyngeal muscle activity. Patients with diminished upper airway protection often exhibit cough impairment resulting in aspiration pneumonia. Computational Fluid Dynamics (CFD) technique was used to simulate airflow and penetrant behavior in the airway geometry reconstructed from Computed Tomography (CT) images acquired from participants. The second study describes Obstructive Sleep Apnea (OSA) and the effects of dilator muscular activation on the human retro-lingual airway in OSA. Computations were performed for the inspiration stage of the breathing cycle, utilizing a fluid-structure interaction (FSI) method to couple structural deformation with airflow dynamics. The spatiotemporal deformation of the structures surrounding the airway wall was predicted and found to be in general agreement with observed changes in luminal opening and the distribution of airflow from upright to supine posture. The third study describes the effects of cancer of the tongue base on tongue motion during swallow. A three-dimensional biomechanical model was developed and used to calculate the spatiotemporal deformation of the tongue under a sequence of movements which simulate the oral stage of swallow

Recent hominim cranial form and function

This thesis aims to assess if biting mechanics drives craniofacial morphology in recent hominins. To that end, a virtual functional morphology toolkit, that includes Finite Element Analysis (FEA) and Geometric Morphometrics (GM), is used to simulate biting, measure bite force and quantify deformations arising due to simulated biting in Homo sapiens and its proposed ancestral species, Homo heidelbergensis. Moreover, the mechanical significance of the frontal sinus and of the brow-ridge is also assessed in Kabwe 1 (a Homo heidelbergensis specimen). The frontal sinus is examined by comparing the mechanical performance in three FE models with varying sinus morphology. A similar approach is applied to the brow-ridge study. This approach relies on the assumption that FEA approximates reality. Thus, a validation study compares the deformations experienced by a real cranium under experimental loading with those experienced by an FE model under equivalent virtual loading to verify this assumption. A sensitivity analysis examines how simplifications in segmentation impact on FEA results. Lastly, the virtual reconstruction of Kabwe 1 is described.Results show that prediction of absolute strain magnitudes is not precise, but the distribution of regions of larger and smaller (i.e. pattern of) deformations experienced by the real cranium is reasonably approximated by FEA, despite discrepancies in the alveolus. Simplification of segmentation stiffens the model but has no impact on the pattern of deformations, with the exception of the alveolus. Comparison of the biting performance of Kabwe 1 and H. sapiens suggests that morphological differences between the two species are likely not driven by selection of the masticatory system. Frontal sinus morphogenesis and morphology are possibly impacted by biting mechanics in the sense that very low strains are experienced by this region. Because bone adapts to strains, the frontal sinus is possibly impacted by this mechanism. Lastly, biting mechanics has limited impact on brow-ridge morphology and does not explain fully the enormous brow-ridge of Kabwe 1. Hence, other explanations are necessary to explain this prominent feature

Mechanics of airflow in human inhalation

The mechanics of airflow in the large airways during inspiration affects important physiological functions such as ventilation, olfaction, heat exchange and mass transfer. The behaviour of the airflow is important not only for healthcare applications including diagnosis, intervention planning and assessment, but for inhalation toxicology. This research aims to further the understanding of human nasal physiology through computational modelling. Specifically, the effects of transient inhalation conditions on flow dynamics and transport were characterised and the changes in flow behaviour in response to certain pathologies quantified. The key findings can be summarised as follows: Firstly, the time scales for airflow in the large airways have been identified and the initial flow patterns revealed. Three phases in the temporal behaviour of the flow were identified (flow initiation, quasi-equilibrium and decay). The duration of each phase differs depending on the quantity of interest. Flow in the nose was characterised as transitional, whilst in parts of the descending airways it is turbulent, particularly in the faster moving regions around the jets which may occur in the pharynx, larynx and at the superior end of the trachea. The bulk of the flow is biased to fill only certain regions of the airways, whilst other regions carry little flow, due to features upstream. Analysis of cross-sectional images provided by medical imaging does not necessarily provide a representative view of the area available to the flow. Various scalar species were employed to represent the fate of nanoparticles and gaseous species within the airways. Only species with high diffusion rates exhibited significant absorption at the airway walls. Airway pathologies often cause changes to the geometry of the airway. One such pathology, the goitre, was found to curve the trachea and in some cases cause constriction. Both these geometric changes were found to increase the pressure loss and energy required to drive flow through the trachea. Furthermore, the flow in pathological cases was more disturbed. High resolution simulations have been used to address these topics and the scales simulated have been analysed in terms of the smallest features possible in the flow to determine their fidelity.Open Acces

Unveiling advanced mechanisms of inhalable drug aerosol dynamics using computational fluid dynamics and discrete element method

Capsule-based dry powder inhalers (DPIs) are widely used to treat chronic obstructive pulmonary disease (COPD) by delivering active pharmaceutical ingredients (APIs) via inhalation into human respiratory systems. Previous research has shown that the actuation flow rate, aerodynamic particle size distribution (APSD), and particle shape of lactose carriers are factors that can influence the particle deposition patterns in human respiratory systems. Understanding the dynamics of APIs transport in DPIs and airways can provide significant value for the design optimization of DPIs and particle shapes to enhance the delivery of APIs to the designated lung sites, i.e., small airways. Thus, it is necessary to investigate how to modulate the above-mentioned factors to increase the delivery efficacy to small airways and enhance the therapeutic effect to treat COPD. Compared with in vitro and in vivo methods, computational fluid-particle dynamics (CFPD) models allow researchers to study the transport dynamics of inhalable therapeutic dry powders in both DPI and human respiratory systems. However, existing CFPD models neglect particle-particle interactions, and most existing airway models lack peripheral lung airway and neglect the airway deformation kinematics. Such deficiencies can lead to inaccurate predictions of particle transport and deposition. This study developed a one-way coupled computational fluid dynamics (CFD) and discrete element method (DEM) model to simulate the particle-particle and particle-device interactions, and the transport of API-carrier dry powder mixtures with different shapes of carriers in a DPI flow channel. The influence of actuation flow rate (30 to 90 L/min) and particle shape (aspect ratio equals 1, 5, and 10) on lactose carrier dynamics in a representative DPI, i.e., SpirivaTM HandihalerTM, has been investigated. Subsequently, an elastic truncated whole-lung model has also been developed to predict particle transport and deposition from mouth to alveoli, with disease-specific airway deformation kinematics. Numerical results indicate that 90 L/min actuation flow rate generates the highest delivery efficiency of Handihaler, as approximately 26% API reaches the deep lung region. The elastic truncated whole-lung modeling results show that noticeable differences of predictions between static and elastic lung models can be found, which demonstrates the necessity to model airway deformation kinematics in virtual lung models

Numerical modelling of fluid-structure interactions for fluid-induced instability in the upper airway

This study is concerned with fluid-structure interactions (FSI) involved in the human upper airway, in particular, those associated with snoring and obstructive sleep apnoea/hypopnoea syndrome (OSAHS). Further examining this area of interest, the goal of the current research is to contribute further understanding and enhance development of computational modelling, for retroglossal obstruction and palatal snoring. To that end, the investigation was divided into three major parts. Firstly, extending previous laminar, 2-D reduced Navier-Stokes model, an idealised 3-D computational model was constructed for studying retroglossal obstruction. A full Navier-Stokes solver in an Arbitrary Lagrangian-Eulerian (ALE) framework was coupled to a linear thin shell, where both laminar and turbulent flow was investigated. Numerical results showed increase flow-induced tongue replica deflection under turbulent conditions and demonstrated cross-flow pressures that may encourage side wall collapse. In the second part of the thesis, palatal snoring was further examined and its potential to detect retroglossal obstruction was proposed. In order to investigate this, flow-induced instability of a cantilever plate in an obstructed channel was modeled and a relationship between critical velocity and obstruction depth was established. Correlating the critical velocity with typical breathing flow curve, a time difference between palatal snoring episodes or onset of palatal snoring, may represent a key variable for non-invasive measurement of retroglossal obstruction severity. Further enhancement of the 2-D computational model by including contact was proposed using an immersed boundary method (IBM). This may represent a more complete model of palatal snoring by modelling pre- and post-contact response of unstable cantilever plate, which showed potential to capture more complicated palatal snoring signals. Finally, the third part of this thesis examined flow-induced instability of a soft palate in a 3-D realistic upper airway. To model this, a full 3-D Navier-Stokes solver under an ALE framework was coupled to a non-linear soft palate model. Appropriate soft palate properties were applied giving palatal snoring frequency within range of clinically measured values. Palatal flutter was observed at high flow rates, demonstrating irreversible transfer of flow energy to soft palate. This computational model may perhaps be exploited for future investigation of more accurate palatal snoring, necessary for developing non-invasive snoring signals for measurement or diagnosis of retroglossal obstruction