Combined 3D thinning and greedy algorithm to approximate realistic particles with corrected mechanical properties

The shape of irregular particles has significant influence on micro- and macro-scopic behavior of granular systems. This paper presents a combined 3D thinning and greedy set-covering algorithm to approximate realistic particles with a clump of overlapping spheres for discrete element method (DEM) simulations. First, the particle medial surface (or surface skeleton), from which all candidate (maximal inscribed) spheres can be generated, is computed by the topological 3D thinning. Then, the clump generation procedure is converted into a greedy set-covering (SCP) problem. To correct the mass distribution due to highly overlapped spheres inside the clump, linear programming (LP) is used to adjust the density of each component sphere, such that the aggregate properties mass, center of mass and inertia tensor are identical or close enough to the prototypical particle. In order to find the optimal approximation accuracy (volume coverage: ratio of clump's volume to the original particle's volume), particle flow of 3 different shapes in a rotating drum are conducted. It was observed that the dynamic angle of repose starts to converge for all particle shapes at 85% volume coverage (spheres per clump < 30), which implies the possible optimal resolution to capture the mechanical behavior of the system.Comment: 34 pages, 13 figure

Generating Second Order (Co)homological Information within AT-Model Context

In this paper we design a new family of relations between (co)homology classes, working with coefficients in a field and starting from an AT-model (Algebraic Topological Model) AT(C) of a finite cell complex C These relations are induced by elementary relations of type “to be in the (co)boundary of” between cells. This high-order connectivity information is embedded into a graph-based representation model, called Second Order AT-Region-Incidence Graph (or AT-RIG) of C. This graph, having as nodes the different homology classes of C, is in turn, computed from two generalized abstract cell complexes, called primal and dual AT-segmentations of C. The respective cells of these two complexes are connected regions (set of cells) of the original cell complex C, which are specified by the integral operator of AT(C). In this work in progress, we successfully use this model (a) in experiments for discriminating topologically different 3D digital objects, having the same Euler characteristic and (b) in designing a parallel algorithm for computing potentially significant (co)homological information of 3D digital objects.Ministerio de Economía y Competitividad MTM2016-81030-PMinisterio de Economía y Competitividad TEC2012-37868-C04-0

Flow and retreat of the Late Quaternary Pine Island-Thwaites palaeo-ice stream, West Antarctica

Multibeam swath bathymetry and sub-bottom profiler data are used to establish constraints on the flow and retreat history of a major palaeo-ice stream that carried the combined discharge from the parts of the West Antarctic Ice Sheet now occupied by the Pine Island and Thwaites glacier basins. Sets of highly elongated bedforms show that, at the last glacial maximum, the route of the Pine Island-Thwaites palaeo-ice stream arced north-northeast following a prominent cross-shelf trough. In this area, the grounding line advanced to within similar to 68 km of, and probably reached, the shelf edge. Minimum ice thickness is estimated at 715 m on the outer shelf, and we estimate a minimum ice discharge of similar to 108 km(3) yr(-1) assuming velocities similar to today's Pine Island glacier (similar to 2.5 km yr(-1)). Additional bed forms observed in a trough northwest of Pine Island Bay likely formed via diachronous ice flows across the outer shelf and demonstrate switching ice stream behavior. The "style" of ice retreat is also evident in five grounding zone wedges, which suggest episodic deglaciation characterized by halts in grounding line migration up-trough. Stillstands occurred in association with changes in ice bed gradient, and phases of inferred rapid retreat correlate to higher bed slopes, supporting theoretical studies that show bed geometry as a control on ice margin recession. However, estimates that individual wedges could have formed within several centuries still imply a relatively rapid overall retreat. Our findings show that the ice stream channeled a substantial fraction of West Antarctica's discharge in the past, just as the Pine Island and Thwaites glaciers do today

Effect of the contraction ratio upon viscoelastic fluid flow in three-dimensional square-square contractions

In this work we investigate the laminar flow through square–square sudden contractions with various contraction ratios (CR¼2.4, 4,8and12), using a Newtonian fluid and a shear-thinning viscoelastic fluid. Visualizations of the flow patterns were carried out using streakline photography and detailed velocity field measurements were performed using particle image velocimetry. The experimental results are compared with numerical predictions obtained using a finite-volume method. For the Newtonian fluid, a corner vortex is found upstream of the contraction and increasing flow inertia leads to a reduction of the vortex size. Good agreement is observed between experiments and numerical simulations. For the shear-thinning fluid flow a corner vortex is also observed upstream of the contraction independently of the contraction ratio. Increasing the elasticity of the flow, while still maintaining low inertia flow conditions, leads to a strong increase of the vortex size, until an elastic instability sets in and the flow becomes time-dependent at DeE200, 300, 70 and 450 for CR¼2.4, 4, 8 and 12, respectively. At low contraction ratios, viscoelasticity brings out an anomalous divergent flow upstream of the contraction. For both fluids studied the flow presents a complex three-dimensional helical vortex structure which is well predicted by numerical simulations. However, for the viscoelastic fluid flow the maximum Deborah number achieved in the numerical simulations is about one order of magnitude lower than the critical Deborah number for the onset of the elastic instability found in the experiments