An improved diffuse interface method for three-dimensional multiphase flows with complex interface deformation

An improved diffuse interface (DI) method is proposed for accurately capturing complex interface deformation in simulations of three-dimensional (3D) multiphase flows. In original DI methods, the unphysical phenomenon of interface thickening or blurring can affect the accuracy of numerical simulations, especially for flows with large density ratio and high Reynolds number. To remove this drawback, in this article, an interface-compression term is introduced into the Cahn-Hilliard equation to suppress the interface dispersion. The additional term only takes effect in the interface region and works normal to the interface. The difference of the current method from the previous work is that the compression rate can be adjusted synchronously according to the magnitude of local vorticity, which is strongly correlated to the interface dispersion and changes with the computational time and interface position. Numerical validations of the proposed method are implemented by simulating problems of Laplace law, Rayleigh-Taylor instability, bubble rising in a channel, and binary droplet collision. The obtained results agree well with the analytical solutions and published data. The numerical results show that the phenomenon of interface dispersion is suppressed effectively and the tiny interfacial structures in flow field can be captured accurately.     

Motion of deformable drops in pipes and channels using Navier–Stokes equations

The motion of deformable drops in pipes and channels is studied using a level set approach in order to capture the interface of two fluids. The interface is described as the zero level set of a smooth function, which is defined to be the signed normal distance from the interface. In order to solve the Navier–Stokes equations, a second‐order projection method is used. The dimensionless parameters of the problem are the relative size of the drop to the size of the pipe or channel cross‐section, the ratio of the drop viscosity to the viscosity of the suspending fluid and the relative magnitude of viscous forces to the surface tension forces. The shape of the drop, the velocity field and the additional pressure loss due to the presence of the drop, varying systematically with the above‐mentioned dimensionless parameters, are computed. Copyright © 2000 John Wiley & Sons, Ltd.     

液滴撞击结构液靶成坑特点分析

本文介绍了液滴对多层液靶的撞击实验研究,实验中液体从不同高度滴落,同时改变液靶的液层厚度,通过高速摄影机记录撞击成坑的过程,将胶片扫描数据输入计算机,并通过图像软件处理获得最终成坑数据。由实验结果分析了液滴撞击多层液靶成坑的特点,在一定的撞击能量情况下,撞击成坑现象发生于第一层液靶内,界面能在多层液靶不互溶时成为影响成坑状态的一个重要参数。根据能量守恒定律,并通过对撞击能量转化过程的分析,得出了撞击成坑最大坑径公式。     

Simulation of incompressible multiphase flows with complex geometry using etching multiblock method

The incompressible two-phase flows are simulated using combination of an etching multiblock method and a diffuse interface (DI) model, particularly in the complex domain that can be decomposed into multiple rectangular subdomains. The etching multiblock method allows natural communications between the connected subdomains and the efficient parallel computation. The DI model can consider two-phase flows with a large density ratio, and simulate the flows with the moving contact line (MCL) when a geometric formulation of the MCL model is included. Therefore, combination of the etching method and the DI model has potential to deal with a variety of two-phase flows in industrial applications. The performance is examined through a series of numerical experiments. The convergence of the etching method is firstly tested by simulating single-phase flows past a square cylinder, and the method for the multiphase flow simulation is validated by investing drops dripping from a pore. The numerical results are compared with either those from other researchers or experimental data. Good agreement is achieved. The method is also used to investigate the impact of a droplet on a grooved substrate and droplet generation in flow focusing devices.     

Acoustic emission assessment of interface cracking in thermal barrier coatings

In this paper, acoustic emission (AE) and digital image correlation methods were applied to monitor interface cracking in thermal barrier coatings under compression. The interface failure process can be identified via its AE features, including buckling, delamination incubation and spallation. According to the Fourier transformation of AE signals, there are four different failure modes: surface vertical cracks, opening and sliding interface cracks, and substrate deformation. The characteristic frequency of AE signals from surface vertical cracks is 0.21 MHz, whilst that of the two types of interface cracks are 0.43 and 0.29 MHz, respectively. The energy released of the two types of interface cracks are 0.43 and 0.29 MHz, respectively. Based on the energy released from cracking and the AE signals, a relationship is established between the interface crack length and AE parameters, which is in good agreement with experimental results.     

AN EXPERIMENTAL METHOD BASED ON ENERGY RELEASE RATE OF INTERFACIAL FRACTURE

基于对断裂力学常用实验方法的研究,结合界面断裂问题的特殊性,以断裂力学为理论基础,通过能量释放率建立了界面断裂测量的实验分析方法,并且利用文献中的实验数据进行了验证,取得了良好的一致性.该方法通过测量试验件的载荷位移关系,利用裂纹扩展过程中的能量变化关系得到该裂纹长度下的临界能量释放率;在此基础上,根据试件的阻抗能量曲线预测结构的最大承载能力.该方法以能量释放率为理论基础,为界面裂纹的强度分析提供了合理的手段,基于能量角度建立的实验分析方法也具有良好的实用性和适用性.     

Two-phase pressure drop in multiple thick- and thin-orifice plates

Two-phase flow pressure changes through singularities such as sudden expansion and sudden contraction of thick- and thin-orifice plates were modeled. The modeling was based on the reversible and irreversible losses through contractions and expansions. The volume-averaged momentum equation and the reversible mechanical energy equation were used to evaluate the irreversibilities. Local slip ratios, which were necessary for the prediction of pressure drop through these singularities, were correlated from a large number of experimental data. To check the validity of the analytical predictions, an experimental test section was designed, and experiments were performed to produce benchmark pressure-drop data for thin- and thick-plate orifices. The working fluid used for these experiments was R-113. The predictive methods developed agree well with the experiments by the authors and with a wide range of two-phase flow test results obtained by others for steam-water systems. A parametric study shows the relative importance of geometry and of the flow variables such as quality, liquid-to-vapor density, and viscosity ratios on the pressure drop multipliers conventionally used in two-phase flows.     

Fluid–structure interaction of a splitter plate in a convergent channel

The fluid–structure interaction (FSI) of a splitter plate in a convergent channel flow is studied by measuring both the flow field and the plate vibration. Particle Image Velocimetry (PIV) measurements show that the wake generated by the plate is characterized by cellular vortex shedding. Mean and RMS velocities presented in the plane normal to the main flow direction visualize clearly the cellular structure and related secondary flows. To evaluate the energy and spatial organization of the vortex shedding, spectral and correlation estimation methods are adapted to the PIV data. By presenting the spanwise variation of the streamwise spectra along the trailing edge, the nature of the cellular vortex shedding becomes evident. 2D space-correlation function reveals that the shedding in two neighboring cells occurs in a 180-degree phase shift. The vibration of the plate is studied with Digital Imaging (DI) and Laser Vibrometer (LV). The DI is based on images measured by the PIV system. An image-processing algorithm is used to detect the plate tip location and velocity simultaneously with the estimation of the fluid velocity field. The LV is used for the time-resolved measurement of the plate vibration. The results show that the plate vibrates in a very distinct mode characterized by a spanwise standing wave along the plate-trailing edge. This mode, in turn, causes the cellular vortex shedding.     

共混/填充高分子复合材料界面力学行为实验研究进展

本文首先介绍了高聚物复合材料界面微观结构的特点和界面破坏的微观形式, 在此基础上,详细回顾了粒子填充高分子复合材料和共混高分子合金两类材料的界面力学性能实验研究的进展,总结了研究粒子/基体界面力学性能的实验方法和用这些方法获得的实验结果。介绍的实验方法从宏观拉伸实验到细观在位拉伸实验, 从定性的断口形貌分析到定量的界面强度和粘结能计算, 从不同层次和角度阐述了研究界面性能各种方法的优缺点。最后, 作者根据综合分析的结果, 指出了实验研究界面力学性能遇到的主要困难来自实验仪器的局限性和实验方法的匮乏性。为了获得可靠的界面性能, 一方面应该设计出利于表征的界面, 另一方面需要在实验技术上有所创新。相信该文对于深入了解和研究共混/填充高分子复合材料的界面问题, 以及其它类型材料的界面问题具有一定的借鉴意义。      

Equivalent inclusions in micromechanics with interface energy effect

In order to apply classical micromechanics in predicting the effective properties of nanocomposites incorporating interface energy, a concept of equivalent inclusion(EI) is usually adopted. The properties of EI are obtained by embedding a single inclusion with the interface into an infinite matrix. However, whether such an EI is universal for different micromechanics-based methods is rarely discussed in the literature. In this paper, the interface energy theory is used to study the applicability of the above mentioned EI. It is found that some elastic properties of the EI are related only to the properties of the inclusion and the interface, whereas others are also related to the properties of the matrix. The former properties of the EI can be applied to both the classical Mori-Tanaka method(MTM) and the generalized self-consistent method(GSCM). However, the latter can be applied only to the MTM. Two kinds of new EIs are proposed for the GSCM and used to estimate the effective mechanical properties of nanocomposites.     

双层金属纳米板界面能密度的尺寸效应

界面能密度是表征纳米复合材料与结构界面力学性质的重要物理量.采用分子动力学方法计算了不同面心立方金属晶体构成的双材料纳米薄板结构的界面能密度,分析了界面晶格结构形貌变化及界面效应对原子势能的影响.结果表明:双材料纳米薄板界面具有周期性褶皱状疏密相间的晶格结构形貌,界面上原子势能亦呈现周期性分布特性,而靠近界面的两侧原子势能与板内原子势能具有明显差异.拉格朗日界面能密度和欧拉界面能密度均随双层薄板厚度的增加而增加,最终趋向于块体双材料结构的界面能密度.     

基于分维理论的爆炸焊接界面形貌的定量表征

以爆炸焊接界面形貌为研究对象,利用分形理论研究其界面特征。运用三维超景深显微镜获得界面形貌图像,利用Matlab图像分析技术对界面形貌进行二值化处理,基于分形理论计算图像分维值以及多重分维谱。由分维值及多重分形谱的物理意义可知,轮廓分维值是对界面不规则程度宏观度量,而多重分形谱能定量表征界面的起伏程度及高度分布。因而,通过对界面结构进行分形几何分析,可实现界面形貌定量表征,弥补当前定性分析的不足。     

Determinations of both length scale and surface elastic parameters for fcc metals

In the present research, a systematical study of trans-scale mechanics theory is performed. The surface/interface energy density varying with material deformation is considered, and the general surface/interface elastic constitutive equations are derived. New methods to determine the material length scale parameter and the surface elastic parameters based on a simple quasi-continuum method, i.e. the Cauchy–Born rule, are developed and applied to typical fcc metals. In the present research, the material length parameters will be determined through an equivalent condition of the strain energy density calculated by adopting the strain gradient theory and by adopting the Cauchy–Born rule, respectively. Based on the surface constitutive equations obtained in the present research, the surface elastic parameters are calculated by using the Gibbs definition of surface energy density and the Cauchy–Born rule method.     

Flow of linear molecules through a 4:1:4 contraction–expansion using non-equilibrium molecular dynamics: Extensional rheology and pressure drop

In this work, non-equilibrium molecular dynamics simulations are used to generate the flow of linear polymer chains (monomer-springs with FENE potential) and a Lennard–Jones fluid (Newtonian fluid) through a contraction–expansion (4:1:4) geometry. An external force field simulating a constant pressure gradient upstream the contraction region induces the flow, where the confining action of the walls is represented by a Lennard–Jones potential. The equations of motion are solved through a multiple-step integration algorithm coupled to a Nosé-Hoover dynamics [S. Nose, A unified formulation of the constant temperature molecular dynamics methods, J. Chem. Phys. 81 (1984) 511–519], i.e., to simulate a thermostat, which maintains a constant temperature. In this investigation, we assume that the energy removed by the thermostat is related to the viscous dissipation along the contraction–expansion geometry. A non-linear increasing function between the pressure drop and the mean velocity along the contraction for the linear molecules is found, being an order of magnitude larger than that predicted for the Lennard–Jones fluid. The pressure drop of both systems (the linear molecules and Lennard–Jones fluid) is related to the dissipated energy at the contraction entry. The large deformation that the linear molecules experience and the evolution of the normal stress at the contraction entry follow a different trajectory in the relaxation process past the contraction, generating large hysteresis loops. The area enclosed by these cycles is related to the dissipated energy. Large shear stresses developed near the re-entrant corners as well as the vortex formation, dependent on the Deborah number, are also predicted at the exit of the contraction. To our knowledge, for the first time, the excessive pressure losses found in experimental contraction flows can be explained theoretically.     

Deflection criterion of a crack approaching a bi-material interface obliquely

Deflection criterion for oblique cracks terminating at a bi-material interface was established numerically based on a remote loading condition where the crack deflection event took place well within the K-dominant stress field. The criterion was described in terms of the ratio of the energy release rate of a deflected crack (G _d) to the maximum energy release rate of a penetrated crack (). The criterion was markedly more conservative than the existing solution based on wedge loading which did not converge with respect to limited number of a/L ratios in the literature (a is the size of the putative crack; L is the loading distance). Further, the criterion established herein for the cracks slightly oblique from the normal direction to the interface was more conservative than the crack normal to the interface.     

A PLIC volume tracking method for the simulation of two‐fluid flows

Numerical methodologies for computer simulations of two‐fluid flows are presented. These methodologies fall into the category of volume tracking methods with piecewise‐linear interface calculation (PLIC). The scope of this work is limited to laminar flows of immiscible, non‐reacting, incompressible Newtonian fluids, without phase change, in planar two‐dimensional geometries. The following novel or enhanced procedures are proposed: a parallelogram scheme for multidimensional advection of the volume‐fraction field; a circle‐fit technique for the orientation of the interface segments and the calculation of curvature; a novel contact angle treatment; and a staggered formulation for volumetric body forces that can accurately balance pressure forces in the vicinity of the interface. In addition, surface‐tension‐derived and hydrostatic‐derived pressure adjustments are introduced as a means of accurately calculating pressure forces in cells that contain the interface, so as to minimize the non‐physical flows that afflict many available volume tracking methods. The proposed method is validated using four test problems that involve simulations of pure advection, a static drop, an oscillating bubble, and a static meniscus. Copyright © 2006 John Wiley & Sons, Ltd.     

Application of the Density-Functional Method for Numerical Simulation of Flows of Multispecies Multiphase Mixtures

Numerical examples of application of the density functional used to describe isothermal flows of two-phase two-species mixtures are given. The following flows are calculated in a two-dimensional formulation: impact of a drop on a liquid layer, breakdown of a drop in the velocity field of the Couette flow, formation of the wetting angle of a drop on a solid surface, and development of the Rayleigh–Taylor and Kelvin–Helmholtz instabilities at the gas–liquid interface.     

爆炸驱动液体分层现象和液滴尺寸模拟

爆炸驱动液体介质外界面的分散和破碎是气溶胶云团形成的重要过程。采用基于维数分裂的欧拉程序和Youngs混合界面处理方法,对中心药爆炸驱动甘油和水介质流场的液体分层现象进行了数值模拟。结合试验结果推断提出了液滴形成过程的三种并存机制:外层射流破碎、内层R-T失稳和中间液层"空化"破碎,分别建立了不同液层破碎液滴的尺寸模拟方法。对比给出抛撒甘油和水装置初级液滴的尺寸分布及最外层理论射流量。     

Characterization of interfacial failure using a reflected light technique

A simple optical method of observing and tracking interface failure in transparent matrix composites is developed to extend the range of material systems that can be tested with the cruciform test for measuring the normal strength of a bi-material interface. This technique of detecting debonds using reflected light is demonstrated on two model material systems with fibers 140 and 15 μm in diameter. The reflected light technique is more reliable and provides information about the interface that cannot be obtained with previous methods of debond detection, such as photoelasticity, surface strain gages, and acoustic emission. Using this technique, debond initiation, location, length, and shape are measured as a function of applied load; thus, all the required parameters for calculating the normal strength of the interface, as well as the interfacial energy release rate are acquired.