Viscoelastic theory for nematic interfaces

A complete macroscopic theory for compressible nematic-viscous fluid interfaces is developed and used to characterize the interfacial elastic, viscous, and viscoelastic material properties. The derived expression for the interfacial stress tensor includes elastic and viscous components. Surface gradients of the interfacial elastic stress tensor generates tangential Marangoni forces as well as normal forces. The latter may be present even in planar surfaces, implying that in principle static planar interfaces may accommodate pressure jumps. The asymmetric interfacial viscous stress tensor takes into account the surface nematic ordering and is given in terms of the interfacial rate of deformation and interfacial Jaumann derivative. The material function that describes the anisotropic viscoelasticity is the dynamic interfacial tension, which includes the interfacial tension and dilational viscosities. Viscous dissipation due to interfacial compressibility is described by the anisotropic dilational viscosity, and it is shown to describe the Boussinesq surface fluid appropriate for Newtonian interfaces when the director is homeotropic. Three characteristic interfacial shear viscosities are defined according to whether the surface orientation is along the velocity direction, the velocity gradient, or the unit normal. In the last case the expression reduces to the interfacial shear viscosity of the Boussinesq surface fluid. The theory provides a theoretical framework to study interfacial stability, thin liquid film stability and hydrodynamics, and any other interfacial rheology phenomena.     

Steady-profile fingering flows in Marangoni driven thin films

We present experimental and computational results indicating the existence of finite-amplitude fingering solutions in a flow of a thin film of a viscous fluid driven by thermally induced Marangoni stresses. Using carefully controlled experiments, spatially periodic perturbations to the contact line of an initially uniform thin film flow are shown to lead to the development of steady-profile two-dimensional traveling wave fingers. Using an infrared laser and scanning mirror, we impose thermal perturbations with a known wavelength to an initially uniform advancing fluid front. As the front advances in the experiment, we observe convergence to fingers with the initially prescribed wavelength. Experiments and numerical computations show that this family of solutions arises from a subcritical bifurcation.     

Formation of vortex flows in thin nematic cells

The director field reorientation and the formation of steady vortex flows in doubly and obliquely oriented thin liquid-crystal cells under the temperature gradient have been theoretically described within the nonlinear generalization of the classical Ericksen-Leslie theory for describing liquid crystal hydrodynamics. This generalization allowing the consideration of thermomechanical contributions to the shear stress and to the entropy balance equation has made it possible to describe the formation of a steady two-vortex flow in doubly oriented liquid-crystal cells.     

Chaotic boundaries of nematic polymers in mixed shear and extensional flows

Chaotic orientational dynamics of sheared nematic polymers is documented in laboratory experiments and predicted by Doi-Hess kinetic theory for infinitely thin rods. We address robustness of rheochaos when simple shear is modified by a planar straining flow, and the macromolecules have finite aspect ratio. We predict persistence of sheared chaotic response up to a threshold straining flow strength and minimum aspect ratio, beyond which chaotic behavior is arrested. More intriguing, a straining component can induce chaos from periodic shear responses.     

Simulations of MHD flows with moving interfaces

We report on the numerical simulation of a two-fluid magnetohydrodynamics problem arising in the industrial production of aluminium. The motion of the two non-miscible fluids is modeled through the incompressible Navier–Stokes equations coupled with the Maxwell equations. Stabilized finite elements techniques and an arbitrary Lagrangian–Eulerian formulation (for the motion of the interface separating the two fluids) are used in the numerical simulation. With a view to justifying our strategy, details on the numerical analysis of the problem, with a special emphasis on conservation and stability properties and on the surface tension discretization, as well as results on tests cases are provided. Examples of numerical simulations of the industrial case are eventually presented.     

A front-tracking/ghost-fluid method for fluid interfaces in compressible flows

A front-tracking/ghost-fluid method is introduced for simulations of fluid interfaces in compressible flows. The new method captures fluid interfaces using explicit front-tracking and defines interface conditions with the ghost-fluid method. Several examples of multiphase flow simulations, including a shock–bubble interaction, the Richtmyer–Meshkov instability, the Rayleigh–Taylor instability, the collapse of an air bubble in water and the breakup of a water drop in air, using the Euler or the Navier–Stokes equations, are performed in order to demonstrate the accuracy and capability of the new method. The computational results are compared with experiments and earlier computational studies. The results show that the new method can simulate interface dynamics accurately, including the effect of surface tension. Results for compressible gas–water systems show that the new method can be used for simulations of fluid interface with large density differences.     

Theory of surface excess Miesowicz viscosities of planar nematic liquid crystal-isotropic fluid interfaces

An expression for the surface excess stress tensor for planar compressible interfaces between rod-like nematic liquid crystals and isotropic viscous fluids is derived using the classical surface excess theory formalism, adapted to capture the intrinsic anisotropy of the nematic orientational ordering. A required step in the theory is to find the actual stress tensor in the three-dimensional interfacial region, which is obtained by a decomposition of the kinematic fields (rate of deformation tensor and director Jaumann derivative) into tangential, normal, and mixed components with respect to the interface. The viscosity coefficients appearing in the surface excess stress tensor are expressed in terms of interfacial and bulk viscosities for planar, constant orientation, flows. The expressions are used to define the three fundamental surface excess Miesowicz shear viscosities, in analogy with the three bulk Miesowicz shear viscosities. The ordering in the magnitudes of the surface excess Miesowicz shear viscosities is shown to depend on the magnitude of the surface scalar nematic order parameter relative to that of the adjoining bulk nematic phase. When the surface scalar order parameter is greater than in the bulk, the classical ordering in terms of magnitudes of the three bulk Miesowicz shear viscosities is recovered. On the other hand, when the surface scalar order parameter is smaller than in the bulk, the classical ordering in terms of magnitudes of the three viscosities does not hold, and inequality transitions are predicted as the surface scalar order parameter increases towards the bulk value. Received 5 July 1999 and Received in final form 16 November 1999     

Electric field generated solitons,disclinations and vortical flows in freely suspended nematic 8CB

We present the results of optical studies on the instabilities in a substrate-free nematic 8CB film, subject to an in-plane electric field. The initial director field involves a −1/2 strength disclination loop, separating the central pseudoisotropic zone from the splay-bend (SB) birefringent boundary. Three regimes of sample thickness are distinguishable on the basis of field-induced instabilities. Thick (∼75 μm) films display growth of SB zones (independently of disclination movement), wall-formation, and reversible transition between walls and disclinations. Moderately thick films, a few μm in the central part, exhibit distinctive undulations at the border of advancing SB layers. Submicron thin films, with smectic-like homeotropic central plateau, show spectacular isotropic vortex-pairs at either end of this plateau. Further, the end regions of the birefringent zone exhibit both electro-convective flows and reorientational effects. The latter are associated with the formation of open and closed walls, and loop-wall emission. The final high field instability involves jet-like flows at the two ends of the film.     

Finite element approximation of nematic liquid crystal flows using a saddle-point structure

In this work, we propose finite element schemes for the numerical approximation of nematic liquid crystal flows, based on a saddle-point formulation of the director vector sub-problem. It introduces a Lagrange multiplier that allows to enforce the sphere condition. In this setting, we can consider the limit problem (without penalty) and the penalized problem (using a Ginzburg–Landau penalty function) in a unified way. Further, the resulting schemes have a stable behavior with respect to the value of the penalty parameter, a key difference with respect to the existing schemes. Two different methods have been considered for the time integration. First, we have considered an implicit algorithm that is unconditionally stable and energy preserving. The linearization of the problem at every time step value can be performed using a quasi-Newton method that allows to decouple fluid velocity and director vector computations for every tangent problem. Then, we have designed a linear semi-implicit algorithm (i.e. it does not involve nonlinear iterations) and proved that it is unconditionally stable, verifying a discrete energy inequality. Finally, some numerical simulations are provided.     

Influence of barrier effects at interfaces on dynamic scattering of light in a nematic liquid crystal

The intensity of light dynamic scattering by a nematic liquid crystal with negative dielectric anisotropy (ZhK-440) as a function of the constant electric field strength has been studied under different interfacial conditions. It has been shown that potential barriers that arise at the interfaces not only influence the scattering intensity, but may also radically change the form of the dependence; i.e., the curves of optical transmission in the direction of the incident beam may have a minimum under certain field strengths. At higher strengths, the cells become totally transparent again. This anomalous behavior of the transmission curve is associated with the fact that the conductivity of the cells drops below a critical value with growing field strength, as the resistance of the space charge region at the inversely biased junction of the extraction electrode rises. In addition, it has been shown that the high resistance of the cells at low voltages may be attributed not only to a low concentration of impurities in the liquid crystal, but also to a low emissivity of the injecting electrode and a weak electron affinity of the extraction electrode.     

An immersed interface method for Stokes flows with fixed/moving interfaces and rigid boundaries

We present an immersed interface method for solving the incompressible steady Stokes equations involving fixed/moving interfaces and rigid boundaries (irregular domains). The fixed/moving interfaces and rigid boundaries are represented by a number of Lagrangian control points. In order to enforce the prescribed velocity at the rigid boundaries, singular forces are applied on the fluid at these boundaries. The strength of singular forces at the rigid boundary is determined by solving a small system of equations. For the deformable interfaces, the forces that the interface exerts on the fluid are calculated from the configuration (position) of the deformed interface. The jumps in the pressure and the jumps in the derivatives of both pressure and velocity are related to the forces at the fixed/moving interfaces and rigid boundaries. These forces are interpolated using cubic splines and applied to the fluid through the jump conditions. The positions of the deformable interfaces are updated implicitly using a quasi-Newton method (BFGS) within each time step. In the proposed method, the Stokes equations are discretized via the finite difference method on a staggered Cartesian grid with the incorporation of jump contributions and solved by the conjugate gradient Uzawa-type method. Numerical results demonstrate the accuracy and ability of the proposed method to simulate incompressible Stokes flows with fixed/moving interfaces on irregular domains.     

Stability considerations associated with the meniscoid particle band at advancing interfaces in hele-shaw suspension flows

We present a new convergent strong-coupling expansion for two-level atoms in external periodic fields, free of secular terms. As a first application, we show that the coherent destruction of tunneling is a third-order effect. We also present an exact treatment of the high-frequency region, and compare it with the theory of averaging. The qualitative frequency spectrum of the transition probability amplitude contains an effective Rabi frequency.     

Simple and efficient relaxation methods for interfaces separating compressible fluids,cavitating flows and shocks in multiphase mixtures

Numerical approximation of the five-equation two-phase flow of Kapila et al. [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10) (2001) 3002–3024] is examined. This model has shown excellent capabilities for the numerical resolution of interfaces separating compressible fluids as well as wave propagation in compressible mixtures [A. Murrone, H. Guillard, A five equation reduced model for compressible two phase flow problems, Journal of Computational Physics 202(2) (2005) 664–698; R. Abgrall, V. Perrier, Asymptotic expansion of a multiscale numerical scheme for compressible multiphase flows, SIAM Journal of Multiscale and Modeling and Simulation (5) (2006) 84–115; F. Petitpas, E. Franquet, R. Saurel, O. Le Metayer, A relaxation-projection method for compressible flows. Part II. The artificial heat exchange for multiphase shocks, Journal of Computational Physics 225(2) (2007) 2214–2248]. However, its numerical approximation poses some serious difficulties. Among them, the non-monotonic behavior of the sound speed causes inaccuracies in wave’s transmission across interfaces. Moreover, volume fraction variation across acoustic waves results in difficulties for the Riemann problem resolution, and in particular for the derivation of approximate solvers. Volume fraction positivity in the presence of shocks or strong expansion waves is another issue resulting in lack of robustness. To circumvent these difficulties, the pressure equilibrium assumption is relaxed and a pressure non-equilibrium model is developed. It results in a single velocity, non-conservative hyperbolic model with two energy equations involving relaxation terms. It fulfills the equation of state and energy conservation on both sides of interfaces and guarantees correct transmission of shocks across them. This formulation considerably simplifies numerical resolution. Following a strategy developed previously for another flow model [R. Saurel, R. Abgrall, A multiphase Godunov method for multifluid and multiphase flows, Journal of Computational Physics 150 (1999) 425–467], the hyperbolic part is first solved without relaxation terms with a simple, fast and robust algorithm, valid for unstructured meshes. Second, stiff relaxation terms are solved with a Newton method that also guarantees positivity and robustness. The algorithm and model are compared to exact solutions of the Euler equations as well as solutions of the five-equation model under extreme flow conditions, for interface computation and cavitating flows involving dynamics appearance of interfaces. In order to deal with correct dynamic of shock waves propagating through multiphase mixtures, the artificial heat exchange method of Petitpas et al. [F. Petitpas, E. Franquet, R. Saurel, O. Le Metayer, A relaxation-projection method for compressible flows. Part II. The artificial heat exchange for multiphase shocks, Journal of Computational Physics 225(2) (2007) 2214–2248] is adapted to the present formulation.     

Experimental investigation of Marangoni convection in nanofluids

Onset of regular convection in nanoparticle doped isotropic and anisotropic fluids with one free surface under absorption of light with Gaussian distribution of intensity due to temperature dependence of surface tension coefficient is experimentally studied. It is shown that nanoparticles essentially increase the thermal conductivity of the mixture, which leads to decrease of the velocity of convective motion. Dependence of velocity of hydrodynamic motions on intensity of laser beam and concentration of nanoparticles is also studied.     

Threshold onset of Marangoni convection in narrow channels

There are several experimental studies where the Marangoni convection begins only at a certain difference in the surface tension, i.e. in a threshold way. This effect contradicts a traditional point of view according to which the surface flow in Newtonian fluids should begin at an arbitrary small difference in surface tension. To explore this phenomenon in detail we investigated the initiation of the Marangoni convection at a free liquid surface caused by injection of a droplet of surfactant. It was found that the surface motion occurs in a threshold manner, i.e. when a surfactant concentration in the droplet approaches a certain critical value. The described phenomenon is more important in narrow channels and essentially depends both on the purity of the basic liquid and on the surfactant used. Based on the experimental results, a hypothesis about an important role of residual impurities contained in basic liquids which can thoroughly change a surface rheology was suggested. The theoretical model taking into account special rheological properties in the free surface is considered. The results of the numerical simulation are in a good agreement with the experimental observations.