A numerical model for inert gas transport in the lung based on fractal airway morphology

Patients suffering from cystic fibrosis (CF) show thick secretions, mucus plugging and bronchiectasis in bronchial and alveolar ducts. This results in substantial structural changes of the airway morphology and heterogeneous ventilation. Disease progression and treatment effects are monitored by so-called gas washout tests, where the change in concentration of an inert gas is measured over a single or multiple breaths. The result of the tests based on the profile of the measured concentration is a marker for the severity of the ventilation inhomogeneity strongly affected by the airway morphology. However, it is hard to localize underlying obstructions to specific parts of the airways, especially if occurring in the lung periphery. In order to support the analysis of lung function tests (e.g. multi-breath washout), we developed a numerical model of the entire airway tree, coupling a lumped parameter model for the lung ventilation with a 4th-order accurate finite difference model of a 1D advection-diffusion equation for the transport of an inert gas. The boundary conditions for the flow problem comprise the pressure and flow profile at the mouth, which is typically known from clinical washout tests. The natural asymmetry of the lung morphology is approximated by a generic, fractal, asymmetric branching scheme which we applied for the conducting airways. A conducting airway ends when its dimension falls below a predefined limit. A model acinus is then connected to each terminal airway. The morphology of an acinus unit comprises a network of expandable cells. A regional, linear constitutive law describes the pressure-volume relation between the pleural gap and the acinus. The cyclic expansion (breathing) of each acinus unit depends on the resistance of the feeding airway and on the flow resistance and stiffness of the cells themselves. Special care was taken in the development of a conservative numerical scheme for the gas transport across bifurcations, handling spatially and temporally varying advective and diffusive fluxes over a wide range of scales. Implicit time integration was applied to account for the numerical stiffness resulting from the discretized transport equation. Local or regional modification of the airway dimension, resistance or tissue stiffness are introduced to mimic pathological airway restrictions typical for CF. This leads to a more heterogeneous ventilation of the model lung. As a result the concentration in some distal parts of the lung model remains increased for a longer duration. The inert gas concentration at the mouth towards the end of the expirations is composed of gas from regions with very different washout efficiency. This results in a steeper slope of the corresponding part of the washout profile.

Transport of volatile chlorinated hydrocarbons in unsaturated aggregated media

Transport of volatile hydrocarbons in soils is largely controlled by interactions of vapours with the liquid and solid phase. Sorption on solids of gaseous or dissolved comPounds may be important. Since the contact time between a chemical and a specific sorption site can be rather short, kinetic or mass-transfer resistance effects may be relevant. An existing mathematical model describing advection and diffusion in the gas phase and diffusional transport from the gaseous phase into an intra-aggregate water phase is modified to include linear kinetic sorption on ps-solid and water-solid interfaces. The model accounts for kinetic mass transfer between all three phases in a soil. The solution of the Laplace-transformed equations is inverted numerically. We performed transient column experiments with 1,1,2-Trichloroethane, Trichloroethylene, and Tetrachloroethylene using air-dry solid and water-saturated porous glass beads. The breakthrough curves were calculated based on independently estimated parameters. The model calculations agree well with experimental data. The different transport behaviour of the three compounds in our system primarily depends on Henry's constants.

Flux and resident injection in gaseous advection experiments

The occurrence of gaseous pollutants in soils has stimulated many experimental activities, including forced ventilation in the field as well as laboratory transport experiments with gases. The dispersion coefficient in advective-dispersive gas phase transport is often dominated by molecular diffusion, which leads to a large overall dispersivity gamma. Under such conditions it is important to distinguish between flux and resident modes of solute injection and detection. The influence of the inlet type oil the macroscopic injection mode was tested in two series of column experiments with gases at different mean flow velocities nu. First we compared infinite resident and flux injections, and second, semi-infinite resident and flux injections. It is shown that the macroscopically apparent injection condition depends on the geometry of the inlet section. A reduction of the cross-sectional area of the inlet relative to that of the column is very effective in excluding the diffusive solute input, thus allowing us to use the solutions for a flux Injection also at rather low mean flow velocities nu. If the whole cross section of a column is exposed to a large reservoir like that of ambient air, a semi-infinite resident injection is established, which can be distinguished from a flux injection even at relatively high velocities nu, depending on the mechanical dispersivity of the porous medium.

Optimized CUDA-Based PDE Solver for Reaction Diffusion Systems on Arbitrary Surfaces

Partial differential equation (PDE) solvers are commonly employed to study and characterize the parameter space for reaction-diffusion (RD) systems while investigating biological pattern formation. Increasingly, biologists wish to perform such studies with arbitrary surfaces representing ‘real’ 3D geometries for better insights. In this paper, we present a highly optimized CUDA-based solver for RD equations on triangulated meshes in 3D. We demonstrate our solver using a chemotactic model that can be used to study snakeskin pigmentation, for example. We employ a finite element based approach to perform explicit Euler time integrations. We compare our approach to a naive GPU implementation and provide an in-depth performance analysis, demonstrating the significant speedup afforded by our optimizations. The optimization strategies that we exploit could be generalized to other mesh based processing applications with PDE simulations.