A model for crack formation during active solid pyrolysis of a char-forming solid: crack patterns; surface area generation; volatile mass efflux

A model is developed for the formation and propagation of cracks in a material sample that is heated at its top surface, pyrolyses, and then thermally degrades to form char. In this work the sample is heated uniformly over its entire top surface by a hypothetical flame (a heat source). The pyrolysis mechanism is described by a one-step overall reaction that is dependent nonlinearly on the temperature (Arrhenius form). Stresses develop in response to the thermal degradation of the material by means of a shrinkage strain caused by local mass loss during pyrolysis. When the principal stress exceeds a prescribed threshold value, the material forms a local crack. Cracks are found to generally originate at the surface in response to heating, but occasionally they form in the bulk, away from ever-changing material boundaries. The resulting cracks evolve and form patterns whose characteristics are described. Quantities examined in detail are: the crack spacing in the pyrolysis zone; the crack length evolution; the formation and nature of crack loops which are defined as individual cracks that have joined to form loops that are disconnected from the remaining material; the formation of enhanced pyrolysis area; and the impact of all of the former quantities on mass flux. It is determined that the mass flux from the sample can be greatly enhanced over its nominal (non-cracking) counterpart. The mass efflux profile qualitatively resembles those observed in Cone Calorimeter tests.     

Second-order random wave solutions for interfacial internal waves in N-layer density-stratified fluid

This paper studies the random internal wave equations describing the density interface displacements and the velocity potentials of N-layer stratified fluid contained between two rigid walls at the top and bottom. The density interface displacements and the velocity potentials were solved to the second-order by an expansion approach used by Longuet-Higgins （1963） and Dean （1979） in the study of random surface waves and by Song （2004） in the study of second- order random wave solutions for internal waves in a two-layer fluid. The obtained results indicate that the first-order solutions are a linear superposition of many wave components with different amplitudes, wave numbers and frequencies, and that the amplitudes of first-order wave components with the same wave numbers and frequencies between the adjacent density interfaces are modulated by each other. They also show that the second-order solutions consist of two parts： the first one is the first-order solutions, and the second one is the solutions of the second-order asymptotic equations, which describe the second-order nonlinear modification and the second-order wave-wave interactions not only among the wave components on same density interfaces but also among the wave components between the adjacent density interfaces. Both the first-order and second-order solutions depend on the density and depth of each layer. It is also deduced that the results of the present work include those derived by Song （2004） for second-order random wave solutions for internal waves in a two-layer fluid as a particular case.     

Theory of interfacial phase transitions in surfactant systems

The spin-1 Ising model, which is equivalent to the three-component lattice gas model, is used to study wetting transitions in three-component surfactant systems consisting of an oil, water, and a nonionic surfactant. Phase equilibria, interfacial profiles, and interfacial tensions for three-phase equilibrium are determined in mean field approximation, for a wide range of temperature and interaction parameters. Surfactant interaction parameters are found to strongly influence interfacial tensions, reducing them in some cases to ultralow values. Interfacial tensions are used to determine whether the middle phase, rich in surfactant, wets or does not wet the interface between the oil-rich and water-rich phases. By varying temperature and interaction parameters, a wetting transition is located and found to be of the first order. Comparison is made with recent experimental results on wetting transitions in ternary surfactant systems.This paper is dedicated to J. K. Percus in honor of his 65th birthday.     

Hydrophobicity-induced drying transition in alkanethiol self-assembled monolayer—water interface

During the course of our investigation of the electron transfer properties of some redox species through highly hydrophobic long chain alkanethiol molecules on gold in aqueous and non-aqueous solvents, we obtained some intriguing results such as unusually low interfacial capacitance, very high values of impedance and film resistance, all of which pointed to the possible existence of a nanometer size interfacial gap between the hydrophobic monolayer and aqueous electrolyte. We explain this phenomenon by a model for the alkanethiol monolayer—aqueous electrolyte interface, in which the extremely hydrophobic alkanethiol film repels water molecules adjacent to it and in the process creates a shield between the monolayer film and water. This effectively increases the overall thickness of the dielectric layer that is manifested as an abnormally low value of interfacial capacitance. This behaviour is very much akin to the ‘drying transition’ proposed by Lum, Chandler and Weeks in their theory of length scale dependent hydrophobicity. For small hydrophobic units consisting of apolar solutes, the water molecules can reorganize around them without sacrificing their hydrogen bonds. Since for an extended hydrophobic unit, the existence of hydrogen bonded water structure close to it is geometrically unfavourable, there is a net depletion of water molecules in the vicinity leading to the possible creation of a hydrophobic interfacial gap.     

Closed-form solutions for a rectangular layer with central crack under anti-plane shear stress

We consider the problem of determining the stress distributionin a finite rectangular elastic layer containing a Griffithcrack which is opened by internal shear stress acting alongthe length of the crack. The mode III crack is assumed to belocated in the middle plane of the rectangular layer. The followingtwo problems are considered: (A) the central crack is perpendicularto the two fixed lateral surfaces and parallel to the othertwo stress-free surfaces; (B) all the lateral surfaces of therectangular layer are clamped and the central crack is parallelto the two lateral surfaces. By using Fourier transformations,we reduce the solution of each problem to the solution of dualintegral equations with sine kernels and a weight function whichare solved exactly. Finally, we derive closed-form expressionsfor the stress intensity factor at the tip of the crack andthe numerical values for the stress intensity factor at theedges of the cracks are presented in the form of tables.     

In-Situ Intrinsic Interface Strength Measurement at Elevated Temperatures and its Relationship to Interfacial Structure

A previously developed laser spallation technique has been modified to measure the tensile strength of thin film interfaces in-situ at temperatures up to 1100°C. Tensile strengths of Nb/A-plane sapphire, FeCrAl/A-plane sapphire and FeCrAlY/A-plane sapphire were measured up to 950°C. The measured strengths at high temperatures were substantially lower compared with their corresponding strengths at ambient temperature. For example, at 850°C, the interface tensile strength for the Nb/sapphire (151 ± 17 MPa), FeCrAl/sapphire (62 ± 8 MPa) and FeCrAlY/sapphire (82 ± 11 MPa) interface systems were lower by factors of approximately, 3, 5, and 8, respectively, over their corresponding ambient values. These results underscore the importance of using such in-situ measured values under operating conditions as the failure criterion in any life prediction or reliability models of such coated systems where local interface temperature excursions are expected. The results on alloy film interfaces also demonstrate that the presence of Y increases the strength of FeCrAl/Al₂O₃ interfaces.     

Stress State of Multilayer Thin-Walled Elements with Interfacial Defects of Structure

Based on the discrete-structural theory of thin plates and shells, a calculation model for thin-walled elements consisting of a number of rigid anisotropic layers is put forward. It is assumed that the transverse shear and compression stresses are equal on the interfaces. Elastic slippage is allowed over the interfaces between adjacent layers. The solution to the problem is obtained in a geometrically nonlinear statement with account of the influence of transverse shear and compression strains. The stress-strain state of circular two-layer transversely isotropic plates, both without defects and with a local area of adhesion failure at their center, is investigated numerically and experimentally. It is found that the kinematic and static contact conditions on the interfaces of layered thin-walled structural members greatly affect the magnitude of stresses and strains. With the use of three variants of calculation models, in the cases of perfect and weakened contact conditions between layers, the calculation results for circular plates are compared. It is revealed that the variant suggested in this paper adequately reflects the behavior of layered thin-walled structural elements under large deformations. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 41, No. 6, pp. 761–772, November–December, 2005.     

Shape sensitivity of a plane crack front

The 3D‐elasticity model of a solid with a plane crack under the stress‐free boundary conditions at the crack is considered. We investigate variations of a solution and of energy functionals with respect to perturbations of the crack front in the plane. The corresponding expansions at least up to the second‐order terms are obtained. The strong derivatives of the solution are constructed as an iterative solution of the same elasticity problem with specified right‐hand sides. Using the expansion of the potential and surface energy, we consider an approximate quadratic form for local shape optimization of the crack front defined by the Griffith criterion. To specify its properties, a procedure of discrete optimization is proposed, which reduces to a matrix variational inequality. At least for a small load we prove its solvability and find a quasi‐static model of the crack growth depending on the loading parameter. Copyright © 2003 John Wiley & Sons, Ltd.     

Evaluation of the Lazarus-Leblond constants in the asymptotic model of the interfacial wavy crack

The paper addresses the problem of a semi-infinite plane crack along the interface between two isotropic half-spaces. Two methods of solution have been considered in the past: Lazarus and Leblond [1998a. Three-dimensional crack-face weight functions for the semi-infinite interface crack-I: variation of the stress intensity factors due to some small perturbation of the crack front. J. Mech. Phys. Solids 46, 489-511, 1998b. Three-dimensional crack-face weight functions for the semi-infinite interface crack-II: integrodifferential equations on the weight functions and resolution J. Mech. Phys. Solids 46, 513-536] applied the “special” method by Bueckner [1987. Weight functions and fundamental fields for the penny-shaped and the half-plane crack in three space. Int. J. Solids Struct. 23, 57-93] and found the expression of the variation of the stress intensity factors for a wavy crack without solving the complete elasticity problem; their solution is expressed in terms of the physical variables, and it involves five constants whose analytical representation was unknown; on the other hand, the “general” solution to the problem has been recently addressed by Bercial-Velez et al. [2005. High-order asymptotics and perturbation problems for 3D interfacial cracks. J. Mech. Phys. Solids 53, 1128-1162], using a Wiener-Hopf analysis and singular asymptotics near the crack front.The main goal of the present paper is to complete the solution to the problem by providing the connection between the two methods. This is done by constructing an integral representation for Lazarus-Leblond's weight functions and by deriving the closed form representations of Lazarus-Leblond's constants.     

Front roughening in three-dimensional imbibition

We investigate the structure and dynamics of the interface between two immiscible liquids in a three-dimensional disordered porous medium. We apply a phase-field model that includes explicitly disorder and discuss both spontaneous and forced imbibition. The structure of the interface is dominated by a length scale ξ_× which arises from liquid conservation. We further show that disorder in the capillary and permeability act on different length scales and give rise to different scalings and structures of the interface properties. We conclude with a range of applications.