Search

We investigate the dynamics of two-dimensional assemblies of rod-shaped magnetic colloids under the influence of an external rotating magnetic field. Using Molecular Dynamics, we simulate the formation of patterns that emerge based on the synchronization degree between the magnetic rods and the rotating field. We then explore the structural and dynamic characteristics of the resulting steady states, examining their evolution as a function of changes in the rods' aspect ratio, the strength of the external magnetic field, and its rotation frequency. Three distinct synchronization regimes of the rods with the magnetic field are clearly observed. A detailed set of phase diagrams illustrates the complex relationship between the magnitude of the external magnetic field and its rotation frequency and how these parameters govern the formation of unique self-organized structures., Comment: 10 pages, 9 figures

Stressed soft materials commonly present viscoelastic signatures in the form of power-law or exponential decay. Understanding the origins of such rheologic behaviors is crucial to find proper technological applications. Using an elastic network model of macromolecules immersed in a viscous fluid, we numerically reproduce those characteristic viscoelastic relaxations and show how the microscopic interactions determine the rheologic response. We find that exponential relaxations are indeed the most common behavior. However, power laws may arise when drag forces between the macromolecules and the fluid are sublinear, which is related to micro-deformations of the macromolecules.

Rheological responses are the most relevant features to describe soft matter. So far, such constitutive relations are still not well understood in terms of small scale properties, although this knowledge would help the design of synthetic and bio-materials. Here, we investigate, computational and analytically, how mesoscopic-scale interactions influence the macroscopic behavior of viscoelastic materials. We design a coarse-grained approach where the local elastic and viscous contributions can be controlled. Applying molecular dynamics simulations, we mimic real indentation assays. When elastic forces are dominant, our model reproduces the hertzian behavior. However, when friction increases, it restores the Standard Linear Solid model. We show how the response parameters depend on the microscopic elastic and viscous contributions. Moreover, our findings also suggest that the relaxation times, obtained in relaxation and oscillatory experiments, obey a universal behavior in viscoelastic materials.

We study the contact process on spatially embedded networks, consisting of a regular square lattice with long-range connections. To generate the networks, a long-range connection is randomly added to each node $i$ of a square lattice, following the probability, $P_{ij}\sim{r_{ij}^{-\alpha}}$ , where $r_{ij}$ is the Manhattan distance between nodes $i$ and $j$, and the exponent $\alpha$ is a tunable parameter. Extensive Monte Carlo simulations and a finite-size scaling analysis for different values of $\alpha$ reveal a crossover from the mean-field to $2d$ Directed Percolation universality class with increasing $\alpha$, in the range $3<\alpha<4$., Comment: 5 pages, 7 figures

We propose a modified Invasion Percolation (IP) model to simulate the infiltration of glue into a porous medium under gravity in 2D. Initially, the medium is saturated with air and then invaded by a fluid that has a hardening effect taking place from the interface towards the interior by contact with the air. To take into account that interfacial hardening, we use an IP model where capillary pressures of the growth sites are increased with time. In our model, if a site stays for a certain time at interface, it becomes a hard site and cannot be invaded anymore. That represents the glue interface becoming hard due to exposition with the air. Buoyancy forces are included in this system through the Bond number which represents the competition between the hydrostatic and capillary forces. We then compare our results with results from literature of non-hardening fluids in each regime of Bond number. We see that the invasion patterns change strongly with hardening while the non-hardening behavior remains basically not affected.

We investigate the role of disorder on the fracturing process of heterogeneous materials by means of a two-dimensional fuse network model. Our results in the extreme disorder limit reveal that the backbone of the fracture at collapse, namely the subset of the largest fracture that effectively halts the global current, has a fractal dimension of $1.22 \pm 0.01$. This exponent value is compatible with the universality class of several other physical models, including optimal paths under strong disorder, disordered polymers, watersheds and optimal path cracks on uncorrelated substrates, hulls of explosive percolation clusters, and strands of invasion percolation fronts. Moreover, we find that the fractal dimension of the largest fracture under extreme disorder, $d_f=1.86 \pm 0.01$, is outside the statistical error bar of standard percolation. This discrepancy is due to the appearance of trapped regions or cavities of all sizes that remain intact till the entire collapse of the fuse network, but are always accessible in the case of standard percolation. Finally, we quantify the role of disorder on the structure of the largest cluster, as well as on the backbone of the fracture, in terms of a distinctive transition from weak to strong disorder characterized by a new crossover exponent., Comment: 5 pages, 4 figures

We investigate the effect of a temperature gradient on oil recovery in a two-dimensional pore-network model. The oil viscosity depends on temperature as, $\mu_o=exp(B/T)$, where $B$ is a physico-chemical parameter depending on the type of oil, and $T$ is the temperature. A temperature gradient is applied across the medium in the flow direction. Initially, the porous medium is saturated with oil and, then, another fluid is injected. We have considered two cases representing different injection strategies. In the first case, the invading fluid viscosity is constant (finite viscosity ratio) while in the second one, the invading fluid is inviscid (infinite viscosity ratio). Our results show that, for the case of finite viscosity ratio, recovery increases with $\Delta T$ independently on strength or sign of the gradient. For an infinite viscosity ratio, a positive temperature gradient is necessary to enhance recovery. Moreover, we show that, for $\Delta T>0$, the percentage of oil recovery generally decreases (increases) with $B$ for a finite (infinite) viscosity ratio. Finally, we also extend our results for infinite viscosity ratio to a three-dimensional porous media geometry.