Ladder operator for the one-dimensional Hubbard model

The one-dimensional Hubbard model is integrable in the sense that it has an infinite family of conserved currents. We explicitly construct a ladder operator which can be used to iteratively generate all of the conserved current operators. This construction is different from that used for Lorentz invariant systems such as the Heisenberg model. The Hubbard model is not Lorentz invariant, due to the separation of spin and charge excitations. The ladder operator is obtained by a very general formalism which is applicable to any model that can be derived from a solution of the Yang-Baxter equation.

Integrability of the Russian doll BCS model

We show that integrability of the BCS model extends beyond Richardson's model (where all Cooper pair scatterings have equal coupling) to that of the Russian doll BCS model for which the couplings have a particular phase dependence that breaks time-reversal symmetry. This model is shown to be integrable using the quantum inverse scattering method, and the exact solution is obtained by means of the algebraic Bethe ansatz. The inverse problem of expressing local operators in terms of the global operators of the monodromy matrix is solved. This result is used to find a determinant formulation of a correlation function for fluctuations in the Cooper pair occupation numbers. These results are used to undertake exact numerical analysis for small systems at half-filling.

Turbulence in Hydraulic Jumps: Experimental Measurements

A hydraulic jump is the transition from a supercritical open channel flow to a subcritical regime. It is characterised by a highly turbulent flow with macro-scale vortices, some kinetic energy dissipation and a bubbly two-phase flow structure. New air-water flow measurements were performed in hydraulic jump flows for a range of inflow Froude numbers. The experiments were conducted in a large-size facility using two types of phase-detection intrusive probes: i.e., single-tip and double-tip conductivity probes. These were complemented by some measurements of free-surface fluctuations using ultrasonic displacement meters. The present study was focused on the turbulence characteristics of hydraulic jumps with partially-developed inflow conditions. The void fraction measurements showed the presence of an advective diffusion shear layer in which the void fractions profiles matched closely an analytical solution of the advective diffusion equation for air bubbles. The present results highlighted some influence of the inflow Froude number onto the air bubble entrainment process. At the largest Froude numbers, the advected air bubbles were more thoroughly dispersed vertically, and larger amount of air bubbles were detected in the turbulent shear layer. In the air-water mixing layer, the maximum void fraction and bubble count rate data showed some longitudinal decay function in the flow direction. Such trends were previously reported in the literature. The measurements of interfacial velocity and turbulence level distributions provided new information on the turbulent velocity field in the highly-aerated shear region. The present data suggested some longitudinal decay of the turbulence intensity. The velocity profiles tended to follow a wall jet flow pattern. The air–water turbulent time and length scales were deduced from some auto- and cross-correlation analyses based upon the method of CHANSON (2006,2007). The results provided the integral turbulent time and length scales of the eddy structures advecting the air bubbles in the developing shear layer. The experimental data showed that the auto-correlation time scale Txx was larger than the transverse cross-correlation time scale Txz. The integral turbulence length scale Lxz was a function of the inflow conditions, of the streamwise position (x-x1)/d1 and vertical elevation y/d1. Herein the dimensionless integral turbulent length scale Lxz/d1 was closely related to the inflow depth: i.e., Lxz/d1 = 0.2 to 0.8, with Lxz increasing towards the free-surface. The free-surface fluctuations measurements showed large turbulent fluctuations that reflected the dynamic, unsteady structure of the hydraulic jumps. A linear relationship was found between the normalized maximum free-surface fluctuation and the inflow Froude number.

Exact solution of the A(n-1)((1)) trigonometric vertex model with non-diagonal open boundaries

The A(n-1)((1)) trigonometric vertex model with generic non-diagonal boundaries is studied. The double-row transfer matrix of the model is diagonalized by algebraic Bethe ansatz method in terms of the intertwiner and the corresponding face-vertex relation. The eigenvalues and the corresponding Bethe ansatz equations are obtained.

Superconducting Correlations in Metallic Nanoparticles: Exact Solution of the BCS Model by the Algebraic Bethe Ansatz

Superconducting pairing of electrons in nanoscale metallic particles with discrete energy levels and a fixed number of electrons is described by the reduced Bardeen, Cooper, and Schrieffer model Hamiltonian. We show that this model is integrable by the algebraic Bethe ansatz. The eigenstates, spectrum, conserved operators, integrals of motion, and norms of wave functions are obtained. Furthermore, the quantum inverse problem is solved, meaning that form factors and correlation functions can be explicitly evaluated. Closed form expressions are given for the form factors and correlation functions that describe superconducting pairing.

A master equation model for bimolecular reaction via multi-well isomerization intermediates

A reversible linear master equation model is presented for pressure- and temperature-dependent bimolecular reactions proceeding via multiple long-lived intermediates. This kinetic treatment, which applies when the reactions are measured under pseudo-first-order conditions, facilitates accurate and efficient simulation of the time dependence of the populations of reactants, intermediate species and products. Detailed exploratory calculations have been carried out to demonstrate the capabilities of the approach, with applications to the bimolecular association reaction C3H6 + H reversible arrow C3H7 and the bimolecular chemical activation reaction C2H2 +(CH2)-C-1--> C3H3+H. The efficiency of the method can be dramatically enhanced through use of a diffusion approximation to the master equation, and a methodology for exploiting the sparse structure of the resulting rate matrix is established.

On the validity of the shrinking core model in noncatalytic gas solid reaction

By using a matched asymptotic expansion technique, the shrinking core model (SCM) used in non-catalytic gas solid reactions with general kinetic expression is rigorously justified in this paper as a special case of the homogeneous model when the reaction rate is much faster than that of diffusion. The time-pendent velocity of the moving reacted-unreacted interface is found to be proportional to the gas flux at that interface for all geometries of solid particles, and the thickness order of the reaction zone and also the degree of chemical reaction at the interface is discussed in this paper.

High-Frequency Turbulence and Sediment Flux Measurements in an Upper Estuarine Zone

In natural estuaries, the predictions of scalar dispersion are rarely predicted accurately because of a lack of fundamental understanding of the turbulence structure in estuaries. Herein detailed turbulence field measurements were conducted continuously at high frequency for 50 hours in the upper zone of a small subtropical estuary with semi-diurnal tides. Acoustic Doppler velocimetry was deemed the most appropriate measurement technique for such shallow water depths (less than 0.4 m at low tides), and a thorough post-processing technique was applied. In addition, some experiments were conducted in laboratory under controlled conditions using water and soil samples collected in the estuary to test the relationship between acoustic backscatter strength and suspended sediment load. A striking feature of the field data set was the large fluctuations in all turbulence characteristics during the tidal cycle, including the suspended sediment flux. This feature was rarely documented.

Turbulence Characterisation of High-Velocity Free-Surface Flows

In high-velocity free-surface flows, air is continuously being trapped and released through the free-surface. Such high-velocity highly-aerated flows cannot be studied numerically because of the large number of relevant equations and parameters. Herein an advanced signal processing of traditional single- and dual-tip conductivity probes provides some new information on the air-water turbulent time and length scales. The technique is applied to turbulent open channel flows in a large-size facility. The auto- and cross-correlation analyses yield some characterisation of the large eddies advecting the bubbles. The transverse integral turbulent length and time scales are related to the step height: i.e., Lxy/h ~ 0.02 to 0.2, and T.sqrt(g/h) ~ 0.004 to 0.04. The results are irrespective of the Reynolds numbers. The present findings emphasise that turbulent dissipation by large-scale vortices is a significant process in the intermediate zone between the spray and bubbly flow regions (0.3 < C < 0.7). Some self-similar relationships were observed systematically at both macroscopic and microscopic levels. The results are significant because they provide a picture general enough to be used to characterise the air-water flow field in prototype spillways.

Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.

Combining GIS and Expert Opinion to Model Landscapes for a Smallholder Forest Extension Program in Leyte, the Philippines

A Geographic Information System (GIS) was used to model datasets of Leyte Island, the Philippines, to identify land which was suitable for a forest extension program on the island. The datasets were modelled to provide maps of the distance of land from cities and towns, land which was a suitable elevation and slope for smallholder forestry and land of various soil types. An expert group was used to assign numeric site suitabilities to the soil types and maps of site suitability were used to assist the selection of municipalities for the provision of extension assistance to smallholders. Modelling of the datasets was facilitated by recent developments of the ArcGIS® suite of computer programs and derivation of elevation and slope was assisted by the availability of digital elevation models (DEM) produced by the Shuttle Radar Topography (SRTM) mission. The usefulness of GIS software as a decision support tool for small-scale forestry extension programs is discussed.

Parliament House Competition Model, Canberra

View of model for competition entry.