Shear cell rupture of nematic liquid crystal droplets in viscous fluids

We model the hydrodynamics of a shear cell experiment with an immiscible nematic liquid crystal droplet in a viscous fluid using an energetic variational approach and phase-field methods [86]. The model includes the coupled system for the flow field for each phase, a phase-field function for the diffuse interface and the orientational director field of the liquid crystal phase. An efficient numerical scheme is implemented for the two-dimensional evolution of the shear cell experiment for this initial data. The same model reduces to an immiscible viscous droplet in a viscous fluid, which we simulate first to compare with other numerical and experimental behavior. Then we simulate drop deformation by varying capillary number (independent of liquid crystal physics), liquid crystal interfacial anchoring energy and Oseen–Frank distortional elastic energy. We show the number of eventual droplets (one to several) and “beads on a string” behavior are tunable with these three physical parameters. All stable droplets possess signature quadrupolar shear and normal stress distributions. The liquid crystal droplets always possess a global surface defect structure, called a boojum, when tangential surface anchoring is imposed. Boojums [79], [32] consist of degree +1/2 and ?1/2 surface defects within a bipolar global orientational structure.     

Experimental techniques in natural convection

Experimental techniques in natural convection heat transfer employed in the author's laboratory are introduced. The techniques are mostly related to visualization of flow, temperature field, and heat flux distribution in fluids. Three topics are presented, the first being natural convection in a horizontal rectangular liquid layer driven by surface tension and buoyancy. The patterns of flow were visualized by suspending fine aluminum flakes in the liquid. At the same time, the distribution of the temperature gradient in the liquid was visualized by an optical method making use of the refraction of light. The second topic is the onset of oscillatory convection in the Czochralski growth melt. In this case a forced flow due to rotation of the crystal and the vessel is superimposed on the buoyancy convection, resulting in an oscillatory flow under certain circumstances. The flow pattern and the temperature distribution in the liquid were visualized simultaneously by suspending in the liquid a microencapsulated temperature-sensitive liquid crystal. Periodical oscillation of the flow and the temperature was clearly recognized. The third topic is the rollover of double liquid layers that were stratified stably due to a density difference. A small-scale experiment was carried out to clarify the basic mechanism of rollover. The tracer method was used to visualize boundary layer flow along the vertical side wall and the shadowgraph technique to visualize the density distribution in the liquid layers. The article emphasizes the importance of visual observation in the investigation of natural convection phenomena.     

Optical response of a layer of a nematic liquid crystal to an air flow

A liquid crystal coating is used to measure the surface friction created when a flat plate is subjected to an air flow. Surface friction is determined from the optical response of the nematic liquid crystal coating to the flow. The proposed method does not require precise monitoring of the thickness of the coating or the angles of illumination and observation. This makes it possible to eventually progress to panoramic measurements of surface aerodynamic characteristics. Institute of Theoretical and Applied Mechanics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 4, pp. 102–109, July–August, 1998.     

半浮区液桥热毛细振荡流

采用非定常、三维直接数值模拟方法研究大Pr数半浮区液桥热毛细对流从定常流向振荡流的过渡过程．文中详细描述了热毛细振荡流的起振和振荡特征,给出了液桥横截面上振荡流的流场和温度分布．在地面引力场条件下计算的结果与地面实验的结果进行比较,得出液桥水平截面上的流场和温度分布图样以一定的速度旋转,自由表面固定点处流体的环向流速正、负交替变化的一致结论．     

Surface oil flow technique and liquid crystal thermography for flow visualization in impulse,wind tunnels

This paper describes flow visualization techniques employing surface oil flow and liquid crystal thermography suitable for use in impulse wind tunnels. High spatial resolution photographs of oil flow patterns and liquid crystal thermograms have been obtained within test times ranging from 7 to 500 ms and have been shown to be very useful for revealing the detailed features of 3-D separated flow. The results from oil flow patterns, liquid crystal thermograms, schlieren photographs and heat flux measurements are shown to be in good agreement.     

Detection of flow separation and reattachment using shear-sensitive liquid crystals

Coatings of pure chiral nematic liquid crystals are known to change colour under different levels of surface shear stress. In this study, the liquid crystal was used to provide information about flow separation and reattachment on both a two-dimensional aerofoil and a delta wing. The tests were carried out at a free-stream velocity of 28 m/s and a number of incidence angles. The Reynolds numbers based on the central chord length of the models were 200,000 and 270,000 for the aerofoil and delta wing models, respectively. The study showed that locations of boundary layer separation and reattachment can be identified from spatial variations in the surface colour; the agreement between the results and those obtained using surface oil flow was good. Issues relating to interpretation of the crystal colour pattern and the limitation of this technique in detection of flow separation were also discussed.     

聚变堆相关的液态金属膜流初步实验研究

在磁约束核聚变堆的面对等离子部件设计中,液态金属锂膜流因具有带走杂质、保护面对等离子固壁等优点而被认为是优选方案之一. 然而,如何克服聚变堆中强磁场环境下产生的磁流体力学效应并形成大面积均匀铺展锂膜流动是目前亟需解决的问题.本文通过搭建室温液 态镓铟锡回路和高温液态锂回路,开展了两种不同特性的液态金属膜流实验, 并采用传统可视化方法获得了展向磁场存在时镓铟锡和锂在导电底板形成的液膜流动表面特征.实验结果 表明: 无磁场时,两种液态金属膜流流动表面波动特性与常规流体膜流均一致, 即随着流动雷诺数的增加表面波动变得更为混乱; 而展向磁场存在时,镓铟锡膜流表面波动变得更为规则, 且沿着磁场方向平行排列,表现为拟二维波动的特征; 而锂膜流却产生了明显的磁流体 力学阻力效应,表现为在流动方向局部产生锂滞留现象, 且滞留点随雷诺数增大向下游移动. 最后通过膜流受力分析,进一步阐述了锂膜流受到比镓铟锡膜流更为严重磁流体力学效应影响的原因.      

Abnormal rotation of a deformed liquid crystal droplet immersed in an isotropic fluid after transient flow

The morphology evolution of liquid crystal droplets immersed in an isotropic fluid in flow field is found to be different from flexible polymer droplets. In this paper, we investigated the retraction of a liquid crystal droplet after transient flow. It is found that the liquid crystal droplet will rotate during the shape recovery, which has never been observed for an isotropic droplet. The factors that influence the rotational angle of a single liquid crystal droplet during retraction progress were studied, including the temperature, the dimension of the droplets, the time of shear flow, the shear rate, the flow type, and the properties of liquid crystal molecules. The rotation of liquid crystal droplet during shape recovery is ascribed to both the bulk elasticity of liquid crystal droplets and the anisotropic properties of the interface between liquid crystal and isotropic fluid.     

Flow and temperature field measurements of thermal convection in a small vertical gap using liquid crystals

Thermal convection in a small vertical gap is studied experimentally applying digital particle image velocimetry/thermometry. This optical method enables the simultaneous measurement of two-dimensional flow and temperature fields in a liquid. The principle is based on seeding the liquid flow medium with thermochromic liquid crystal particles. The temperature is measured by the crystal particles which change their reflected colour as function of temperature. The flow velocity is measured by using the same particles as flow tracers. The investigation shall contribute to the understanding of the fluid mechanical behaviour of biological liquids within micro reactor systems. However, the problem is also of fundamental interest as far as heat and mass transfer is concerned. Measured temperature and flow velocity fields are presented and discussed. Presented in part at the 4th Chemnitz/Hamburger Colloquium (CHC) on Microflows, Hamburg, Germany, November 2004.     

Viscous effects in liquid encapsulated liquid bridges

An analytical derivation of the surface deflections and the streamfunctions for the flow inside a liquid encapsulated liquid bridge has been derived using an asymptotic expansion about a small capillary number. The model assumes an initially flat and cylindrical interface under the assumption that the densities of both fluids are equal. To simplify the analysis, the top and bottom walls are assumed to be stress-free and the Reynolds number is assumed to be negligible. Flow is generated either by a moving outer wall (shear-driven flow) or by applying a temperature difference across the top and bottom walls (Marangoni-driven flow). The resulting equations show that for the shear-driven flow, as the viscosity ratio increases, the surface deflections increase monotonically. For the Marangoni-driven flow there exist values of the viscosity ratio where the surface deflections reach a minimum and then switch signs. This investigation shows that it may be possible in more realistic systems to use an outer encapsulating liquid of the proper viscosity ratio to stabilize the liquid–liquid interface during float zone crystal growth.     

Skin friction measurements on reflective surfaces using nematic liquid crystal

 A new liquid crystal technique for full surface skin friction measurements is introduced. With the new technique, the transmission of light through nematic liquid crystal viewed in reflection provides a quantitative measurement of the skin friction. The measurement technique is discussed with the aid of a model which describes the rotation dynamics of the liquid crystal molecules. Calibration experiments performed in a laminar flow duct demonstrate that the model captures the essential physics of the new technique. The measurement of skin friction downstream of a three dimensional roughness element in an incompressible laminar boundary layer is then presented as a demonstration of the utility of the technique. Received: 5 June 1998/Accepted: 13 November 1998     

Equalizing calculus in Trefftz method for solving two-dimensional temperature field of FC-72 flowing along the minichannel

The paper focuses on the numerical solution to two-dimensional temperature field of boiling liquid flowing along a vertical, asymmetrically heated minichannel with a rectangular cross-section. One of the walls of a minichannel is DC supplied single-sided enhanced foil with mini-recesses distributed unevenly in the selected area. The parallel walls are made of glass panes for thermal insulation and they are intended for observation of the two-phase flow and the void fraction. The thin layer of thermosensitive liquid crystal paint on the outer side of the foil enabled to record two-dimensional temperature distribution of outer foil surface. The paper described computations based on Trefftz method for finding two-dimensional temperature field of boiling liquid flowing along the minichannel. The presented research is limited only to the liquid phase of the two-phase mixture observed in the minichannel. The velocity of liquid flowing through the minichannel is represented by a piecewise linear approximating function. To solve energy equation for liquid phase, Trefftz functions specially generated for this purpose were employed. Temperature field in the fluid was approximated by a linear combination of Trefftz functions. Equalizing calculus was applied to the Trefftz method to smooth temperature measurements and reduce measurement errors. Temperature at the interface between working fluid and foil amounts to the saturation temperature. Temperature distribution in the foil and the glass pane was also computed using proper Trefftz functions.     

Dispersed multiphase flow with air-driven runback of a liquid layer at a moving boundary

Multiphase flow with impinging droplets on an icing surface with a flowing supercooled surface layer is investigated. The air-assisted flowing layer is modelled with the cross-phase shear stresses imparted at the moving liquid/air interface. Runback and runoff of the surface layer are predicted by mass flow across the boundaries between adjacent elements in the numerical formulation. This liquid runoff is determined by coupled heat and momentum balances for the unfrozen water layer. The numerical analysis is developed with a control-volume-based finite element method (CVFEM). An Eulerian formulation with volume averaging is developed to accommodate the near-wall elements containing both dispersed and continuous phases. The predicted results are successfully validated through comparisons with analytical solutions and measured data.     

Adhesion of oriented nematic droplets

We propose and analyze a two-dimensional model for the equilibrium of oriented droplets of nematic liquid crystals that may adhere to a rigid substrate, while surrounded by an isotropic environment. We obtain the contact condition at the edge where the liquid crystal, the substrate, and the environment come together. We further develop a fairly general method to arrive at the equilibrium shapes of a drop, which is then applied to the case where the surface tension at the liquid crystal interface is given by Rapini and Papoular's expression. In this case, we also predict the existence of concave equilibrium shapes. Here is indeed the main difference between this method and Wulff's construction, which always yields convex equilibrium shapes for a drop free from adhesion. Received February 22, 2000     

Numerical simulation of 3D free surface flows,with multiple incompressible immiscible phases. Applications to impulse waves

A numerical method for the solution to the density‐dependent incompressible Navier–Stokes equations modeling the flow of N immiscible incompressible liquid phases with a free surface is proposed. It allows to model the flow of an arbitrary number of liquid phases together with an additional vacuum phase separated with a free surface. It is based on a volume‐of‐fluid approach involving N indicator functions (one per phase, identified by its density) that guarantees mass conservation within each phase. An additional indicator function for the whole liquid domain allows to treat boundary conditions at the interface between the liquid domain and a vacuum. The system of partial differential equations is solved by implicit operator splitting at each time step: first, transport equations are solved by a forward characteristics method on a fine Cartesian grid to predict the new location of each liquid phase; second, a generalized Stokes problem with a density‐dependent viscosity is solved with a FEM on a coarser mesh of the liquid domain. A novel algorithm ensuring the maximum principle and limiting the numerical diffusion for the transport of the N phases is validated on benchmark flows. Then, we focus on a novel application and compare the numerical and physical simulations of impulse waves, that is, waves generated at the free surface of a water basin initially at rest after the impact of a denser phase. A particularly useful application in hydraulic engineering is to predict the effects of a landslide‐generated impulse wave in a reservoir. Copyright © 2014 John Wiley & Sons, Ltd.     

Morphological stability analysis of vesicles with mechanical–electrical coupling effects

Using a recently established liquid crystal model for vesicles, we present a theoretical method to analyze the morphological stability of liquid crystal vesicles in an electric field. The coupled mechanical-electrical effects associated with elastic bending, osmotic pressure, surface tension, Max- well pressure, as well as flexoelectric and dielectric proper- ties of the membrane are taken into account. The first and second variations of the free energy are derived in a com- pact form by virtue of the surface variational principle. The former leads to the shape equation of a vesicle embedded in an electric field, and the latter allows us to examine the stabil- ity of a given vesicle morphology. As an illustrative exam- ple, we analyze the stability of a spherical vesicle under a uniform electric field. This study is helpful for understanding and revealing the morphological evolution mechanisms of vesicles in electric fields and some associated phenomena of cells.     

The effect of an insoluble surfactant on the skin friction of a bubble

In this paper, we develop a novel moving mesh method suitable for solving axisymmetric free-boundary problems, including the Marangoni effect induced by surfactant or temperature variation. This method employs a body-fitted grid system where the gas–liquid interface is one line of the grid system. We model the surfactant equation of state with a non-linear Langmuir law, and, for simplicity, we limit ourselves to the situation of an insoluble surfactant. We solve complicated dynamic boundary conditions accurately on the gas–liquid interface in the framework of finite-volume methods. Our method is used to study the effect of a surfactant on the skin friction of a bubble in a uniaxial flow. For the limiting case where the surface diffusivity is zero, the effect of a tangential stress generated by the surface tension gradient, allows us to explain a new phenomenon in high concentration regimes: larger surface tension, but also larger deformation. Furthermore, this condition leads to the formation of boundary layers and flow separation at high Reynolds numbers. The influence of these complex flow patterns is examined.     

On the influence of buoyancy on the surface tension driven flow around a bubble on a heated wall

The surface tension driven flow in the liquid vicinity of gas bubbles on a heated solid wall has been investigated both, in a reduced gravity environment aboard a sounding rocket, and in an earth-bound experiment. Both experiments deal with temperature gradients within the liquid surrounding of a bubble which cause variations of the surface tension. These, in turn, lead to a liquid flow around the bubble periphery termed thermocapillary or thermal Marangoni-convection. On Earth, this phenomenon is widely masked by buoyancy. We therefore carried out an experiment under reduced gravitational acceleration. In order to simultaneously observe and record the flow field and the temperature field liquid crystal tracers have been applied. These particles offer the feature of selectively reflecting certain wavelengths of incident white light depending on the crystals temperature. Although the bubble injection system did not perform nominally during the flight experiment, some interesting flow characteristics could be observed. Comparison of results obtained in microgravity to data measured on Earth reveal that due to the interaction of thermocapillarity and buoyancy a very compact vortex flow results on ground, while in microgravity the influence on the surface tension driven flow penetrates much deeper into the bulk. This result is of special interest regarding the production of materials in space. Dedicated to Professor Dr. Julius Siekmann on the occasion of his 70th birthday The work described herein was supported by the German space agency DARA (Deutsche Agentur für Raumfahrtangelegenheiten GmbH) through DARA Grant 50 WM 9434. The authors thank the European Space Agency (ESA) for the opportunity to conduct the TEXUS 33 sounding rocket experiment. The flight hardware has been partly built by Daimler-Benz-Aerospace which is gratefully acknowledged. Also, the authors are indebted to Mr. H.-H. Wolf for his careful evaluation of the particle images     

On the stability of shear flows in nematics

A theoretical study of the effect of an applied magnetic field on the stability of the flow of nematic slabs subjected to an arbitrary shear is presented. Homeotropic boundary conditions with strong anchoring and a constant magnetic field applied perpendicular to the plates are considered. We discuss the general conditions on the control parameters under which the flow is stable, for a low molecular weight liquid crystal and for a polymer liquid crystal, and obtain estimations of the critical values.     

高超声速边界层液膜演化过程和冷却机理研究

高超声速液膜冷却技术是通过一系列狭缝或孔洞压出冷却工质,在飞行器表面边界层形成一层低温冷却膜,阻止高超声速气流对飞行器的气动加热.其作为一种主动冷却方式在高超声速飞行器表面热防护有着巨大的应用潜力.文章采用数值方法,结合VOF模型,研究25 km飞行高度和Ma=5气流条件下的液膜铺展情况,并通过不同冷却工质的入射速度、角度、表面张力和黏性系数条件,讨论了液膜在平板上的演化过程和冷却机理.结果表明,在气流作用下,液膜向壁面下游发展,液膜的存在导致边界层分离,连续液膜会在一定位置断裂为液块,然后进一步破碎为液滴.入射条件和液体性质的改变,会影响液膜沿流向的发展,具体表现在连续液膜断裂点的位置和连续液膜的厚度.在所设定的计算域内,壁面热流降低了80%～95%,液膜对壁面的冷却效率随着液膜形态的变化而变化.