Coarsening dynamics of ternary amphiphilic fluids and the self-assembly of the gyroid and sponge mesophases: lattice-Boltzmann simulations

By means of a three-dimensional amphiphilic lattice-Boltzmann model with short-range interactions for the description of ternary amphiphilic fluids, we study how the phase separation kinetics of a symmetric binary immiscible fluid is altered by the presence of the amphiphilic species. We find that a gradual increase in amphiphile concentration slows down domain growth, initially from algebraic, to logarithmic temporal dependence, and, at higher concentrations, from logarithmic to stretched-exponential form. In growth-arrested stretched-exponential regimes, at late times we observe the self-assembly of sponge mesophases and gyroid liquid crystalline cubic mesophases, hence confirming that (a) amphiphile-amphiphile interactions need not be long-ranged in order for periodically modulated structures to arise in a dynamics of competing interactions, and (b) a chemically-specific model of the amphiphile is not required for the self-assembly of cubic mesophases, contradicting claims in the literature. We also observe a structural order-disorder transition between sponge and gyroid phases driven by amphiphile concentration alone or, independently, by the amphiphile-amphiphile and the amphiphile-binary fluid coupling parameters. For the growth-arrested mesophases, we also observe temporal oscillations in the structure function at all length scales; most of the wavenumbers show slow decay, and long-term stationarity or growth for the others. We ascribe this behaviour to a combination of complex amphiphile dynamics leading to Marangoni flows.Comment: 16 pages, 13 figures. Accepted for publication in Phys. Rev. E. (Replaced for the latest version, in press.) Higher-quality figures can be sent upon reques

Stress response and structural transitions in sheared gyroidal and lamellar amphiphilic mesophases: lattice-Boltzmann simulations

We report on the stress response of gyroidal and lamellar amphiphilic mesophases to steady shear simulated using a bottom-up lattice-Boltzmann model for amphiphilic fluids and sliding periodic (Lees-Edwards) boundary conditions. We study the gyroid per se (above the sponge-gyroid transition, of high crystallinity) and the molten gyroid (within such a transition, of shorter-range order). We find that both mesophases exhibit shear-thinning, more pronounced and at lower strain rates for the molten gyroid. At late times after the onset of shear, the skeleton of the crystalline gyroid becomes a structure of interconnected irregular tubes and toroidal rings, mostly oriented along the velocity ramp imposed by the shear, in contradistinction with free-energy Langevin-diffusion studies which yield a much simpler structure of disentangled tubes. We also compare the shear stress and deformation of lamellar mesophases with and without amphiphile when subjected to the same shear flow applied normal to the lamellae. We find that the presence of amphiphile allows (a) the shear stress at late times to be higher than in the case without amphiphile, and (b) the formation of rich patterns on the sheared interface, characterised by alternating regions of high and low curvature.Comment: 15 pages, 10 figures, Physical Review E, in pres

Three-dimensional lattice-Boltzmann simulations of critical spinodal decomposition in binary immiscible fluids

We use a modified Shan-Chen, noiseless lattice-BGK model for binary immiscible, incompressible, athermal fluids in three dimensions to simulate the coarsening of domains following a deep quench below the spinodal point from a symmetric and homogeneous mixture into a two-phase configuration. We find the average domain size growing with time as $t^\gamma$, where $\gamma$ increases in the range $0.545 < \gamma < 0.717$, consistent with a crossover between diffusive $t^{1/3}$ and hydrodynamic viscous, $t^{1.0}$, behaviour. We find good collapse onto a single scaling function, yet the domain growth exponents differ from others' works' for similar values of the unique characteristic length and time that can be constructed out of the fluid's parameters. This rebuts claims of universality for the dynamical scaling hypothesis. At early times, we also find a crossover from $q^2$ to $q^4$ in the scaled structure function, which disappears when the dynamical scaling reasonably improves at later times. This excludes noise as the cause for a $q^2$ behaviour, as proposed by others. We also observe exponential temporal growth of the structure function during the initial stages of the dynamics and for wavenumbers less than a threshold value.Comment: 45 pages, 18 figures. Accepted for publication in Physical Review

Numerical simulations of complex fluid-fluid interface dynamics

Interfaces between two fluids are ubiquitous and of special importance for industrial applications, e.g., stabilisation of emulsions. The dynamics of fluid-fluid interfaces is difficult to study because these interfaces are usually deformable and their shapes are not known a priori. Since experiments do not provide access to all observables of interest, computer simulations pose attractive alternatives to gain insight into the physics of interfaces. In the present article, we restrict ourselves to systems with dimensions comparable to the lateral interface extensions. We provide a critical discussion of three numerical schemes coupled to the lattice Boltzmann method as a solver for the hydrodynamics of the problem: (a) the immersed boundary method for the simulation of vesicles and capsules, the Shan-Chen pseudopotential approach for multi-component fluids in combination with (b) an additional advection-diffusion component for surfactant modelling and (c) a molecular dynamics algorithm for the simulation of nanoparticles acting as emulsifiers.Comment: 24 pages, 12 figure

Optimization of chaotic micromixers using finite time Lyapunov exponents

In microfluidics mixing of different fluids is a highly non-trivial task due to the absence of turbulence. The dominant process allowing mixing at low Reynolds number is therefore diffusion, thus rendering mixing in plain channels very inefficient. Recently, passive chaotic micromixers such as the staggered herringbone mixer were developed, allowing efficient mixing of fluids by repeated stretching and folding of the fluid interfaces. The optimization of the geometrical parameters of such mixer devices is often performed by time consuming and expensive trial and error experiments. We demonstrate that the application of the lattice Boltzmann method to fluid flow in highly complex mixer geometries together with standard techniques from statistical physics and dynamical systems theory can lead to a highly efficient way to optimize micromixer geometries. The strategy applies massively parallel fluid flow simulations inside a mixer, where massless and noninteracting tracer particles are introduced. By following their trajectories we can calculate finite time Lyapunov exponents in order to quantify the degree of chaotic advection inside the mixer. The current report provides a review of our results published in [1] together with additional details on the simulation methodology

Acute quadriplegia caused by necrotizing myopathy in a renal transplant recipient with severe pneumonia: acute onset and complete recovery

Critical illness polyneuropathy and myopathy are multifaceted complications that follow severe illnesses involving the sensorimotor axons and proximal skeletal muscles. These syndromes have rarely been reported among renal transplant recipients. In this paper, we report a case of acute quadriplegia caused by necrotizing myopathy in a renal transplant recipient with severe pneumonia. The muscle strength in the patient’s extremities improved gradually after four weeks of comprehensive treatment, and his daily life activities were normal a year after being discharged

A Markov computer simulation model of the economics of neuromuscular blockade in patients with acute respiratory distress syndrome

BACKGROUND: Management of acute respiratory distress syndrome (ARDS) in the intensive care unit (ICU) is clinically challenging and costly. Neuromuscular blocking agents may facilitate mechanical ventilation and improve oxygenation, but may result in prolonged recovery of neuromuscular function and acute quadriplegic myopathy syndrome (AQMS). The goal of this study was to address a hypothetical question via computer modeling: Would a reduction in intubation time of 6 hours and/or a reduction in the incidence of AQMS from 25% to 21%, provide enough benefit to justify a drug with an additional expenditure of $267 (the difference in acquisition cost between a generic and brand name neuromuscular blocker)? METHODS: The base case was a 55 year-old man in the ICU with ARDS who receives neuromuscular blockade for 3.5 days. A Markov model was designed with hypothetical patients in 1 of 6 mutually exclusive health states: ICU-intubated, ICU-extubated, hospital ward, long-term care, home, or death, over a period of 6 months. The net monetary benefit was computed. RESULTS: Our computer simulation modeling predicted the mean cost for ARDS patients receiving standard care for 6 months to be$62,238 (5% – 95% percentiles $42,259 –$83,766), with an overall 6-month mortality of 39%. Assuming a ceiling ratio of $35,000, even if a drug (that cost$267 more) hypothetically reduced AQMS from 25% to 21% and decreased intubation time by 6 hours, the net monetary benefit would only equal $137. CONCLUSION: ARDS patients receiving a neuromuscular blocker have a high mortality, and unpredictable outcome, which results in large variability in costs per case. If a patient dies, there is no benefit to any drug that reduces ventilation time or AQMS incidence. A prospective, randomized pharmacoeconomic study of neuromuscular blockers in the ICU to asses AQMS or intubation times is impractical because of the highly variable clinical course of patients with ARDS

Self-assembly of the gyroid cubic mesophase: Lattice-Boltzmann simulations

We present the first simulations of the self-assembly kinetics of the gyroid cubic mesophase using a Boltzmann transport method. No macroscopic parameters are included in the model and three-dimensional hydrodynamics is emergent from the microscopic conservation laws. The self-assembly arises from local inter-particle interactions in an initially homogeneous, phase-segregating binary fluid with dispersed amphiphile. The mixture evolves in discrete time according to the dynamics of a set of coupled Boltzmann-BGK equations on a lattice. We observe a transient microemulsion phase during self-assembly, the structure function peaks and direct-space imaging unequivocally identifying the gyroid at later times. For larger lattices, highly ordered subdomains are separated by grain boundaries. Relaxation towards the ordered equilibrium structure is very slow compared to the diffusive and microemulsion-assembling transients, the structure function oscillating in time due to a combination of Marangoni effects and long-time-scale defect dynamics