Lattice Boltzmann method investigation of a reactive electro-kinetic flow in porous media: towards a phenomenological model.

A model based on the Lattice Boltzmann method is developed to study the flow of reactive electro-kinetic fluids in porous media. The momentum, concentration and electric/potential fields are simulated via the Navier-Stokes, advection-diffusion/Nernst-Planck and Poisson equations, respectively. With this model, the total density and velocity fields, the concentration of reactants and reaction products, including neutral and ionized species, the electric potential and the interaction forces between the fields can be studied, and thus we provide an insight into the interplay between chemistry, flow and the geometry of the porous medium. The results show that the conversion efficiency of the reaction can be strongly influenced by the fluid velocity, reactant concentration and by porosity of the porous medium. The fluid velocity determines how long the reactants stay in the reaction areas, the reactant concentration controls the amount of the reaction material and with different dielectric constant, the porous medium can distort the electric field differently. All these factors make the reaction conversion efficiency display a non-trivial and non-monotonic behaviour as a function of the flow and reaction parameters. To better illustrate the dependence of the reaction conversion efficiency on the control parameters, based on the input from a number of numerical investigations, we developed a phenomenological model of the reactor. This model is capable of capturing the main features of the causal relationship between the performance of the reactor and the main test parameters. Using this model, one could optimize the choice of reaction and flow parameters in order to improve the performance of the reactor and achieve higher production rates. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

Competition between Ekman Plumes and Vortex Condensates in Rapidly Rotating Thermal Convection.

We perform direct numerical simulations of rotating Rayleigh-Bénard convection (RRBC) of fluids with low (Pr=0.1) and high (Pr≈5) Prandtl numbers in a horizontally periodic layer with no-slip bottom and top boundaries. No-slip boundaries are known to actively promote the formation of plumelike vertical disturbances, through so-called Ekman pumping, that control the ambient flow at sufficiently high rotation rates. At both Prandtl numbers, we demonstrate the presence of competing large-scale vortices (LSVs) in the bulk. Strong buoyant forcing and rotation foster the quasi-two-dimensional turbulent state of the flow that leads to the upscale transfer of kinetic energy that forms the domain-filling LSV condensate. The Ekman plumes from the boundary layers are sheared apart by the large-scale flow, yet we find that their energy feeds the upscale transfer. Our results of RRBC simulations substantiate the emergence of large-scale flows in nature regardless of the specific details of the boundary conditions.

Statistical properties of thermally expandable particles in soft-turbulence Rayleigh-Bénard convection.

The dynamics of inertial particles in Rayleigh-Bénard convection, where both particles and fluid exhibit thermal expansion, is studied using direct numerical simulations (DNS) in the soft-turbulence regime. We consider the effect of particles with a thermal expansion coefficient larger than that of the fluid, causing particles to become lighter than the fluid near the hot bottom plate and heavier than the fluid near the cold top plate. Because of the opposite directions of the net Archimedes' force on particles and fluid, particles deposited at the plate now experience a relative force towards the bulk. The characteristic time for this motion towards the bulk to happen, quantified as the time particles spend inside the thermal boundary layers (BLs) at the plates, is shown to depend on the thermal response time, [Formula: see text], and the thermal expansion coefficient of particles relative to that of the fluid, [Formula: see text]. In particular, the residence time is constant for small thermal response times, [Formula: see text], and increasing with [Formula: see text] for larger thermal response times, [Formula: see text]. Also, the thermal BL residence time is increasing with decreasing K. A one-dimensional (1D) model is developed, where particles experience thermal inertia and their motion is purely dependent on the buoyancy force. Although the values do not match one-to-one, this highly simplified 1D model does predict a regime of a constant thermal BL residence time for smaller thermal response times and a regime of increasing residence time with [Formula: see text] for larger response times, thus explaining the trends in the DNS data well.

Effective Propulsion in Swimming: Grasping the Hydrodynamics of Hand and Arm Movements.

In this paper, a literature review is presented regarding the hydrodynamic effects of different hand and arm movements during swimming with the aim to identify lacunae in current methods and knowledge, and to distil practical guidelines for coaches and swimmers seeking to increase swimming speed. Experimental and numerical studies are discussed, examining the effects of hand orientation, thumb position, finger spread, sculling movements, and hand accelerations during swimming, as well as unsteady properties of vortices due to changes in hand orientation. Collectively, the findings indicate that swimming speed may be increased by avoiding excessive sculling movements and by spreading the fingers slightly. In addition, it appears that accelerating the hands rather than moving them at constant speed may be beneficial, and that (in front crawl swimming) the thumb should be abducted during entry, catch, and upsweep, and adducted during the pull phase. Further experimental and numerical research is required to confirm these suggestions and to elucidate their hydrodynamic underpinnings and identify optimal propulsion techniques. To this end, it is necessary that the dynamical motion and resulting unsteady effects are accounted for, and that flow visualization techniques, force measurements, and simulations are combined in studying those effects.

Spinodal decomposition in homogeneous and isotropic turbulence.

We study the competition between domain coarsening in a symmetric binary mixture below critical temperature and turbulent fluctuations. We find that the coarsening process is arrested in the presence of turbulence. The physics of the process shares remarkable similarities with the behavior of diluted turbulent emulsions and the arrest length scale can be estimated with an argument similar to the one proposed by Kolmogorov and Hinze for the maximal stability diameter of droplets in turbulence. Although, in the absence of flow, the microscopic diffusion constant is negative, turbulence does effectively arrest the inverse cascade of concentration fluctuations by making the low wavelength diffusion constant positive for scales above the Hinze length.

Reynolds number scaling and energy spectra in geostrophic convection.

We report flow measurements in rotating Rayleigh-Bénard convection in the rotationally-constrained geostrophic regime. We apply stereoscopic particle image velocimetry to measure the three components of velocity in a horizontal cross-section of a water-filled cylindrical convection vessel. At a constant, small Ekman number Ek = 5 × 10-8 we vary the Rayleigh number Ra between 1011 and 4 × 1012 to cover various subregimes observed in geostrophic convection. We also include one nonrotating experiment. The scaling of the velocity fluctuations (expressed as the Reynolds number Re) is compared to theoretical relations expressing balances of viscous-Archimedean-Coriolis (VAC) and Coriolis-inertial-Archimedean (CIA) forces. Based on our results we cannot decide which balance is most applicable here; both scaling relations match equally well. A comparison of the current data with several other literature datasets indicates a convergence towards diffusion-free scaling of velocity as Ek decreases. However, the use of confined domains leads at lower Ra to prominent convection in the wall mode near the sidewall. Kinetic energy spectra point at an overall flow organisation into a quadrupolar vortex filling the cross-section. This quadrupolar vortex is a quasi-two-dimensional feature; it only manifests in energy spectra based on the horizontal velocity components. At larger Ra the spectra reveal the development of a scaling range with exponent close to -5/3, the classical exponent for inertial-range scaling in three-dimensional turbulence. The steeper Re(Ra) scaling at low Ek and development of a scaling range in the energy spectra are distinct indicators that a fully developed, diffusion-free turbulent bulk flow state is approached, sketching clear perspectives for further investigation.

Lagrangian acceleration of passive tracers in statistically steady rotating turbulence.

The statistical properties of the Lagrangian acceleration vector of passive tracers in statistically steady rotating turbulence is studied by particle tracking velocimetry. Direct effects of the background rotation are the suppression of high-acceleration events parallel to the (vertical) rotation axis, the enhancement of high-acceleration events for the horizontal acceleration, and the strong amplification of the autocorrelation of the acceleration component perpendicular to both the rotation vector Ω and local velocity vector u. The autocorrelation of the acceleration component in the plane set up by Ω and u is only mildly enhanced.

Force balance in rapidly rotating Rayleigh-Bénard convection.

The force balance of rotating Rayleigh-Bénard convection regimes is investigated using direct numerical simulation on a laterally periodic domain, vertically bounded by no-slip walls. We provide a comprehensive view of the interplay between governing forces both in the bulk and near the walls. We observe, as in other prior studies, regimes of cells, convective Taylor columns, plumes, large-scale vortices (LSVs) and rotation-affected convection. Regimes of rapidly rotating convection are dominated by geostrophy, the balance between Coriolis and pressure-gradient forces. The higher-order interplay between inertial, viscous and buoyancy forces defines a subdominant balance that distinguishes the geostrophic states. It consists of viscous and buoyancy forces for cells and columns, inertial, viscous and buoyancy forces for plumes, and inertial forces for LSVs. In rotation-affected convection, inertial and pressure-gradient forces constitute the dominant balance; Coriolis, viscous and buoyancy forces form the subdominant balance. Near the walls, in geostrophic regimes, force magnitudes are larger than in the bulk; buoyancy contributes little to the subdominant balance of cells, columns and plumes. Increased force magnitudes denote increased ageostrophy near the walls. Nonetheless, the flow is geostrophic as the bulk. Inertia becomes increasingly more important compared to the bulk, and enters the subdominant balance of columns. As the bulk, the near-wall flow loses rotational constraint in rotation-affected convection. Consequently, kinetic boundary layers deviate from the expected behaviour from linear Ekman boundary layer theory. Our findings elucidate the dynamical balances of rotating thermal convection under realistic top/bottom boundary conditions, relevant to laboratory settings and large-scale natural flows.

Finite-size effects lead to supercritical bifurcations in turbulent rotating Rayleigh-Bénard convection.

In turbulent thermal convection in cylindrical samples with an aspect ratio Γ≡D/L (D is the diameter and L the height), the Nusselt number Nu is enhanced when the sample is rotated about its vertical axis because of the formation of Ekman vortices that extract additional fluid out of thermal boundary layers at the top and bottom. We show from experiments and direct numerical simulations that the enhancement occurs only above a bifurcation point at a critical inverse Rossby number 1/Ro(c), with 1/Ro(c)∝1/Γ. We present a Ginzburg-Landau-like model that explains the existence of a bifurcation at finite 1/Ro(c) as a finite-size effect. The model yields the proportionality between 1/Ro(c) and 1/Γ and is consistent with several other measured or computed system properties.

Modifying oculomotor activity in awake subjects increases the amplitude of eye movements during REM sleep.

The eye movements of human subjects were experimentally modified while they were awake to determine the effect of waking experience on electroculographic activity during rapid eye movement (REM) sleep. After normal eye movements were monitored under controlled conditions, subjects spent 5 days wearing goggles that contained minification lenses and that curtailed vision to a 5 degree field. The amplitude and frequency of eye movements decreased when subjects were awake and increased during REM sleep; sleep stage durations and distributions were unaffected. Values returned to normal in the first 24 hours of recovery.

Directional change of tracer trajectories in rotating Rayleigh-Bénard convection.

The angle of directional change of tracer trajectories in rotating Rayleigh-Bénard convection is studied as a function of the time increment τ between two instants of time along the trajectories, both experimentally and with direct numerical simulations. Our aim is to explore the geometrical characterization of flow structures in turbulent convection in a wide range of timescales and how it is affected by background rotation. We find that probability density functions (PDFs) of the angle of directional change θ(t,τ) show similar behavior as found in homogeneous isotropic turbulence, up to the timescale of the large-scale coherent flow structures. The scaling of the averaged (over particles and time) angle of directional change Θ(τ)=〈|θ(t,τ)|〉 with τ shows a transition from the ballistic regime [Θ(τ)∼τ^{c} with c=1] for τ≲τ_{η}, with τ_{η} the Kolmogorov timescale, to a scaling with smaller exponent c for τ_{η}≲τ≲T_{L}, with T_{L} the Lagrangian integral timescale. This scaling exponent is approximately constant in the weakly rotating regime (Rossby number Ro≳2.5) and is decreasing for increasing rotation rates when Ro≲2.5. We show that this trend in the scaling exponent is related with the large-scale coherent structures in the flow; the large-scale circulation for Ro≳2.5 and vertically aligned vortices emerging from the boundary layers (BLs) near the top and bottom plates and penetrating into the bulk for Ro≲2.5. In the viscous BLs, the PDFs of θ(t,τ) and scaling properties of Θ(τ) are in general different from those measured in the bulk and depend on the type of boundary layer, in particular whether the BL is of Prandtl-Blasius type (Ro≳2.5) or of Ekman type (Ro≲2.5). When it is of Ekman type, a stronger dynamic coupling exists between the BL and the bulk of the flow, resulting in similar scaling exponents in BL and bulk.

Velocity and acceleration statistics in rapidly rotating Rayleigh-Bénard convection.

Background rotation causes different flow structures and heat transfer efficiencies in Rayleigh-Bénard convection (RBC). Three main regimes are known: rotation-unaffected, rotation-affected and rotation-dominated. It has been shown that the transition between rotation-unaffected and rotation-affected regimes is driven by the boundary layers. However, the physics behind the transition between rotation-affected and rotation-dominated regimes are still unresolved. In this study, we employ the experimentally obtained Lagrangian velocity and acceleration statistics of neutrally buoyant immersed particles to study the rotation-affected and rotation-dominated regimes and the transition between them. We have found that the transition to the rotation-dominated regime coincides with three phenomena; suppressed vertical motions, strong penetration of vortical plumes deep into the bulk and reduced interaction of vortical plumes with their surroundings. The first two phenomena are used as confirmations for the available hypotheses on the transition to the rotation-dominated regime while the last phenomenon is a new argument to describe the regime transition. These findings allow us to better understand the rotation-dominated regime and the transition to this regime.

Quantifying Photothermal and Hot Charge Carrier Effects in Plasmon-Driven Nanoparticle Syntheses.

The excitation of localized surface plasmon resonances in Au and Ag colloids can be used to drive the synthesis of complex nanostructures, such as anisotropic prisms, bipyramids, and core@shell nanoparticles. Yet, after two decades of research, it is challenging to paint a complete picture of the mechanisms driving such light-induced chemical transformations. In particular, whereas the injection of hot charge carriers from the metal nanoparticles is usually proposed as the dominant mechanism, the contribution of plasmon-induced heating can often not be neglected. Here, we tackle this uncertainty and quantify the contribution of different activation mechanisms using a temperature-sensitive synthesis of Au@Ag core@shell nanoparticles. We compare the rate of Ag shell growth in the dark at different temperatures with the one under plasmon excitation with varying laser intensities. Our controlled illumination geometry, coupled to numerical modeling of light propagation and heat diffusion in the reaction volume, allows us to quantify both localized and collective heating effects and determine their contribution to the total growth rate of the nanoparticles. We find that nonthermal effects can be dominant, and their relative contribution depends on the fraction of nanoparticle suspension under irradiation. Understanding the mechanism of plasmon-activated chemistry at the surface of metal nanoparticles is of paramount importance for a wide range of applications, from the rational design of novel light-assisted nanoparticle syntheses to the development of plasmonic nanostructures for catalytic and therapeutic purposes.

Immunopathologic studies in relapsing polychondritis.

Serial studies have been performed on three patients with relapsing polychondritis in an attempt to define a potential immunopathologic role for degradation constituents of cartilage in the causation and/or perpetuation of the inflammation observed. Crude proteoglycan preparations derived by disruptive and differential centrifugation techniques from human costal cartilage, intact chondrocytes grown as monolayers, their homogenates and products of synthesis provided antigenic material for investigation. Circulating antibody to such antigens could not be detected by immunodiffusion, hemagglutination, immunofluorescence or complement mediated chondrocyte cytotoxicity as assessed by (51)Cr release. Similarly, radiolabeled incorporation studies attempting to detect de novo synthesis of such antibody by circulating peripheral blood lymphocytes as assessed by radioimmunodiffusion, immune absorption to neuraminidase treated and untreated chondrocytes and immune coprecipitation were negative. Delayed hypersensitivity to cartilage constituents was studied by peripheral lymphocyte transformation employing [(3)H]thymidine incorporation and the release of macrophage aggregation factor. Positive results were obtained which correlated with periods of overt disease activity. Similar results were observed in patients with classical rheumatoid arthritis manifesting destructive articular changes. This study suggests that cartilage antigenic components may facilitate perpetuation of cartilage inflammation by cellular immune mechanisms.

Response of the rheumatoid synovial membrane to exogenous immunization.

The secondary immune response to tetanus toxoid in 14 patients with rheumatoid arthritis (RA) has been studied in suspension cultures of peripheral blood lymphocytes (PBL) and synovial membrane obtained at synovectomy. Sequential cultures of PBL from three normal subjects established the optimal time of antibody response at 5 days. At this time, the antitetanus antibody produced was predominantly IgG, comprising half of this immunoglobulin fraction. Rheumatoid synovium synthesized 5-9 times more IgG than PBL, expressed as per cent of total protein synthesis, but only negligible amounts of tetanus antibody. The same results were observed in synovial cultures following repeated immunization and after the additional intra-articular injection of tetanus antigen. This marked limitation of the synovium to respond to exogenous antigen in spite of its large immunoglobulin production was considered consistent with a prior commitment of the synovial lymphoid infiltrate to other antigen.

Finite-size effects on bacterial population expansion under controlled flow conditions.

The expansion of biological species in natural environments is usually described as the combined effect of individual spatial dispersal and growth. In the case of aquatic ecosystems flow transport can also be extremely relevant as an extra, advection induced, dispersal factor. We designed and assembled a dedicated microfluidic device to control and quantify the expansion of populations of E. coli bacteria under both co-flowing and counter-flowing conditions, measuring the front speed at varying intensity of the imposed flow. At variance with respect to the case of classic advective-reactive-diffusive chemical fronts, we measure that almost irrespective of the counter-flow velocity, the front speed remains finite at a constant positive value. A simple model incorporating growth, dispersion and drift on finite-size hard beads allows to explain this finding as due to a finite volume effect of the bacteria. This indicates that models based on the Fisher-Kolmogorov-Petrovsky-Piscounov equation (FKPP) that ignore the finite size of organisms may be inaccurate to describe the physics of spatial growth dynamics of bacteria.

The effect of finger spreading on drag of the hand in human swimming.

The effect of finger spread on overall drag on a swimmer's hand is relatively small, but could be relevant for elite swimmers. There are many sensitivities in measuring this effect. A comparison between numerical simulations, experiments and theory is urgently required to observe whether the effect is significant. In this study, the beneficial effect of a small finger spread in swimming is confirmed using three different but complementary methods. For the first time numerical simulations and laboratory experiments are conducted on the exact same 3D model of the hand with attached forearm. The virtual version of the hand with forearm was implemented in a numerical code by means of an immersed boundary method and the 3D printed physical version was studied in a wind tunnel experiment. An enhancement of the drag coefficient of 2% and 5% compared to the case with closed fingers was found for the numerical simulation and experiment, respectively. A 5% and 8% favorable effect on the (dimensionless) force moment at an optimal finger spreading of 10° was found, which indicates that the difference is more outspoken in the force moment. Moreover, an analytical model is proposed, using scaling arguments similar to the Betz actuator disk model, to explain the drag coefficient as a function of finger spacing.

Lack of effect of donor-recipient ABO mismatching on outcome following allogeneic hematopoietic stem cell transplantation.

Several recently published studies have suggested that patients who undergo ABO mismatched hematopoietic stem cell transplantation may be at increased risk for relapse, graft-versus-host disease, transplant-related mortality, and/or all-cause mortality. To investigate this issue further, we analyzed potential associations between the donor-recipient ABO mismatch pattern and the above outcome measures among 240 consecutive patients who underwent allogeneic hematopoietic stem cell transplantation at our institution. Our analyses uncovered no significant associations between donor-recipient ABO mismatch pattern and overall survival, event-free survival, transplant-related mortality, incidence of acute graft-versus-host disease (GVHD), or incidence of chronic GVHD. Our data do not support recent assertions that donor-recipient ABO mismatching is a major risk factor for patients undergoing allogeneic transplant, nor do they support recent assertions that ABO matching should be an important consideration in selecting allogeneic hematopoietic stem cell donors.

Transitions in turbulent rotating convection: A Lagrangian perspective.

Using measurements of Lagrangian acceleration in turbulent rotating convection and accompanying direct numerical simulations, we show that the transition between turbulent states reported earlier [e.g., S. Weiss et al., Phys. Rev. Lett. 105, 224501 (2010)PRLTAO0031-900710.1103/PhysRevLett.105.224501] is a boundary-layer transition between the Prandtl-Blasius type (typical of nonrotating convection) and Ekman type.

Clustering of vertically constrained passive particles in homogeneous isotropic turbulence.

We analyze the dynamics of small particles vertically confined, by means of a linear restoring force, to move within a horizontal fluid slab in a three-dimensional (3D) homogeneous isotropic turbulent velocity field. The model that we introduce and study is possibly the simplest description for the dynamics of small aquatic organisms that, due to swimming, active regulation of their buoyancy, or any other mechanism, maintain themselves in a shallow horizontal layer below the free surface of oceans or lakes. By varying the strength of the restoring force, we are able to control the thickness of the fluid slab in which the particles can move. This allows us to analyze the statistical features of the system over a wide range of conditions going from a fully 3D incompressible flow (corresponding to the case of no confinement) to the extremely confined case corresponding to a two-dimensional slice. The background 3D turbulent velocity field is evolved by means of fully resolved direct numerical simulations. Whenever some level of vertical confinement is present, the particle trajectories deviate from that of fluid tracers and the particles experience an effectively compressible velocity field. Here, we have quantified the compressibility, the preferential concentration of the particles, and the correlation dimension by changing the strength of the restoring force. The main result is that there exists a particular value of the force constant, corresponding to a mean slab depth approximately equal to a few times the Kolmogorov length scale η, that maximizes the clustering of the particles.