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.

Influência da densidade de estocagem sobre o crescimento, ciclo de muda e o comportamento em juvenis do camarão marinho Litopenaeus vannamei

One of the factors that may interfere with the cultivation of Litopenaeus vannamei is the population density. This study aimed to assess the effect of density on growth, mortality, physical integrity and behavior of shrimp. The study was divided into two stages. At first, the shrimp were placed in tanks at densities of 50, 75 and 100 shrimp m-2. The animals were monitored in relation to the degree of proventricular filling, the stage of the molt cycle and physical integrity three times a week and in relation to the weight and length once a week. Mortality, growth and proventricular filling were not influenced by the density; frequency of records in specific stages of the molt cycle varied according to the density. The lower proportion of broken appendages and higher frequency of necrotic lesions occurred in lower density. The second stage of the research, conducted in aquaria, was divided into two parts. The first described social or feeding behavioral categories: slow displacement by contact, slow displacement by approximation, abrupt displacement by contact, abrupt displacement by approximation, reactivity, cannibalism, occupying the tray, get feed in the tray and get feed outside the tray. In the second part, these and other behavioral categories, described in the literature, were recorded in densities of 50, 75 and 100 shrimp m-2. Mortality was more frequent in higher density. The frequency of most behaviors mentioned above was very low, not differing between densities or being too low to determine differences between them. The behavioral profile of animals in different densities was, in general, very similar, with no difference in exploration, digging and cleaning between the densities. Even so, inactivity, feeding, crawling, burrowing, swimming, and proximity between animals were influenced by the density. These results suggest that some behaviors suffer greater interference from population density. However, the density may not have a broader influence on the animals when other factors, such as physico-chemical parameters of water and feed offering, are adequate

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.

Finite-volume form factors in semi-classical approximation

A semi-classical approach is used to obtain Lorentz covariant expressions for the form factors between the kink states of a quantum field theory with degenerate vacua. Implemented on a cylinder geometry it provides an estimate of the spectral representation of correlation functions in a finite volume. Illustrative examples of the applicability of the method are provided by the sine-Gordon and the broken phi(4) theories in 1 + 1 dimensions. (C) 2003 Elsevier B.V. All rights reserved.

Quark condensates and high density quark matter

We consider a Coulomb gauge quark model which includes an explicit construct for a nontrivial vacuum structure in QCD. The dynamics is described by a Hamiltonain that contains a linearly rising confining potential and longitudinal and transverse Coulomb-type interactions. The Coulomb potential gives rise to ultraviolate divergences which are non-perturbatively renormalized by adding appropriate counter terms to the Hamiltonian. The equation of state for u and d quark matter at zero temperature is derived in the Hartree-Fock approximation.

A dynamical gluon mass solution in Mandelstam's approximation

We discuss the pure gauge Schwinger-Dyson equation for the gluon propagator in the Landau gauge within an approximation proposed by Mandelstam many years ago. We show that a dynamical gluon mass arises as a solution. This solution is obtained numerically in the full range of momenta that we have considered without the introduction of any ansatz or asymptotic expression in the infrared region. The vertex function that we use follows a prescription formulated by Cornwall to determine the existence of a dynamical gluon mass in the light cone gauge. The renormalization procedure differs from the one proposed by Mandelstam and allows for the possibility of a dynamical gluon mass. Some of the properties of this solution, such as its dependence on A(QCD) and its perturbative scaling behavior are also discussed.