A new membrane formulation for modelling the flow of stomatocyte, discocyte, and echinocyte red blood cells.

In this work, a numerical model that enables simulation of the deformation and flow behaviour of differently aged Red Blood Cells (RBCs) is developed. Such cells change shape and decrease in deformability as they age, thus impacting their ability to pass through the narrow capillaries in the body. While the body filters unviable cells from the blood naturally, cell aging poses key challenges for blood stored for transfusions. Therefore, understanding the influence RBC morphology and deformability have on their flow is vital. While several existing models represent young Discocyte RBC shapes well, a limited number of numerical models are developed to model aged RBC morphologies like Stomatocytes and Echinocytes. The existing models are also limited to shear and stretching simulations. Flow characteristics of these morphologies are yet to be investigated. This paper aims to develop a new membrane formulation for the numerical modelling of Stomatocyte, Discocytes and Echinocyte RBC morphologies to investigate their deformation and flow behaviour. The model used represents blood plasma using the Lattice Boltzmann Method (LBM) and the RBC membrane using the discrete element method (DEM). The membrane and the plasma are coupled by the Immersed Boundary Method (IBM). Previous LBM-IBM-DEM formulations represent RBC membrane response based on forces generated from changes in the local area, local length, local bending, and cell volume. In this new model, two new force terms are added: the local area difference force and the local curvature force, which are specially incorporated to model the flow and deformation behaviour of Stomatocytes and Echinocytes. To verify the developed model, the deformation behaviour of the three types of RBC morphologies are compared to well-characterised stretching and shear experiments. The flow modelling capabilities of the method are then demonstrated by modelling the flow of each cell through a narrow capillary. The developed model is found to be as accurate as benchmark Smoothed Particle Hydrodynamics (SPH) approaches while being significantly more computationally efficient.

A three-dimensional (3-D) meshfree-based computational model to investigate stress-strain-time relationships of plant cells during drying.

A better understanding of plant cell micromechanics would enhance the current opinion on "how things are happening" inside a plant cell, enabling more detailed insights into plant physiology as well as processing plant biomaterials. However, with the contemporary laboratory equipment, the experimental investigation of cell micromechanics has been a challenging task due to diminutive spatial and time scales involved. In this investigation, a three-dimensional (3-D) coupled Smoothed Particle Hydrodynamics (SPH) and Coarse-Grained (CG) computational approach has been employed to model micromechanics of single plant cells going through drying or dehydration. This meshfree-based computational model has conclusively demonstrated that it can effectively simulate the behaviour of stress and strain in a plant cell being compressed at different levels of dryness: ranging from a fresh state to an extremely dried state. In addition, different biological and physical circumstances have been approximated through the proposed novel computational framework in the form of different turgor pressures, strain rates, mechanical properties and cell sizes. The proposed computational framework has potential not only to study the micromechanical characteristics of plant cellular structure during drying, but also other equivalent, biological structures and processes with relevant modifications. There are no underlying difficulties in adopting the model to replicate other types of cells and more sophisticated micromechanical phenomena of the cells under different external loading conditions.

Application of high-order lattice Boltzmann pseudopotential models.

Higher-order lattice Boltzmann (LB) pseudopotential models have great potential for solving complex fluid dynamics in various areas of modern science. The discreteness of the lattice discretization makes these models an attractive choice due to their flexibility, capacity to capture hydrodynamic details, and inherent adaptability to parallel computations. Despite those advantages, the discreteness makes high-order LB models difficult to apply due to the larger lattice structure, for which basic fundamental properties, namely diffusion coefficient and contact angle, remain unknown. This work addresses this by providing general continuum solutions for those two basic properties and demonstrating these solutions to compare favorably against known theory. Various high-order LB models are shown to reproduce the sinusoidal decay of a binary miscible mixture accurately and consistently. Furthermore, these models are shown to reproduce neutral, hydrophobic, and hydrophilic contact angles. Discrete differences are shown to exist, which are captured at the discrete level and confirmed through droplet shape analysis. This work provides practical tools that allow for high-order LB pseudopotential models to be used to simulate multicomponent flows.

Deformation behaviour of stomatocyte, discocyte and echinocyte red blood cell morphologies during optical tweezers stretching.

The red blood cell (RBC) deformability is a critical aspect, and assessing the cell deformation characteristics is essential for better diagnostics of healthy and deteriorating RBCs. There is a need to explore the connection between the cell deformation characteristics, cell morphology, disease states, storage lesion and cell shape-transformation conditions for better diagnostics and treatments. A numerical approach inspired from the previous research for RBC morphology predictions and for analysis of RBC deformations is proposed for the first time, to investigate the deformation characteristics of different RBC morphologies. The present study investigates the deformability characteristics of stomatocyte, discocyte and echinocyte morphologies during optical tweezers stretching and provides the opportunity to study the combined contribution of cytoskeletal spectrin network and the lipid-bilayer during RBC deformation. The proposed numerical approach predicts agreeable deformation characteristics of the healthy discocyte with the analogous experimental observations and is extended to further investigate the deformation characteristics of stomatocyte and echinocyte morphologies. In particular, the computer simulations are performed to investigate the influence of direct stretching forces on different equilibrium cell morphologies on cell spectrin link extensions and cell elongation index, along with a parametric analysis on membrane shear modulus, spectrin link extensibility, bending modulus and RBC membrane-bead contact diameter. The results agree with the experimentally observed stiffer nature of stomatocyte and echinocyte with respect to a healthy discocyte at experimentally determined membrane characteristics and suggest the preservation of relevant morphological characteristics, changes in spectrin link densities and the primary contribution of cytoskeletal spectrin network on deformation behaviour of stomatocyte, discocyte and echinocyte morphologies during optical tweezers stretching deformation. The numerical approach presented here forms the foundation for investigations into deformation characteristics and recoverability of RBCs undergoing storage lesion.

A Review of Respiratory Anatomical Development, Air Flow Characterization and Particle Deposition.

The understanding of complex inhalation and transport processes of pollutant particles through the human respiratory system is important for investigations into dosimetry and respiratory health effects in various settings, such as environmental or occupational health. The studies over the last few decades for micro- and nanoparticle transport and deposition have advanced the understanding of drug-aerosol impacts in the mouth-throat and the upper airways. However, most of the Lagrangian and Eulerian studies have utilized the non-realistic symmetric anatomical model for airflow and particle deposition predictions. Recent improvements to visualization techniques using high-resolution computed tomography (CT) data and the resultant development of three dimensional (3-D) anatomical models support the realistic representation of lung geometry. Yet, the selection of different modelling approaches to analyze the transitional flow behavior and the use of different inlet and outlet conditions provide a dissimilar prediction of particle deposition in the human lung. Moreover, incorporation of relevant physical and appropriate boundary conditions are important factors to consider for the more accurate prediction of transitional flow and particle transport in human lung. This review critically appraises currently available literature on airflow and particle transport mechanism in the lungs, as well as numerical simulations with the aim to explore processes involved. Numerical studies found that both the Euler-Lagrange (E-L) and Euler-Euler methods do not influence nanoparticle (particle diameter ≤50 nm) deposition patterns at a flow rate ≤25 L/min. Furthermore, numerical studies demonstrated that turbulence dispersion does not significantly affect nanoparticle deposition patterns. This critical review aims to develop the field and increase the state-of-the-art in human lung modelling.

Euler-Lagrange Prediction of Diesel-Exhaust Polydisperse Particle Transport and Deposition in Lung: Anatomy and Turbulence Effects.

In clinical assessments, the correlation between atmospheric air pollution and respiratory damage is highly complicated. Epidemiological studies show that atmospheric air pollution is largely responsible for the global proliferation of pulmonary disease. This is particularly significant, since most Computational Fluid Dynamics (CFD) studies to date have used monodisperse particles, which may not accurately reflect realistic inhalation patterns, since atmospheric aerosols are mostly polydisperse. The aim of this study is to investigate the anatomy and turbulent effects on polydisperse particle transport and deposition (TD) in the upper airways. The Euler-Lagrange approach is used for polydisperse particle TD prediction in both laminar and turbulent conditions. Various anatomical models are adopted to investigate the polydisperse particle TD under different flow conditions. Rossin-Rammler diameter distribution is used for the distribution of the initial particle diameter. The numerical results illustrate that airflow rate distribution at the right lung of a realistic model is higher than a non-realistic model. The CFD study also shows that turbulence effects on deposition are higher for larger diameter particles than with particles of smaller diameter. A significant amount of polydisperse particles are also shown to be deposited at the tracheal wall for CT-based model, whereas particles are mostly deposited at the carinal angle for the non-realistic model. A comprehensive, polydisperse particle TD analysis would enhance understanding of the realistic deposition pattern and decrease unwanted therapeutic aerosol deposition at the extrathoracic airways.

Interaction pressure tensor on high-order lattice Boltzmann models for nonideal fluids.

In this work we address the application of pseudopotentials directly on high-order lattice Boltzmann models. We derive a general expression for the pressure tensor on high-order lattices considering all nonideal interactions, including intra- and intermolecular interactions, following the discrete lattice theory introduced by X. Shan [Phys. Rev. E 77, 066702 (2008)PLEEE81539-375510.1103/PhysRevE.77.066702]. From the derived expression, a generalized continuum approximation, truncated at fourth-order isotropy, is obtained that is readily applicable to high-order lattices. With this, we demonstrate that high-order lattice models with pseudopotentials can satisfy thermodynamic consistency. The derived generalized expression and continuum approximation are validated for the case of a flat interface and compared against the standard definition available from the literature. The generalized expression is also shown to accurately reproduce the Laplace experiment for a variety of high-order lattice structures. This work sets the preliminary steps towards the application of high-order lattice models for simulating nonideal fluid mixtures.

Molecular insights on the interference of simplified lung surfactant models by gold nanoparticle pollutants.

Inhaled nanoparticles (NPs) are experienced by the first biological barrier inside the alveolus known as lung surfactant (LS), a surface tension reducing agent, consisting of phospholipids and proteins in the form of the monolayer at the air-water interface. The monolayer surface tension is continuously regulated by the alveolus compression and expansion and protects the alveoli from collapsing. Inhaled NPs can reach deep into the lungs and interfere with the biophysical properties of the lung components. The interaction mechanisms of bare gold nanoparticles (AuNPs) with the LS monolayer and the consequences of the interactions on lung function are not well understood. Coarse-grained molecular dynamics simulations were carried out to elucidate the interactions of AuNPs with simplified LS monolayers at the nanoscale. It was observed that the interactions of AuNPs and LS components deform the monolayer structure, change the biophysical properties of LS and create pores in the monolayer, which all interfere with the normal lungs function. The results also indicate that AuNP concentrations >0.1 mol% (of AuNPs/lipids) hinder the lowering of the LS surface tension, a prerequisite of the normal breathing process. Overall, these findings could help to identify the possible consequences of airborne NPs inhalation and their contribution to the potential development of various lung diseases.

A coarse-grained multiscale model to simulate morphological changes of food-plant tissues undergoing drying.

Numerical modelling has emerged as a powerful and effective tool to study various dynamic behaviours of biological matter. Such numerical modelling tools have contributed to the optimisations of food drying parameters leading to higher quality end-products in the field of food engineering. In this context, one of the most recent developments is the meshfree-based numerical models, which have demonstrated enhanced capabilities to model cellular deformations during drying, providing many benefits compared to conventional grid-based modelling approaches. However, the potential extension of this method for simulating bulk level tissues has been a challenge due to the increased requirement for higher computational time and resources. As a solution for this, by incorporating meshfree features, a novel coarse-grained multiscale numerical model is proposed in this work to predict bulk level (macroscale) deformations of food-plant tissues during drying. Accordingly, realistic simulation of morphological changes of apple tissues can now be performed with just 2% of the previous computational time in microscale and macroscale simulations can also be conducted. Compared to contemporary multiscale models, this modelling approach provides more convenient computational implementation as well. Thus, this novel approach can be recommended for efficiently and accurately simulating morphological changes of cellular materials undergoing drying processes, while confirming its potential future expansion to efficiently and accurately predict morphological changes of heterogeneous plant tissues in different spatial scales.

Polydisperse Microparticle Transport and Deposition to the Terminal Bronchioles in a Heterogeneous Vasculature Tree.

The atmospheric particles from different sources, and the therapeutic particles from various drug delivery devices, exhibit a complex size distribution, and the particles are mostly polydisperse. The limited available in vitro, and the wide range of in silico models have improved understanding of the relationship between monodisperse particle deposition and therapeutic aerosol transport. However, comprehensive polydisperse transport and deposition (TD) data for the terminal airways is still unavailable. Therefore, to benefit future drug therapeutics, the present numerical model illustrates detailed polydisperse particle TD in the terminal bronchioles for the first time. Euler-Lagrange approach and Rosin-Rammler diameter distribution is used for polydisperse particles. The numerical results show higher deposition efficiency (DE) in the right lung. Specifically, the larger the particle diameter (dp > 5 µm), the higher the DE at the bifurcation area of the upper airways is, whereas for the smaller particle (dp < 5 µm), the DE is higher at the bifurcation wall. The overall deposition pattern shows a different deposition hot spot for different diameter particle. These comprehensive lobe-specific polydisperse particle deposition studies will increase understanding of actual inhalation for particle TD, which could potentially increase the efficiency of pharmaceutical aerosol delivery at the targeted position of the terminal airways.

Application of a coupled smoothed particle hydrodynamics (SPH) and coarse-grained (CG) numerical modelling approach to study three-dimensional (3-D) deformations of single cells of different food-plant materials during drying.

Numerical modelling has gained popularity in many science and engineering streams due to the economic feasibility and advanced analytical features compared to conventional experimental and theoretical models. Food drying is one of the areas where numerical modelling is increasingly applied to improve drying process performance and product quality. This investigation applies a three dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) and Coarse-Grained (CG) numerical approach to predict the morphological changes of different categories of food-plant cells such as apple, grape, potato and carrot during drying. To validate the model predictions, experimental findings from in-house experimental procedures (for apple) and sources of literature (for grape, potato and carrot) have been utilised. The subsequent comaprison indicate that the model predictions demonstrate a reasonable agreement with the experimental findings, both qualitatively and quantitatively. In this numerical model, a higher computational accuracy has been maintained by limiting the consistency error below 1% for all four cell types. The proposed meshfree-based approach is well-equipped to predict the morphological changes of plant cellular structure over a wide range of moisture contents (10% to 100% dry basis). Compared to the previous 2-D meshfree-based models developed for plant cell drying, the proposed model can draw more useful insights on the morphological behaviour due to the 3-D nature of the model. In addition, the proposed computational modelling approach has a high potential to be used as a comprehensive tool in many other tissue morphology related investigations.

Ultrafine particle transport and deposition in a large scale 17-generation lung model.

To understand how to assess optimally the risks of inhaled particles on respiratory health, it is necessary to comprehend the uptake of ultrafine particulate matter by inhalation during the complex transport process through a non-dichotomously bifurcating network of conduit airways. It is evident that the highly toxic ultrafine particles damage the respiratory epithelium in the terminal bronchioles. The wide range of in silico available and the limited realistic model for the extrathoracic region of the lung have improved understanding of the ultrafine particle transport and deposition (TD) in the upper airways. However, comprehensive ultrafine particle TD data for the real and entire lung model are still unavailable in the literature. Therefore, this study is aimed to provide an understanding of the ultrafine particle TD in the terminal bronchioles for the development of future therapeutics. The Euler-Lagrange (E-L) approach and ANSYS fluent (17.2) solver were used to investigate ultrafine particle TD. The physical conditions of sleeping, resting, and light activity were considered in this modelling study. A comprehensive pressure-drop along five selected path lines in different lobes was calculated. The non-linear behaviour of pressure-drops is observed, which could aid the health risk assessment system for patients with respiratory diseases. Numerical results also showed that ultrafine particle-deposition efficiency (DE) in different lobes is different for various physical activities. Moreover, the numerical results showed hot spots in various locations among the different lobes for different flow rates, which could be helpful for targeted therapeutical aerosol transport to terminal bronchioles and the alveolar region.

Angiotensin II type 1 receptor gene A1166C polymorphism and breast cancer susceptibility.

Numerous studies have evaluated the association between the angiotensin II type-1 receptor (AGTR1) gene A1166C polymorphism and breast cancer risk. However, the specific association is controversial. The aim of the present study was to derive a more precise estimation of the relationship. A comprehensive research was conducted of the PubMed and the Google Scholar databases through February 2015. Data were assessed using STATA version 12.0. Pooled odds ratios with 95%CIs were derived from the fixed-effect or random-effect models. A total of 911 patients with breast cancer and 1284 controls from 5 case-control studies were included in this meta-analysis. The meta-analysis results showed no significant association between the AGTR1 gene A1166C polymorphism and breast cancer risk. Similarly, in the subgroup analysis regarding ethnicity, no associations were observed. Heterogeneity and publication bias were not observed in this meta-analysis. The A1166C polymorphism in the AGTR1 gene may not be a risk factor for breast cancer. Further, large, and well-designed studies are needed to confirm this conclusion.

Effect of joint mechanism on vehicle redirectional capability of water-filled road safety barrier systems.

Portable water-filled barriers (PWFBs) are roadside appurtenances that prevent vehicles from penetrating into temporary construction zones on roadways. PWFBs are required to satisfy the strict regulations for vehicle re-direction in tests. However, many of the current PWFBs fail to re-direct the vehicle at high speeds due to the inability of the joints to provide appropriate stiffness. The joint mechanism hence plays a crucial role in the performance of a PWFB system at high speed impacts. This paper investigates the desired features of the joint mechanism in a PWFB system that can re-direct vehicles at high speeds, while limiting the lateral displacement to acceptable limits. A rectangular "wall" representative of a 30m long barrier system was modeled and a novel method of joining adjacent road barriers was introduced through appropriate pin-joint connections. The impact response of the barrier "wall" and the vehicle was obtained and the results show that a rotational stiffness of 3000kNm/rad at the joints seems to provide the desired features of the PWFB system to re-direct impacting vehicles and restrict the lateral deflection. These research findings will be useful to safety engineers and road barrier designers in developing a new generation of PWFBs for increased road safety.

A particle based model to simulate microscale morphological changes of plant tissues during drying.

Fundamental understanding on microscopic physical changes of plant materials is vital to optimize product quality and processing techniques, particularly in food engineering. Although grid-based numerical modelling can assist in this regard, it becomes quite challenging to overcome the inherited complexities of these biological materials especially when such materials undergo critical processing conditions such as drying, where the cellular structure undergoes extreme deformations. In this context, a meshfree particle based model was developed which is fundamentally capable of handling extreme deformations of plant tissues during drying. The model is built by coupling a particle based meshfree technique: Smoothed Particle Hydrodynamics (SPH) and a Discrete Element Method (DEM). Plant cells were initiated as hexagons and aggregated to form a tissue which also accounts for the characteristics of the middle lamella. In each cell, SPH was used to model cell protoplasm and DEM was used to model the cell wall. Drying was incorporated by varying the moisture content, the turgor pressure, and cell wall contraction effects. Compared to the state of the art grid-based microscale plant tissue drying models, the proposed model can be used to simulate tissues under excessive moisture content reductions incorporating cell wall wrinkling. Also, compared to the state of the art SPH-DEM tissue models, the proposed model better replicates real tissues and the cell-cell interactions used ensure efficient computations. Model predictions showed good agreement both qualitatively and quantitatively with experimental findings on dried plant tissues. The proposed modelling approach is fundamentally flexible to study different cellular structures for their microscale morphological changes at dehydration.

Molecular investigation of the mechanical properties of single actin filaments based on vibration analyses.

The mechanical vibration properties of single actin filaments from 50 to 288 nm are investigated by the molecular dynamics simulation in this study. The natural frequencies obtained from the molecular simulations agree with those obtained from the analytical solution of the equivalent Euler-Bernoulli beam model. Through the convergence study of the mechanical properties with respect to the filament length, it was found that the Euler-Bernoulli beam model can only be reliably used when the single actin filament is of the order of hundreds of nanometre scale. This molecular investigation not only provides the evidence for the use of the continuum beam model in characterising the mechanical properties of single actin filaments, but also clarifies the criteria for the effective use of the Euler-Bernoulli beam model.

Measurements of h(c)(1P(1)) in psi' decays.

We present measurements of the charmonium state h(c)(1P(1)) made with 106x10(6) psi' events collected by BESIII at BEPCII. Clear signals are observed for psi'-->pi0 h(c) with and without the subsequent radiative decay h(c)-->gamma eta(c). First measurements of the absolute branching ratios B(psi'-->pi0 h(c)) = (8.4+/-1.3+/-1.0) x 10(-4) and B(h(c)-->gamma eta(c)) = (54.3+/-6.7+/-5.2)% are presented. A statistics-limited determination of the previously unmeasured h(c) width leads to an upper limit Gamma(h(c))<1.44 MeV (90% confidence). Measurements of M(h(c)) = 3525.40+/-0.13+/-0.18 MeV/c2 and B(psi'-->pi0 h(c)) x B(h(c)-->gamma eta(c)) = (4.58+/-0.40+/-0.50) x 10(-4) are consistent with previous results.

Anomalous line shape of the cross section for e{+}e{-}--> hadrons in the center-of-mass energy region between 3.650 and 3.872 GeV.

We observe an obvious anomalous line shape of the e;{+}e;{-}--> hadrons total cross sections in the energy region between 3.700 and 3.872 GeV. It is inconsistent with the explanation for only one simple psi(3770) resonance with a statistical significance of 7sigma. The anomalous line shape may be explained by two possible enhancements of the inclusive hadron production near the center-of-mass energies of 3.764 and 3.779 GeV, indicating that either there is likely a new structure in addition to the psi(3770) resonance around 3.773 GeV, or there are some physics effects reflecting the DD[over ] production dynamics.

Search for the invisible decay of J/psi in psi(2S) --> pi(+)pi(-) J/psi.

Using psi(2S) --> pi(+)pi(-) J/psi events in a sample of 14.0 x 10(6) psi(2S) decays collected with the BES-II detector, a search for the decay of the J/psi to invisible final states is performed. No signal is found, and an upper limit at the 90% confidence level is determined to be 1.2 x 10(-2) for the ratio B(J/psi --> invisible)/B(J/psi-->mu(+)mu(-)). This is the first search for J/psi decays to invisible final states.