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.     

Vibration Properties of Thin Films of Transition 3d-Metals on the Cu (111) Surface

Russian Physics Journal - Using the potentials, generated by a tight binding approximation, the equilibrium atomic structure and the phonon spectra of Cr and Ni thin films (one monolayer) deposited...     

A new study of associating inhomogeneous fluids with classical density functional theory

ABSTRACT

A new density functional for the study of associating inhomogeneous fluids based on Wertheim's first-order thermodynamic perturbation theory is presented and compared to the most currently used associating density functionals. This functional is developed using the weighted density approximation in the range of association of hard spheres. We implement this functional within the framework of classical density functional theory together with modified fundamental measure theory to account for volume exclusion of hard spheres. This approach is tested against molecular simulations from literature of pure associating hard spheres and mixtures of non-associationg and associating hard spheres with different number of bonding sites close to a hard uniform wall. Furthermore, we compare and review our results with the performance of associating functionals from literature, one based on fundamental measure theory and the inhomogeneous version of Wertheim's perturbation theory. Results obtained with classical DFT and the three functionals show excellent agreement with molecular simulations in systems with one hard wall. For the cases of small pores where only one or two layers of fluid are allowed discrepancies between results with classical DFT and molecular simulations were found.     

Spectral estimates for saddle point matrices arising in weak constraint four‐dimensional variational data assimilation

We consider the large sparse symmetric linear systems of equations that arise in the solution of weak constraint four‐dimensional variational data assimilation, a method of high interest for numerical weather prediction. These systems can be written as saddle point systems with a 3 × 3 block structure but block eliminations can be performed to reduce them to saddle point systems with a 2 × 2 block structure, or further to symmetric positive definite systems. In this article, we analyse how sensitive the spectra of these matrices are to the number of observations of the underlying dynamical system. We also obtain bounds on the eigenvalues of the matrices. Numerical experiments are used to confirm the theoretical analysis and bounds.