Instabilities, Nonequilibrium and Active Behavior in Complex Fluids
by H. Pleiner
(external collaboration with W. Zimmermann, H.R. Brand, B. Huke, M. Lücke, D. Laroze, and D. Svenšek )
The investigation of instabilities of complex fluids under external fields is motivated by several important aspects. First, detection of the onset allows very precise measurements of relevant material parameters, which very often cannot be measured directly. Second, depending on the nature of the complex fluid, a host of different external fields, like electric and magnetic fields, temperature and concentration gradients, flow and mechanical fields, can be used to drive the system out of equilibrium. Third, in complex fluids even first instabilities can show non-trivial features, thus allowing a rather detailed theoretical description of complicated instability types as well as a quantitative comparison with experiments. Advanced topics are concerned with the spatial structures above onset, their defects and dynamics.
Topic A: smectic films under external electric fields are considered. The advantage of using films lies in their superior visualization properties and in the possibility to explore new geometries in the experimental setup. Generally, one can get instabilities within the layers (e.g. of nematic-like in-plane structures) as well as those of the layers (undulation or buckling).
Topic B is related to shear induced reorientation of layered systems, like amphiphilic L_{α} phases, ordered block copolymers and smectic liquid crystals. We describe the onset of reorientation by a shear induced Helfrich-Hurault type instability using generalized hydrodynamic equations suitable for layered systems. Shear flow (in the x-direction, as shown below) tilts the director ˆn somewhat away from the layer normal ˆp, thus reducing the natural layer spacing. If the layer spacing is kept constant due to the experimental setup, this situation is equivalent to a dilation of the layers. Above a certain threshold undulational instability starts in the plane perpendicular to the shear plane, which eventually leads to a complete reorientation of the layers.
Topic C deals with (mostly thermal) instabilities in ferrofluids. Ferrofluids are dispersions of heavy solid ferromagnetic grains suspended in a carrier liquid. With a typical diameter of 10 nm the particles are pretty large on molecular length scales, resulting in an extremely small particle mobility (very small Lewis number L). This leads to a situation where de-mixing effects take place on time scales far beyond any reasonable observation time. On the other hand ferrofluids are characterized by a large thermo-diffusive or Soret effect (large separation ratio ψ). As a result of the joint action of thermal and solutal buoyancy forces (ψ>0) the critical Rayleigh number Ra_c for the onset of convection is dramatically reduced as compared to the pure-fluid reference value Ra^{0}_c (for ψ=0). However, Ra_c is experimentally inaccessible due to the extremely slow growth of convection patterns, requiring extremely large observation times. For Rayleigh numbers of the order of Ra^{0}_c (ε=0 in the Fig. below) an initially motionless configuration with a uniform concentration distribution grows to a finite amplitude convection pattern. This happens within a time, small compared to the creeping solutal diffusion time, but almost as fast as pure-fluid convection does above Ra^{0}_c [A].
The influence of an external magnetic field on this bifurcation scenario is of experimental interest, since the field is an additional control parameter that is easily manipulated. The full system of hydrodynamic equations for magnetic binary mixtures is derived. Numerical calculation show that only two effects are important, Kelvin force and magnetophoresis. Another striking feature is the enhanced significance of the boundary layer due to the magnetic boundary condition. An approximate analytical solution of the boundary layer problem for concentration and magnetic potential is given. The bifurcation scenario (the saturation amplitude A_{sat} as a function of the external field) is shown below, where M_{1} is proportional to the field intensity reflecting the Kelvin force. The dashed and solid curves are the analytical and numerical results, respectively; ψ_{m} is the magnetic separation ratio based on the concentration dependence of the magnetic susceptibility [B].
The case when the separation ratio (the Soret coefficient) is negative reveals a much richer variety of phenomena than that of positive separation ratio. In particular, for heating from below we find a linear oscillatory instability, whose amplitude, however, relaxes to zero on the long turn and is, thus transient only and, at higher Rayleigh numbers, a finite amplitude stationary instability coexistent with the linearly stable convection-free state (cf. below). By heating from above short-length-scale convective structures occur, whose wavelength depends on the Rayleigh number [C].
We investigate the bifurcation properties and pattern selection in Rayleigh-Bénard convecting binary mixtures for parameter combinations that are typical for colloidal fluids. We study stationary rolls, square and crossroll patterns for positive separation ratios. For colloids with positive separation ratios we found essentially the same sequence of stable structures as for molecular liquids with squares being stable at small reduced Rayleigh number r, rolls being stable at large r, and oscillatory and stationary crossroll patterns in between. For small Lewis number L the intervall where crossrolls exist is very small as long as ψ is smaller than 30 (cf. below). For negative separation ratios and Lewis numbers typical for colloidal solutions neither stable roll nor travelling wave convection is possible. [D].
We describe a convective instability of a stratified ferrofluid subject to a magnetic field: In earth's gravity field, the density difference between particles and solvent is responsible for a sedimentation effect. Waiting long enough a concentration gradient evolves. If a strong enough external magnetic field is applied along the concentration gradient, instability sets in. The ferrofluid is described as a binary mixture featuring two important dimensionless numbers, B the barometric and B_m the magnetic barometric one. At the field strength, where the latter exceeds the former, instability sets in. We discuss [E] the growth rate of linear perturbations and the amplitude of nonlinear stationary convection rolls. Below, the amplitude of the velocity is shown as a function of time with B=100, and B_m= 150, 200, 250, 300 (curves 1-4, respectively). The inset magnifies the long time behavior. The latter is governed by the balance of concentration advection and sedimentation.
For a discussion of instabilities in ferronematics (thermal convection and Saffman-Taylor) and ferrogels and ferrofluids (Rosensweig), cf. also Unconventional Ferrofluids
Topic F: Active Systems:
Active systems are (predominantly) biological systems that are driven out of equilibrium, not by externally applied driving fields, but internally - typically by chemical reactions. As long as the internal driving forces are operating, the units of the active system move, and due to this motion an ordered spatial structure is created. Examples are schools of fish or flocks of birds, pattern forming growing bacteria, biological motors (myosin and actin), or general suspensions of active particles. If the internal drive stops, e.g. due to lack of nutrients, the system becomes passive, with no motion and no remaining localized collective spatial structures, anymore. Orientational order can be present through spontaneous breaking of a continuous symmetry in equilibrium systems, like in nematic and polar nematic liquid crystals. In these equilibrium systems, the spatial structure is static in nature. Similarly, static orientational order can be induced in some systems by strong external (magnetic or electric) fields or mechanical stresses. In both cases the spatial structure and its dynamics is described by a second rank orientational tensor (which can be used to define a director) or a polar vector. These quantities are static in the sense that they do not change under time reversal. It us rather popular to describe also active systems by those static quantities (and adding a so-called "active" term). However, we take a different point of view and describe the active structure as a dynamic one, created by the motion of active units. The natural candidate for the variable to describe the dynamic structure is then (in the polar case [F] ) the velocity and (in the non-polar case [G] ) the mean angular velocity of the active entities. Since there is generally also a non-active background motion (flowing water or solvent, windy air, etc.), there are two different velocities in the polar case. As a result, the hydrodynamics of these (polar) dynamically ordered systems shows, e.g., a normal mode structure quite different compared to the static descriptions, in particular an anisotropic first sound (impossible for equilibrium systems) and a propagating second sound (in equilibrium systems only possible for superfluids and for elastic media). In addition, there are linear couplings between (active) flow and e.g. densities and concentrations due to the genuine two-fluid transport derivatives. The very non-equilibrium nature of these systems is manifest by the existence of a non-zero (active) velocity even in the stationary limit. If this stationary active velocity goes to zero, the passive (equilibrium) state is reached - a transition that is described by well-known bifurcation scenarios. A simplified model on the basis of this dynamic order description has recently be used to successfully explain the stop-and-go dynamics of the growth front of Proteus mirabilis bacteria colonies [H]. Fig.1 shows stroboscopically the time evolution of the growth front, while Fig. 2 demonstrates the stop-and-go behavior.
Fig.1:
Fig.2:
Recent Publications:
Topic A (Smectic Instabilities):
W. Zimmermann, S. Ried, H. Pleiner, and H.R. Brand, "Polarization Frederiks-Transition in SmC*
Liquid Crystal Films", Europhys. Lett. 33, 521 (1996), [839 kB pdf-file] and
Proceedings Freiburger Arbeitstagung Flüssigkristalle 25, 30-3 (1996), [194 kB pdf-file]
S. Ried, H. Pleiner, W. Zimmermann, and H.R. Brand, "Electroconvective Instabilities in SmC*
Liquid Crystal Films", Phys. Rev. E 53, 6101 (1996), [1.17 MB pdf-file]
H. Pleiner and S. Ried, "Instabilities in smectic films", in Nonlinear Physics of Complex Systems
-- Current Status and Future Trends, eds. J. Parisi, S.C. Müller, and W. Zimmermann, Springer Berlin,
p.73 (1996), [870 kB pdf-file]
H. Pleiner, R. Stannarius, and W. Zimmermann, "Electrically driven instabilities in smectic liquid crystal films",
in Evolution of Spontaneous Structures in Dissipative Continuous Systems , eds. F. Busse and S.C. Müller,
Springer Berlin, p. 295 (1998), [782 kB pdf-file]
S. Ried, H. Pleiner and W. Zimmermann, "Electrically driven hybrid instabilities in smectic C liquid crystal films",
Europhys. Lett. 52, 525 (2000), [981 kB pdf-file]
Topic B (Layer Reorientation):
G.K. Auernhammer, H.R. Brand, and H. Pleiner, "The Undulation Instability in Layered Systems under Shear Flow - A Simple Model",
Rheolog. Acta 39, 215 (2000), [865 kB pdf-file] and
"Destabilization of a Layered System by Shear Flow",
Proceedings Freiburger Arbeitstagung Flüssigkristalle 29, 19 (2000), [178 kB pdf-file]
G.K. Auernhammer, H.R. Brand, and H. Pleiner, "Shear induced instabilities in layered liquids",
Phys. Rev. E 66, 061707 (2002), [513 kB pdf-file] DOI: 10.1103/PhysRevE.66.061707
Erratum: Phys. Rev. E 71, 049901(E) (2005) [282 kB pdf-file] DOI: 10.1103/PhysRevE.71.049901
For a comparison with simulation results cf.
T. Soddemann, G.K. Auernhammer, H.X. Guo, B. Dünweg and K. Kremer,
"Shear-induced undulation of smectic-A: Molecular dynamics simulations vs analytical theory",
Eur. Phys. J. E 13, 141 (2004) [417 kB pdf-file] DOI: 10.1140epje/e2004-00045-0
Topic C (Convective Ferrofluid Instabilities):
A. Ryskin, H.-W. Müller, and H. Pleiner, "Thermal convection in binary fluid mixtures with a weak concentration diffusivity but strong solutal buoyancy forces",
Phys. Rev. E 67, 046302 (2003), [150 kB pdf-file] DOI: 10.1103/PhysRevE.67.046302 [A] and
"Thermodiffusive effects in convection of ferrofluids",
Magnetohydrodynamics, 39, 51 (2003) [90 kB pdf-file]
A. Ryskin and H. Pleiner, "The influence of a magnetic field on the Soret-dominated thermal convection in ferrofluids",
Phys. Rev. E 69, 046301 (2004) [210 kB pdf-file] DOI: 10.1103/PhysRevE.69.046301 [B]
A. Ryskin and H. Pleiner, "Thermal convection in colloidal suspensions with negative separation ratio",
Phys. Rev. E 71, 056303 (2005) [432 kB pdf-file] DOI: 10.1103/PhysRevE.71.056303 [C]
B. Huke, H. Pleiner, and M. Lücke, "Convection patterns in colloidal solutions",
Phys. Rev. E 75, 036203 (2007) [694 kB pdf-file] DOI: 10.1103/PhysRevE.75.036203 [D]
A. Ryskin and H. Pleiner, "Magnetic Field Driven Instability in Stratified Ferrofluids",
Phys. Rev. E 75, 056303 (2007) [629 kB pdf-file] DOI: 10.1103/PhysRevE.75.056303 [E]
S. Bohlius, H. Pleiner, and H.R. Brand, "Solution of the Adjoint Problem for Instabilities with a Deformable Surface: Rosensweig and Marangoni Instability",
Phys. Fluids 19, 094103 (2007) [241 kB pdf-file] DOI: 10.1063/1.2757709
B. Huke, H. Pleiner, and M. Lücke, "Thermal convection in mixtures with an inversion of the separation ratio",
Phys. Rev. E 78, 046315 (2008) [5.2 MB pdf-file] or the smaller version [786 kB pdf-file] DOI: 10.1103/PhysRevE.78.046315
A. Ryskin and H. Pleiner, "Influence of sedimentation on convective instabilities in colloidal suspensions",
Int. J. Bif. Chaos 20, 225-34 (2010) [897 kB pdf-file] DOI: 10.1142/S0218127410025661; also at arXiv:0905.0825v1 [nlin.PS]
L.M. Pérez, J. Bragard, D. Laroze, J. Martinez-Mardones, and H. Pleiner, "Thermal convection thresholds in a Oldroyd magnetic fluid",
J. Magn. Magn. Mater. 323, 691 (2011) [856 kB pdf-file] DOI:10.1016/j.jmmm.2010.10.022
D. Laroze, J. Martinez-Mardones, and H. Pleiner, "Bénard-Marangoni instability in a linear viscoelastic magnetic fluid",
Ann. Trans. Nordic Rheol. Soc. 20, 223 (2012) [553 kB pdf-file]
and
"Bénard-Marangoni instability in a viscoelastic ferrofluid",
Eur. Phys. J. ST 219, 71 (2013) [670 kB pdf-file] DOI: 10.1140/epjst/e2013-01782-6
Y. Hayase, M. Kobayashi, D. Vollmer, H. Pleiner, and G.K. Auernhammer, "Asymmetric oscillations during phase separation under continuous cooling - A simple model",
J. Chem. Phys. 129, 184109 (2008) [748 kB pdf-file] DOI: 10.1063/1.3009867
Topic E (Chaos and Localized States):
J. Bragard, H. Pleiner, O.J. Suarez, P. Vargas, J.A.C. Gallas, and D. Laroze, "Chaotic dynamics of a magnetic nanoparticle",
Phys. Rev. E 84, 037202 (2011) [8.1 MB pdf-file] DOI: 10.1103/PhysRevE.84.037202
D. Laroze, J. Bragard, O.J. Suarez, and H. Pleiner, "Characterization of the Chaotic Magnetic Particle Dynamics",
IEEE Trans. Mag. 47, 3032 (2011) [2.5 MB pdf-file] DOI: 10.1109/TMAG.2011.2158072
D. Urzagasti, D. Laroze, M.G. Clerc, S. Coulibaly, and H. Pleiner, "Two-soliton precession state in a parametrically driven magnetic wire",
J. Appl. Phys. 111, 07D111 (2012) [3.1 MB pdf-file] DOI: 10.1063/1.3672872
D. Laroze, D. Becerra-Alonso, J.A.C. Gallas, and H. Pleiner, "Magnetization dynamics under a quasiperiodic magnetic field",
IEEE Trans. Mag. 48, 3567 (2012) [390 kB pdf-file] DOI: 10.1109/TMAG.2012.2207378
D. Laroze, P. G. Siddheshwar, and H. Pleiner, "Chaotic convection in a ferrofluid",
Commun. Nonlinear Sci. Numer. Simulat. 18, 2436-47 (2013) [14 MB pdf-file] DOI: 10.1016/j.cnsns.2013.01.016
D. Urzagasti, D. Laroze, M. G. Clerc, and H. Pleiner, "Breather Soliton Solutions in a Parametrically Driven Magnetic Wire",
Europhys. Lett. 104, 40001 (2013) [6 MB pdf-file] DOI: 10.1209/0295-5075/104/40001
D. Urzagasti, D. Laroze, and H. Pleiner, "Localized chaotic patterns in weakly dissipative systems",
Eur. Phys. J. ST 223, 141-54 (2014) [4.5 MB pdf-file] DOI: 10.1140/epjst/e2014-02089-x
D. Laroze, and H. Pleiner, "Thermal convection in a nonlinear non-Newtonian magnetic fluid",
Ann. Trans. Nordic Rheol. Soc. 21, 35-42 (2013) [11.8 MB pdf-file]
and
Commun. Nonlinear Sci. Numer. Simulat. 26, 167 (2015) [4 MB pdf-file] DOI: 10.1016/j.cnsns.2015.01.002
Topic F (Active Systems):
H.R. Brand, H. Pleiner, and D. Svenšek, "Macroscopic behavior of systems with an axial dynamic preferred direction",
Eur. Phys. J. E 34, 128 (2011), [262 kB pdf-file] DOI: 10.1140/epje/i2011-11128-2 [G]
D. Svenšek, H. Pleiner, and H.R. Brand, "Collective stop-and-go dynamics of active bacteria swarms",
Phys. Rev. Lett. 111, 228101 (2013) [8.1 MB pdf-file] DOI: 10.1103/PhysRevLett.111.228101 [H]
H. Pleiner, D. Svenšek, and H.R. Brand, "Active Polar Two-Fluid Macroscopic Dynamics",
Eur. Phys. J. E 36, 135 (2013) [307 kB pdf-file] DOI: 10.1140/epje/i2013-13135-7[F]
H.R. Brand, H. Pleiner, and D. Svenšek, "Reversible and dissipative macroscopic contributions to the stress tensor: Active or passive?",
Eur. Phys. J. E 37, 83 (2014), [254 kB pdf-file] DOI: 10.1140/epje/i2014-14083-4 [I]
H. Pleiner, D. Svenšek, and H.R. Brand, "Hydrodynamics of Active Polar Systems in a (Visco)Elastic Background"
Rheologica Acta 55, 857 (2016), [376 kB pdf-file] DOI: 10.1007/s00397-016-0957-0 [I]
preprints of the theory group
Other research topics:
Banana and Tetrahedral Phases
Rheology of Complex Fluids
Liquid Crystalline Polymers and Elastomers
Membranes, Films and Surface Waves
General Mesophases
Unconventional Ferrofluids
Back to:
HARALD PLEINER
Last modified: November 12^{th}, 2013