Mountain View
biomedical and healthCAre Data Discovery Index Ecosystem
help Advanced Search
Title: Contact changes in shear-stabilized jammed packings      
dateReleased:
09-03-2016
privacy:
information not avaiable
aggregation:
instance of dataset
dateCreated:
09-03-2016
refinement:
raw
ID:
doi:10.5281/ZENODO.59217
creators:
Dagois-Bohy, Simon
Henkes, Silke
Simon, Johannes
Tighe, Brian P.
Deen, Merlijn S. Van
Hecke, Martin Van
Zeravcic, Zorana
availability:
available
types:
other
description:
Authors are listed in alphabetical order. This data set contains the result of small simple shear deformations applied to approximately 140,000 shear-stabilized jammed packings (see [...]), focusing on contact changes, as described in [1,2,3,4]. These packings contain \(N = 16 \ldots 4096\) particles with harmonic interactions, under a confining pressure  \(p=10^{-7}\ldots10^{-2}\). Ensemble sizes range from 10 (N=4096) to 5000 (N=16).  In addition, a data file summarizing properties of the first contact change for each packing is provided. Particle interactions The simulation code minimizes the enthalpy \(H = \sum_{} \frac{k}{2} \delta_{ij}^2 + pL^2\) where L² is the simulation box area, p the externally applied pressure, k=1 the spring constant and  \(\delta_{ij} = \left\{ \begin{array}{ll} R_i + R_j - |\vec{r_{ij}}| & \textrm{if } |\vec{r_{ij}}| < R_i + R_j, \\ 0 & \textrm{otherwise.} \end{array}\right.\)   During shear, the boundary conditions are changed, and the system is relaxed to the new state. The simulation uses a bisection algorithm to efficiently step towards each subsequent contact change; see [3,4] for details. Data files The packings are stored in HDF5 data files. For each ensemble, we provide two files: one with and one without particle positions: N1024~P3162e-3_shear_noparticles.h5 includes all simulation data, but omits particle positions (see below for which data is included). For small data sets, N1024~P3162e-3_shear.h5 contains all simulations and all particle positions For large data sets, N1024~P3162e-3_shear_partial.h5 contains a subset of all simulations, but with all particle positions. Full particle positions for all simulations are available upon request to the authors. Please contact Martin van Hecke . All files follow the same HDF5 layout: Example name: N1024~P3162e-3_shear.h5 and N1024~P3162e-3_shear_noparticles.h5 Packings with \(N=1024\) particles Pressure \(p = 3.162\cdot 10^{-3}\) Packings are stored in a directory structure, e.g. /N1024/P3.1620e-03/0090/SR for the packing with id 0090. ​This directory contains a table 'data' indicating system parameters for each simulation step: boundary conditions L1 and L2 (also as L, alpha, delta) pressure P, strain gamma, stresses s_xy (simple shear), s_xx and s_yy,  number of contacts Ncontacts, contact number Z and number of rattlers #rattler number of changed contacts for this contact change bisection Nchanges, N+ (created), N- (broken) contact number Z energy U, enthalpy H and their change in the last relaxation step (dU, dH) step runtime t_run (seconds), #CG, #FIRE (number of conjugate gradient and FIRE iterations) path to the packing corresponding to this state (does not always exist for each state for older simulations) Each state is saved in /N1024/P3.1620e-03/0090/SR/0000 (initial), /N1024/P3.1620e-03/0090/SR/0001, ...etc. States are not included in the _noparticles.h5 files Some files omit intermediate positions, and only store positions just before and just after a contact change. The format of these directories is the same as in https://dx.doi.org/10.5281/zenodo.59216. Finally, we provide a summary file (shear_summary_cache.h5) which contains one table ('data') with properties of the first contact change of all packings. We provide the following columns: The variable postfix determines whether the value was calculated in the initial state (_base), just before the first contact change (_min) or just after the first contact change (_plus).   General/simulation properties Number of particles 'N' Random seed ['num', 'PackingNumber_base'] External pressure 'P0_base' Simulation step ['i_min', 'i_plus'] Relaxation statistics Last change in enthalpy during relaxation ['dH_base', 'dH_plus', 'dH_min'] Last change in energy during relaxation ['dU_base', 'dU_plus', 'dU_min'] Maximum gradient ['maxGrad_base', 'gg_min', 'gg_plus'] Initial simulation runtime ['runtime (s)_base'] State properties Number of rattlers 'N - Ncorrected_base' Number of non-rattler particles ['Neff_min', 'Neff_plus'] Number of contacts ['Ncontacts_plus', 'Ncontacts_min'] Contact number z ['Z_base', 'Z_min', 'Z_plus'] Internal pressure ['P', 'P_base', 'P_min', 'P_plus'] Mean overlap δ ['mean_delta_base'] Packing fraction ['phi_base', 'phi_min', 'phi_plus'] Enthalpy ['H_base', 'H_plus', 'H_min'] Energy ['Uhelper_base', 'U_min', 'U_plus'] Simple shear parameter alpha ['alpha_base', 'alpha_min', 'alpha_plus'] Pure shear parameter delta ['delta_base', 'delta_plus', 'delta_min'] Square root of area ['L_base', 'L_min', 'L_plus'] Stresses on boundaries: xx  ['sxx_base', 's_xx_min', 's_xx_plus',] yy  [ 'syy_base', 's_yy_plus', 's_yy_min',] xy ['sxy_base', 's_xy_plus', 's_xy_min'] Elastic moduli:  ['c1_base', 'c1_min', 'c1_plus', 'c2_base', 'c2_min', 'c2_plus', 'c3_base', 'c3_min', 'c3_plus', 'c4_base', 'c4_plus', 'c4_min', 'c5_base', 'c5_plus', 'c5_min', 'c6_base', 'c6_min', 'c6_plus', 'Dac_base', 'Dac_plus', 'Dac_min', AC component of G(θ) ['Gac_base', 'Gac_min', 'Gac_plus'] DC component of G(θ) ['Gdc_base', 'Gdc_min', 'Gdc_plus',] AC component of U(θ) ['Uac_base', 'Uac_plus', 'Uac_min'] DC component of U(θ) ['Udc_base', 'Udc_plus', 'Udc_min'] Simple shear ['Galpha_base', 'Galpha_plus', 'Galpha_min' ] Contact change properties Applied strain gamma ['gamma_plus', 'gamma_min'] gamma_min is used as contact change strain Number of created/broken contacts ['N+_plus', 'N+_min', 'N-_plus', 'N-_min'] Number of changed contacts (=N+ + N-) ['Nchanges_plus', 'Nchanges_min'] Making & breaking strain from upar and uperp: simple linear (SL) solution: ['gmk_SL_base', 'gbk_SL_base'] full quadratic (FQ) solution: ['gmk_FQ_base' 'gbk_FQ_base'] G up to CC from fit σ=Gγ & error bar ['Glin', 'Glinerr'] G up to CC from fit σ=Gγ + λγ² & error bar ['Gquad', 'Gquaderr'] λ up to CC from fit σ=Gγ + λγ² & error bar ['lambdaquad', 'lambdaquaderr'] [1] Simon Dagois-Bohy, Brian P. Tighe, Johannes Simon, Silke Henkes, and Martin van Hecke. Soft-Sphere Packings at Finite Pressure but Unstable to Shear. Phys. Rev. Lett. 109, 095703. http://dx.doi.org/10.1103/PhysRevLett.109.095703 [2] Merlijn S. van Deen, Johannes Simon, Zorana Zeravcic, Simon Dagois-Bohy, Brian P. Tighe, and Martin van Hecke. Contact changes near jamming. Phys. Rev. E 90 020202(R). http://dx.doi.org/10.1103/PhysRevE.90.020202 [3] Merlijn S. van Deen, Brian P. Tighe, and Martin van Hecke. Contact Changes of Sheared Systems: Scaling, Correlations, and Mechanisms. arXiv:1606.04799. https://arxiv.org/abs/1606.04799 [4] Merlijn S. van Deen. Mechanical Response of Foams: Elasticity, Plasticity, and Rearrangements. PhD Thesis, Leiden University, 2016. https://openaccess.leidenuniv.nl/handle/1887/40902
accessURL: https://doi.org/10.5281/ZENODO.59217
storedIn:
Zenodo
qualifier:
not compressed
format:
HTML
accessType:
landing page
authentication:
none
authorization:
none
abbreviation:
ZENODO
homePage: https://zenodo.org/
ID:
SCR:004129
name:
ZENODO

Feedback?

If you are having problems using our tools, or if you would just like to send us some feedback, please post your questions on GitHub.