Solution of the dual reflection equation for An−1(1)A^{(1)}_{n-1} SOS model

We obtain a diagonal solution of the dual reflection equation for elliptic $A^{(1)}_{n-1}$ SOS model. The isomorphism between the solutions of the reflection equation and its dual is studied.Comment: Latex file 12 pages, added reference

Random walks on finite lattice tubes

Exact results are obtained for random walks on finite lattice tubes with a single source and absorbing lattice sites at the ends. Explicit formulae are derived for the absorption probabilities at the ends and for the expectations that a random walk will visit a particular lattice site before being absorbed. Results are obtained for lattice tubes of arbitrary size and each of the regular lattice types; square, triangular and honeycomb. The results include an adjustable parameter to model the effects of strain, such as surface curvature, on the surface diffusion. Results for the triangular lattice tubes and the honeycomb lattice tubes model diffusion of adatoms on single walled zig-zag carbon nano-tubes with open ends.Comment: 22 pages, 4 figure

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.

Transport in a highly asymmetric binary fluid mixture

We present molecular dynamics calculations of the thermal conductivity and viscosities of a model colloidal suspension with colloidal particles roughly one order of magnitude larger than the suspending liquid molecules. The results are compared with estimates based on the Enskog transport theory and effective medium theories (EMT) for thermal and viscous transport. We find, in particular, that EMT remains well applicable for predicting both the shear viscosity and thermal conductivity of such suspensions when the colloidal particles have a ``typical'' mass, i.e. much larger than the liquid molecules. Very light colloidal particles on the other hand yield higher thermal conductivities, in disagreement with EMT. We also discuss the consequences of these results to some proposed mechanisms for thermal conduction in nanocolloidal suspensions.Comment: 13 pages, 6 figures, to appear in Physical Review E (2007

Locating the source of projectile fluid droplets

The ill-posed projectile problem of finding the source height from spattered droplets of viscous fluid is a longstanding obstacle to accident reconstruction and crime scene analysis. It is widely known how to infer the impact angle of droplets on a surface from the elongation of their impact profiles. However, the lack of velocity information makes finding the height of the origin from the impact position and angle of individual drops not possible. From aggregate statistics of the spatter and basic equations of projectile motion, we introduce a reciprocal correlation plot that is effective when the polar launch angle is concentrated in a narrow range. The vertical coordinate depends on the orientation of the spattered surface, and equals the tangent of the impact angle for a level surface. When the horizontal plot coordinate is twice the reciprocal of the impact distance, we can infer the source height as the slope of the data points in the reciprocal correlation plot. If the distribution of launch angles is not narrow, failure of the method is evident in the lack of linear correlation. We perform a number of experimental trials, as well as numerical calculations and show that the height estimate is insensitive to aerodynamic drag. Besides its possible relevance for crime investigation, reciprocal-plot analysis of spatter may find application to volcanism and other topics and is most immediately applicable for undergraduate science and engineering students in the context of crime-scene analysis.Comment: To appear in the American Journal of Physics (ms 23338). Improved readability and organization in this versio

Velocity, energy and helicity of vortex knots and unknots

In this paper we determine the velocity, the energy and estimate writhe and twist helicity contributions of vortex filaments in the shape of torus knots and unknots (toroidal and poloidal coils) in a perfect fluid. Calculations are performed by numerical integration of the Biot-Savart law. Vortex complexity is parametrized by the winding number $w$, given by the ratio of the number of meridian wraps to that of the longitudinal wraps. We find that for $w<1$ vortex knots and toroidal coils move faster and carry more energy than a reference vortex ring of same size and circulation, whereas for $w>1$ knots and poloidal coils have approximately same speed and energy of the reference vortex ring. Helicity is dominated by the writhe contribution. Finally, we confirm the stabilizing effect of the Biot-Savart law for all knots and unknots tested, that are found to be structurally stable over a distance of several diameters. Our results also apply to quantized vortices in superfluid $^4$He.Comment: 17 pages, 8 figures, 2 table

Turbulent Pair Diffusion

Kinematic Simulations of turbulent pair diffusion in planar turbulence with a -5/3 energy spectrum reproduce the results of the laboratory measurements of Jullien Phys. Rev. Lett. 82, 2872 (1999), in particular the stretched exponential form of the PDF of pair separations and their correlation functions. The root mean square separation is found to be strongly dependent on initial conditions for very long stretches of times. This dependence is consistent with the topological picture of turbulent pair diffusion where pairs initially close enough travel together for long stretches of time and separate violently when they meet straining regions around hyperbolic points. A new argument based on the divergence of accelerations is given to support this picture

Exact solution and surface critical behaviour of open cyclic SOS lattice models

We consider the $L$-state cyclic solid-on-solid lattice models under a class of open boundary conditions. The integrable boundary face weights are obtained by solving the reflection equations. Functional relations for the fused transfer matrices are presented for both periodic and open boundary conditions. The eigen-spectra of the unfused transfer matrix is obtained from the functional relations using the analytic Bethe ansatz. For a special case of crossing parameter $\lambda=\pi/L$, the finite-size corrections to the eigen-spectra of the critical models are obtained, from which the corresponding conformal dimensions follow. The calculation of the surface free energy away from criticality yields two surface specific heat exponents, $\alpha_s=2-L/2\ell$ and $\alpha_1=1-L/\ell$, where $\ell=1,2,\cdots,L-1$ coprime to $L$. These results are in agreement with the scaling relations $\alpha_s=\alpha_b+\nu$ and $\alpha_1=\alpha_b-1$.Comment: 13 pages, LaTeX, to appear in J. Phys.

The packing of two species of polygons on the square lattice

We decorate the square lattice with two species of polygons under the constraint that every lattice edge is covered by only one polygon and every vertex is visited by both types of polygons. We end up with a 24 vertex model which is known in the literature as the fully packed double loop model. In the particular case in which the fugacities of the polygons are the same, the model admits an exact solution. The solution is obtained using coordinate Bethe ansatz and provides a closed expression for the free energy. In particular we find the free energy of the four colorings model and the double Hamiltonian walk and recover the known entropy of the Ice model. When both fugacities are set equal to two the model undergoes an infinite order phase transition.Comment: 21 pages, 4 figure

Inertial range scaling in numerical turbulence with hyperviscosity

Numerical turbulence with hyperviscosity is studied and compared with direct simulations using ordinary viscosity and data from wind tunnel experiments. It is shown that the inertial range scaling is similar in all three cases. Furthermore, the bottleneck effect is approximately equally broad (about one order of magnitude) in these cases and only its height is increased in the hyperviscous case--presumably as a consequence of the steeper decent of the spectrum in the hyperviscous subrange. The mean normalized dissipation rate is found to be in agreement with both wind tunnel experiments and direct simulations. The structure function exponents agree with the She-Leveque model. Decaying turbulence with hyperviscosity still gives the usual t^{-1.25} decay law for the kinetic energy, and also the bottleneck effect is still present and about equally strong.Comment: Final version (7 pages