Bounds on integrals of the Wigner function: The hyperbolic case

Wigner functions play a central role in the phase space formulation of quantum mechanics. Although closely related to classical Liouville densities, Wigner functions are not positive definite and may take negative values on subregions of phase space. We investigate the accumulation of these negative values by studying bounds on the integral of an arbitrary Wigner function over noncompact subregions of the phase plane with hyperbolic boundaries. We show using symmetry techniques that this problem reduces to computing the bounds on the spectrum associated with an exactly solvable eigenvalue problem and that the bounds differ from those on classical Liouville distributions. In particular, we show that the total "quasiprobability" on such a region can be greater than 1 or less than zero. (C) 2005 American Institute of Physics.

On the chaotic instability of a nonsliding liquid-filled top with a small spheroidal base via Melnikov-Holmes-Marsden integrals

Chaotic orientations of a top containing a fluid filled cavity are investigated analytically and numerically under small perturbations. The top spins and rolls in nonsliding contact with a rough horizontal plane and the fluid in the ellipsoidal shaped cavity is considered to be ideal and describable by finite degrees of freedom. A Hamiltonian structure is established to facilitate the application of Melnikov-Holmes-Marsden (MHM) integrals. In particular, chaotic motion of the liquid-filled top is identified to be arisen from the transversal intersections between the stable and unstable manifolds of an approximated, disturbed flow of the liquid-filled top via the MHM integrals. The developed analytical criteria are crosschecked with numerical simulations via the 4th Runge-Kutta algorithms with adaptive time steps.

Numerical methods for strong solutions of stochastic differential equations: an overview

This paper gives a review of recent progress in the design of numerical methods for computing the trajectories (sample paths) of solutions to stochastic differential equations. We give a brief survey of the area focusing on a number of application areas where approximations to strong solutions are important, with a particular focus on computational biology applications, and give the necessary analytical tools for understanding some of the important concepts associated with stochastic processes. We present the stochastic Taylor series expansion as the fundamental mechanism for constructing effective numerical methods, give general results that relate local and global order of convergence and mention the Magnus expansion as a mechanism for designing methods that preserve the underlying structure of the problem. We also present various classes of explicit and implicit methods for strong solutions, based on the underlying structure of the problem. Finally, we discuss implementation issues relating to maintaining the Brownian path, efficient simulation of stochastic integrals and variable-step-size implementations based on various types of control.

Application of X-ray microtomography to numerical simulations of agglomerate breakage by distinct element method

Computer-aided tomography has been used for many years to provide significant information about the internal properties of an object, particularly in the medical fraternity. By reconstructing one-dimensional (ID) X-ray images, 2D cross-sections and 3D renders can provide a wealth of information about an object's internal structure. An extension of the methodology is reported here to enable the characterization of a model agglomerate structure. It is demonstrated that methods based on X-ray microtomography offer considerable potential in the validation and utilization of distinct element method simulations also examined.

On the chaotic rotation of a liquid-filled gyrostat via the Melnikov-Holmes-Marsden integral

The non-linear motions of a gyrostat with an axisymmetrical, fluid-filled cavity are investigated. The cavity is considered to be completely filled with an ideal incompressible liquid performing uniform rotational motion. Helmholtz theorem, Euler's angular momentum theorem and Poisson equations are used to develop the disturbed Hamiltonian equations of the motions of the liquid-filled gyrostat subjected to small perturbing moments. The equations are established in terms of a set of canonical variables comprised of Euler angles and the conjugate angular momenta in order to facilitate the application of the Melnikov-Holmes-Marsden (MHM) method to investigate homoclinic/heteroclinic transversal intersections. In such a way, a criterion for the onset of chaotic oscillations is formulated for liquid-filled gyrostats with ellipsoidal and torus-shaped cavities and the results are confirmed via numerical simulations. (c) 2006 Elsevier Ltd. All rights reserved.

Sigmoidal cosine series on the interval

We construct a set of functions, say, psi([r])(n) composed of a cosine function and a sigmoidal transformation gamma(r) of order r > 0. The present functions are orthonormal with respect to a proper weight function on the interval [-1, 1]. It is proven that if a function f is continuous and piecewise smooth on [-1, 1] then its series expansion based on psi([r])(n) converges uniformly to f so long as the order of the sigmoidal transformation employed is 0 < r

Vibration of symmetrically laminated thick super elliptical plates

This paper reports a free vibration analysis of thick plates with rounded corners subject to a free, simply-supported or clamped boundary condition. The plate perimeter is defined by a super elliptic function with a power defining the shape ranging from an ellipse to a rectangle. To incorporate transverse shear deformation, the Reddy third-order plate theory is employed. The energy integrals incorporating shear deformation and rotary inertia are formulated and the p-Ritz procedures are used to derive the governing eigenvalue equation. Numerical examples for plates with different shapes and boundary conditions are solved and their frequency parameters, where possible, are compared with known results. Parametric studies are carried out to show the sensitivities of frequency parameters by varying the geometry, fibre stacking sequence, and boundary condition. (C) 1999 Academic Press.

Cracks of higher modes in Cosserat continua

Rotational degrees of freedom in Cosserat continua give rise to higher fracture modes. Three new fracture modes correspond to the cracks that are surfaces of discontinuities in the corresponding components of independent Cosserat rotations. We develop a generalisation of J- integral that includes these additional degrees of freedom. The obtained path-independent integrals are used to develop a criterion of crack propagation for a special type of failure in layered materials with sliding layers. This fracture propagates as a progressive bending failure of layers – a “bending crack that is, a crack that can be represented as a distribution of discontinuities in the layer bending. This situation is analysed using a 2D Cosserat continuum model. Semi-infinite bending crack normal to layering is considered. The moment stress concentrates along the line that is a continuation of the crack and has a singularity of the power − 1/4. A model of process zone is proposed for the case when the breakage of layers in the process of bending crack propagation is caused by a crack (microcrack in our description) growing across the layer adjacent to the crack tip. This growth is unstable (in the moment-controlled loading), which results in a typical descending branch of moment stress – rotation discontinuity relationship and hence in emergence of a Barenblatt-type process zone at the tip of the bending crack.