Comparative study of two lattice Boltzmann multiphase models for simulating wetting phenomena: implementing static contact angles based on the geometric formulation

Wetting phenomena are widespread in nature and industrial applications. In general, systems concerning wetting phenomena are typical multicomponent/multiphase complex fluid systems. Simulating the behavior of such systems is important to both scientific research and practical applications. It is challenging due to the complexity of the phenomena and difficulties in choosing an appropriate numerical method. To provide some detailed guidelines for selecting a suitable multiphase lattice Boltzmann model, two kinds of lattice Boltzmann multiphase models, the modified S-C model and the H-C-Z model, are used in this paper to investigate the static contact angle on solid surfaces with different wettability combined with the geometric formulation (Ding, H. and Spelt, P. D. M. Wetting condition in diffuse interface simulations of contact line motion. Physical Review E, 75(4), 046708 (2007)). The specific characteristics and computational performance of these two lattice Boltzmann method (LBM) multiphase models are analyzed including relationship between surface tension and the control parameters, the achievable range of the static contact angle, the maximum magnitude of the spurious currents (MMSC), and most importantly, the convergence rate of the two models on simulating the static contact angle. The results show that a wide range of static contact angles from wetting to non-wetting can be realized for both models. MMSC mainly depends on the surface tension. With the numerical parameters used in this work, the maximum magnitudes of the spurious currents of the two models are on the same order of magnitude. MMSC of the S-C model is universally larger than that of the H-C-Z model. The convergence rate of the S-C model is much faster than that of the H-C-Z model. The major foci in this work are the frequently-omitted important details in simulating wetting phenomena. Thus, the major findings in this work can provide suggestions for simulating wetting phenomena with LBM multiphase models along with the geometric formulation.     

Simulations of free surface flows with implementation of surface tension and interface sharpening in the two-fluid model

Two-fluid model used for free surface flows with large characteristic scales is improved; the smeared interface is sharpened with conservative level set method and the surface tension force with wetting angle is implemented. Surface tension force is split between two phases with several models. Detailed analysis showed the splitting of surface tension force with volume averaging as the most appropriate. The improved two-fluid model with interface sharpening and implemented surface tension is validated on several test cases. The pressure jump over a droplet interface test case showed that the pressure jump in simulation converges with grid refinement to the analytical one. The parasitic currents in simulation are one order of magnitude larger than in simulation with volume of fluid model. In the oscillating droplet test case the time period of oscillating droplet with initially ellipsoid or square shape is similar to the analytical time period. In the rising bubble test case, the rising bubble position, terminal velocity, and circularity are similar to the one observed in simulations with level set model. The wetting angle is implemented in the two-fluid model with interface sharpening and surface tension force. Model is tested in the simulation of droplet in contact with wall with different wetting angles.     

An approach to control the spurious currents in a multiphase lattice Boltzmann method and to improve the implementation of initial condition

Multiphase lattice Boltzmann methods are known to generate spurious or parasitic currents at the fluid–fluid interfaces. This nonphysical phenomenon has to be avoided, or at least controlled, in order to achieve reliable solutions. In this article, a method to control these fictitious velocities via lattice refinement is proposed, which is based on interface thickness control for which both the spurious currents and the physical fluid–fluid interface thickness vanishes as the spatial resolution increases. It has been found that a proper interface thickness adjustment is required as the lattice refinement is applied, or an increase in spurious currents, instead of a reduction, can occur. By combining the new method with an appropriate multiphase flow initialization, the overall stability for high density O(1000) and viscosity O(100) ratios is greatly improved. Although this research has been conducted with a Rothman and Keller type lattice Boltzmann model, it is believed that other types of multiphase lattice Boltzmann models could benefit from the basic ideas underlying this research. Copyright © 2015 John Wiley & Sons, Ltd.     

Issue Information

An integrated finite element method (FEM) is proposed to simulate incompressible two‐phase flows with surface tension effects, and three different surface tension models are applied to the FEM to investigate spurious currents and temporal stability. A Q2Q1 element is adopted to solve the continuity and Navier–Stokes equations and a Q2‐iso‐Q1 to solve the level set equation. The integrated FEM solves pressure and velocity simultaneously in a strongly coupled manner; the level set function is reinitialized by adopting a direct approach using interfacial geometry information instead of solving a conventional hyperbolic‐type equation. In addition, a consistent continuum surface force (consistent CSF) model is utilized by employing the same basis function for both surface tension and pressure variables to damp out spurious currents and to estimate the accurate pressure distribution. The model is further represented as a semi‐implicit manner to improve temporal stability with an increased time step. In order to verify the accuracy and robustness of the code, the present method is applied to a few benchmark problems of the static bubble and rising bubble with large density and viscosity ratios. The Q2Q1‐integrated FEM coupled with the semi‐implicit consistent CSF demonstrates the significantly reduced spurious currents and improved temporal stability. The numerical results are in good qualitative and quantitative agreements with those of the existing studies. Copyright © 2015 John Wiley & Sons, Ltd.     

The Zero Surface Tension Limit of Two-Dimensional Interfacial Darcy Flow

We perform energy estimates for a sharp-interface model of two-dimensional, two-phase Darcy flow with surface tension. A proof of well-posedness of the initial value problem follows from these estimates. In general, the time of existence of these solutions will go to zero as the surface tension parameter vanishes. We then make two additional estimates, in the case that a stability condition is satisfied by the initial data: we make an additional energy estimate which is uniform in the surface tension parameter, and we make an estimate for the difference of two solutions with different values of the surface tension parameter. These additional estimates allow the zero surface tension limit to be taken, showing that solutions of the initial value problem in the absence of surface tension are the limit of solutions of the initial value problem with surface tension as surface tension vanishes.     

Plotnikov——Toland板模型中水弹性孤立波的迎撞

通过奇异摄动方法研究了在薄冰层覆盖的不可压缩理想流体表面上传播的两个水弹性孤立波之间的迎面碰撞.借助特殊的 Cosserat 超弹性壳 理论以及Kirchhoff--Love 板理论,冰层由 Plotnikov--Toland板模型描述.流体运动采用浅水假设和Boussinesq 近似. 应用Poincaré--Lighthill--Kuo 方法进行坐标变形,进而渐近求解控制方程及边界条件, 给出了三阶解的显式表达. 可以观察到碰撞后的孤立波不会改变它们的形状和振幅. 波浪轮廓在碰撞之前是对称的, 而在碰撞之后变成不对称的并且在波传播方向上向后倾斜. 弹性板和流体表面张力减小了波幅. 图示比 较了本文与已有结果可知线性板模型可作为本文的一个特例.      

Stress concentration due to an array or hemispherical cavities at the surface of an elastic half-space

The paper investigates the perturbation in an otherwise uniform stress field in an elastic half-space due to a doubly-periodic array of small hemispherical holes at the free surface. The solution is obtained using three potential functions of double Fourier series form in Galerkin's strain potential solution, the coefficients of which are determined using the collocation method. The unperturbed field is taken to be one of uniform plane stress parallel to the free surface. Two special cases are studied—uniform tension and uniform shear stress. Numerical results for these cases can be generalized by superposition to give solutions for a general state of biaxial plane stress. It is found that, for both tension and shear, the maximum stress concentration occurs at the bottom of the holes. The stress concentration factor increases with the ratio of hole spacing to radius, approaching the known solution for a single hemispherical hole at large ratios.     

A marker‐and‐cell approach to free surface 2‐D multiphase flows

This work describes a methodology to simulate free surface incompressible multiphase flows. This novel methodology allows the simulation of multiphase flows with an arbitrary number of phases, each of them having different densities and viscosities. Surface and interfacial tension effects are also included. The numerical technique is based on the GENSMAC front‐tracking method. The velocity field is computed using a finite‐difference discretization of a modification of the Navier–Stokes equations. These equations together with the continuity equation are solved for the two‐dimensional multiphase flows, with different densities and viscosities in the different phases. The governing equations are solved on a regular Eulerian grid, and a Lagrangian mesh is employed to track free surfaces and interfaces. The method is validated by comparing numerical with analytic results for a number of simple problems; it was also employed to simulate complex problems for which no analytic solutions are available. The method presented in this paper has been shown to be robust and computationally efficient. Copyright © 2012 John Wiley & Sons, Ltd.     

Ratchetting under tension–torsion loadings: experiments and modelling

This paper is concerned with the mechanical behaviour of 316 austenitic stainless steel under multiaxial loadings and particular attention is paid to ratchetting under tension–torsion non-proportional loadings. First, a series of uniaxial tests and biaxial tests has been carried out in order to calibrate five different cyclic plasticity models based on an isotropic hardening rule and a non-linear kinematic hardening rule. It is shown that this class of models gives quite good agreement between the experimental and numerical results. Second, another series of ratchetting tests has been carried out under tension–torsion loadings in order to test the prediction capacities of the previous models. It is shown that whereas the models have been calibrated with similar loading paths, four of the five selected models give poor predictions.     

Spurious currents in finite element based level set methods for two‐phase flow

A study of spurious currents in continuous finite element based simulations of the incompressible Navier–Stokes equations for two‐phase flows is presented on the basis of computations on a circular drop in equilibrium. The conservative and the standard level set methods are used. It is shown that a sharp surface tension force, expressed as a line integral along the interface, can give rise to large spurious currents and oscillations in the pressure that do not decrease with mesh refinement. If instead a regularized surface tension representation is used, exact force balance at the interface is possible, both for a fully coupled discretization approach and for a fractional step projection method. However, the numerical curvature calculation introduces errors that cause spurious currents. Different ways to extend the curvature from the interface to the whole domain are discussed and investigated. The impact of using different finite element spaces and stabilization methods is also considered. Copyright © 2011 John Wiley & Sons, Ltd.     

Prediction of resonantly forced vibrations for suspended footbridges

Abstract

The aim of the present paper is to investigate the performance of suspended footbridges under pedestrian loads. Indeed, several Authors have underlined the possible activation of large amplitude oscillations into suspended footbridges due to the nonlinear behavior of the hangers. In fact, the last ones act as linear elastic springs in tension and do not react in compression. Consequently, if the whole suspended footbridge or parts of it undergo large amplitude oscillations, the initial hangers’ pretension stress may become zero and the slackening may start. In these cases, the stiffness of the footbridge deck decreases drastically, and a complex dynamic response may occur. Hence, the footbridge may show unexpected vertical and torsional oscillations that the “classical” models cannot predict; these models, in fact, assume a bilateral behavior for the suspended system and, consequently, no variations of the global stiffness during the motion. Here, the response of suspended footbridges is evaluated by using a continuous model obtained by adopting the nonlinear equivalent regularization technique proposed for long span suspended bridges. The dynamic analysis, performed by means of a perturbation method, shows the possibility of the coexistence of multiple solutions, some of which are characterized by high amplitudes and by the activation of the afore-described slackening phenomenon. The response is evaluated for several values of loading, mechanical and geometrical parameters, with the main aim of highlighting the characteristics and the stability of the investigated oscillations and obtaining information and/or design indications to prevent such phenomena.

Communicated by Dumitru Caruntu.     

On surface tension modelling using the level set method

The paper describes and compares the performance of two options for numerically representing the surface tension force in combination with the level set interface‐tracking method. In both models, the surface tension is represented as a body force, concentrated near the interface, but the technical implementation is different: the first model is based on a traditional level set approach in which the force is distributed in a band around the interface using a regularized delta function, whereas in the second, the force is partly distributed in a band around the interface and partly localized to the actual computational cells containing the interface. A comparative study, involving analysis of several two‐phase flows with moving interfaces, shows that in general the two surface tension models produce results of similar accuracy. However, in the particular case of merging and pinching‐off of interfaces, the traditional level set model of surface tension produces an error that results in non‐converging solutions for film‐like interfaces (i.e. ones involving large contact areas). In contrast, the second model, based on the localized representation of the surface tension force, displays consistent first‐order convergence. Copyright © 2008 John Wiley & Sons, Ltd.     

Effects of a transverse flows with variable viscosity in micropolar fluids

A regular perturbation analysis is presented for three laminar natural convection flows in micropolar fluids in liquids with temperature dependent viscosity: a freely-rising plane plume, the flow above a horizontal line source on an adiabatic surface (a plane wall plume) and the flow adjacent to a vertical uniform flux surface. While these flows have well-known power-low similarity solutions when the fluid viscosity is taken to be constant, they are non-similar when the viscosity is considered to a function of temperature. A single similar flow, that adjacent to a vertical isothermal surface, is also analysed for comparison in order to estimate the extent of validity of perturbation analysis. The formulation used here provides a unified treatment of variable viscosity effects on those four flows. Computed first-order perturbation quantities are presented for all four flows. Numerical results for velocity, angular velocity and thermal functions has been shown graphically or tabulated for different values of micropolar parameters. Received on 20 October 1997     

Investigation of droplet collisions for solutions with different solids content

The collision behaviour of droplets and the collision outcome are investigated for high viscous polymer solutions. For that purpose, two droplet chains produced by piezoelectric droplet generators are directed towards each other at a certain angle so that individual droplet pairs collide. For recording the collision event, one double-image and one high-speed CCD camera were used. One camera is positioned perpendicular to the collision plane recording the outcome of the collision, and the second camera is aligned parallel to the collision plane to assure that the droplet chains are exactly in one plane. A new approach for tracking droplets in combination with an extended particle tracking velocimetry algorithm has been developed. Time-resolved series of pictures were used to analyse the dynamics of droplet collisions. The three different water soluble substances were saccharose and 1-Ethenyl-2-pyrrolidone (PVP) with different molecular weights (K17, K30). The solvent was demineralised water. The solids contents ranged from 20 to 60 %, 5 to 25 % and 5 to 35 %, yielding dynamic viscosities in the range of 2–60 mPa s. Results were collected for different pairs of impact angles and Weber numbers in order to establish common collision maps for characterising the outcomes. Here, relative velocities between 0.5 and 4 m/s and impact parameters in the interval from 0 to 1 for equal-sized droplets (Δ = 1) have been investigated. Additionally, satellite formation will be discussed exemplarily for K30. A comparison with common models of different authors (Ashgriz and Poo in J Fluid Mech 221:183–204, 1990; Estrade et al. in Int J Heat Fluid Flow 20:486–491, 1999) mainly derived for low viscous droplets revealed that the upper limit of their validity is given by an Ohnesorge number of Oh = 0.115 and a capillary number of Ca = 0.577. For higher values of these non-dimensional parameters and hence higher dynamic viscosities, these models are unable to predict correctly the boundaries between collision scenarios. The model proposed by Jiang et al. (J Fluid Mech 234:171–190, 1992), which includes viscous dissipation, is able to predict the boundary between coalescence and stretching separation for higher viscosities (i.e. Oh > 0.115 and Ca > 0.577). However, the model constants are not identical for different solution properties. As a conclusion, an alteration of the collision appearance takes place because of the relative importance between surface tension and viscosity.     

Convective magnetohydrodynamic two fluid flow and heat transfer in an inclined channel

 The problem of fully developed free convection two fluid magnetohydrodynamic flow in an inclined channel is investigated. The governing momentum and energy equations are coupled and highly nonlinear due to dissipation terms, solutions are found employing perturbation technique for small values of Pr · Ec (=ɛ) the product of Prandtl number and Eckert number. Effects of Grashof number, Hartmann number, inclination angle, the ratios of electrical conductivities, viscosities and heights of two fluids on the flow are explored. It is observed that the flow can be controlled effectively by suitable adjustment of the values for the ratios of heights, electrical conductivities and viscosities of the two fluids. Received on 10 December 1999     

横观各向同性压电材料中裂纹问题的边界元法

采用Somigiliana公式给出了三维横观各向同性压电材料中的非渗漏裂纹问题的一般解和超奇异积分方程,其中未知函数为裂纹面上的位移间断和电势间断．在此基础上,使用有限部积分和边界元结合的方法,建立了超奇异积分方程的数值求解方法,并给出了一些典型数值算例的应力强度因子和电位移强度因子的数值结果,结果令人满意．     

Particle-driven gravity currents: asymptotic and box model solutions

We employ shallow water analysis to model the flow of particle-driven gravity currents above a horizontal boundary. While there exist similarity solutions for the propagation of a homogeneous gravity current, in which the density difference between the current and ambient is constant, there are no such similarity solutions for particle-driven currents. However, because the settling velocity of the particles is often much less than the initial velocity of propagation of these currents, we can develop an asymptotic series to obtain the deviations from the similarity solutions for homogeneous currents which describe particle-driven currents. The asymptotic results render significant insight into the dynamics of these flows and their domain of validity is determined by comparison with numerical integration of the governing equations and also with experimental measurements. An often used simplification of the governing equations leads to `box' models wherein horizontal variations within the flow are neglected. We show how to derive these models rigorously by taking horizontal averages of the governing equations. The asymptotic series are then used to explain the origin of the scaling of these `box' models and to assess their accuracy.     

Stability and capillary dynamics of circular vortex sheets

We investigate the motion of circular vortex sheets with surface tension. A linear stability analysis shows that high modes of the circular vortex sheet are stabilized by surface tension, and the sheet is stable if surface tension is larger than a critical value. The modes of perturbations, n = 1 and 2, are always stable, regardless of surface tension, and the mode n = 3 is also stable for large surface tension. The numerical results show that a stable vortex sheet rotates and oscillates weakly. The oscillations of each mode of the interface mainly consist of two travelling waves of different frequencies in time. The amplitude and the period of the oscillation depend on the mode of the perturbation and surface tension. We also perform long-time computations for the unstable evolution of circular sheets. For a high Weber number, ripples are produced on the sheets, as well as pinching and self-intersection. It is found that the appearance of ripples is associated with the growth of noise. For an intermediate Weber number, the sheet evolves to an exotic structure with small spikes on the fingers, while for a low Weber number, it is nonlinearly stable.     

The evolution and interaction of Marangoni-Bénard solitary waves

Marangoni-Bénard convection is the process by which oscillatory waves are generated on an interface due to a change in surface tension. This process, which can be mass or temperature driven is described by a perturbed Korteweg-de Vries (KdV) equation. The evolution and interaction of solitary waves generated by Marangoni-Bénard convection is examined. The solitary wave with steady-state amplitude, which occurs when the excitation and friction terms of the perturbed KdV equation are in balance is found to second-order in the perturbation parameter. This solitary wave has a fixed amplitude, which depends on the coefficients of the perturbation terms in the governing equation. The evolution of a solitary wave of arbitrary amplitude to the steady-state amplitude is also found, to first-order in the perturbation parameter. In addition, by using a perturbation method based on inverse scattering, it is shown that the interaction of two solitary waves is not elastic with the change in wave amplitude determined. Numerical solutions of the perturbed KdV equation are presented and compared to the asymptotic solutions.     

被均匀流缓慢调制的有限水深毛细重力波

非线性的存在会产生高次谐波,这些谐波又反作用于原来的低次谐波,使波幅发生缓慢变化,从而产生缓慢调制现象.这里从考虑均匀流作用下的毛细重力水波基本方程出发,在不可压缩、无旋、无黏条件假设下,使用多重尺度分析方法推导出了在均匀流影响下有限深水毛细重力波振幅所满足的非线性Schr¨odinger方程(NLSE).分析了NLSE解的调制不稳定性.给出了毛细重力波调制不稳定的条件和钟型孤立波产生的条件.分析了无量纲最大不稳定增长率随无量纲水深和表面张力的变化趋势.同时给出了无量纲不稳定增长率随无量纲微扰动波数变化的曲线,呈现出了先增大后减小的趋势.最后指出均匀顺流减小了无量纲不稳定增长率及最大增长率,逆流增大了它们.由表面张力作用产生的毛细波及重力与表面张力共同作用产生的毛细重力波,与流的相互作用可以改变海表粗糙度和海洋上层流场结构,进而影响海气界面动量、热量及水汽的交换.了解海表这些短波动力机制,对卫星遥感的精确测量、海气相互作用的研究及海气耦合模式的改进等有重要意义.