Linking snow microstructure to its macroscopic elastic stiffness tensor: A numerical homogenization method and its application to 3‐D images from X‐ray tomography

The full 3D macroscopic mechanical behavior of snow is investigated by solving Kinematically Uniform Boundary Condition (KUBC) problems derived from homogenization theories over 3D images obtained by X-ray tomography. Snow is modeled as a porous cohesive material and its mechanical stiffness tensor is computed within the framework of the elastic behavior of ice. The size of the optimal representative elementary volume, expressed in terms of correlation lengths, is determined through a convergence analysis of the computed effective properties. A wide range of snow densities is explored and power laws with high regression coefficients are proposed to link the Young's and shear moduli of snow to its density. The degree of anisotropy of these properties is quantified and Poisson's ratios are also provided. Finally, the influence of the main types of metamorphism (isothermal, temperature gradient and wet snow metamorphism) on the elastic properties of snow and on their anisotropy is reported.

[1]  M. Schneebeli.  Numerical simulation of elastic stress in the microstructure of snow , 2004, Annals of Glaciology.

[2]  I. Baker,et al.  Evolution of individual snowflakes during metamorphism , 2010 .

[3]  Christoph H. Arns,et al.  Elastic and transport properties of cellular solids derived from three-dimensional tomographic images , 2006, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[4]  P. K. Srivastava,et al.  Micromechanical analysis of deformation of snow using X-ray tomography , 2014 .

[5]  C. Lantuéjoul,et al.  Ergodicity and integral range , 1991 .

[6]  R. Zimmerman.  Elastic moduli of a solid containing spherical inclusions , 1991 .

[7]  Cécile Coléou,et al.  From snow X-ray microtomograph raw volume data to micromechanics modeling: first results , 2004, Annals of Glaciology.

[8]  S. Morin,et al.  Numerical and experimental investigations of the effective thermal conductivity of snow , 2011 .

[9]  J. Schweizer,et al.  The temperature dependence of the effective elastic shear modulus of snow , 2002 .

[10]  Frédéric Flin,et al.  The temperature-gradient metamorphism of snow: vapour diffusion model and application to tomographic images , 2008, Annals of Glaciology.

[11]  S. Lejeunes,et al.  Une Toolbox Abaqus pour le calcul de propriétés effectives de milieux hétérogènes , 2011 .

[12]  E. Brun,et al.  A numerical model to simulate snow-cover stratigraphy for operational avalanche forecasting , 1992, Journal of Glaciology.

[13]  R. Hill.  Elastic properties of reinforced solids: some theoretical principles , 1963 .

[14]  H. Löwe,et al.  Hot-pressure sintering of low-density snow analyzed by X-ray microtomography and in situ microcompression , 2014 .

[15]  M. Schneebeli,et al.  Numerical simulation of microstructural damage and tensile strength of snow , 2014 .

[16]  S. R. D. Roscoat,et al.  Study of a temperature gradient metamorphism of snow from 3D images: time evolution of microstructures, physical properties and their associated anisotropy , 2013 .

[17]  D. Jeulin,et al.  Determination of the size of the representative volume element for random composites: statistical and numerical approach , 2003 .

[18]  Christian Geindreau,et al.  3-D image-based numerical computations of snow permeability: links to specific surface area, density, and microstructural anisotropy , 2012 .

[19]  M. Schneebeli,et al.  Three-dimensional microstructure and numerical calculation of elastic properties of alpine snow with a focus on weak layers , 2014 .

[20]  Martin Schneebeli,et al.  Vapor flux and recrystallization during dry snow metamorphism under a steady temperature gradient as observed by time-lapse micro-tomography , 2012 .

[21]  K. Tanaka,et al.  Average stress in matrix and average elastic energy of materials with misfitting inclusions , 1973 .

[22]  P. K. Srivastava,et al.  Observation of temperature gradient metamorphism in snow by X-ray computed microtomography: measurement of microstructure parameters and simulation of linear elastic properties , 2010, Annals of Glaciology.

[23]  B. Lesaffre,et al.  CellDyM: A room temperature operating cryogenic cell for the dynamic monitoring of snow metamorphism by time‐lapse X‐ray microtomography , 2015 .

[24]  P. Bartelt,et al.  A physical SNOWPACK model for the Swiss avalanche warning: Part I: numerical model , 2002 .

[25]  Cécile Coléou,et al.  Three-dimensional geometric measurements of snow microstructural evolution under isothermal conditions , 2004, Annals of Glaciology.

[26]  David A. Boas,et al.  Tetrahedral mesh generation from volumetric binary and grayscale images , 2009, 2009 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[27]  Martin Schneebeli,et al.  A general treatment of snow microstructure exemplified by an improved relation for the thermal conductivity , 2012 .

[28]  P. Gumbsch,et al.  Anticrack Nucleation as Triggering Mechanism for Snow Slab Avalanches , 2008, Science.

[29]  R. Huiskes,et al.  Direct mechanics assessment of elastic symmetries and properties of trabecular bone architecture. , 1996, Journal of biomechanics.

[30]  S. Shtrikman,et al.  A variational approach to the theory of the elastic behaviour of multiphase materials , 1963 .

[31]  Sergey A. Sokratov,et al.  A microstructural approach to model heat transfer in snow , 2005 .

[32]  F. Flin,et al.  Macroscopic modeling for heat and water vapor transfer in dry snow by homogenization. , 2014, The journal of physical chemistry. B.

[33]  Edward E. Adams,et al.  FINE STRUCTURE LAYERING IN RADIATION RECRYSTALLIZED SNOW , 2014 .

[34]  J. Schweizer,et al.  Influence of weak layer heterogeneity and slab properties on slab tensile failure propensity and avalanche release area , 2014 .