Fourth-order method for solving the Navier-Stokes equations in a constricting channel

A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.

Some results on stability of retarded functional differential equations using dichotomic map techniques

This paper deals with the study of the stability of nonautonomous retarded functional differential equations using the theory of dichotomic maps. After some preliminaries, we prove the theorems on simple and asymptotic stability. Some examples are given to illustrate the application of the method. Main results about asymptotic stability of the equation x′(t) = -b(t)x(t - r) and of its nonlinear generalization x′(t) = b(t) f (x(t - r)) are established. © 1998 Kluwer Academic Publishers.

Solid mass falling into the reservoirs shallow water equations validity from the engineering point of view

Water waves generated by a solid mass is a complex phenomenon discussed in this paper by numerical and experimental approaches. A model based on shallow water equations with shocks (Saint Venant) has developed. It can reproduce the amplitude and the energy of the wave quite well, but because it consistently generates a hydraulic jump, it is able to reproduce the profile, in the case of high relative thickness of slide, but in the case of small relative thickness it is unable to reproduce the amplitude of the wave. As the momentum conservation is not verified during the phase of wave creation, a second technique based on discharge transfer coefficient α, is introduced at the zone of impact. Numerical tests have been performed and validated this technique from the experimental results of the wave's height obtained in a flume.

Unitary Operator Bases and Q-deformed Algebras

Starting from the Schwinger unitary operator bases formalism constructed out of a finite dimensional state space, the well-known q-deformed commutation relation is shown to emerge in a natural way, when the deformation parameter is a root of unity.

A note on the superstring BRST operator

We write the BRST operator of the N = 1 superstring as, Q = e-R(1/2πiφdzγ2b)eR where y and b are super-reparameterization ghosts. This provides a trivial proof that Q is nilpotent. © 1999 Published by Elsevier Science B.V. All rights reserved.

Equivalence of the Duffin-Kemmer-Petiau and Klein-Gordon-Fock equations

A strict proof of the equivalence of the Duffin-Kemmer-Petiau and Klein-Gordon Fock theories is presented for physical S-matrix elements in the case of charged scalar particles minimally interacting with an external or quantized electromagnetic field. The Hamiltonian canonical approach to the Duffin - Kemmer Petiau theory is first developed in both the component and the matrix form. The theory is then quantized through the construction of the generating functional for the Green's functions, and the physical matrix elements of the S-matrix are proved to be relativistic invariants. The equivalence of the two theories is then proved for the matrix elements of the scattered scalar particles using the reduction formulas of Lehmann, Symanzik, and Zimmermann and for the many-photon Green's functions.

Lamé differential equations and electrostatics

The problem of existence and uniqueness of polynomial solutions of the Lamé differential equation A(x)y″ + 2B(x)y′ + C(x)y = 0, where A(x),B(x) and C(x) are polynomials of degree p + 1,p and p - 1, is under discussion. We concentrate on the case when A(x) has only real zeros aj and, in contrast to a classical result of Heine and Stieltjes which concerns the case of positive coefficients rj in the partial fraction decomposition B(x)/A(x) = ∑j p=0 rj/(x - aj), we allow the presence of both positive and negative coefficients rj. The corresponding electrostatic interpretation of the zeros of the solution y(x) as points of equilibrium in an electrostatic field generated by charges rj at aj is given. As an application we prove that the zeros of the Gegenbauer-Laurent polynomials are the points of unique equilibrium in a field generated by two positive and two negative charges. © 2000 American Mathematical Society.

Galilean covariance and non-relativistic Bhabha equations

We apply a five-dimensional formulation of Galilean covariance to construct non-relativistic Bhabha first-order wave equations which, depending on the representation, correspond either to the well known Dirac equation (for particles with spin 1/2) or the Duffin-Kemmer-Petiau equation (for spinless and spin 1 particles). Here the irreducible representations belong to the Lie algebra of the 'de Sitter group' in 4 + 1 dimensions, SO(5, 1). Using this approach, the non-relativistic limits of the corresponding equations are obtained directly, without taking any low-velocity approximation. As a simple illustration, we discuss the harmonic oscillator.

Massive superstring vertex operator in D = 10 superspace

Using the pure spinor formalism for the superstring, the vertex operator for the first massive states of the open superstring is constructed in a manifestly super-Poincaré covariant manner. This vertex operator describes a massive spin-two multiplet in terms of ten-dimensional superfields. © SISSA/ISAS 2002.

Solutions for a class of integro-differential equations with time periodic coefficients

In this work, a series solution is found for the integro-differential equation y″ (t) = -(ω2 c + ω2 f sin2 ωpt)y(t) + ωf (sin ωpt) z′ (0) + ω2 fωp sin ωpt ∫t 0 (cos ωps) y(s)ds, which describes the charged particle motion for certain configurations of oscillating magnetic fields. As an interesting feature, the terms of the solution are related to distinct sequences of prime numbers.

A FEM approach to the equations of magneto-aerodynamics

In this work, a Finite Element Method treatment is outlined for the equations of Magnetoaerodynamics. In order to provide a good basis for numerical treatment of Magneto-aerodynamics, a full version of the complete equations is presented and FEM contribution matrices are deduced, as well as further terms of stabilization for the compressible flow case.

Supersymmetry and the KdV equations for integrable hierarchies with a half-integer gradation

Supersymmetry is formulated for integrable models based on the sl(2 1) loop algebra endowed with a principal gradation. The symmetry transformations which have half-integer grades generate supersymmetry. The sl(2 1) loop algebra leads to N=2 supersymmetric mKdV and sinh-Gordon equations. The corresponding N=1 mKdV and sinh-Gordon equations are obtained via reduction induced by twisted automorphism. Our method allows for a description of a non-local symmetry structure of supersymmetric integrable models. © 2003 Elsevier B.V. All rights reserved.

A dynamical gluon mass solution in a coupled system of the Schwinger-Dyson equations

We study numerically the Schwinger-Dyson equations for the coupled system of gluon and ghost propagators in the Landau gauge and in the case of pure gauge QCD. We show that a dynamical mass for the gluon propagator arises as a solution while the ghost propagator develops an enhanced behavior in the infrared regime of QCD. Simple analytical expressions are proposed for the propagators, and the mass dependency on the ΛQCD scale and its perturbative scaling are studied. We discuss the implications of our results for the infrared behavior of the coupling constant, which, according to fits for the propagators infrared behavior, seems to indicate that α s(q2) → 0 as q2 → 0. © SISSA/ISAS 2004.

Schwinger-dyson equations and dynamical gluon mass generation

We discuss the solutions obtained for the gluon propagador in Landau gauge within two distinct approximations for the Schwinger-Dyson equations (SDE). The first, named Mandelstam's approximation, consist in neglecting all contributions that come from fermions and ghosts fields while in the second, the ghosts fields are taken into account leading to a coupled system of integral equations. In both cases we show that a dynamical mass for the gluon propagator can arise as a solution. © 2005 American Institute of Physics.

Substitution-newton-Raphson method for the solution of electric network equations

The power flow problem, in transmission networks, has been well solved, for most cases, using Newton-Raphson method (NR) and its decoupled versions. Generally speaking, the solution of a non-linear system of equations refers to two methods: NR and Successive Substitution. The proposal of this paper is to evaluate the potential of the Substitution-Newton-Raphson Method (SNR), which combines both methods, on the solution of the power flow problem. Simulations were performed using a two-bus test network in order to observe the characteristics of these methods. It was verified that the NR is faster than SNR, in terms of convergence, considering non-stressed scenarios. For those cases where the power flow in the network is closed to the limits (stressed system), the SNR converges faster. This paper presents the power flow formulation of the SNR and describes its potential for its application in special cases such as stressed scenarios. © 2006 IEEE.