Towards realistic simulations of lava dome growth using the level set method

The level set method has been implemented in a computational volcanology context. New techniques are presented to solve the advection equation and the reinitialisation equation. These techniques are based upon an algorithm developed in the finite difference context, but are modified to take advantage of the robustness of the finite element method. The resulting algorithm is tested on a well documented Rayleigh–Taylor instability benchmark [19], and on an axisymmetric problem where the analytical solution is known. Finally, the algorithm is applied to a basic study of lava dome growth.

Determination of coalescence kernels for high-shear granulation using DEM simulations

Discrete element method (DEM) modeling is used in parallel with a model for coalescence of deformable surface wet granules. This produces a method capable of predicting both collision rates and coalescence efficiencies for use in derivation of an overall coalescence kernel. These coalescence kernels can then be used in computationally efficient meso-scale models such as population balance equation (PBE) models. A soft-sphere DEM model using periodic boundary conditions and a unique boxing scheme was utilized to simulate particle flow inside a high-shear mixer. Analysis of the simulation results provided collision frequency, aggregation frequency, kinetic energy, coalescence efficiency and compaction rates for the granulation process. This information can be used to bridge the gap in multi-scale modeling of granulation processes between the micro-scale DEM/coalescence modeling approach and a meso-scale PBE modeling approach.

Gauge P-representations for quantum-dynamical problems: Removal of boundary terms

P-representation techniques, which have been very successful in quantum optics and in other fields, are also useful for general bosonic quantum-dynamical many-body calculations such as Bose-Einstein condensation. We introduce a representation called the gauge P representation, which greatly widens the range of tractable problems. Our treatment results in an infinite set of possible time evolution equations, depending on arbitrary gauge functions that can be optimized for a given quantum system. In some cases, previous methods can give erroneous results, due to the usual assumption of vanishing boundary conditions being invalid for those particular systems. Solutions are given to this boundary-term problem for all the cases where it is known to occur: two-photon absorption and the single-mode laser. We also provide some brief guidelines on how to apply the stochastic gauge method to other systems in general, quantify the freedom of choice in the resulting equations, and make a comparison to related recent developments.

On a semilinear Schrodinger equation with critical Sobolev exponent

We consider the semilinear Schrodinger equation -Deltau+V(x)u= K(x) \u \ (2*-2 u) + g(x; u), u is an element of W-1,W-2 (R-N), where N greater than or equal to4, V, K, g are periodic in x(j) for 1 less than or equal toj less than or equal toN, K>0, g is of subcritical growth and 0 is in a gap of the spectrum of -Delta +V. We show that under suitable hypotheses this equation has a solution u not equal 0. In particular, such a solution exists if K equivalent to 1 and g equivalent to 0.

Existence of Multiple Solutions for Second Order Boundary Value Problems

We give conditions on f involving pairs of lower and upper solutions which lead to the existence of at least three solutions of the two point boundary value problem y" + f(x, y, y') = 0, x epsilon [0, 1], y(0) = 0 = y(1). In the special case f(x, y, y') = f(y) greater than or equal to 0 we give growth conditions on f and apply our general result to show the existence of three positive solutions. We give an example showing this latter result is sharp. Our results extend those of Avery and of Lakshmikantham et al.

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.

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.

Integrable eight-state supersymmetric U model with boundary terms and its Bethe ansatz solution

A class of integrable boundary terms for the eight-state supersymmetric U model are presented by solving the graded reflection equations. The boundary model is solved by using the coordinate Bethe ansatz method and the Bethe ansatz equations are obtained. (C) 1998 Elsevier Science B.V.

An open-boundary integrable model of three coupled XY spin chains

The integrable open-boundary conditions for the model of three coupled one-dimensional XY spin chains are considered in the framework of the quantum inverse scattering method. The diagonal boundary K-matrices are found and a class of integrable boundary terms is determined. The boundary model Hamiltonian is solved by using the coordinate space Bethe ansatz technique and Bethe ansatz equations are derived. (C) 1998 Elsevier Science B.V.

Quantum integrability and exact solution of the supersymmetric U model with boundary terms

Quantum integrability is established for the one-dimensional supersymmetric U model with boundary terms by means of the quantum inverse-scattering method. The boundary supersymmetric U chain is solved by using the coordinate-space Bethe-ansatz technique and Bethe-ansatz equations are derived. This provides us with a basis for computing the finite-size corrections to the low-lying energies in the system. [S0163-1829(98)00425-1].

Nine classes of integrable boundary conditions for the eight-state supersymmetric fermion model

Nine classes of integrable boundary conditions for the eight-state supersymmetric model of strongly correlated fermions are presented. The boundary systems are solved by using the coordinate Bethe ansatz method and the Bethe ansatz equations for all nine cases are given.

Algebraic Bethe ansatz for integrable Kondo impurities in the one-dimensional supersymmetric t-J model

An integrable Kondo problem in the one-dimensional supersymmetric t-J model is studied by means of the boundary supersymmetric quantum inverse scattering method. The boundary K matrices depending on the local moments of the impurities are presented as a nontrivial realization of the graded reflection equation algebras in a two-dimensional impurity Hilbert space. Further, the model is solved by using the algebraic Bethe ansatz method and the Bethe ansatz equations are obtained. (C) 1999 Elsevier Science B.V.

Integrable Kondo impurities in the one-dimensional supersymmetric U model of strongly correlated electrons

Integrable Kondo impurities in the one-dimensional supersymmetric U model of strongly correlated electrons are studied by means of the boundary graded quantum inverse scattering method. The boundary K-matrices depending on the local magnetic moments of the impurities are presented as non-trivial realizations of the reflection equation algebras in an impurity Hilbert space. Furthermore, the model Hamiltonian is diagonalized and the Bethe ansatz equations are derived. It is interesting to note that our model exhibits a free parameter in the bulk Hamiltonian but no free parameter exists on the boundaries. This is in sharp contrast to the impurity models arising from the supersymmetric t-J and extended Hubbard models where there is no free parameter in the bulk but there is a free parameter on each boundary.

Integrable Kondo impurities in the one-dimensional supersymmetric extended Hubbard model

An integrable Kondo problem in the one-dimensional supersymmetric extended Hubbard model is studied by means of the boundary graded quantum inverse scattering method. The boundary K-matrices depending on the local moments of the impurities are presented as a non-trivial realization of the graded reflection equation algebras in a two-dimensional impurity Hilbert space. Further, the model is solved by using the algebraic Bethe ansatz method and the Bethe ansatz equations are obtained.