Boundary-induced anisotropy of the avalanches in the sandpile automaton

We study numerically the avalanches in a two--dimensional critical height sandpile model with sand grains added at the center of the system. Smaller avalanches near the center of the system are isotropic. Larger avalanches are, however, affected by the boundary of the system, to a degree that increases with the avalanche size. Up to linear system size $L=1001$, we did not find an obvious indication for lattice--induced anisotropy.Comment: 7 pages, LaTeX, preprint HLRZ 38/9

Monte Carlo simulations of random copolymers at a selective interface

We investigate numerically using the bond--fluctuation model the adsorption of a random AB--copolymer at the interface between two solvents. From our results we infer several scaling relations: the radius of gyration of the copolymer in the direction perpendicular to the interface ($R_{gz}$) scales with $\chi$, the interfacial selectivity strength, as $R_{gz}=N^{\nu}f(\sqrt{N}\chi)$ where $\nu$ is the usual Flory exponent and $N$ is the copolymer's length; furthermore the monomer density at the interface scales as $\chi^{2\nu}$ for small $\chi$. We also determine numerically the monomer densities in the two solvents and discuss their dependence on the distance from the interface.Comment: Latex text file appended with figures.tar.g

Steady state properties of a driven granular medium

We study a two-dimensional granular system where external driving force is applied to each particle in the system in such a way that the system is driven into a steady state by balancing the energy input and the dissipation due to inelastic collision between particles. The velocities of the particles in the steady state satisfy the Maxwellian distribution. We measure the density-density correlation and the velocity-velocity correlation functions in the steady state and find that they are of power-law scaling forms. The locations of collision events are observed to be time-correlated and such a correlation is described by another power-law form. We also find that the dissipated energy obeys a power-law distribution. These results indicate that the system evolves into a critical state where there are neither characteristic spatial nor temporal scales in the correlation functions. A test particle exhibits an anomalous diffusion which is apparently similar to the Richardson law in a three-dimensional turbulent flow.Comment: REVTEX, submitted to Phys. Rev.

Density waves and 1/f1/f density fluctuations in granular flow

We simulate the granular flow in a narrow pipe with a lattice-gas automaton model. We find that the density in the system is characterized by two features. One is that spontaneous density waves propagate through the system with well-defined shapes and velocities. The other is that density waves are so distributed to make the power spectra of density fluctuations as $1/f^{\alpha}$ noise. Three important parameters make these features observable and they are energy dissipation, average density and the rougness of the pipe walls.Comment: Latex (with ps files appended

Velocity and density profiles of granular flow in channels using lattice gas automaton

We have performed two-dimensional lattice-gas-automaton simulations of granular flow between two parallel planes. We find that the velocity profiles have non-parabolic distributions while simultaneously the density profiles are non-uniform. Under non-slip boundary conditions, deviation of velocity profiles from the parabolic form of newtonian fluids is found to be characterized solely by ratio of maximal velocity at the center to the average velocity, though the ratio depends on the model parameters in a complex manner. We also find that the maximal velocity ($u_{max}$) at the center is a linear function of the driving force (g) as $u_{max} = \alpha g - \delta$ with non-zero $\delta$ in contrast with newtonian fluids. Regarding density profiles, we observe that densities near the boundaries are higher than those in the center. The width of higher densities (above the average density) relative to the channel width is a decreasing function of a variable which scales with the driving force (g), energy dissipation parameter ($\epsilon$) and the width of the system (L) as $g^{\mu} L^{\nu}/\epsilon$ with exponents $\mu = 1.4 \pm 0.1$ and $\nu = 0.5 \pm 0.1$. A phenomenological theory based on a scaling argument is presented to interpret these findings.Comment: Latex, 15 figures, to appear in PR