RELATIVISTIC 3-PARTICLE SCATTERING EQUATIONS

We derive a set of relativistic three-particle scattering equations in the three-particle c.m. frame employing a relativistic three-particle propagator suggested long ago by Ahmadzadeh and Tjon in the c.m. frame of a two-particle subsystem. We make the coordinate transformation of this propagator from the c.m. frame of the two-particle subsystem to the three-particle c.m. frame. We also point out that some numerical applications of the Ahmadzadeh and Tjon propagator to the three-nucleon problem use unnecessary nonrelativistic approximations which do not simplify the computational task, but violate constraints of relativistic unitarity and/or covariance.

A W-matrix methodology for solving sparse network equations on multiprocessor computers

This paper describes a methodology for solving efficiently the sparse network equations on multiprocessor computers. The methodology is based on the matrix inverse factors (W-matrix) approach to the direct solution phase of A(x) = b systems. A partitioning scheme of W-matrix , based on the leaf-nodes of the factorization path tree, is proposed. The methodology allows the performance of all the updating operations on vector b in parallel, within each partition, using a row-oriented processing. The approach takes advantage of the processing power of the individual processors. Performance results are presented and discussed.

Analytical attitude propagation with non-singular variables and gravity gradient torque for spin stabilized satellite

An analytical approach for the spin stabilized satellite attitude propagation is presented using the non-singular canonical variables to describe the rotational motion. Two sets of variables were introduced for Fukushima in 1994 by a canonical transformation and they are useful when the angle between z-satellite axis of a coordinate system fixed in artificial satellite and the rotational angular momentum vector is zero or when the angle between Z-equatorial axis and rotation angular momentum vector is zero. Analytical solutions for rotational motion equations and torque-free motion are discussed in terms of the elliptic functions and by the application of some simplification to get an approximated solution. These solutions are compared with a numerical solution and the results show a good agreement for many rotation periods. When the mean Hamiltonian associated with the gravity gradient torque is included, an analytical solution is obtained by the application of the successive approximations' method for the satellite in an elliptical orbit. These solutions show that the magnitude of the rotation angular moment is not affected by the gravity gradient torque but this torque causes linear and periodic variations in the angular variables, long and short periodic variations in Z-equatorial component of the rotation angular moment and short periodic variations in x-satellite component of the rotation angular moment. The goal of this analysis is to emphasize the geometrical and physical meaning of the non-singular variables and to validate the approximated analytical solution for the rotational motion without elliptic functions for a non-symmetrical satellite. The analysis can be applied for spin stabilized satellite and in this case the general solution and the approximated solution are coincidence. Then the results can be used in analysis of the space mission of the Brazilian Satellites. (C) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.

NUMERICAL STUDY OF A PERTURBED EINSTEIN-DESITTER COSMOLOGICAL MODEL

Given the ever-increasing scale of structures discovered in the universe, we solve Einstein's equations numerically, under simplifying assumptions, to examine how this lack of uniformity affects the metric of Einstein-de Sitter cosmology. The results confirm the qualitative conclusion of Barrow, that a large density contrast is compatible with much smaller metric perturbations. The contribution of this peculiar gravity to the redshift might complicate studies of peculiar motions of galaxies, although it appears that the distortion is small for nearby clusters.

On some classes of exactly-solvable Klein-Gordon equations

In this work we discuss some exactly solvable Klein-Gordon equations. We basically discuss the existence of classes of potentials with different nonrelativistic limits, but which shares the intermediate effective Schroedinger differential equation. We comment about the possible use of relativistic exact solutions as approximations for nonrelativistic inexact potentials. (c) 2005 Elsevier B.V. All rights reserved.

MIXMASTER NUMERICAL BEHAVIOR AND GENERALIZATIONS

The numerical simulation of the mixmaster universe serves the purpose of suggesting two kinds of results. The intrinsic time evolution, during contraction, will be seen to be nonchaotic. This is a necessary feature of relativistic cosmological models undergoing this kind of motion. The mixmaster model also provides a clue on how to define chaoticity for systems described by nonautonomous sets of differential equations.

Spatial distribution of solutes and water in sucrose solution dehydrated apples

Half-fresh apples were immersed in sucrose solution (50% w/w, 27 degrees C) during different times of exposition (2, 4, and 8 h). Then each fruit was sliced from the transversal exposed surface. Density, water, and sugar content were determined for each slice. A mathematical model was fitted to experimental data of water and sucrose content considering the global flux and the tissue shrinkage. By numerical analysis, the binary effective diffusion coefficients as a function of concentration were calculated, using material coordinates and integrating simultaneously two differential equations (for water and sucrose). The coefficients obtained are one or even two orders of magnitude lower than the ones for pure solutions and present an unusual concentration dependence. This comparison shows the influence of the tissue resistance to the diffusion.

Numerical simulation of viscous three-dimensional flow over aerospace vehicles using Euler/boundary-layer zonal approach

This paper presents a viscous three-dimensional simulations coupling Euler and boundary layer codes for calculating flows over arbitrary surfaces. The governing equations are written in a general non orthogonal coordinate system. The Levy-Lees transformation generalized to three-dimensional flows is utilized. The inviscid properties are obtained from the Euler equations using the Beam and Warming implicit approximate factorization scheme. The resulting equations are discretized and approximated by a two-point fmitedifference numerical scheme. The code developed is validated and applied to the simulation of the flowfield over aerospace vehicle configurations. The results present good correlation with the available data.

Tau-functions and dressing transformations for zero-curvature affine integrable equations

The solutions of a large class of hierarchies of zero-curvature equations that includes Toda- and KdV-type hierarchies are investigated. All these hierarchies are constructed from affine (twisted or untwisted) Kac-Moody algebras g. Their common feature is that they have some special vacuum solutions corresponding to Lax operators lying in some Abelian (up to the central term) subalgebra of g; in some interesting cases such subalgebras are of the Heisenberg type. Using the dressing transformation method, the solutions in the orbit of those vacuum solutions are constructed in a uniform way. Then, the generalized tau-functions for those hierarchies are defined as an alternative set of variables corresponding to certain matrix elements evaluated in the integrable highest-weight representations of g. Such definition of tau-functions applies for any level of the representation, and it is independent of its realization (vertex operator or not). The particular important cases of generalized mKdV and KdV hierarchies as well as the Abelian and non-Abelian affine Toda theories are discussed in detail. © 1997 American Institute of Physics.

Infrared finite solutions for the gluon propagator and the QCD vacuum energy

Nonperturbative infrared finite solutions for the gluon polarization tensor have been found, and the possibility that gluons may have a dynamically generated mass is supported by recent Monte Carlo simulation on the lattice. These solutions differ among themselves, due to different approximations performed when solving the Schwinger-Dyson equations for the gluon polarization tensor. Only approximations that minimize energy are meaningful, and, according to this, we compute an effective potential for composite operators as a function of these solutions in order to distinguish which one is selected by the vacuum. © 1997 Elsevier Science B.V.

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.

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.

Numerical study of the spherically symmetric Gross-Pitaevskii equation in two space dimensions

The Gross-Pitaevskii equation for Bose-Einstein condensation (BEC) in two space dimensions under the action of a harmonic oscillator trap potential for bosonic atoms with attractive and repulsive interparticle interactions was numerically studied by using time-dependent and time-independent approaches. In both cases, numerical difficulty appeared for large nonlinearity. Nonetheless, the solution of the time-dependent approach exhibited intrinsic oscillation with time iteration which is independent of space and time steps used in discretization.

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.

Confinement and soliton solutions in the SL(3) Toda model coupled to matter fields

We consider an integrable conformally invariant two-dimensional model associated to the affine Kac-Moody algebra sl3(ℂ). It possesses four scalar fields and six Dirac spinors. The theory does not possesses a local Lagrangian since the spinor equations of motion present interaction terms which are bilinear in the spinors. There exists a submodel presenting an equivalence between a U(1) vector current and a topological current, which leads to a confinement of the spinors inside the solitons. We calculate the one-soliton and two-soliton solutions using a procedure which is a hybrid of the dressing and Hirota methods. The soliton masses and time delays due to the soliton interactions are also calculated. We give a computer program to calculate the soliton solutions. © 2002 Published by Elsevier Science B.V.