液滴平壁铺展过程中的滞后效应及力学机制研究

研究了液滴平壁铺展过程中的接触角滞后效应,从接触线附近流体的压力、速度以及能量分布等角度考虑滞后效应的成因和变化规律.在此基础之上分析了固体表面粗糙度对滞后效应的影响,并借助部分三维形貌参数(ISO 25178)建立了固体表面形貌与接触角滞后效应之间的量化关系.为了研究以上内容,应用数值仿真软件建立了液滴铺展动力学模型,并结合固体表面上的滞后性实验进行了相关验证.研究结果表明:由于表面粗糙度的存在,液滴铺展至平衡位置时,位于铺展前沿的液体分子被钉扎在固体表面的凹坑或低谷中,使得前沿接触角逐渐增大,后沿接触角逐渐减小,接触角发生滞后;驱动液滴铺展的Laplace压力和自身重力与阻碍液滴铺展的黏性阻力之间的平衡关系,是接触角发生滞后的主要力学机制.另外,实验结果表明接触角滞后效应与固体表面形貌密切相关,具有相同表面粗糙度(Sa)的固体表面,由于表面形貌不同接触角滞后效应可能存在明显的差异.     

Capillary Drops: Contact angle hysteresis and sticking drops

This paper is concerned with the homogenization of a capillary equation for liquid drops lying on an inhomogeneous solid plane. We show in particular that the homogenization of the Young–Laplace law leads to a contact angle condition of the form , which justifies the so-called contact angle hysteresis phenomenon.     

Numerical Solution of Contact Problems using Level Set Methods on Voxel Discretizations

Fast solvers performing on a regular grid are often used for problems in elasticity, in order to avoid expensive mesh generations. However, if overlaps between solid materials occur without any interactions, these might deteriorate the results, especially for the stresses. Therefore, we want to present an approach for developing numerical methods for contact problems on a regular mesh with the help of signed distance data and multi-material voxels. In this contribution we will focus on problems in linear elastostatics with contact between several different bodies. Finally, we present the results from a numerical test for the two dimensional Hertz problem, solved on a triangular regular mesh. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Contact deformation of an elastic composition of a half-plane and a strip of variable width

Asymptotic dependences of the deformation on the contact stresses are derived for a strip of variable width bound to an elastic half-plane. Similar relations were previously obtained in [1] for a strip of constant width.     

A moving grid finite-element method using grid deformation

A new front tracking method based on grid deformation is combined with a streamline upwind Petrov–Galerkin finite element method to produce an adaptive method for application to convection-dominated transport equations. In this work, we examine the practicality of this approach and demonstrate the method on one-dimensional examples. © 1995 John Wiley & Sons, Inc.     

Spatial effects in problems on the deformation of cylindrical shells of anisotropic composites

Mechanics of Composite Materials -     

CFD studies on burner secondary airflow

In many fossil power plants operating today, there is insufficient means to assure the proper balancing of the secondary airflows between the individual burners of wall-fired units. This mismatch leads to decreased boiler efficiency and increased emissions. In this study, a computational fluid dynamics (CFD) modeling of a fossil power plant wind box was performed. The model solved the three-dimensional Reynolds averaged Navier–Stokes equations with the k–ε turbulence model. The CFD results were validated by the experimental data taken from a 1/8th scale model of a wall-fired fossil unit. Simulations under various mass flow rates specified at inlet, various baffle positions and two opening conditions of the burners were obtained to identify the optimum design in terms of the equalization of the secondary airflow through the burners. This study demonstrated that the combination of experimental and CFD approach can be an effective tool in the research of burner secondary airflow balancing.     

Investigation of elastoplastic effects of cables under large spatial deformation

Cables are complex components consisting of a multi-layer structure and various materials. The structural setup includes for example conducting wires, isolating shields and protecting sheaths. This leads to several inelastic effects under large deformations like pull-out of wires, delamination of layers or friction between the constituents. The materials used in cables belong to different material classes and consequently show different behavior under load. Elastoplastic behavior has to be expected for metallic wires, whereas polymer layers behave viscoelastically. The combination of these inelastic effects caused by the structure and constituents of cables motivates the inclusion of inelasticity in the material model on a phenomenological level. Since cables are flexible, slender structures, they can be described physically correctly by the theory of Cosserat rods. In this context, the constitutive equations are formulated in terms of the sectional quantities. The related model parameters have to be determined in suitable experiments. As cables undergo large multiaxial deformations in applications, uniaxial experiments are not sufficient for their characterization. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Progressive necking and orientation effects in the deformation of kapron

The stretching of kapron has been investigated by x-ray structural analysis using an estimate of the degree of orientation based on development of the Debye ring. It is established that with the progressive development of a first, second and third neck there is a decrease in the width of the texture maxima (increase in the degree of preferred orientation). The width of the texture maxima varies with the relative reduction of the cross section approximately in accordance with a linear law.Mekhanika Polimerov, Vol. 1, No. 4, pp. 7–11, 1965     

Deformed gap space using macro-micro FEA model and transferred into a CFD model

Using a cylindrical nozzle and seat of a Pressure Relief Valve (PRV) the surface form and waviness is modelled using actual metrological data i. e. average surface form and waviness (W_a and W_sm) in a symmetry manner. To model the surface waviness the technique used is based on the summing technique created by Tsukizoe & Hisakado [1, 2] for micro contact analysis. Due to the actual surface form measurements being in the micro-meter range, the model is required to incorporate micro and macro-meter dimensions. The material in question is stainless steel. The deformed finite element analysis model is then transferred into a CAD geometry allowing the void space to be meshed and solved using computational fluid dynamics. © 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Patient-specific model of lung deformation using spatially dependent constitutive parameters

Breathing-induced spatially dependent lung deformation is predicted using patient-specific elastic properties with the contact–impact analysis model. The lung geometry is derived from 4D CT scan data of real patients. The spatially varying Young’s modulus for the patient is obtained from a previous study that used inverse deformation of the lung. The compact–impact analysis is implemented using the finite element method. The predicted lung deformation is compared with the results based on linear elasticity. The results are consistent with physiology, indicating large deformations near the diaphragm and smaller values at remote locations on the lobe. The effect of non-linearity of elastic property is most significant at the remote locations where the diaphragm-induced deformation is significantly attenuated.     

CFD methodology for sedimentation tanks: The effect of secondary phase on fluid phase using DPM coupled calculations

Computational Fluid Dynamics (CFD) methods are employed in order to simulate the 3D hydrodynamics and flow behaviour in a sedimentation tank. Unlike most of the previous numerical investigations, in the present paper the momentum exchange between the primary and the secondary phase is taken into account, using a Lagrangian method (discrete phase model) with two-way coupled calculations. By computing particle trajectories the proposed numerical model can track the momentum gained or lost by the particle stream that follows that trajectory and these quantities can be incorporated in the subsequent continuous phase calculations. Thus, while the continuous phase always impacts the discrete phase, the effect of the discrete phase trajectories on the continuum can be incorporated. This interchange affects fluid velocity, especially in the case of large particles sizes, which have a greater relaxation time in relation to the characteristic time of the tank. The present investigation compares a series of numerical simulations for a sedimentation tank with varying particle diameters and volume fractions, in order to identify the influence of the secondary phase to the primary phase and vice-versa and the way that this influence affects the efficiency of the tank.     

Predicting initial erosion during the hole erosion test by using turbulent flow CFD simulation

CFD simulation with enhanced modeling of turbulence and near-wall treatment is used to model water–clay mixtures flowing through a cylindrical pipe domain. Effects on the wall-shear stress resulting from varying water clay content and applied hydraulic gradient are analyzed. Various parametric studies were performed and had shown that the two-dimensional modelling introduced in the present study does not yield a uniform wall-shear stress along the pipe wall and that clay concentration affects significantly the wall-shear stress value. This is in contrast with the common hypothesis used in one-dimensional modeling approaches where this stress is assumed constant and which gives rise to uniform erosion along the pipe wall. The obtained results had enabled predicting more realistically erosion amount and had allowed for understanding the irregular eroded hole wall shape as observed experimentally after performing the standard hole erosion test.     

Contact Elements on Fibered Manifolds

For every product preserving bundle functor T ^μ on fibered manifolds, we describe the underlying functor of any order (r, s, q), s ≥ r ≤ q. We define the bundle K_k,l^r,s,q YK_{k,l}^{r,s,q} Y of (k, l)-dimensional contact elements of the order (r, s, q) on a fibered manifold Y and we characterize its elements geometrically. Then we study the bundle of general contact elements of type μ. We also determine all natural transformations of K_k,l^r,s,q YK_{k,l}^{r,s,q} Y into itself and of T( K_k,l^r,s,q Y )T\left( {K_{k,l}^{r,s,q} Y} \right) into itself and we find all natural operators lifting projectable vector fields and horizontal one-forms from Y to K_k,l^r,s,q YK_{k,l}^{r,s,q} Y .     

Contact processes on scale-free networks

We investigate the contact process on random graphs generated from the configuration model for scale-free complex networks with the power law exponent β E （2, 3]. Using the neighborhood expansion method, we show that, with positive probability, any disease with an infection rate λ 〉 0 can survive for exponential time in the number of vertices of the graph. This strongly supports the view that stochastic scale-free networks are remarkably different from traditional regular graphs, such as, Z^d and classical Erdos-Renyi random graphs.     

Analysis of functionally graded plates using higher order shear deformation theory

This work addresses a static analysis of functionally graded material (FGM) plates using higher order shear deformation theory. In the theory the transverse shear stresses are represented as quadratic through the thickness and hence it requires no shear correction factor. The material property gradient is assumed to vary in the thickness direction. Mori and Tanaka theory (1973) [1] is used to represent the material property of FGM plate at any point. The thermal gradient across the plate thickness is represented accurately by utilizing the thermal properties of the constituent materials. Results have been obtained by employing a C° continuous isoparametric Lagrangian finite element with seven degrees of freedom for each node. The convergence and comparison studies are presented and effects of the different material composition and the plate geometry (side-thickness, side–side) on deflection and temperature are investigated. Effect of skew angle on deflection and axial stress of the plate is also studied. Effects of material constant n on deflection and the temperature distribution are also discussed in detail.     

Optimising the design of fume extraction hoods using a combination of engineering and CFD modelling

Fume and hygiene hoods are widely used to prevent fugitive emissions from charge ports, tap holes and many other openings in mineral processing and smelting vessels. The highly buoyant nature of the fume combined with often complex geometries make the design of these hoods difficult with traditional engineering tools. However, by combining the traditional engineering approach with computational fluid dynamics (CFD) techniques, a clear understanding of the shortfalls of an existing system can be obtained, and an optimised hood design can be achieved. This paper reports on a combined engineering and CFD analysis of a fume extraction system for a zinc slag fumer charge port. The engineering model revealed that the existing plant components (bag house and fan) were not capable of capturing the required amount of fume, and that the original hood design was flawed. The CFD model was then used to predict the fume capture and emission from the existing hood. CFD model predictions showed that increasing the draft flow rate by an order of magnitude would only give a marginal improvement in fume capture. Using findings of both the models enabled a new fume capture hood to be designed. CFD analysis of the new hood revealed that a significant improvement in fume capture is possible. Construction and installation of the hood has been performed and a 65% reduction in fume emission was achieved, thus significantly mitigating a long-standing emission problem.