Magnetic Susceptibility of an integrable anisotropic spin ladder system

We investigate the thermodynamics of a spin ladder model which possesses a free parameter besides the rung and leg couplings. The model is exactly solved by the Bethe Ansatz and exhibits a phase transition between a gapped and a gapless spin excitation spectrum. The magnetic susceptibility is obtained numerically and its dependence on the anisotropy parameter is determined. A connection with the compounds KCuCl3, Cu2(C5H12N2)2Cl4 and (C5H12N)2CuBr4 in the strong coupling regime is made and our results for the magnetic susceptibility fit the experimental data remarkably well.Comment: 12 pages, 12 figures included, submitted to Phys. Rev.

Exact solution for random walks on the triangular lattice with absorbing boundaries

The problem of a random walk on a finite triangular lattice with a single interior source point and zig-zag absorbing boundaries is solved exactly. This problem has been previously considered intractable.Comment: 10 pages, Latex, IOP macro

Motion and homogenization of vortices in anisotropic Type II superconductors

The motion of vortices in an anisotropic superconductor is considered. For a system of well-separated vortices, each vortex is found to obey a law of motion analogous to the local induction approximation, in which velocity of the vortex depends upon the local curvature and orientation. A system of closely packed vortices is then considered, and a mean field model is formulated in which the individual vortex lines are replaced by a vortex density

Optimal model parameters for multi-objective large-eddy simulations

A methodology is proposed for the assessment of error dynamics in large-eddy simulations. It is demonstrated that the optimization of model parameters with respect to one flow property can be obtained at the expense of the accuracy with which other flow properties are predicted. Therefore, an approach is introduced which allows to assess the total errors based on various flow properties simultaneously. We show that parameter settings exist, for which all monitored errors are "near optimal," and refer to such regions as "multi-objective optimal parameter regions." We focus on multi-objective errors that are obtained from weighted spectra, emphasizing both large- as well small-scale errors. These multi-objective optimal parameter regions depend strongly on the simulation Reynolds number and the resolution. At too coarse resolutions, no multi-objective optimal regions might exist as not all error-components might simultaneously be sufficiently small. The identification of multi-objective optimal parameter regions can be adopted to effectively compare different subgrid models. A comparison between large-eddy simulations using the Lilly-Smagorinsky model, the dynamic Smagorinsky model and a new Re-consistent eddy-viscosity model is made, which illustrates this. Based on the new methodology for error assessment the latter model is found to be the most accurate and robust among the selected subgrid models, in combination with the finite volume discretization used in the present study

The Generation of Magnetic Fields Through Driven Turbulence

We have tested the ability of driven turbulence to generate magnetic field structure from a weak uniform field using three dimensional numerical simulations of incompressible turbulence. We used a pseudo-spectral code with a numerical resolution of up to $144^3$ collocation points. We find that the magnetic fields are amplified through field line stretching at a rate proportional to the difference between the velocity and the magnetic field strength times a constant. Equipartition between the kinetic and magnetic energy densities occurs at a scale somewhat smaller than the kinetic energy peak. Above the equipartition scale the velocity structure is, as expected, nearly isotropic. The magnetic field structure at these scales is uncertain, but the field correlation function is very weak. At the equipartition scale the magnetic fields show only a moderate degree of anisotropy, so that the typical radius of curvature of field lines is comparable to the typical perpendicular scale for field reversal. In other words, there are few field reversals within eddies at the equipartition scale, and no fine-grained series of reversals at smaller scales. At scales below the equipartition scale, both velocity and magnetic structures are anisotropic; the eddies are stretched along the local magnetic field lines, and the magnetic energy dominates the kinetic energy on the same scale by a factor which increases at higher wavenumbers. We do not show a scale-free inertial range, but the power spectra are a function of resolution and/or the imposed viscosity and resistivity. Our results are consistent with the emergence of a scale-free inertial range at higher Reynolds numbers.Comment: 14 pages (8 NEW figures), ApJ, in press (July 20, 2000?

A refined Razumov-Stroganov conjecture II

We extend a previous conjecture [cond-mat/0407477] relating the Perron-Frobenius eigenvector of the monodromy matrix of the O(1) loop model to refined numbers of alternating sign matrices. By considering the O(1) loop model on a semi-infinite cylinder with dislocations, we obtain the generating function for alternating sign matrices with prescribed positions of 1's on their top and bottom rows. This seems to indicate a deep correspondence between observables in both models.Comment: 21 pages, 10 figures (3 in text), uses lanlmac, hyperbasics and epsf macro

Pearling instability of nanoscale fluid flow confined to a chemical channel

We investigate the flow of a nano-scale incompressible ridge of low-volatility liquid along a "chemical channel": a long, straight, and completely wetting stripe embedded in a planar substrate, and sandwiched between two extended less wetting solid regions. Molecular dynamics simulations, a simple long-wavelength approximation, and a full stability analysis based on the Stokes equations are used, and give qualitatively consistent results. While thin liquid ridges are stable both statically and during flow, a (linear) pearling instability develops if the thickness of the ridge exceeds half of the width of the channel. In the flowing case periodic bulges propagate along the channel and subsequently merge due to nonlinear effects. However, the ridge does not break up even when the flow is unstable, and the qualitative behavior is unchanged even when the fluid can spill over onto a partially wetting exterior solid region.Comment: 17 pages, 12 figures, submitted to Physics of Fluids, fixed equation numbering after Eq. (17