A combined wavelet-element free Galerkin method for numerical calculations of electromagnetic fields

A combined wavelet-element free Galerkin (EFG) method is proposed for solving electromagnetic EM) field problems. The bridging scales are used to preserve the consistency and linear independence properties of the entire bases. A detailed description of the development of the discrete model and its numerical implementations is given to facilitate the reader to. understand the proposed algorithm. A numerical example to validate the proposed method is also reported.

Reformulation of variational inequalities on a simplex and compactification of complementarity problems

Many variational inequality problems (VIPs) can be reduced, by a compactification procedure, to a VIP on the canonical simplex. Reformulations of this problem are studied, including smooth reformulations with simple constraints and unconstrained reformulations based on the penalized Fischer-Burmeister function. It is proved that bounded level set results hold for these reformulations under quite general assumptions on the operator. Therefore, it can be guaranteed that minimization algorithms generate bounded sequences and, under monotonicity conditions, these algorithms necessarily nd solutions of the original problem. Some numerical experiments are presented.

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.

NUMERICAL STUDY OF TRANSPORT IN A DISSIPATIVE MEDIUM

A numerical study of propagation of a particle through a one-dimensional dissipative medium is presented. The medium is composed of several dissipative sections, which are characterized by their friction coefficients eta. In particular, we have considered two types of friction coefficients distributed orderly or disorderly along the chain. For the same relative proportion of the coefficients, we have found that transport can be enhanced in the disordered distribution in comparison with the ordered one. We also show how this can be considered an approximated way to treat the propagation in a dissipative medium with a position-dependent friction coefficient.

An analytical and numerical study of plane change maneuvers using aerodynamic force

The goal of the present work is to analyze space missions that use the terrestrial atmosphere to accomplish orbital maneuvers that involve a plane change. A set of analytical solutions is presented for the variation of the orbital elements due to a single passage through the atmosphere, assuming that the interval the spacecraft travels through the atmosphere is not too large. The study considers both the lift influence on the spacecraft orbit as well as drag. The final equations are tested with numerical integration and can be considered in accordance with the numerical results whenever the perigee height is larger than a critical value. Next, a numerical study of the ratio between the velocity increment required to correct the semimajor axis decay due to the atmospheric passage and the velocity variation required to obtain the change in the inclination is also presented. This analysis can be used to decide if a maneuver passing through the atmosphere can decrease the fuel consumption of the mission and, in the cases where this technique can be used, if a multiple passage is more efficient than a single passage.

Non-linear boundary element analysis of floor slabs reinforced with rectangular beams

In this work, a numerical model to perform non-linear analysis of building floor structures is proposed. The presented model is derived from the Kirchhoff-s plate bending formulation of the boundary element method (BENI) for zoned domains, in which the plate stiffness is modified by the presence of membrane effects. In this model, no approximation of the generalized forces along the interface is required and the compatibility and equilibrium conditions along interfaces are imposed at the integral equation level. In order to reduce the number of degrees of freedom, the Navier Bernoulli hypothesis is assumed to simplify the strain field for the thin sub-regions (rectangular beams). The non-linear formulation is obtained from the linear formulation by incorporating initial internal force fields, which are approximated by using the well-known cell sub-division. Then, the non-linear solution of algebraic equations is obtained by using the concept of the consistent tangent operator. The Von Mises criterion is adopted to govern the elasto-plastic material behaviour checked at points along the plate thickness and along the rectangular beam element axes. The numerical representations are accurately obtained by either computing analytically the element integrals or performing the numerical integration accurately using an appropriate sub-elementation scheme. (C) 2007 Elsevier Ltd. All rights reserved.

NUMERICAL DETERMINATION OF THE ORDER OF PHASE-TRANSITIONS IN ISING SYSTEMS WITH MULTISPIN INTERACTIONS

The method of the fourth-order cumulant of Challa, Landau, and Binder is used together with the Monte Carlo histogram technique of Ferrenberg and Swendsen to study the order of the phase transitions of two-dimensional Ising systems with multispin interactions in the horizontal direction and two-body interactions in the vertical direction.

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.

PHYLOGENETIC STUDIES OF SOME SPECIES OF THE GENUS COFFEA .2. NUMERICAL-ANALYSIS OF ISOENZYMATIC DATA

Thirteen species of Coffea were studied for five enzymes systems, including alpha and beta esterase, alkaline phosphatase, acid phosphatase, malate dehydrogenase and acid dehydrogenase. Three coefficients of similarity: Simple Matching, Jaccard and Ochiai and three different clustering methods: Single Linkage, Complete Linkage and Unweighted Pair Group, using Arithmetic Averages (UPGMA) were used to analyse the data.The phylogenetic relationships among the twelve diploid species and between them and the tetraploid species C. arabica showed that similarity among species of the same subsection is not always greater than among species of different subsections. In addition, although there are several similarity groups in common, established by isoenzymatic polymorphism, morphological characteristics, chemical data, crossability and geographic distribution, there is no common trend among the phylogenetic relationships as indicated by all these different evaluating procedures.

ITERATIVE NUMERICAL-SOLUTION OF SCATTERING PROBLEMS

An iterative Neumann series method, employing a real auxiliary scattering integral equation, is used to calculate scattering lengths and phase shifts for the atomic Yukawa and exponential potentials. For these potentials the original Neumann series diverges. The present iterative method yields results that are far better, in convergence, stability and precision, than other momentum space methods. Accurate result is obtained in both cases with an estimated error of about 1 in 10(10) after some 8-10 iterations.

Numerical study of low-cost alternative orbits around the Moon

In this paper, we have investigated a region of direct stable orbits around the Moon, whose stability is related to the H2 Family of periodic orbits and to the quasi-periodic orbits that oscillate around them. The stability criteria adopted was that the path did not escape from the Moon during an integration period of 1000 days (remaining with negative two-body Moon-probe orbital energy during this period). Considering the three-dimensional four-body Sun-Earth-Moon-probe problem, we investigated the evolution of the size of the stability region, taking into account the eccentricity of the Earth's orbit, the eccentricity and inclination of the Moon's orbit, and the solar radiation pressure on the probe. We also investigated the evolution of the region's size and its location by varying the inclination of the probe's initial osculating orbit relative to the Moon's orbital plane between 0 degrees and 180 degrees. The size of the stability region diminishes; nevertheless, it remains significant for 0 <= i <= 25 degrees and 35 degrees <= i <= 45 degrees. The orbits of this region could be useful for missions by space vehicles that must remain in orbit around the Moon for periods of up to 1000 days, requiring low maintenance costs. (c) 2005 COSPAR. Published by Elsevier Ltd. All rights reserved.

Numerical study of circulation on the inner Amazon Shelf

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

Integration Methods Used in Numerical Simulations of Electromagnetic Transients

This paper made an analysis of some numerical integration methods that can be used in electromagnetic transient simulations. Among the existing methods, we analyzed the trapezoidal integration method (or Heun formula), Simpson's Rule and Runge-Kutta. These methods were used in simulations of electromagnetic transients in power systems, resulting from switching operations and maneuvers that occur in transmission lines. Analyzed the characteristics such as accuracy, computation time and robustness of the methods of integration.