A study of a numerical solution of the steady two dimensions Navier-Stokes equations in a constricted channel problem by a compact fourth order method

We present a numerical solution for the steady 2D Navier-Stokes equations using a fourth order compact-type method. The geometry of the problem is a constricted symmetric channel, where the boundary can be varied, via a parameter, from a smooth constriction to one possessing a very sharp but smooth corner allowing us to analyse the behaviour of the errors when the solution is smooth or near singular. The set of non-linear equations is solved by the Newton method. Results have been obtained for Reynolds number up to 500. Estimates of the errors incurred have shown that the results are accurate and better than those of the corresponding second order method. (C) 2002 Elsevier B.V. All rights reserved.

De Sitter Relativity: a New Road to Quantum Gravity?

The Poincar, group generalizes the Galilei group for high-velocity kinematics. The de Sitter group is assumed to go one step further, generalizing Poincar, as the group governing high-energy kinematics. In other words, ordinary special relativity is here replaced by de Sitter relativity. In this theory, the cosmological constant I > is no longer a free parameter, and can be determined in terms of other quantities. When applied to the whole universe, it is able to predict the value of I > and to explain the cosmic coincidence. When applied to the propagation of ultra-high energy photons, it gives a good estimate of the time delay observed in extragalactic gamma-ray flares. It can, for this reason, be considered a new paradigm to approach the quantum gravity problem.

On homoclinic chaos and a control chaos strategy applied to a free-joint nonholonomic manipulator

In this work, the occurrence of chaos (homoclinic scene) is verified in a robotic system with two degrees of freedom by using Poincare-Mel'nikov method. The studied problem was based on experimental results of a two-joint planar manipulator-first joint actuated and the second joint free-that resides in a horizontal plane. This is the simplest model of nonholonomic free-joint manipulators. The purpose of the present study is to verify analytically those results and to suggest a control strategy.

Boundary crisis and transient in a dissipative relativistic standard map

Some dynamical properties for a problem concerning the acceleration of particles in a wave packet are studied. The model is described in terms of a two-dimensional nonlinear map obtained from a Hamiltonian which describes the motion of a relativistic standard map. The phase space is mixed in the sense that there are regular and chaotic regions coexisting. When dissipation is introduced, the property of area preservation is broken and attractors emerge. We have shown that a tiny increase of the dissipation causes a change in the phase space. A chaotic attractor as well as its basin of attraction are destroyed thereby leading the system to experience a boundary crisis. We have characterized such a boundary crisis via a collision of the chaotic attractor with the stable manifold of a saddle fixed point. Once the chaotic attractor is destroyed, a chaotic transient described by a power law with exponent 1 is observed. (C) 2011 Elsevier B.V. All rights reserved.

Dynamical properties of a particle in a wave packet: Scaling invariance and boundary crisis

Some dynamical properties present in a problem concerning the acceleration of particles in a wave packet are studied. The dynamics of the model is described in terms of a two-dimensional area preserving map. We show that the phase space is mixed in the sense that there are regular and chaotic regions coexisting. We use a connection with the standard map in order to find the position of the first invariant spanning curve which borders the chaotic sea. We find that the position of the first invariant spanning curve increases as a power of the control parameter with the exponent 2/3. The standard deviation of the kinetic energy of an ensemble of initial conditions obeys a power law as a function of time, and saturates after some crossover. Scaling formalism is used in order to characterise the chaotic region close to the transition from integrability to nonintegrability and a relationship between the power law exponents is derived. The formalism can be applied in many different systems with mixed phase space. Then, dissipation is introduced into the model and therefore the property of area preservation is broken, and consequently attractors are observed. We show that after a small change of the dissipation, the chaotic attractor as well as its basin of attraction are destroyed, thus leading the system to experience a boundary crisis. The transient after the crisis follows a power law with exponent -2. (C) 2011 Elsevier Ltd. All rights reserved.

Analytical study of the nonlinear behavior of a shape memory oscillator: Part I-primary resonance and free response at low temperatures

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Energy gain induced by boundary crisis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Rapidly varying boundaries in equations with nonlinear boundary conditions. The case of a Lipschitz deformation

We analyze the behavior of solutions of nonlinear elliptic equations with nonlinear boundary conditions of type partial derivative u/partial derivative n + g( x, u) = 0 when the boundary of the domain varies very rapidly. We show that the limit boundary condition is given by partial derivative u/partial derivative n+gamma(x) g(x, u) = 0, where gamma(x) is a factor related to the oscillations of the boundary at point x. For the case where we have a Lipschitz deformation of the boundary,. is a bounded function and we show the convergence of the solutions in H-1 and C-alpha norms and the convergence of the eigenvalues and eigenfunctions of the linearization around the solutions. If, moreover, a solution of the limit problem is hyperbolic, then we show that the perturbed equation has one and only one solution nearby.

Free actions of abelian p-groups on the n-torus

In this work we make some contributions to the theory of actions of abelian p-groups on the n-Torus T-n. Set congruent to Z(pk1)(h1) x Z(pk2)(h2) x...x Z(pkr)(hr), r >= 1, k(1) >= k(2) >=...>= k(r) >= 1, p prime. Suppose that the group H acts freely on T-n and the induced representation on pi(1)(T-n) congruent to Z(n) is faithful and has first Betti number b. We show that the numbers n, p, b, k(i) and h(i) (i = 1,..,r) satisfy some relation. In particular, when H congruent to Z(p)(h), the minimum value of n is phi(p) + b when b >= 1. Also when H congruent to Z(pk1) x Z(p) the minimum value of n is phi(p(k1)) + p - 1 + b for b >= 1. Here phi denotes the Euler function.

VERY RAPIDLY VARYING BOUNDARIES IN EQUATIONS WITH NONLINEAR BOUNDARY CONDITIONS. THE CASE of A NON UNIFORMLY LIPSCHITZ DEFORMATION

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Semilinear parabolic problems in thin domains with a highly oscillatory boundary

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

Effect of magnesium addition on the vacancy-Mg binding free enthalpy in Al-Cr alloys

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Control of retention and fatigue-free characteristics in CaBi4Ti4O15 thin films prepared by chemical method

Ferroelectric CaBi4Ti4O15 (CBTi144) thin films were deposited on Pt/Ti/SiO2/Si substrates by the polymeric precursor method. The films present a single phase of layered-structured perovskite with polar axis orientation after annealing at 700 degrees C for 2 h in static air and oxygen atmosphere. The a/b-axis orientation of the ferroelectric film is considered to be associated with the preferred orientation of the Pt bottom electrode. It is noted that the films annealed in static air showed good polarization fatigue characteristics at least up to 10(10) bipolar pulse cycles and excellent retention properties up to 10(4) s. on the other hand, oxygen atmosphere seems to be crucial in the decrease of both, fatigue and retention characteristics of the capacitors. Independently of the applied electric field, the retained switchable polarization approached a nearly steady-state value after a retention time of 10 s. (C) 2006 Elsevier B.V. All rights reserved.

Aggregation in the generalized transportation problem

Aggregation disaggregation is used to reduce the analysis of a large generalized transportation problem to a smaller one. Bounds for the actual difference between the aggregated objective and the original optimal value are used to quantify the error due to aggregation and estimate the quality of the aggregation. The bounds can be calculated either before optimization of the aggregated problem (a priori) or after (a posteriori). Both types of the bounds are derived and numerically compared. A computational experiment was designed to (a) study the correlation between the bounds and the actual error and (b) quantify the difference of the error bounds from the actual error. The experiment shows a significant correlation between some a priori bounds, the a posteriori bounds and the actual error. These preliminary results indicate that calculating the a priori error bound is a useful strategy to select the appropriate aggregation level, since the a priori bound varies in the same way that the actual error does. After the aggregated problem has been selected and optimized, the a posteriori bound provides a good quantitative measure for the error due to aggregation.

Self-consistent equilibrium calculation through a direct variational technique in tokamak plasmas

A self-consistent equilibrium calculation, valid for arbitrary aspect ratio tokamaks, is obtained through a direct variational technique that reduces the equilibrium solution, in general obtained from the 2D Grad-Shafranov equation, to a 1D problem in the radial flux coordinate rho. The plasma current profile is supposed to have contributions of the diamagnetic, Pfirsch-Schluter and the neoclassical ohmic and bootstrap currents. An iterative procedure is introduced into our code until the flux surface averaged toroidal current density (J(T)), converges to within a specified tolerance for a given pressure profile and prescribed boundary conditions. The convergence criterion is applied between the (J(T)) profile used to calculate the equilibrium through the variational procedure and the one that results from the equilibrium and given by the sum of all current components. The ohmic contribution is calculated from the neoclassical conductivity and from the self-consistently determined loop voltage in order to give the prescribed value of the total plasma current. The bootstrap current is estimated through the full matrix Hirshman-Sigmar model with the viscosity coefficients as proposed by Shaing, which are valid in all plasma collisionality regimes and arbitrary aspect ratios. The results of the self-consistent calculation are presented for the low aspect ratio tokamak Experimento Tokamak Esferico. A comparison among different models for the bootstrap current estimate is also performed and their possible Limitations to the self-consistent calculation is analysed.