Problem of Low-frequency Localized Oscillations in a Thin Film with Growing Islands

Influence of localization waves to islands growing on a thin film is investigated. The film is modelled as a fluid layer covered by an inertial surface with the variable density of mass and surface tension. Mathematically, the problem is reduces to analysis of a system of non-linear equations describing the growth of island nuclei and wave propagation in the films. The existence of trapped modes for the corresponding frequency-domain problem is established. We show that for large time wave localization near islands gives some contribution in the increase of the velocity of island growth.     

Atomistic simulation of stress evolution during island growth

We report the results of a series of hybrid molecular dynamics simulations of the growth of islands on a substrate for several different island/substrate interface energies. When the interface energy is small, the islands tend to be thin and broad and the magnitude of the compressive stress-thickness product is relatively large. As the interface energy increases, the islands become taller and thinner and the magnitude of the compressive stress-thickness product decreases. This trend is consistent with experimental observations. The island aspect ratio dependence on interface energy follows from consideration of the equilibrium wetting angle. The effect of interface energy on the stress-thickness product shows that the island shape, surface/interface stresses and island stresses are self-equilibrated. A simple theory is developed that shows that the stress-thickness product is simply proportional to the substrate coverage and the substrate surface stress. The present simulations yield a simple, accurate, validated theory for stress development during the pre-coalescence stage of film growth.     

Kinetic wrinkling of an elastic film on a viscoelastic substrate

A compressed elastic film on a compliant substrate can form wrinkles. On an elastic substrate, equilibrium and energetics set the critical condition and select the wrinkle wavelength and amplitude. On a viscous substrate, wrinkle grows over time and the kinetics selects the fastest growing wavelength. More generally, on a viscoelastic substrate, both energetics and kinetics play important roles in determining the critical condition, the growth rate, and the wavelength. This paper studies the wrinkling process of an elastic film on a viscoelastic layer, which in turn lies on a rigid substrate. The film is elastic and modeled by the nonlinear von Karman plate theory. The substrate is linear viscoelastic with a relaxation modulus typical of a cross-linked polymer. Beyond a critical stress, the film wrinkles by the out-of-plane displacement but remains bonded to the substrate. This study considers plane strain wrinkling and neglects the in-plane displacement. A classification of the wrinkling behavior is made based on the critical conditions at the elastic limits, the glassy and rubbery states of the viscoelastic substrate. Linear perturbation analyses are conducted to reveal the kinetics of wrinkling in films subjected to intermediate and large compressive stresses. It is shown that, depending on the stress level, the growth of wrinkles at the initial stage can be exponential, accelerating, linear, or decelerating. In all cases, the wrinkle amplitude saturates at an equilibrium state after a long time. Subsequently, both amplitude and wavelength of the wrinkle evolve, but the process is kinetically constrained and slow compared to the initial growth.     

Cracking mechanisms in lithiated silicon thin film electrodes

The massive cracking of silicon thin film electrodes in lithium ion batteries is associated with the colossal volume changes that occur during lithiation and delithiation cycles. However, the underlying cracking mechanism or even whether fracture initiates during lithiation or delithiation is still unknown. Here, we model the stress generation in amorphous silicon thin films during lithium insertion, fully accounting for the effects of finite strains, plastic flow, and pressure-gradients on the diffusion of lithium. Our finite element analyses demonstrate that the fracture of lithiated silicon films occurs by a sequential cracking mechanism which is distinct from fracture induced by residual tension in conventional thin films. During early-stage lithiation, the expansion of the lithium-silicon subsurface layer bends the film near the edges, and generates a high tensile stress zone at a critical distance away within the lithium-free silicon. Fracture initiates at this high tension zone and creates new film edges, which in turn bend and generate high tensile stresses a further critical distance away. Under repeated lithiation and delithiation cycles, this sequential cracking mechanism creates silicon islands of uniform diameter, which scales with the film thickness. The predicted island sizes, as well as the abrupt mitigation of fracture below a critical film thickness due to the diminishing tensile stress zone, is quantitatively in good agreement with experiments.     

Instability of an elastic film on a viscous layer

A flat, compressed elastic film on a viscous layer is unstable. The film can form wrinkles to reduce the elastic energy. In this paper, we are interested in the two-dimensional models for thin films bonded to a viscous layer and in particular we focus on generic instabilities evidenced in this context by Suo and coworkers [Huang, Z., Hong, W., Suo, Z., 2005. Non linear analyses of wrinkles in a film bonded to a compliant substrate. J. Mech. Phys. Solids 53, 2101–2118; Lo, Y.H., 1991. New approach to grow pseudomorphic structures over the critical thickness. Appl. Phys. Lett. 59, 2311–2320]. We present a rigorous linear perturbation analysis for anisotropic materials, that allows the prediction of both the orientation of the corrugations of the thin film, and the wavelength that maximize the growth velocity. Finally, we compare our theoretical estimates to experimental results for a In_0.65Ga_0.35As alloy constraint to InP.     

Instability of a compressed elastic film on a viscous layer

A flat, compressed elastic film on a viscous layer is unstable. The film can form wrinkles to reduce the elastic energy. A linear perturbation analysis is performed to determine the critical wave number and the growth rate of the unstable modes. While the viscous layer has no effect on the critical wave number, its viscosity and thickness set the time scale for the growth of the perturbations. The fastest growing wave number and the corresponding growth rate are obtained as functions of the compressive strain and the thickness ratio between the viscous layer and the elastic film. The present analysis is valid for all thickness range of the viscous layer. In the limits of infinitely thick and thin viscous layers, the results reduce to those obtained in the previous studies.     

Equilibrium shapes of strained islands with non-zero contact angles

The equilibrium morphology of a strained island on an elastic substrate is determined. The island is assumed to partially wet the substrate (Volmer-Weber growth) and thus makes a non-zero contact angle with the surface. Both isotropic and anisotropic misfit strain are allowed. Two- and three-dimensional equilibrium island shapes are determined by using expressions for the elastic strain energy in the small-slope approximation. In this limit, the problem can be reduced to a singular integral-differential equation for the island thickness. We find that when there is a non-zero contact angle, all island shapes, for a given ratio of the elastic stress to surface energy, attain a form that is independent of the specific contact angle under an appropriate scaling. We show that for islands with non-zero contact angles, as the island volume increases, the shape approaches the geometry of a completely wetting island. But when the volume decreases, these islands approach a point while islands with a zero contact angle, approach a finite length line segment of zero volume. Multiple-hump equilibrium shapes are found. Single-humped islands are shown to have a lower chemical potential than multiple-humped islands, implying that they are the most stable. This conclusion is shown to be consistent with a stability analysis of the two-dimensional case. The effects of a tetragonal misfit strain on the three-dimensional island shape is investigated.     

Nonlinear analyses of wrinkles in a film bonded to a compliant substrate

Subject to a compressive membrane force, a film bonded to a compliant substrate often forms a pattern of wrinkles. This paper studies such wrinkles in a layered structure used in several recent experiments. The structure comprises a stiff film bonded to a compliant substrate, which in turn is bonded to a rigid support. Two types of analyses are performed. First, for sinusoidal wrinkles, by minimizing energy, we obtain the wavelength and the amplitude of the wrinkles for substrates of various moduli and thicknesses. Second, we develop a method to simultaneously evolve the two-dimensional pattern in the film and the three-dimensional elastic field in the substrate. The simulations show that the wrinkles can evolve into stripes, labyrinths, or herringbones, depending on the anisotropy of the membrane forces. Statistical averages of the amplitude and wavelength of wrinkles of various patterns correlate well with the analytical solution of the sinusoidal wrinkles.     

On the Development of Negative Pressure in Oil Film and the Characteristics of Journal Bearing

The stability and vibration characteristics of turbo-rotors running in journal bearings depend strongly on the characteristics of journal bearing determined by the oil film pressure and negative pressure developed in journal bearings. This paper is concerned with clarifying the development of this negative pressure and its influences by taking into account air bubbles in oil film, especially the surface dilatational viscosity, which decelerates the compression and expansion of air bubbles. Deceleration of the bubble expansion under negative pressure helps the bubbles to withstand greater negative pressure. When the negative pressure is developed in the oil film, the locus of center of journal running in bearing deviates from the well-known semi-circle like shape, and the journal center is pushed outwards and horizontally under vertical static load. This means that the stability of a turbo-rotor will be reduced when negative pressure is developed.     

Fracture mechanics analysis on Smart-Cut^R technology. Part 2： Effect of bonding flaws

In Part 2 of the paper on the Smart-Cut process, the effects of bonding flaws characterized by the size and internal pressure before and after splitting are studied by using fracture mechanics models. It is found that the bonding flaws with large size are prone to cause severe deviation of defect growth, leading to a non-transferred area of thin layer when splitting. In a practical Smart-Cut process where the internal pressure of bonding flaws is very small, large interfacial defects always promote defect growth in the splitting process. Meanwhile, increasing the internal pressure of the bonding flaws decreases the defect growth and its deviation before splitting. The mechanism of relaxation of stiffener constraint is proposed to clarify the effect of bonding flaws. Moreover, the progress of the splitting process is analyzed when bonding flaws are present. After splitting, those bonding flaws with large size and high internal pressure are vulnerable for the blistering of the thin film during high-temperature annealing.     

PROGRESSIVE FRACTURE MODELING OF THE FAILURE WAVE IN IMPACTED GLASS

The failure wave has been observed propagating in glass under impact loading since 1991. It is a continuous fracture zone which may be associated with the damage accumulation process during the propagation of shock waves. A progressive fracture model was proposed to describe the failure wave formation and propagation in shocked glass considering its heterogeneous meso-structures. The original and nucleated microcracks will expand along the pores and other defects with concomitant dilation when shock loading is below the Hugoniot Elastic Limit. The governing equation of the failure wave is characterized by inelastic bulk strain with material damage and fracture. And the inelastic bulk strain consists of dilatant strain from nucleation and expansion of microcracks and condensed strain from the collapse of the original pores. Numerical simulation of the free surface velocity was performed and found in good agreement with planar impact experiments on K9 glass at China Academy of Engineering Physics. And the longitudinal, lateral and shear stress histories upon the arrival of the failure wave were predicted, which present the diminished shear strength and lost spall strength in the failed layer.     

Experimental study of steady concentration fields in turbulent wakes

 The pollutant transportation process in turbulent wakes is studied experimentally using planar laser-induced fluorescence (PLIF). The concentration fields in the very near wake region behind typical bluff bodies are measured for steady flow. The characteristics of the mean and instantaneous concentration fields behind circular and sinusoidal islands and peninsulas are investigated. The results indicate that the pollutant distribution is closely related with the unsteady vortex shedding in the flow field. Compared with that of the circular island, more pollutants enter into the wake generated by the sinusoidal-shaped island. The time needed for pollutants to accumulate in or drain out of the wake region after the peninsula before reaching a relatively constant value is longer than that for the islands, regardless of the island or peninsula shape. The results will facilitate pollutant control behind sea islands and other natural or man-made structures in water. Also the results provide some fundamental data for checking numerical models. Received: 11 November 2000/Accepted: 26 January 2001     

BUCKLING ANALYSIS OF STRETCHABLE FERROELECTRIC THIN FILM ON ELASTOMERIC SUBSTRATES

The thin stiff films on pre-stretched compliant substrates can form wrinkles, which can be controlled in micro and nanoscale systems to generate smart structures. Recently, buck- led piezoelectric/ferroelectrie nanoribbons have been reported to show an enhancement in the piezoelectric effect and stretchability, which can be applied in energy harvesting devices, sensors and memory devices instead of polymeric polyvinylidine fluoride （PVDF）. This paper studies the buckling and post-buckling process of ferroelectric thin films bonded to the pre-stretched soft layer, which in turn lies on a rigid support. Nonlinear electromechanical equations for the buckling of thin piezoelectric plates are deduced and employed to model the ferroelectric film poled in the thickness direction. Two buckling modes are analyzed and discussed： partially de-adhered buck- ling and fully adhered buckling. Transition from one buckling mode to the other is predicted and the effect of piezoelectricity on the critical buckling condition of piezoelectric film is examined.     

高超声速溢流冷却实验研究

高超声速溢流冷却是一种新型的飞行器热防护方法,基本思想为:在高热流区布置溢流孔,控制冷却液以溢流方式流出,之后通过飞行器表面摩阻作用展布为液膜,形成热缓冲层以降低飞行器表面热流. 目前,溢流冷却技术还处于探索阶段,实现工程应用前还需开展大量的实验验证和机理研究工作. 本文首次开展溢流冷却的实验研究工作,采用热流测量、液膜厚度测量及液膜流动特性观测技术,搭建了完善的溢流冷却风洞实验平台,对溢流冷却热防护性能和高超声速条件下液膜流动规律进行了初步研究. 研究表明:(1) 高超声速流场中通过溢流能够在飞行器表面形成液膜并有效隔离外部高温气流,可降低飞行器表面热流率;(2) 楔面上的液膜前缘流动是一个逐渐减速的过程,增加冷却液流量液膜厚度变化不明显,但液膜前缘运动速度增大;(3) 液膜层存在表面波,在时间和空间方向发生演化,导致液膜厚度的微弱扰动;(4) 液膜层存在横向展宽现象,即液膜层宽度大于溢流缝宽度. 原因是液膜层与流场边界层条件不匹配,存在压力梯度,迫使冷却液向低压区流动,从而展宽液膜层,并且流量越高,横向展宽现象越明显.      

Temporal evolution and instability in a viscoelastic dielectric elastomer

Dielectric elastomer transducers are being developed for applications in stretchable electronics, tunable optics, biomedical devices, and soft machines. These transducers exhibit highly nonlinear electromechanical behavior: a dielectric membrane under voltage can form wrinkles, undergo snap-through instability, and suffer electrical breakdown. We investigate temporal evolution and instability by conducting a large set of experiments under various prestretches and loading rates, and by developing a model that allows viscoelastic instability. We use the model to classify types of instability, and map the experimental observations according to prestretches and loading rates. The model describes the entire set of experimental observations. A new type of instability is discovered, which we call wrinkle-to-wrinkle transition. A flat membrane at a critical voltage forms wrinkles and then, at a second critical voltage, snaps into another state of winkles of a shorter wavelength. This study demonstrates that viscoelasticity is essential to the understanding of temporal evolution and instability of dielectric elastomers.     

Ratcheting and wrinkling of tubes due to axial cycling under internal pressure: Part II analysis

Part I presented a set of experiments in which pressurized tubes were cycled axially under stress control about a compressive mean stress. This loading history causes biaxial ratcheting involving compressive axial strain and expansion of the diameter of the tube. The compressive strain in turn induces the initiation and growth of axisymmetric wrinkles. Persistent cycling resulted in localization of the wrinkles and collapse. In Part II the problem is first modeled as a shell with initial axisymmetric imperfections while the biaxial ratcheting of the material is modeled using the Dafalias–Popov two-surface nonlinear kinematic hardening model. It is demonstrated that when suitably calibrated this modeling framework reproduces the prevalent ratcheting deformations and the evolution of wrinkling including the conditions at collapse accurately for all experiments. The calibrated model is then used to evaluate the ratcheting behavior of pipes under thermal-pressure cyclic loading histories experienced by axially restrained pipelines.     

Heat and moisture transfer in gaps between sweating imitation skin and nonwoven cloth: effect of gap space and alignment of skin and clothing on the moisture transfer

This study investigates heat and moisture transfer between a sweating film and a nonwoven sheet both experimentally and numerically. A mathematical model based on heat conduction and moisture diffusion in both the air gap and cloth is presented. The evaporation rate and surface temperature of the sweating film are well predicted under various conditions such as air gap height, heating conditions, and sweating film orientation by evaluating the effective thermal conductivity and diffusion coefficient from the empirical equations of the Nusselt number for a fluid layer, even though the air gap height is sufficiently large to cause natural convections.     

石墨氧化物分子沉积膜的制备及其摩擦学行为研究

采用分子沉积技术制备了石墨氧化物（GO）－聚二丙烯二甲铵（PDDA）多层超薄膜，用紫外和红外光谱对超薄膜的结构进行了分析，用原了力显微镜考察了薄膜的表面形貌及纳米摩擦学行为。结果表明：石墨氧化物分子沉积膜可以降低玻璃表面的摩擦系数；其摩擦系数随载荷增大而降低；该分子沉积薄膜良好的润滑性能归因于其较低的表面剪切强度。     

Heat transfer in bubble layers at high pressures

An original experimental investigation of heat transfer with steam condensation on a surface of a horizontal cooled tube immersed in a bubbling layer was carried out. A copper test section 16 mm in diameter and 285 mm in length was placed in a bubbling column 295 mm in diameter. Experiments were made under a pressure of 0.72-3.8 MPa with volume steam content 0-0.18, steam superficial velocities 0-0.18 m/s, and liquid-wall temperature difference 38–106 K. The heat transfer process in a bubbling layer under high pressures is shown to be of considerably intensity; with moderate values of steam content heat transfer coefficients reach 10–12 kW/(m²·K). The use of the known correlations assumed for the case of air bubbling under atmospheric pressure results in systematically underestimating heat transfer by 30–80%. Data were obtained on heat transfer with film condensation of steam and natural convection of subcooled water at high temperature differences outside the range investigated earlier. Experimental data table is appended.     

Cooling of a viscoelastic film during unsteady extrusion from an annular die

The unsteady extrusion of a viscoelastic film from an annular and axisymmetric die is examined. External, elastic, viscous and inertia forces deform the film, which is simultaneously cooled via forced convection to the ambient air. This moving boundary problem is solved by mapping the liquid/air interfaces onto fixed ones and by employing a regular perturbation expansion for all the dependent variables. The ratio of the film thickness to its inner radius at the exit of the die is used as the small parameter in the perturbation expansion. The fluid mechanical aspects of the process depend on the Stokes, Deborah, Reynolds, and Capillary numbers. The heat transfer in the film and to the environment gives rise to four additional dimensionless groups: the Peclet, Biot and Brinkman numbers and the activation energy, which determines the temperature dependence of fluid viscosity and elasticity. A variable heat transfer coefficient is also considered. For typical fluid properties and process conditions, the Peclet number is very large. In this case it is the ratio of the Biot to the Peclet number, the Stanton number, which arises in the energy conservation equation. It is shown that film cooling becomes important when the Stanton number and/or the activation energy are in the high-end of their typical values. In such cases, the cooling of the parison leads to a more uniform flow and shape for the film. The influence on the process of a variable heat transfer coefficient and the Brinkman number is small. Received: 7 April 1999/Accepted: 10 August 1999