Assessment of a high-order finite difference upwind scheme for the simulation of convection-diffusion problems

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

Numerical approximation of the Ginzburg-Landau equation with memory effects in the dynamics of phase transitions

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

Finite-difference calculation of donor energy levels in a spherical quantum dot subject to a magnetic field

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

Mathematical tests about the existence and applications of PPS-wavelets in function approximation problems

Neural networks and wavelet transform have been recently seen as attractive tools for developing eficient solutions for many real world problems in function approximation. Function approximation is a very important task in environments where computation has to be based on extracting information from data samples in real world processes. So, mathematical model is a very important tool to guarantee the development of the neural network area. In this article we will introduce one series of mathematical demonstrations that guarantee the wavelets properties for the PPS functions. As application, we will show the use of PPS-wavelets in pattern recognition problems of handwritten digit through function approximation techniques.

A fully adaptive multiresolution scheme for shock computations

The scheme is based on Ami Harten's ideas (Harten, 1994), the main tools coming from wavelet theory, in the framework of multiresolution analysis for cell averages. But instead of evolving cell averages on the finest uniform level, we propose to evolve just the cell averages on the grid determined by the significant wavelet coefficients. Typically, there are few cells in each time step, big cells on smooth regions, and smaller ones close to irregularities of the solution. For the numerical flux, we use a simple uniform central finite difference scheme, adapted to the size of each cell. If any of the required neighboring cell averages is not present, it is interpolated from coarser scales. But we switch to ENO scheme in the finest part of the grids. To show the feasibility and efficiency of the method, it is applied to a system arising in polymer-flooding of an oil reservoir. In terms of CPU time and memory requirements, it outperforms Harten's multiresolution algorithm.The proposed method applies to systems of conservation laws in 1Dpartial derivative(t)u(x, t) + partial derivative(x)f(u(x, t)) = 0, u(x, t) is an element of R-m. (1)In the spirit of finite volume methods, we shall consider the explicit schemeupsilon(mu)(n+1) = upsilon(mu)(n) - Deltat/hmu ((f) over bar (mu) - (f) over bar (mu)-) = [Dupsilon(n)](mu), (2)where mu is a point of an irregular grid Gamma, mu(-) is the left neighbor of A in Gamma, upsilon(mu)(n) approximate to 1/mu-mu(-) integral(mu-)(mu) u(x, t(n))dx are approximated cell averages of the solution, (f) over bar (mu) = (f) over bar (mu)(upsilon(n)) are the numerical fluxes, and D is the numerical evolution operator of the scheme.According to the definition of (f) over bar (mu), several schemes of this type have been proposed and successfully applied (LeVeque, 1990). Godunov, Lax-Wendroff, and ENO are some of the popular names. Godunov scheme resolves well the shocks, but accuracy (of first order) is poor in smooth regions. Lax-Wendroff is of second order, but produces dangerous oscillations close to shocks. ENO schemes are good alternatives, with high order and without serious oscillations. But the price is high computational cost.Ami Harten proposed in (Harten, 1994) a simple strategy to save expensive ENO flux calculations. The basic tools come from multiresolution analysis for cell averages on uniform grids, and the principle is that wavelet coefficients can be used for the characterization of local smoothness.. Typically, only few wavelet coefficients are significant. At the finest level, they indicate discontinuity points, where ENO numerical fluxes are computed exactly. Elsewhere, cheaper fluxes can be safely used, or just interpolated from coarser scales. Different applications of this principle have been explored by several authors, see for example (G-Muller and Muller, 1998).Our scheme also uses Ami Harten's ideas. But instead of evolving the cell averages on the finest uniform level, we propose to evolve the cell averages on sparse grids associated with the significant wavelet coefficients. This means that the total number of cells is small, with big cells in smooth regions and smaller ones close to irregularities. This task requires improved new tools, which are described next.

BOUND HYDROGENIC STATES IN A MAGNETIC FIELD: A FINITE-DIFFERENCE APPROACH

A finite-difference scheme is used to calculate bound electronic states of an electron in a hydrogen atom subject to a magnetic field. The numerical results are in good agreement with exact results, in the absence of the magnetic field, and with a two-parameters variational calculation, when the magnetic field is applied.

Cubic approximation for the transverse displacement in BEM for elastic plates analysis

This work presents an application for the plate analysis formulation by BEM where 3 boundary equations are used, written for the transverse displacement w and the normal and tangential derivatives partial derivativew/partial derivativen and partial derivativew/partial derivatives. In this extension, the transverse displacement w is approximated by a cubic polynomial and, as a consequence, partial derivativew/partial derivatives has a quadratic approximation. This alternative BEM formulation improves the analysis of thin plates, when compared to the formulation using the linear approximation for the displacements, mainly in the obtaining of the bending moments at the boundary of the plate. The implementation of this proposal to the computational codes is simple. (C) 2004 Published by Elsevier Ltd.

Stability analysis of power systems described with detailed models by automatic method

The power system stability analysis is approached taking into explicit account the dynamic performance of generators internal voltages and control devices. The proposed method is not a direct method in the usual sense since conclusion for stability or instability is not exclusively based on energy function considerations but it is automatic since the conclusion is achieved without an analyst intervention. The stability test accounts for the nonconservative nature of the system with control devices such as the automatic voltage regulator (AVR) and automatic generation control (AGC) in contrast with the well-known direct methods. An energy function is derived for the system with machines forth-order model, AVR and AGC and it is used to start the analysis procedure and to point out criticalities. The conclusive analysis itself is made by means of a method based on the definition of a region surrounding the equilibrium point where the system net torque is equilibrium restorative. This region is named positive synchronization region (PSR). Since the definition of the PSR boundaries have no dependence on modelling approximation, the PSR test conduces to reliable results. (C) 2008 Elsevier Ltd. All rights reserved.

Predictors of health status do not change over three-year periods and exacerbation makes difference in chronic obstructive pulmonary disease

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

Electroacupuncture analgesia in dogs: is there a difference between uni- and bi-lateral stimulation?

Objective To compare the analgesic effect of uni- and bi-lateral electroacupuncture (EA) in response to thermal and mechanical nociceptive stimuli and to investigate the cardiorespiratory, endocrine, and behavioral changes in dogs submitted to EA.Study design Prospective, randomized cross-over experimental study.Animals Eight adult, clinically healthy, cross-breed dogs, weighing 13 +/- 4 kg.Methods Dogs underwent electrostimulation at false acupoints (T-false); bilateral EA at acupoints, stomach 36, gall bladder 34 and spleen 6 (T-EA/bil); unilateral EA at the same points (T-EA/uni) or were untreated (T-control). All animals received acepromazine (0.05 mg kg(-1)) IV; and heart rate, pulse oximetry, indirect arterial blood pressure, respiratory rate, PECO2, rectal temperature, and plasma cortisol concentration were measured before, during, and after EA. Analgesia was tested using thoracic and abdominal cutaneous thermal and mechanical stimuli, and an interdigital thermal stimulus. Behavior was classified as calm or restless. Analysis of variance for repeated measures followed by Tukey's test was used for analysis of the data.Results There were no cardiorespiratory differences among the treatments. The cutaneous pain threshold was higher after EA, compared with false points. The latency period was shorter and analgesia was more intense in T-EA/bil than T-EA/uni, when both were compared with T-false and T-control. Six out of eight animals treated with EA were calm during treatment, and 5/8 and 4/8 of the T-false and T-control animals, respectively, were restless. Latency to interdigital thermal stimulation increased in T-EA/bil compared with the others. There was no difference in plasma cortisol concentrations among the treatments.Conclusions Bilateral EA produced a shorter latency period, a greater intensity, and longer duration of analgesia than unilateral stimulation, without stimulating a stress response.Clinical relevance Bilateral EA produces a better analgesic effect than unilateral EA.

s-wave approximation for asymmetry in nonmesonic decay of finite hypernuclei

We establish the bridge between the commonly used Nabetani-Ogaito-Sato-Kishimoto (NOSK) formula for the asymmetry parameter a(Lambda) in the Lambda p -> np emission of polarized hypernuclei, and the shell-model (SM) formalism for finite hypernuclei. We demonstrate that the s-wave approximation leads to a SM formula for a(Lambda) that is as simple as the NOSK one and that reproduces the exact results for (5)(Lambda)He and (12)(Lambda)C better than initially expected. The simplicity achieved here is indeed remarkable. The new formalism makes the theoretical evaluation of a(Lambda) more transparent and explains clearly why the one-meson exchange model is unable to account for the experimental data of (5)(Lambda)He.

Low-energy direct muon transfer from H to Ne10+, S16+ and Ar18+ using the two-state close-coupling approximation to the Faddeev-Hahn-type equation

We perform a three-body calculation of direct muon-transfer rates from thermalized muonic hydrogen isotopes to bare nuclei Ne10+, S16+ and Ar18+ employing integro-differential Faddeev-Hahn-type equations in configuration space with a two-state close-coupling approximation scheme. All Coulomb potentials including the strong final-state Coulomb repulsion are treated exactly. A long-range polarization potential is included in the elastic channel to take into account the high polarizability of the muonic hydrogen. The transfer rates so-calculated are in good agreement with recent experiments. We find that the muon is captured predominantly in the n = 6, 9 and 10 states of muonic Ne10+, S16+ and Ar18+, respectively.

Coulomb gauge quark model, vacuum condensates and high density quark matter

We consider here a Coulomb gauge quark model which includes an explicit construct for a nontrivial vacuum structure in QCD at finite density. Non-perturbative renormalization of ultraviolet diverges is performed by adding counterterms. The equation of state for u and d quark matter at zero temperature is calculated in the Hartree-Fock approximation.

Bose-Einstein condensation dynamics in three dimensions by the pseudospectral and finite-difference methods

We suggest a pseudospectral method for solving the three-dimensional time-dependent Gross-Pitaevskii (GP) equation, and use it to study the resonance dynamics of a trapped Bose-Einstein condensate induced by a periodic variation in the atomic scattering length. When the frequency of oscillation of the scattering length is an even multiple of one of the trapping frequencies along the x, y or z direction, the corresponding size of the condensate executes resonant oscillation. Using the concept of the differentiation matrix, the partial-differential GP equation is reduced to a set of coupled ordinary differential equations, which is solved by a fourth-order adaptive step-size control Runge-Kutta method. The pseudospectral method is contrasted with the finite-difference method for the same problem, where the time evolution is performed by the Crank-Nicholson algorithm. The latter method is illustrated to be more suitable for a three-dimensional standing-wave optical-lattice trapping potential.

Combined D0 measurements constraining the CP-violating phase and width difference in the B-s(0) system

We combine the D0 measurement of the width difference between the light and heavy B-s(0) mass eigenstates and of the CP-violating mixing phase determined from the time-dependent angular distributions in the B-s(0)-> J/psi phi decays along with the charge asymmetry in semileptonic decays also measured with the D0 detector. With the additional constraint from the world average of the flavor-specific B-s(0) lifetime, we obtain Delta Gamma(s)equivalent to(Gamma(L)-Gamma(H))=0.13 +/- 0.09 ps(-1) and vertical bar phi(s)vertical bar=0.70(-0.47)(+0.39) or Delta Gamma(s)=-0.13 +/- 0.09 ps(-1) and vertical bar phi(s)vertical bar=2.44(-0.39)(+0.47). The data sample corresponds to an integrated luminosity of 1.1 fb(-1) accumulated with the D0 detector at the Fermilab Tevatron Collider.