Three-dimensional micro-droplet collision simulation using the Lattice Boltzmann method

This paper was presented at the 3rd Micro and Nano Flows Conference (MNF2011), which was held at the Makedonia Palace Hotel, Thessaloniki in Greece. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, Aristotle University of Thessaloniki, University of Thessaly, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute.The modelling of binary droplet collisions has important applications in many engineering problems, including spray coating and fuel injection. The Lattice Boltzmann method (LBM) is a well established technique for modelling multiphase fluids, and does so without the difficulties of explicit interface tracking found in other CFD methods. However, simulating droplet collisions under realistic conditions remains a complex problem. Challenges include reproducing the different collision outcomes observed experimentally (Qian and Law, 1997), and maintaining a stable simulation at sufficiently high Reynolds and Weber numbers, and with a high density ratio between the liquid and gas phases. Although previous studies have achieved these goals individually, they have not been successfully combined to simulate droplet collisions with realistic physical parameters. A number of different methods for extending the LBM for multiphase flow exist, with the Shan-Chen interparticle potential method (Shan and Chen, 1993) being the basic model used here. Many extensions to improve the original Shan-Chen method have been proposed, to improve achievable Reynolds number and density ratio. Using combinations of these, both coalescence and separation of two-dimensional droplets were successfully simulated at density ratios of order 1000, and high Weber numbers (Lycett-Brown et al., 2011). In this study, the developed methodologies in Lycett-Brown et al. (2011) are extended to simulate three dimensional micro-droplet collisions by making use of the LBM’s excellent scalability on massively parallel computers. These high-resolution simulations are also compared with low-resolution three-dimensional simulations using a multiple-relaxation-time LBM approach (Monaco and Luo, 2008).This study is funded by the Engineering and Physical Sciences Research Council for Grant No. EP/I000801/1 and a HEC Studentship

How are adjacent spinal levels affected by vertebral fracture and by vertebroplasty? A biomechanical study on cadaveric spines

Background Context Spinal injuries and surgery may have important effects on neighboring spinal levels, but previous investigations of adjacent-level biomechanics have produced conflicting results. We use “stress profilometry” and noncontact strain measurements to investigate thoroughly this long-standing problem. Purpose This study aimed to determine how vertebral fracture and vertebroplasty affect compressive load-sharing and vertebral deformations at adjacent spinal levels. Study Design We conducted mechanical experiments on cadaver spines. Methods Twenty-eight cadaveric spine specimens, comprising three thoracolumbar vertebrae and the intervening discs and ligaments, were dissected from fourteen cadavers aged 67–92 years. A needle-mounted pressure transducer was used to measure the distribution of compressive stress across the anteroposterior diameter of both intervertebral discs. “Stress profiles” were analyzed to quantify intradiscal pressure (IDP) and concentrations of compressive stress in the anterior and posterior annulus. Summation of stresses over discrete areas yielded the compressive force acting on the anterior and posterior halves of each vertebral body, and the compressive force resisted by the neural arch. Creep deformations of vertebral bodies under load were measured using an optical MacReflex system. All measurements were repeated following compressive injury to one of the three vertebrae, and again after the injury had been treated by vertebroplasty. The study was funded by a grant from Action Medical Research, UK ($143,230). Authors of this study have no conflicts of interest to disclose. Results Injury usually involved endplate fracture, often combined with deformation of the anterior cortex, so that the affected vertebral body developed slight anterior wedging. Injury reduced IDP at the affected level, to an average 47% of pre-fracture values (p<.001), and transferred compressive load-bearing from nucleus to annulus, and also from disc to neural arch. Similar but reduced effects were seen at adjacent (non-fractured) levels, where mean IDP was reduced to 73% of baseline values (p<.001). Vertebroplasty partially reversed these changes, increasing mean IDP to 76% and 81% of baseline values at fractured and adjacent levels, respectively. Injury also increased creep deformation of the vertebral body under load, especially in the anterior region where a 14-fold increase was observed at the fractured level and a threefold increase was observed at the adjacent level. Vertebroplasty also reversed these changes, reducing deformation of the anterior vertebral body (compared with post-fracture values) by 62% at the fractured level, and by 52% at the adjacent level. Conclusions Vertebral fracture adversely affects compressive load-sharing and increases vertebral deformations at both fractured and adjacent levels. All effects can be partially reversed by vertebroplasty

Vertebroplasty reduces progressive ׳creep' deformity of fractured vertebrae

Elderly vertebrae frequently develop an “anterior wedge” deformity as a result of fracture and creep mechanisms. Injecting cement into a damaged vertebral body (vertebroplasty) is known to help restore its shape and stiffness. We now hypothesise that vertebroplasty is also effective in reducing subsequent creep deformations. Twenty-eight spine specimens, comprising three complete vertebrae and the intervening discs, were obtained from cadavers aged 67–92 years. Each specimen was subjected to increasingly-severe compressive loading until one of its vertebrae was fractured, and the damaged vertebral body was then treated by vertebroplasty. Before and after fracture, and again after vertebroplasty, each specimen was subjected to a static compressive force of 1 kN for 1 h while elastic and creep deformations were measured in the anterior, middle and posterior regions of each adjacent vertebral body cortex, using a 2D MacReflex optical tracking system. After fracture, creep in the anterior and central regions of the vertebral body cortex increased from an average 4513 and 885 microstrains, respectively, to 54,107 and 34,378 microstrains (both increases: P<0.001). Elastic strains increased by a comparable amount. Vertebroplasty reduced creep in the anterior and central cortex by 61% (P=0.006) and 66% (P=0.017) respectively. Elastic strains were reduced by less than half this amount. Results suggest that the beneficial effects of vertebroplasty on the vertebral body continue long after the post-operative radiographs. Injected cement not only helps to restore vertebral shape and elastic properties, but also reduces subsequent creep deformation of the damaged vertebra

Improved forcing scheme in pseudopotential lattice Boltzmann methods for multiphase flow at arbitrarily high density ratios

The pseudopotential lattice Boltzmann method has been widely used to simulate many multiphase flow applications. However, there still exist problems with reproducing realistic values of density ratio and surface tension. In this study, a higher-order analysis of a general forcing term is derived. A forcing scheme is then constructed for the pseudopotential method that is able to accurately reproduce the full range of coexistence curves. As a result, multiphase flow of arbitrarily high density ratios independent of the surface tension can be simulated. Furthermore, the interface width can be tuned to allow for grid refinement and systematic error reduction. Numerical results confirm that the proposed scheme enables independent control of density ratio, surface tension, and interface width simultaneously

ELASTOGRAPHY CAN EFFECTIVELY DECREASE THE NUMBER OF FINE-NEEDLE ASPIRATION BIOPSIES IN PATIENTS WITH CALCIFIED THYROID NODULES

When calcification, frequently found in both benign and malignant nodules, is present in thyroid nodules, non-invasive differentiation with ultrasound becomes challenging. The goal of this study was to evaluate the utility of elastography in differentiating calcified thyroid nodules. Consecutive patients (165 patients with 196 nodules) referred for fine-needle aspiration who had undergone both ultrasound elastography and B-mode examinations were analyzed retrospectively. Calcification was present in 45 benign and 20 malignant nodules. On 65 calcified nodules, elastography had 95% sensitivity, 51.1% specificity, 46.3% positive predictive value and 95.8% negative predictive value in detecting malignancy. Twenty-three of 45 benign calcified nodules were correctly diagnosed with elastography compared with 4 of 45 by B-mode ultrasound. Although it is difficult to differentiate benign and malignant calcified thyroid nodules solely with B-mode ultrasound, elastography has the potential to reduce the number of fine-needle aspiration biopsies performed on calcified nodules. (C) 2014 World Federation for Ultrasound in Medicine & Biology.1188Ysciescopu

Magnetically responsive layer-by-layer microcapsules can be retained in cells and under flow conditions to promote local drug release without triggering ROS production.

Nanoengineered vehicles have the potential to deliver cargo drugs directly to disease sites, but can potentially be cleared by immune system cells or lymphatic drainage. In this study we explore the use of magnetism to hold responsive particles at a delivery site, by incorporation of superparamagnetic iron oxide nanoparticles (SPIONs) into layer-by-layer (LbL) microcapsules. Microcapsules with SPIONs were rapidly phagocytosed by cells but did not trigger cellular ROS synthesis within 24 hours of delivery nor affect cell viability. In a non-directional cell migration assay, SPION containing microcapsules significantly inhibited movement of phagocytosing cells when placed in a magnetic field. Similarly, under flow conditions, a magnetic field retained SPION containing microcapsules at a physiologic wall shear stress of 0.751 dyne cm-2. Even when the SPION content was reduced to 20%, the majority of microcapsules were still retained. Dexamethasone microcrystals were synthesised by solvent evaporation and underwent LbL encapsulation with inclusion of a SPION layer. Despite a lower iron to volume content of these structures compared to microcapsules, they were also retained under shear stress conditions and displayed prolonged release of active drug, beyond 30 hours, measured using a glucocorticoid sensitive reporter cell line generated in this study. Our observations suggest use of SPIONs for magnetic retention of LbL structures is both feasible and biocompatible and has potential application for improved local drug delivery

3D Printed Franz cells - update on optimization of manufacture and evaluation

The evaluation of permeation profiles from cosmetic formulations is considered to be a crucial component in both the development and quality assurance of any new product [1, 2]. Data gathered from such studies allow researchers to assess the viability of delivering different materials to and through biological membranes. To date, laboratory in vitro permeation processes require the use of modified Franz type diffusion cells, conventionally fabricated from glass, which are available in different formats that can be customised to experimental requirements [3]

The effects of dynamic compression on the development of cartilage grafts engineered using bone marrow and infrapatellar fat pad derived stem cells.

Bioreactors that subject cell seeded scaffolds or hydrogels to biophysical stimulation have been used to improve the functionality of tissue engineered cartilage and to explore how such constructs might respond to the application of joint specific mechanical loading. Whether a particular cell type responds appropriately to physiological levels of biophysical stimulation could be considered a key determinant of its suitability for cartilage tissue engineering applications. The objective of this study was to determine the effects of dynamic compression on chondrogenesis of stem cells isolated from different tissue sources. Porcine bone marrow (BM) and infrapatellar fat pad (FP) derived stem cells were encapsulated in agarose hydrogels and cultured in a chondrogenic medium in free swelling (FS) conditions for 21 d, after which samples were subjected to dynamic compression (DC) of 10% strain (1 Hz, 1 h d(-1)) for a further 21 d. Both BM derived stem cells (BMSCs) and FP derived stem cells (FPSCs) were capable of generating cartilaginous tissues with near native levels of sulfated glycosaminoglycan (sGAG) content, although the spatial development of the engineered grafts strongly depended on the stem cell source. The mechanical properties of cartilage grafts generated from both stem cell sources also approached that observed in skeletally immature animals. Depending on the stem cell source and the donor, the application of DC either enhanced or had no significant effect on the functional development of cartilaginous grafts engineered using either BMSCs or FPSCs. BMSC seeded constructs subjected to DC stained less intensely for collagen type I. Furthermore, histological and micro-computed tomography analysis showed mineral deposition within BMSC seeded constructs was suppressed by the application of DC. Therefore, while the application of DC in vitro may only lead to modest improvements in the mechanical functionality of cartilaginous grafts, it may play an important role in the development of phenotypically stable constructs.Funding was provided by the European Research Council Starter Grant (StemRepair—Project number 258463) and a SFI President of Ireland Young Researcher Award (08/Y15/B1336)

Wave attenuation at a salt marsh margin: A case study of an exposed coast on the Yangtze estuary

To quantify wave attenuation by (introduced) Spartina alterniflora vegetation at an exposed macrotidal coast in the Yangtze Estuary, China, wave parameters and water depth were measured during 13 consecutive tides at nine locations ranging from 10 m seaward to 50 m landward of the low marsh edge. During this period, the incident wave height ranged from <0.1 to 1.5 m, the maximum of which is much higher than observed in other marsh areas around the world. Our measurements and calculations showed that the wave attenuation rate per unit distance was 1 to 2 magnitudes higher over the marsh than over an adjacent mudflat. Although the elevation gradient of the marsh margin was significantly higher than that of the adjacent mudflat, more than 80% of wave attenuation was ascribed to the presence of vegetation, suggesting that shoaling effects were of minor importance. On average, waves reaching the marsh were eliminated over a distance of similar to 80 m, although a marsh distance of >= 100 m was needed before the maximum height waves were fully attenuated during high tides. These attenuation distances were longer than those previously found in American salt marshes, mainly due to the macrotidal and exposed conditions at the present site. The ratio of water depth to plant height showed an inverse correlation with wave attenuation rate, indicating that plant height is a crucial factor determining the efficiency of wave attenuation. Consequently, the tall shoots of the introduced S. alterniflora makes this species much more efficient at attenuating waves than the shorter, native pioneer species in the Yangtze Estuary, and should therefore be considered as a factor in coastal management during the present era of sea-level rise and global change. We also found that wave attenuation across the salt marsh can be predicted using published models when a suitable coefficient is incorporated to account for drag, which varies in place and time due to differences in plant characteristics and abiotic conditions (i.e., bed gradient, initial water depth, and wave action).

Robustness and fault tolerance make brains harder to study

Brains increase the survival value of organisms by being robust and fault tolerant. That is, brain circuits continue to operate as the organism needs, even when the circuit properties are significantly perturbed. Kispersky and colleagues, in a recent paper in Neural Systems & Circuits, have found that Granger Causality analysis, an important method used to infer circuit connections from the behavior of neurons within the circuit, is defeated by the mechanisms that give rise to this robustness and fault tolerance