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.     

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.     

Wave flow of a liquid film together with a flow of gas

The wave flow of a thin layer of viscous liquid in conjunction with a flow of gas was considered in a linear formulation earlier [1, 2]. In this paper the problem of the wave flow of a liquid film together with a gas flow is solved in a nonlinear setting. On this basis relationships are derived for calculating the parameters of the film and the hydrodynamic quantities.Ivanovo. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 12–18, January–February, 1972.     

Waves on the surface of a thin layer of viscous liquid

Under the assumption that the boundary layer approximation for the original equations is valid, we show the possibility of the existence of progressive waves on the surface of a vertically flowing film when surface tension is neglected. From the system of equations obtained for a thin layer of viscous liquid flowing down an inclined plane, one equation for perturbations of a thin film follows. Steady solutions of this equation allow periodic discontinuous solutions of the roll-wave type.Translated from Zhurnal Prikladnoi Mekhaniki i Teknicheskoi Fiziki, No. 2, pp. 109–113-March–April, 1973.     

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.     

Stability of flow of a thin layer of viscous magnetic liquid

The laminar flow of a thin layer of heavy viscous magnetic liquid down an inclined wall is examined. The stability and control of the flow of an ordinary liquid are affected only by alteration of the angle of inclination of the solid wall and the velocity of the adjacent gas flow. When magnetic liquids are used [1, 2], an effective method of flow control may be control of the magnetic field. By using magnetic fields of various configurations it is possible to control the flow of a thin film of viscous liquid, modify the stability of laminar film flow, and change the shape of the free surface of the laminarly flowing thin film, a factor which plays a role in mass transfer, whose rate depends on the phase contact surface area. The magnetic field significantly affects the shape of the free surface of a magnetic liquid [3, 4]. In this paper the velocity profile of a layer of viscous magnetic liquid adjoining a gas flow and flowing down an inclined solid wall in a uniform magnetic field is found. It is shown that the flow can be controlled by the magnetic field. The problem of stability of the flow is solved in a linear formulation in which perturbations of the magnetic field are taken into account. The stability condition is found. The flow stability is affected by the nonuniform nature of the field and also by its direction.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 59–65, September–October, 1977.     

Wrinkling of a charged elastic film on a viscous layer

A thin metallic film deposited on a compliant polymeric substrate begins to wrinkle under compression induced in curing process and afterwards cooling of the system. The wrinkle mode depends upon the thin film elasticity, thickness, compressive strain, as well as mechanical properties of the compliant substrate. This paper presents a simple model to study the modulation of the wrinkle mode of thin metallic films bonded on viscous layers in external electric field. During the procedure, linear perturbation analysis was performed for determining the characteristic relation that governs the evolution of the plane-strain wrinkle of the thin films under varying conditions, i.e., the maximally unstable wrinkle mode as a function of the film surface charge, film elasticity and thickness, misfit strain, as well as thickness and viscosity of the viscous layer. It shows that, in proper electric field, thin film may wrinkle subjected to either compression or tension. Therefore, external electric field can be employed to modulate the wrinkle mode of thin films. The present results can be used as the theoretical basis for wrinkling analysis and mode modulation in surface metallic coatings, drying adhesives and paints, and microelectromechanical systems (MEMS), etc.     

Three-dimensional chemically nonequilibrium viscous shock layer on a catalytic surface with allowance for associated heat transfer

A three-dimensional flow of dissociating air past blunt bodies is investigated in the framework of the thin viscous shock layer theory. Multicomponent diffusion and homogeneous chemical reactions, including dissociation, recombination, and exchange reactions, are taken into account. The generalized Rankine-Kugoniot conditions are specified on the shock wave and the conditions which take into account the heterogeneous catalytic reactions, on the surface of the body. The viscous shock layer equations are solved together with the heat equations inside the coating, which is carbon with a deposited thin film of SiO₂, or quartz. The case of a thermally insulated surface is also considered. The problem for the case of the motion of a body along the re-entry trajectory into Earth's atmosphere is investigated numerically. The temperature of the surface and the heat flux toward it are given as a dependence on the height (tine) of the flight for different cases of the specification of the catalytic reactions. It is shown that the difference between the heat fluxes towards the thermally insulated surface and the fluxes toward the heat-conducting surface in the neighborhood of the stagnation point is of the order of 6–12% for all the cases considered. This makes it possible to decouple the solution of the problem of heat conduction in the body.Translated fron Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 140–146, November–December, 1985.deceased     

Investigation of the hypersonic viscous shock layer around blunt bodies in a nonuniform flow

Unseparated viscous gas flow past a body is numerically investigated within the framework of the theory of a thin viscous shock layer [13–15]. The equations of the hypersonic viscous shock layer with generalized Rankine-Hugoniot conditions at the shock wave are solved by a finite-difference method [16] over a broad interval of Reynolds numbers and values of the temperature factor and nonuniformity parameters. Calculation results characterizing the effect of free-stream nonuniformity on the velocity and temperature profiles across the shock layer, the friction and heat transfer coefficients and the shock wave standoff distance are presented. The unseparated flow conditions are investigated and the critical values of the nonuniformity parameter a_k [10] at which reverse-circulatory zones develop on the front of the body are obtained as a function of the Reynolds number. The calculations are compared with the asymptotic solutions [10, 12].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 154–159, May–June, 1987.     

Mechanics of relaxing SiGe islands on a viscous glass

A process has been developed recently to fabricate a structure comprising, from top to bottom, a SiGe thin film, a glass layer, and a Si wafer. The SiGe film is a perfect crystal, and is under biaxial compression. The SiGe film is patterned into islands. On annealing, the glass flows and the islands relax. The resulting strain-free islands are used as substrates, to grow epitaxial optoelectronic devices. This article describes a series of studies on the annealing process, combining experiment and theory. A small island relaxes by expansion, starting at the edges and diffusing to the center. A large island wrinkles before the expansion reaches the center. After some time, the wrinkles either disappear, or cause the island to fracture. We model the island as an elastic plate, and the glass layer as a viscous liquid. The strains in the islands are measured by X-ray diffraction and Raman spectroscopy, and the wrinkle amplitudes by atomic force microscope. The data are compared with the theoretical predictions. We determine the conditions under which the islands relax by expansion without significant wrinking, and demonstrate that a cap layer suppresses wrinkles, relaxing a large island crack-free. The project supported by NSF (CMS-9820713), DARPA (N66001-00-1-8957), ARO (DAA655-98-1-0270), and New Jersey Science and Technology Commission     

Forward roll coating flows of viscoelastic liquids

Roll coating is distinguished by the use of one or more gaps between rotating cylinders to meter and apply a liquid layer to a substrate. Except at low speed, the two-dimensional film splitting flow that occurs in forward roll coating is unstable; a three-dimensional steady flow sets in, resulting in more or less regular stripes in the machine direction. For Newtonian liquids, the stability of the two-dimensional flow is determined by the competition of capillary and viscous forces: the onset of meniscus nonuniformity is marked by a critical value of the capillary number. Although most of the liquids coated industrially are non-Newtonian polymeric solutions and dispersions, most of the theoretical analyses of film splitting flows relied on the Newtonian model. Non-Newtonian behavior can drastically change the nature of the flow near the free surface; when minute amounts of flexible polymer are present, the onset of the three-dimensional instability occurs at much lower speeds than in the Newtonian case.Forward roll coating flow is analyzed here with two differential constitutive models, the Oldroyd-B and the FENE-P equations. The results show that the elastic stresses change the flow near the film splitting meniscus by reducing and eventually eliminating the recirculation present at low capillary number. When the recirculation disappears, the difference of the tangential and normal stresses (i.e., the hoop stress) at the free surface becomes positive and grows dramatically with fluid elasticity, which explains how viscoelasticity destabilizes the flow in terms of the analysis of Graham [M.D. Graham, Interfacial hoop stress and instability of viscoelastic free surface flows, Phys. Fluids 15 (2003) 1702–1710].     

Instability of the motion of a liquid between two rotating spherical surfaces

A study has been made of the stability of the steady-state motion of a viscous incompressible liquid, arising in a thin spherical layer, when both spheres are rotating in the same direction at different angular velocities. For a ratio of the radii of the spheres r₂/r₁=1.10, 1.07, a stability curve is obtained which is analogous to the stability curve for the motion of a liquid between rotating cylinders.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 155–156, July–August, 1970.     

Buckling and postbuckling of a compressed thin film bonded on a soft elastic layer: a three-dimensional analysis

The wrinkling of a stiff thin film bonded on a soft elastic layer and subjected to an applied or residual compressive stress is investigated in the present paper. A three-dimensional theoretical model is presented to predict the buckling and postbuckling behavior of the film. We obtained the analytical solutions for the critical buckling condition and the postbuckling morphology of the film. The effects of the thicknesses and elastic properties of the film and the soft layer on the characteristic wrinkling wavelength are examined. It is found that the critical wrinkling condition of the thin film is sensitive to the compressibility and thickness of the soft layer, and its wrinkling amplitude depends on the magnitude of the applied or residual in-plane stress. The bonding condition between the soft layer and the rigid substrate has a considerable influence on the buckling of the thin film, and the relative sliding at the interface tends to destabilize the system.     

PROPAGATION OF LOVE WAVES IN PRESTRESSED PIEZOELECTRIC LAYERED STRUCTURES LOADED WITH VISCOUS LIQUID

We investigate analytically the effect of initial stress in piezoelectric layered structures loaded with viscous liquid on the dispersive and attenuated characteristics of Love waves, which involves a thin piezoelectric layer bonded perfectly to an unbounded elastic substrate. The effects of initial stress in the piezoelectric layer and the viscous coefficient of the liquid on the phase velocity of Love waves are analyzed. Numerical results are presented and discussed. The analytical method and the results can be useful for the design of chemical and biosensing liquid sensors.     

Effect of dopants and interlayers on the growth of thin insulating films

Surface morphologies of thin dielectric films deposited on gallium arsenide substrates are studied by atomic force microscopy (AFM). The quasi-periodic mesostructure with a corrugated configuration is found to form during the deposition process. A special dopant and thin interlayer at the film–substrate interface are used to decrease the surface roughness. The corrugated Si_xN_yO_z–SiO₂ film surface disappears by introducing Se atoms into the subsurface layer of the semiconductor. The root-mean-square roughness and the fractal dimension techniques are used for the numerical characterization of the surface morphologies of thin insulator films.     

Buckling of a stiff film bound to a compliant substrate—Part I:: Formulation, linear stability of cylindrical patterns, secondary bifurcations

The buckling of a thin elastic film bound to a compliant substrate is studied: we analyze the different patterns that arise as a function of the biaxial residual compressive stress in the film. We first clarify the boundary conditions to be used at the interface between film and substrate. We carry out the linear stability analysis of the classical pattern made of straight stripes, and point out secondary instabilities leading to the formation of undulating stripes, varicose, checkerboard or hexagonal patterns. Straight stripes are found to be stable in a narrow window of load parameters only. We present a weakly nonlinear post-buckling analysis of these patterns: for equi-biaxial residual compression, straight wrinkles are never stable and square checkerboard patterns are found to be optimal just above threshold; for anisotropic residual compression, straight wrinkles are present above a primary threshold and soon become unstable with respect to undulating stripes. These results account for many of the previously published experimental or numerical results on this geometry.     

Thin film/substrate systems featuring arbitrary film thickness and misfit strain distributions. Part I: Analysis for obtaining film stress from non-local curvature information

Current methodologies used for the inference of thin film stress through curvature measurements are strictly restricted to stress and curvature states which are assumed to remain uniform over the entire film/substrate system. Recently Huang, Rosakis and co-workers [Huang, Y., Ngo, D., Rosakis, A.J., 2005. Non-uniform, axisymmetric misfit strain: in thin films bonded on plate substrates/substrate systems: the relation between non-uniform film stresses and system curvatures. Acta Mech. Sin. 21, 362–370; Huang, Y., Rosakis A.J., 2005. Extension of Stoney’s Formula to non-uniform temperature distributions in thin film/substrate systems. The case of radial symmetry. J. Mech. Phys. Solids 53, 2483–2500; Ngo, D., Huang, Y., Rosakis, A. J., Feng, X. 2006. Spatially non-uniform, isotropic misfit strain in thin films bonded on plate substrates: the relation between non-uniform film stresses and system curvatures. Thin Solid Films (in press)] established methods for film/substrate system subject to non-uniform misfit strain and temperature changes. The film stresses were found to depend non-locally on system curvatures (i.e., depend on the full-field curvatures). The existing methods, however, all assume uniform film thickness which is often violated in the thin film/substrate system. We extend these methods to arbitrarily non-uniform film thickness for the thin film/substrate system subject to non-uniform misfit strain. Remarkably the stress-curvature relation for uniform film thickness still holds if the film thickness is replaced by its local value at the point where the stress is evaluated. This result has been experimentally validated in Part II of this paper.     

Unsteady thermocapillary convection in a nonuniformly heated fluid layer

In many technological processes, thin extended layers of nonuniformly heated fluid are used [1–3]. If they are sufficiently thin, thermocapillary forces have a decisive influence on the occurrence and development of motion of the fluid [4–6]. Investigation of convective motion in such a layer is of great interest for estimating the intensity of heat and mass transfer in technological processes. This paper is a study of unsteady thermocapillary motion in a layer of viscous incompressible fluid with free surface in which a thermal inhomogeneity is created at the initial time. Approximate expressions are obtained for the fields of the velocity, temperature, and pressure in the fluid, and also for the shape of the free surface.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 17–25, May–June, 1991.     

Interaction of an elastic boundary with a viscous sublayer of a turbulent boundary layer

The boundary region of a turbulent boundary layer contributes greatly to the drag. Intense turbulence is generated in this region. Below we investigate the interaction of an elastic boundary with a viscous sublayer for a decrease in the Reynolds stresses, and for a corresponding decrease in the drag. It does not seem possible to investigate the general case. Therefore, the problem is solved within the framework of the limitations made by Sternberg [1] for the theory of a viscous sublayer in a turbulent flow near a solid smooth wall.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 58–62, May–June, 1971.The authors thanks G. S. Migirenko for advice and remarks given during a discussion of the work.     

Model of a viscous nonequilibrium shock layer with a thin shock wave

It is shown that the concept of a viscous shock layer with boundary conditions specified in a thin shock wave is unsuitable for analyzing the flow of a chemically reacting gas, even in the case of high Reynolds numbers; it may produce a finite error when determining the parameters of the shock layer.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 175–178, September–October, 1973.