Shallow water flow with cylindrical symmetry

A numerical scheme is presented tor the solution of the shallow water equations in a single radial coordinate. This can prove useful when testing codes for the two-dimensional shallow water equations. The scheme is applied with success to problems involving converging and diverging bores.

Optimal attitude corrections for cylindrical spacecraft

A first order analytical model for optimal small amplitude attitude maneuvers of spacecraft with cylindrical symmetry in an elliptical orbits is presented. The optimization problem is formulated as a Mayer problem with the control torques provided by a power limited propulsion system. The state is defined by Seffet-Andoyer's variables and the control by the components of the propulsive torques. The Pontryagin Maximum Principle is applied to the problem and the optimal torques are given explicitly in Serret-Andoyer's variables and their adjoints. For small amplitude attitude maneuvers, the optimal Hamiltonian function is linearized around a reference attitude. A complete first order analytical solution is obtained by simple quadrature and is expressed through a linear algebraic system involving the initial values of the adjoint variables. A numerical solution is obtained by taking the Euler angles formulation of the problem, solving the two-point boundary problem through the shooting method, and, then, determining the Serret-Andoyer variables through Serret-Andoyer transformation. Numerical results show that the first order solution provides a good approximation to the optimal control law and also that is possible to establish an optimal control law for the artificial satellite's attitude. (C) 2003 COSPAR. Published by Elsevier B.V. Ltd. All rights reserved.

Spherical, cylindrical and tetrahedral symmetries; hydrogenic states at high magnetic field in Si:P

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Equation of state for the magnetic-color-flavor-locked phase and its implications for compact star models

Using the solutions of the gap equations of the magnetic-color-flavor-locked (MCFL) phase of paired quark matter in a magnetic field, and taking into consideration the separation between the longitudinal and transverse pressures due to the field-induced breaking of the spatial rotational symmetry, the equation of state of the MCFL phase is self-consistently determined. This result is then used to investigate the possibility of absolute stability, which turns out to require a field-dependent ""bag constant"" to hold. That is, only if the bag constant varies with the magnetic field, there exists a window in the magnetic field vs bag constant plane for absolute stability of strange matter. Implications for stellar models of magnetized (self-bound) strange stars and hybrid (MCFL core) stars are calculated and discussed.

Nonaxisymmetric magnetorotational instability in ideal and viscous plasmas

The excitation of magnetorotational instability (MRI) in rotating laboratory plasmas is investigated. In contrast to astrophysical plasmas, in which gravitation plays an important role, in laboratory plasmas it can be neglected and the plasma rotation is equilibrated by the pressure gradient. The analysis is restricted to the simple model of a magnetic confinement configuration with cylindrical symmetry, in which nonaxisymmetric perturbations are investigated using the local approximation. Starting from the simplest case of an ideal plasma, the corresponding dispersion relations are derived for more complicated models including the physical effects of parallel and perpendicular viscosities. The Friemann-Rotenberg approach used for ideal plasmas is generalized for the viscous model and an analytical expression for the instability boundary is obtained. It is shown that, in addition to the standard effect of radial derivative of the rotation frequency (the Velikhov effect), which can be destabilizing or stabilizing depending on the sign of this derivative in the ideal plasma, there is a destabilizing effect proportional to the fourth power of the rotation frequency, or, what is the same, to the square of the plasma pressure gradient, and to the square of the azimuthal mode number of the perturbations. It is shown that the instability boundary also depends on the product of the plasma pressure and density gradients, which has a destabilizing effect when it is negative. In the case of parallel viscosity, the MRI looks like an ideal instability independent of viscosity, while, in the case of strong perpendicular viscosity, it is a dissipative instability with the growth rate inversely proportional to the characteristic viscous decay rate. We point out, however, that the modes of the continuous range of the magnetohydrodynamics spectrum are not taken into account in this paper, and they can be more dangerous than those that are considered. (c) 2008 American Institute of Physics.

analytical model and measurements of the target erosion depth profile of balanced and unbalanced planar magnetron cathodes

The erosion depth profile of planar targets in balanced and unbalanced magnetron cathodes with cylindrical symmetry is measured along the target radius. The magnetic fields have rotational symmetry. The horizontal and vertical components of the magnetic field B are measured at points above the cathode target with z = 2 x 10(-3) m. The experimental data reveal that the target erosion depth profile is a function of the angle. made by B with a horizontal line defined by z = 2 x 10(-3) m. To explain this dependence a simplified model of the discharge is developed. In the scope of the model, the pathway lengths of the secondary electrons in the pre-sheath region are calculated by analytical integration of the Lorentz differential equations. Weighting these lengths by using the distribution law of the mean free path of the secondary electrons, we estimate the densities of the ionizing events over the cathode and the relative flux of the sputtered atoms. The expression so deduced correlates for the first time the erosion depth profile of the target with the angle theta. The model shows reasonably good fittings to the experimental target erosion depth profiles confirming that ionization occurs mainly in the pre-sheath zone.

Cylindrically symmetric models of gravitational collapse to black holes: a short review

We survey results about exact cylindrically symmetric models of gravitational collapse in General Relativity. We focus on models which result from the matching of two spacetimes having collapsing interiors which develop trapped surfaces and vacuum exteriors containing gravitational waves. We collect some theorems from the literature which help to decide a priori about eventual spacetime matchings. We revise, in more detail, some toy models which include some of the main mathematical and physical issues that arise in this context, and compute the gravitational energy flux through the matching boundary of a particular collapsing region. Along the way, we point out several interesting open problems.

General clas of inhomogeneous perfect-fluid solutions

We present a general class of solutions to Einstein's field equations with two spacelike commuting Killing vectors by assuming the separation of variables of the metric components. The solutions can be interpreted as inhomogeneous cosmological models. We show that the singularity structure of the solutions varies depending on the different particular choices of the parameters and metric functions. There exist solutions with a universal big-bang singularity, solutions with timelike singularities in the Weyl tensor only, solutions with singularities in both the Ricci and the Weyl tensors, and also singularity-free solutions. We prove that the singularity-free solutions have a well-defined cylindrical symmetry and that they are generalizations of other singularity-free solutions obtained recently.

Dynamics of laser produced silver plasma under film deposition conditions studied using optical emission spectroscopy

Laser produced plasma from silver is generated using a Q-switched Nd:YAG laser. Optical emission spectroscopy is used to carry out time of flight (TOF) analysis of atomic particles. An anomalous double peak profile in the TOF distribution is observed at low pressure. A collection of slower species emerge at reduced pressure below 4 X lO-3 mbar and this species has a greater velocity spread. At high pressure the plasma expansion follows the shockwave model with cylindrical symmetry whereas at reduced pressure it shows unsteady adiabatic expansion (UAE). During UAE the species show a parabolic increases in the expansion time with radial distance whereas during shock wave expansion the exponent is less than one. The angular distribution of the ablated species in the plume is obtained from the measurement of optical density of thin films deposited on to glass substrates kept perpendicular to the plume. There is a sharp variation in the film thickness away from the film centre due to asymmetries in the plume.

A kinematically distorted flux rope model for magnetic clouds

Constant-α force-free magnetic flux rope models have proven to be a valuable first step toward understanding the global context of in situ observations of magnetic clouds. However, cylindrical symmetry is necessarily assumed when using such models, and it is apparent from both observations and modeling that magnetic clouds have highly noncircular cross sections. A number of approaches have been adopted to relax the circular cross section approximation: frequently, the cross-sectional shape is allowed to take an arbitrarily chosen shape (usually elliptical), increasing the number of free parameters that are fit between data and model. While a better “fit” may be achieved in terms of reducing the mean square error between the model and observed magnetic field time series, it is not always clear that this translates to a more accurate reconstruction of the global structure of the magnetic cloud. We develop a new, noncircular cross section flux rope model that is constrained by observations of CMEs/ICMEs and knowledge of the physical processes acting on the magnetic cloud: The magnetic cloud is assumed to initially take the form of a force-free flux rope in the low corona but to be subsequently deformed by a combination of axis-centered self-expansion and heliocentric radial expansion. The resulting analytical solution is validated by fitting to artificial time series produced by numerical MHD simulations of magnetic clouds and shown to accurately reproduce the global structure.

Approximate Riemann solutions of the two-dimensional shallow-water equations

A finite-difference scheme based on flux difference splitting is presented for the solution of the two-dimensional shallow-water equations of ideal fluid flow. A linearised problem, analogous to that of Riemann for gasdynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearised problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. An extension to the two-dimensional equations with source terms, is included. The scheme is applied to a dam-break problem with cylindrical symmetry.

An efficient flux difference splitting algorithm for unsteady duct flows of a real gas

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a real gas in a single spatial coordinate. This includes flow in a duct of variable cross-section, as well as flow with slab, cylindrical or spherical symmetry, as well as the case of an ideal gas, and can be useful when testing codes for the two-dimensional equations governing compressible flow of a real gas. The resulting scheme requires an average of the flow variables across the interface between cells, and this average is chosen to be the arithmetic mean for computational efficiency, which is in contrast to the usual “square root” averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and for a number of equations of state. The results compare favourably with the results from other schemes.

An efficient shock capturing algorithm for compressible flows in a duct of variable cross-section

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a gas in a single spatial co-ordinate. This includes flow in a duct of variable cross-section as well as flow with slab, cylindrical or spherical symmetry and can prove useful when testing codes for the two-dimensional equations governing compressible flow of a gas. The resulting scheme requires an average of the flow variables across the interface between cells and for computational efficiency this average is chosen to be the arithmetic mean, which is in contrast to the usual ‘square root’ averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and a comparison is made in the cylindrical case with results from a two-dimensional problem with no sources.

An efficient numerical scheme for the two-dimensional shallow water equations using arithmetic averaging

A finite difference scheme based on flux difference splitting is presented for the solution of the two-dimensional shallow water equations of ideal fluid flow. A linearised problem, analogous to that of Riemann for gas dynamics is defined, and a scheme, based on numerical characteristic decomposition is presented for obtaining approximate solutions to the linearised problem, and incorporates the technique of operator splitting. An average of the flow variables across the interface between cells is required, and this average is chosen to be the arithmetic mean for computational efficiency leading to arithmetic averaging. This is in contrast to usual ‘square root’ averages found in this type of Riemann solver, where the computational expense can be prohibitive. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids nonphysical, spurious oscillations. An extension to the two-dimensional equations with source terms is included. The scheme is applied to the one-dimensional problems of a breaking dam and reflection of a bore, and in each case the approximate solution is compared to the exact solution of ideal fluid flow. The scheme is also applied to a problem of stationary bore generation in a channel of variable cross-section. Finally, the scheme is applied to two other dam-break problems, this time in two dimensions with one having cylindrical symmetry. Each approximate solution compares well with those given by other authors.