Gravitational aggregation regimes: critical dissipation threshold, optimal rigidity and fractal transition

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autore principale: Trivino, Yohann
Natura: Preprint
Pubblicazione: 2025
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866908694403350528
author Trivino, Yohann
author_facet Trivino, Yohann
contents I present a three-dimensional Discrete Element Method study of self-gravitation and contact mechanics in cold granular assemblies. The model couples direct Newtonian attraction between every particle pair with a linear visco-elastic normal contact law. Particles are treated as non-cohesive spheres; the normal force is parameterized to reproduce a prescribed restitution coefficient. Rotations are integrated using quaternions to avoid singularities. By normalizing the stiffness kn by kstar = G*m^2/R^3 and time by the free-fall time t_ff, I perform systematic parameter campaigns over dissipation (gamma) and normalized stiffness ktilde = kn/kstar. Results reveal three aggregation regimes. For low gamma the particles remain largely dispersive; above a critical gamma of about 5e2 aggregation accelerates until plateaus are reached in the aggregation time T_agg divided by t_ff. For stiffness ktilde on the order of 1e6 the aggregation time reaches a clear minimum. The cluster fraction C/Ntot shows a non-monotonic dependence on ktilde, with optimal cohesion at intermediate rigidity and peripheral isolation at extreme stiffness. Mapping the fractal dimension F across (gamma, ktilde) demonstrates transitions from compact structures (F about 3) to ramified structures (F below 2). These findings quantify how microscopic contact laws govern both the kinetics and microstructure of gravity-driven aggregation, providing a predictive framework for planetesimal formation and for calibrating DEM models against laboratory and micro-gravity experiments.
format Preprint
id arxiv_https___arxiv_org_abs_2512_05130
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Gravitational aggregation regimes: critical dissipation threshold, optimal rigidity and fractal transition
Trivino, Yohann
Soft Condensed Matter
Computational Physics
I present a three-dimensional Discrete Element Method study of self-gravitation and contact mechanics in cold granular assemblies. The model couples direct Newtonian attraction between every particle pair with a linear visco-elastic normal contact law. Particles are treated as non-cohesive spheres; the normal force is parameterized to reproduce a prescribed restitution coefficient. Rotations are integrated using quaternions to avoid singularities. By normalizing the stiffness kn by kstar = G*m^2/R^3 and time by the free-fall time t_ff, I perform systematic parameter campaigns over dissipation (gamma) and normalized stiffness ktilde = kn/kstar. Results reveal three aggregation regimes. For low gamma the particles remain largely dispersive; above a critical gamma of about 5e2 aggregation accelerates until plateaus are reached in the aggregation time T_agg divided by t_ff. For stiffness ktilde on the order of 1e6 the aggregation time reaches a clear minimum. The cluster fraction C/Ntot shows a non-monotonic dependence on ktilde, with optimal cohesion at intermediate rigidity and peripheral isolation at extreme stiffness. Mapping the fractal dimension F across (gamma, ktilde) demonstrates transitions from compact structures (F about 3) to ramified structures (F below 2). These findings quantify how microscopic contact laws govern both the kinetics and microstructure of gravity-driven aggregation, providing a predictive framework for planetesimal formation and for calibrating DEM models against laboratory and micro-gravity experiments.
title Gravitational aggregation regimes: critical dissipation threshold, optimal rigidity and fractal transition
topic Soft Condensed Matter
Computational Physics
url https://arxiv.org/abs/2512.05130