International Journal of Modern Physics C

In this work, an explicit Taylor series expansion- and least square-based lattice Boltzmann method (LBM) is used to simulate two-dimensional unsteady incompressible viscous flows. The new method is based on the standard LBM with introduction of the Taylor series expansion and the least squares approach. The final equation is an explicit form and essentially has no limitation on mesh structure and lattice model. Since the Taylor series expansion is only applied in the spatial direction, the time accuracy of the new method is kept the same as the standard LBM, which seems to benefit for unsteady flow simulation. To validate the new method, two test problems, that is, the vortex shedding behind a circular cylinder at low Reynolds numbers and the oscillating flow in a lid driven cavity, were considered in this work. Numerical results obtained by the new method agree very well with available data in the literature.

We study the effects of on-ramp and off-ramp to the traffic flow in a cellular automaton. The stochastic noise and the extended hopping are included. The phase diagram is obtained. Three distinct phases are observed. The density profiles are analyzed. In the maximum current phase, standing traffic jams emerge in the upstream of the on-ramp. A region of free flow is observed in the downstream of the off-ramp. In between the ramps, a scaled density profile is observed in one branch and a phase separation is observed in the other branch.

Front propagation in a random environment is studied close to the depinning threshold. At zero temperature we show that the depinning force distribution exhibits a universal behavior. This property is used to estimate the velocity of the front at very low temperature. We obtain a Arrhenius behavior with a prefactor depending on the temperature as a power law. These results are supported by numerical simulations.

We have simulated the gold deposition on arsenic and gallium terminated GaAs(001) surfaces at low and room temperatures. It has been found that gallium terminated surface is relatively less stable in comparison to the arsenic terminated surface. On the other hand, a single gold adatom on the surface has different characteristics than full coverage gold atoms on the surface; a single gold atom diffuses into the surface at room temperature. Simulations have been performed by considering classical molecular-dynamics technique using an empirical many-body potential energy function comprising two- and three-body atomic interactions.

We study the behavior of one-dimensional chaotic signals and filtering using the discrete Haar wavelet and Daubechies wavelets.

In order to study fluctuations in percolating systems, lattices for sizes up to L = 100 000 have been simulated several thousand times using the Hoshen–Kopelman algorithm. Distributions of cluster numbers are Gaussians for small clusters and half-sided quasi-Gaussians for large clusters. The variance of cluster numbers is proportional to the mean, with power-law deviations for small clusters. Higher moments like skewness and kurtosis were also studied.

Complementarity is one of the main features underlying the interactions in biological and biochemical systems. Inspired by those systems we propose a model for the dynamical evolution of a system composed of agents that interact due to their complementary attributes rather than their similarities. Each agent is represented by a bit-string and has an activity associated to it; the coupling among complementary peers depends on their activity. The connectivity of the system changes in time respecting the constraint of complementarity. We observe the formation of a network of active agents whose stability depends on the rate at which activity of agents diffuses in the system. The model exhibits a nonequilibrium phase transition between the ordered phase, where a stable network is generated, and a disordered phase characterized by the absence of correlation among the agents. The ordered phase exhibits multi-modal distributions of connectivity and activity, indicating a hierarchy of interaction among different populations characterized by different degrees of activity. This model may be used to study the hierarchy observed in social organizations.

A melt of nonentangled polymer chains confined between two smooth and purely repulsive walls is studied for various film thicknesses D and temperatures. The dynamics of the supercooled films is qualitatively identical to that of the bulk, but the walls lead to faster relaxation. To quantify this observation we analyze the data by the mode-coupling theory (MCT) of the glass transition. We find that the critical temperature of MCT, Tc(D), decreases with D and that T - Tc(D) is a relevant temperature scale. The static structure factor and dynamic correlation functions at intermediate times coincide with bulk behavior when compared to the same T - Tc(D).

The article presents a series of numerical simulations of exact solutions of the Einstein equations performed using the Cactus code, a complete three-dimensional machinery for numerical relativity. We describe an application ("thorn") for the Cactus code that can be used for evolving a variety of exact solutions, with and without matter, including solutions used in modern cosmology for modeling the early stages of the universe. Our main purpose has been to test the Cactus code on these well-known examples, focusing mainly on the stability and convergence of the code.

The photon anticorrelation experiment is simulated with the unquantized Maxwell field by implementing realistic instrument bandwidth and signal to noise ratio. The results strongly violate an inequality believed to be satisfied by all semiclassical theories, which opens up a classical interpretation that rivals quantum optics and casts doubt on potential quantum computer applications that involve optical fields.

The Duke–Rubinstein model of gel electrophoresis is applied to calculate the velocity of DNA molecules. We have found that the velocity distribution becomes flat at high electric fields. Simultaneously, the percentage of immobile molecules increases. Effectively, the mean velocity starts to decrease at high fields. The field value, where the mean velocity is maximal, decreases with the molecule length. The results are compared with those from similar calculations obtained by Heukelum and Beljaars within the cage model.

The problem of knowing and characterizing the transitive behavior of a given cellular automaton is a very interesting topic. This paper provides a matrix representation of the global dynamics in reversible one-dimensional cellular automata with a Welch index 1, i.e., those where the ancestors differ just at one end. We prove that the transitive closure of this matrix shows diverse types of transitive behaviors in these systems. Part of the theorems in this paper are reductions of well-known results in symbolic dynamics. This matrix and its transitive closure were computationally implemented, and some examples are presented.