Bounding the stability and rupture condition of emulsion and foam films

A scaling law is presented that provides a complete solution to the equations bounding the stability and rupture of thin films. The scaling law depends on the fundamental physicochemical properties of the film and interface to calculate bounds for the critical thickness and other key film thicknesses, the relevant waveforms associated with instability and rupture, and film lifetimes. Critical thicknesses calculated from the scaling law are shown to bound the values reported in the literature for numerous emulsion and foam films. The majority of critical thickness values are between 15 to 40% lower than the upper bound critical thickness provided by the scaling law.

Scaling laws for the critical rupture thickness or common thin films

Despite decades of experimental and theoretical investigation on thin films, considerable uncertainty exists in the prediction of their critical rupture thickness. According to the spontaneous rupture mechanism, common thin films become unstable when capillary waves. at the interfaces begin to grow. In a horizontal film with symmetry at the midplane. unstable waves from adjacent interfaces grow towards the center of the film. As the film drains and becomes thinner, unstable waves osculate and cause the film to rupture, Uncertainty sterns from a number of sources including the theories used to predict film drainage and corrugation growth dynamics. In the early studies, (lie linear stability of small amplitude waves was investigated in the Context of the quasi-static approximation in which the dynamics of wave growth and film thinning are separated. The zeroth order wave growth equation of Vrij predicts faster wave growth rates than the first order equation derived by Sharma and Ruckenstein. It has been demonstrated in an accompanying paper that film drainage rates and times measured by numerous investigations are bounded by the predictions of the Reynolds equation and the more recent theory of Manev, Tsekov, and Radoev. Solutions to combinations of these equations yield simple scaling laws which should bound the critical rupture thickness of foam and emulsion films, In this paper, critical thickness measurements reported in the literature are compared to predictions from the bounding scaling equations and it is shown that the retarded Hamaker constants derived from approximate Lifshitz theory underestimate the critical thickness of foam and emulsion films, The non-retarded Hamaker constant more adequately bounds the critical thickness measurements over the entire range of film radii reported in the literature. This result reinforces observations made by other independent researchers that interfacial interactions in flexible liquid films are not adequately represented by the retarded Hamaker constant obtained from Lifshitz theory and that the interactions become significant at much greater separations than previously thought. (c) 2005 Elsevier B.V. All rights reserved.

Dynamic rupture in a 3-D particle-based simulation of a rough planar fault

An appreciation of the physical mechanisms which cause observed seismicity complexity is fundamental to the understanding of the temporal behaviour of faults and single slip events. Numerical simulation of fault slip can provide insights into fault processes by allowing exploration of parameter spaces which influence microscopic and macroscopic physics of processes which may lead towards an answer to those questions. Particle-based models such as the Lattice Solid Model have been used previously for the simulation of stick-slip dynamics of faults, although mainly in two dimensions. Recent increases in the power of computers and the ability to use the power of parallel computer systems have made it possible to extend particle-based fault simulations to three dimensions. In this paper a particle-based numerical model of a rough planar fault embedded between two elastic blocks in three dimensions is presented. A very simple friction law without any rate dependency and no spatial heterogeneity in the intrinsic coefficient of friction is used in the model. To simulate earthquake dynamics the model is sheared in a direction parallel to the fault plane with a constant velocity at the driving edges. Spontaneous slip occurs on the fault when the shear stress is large enough to overcome the frictional forces on the fault. Slip events with a wide range of event sizes are observed. Investigation of the temporal evolution and spatial distribution of slip during each event shows a high degree of variability between the events. In some of the larger events highly complex slip patterns are observed.