Effective, anisotropic elasticity tensor of snow, firn, and bubbly ice

HS bounds have the advantage of the correct limiting behavior for low to high-ice volume fractions. We derived parameterization for five transversely isotropic elasticity tensor components, requiring only two free parameters. This parameterization was valid for ice volume fractions ranging from 0.06 to 0.93.

The analysis employing the Thomsen parameter highlighted the dominance of geometrical anisotropy until an ice volume fraction of 0.7. However, to fully comprehend the elasticity of bubbly ice, a comprehensive approach is necessary to integrate coupled elastic theories that account for both geometrical and crystallographic anisotropy. This dataset includes a Jupyter notebook with all the necessary functions required to predict the elasticity tensor of snow for the given ice volume fraction and anisotropy.

Also, the code contains the least squares optimization function to compute the elasticity tensor for the six components of stress and strain. For example, we consider our dataset to calculate the samples' elasticity tensor and reproduce Fig. 7 from the paper. We take the stress and strain values obtained from load states as input for this example.

Also, a .csv file contains all the microstructural information: ice volume fraction, anisotropy, correlation functions, voxels size, and no. of voxels of the samples and the elasticity tensor obtained from finite element simulations and from present work parameterization.