Search on a hypercubic lattice using a quantum random walk. I. d > 2

Random walks describe diffusion processes, where movement at every time step is restricted to only the neighboring locations. We construct a quantum random walk algorithm, based on discretization of the Dirac evolution operator inspired by staggered lattice fermions. We use it to investigate the spatial search problem, that is, to find a marked vertex on a d-dimensional hypercubic lattice. The restriction on movement hardly matters for d > 2, and scaling behavior close to Grover's optimal algorithm (which has no restriction on movement) can be achieved. Using numerical simulations, we optimize the proportionality constants of the scaling behavior, and demonstrate the approach to that for Grover's algorithm (equivalent to the mean-field theory or the d -> infinity limit). In particular, the scaling behavior for d = 3 is only about 25% higher than the optimal d -> infinity value.

On the relation between rotation increments in different tangent spaces

In computational mechanics, finite rotations are often represented by rotation vectors. Rotation vector increments corresponding to different tangent: spaces are generally related by a linear operator, known as the tangential transformation T. In this note, we derive the higher order terms that are usually left out in linear relation. The exact nonlinear relation is also presented. Errors via the linearized T are numerically estimated. While the concept of T arises out of the nonlinear characteristics of the rotation manifold, it has been derived via tensor analysis in the context of computational mechanics (Cardona and Geradin, 1988). We investigate the operator T from a Lie group perspective, which provides a better insight and a 1-1 correspondence between approaches based on tensor analysis and the standard matrix Lie group theory. (C) 2010 Elsevier Ltd. All rights reserved.

Effect of wall thickness on the end corrections of the extended inlet and outlet of a double-tuned expansion chamber

Simple expansion chambers, the simplest of the muffler configurations, have very limited practical application due to the presence of periodic troughs in the transmission loss spectrum which drastically lower the overall transmission loss of the muffler. Tuned extended inlet and outlet can be designed to nullify three-fourths of these troughs, making use of the plane wave theory. These cancellations would not occur unless one altered the geometric lengths for the extended tube in order to incorporate the effect of evanescent higher-order modes (multidimensional effect) through end corrections or lumped inertance approximation at the area discontinuities or junctions. End corrections of the extended inlet and outlet have been studied by several researchers. However the effect of wall thickness of the inlet/outlet duct on end correction has not been studied explicitly. This has significant effect on the tuning of an extended inlet/outlet expansion chamber. It is investigated here experimentally as well as numerically (through use of 3-D FEM software) for stationary medium. Crown Copyright (C) 2010 Published by Elsevier Ltd. All rights reserved.

Q-Band beam wave-guide

Consideration is given to a 25-foot long Q-band (8 mm) confocal, zoned dielectric lens beam waveguide. Numerical expressions for the axial and radial fields are presented. The experimental set-up consisted of uniformly spaced zoned dielectric lenses, a transmitting horn and a receiving horn. It was found that: (1) the wave beam is reiterated when confocal, zoned dielectric lenses act as phase transformers in place of smooth surfaced transformers in beam waveguides; (2) the axial field is oscillatory near the source and the oscillation persists for about 25 cm from the source; (3) the oscillation disappears after one lens is used; (4) higher order modes with higher attenuation rates die out faster than fundamental modes; (5) phase transformers do not alter beam modes; (6) without any lens the beam cross-section broadens significantly in the Z-direction; (7) with one lens the beam exhibits the reiteration phenomenon; and (8) inserting a second lens on the axial and cross-sectional field distribution shows further the reiteration principle.

Path integral quantisation of topological field theories

Conditions for quantum topological invariance of classically topological field theories in the path integral formulation are discussed. Both the three-dimensional Chern-Simons system and a Witten-type topological field theory are shown to satisfy these conditions.

A simple method for potential flow simulation of cascades

A simple method using a combination of conformal mapping and vortex panel method to simulate potential flow in cascades is presented. The cascade is first transformed to a single body using a conformal mapping, and the potential flow over this body is solved using a simple higher order vortex panel method. The advantage of this method over existing methodologies is that it enables the use of higher order panel methods, as are used to solve flow past an isolated airfoil, to solve the cascade problem without the need for any numerical integrations or iterations. The fluid loading on the blades, such as the normal force and pitching moment, may be easily calculated from the resultant velocity field. The coefficient of pressure on cascade blades calculated with this methodology shows good agreement with previous numerical and experimental results.

Stability of spatially developing boundary layers in pressure gradients

A new formulation of the stability of boundary-layer flows in pressure gradients is presented, taking into account the spatial development of the flow and utilizing a special coordinate transformation. The formulation assumes that disturbance wavelength and eigenfunction vary downstream no more rapidly than the boundary-layer thickness, and includes all terms nominally of order R(-1) in the boundary-layer Reynolds number R. In Blasius flow, the present approach is consistent with that of Bertolotti et al. (1992) to O(R(-1)) but simpler (i.e. has fewer terms), and may best be seen as providing a parametric differential equation which can be solved without having to march in space. The computed neutral boundaries depend strongly on distance from the surface, but the one corresponding to the inner maximum of the streamwise velocity perturbation happens to be close to the parallel flow (Orr-Sommerfeld) boundary. For this quantity, solutions for the Falkner-Skan flows show the effects of spatial growth to be striking only in the presence of strong adverse pressure gradients. As a rational analysis to O(R(-1)) demands inclusion of higher-order corrections on the mean flow, an illustrative calculation of one such correction, due to the displacement effect of the boundary layer, is made, and shown to have a significant destabilizing influence on the stability boundary in strong adverse pressure gradients. The effect of non-parallelism on the growth of relatively high frequencies can be significant at low Reynolds numbers, but is marginal in other cases. As an extension of the present approach, a method of dealing with non-similar flows is also presented and illustrated. However, inherent in the transformation underlying the present approach is a lower-order non-parallel theory, which is obtained by dropping all terms of nominal order R(-1) except those required for obtaining the lowest-order solution in the critical and wall layers. It is shown that a reduced Orr-Sommerfeld equation (in transformed coordinates) already contains the major effects of non-parallelism.

Weakly non-linear stability of a hydrodynamic accretion disc

When the cold accretion disc coupling between neutral gas and a magnetic field is so weak that the magnetorotational instability is less effective or even stops working, it is of prime interest to investigate the pure hydrodynamic origin of turbulence and transport phenomena. As the Reynolds number increases, the relative importance of the non-linear term in the hydrodynamic equation increases. In an accretion disc where the molecular viscosity is too small, the Reynolds number is large enough for the non-linear term to have new effects. We investigate the scenario of the `weakly non-linear' evolution of the amplitude of the linear mode when the flow is bounded by two parallel walls. The unperturbed flow is similar to the plane Couette flow, but with the Coriolis force included in the hydrodynamic equation. Although there is no exponentially growing eigenmode, because of the self-interaction, the least stable eigenmode will grow in an intermediate phase. Later, this will lead to higher-order non-linearity and plausible turbulence. Although the non-linear term in the hydrodynamic equation is energy-conserving, within the weakly non-linear analysis it is possible to define a lower bound of the energy (alpha A(c)(2), where A(c) is the threshold amplitude) needed for the flow to transform to the turbulent phase. Such an unstable phase is possible only if the Reynolds number >= 10(3-4). The numerical difficulties in obtaining such a large Reynolds number might be the reason for the negative result of numerical simulations on a pure hydrodynamic Keplerian accretion disc.

Toward a Mean Field Formulation of the Babcock-Leighton Type Solar Dynamo. I. α-Coefficient versus Durney's Double-Ring Approach

We develop a model of the solar dynamo in which, on the one hand, we follow the Babcock-Leighton approach to include surface processes, such as the production of poloidal field from the decay of active regions, and, on the other hand, we attempt to develop a mean field theory that can be studied in quantitative detail. One of the main challenges in developing such models is to treat the buoyant rise of the toroidal field and the production of poloidal field from it near the surface. A previous paper by Choudhuri, Schüssler, & Dikpati in 1995 did not incorporate buoyancy. We extend this model by two contrasting methods. In one method, we incorporate the generation of the poloidal field near the solar surface by Durney's procedure of double-ring eruption. In the second method, the poloidal field generation is treated by a positive α-effect concentrated near the solar surface coupled with an algorithm for handling buoyancy. The two methods are found to give qualitatively similar results.

Confinement-deconfinement transition and spin correlations in a generalized Kitaev model

We present a spin model, namely, the Kitaev model augmented by a loop term and perturbed by an Ising Hamiltonian, and show that it exhibits both confinement-deconfinement transitions from spin liquid to antiferromagnetic/spin-chain/ferromagnetic phases and topological quantum phase transitions between gapped and gapless spin-liquid phases. We develop a fermionic resonating-valence-bonds (RVB) mean-field theory to chart out the phase diagram of the model and estimate the stability of its spin-liquid phases, which might be relevant for attempts to realize the model in optical lattices and other spin systems. We present an analytical mean-field theory to study the confinement-deconfinement transition for large coefficient of the loop term and show that this transition is first order within such mean-field analysis in this limit. We also conjecture that in some other regimes, the confinement-deconfinement transitions in the model, predicted to be first order within the mean-field theory, may become second order via a defect condensation mechanism. Finally, we present a general classification of the perturbations to the Kitaev model on the basis of their effect on it's spin correlation functions and derive a necessary and sufficient condition, within the regime of validity of perturbation theory, for the spin correlators to exhibit a long-ranged power-law behavior in the presence of such perturbations. Our results reproduce those of Tikhonov et al. [Phys. Rev. Lett. 106, 067203 (2011)] as a special case.

Melting and mechanical properties of polymer grafted lipid bilayer membranes

The influence of polymer grafting on the phase behavior and elastic properties of two tail lipid bilayers have been investigated using dissipative particle dynamics simulations. For the range of polymer lengths studied, the L(c) to L(alpha) transition temperature is not significantly affected for grafting fractions, G(f) between 0.16 and 0.25. A decrease in the transition temperature is observed at a relatively high grafting fraction, G(f) = 0.36. At low temperatures, a small increase in the area per head group, a(h), at high G(f) leads to an increase in the chain tilt, inducing order in the bilayer and the solvent. The onset of the phase transition occurs with the nucleation of small patches of thinned membrane which grow and form continuous domains as the temperature increases. This region is the co-existence region between the L(beta)(thick) and the L(alpha)(thin) phases. The simulation results for the membrane area expansion as a function of the grafting density conform extremely well to the scalings predicted by self-consistent mean field theories. We find that the bending modulus shows a small decrease for short polymers (number of beads, N(p) = 10) and low G(f), where the influence of polymer is reduced when compared to the effect of the increased a(h). For longer polymers (N(p) > 15), the bending modulus increases monotonically with increase in grafted polymer. Using the results from mean field theory, we partition the contributions to the bending modulus from the membrane and the polymer and show that the dominant contribution to the increased bending modulus arises from the grafted polymer. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3631940]

A simple numerical method for three-dimensional analysis of simple expansion chamber mufflers of rectangular as well as circular cross-section with a stationary medium

Three-dimensional effects are a primary source of discrepancy between the measured values of automotive muffler performance and those predicted by the plane wave theory at higher frequencies. The basically exact method of (truncated) eigenfunction expansions for simple expansion chambers involves very complicated algebra, and the numerical finite element method requires large computation time and core storage. A simple numerical method is presented in this paper. It makes use of compatibility conditions for acoustic pressure and particle velocity at a number of equally spaced points in the planes of the junctions (or area discontinuities) to generate the required number of algebraic equations for evaluation of the relative amplitudes of the various modes (eigenfunctions), the total number of which is proportional to the area ratio. The method is demonstrated for evaluation of the four-pole parameters of rigid-walled, simple expansion chambers of rectangular as well as circular cross-section for the case of a stationary medium. Computed values of transmission loss are compared with those computed by means of the plane wave theory, in order to highlight the onset (cutting-on) of various higher order modes and the effect thereof on transmission loss of the muffler. These are also compared with predictions of the finite element methods (FEM) and the exact methods involving eigenfunction expansions, in order to demonstrate the accuracy of the simple method presented here.

Microwave cavity resonators. Some perturbation effects and their applications

Lagrange's equation is utilized to show the analogy of a lossless microwave cavity resonator with the conventional LC network. A brief discussion on the resonant frequencies of a microwave cavity resonator and the two degenerate companion modes H01 and E11 appearing in a cavity is given. The first order perturbation theory of a small deformation of the wall of a cavity is discussed. The effects of perturbation, such as the change in the resonant frequency and the Q of a cavity, the change in the electromagnetic field configurations and hence mixing of modes are also discussed. An expression for the coupling coefficient between the two degenerate modes H01 and E11 is derived with the help of the field equations. Results indicate that in the absence of perturbation the above two degenerate modes can co-exist without losing their individual identities. Several applications of the perturbation theory, such as the measurement of the dielectric properties of matter, study of ferromagnetic resonance, etc., are described.

Smooth DMS-FEM: A new approach to solving nearly incompressible nonlinear elasto-static problems

The DMS-FEM, which enables functional approximations with C(1) or still higher inter-element continuity within an FEM-based meshing of the domain, has recently been proposed by Sunilkumar and Roy [39,40]. Through numerical explorations on linear elasto-static problems, the method was found to have conspicuously superior convergence characteristics as well as higher numerical stability against locking. These observations motivate the present study, which aims at extending and exploring the DMS-FEM to (geometrically) nonlinear elasto-static problems of interest in solid mechanics and assessing its numerical performance vis-a-vis the FEM. In particular, the DMS-FEM is shown to vastly outperform the FEM (presently implemented through the commercial software ANSYS (R)) as the former requires fewer linearization and load steps to achieve convergence. In addition, in the context of nearly incompressible nonlinear systems prone to volumetric locking and with no special numerical artefacts (e.g. stabilized or mixed weak forms) employed to arrest locking, the DMS-FEM is shown to approach the incompressibility limit much more closely and with significantly fewer iterations than the FEM. The numerical findings are suggestive of the important role that higher order (uniform) continuity of the approximated field variables play in overcoming volumetric locking and the great promise that the method holds for a range of other numerically ill-conditioned problems of interest in computational structural mechanics. (C) 2011 Elsevier Ltd. All rights reserved.

Coupled Wavenumbers in an Infinite Flexible Fluid-Filled Circular Cylindrical Shell: Comparison between Different Shell Theories

Analytical expressions are found for the wavenumbers in an infinite flexible in vacuo I fluid-filled circular cylindrical shell based on different shell-theories using asymptotic methods. Donnell-Mushtari theory (the simplest shell theory) and four higher order theories, namely Love-Timoshenko, Goldenveizer-Novozhilov, Flugge and Kennard-simplified are considered. Initially, in vacuo and fluid-coupled wavenumber expressions are presented using the Donnell-Mushtari theory. Subsequently, the wavenumbers using the higher order theories are presented as perturbations on the Donnell-Mushtari wavenumbers. Similarly, expressions for the resonance frequencies in a finite shell are also presented, using each shell theory. The basic differences between the theories being what they are, the analytical expressions obtained from the five theories allow one to see how these differences propagate into the asymptotic expansions. Also, they help to quantify the difference between the theories for a wide range of parameter values such as the frequency range, circumferential order, thickness ratio of the shell, etc.