Search

Viscoelastic fluids exhibit elastic instabilities in simple shear flow and flow through curved streamlines. Surprisingly, we found in a porous medium such fluids show strikingly different hydrodynamic instabilities depicted by very large sideways excursions and presence of fast and slow moving lanes which have not been reported before. Particle image velocimetry (PIV) measurements through a pillared microchannel, provide experimental evidence of such instabilities at very low Reynolds number (< 0.01). We observe a transition from a symmetric laminar to an asymmetric flow, which finally transforms to a nonlinear aperiodic flow with strong lateral movements. The instability is characterized by a rapid increase in spatial and temporal fluctuations of velocity components and pressure at a critical Deborah number (De). Our experiments reveal the presence of a fascinating interplay between pore space and fluid rheology., Comment: 11 pages, 8 figures

We describe the implementation of output code optimization in the open source computer algebra system FORM. This implementation is based on recently discovered techniques of Monte Carlo tree search to find efficient multivariate Horner schemes, in combination with other optimization algorithms, such as common subexpression elimination. For systems for which no specific knowledge is provided it performs significantly better than other methods we could compare with. Because the method has a number of free parameters, we also show some methods by which to tune them to different types of problems., Comment: 31 pages

Electronic transport through chaotic quantum dots exhibits universal behaviour which can be understood through the semiclassical approximation. Within the approximation, transport moments reduce to codifying classical correlations between scattering trajectories. These can be represented as ribbon graphs and we develop an algorithmic combinatorial method to generate all such graphs with a given genus. This provides an expansion of the linear transport moments for systems both with and without time reversal symmetry. The computational implementation is then able to progress several orders higher than previous semiclassical formulae as well as those derived from an asymptotic expansion of random matrix results. The patterns observed also suggest a general form for the higher orders., Comment: 37 pages; 19 figures

To study electronic transport through chaotic quantum dots, there are two main theoretical approachs. One involves substituting the quantum system with a random scattering matrix and performing appropriate ensemble averaging. The other treats the transport in the semiclassical approximation and studies correlations among sets of classical trajectories. There are established evaluation procedures within the semiclassical evaluation that, for several linear and non-linear transport moments to which they were applied, have always resulted in the agreement with random matrix predictions. We prove that this agreement is universal: any semiclassical evaluation within the accepted procedures is equivalent to the evaluation within random matrix theory. The equivalence is shown by developing a combinatorial interpretation of the trajectory sets as ribbon graphs (maps) with certain properties and exhibiting systematic cancellations among their contributions. Remaining trajectory sets can be identified with primitive (palindromic) factorisations whose number gives the coefficients in the corresponding expansion of the moments of random matrices. The equivalence is proved for systems with and without time reversal symmetry., Comment: 30 pages, 13 figures; added some examples

Optimizing the cost of evaluating a polynomial is a classic problem in computer science. For polynomials in one variable, Horner's method provides a scheme for producing a computationally efficient form. For multivariate polynomials it is possible to generalize Horner's method, but this leaves freedom in the order of the variables. Traditionally, greedy schemes like most-occurring variable first are used. This simple textbook algorithm has given remarkably efficient results. Finding better algorithms has proved difficult. In trying to improve upon the greedy scheme we have implemented Monte Carlo tree search, a recent search method from the field of artificial intelligence. This results in better Horner schemes and reduces the cost of evaluating polynomials, sometimes by factors up to two., Comment: 5 pages

We present version 4.0 of the symbolic manipulation system FORM. The most important new features are manipulation of rational polynomials and the factorization of expressions. Many other new functions and commands are also added; some of them are very general, while others are designed for building specific high level packages, such as one for Groebner bases. New is also the checkpoint facility, that allows for periodic backups during long calculations. Lastly, FORM 4.0 has become available as open source under the GNU General Public License version 3., Comment: 26 pages. Uses axodraw

We confirm a conjecture about the construction of basis elements for the multiple zeta values (MZVs) at weight 27 and weight 28. Both show as expected one element that is twofold extended. This is done with some lengthy computer algebra calculations using TFORM to determine explicit bases for the MZVs at these weights., Comment: 6 pages

Most theories of homogeneous nucleation are based on a Fokker-Planck-like description of the behavior of the mass of clusters. Here we will show that these approaches are incomplete for a large class of nucleating systems, as they assume the effective dynamics of the clusters to be Markovian, i.e., memoryless. We characterize these non-Markovian dynamics and show how this influences the dynamics of clusters during nucleation. Our results are validated by simulations of a three-dimensional Ising model with locally conserved magnetization., Comment: 4 pages, 4 figures

In many systems, the time scales of the microscopic dynamics and macroscopic dynamics of interest are separated by many orders of magnitude. Examples abound, for instance nucleation, protein folding, and chemical reactions. For these systems, direct simulation of phase space trajectories does not efficiently determine most physical quantities of interest. The last decade has seen the advent of methods circumventing brute force simulation. For most dynamical quantities, these methods all share the drawback of systematical errors. We present a novel method for generating ensembles of phase space trajectories. By sampling small pieces of these trajectories in different phase space domains and piecing them together in a smart way using equilibrium properties, we obtain physical quantities such as transition times. This method does not have any systematic error and is very efficient; the computational effort to calculate the first passage time across a free energy barrier does not increase with the height of the barrier. The strength of the method is shown in the Ising model. Accurate measurements of nucleation times span almost ten orders of magnitude and reveal corrections to classical nucleation theory., Comment: 8 pages, 4 figures

The time evolution of the strength of the Earth's virtual axial dipole moment (VADM) is analyzed by relating it to the Fokker-Planck equation, which describes a random walk with VADM-dependent drift and diffusion coefficients. We demonstrate first that our method is able to retrieve the correct shape of the drift and diffusion coefficients from a time series generated by a test model. Analysis of the Sint-2000 data shows that the geomagnetic dipole mode has a linear growth time of 13 to 33 kyr, and that the nonlinear quenching of the growth rate follows a quadratic function of the type [1-(x/x0)^2]. On theoretical grounds, the diffusive motion of the VADM is expected to be driven by multiplicative noise, and the corresponding diffusion coefficient to scale quadratically with dipole strength. However, analysis of the Sint-2000 VADM data reveals a diffusion which depends only very weakly on the dipole strength. This may indicate that the magnetic field quenches the amplitude of the turbulent velocity in the Earth's outer core., Comment: 11 pages, 6 figures