An enhanced respiratory mechanics model based on double-exponential and fractional calculus

We address mathematical modelling of respiratory mechanics and put forward a model based on double-exponential and fractional calculus for parameter estimation, model simulation, and evaluation based on actual data. Our model has been implemented on a publicly available executable code with adjustable parameters, making it suitable for different applications. Our analysis represents the first application of fractional calculus and double-exponential modelling to respiratory mechanics, and allows us to propose a hybrid model fitting experimental data in different ventilation modes. Furthermore, our model can be used to study the mechanical features of the respiratory system, improve the safety of ventilation techniques, reduce ventilation damages, and provide strong support for fast and adaptive determination of ventilation parameters

Biological Fluid Mechanics Under the Microscope: A Tribute to John Blake

John Blake (1947--2016) was a leader in fluid mechanics, his two principal areas of expertise being biological fluid mechanics on microscopic scales and bubble dynamics. He produced leading research and mentored others in both Australia, his home country, and the UK, his adopted home. This article reviews John Blake's contributions in biological fluid mechanics, as well as giving the author's personal viewpoint as one of the many graduate students and researchers who benefitted from his supervision, guidance and inspiration. The key topics from biological mechanics discussed are: `squirmer' models of protozoa, the method of images in Stokes flow and the `blakelet' solution, discrete cilia modelling via slender body theory, physiological flows in respiration and reproduction, blinking stokeslets in microorganism feeding, human sperm motility, and embryonic nodal cilia.Comment: 23 pages, 11 figures. Submitted versio

Modelling mucociliary clearance

Mathematical modelling of the fluid mechanics of mucociliary clearance (MCC) is reviewed and future challenges for researchers are discussed. The morphology of the bronchial and tracheal airway surface liquid (ASL) and ciliated epithelium are briefly introduced. The cilia beat cycle, beat frequency and metachronal coordination are described, along with the rheology of the mucous layer. Theoretical modelling of MCC from the late 1960s onwards is reviewed, and distinctions between ‘phenomenological’, ‘slender body theory’ and recent ‘fluid–structure interaction’ models are explained.\ud \ud The ASL consists of two layers, an overlying mucous layer and underlying watery periciliary layer (PCL) which bathes the cilia. Previous models have predicted very little transport of fluid in the PCL compared with the mucous layer. Fluorescent tracer transport experiments on human airway cultures conducted by Matsui et al. [Matsui, H., Randell, S.H., Peretti, S.W., Davis, C.W., Boucher, R.C., 1998. Coordinated clearance of periciliary liquid and mucus from airway surfaces. J. Clin. Invest. 102 (6), 1125–1131] apparently showed equal transport in both the PCL and mucous layer. Recent attempts to resolve this discrepancy by the present authors are reviewed, along with associated modelling findings. These findings have suggested new insights into the interaction of cilia with mucus due to pressure gradients associated with the flat PCL/mucus interface. This phenomenon complements previously known mechanisms for ciliary propulsion. Modelling results are related to clinical findings, in particular the increased MCC observed in patients with pseudohypoaldosteronism. Recent important advances by several groups in modelling the fluid–structure interaction by which the cilia movement and fluid transport emerge from specification of internal mechanics, viscous and elastic forces are reviewed. Finally, we discuss the limitations of existing work, and the challenges for the next generation of models, which may provide further insight into this complex and vital system

Outflow boundary conditions for 3D simulations of non-periodic blood flow and pressure fields in deformable arteries

The simulation of blood flow and pressure in arteries requires outflow boundary conditions that incorporate models of downstream domains. We previously described a coupled multidomain method to couple analytical models of the downstream domains with 3D numerical models of the upstream vasculature. This prior work either included pure resistance boundary conditions or impedance boundary conditions based on assumed periodicity of the solution. However, flow and pressure in arteries are not necessarily periodic in time due to heart rate variability, respiration, complex transitional flow or acute physiological changes. We present herein an approach for prescribing lumped parameter outflow boundary conditions that accommodate transient phenomena. We have applied this method to compute haemodynamic quantities in different physiologically relevant cardiovascular models, including patient-specific examples, to study non-periodic flow phenomena often observed in normal subjects and in patients with acquired or congenital cardiovascular disease. The relevance of using boundary conditions that accommodate transient phenomena compared with boundary conditions that assume periodicity of the solution is discussed

Optimization techniques in respiratory control system models

One of the most complex physiological systems whose modeling is still an open study is the respiratory control system where different models have been proposed based on the criterion of minimizing the work of breathing (WOB). The aim of this study is twofold: to compare two known models of the respiratory control system which set the breathing pattern based on quantifying the respiratory work; and to assess the influence of using direct-search or evolutionary optimization algorithms on adjustment of model parameters. This study was carried out using experimental data from a group of healthy volunteers under CO2 incremental inhalation, which were used to adjust the model parameters and to evaluate how much the equations of WOB follow a real breathing pattern. This breathing pattern was characterized by the following variables: tidal volume, inspiratory and expiratory time duration and total minute ventilation. Different optimization algorithms were considered to determine the most appropriate model from physiological viewpoint. Algorithms were used for a double optimization: firstly, to minimize the WOB and secondly to adjust model parameters. The performance of optimization algorithms was also evaluated in terms of convergence rate, solution accuracy and precision. Results showed strong differences in the performance of optimization algorithms according to constraints and topological features of the function to be optimized. In breathing pattern optimization, the sequential quadratic programming technique (SQP) showed the best performance and convergence speed when respiratory work was low. In addition, SQP allowed to implement multiple non-linear constraints through mathematical expressions in the easiest way. Regarding parameter adjustment of the model to experimental data, the evolutionary strategy with covariance matrix and adaptation (CMA-ES) provided the best quality solutions with fast convergence and the best accuracy and precision in both models. CMAES reached the best adjustment because of its good performance on noise and multi-peaked fitness functions. Although one of the studied models has been much more commonly used to simulate respiratory response to CO2 inhalation, results showed that an alternative model has a more appropriate cost function to minimize WOB from a physiological viewpoint according to experimental data.Postprint (author's final draft

Advances in computational modelling for personalised medicine after myocardial infarction

Myocardial infarction (MI) is a leading cause of premature morbidity and mortality worldwide. Determining which patients will experience heart failure and sudden cardiac death after an acute MI is notoriously difficult for clinicians. The extent of heart damage after an acute MI is informed by cardiac imaging, typically using echocardiography or sometimes, cardiac magnetic resonance (CMR). These scans provide complex data sets that are only partially exploited by clinicians in daily practice, implying potential for improved risk assessment. Computational modelling of left ventricular (LV) function can bridge the gap towards personalised medicine using cardiac imaging in patients with post-MI. Several novel biomechanical parameters have theoretical prognostic value and may be useful to reflect the biomechanical effects of novel preventive therapy for adverse remodelling post-MI. These parameters include myocardial contractility (regional and global), stiffness and stress. Further, the parameters can be delineated spatially to correspond with infarct pathology and the remote zone. While these parameters hold promise, there are challenges for translating MI modelling into clinical practice, including model uncertainty, validation and verification, as well as time-efficient processing. More research is needed to (1) simplify imaging with CMR in patients with post-MI, while preserving diagnostic accuracy and patient tolerance (2) to assess and validate novel biomechanical parameters against established prognostic biomarkers, such as LV ejection fraction and infarct size. Accessible software packages with minimal user interaction are also needed. Translating benefits to patients will be achieved through a multidisciplinary approach including clinicians, mathematicians, statisticians and industry partners