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.

The peremptory influence of a uniform background for trapping neutral fermions with an inversely linear potential

The problem of neutral fermions subject to an inversely linear potential is revisited. It is shown that an infinite set of bound-state solutions can be found on the condition that the fermion is embedded in an additional uniform background potential. An apparent paradox concerning the uncertainty principle is solved by introducing the concept of effective Compton wavelength.

New laws for author's rights: adequate for digital preservation?

Introduction. Besides technical and economic-organisational problems, digital preservation also faces legal issues, especially regarding copyright legislation, since all digital preservation strategies involve copying materials and/or using software which is typically copyrighted. The purpose of this paper is to ascertain the extent to which current copyright laws meet the preservation requirements of library materials..Methodology. A cross-sectional analysis of recently updated national copyright laws as well as the impact of the other two protection methods: contractual and technological.Results. Even after the latest updates current copyright legislation is almost useless for digital preservation activities since the opportunities provided by WIPO Copyright Treaty to adapt and extend copyright exceptions and limitations have been used to the full.Conclusion. We need a legislation reform that will make it possible to carry out all required copying and communication activities and software use, even if circumventing technological protection is needed. But that is not enough for licensed works. The best solution for this kind of work is to include specific clauses in the licences that facilitate preservation activities. Thus, cooperation between both parties, libraries and rights' holders, is essential.

Library exceptions in the copyright laws of Ibero-American countries

The development of the digital setting has made it absolutely necessary to revise copyright legislation as a whole, including the exceptions that benefit libraries and similar institutions. Unfortunately, adaptation to the new technological reality is not taking place satisfactorily in most countries - the predominating trend is a refortification of copyright as opposed to user interests, and the maintenance of a certain pre-digital philosophy. In the case of Ibero-America the problem is twofold: aside from obsolescence or a lack of adaptation to the new technological setting, there are countries that have not yet included library-related exceptions in the national laws. Moreover, these happen to be developing countries, whose needs and interests do not coincide with those of the richer nations who paved the path to be followed by international treaties and copy- or copyright agreements. This study looks into the situation of exceptions to copyright to benefit libraries in the countries constituting Ibero-America, with a comparative analysis of the most significant characteristics of their national laws. It is concluded that it is crucial for these countries to take advantage of the options offered through the WIPO Copyright Treaty of 1996 and the results of the WIPO Development Agenda to update their legislation, in order that copyright will be respected, while at the same time making it easier for libraries to continue carrying out their social function in an adequate manner, always with the understanding of the developmental context of these countries.

Some experiments with high order compact methods using a computer algebra software - Part II (non-uniform grid)

We generalize a procedure proposed by Mancera and Hunt [P.F.A. Mancera, R. Hunt, Some experiments with high order compact methods using a computer algebra software-Part 1, Appl. Math. Comput., in press, doi: 10.1016/j.amc.2005.05.015] for obtaining a compact fourth-order method to the steady 2D Navier-Stokes equations in the streamfunction formulation-vorticity using the computer algebra system Maple, which includes conformal mappings and non-uniform grids. To analyse the procedure we have solved a constricted stepped channel problem, where a fine grid is placed near the re-entrant corner by transformation of the independent variables. (c) 2006 Elsevier B.V. All rights reserved.

Phosphorylation of Histone H3S10 in Animal Chromosomes: Is There a Uniform Pattern?

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

Bicomplexes and conservation laws in non-Abelian Toda models

A bicomplex structure is associated with the Leznov-Saveliev equation of integrable models. The linear problem associated with the zero-curvature condition is derived in terms of the bicomplex linear equation. The explicit example of a non-Abelian conformal affine Toda model is discussed in detail and its conservation laws are derived from the zero-curvature representation of its equation of motion.

Evolution of spherical non-uniform perturbations in the universe

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

Soliton dynamics at an interface between a uniform medium and a nonlinear optical lattice

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

Non-uniform drag force on the Fermi accelerator model

Some dynamical properties of a particle suffering the action of a generic drag force are obtained for a dissipative Fermi Acceleration model. The dissipation is introduced via a viscous drag force, like a gas, and is assumed to be proportional to a power of the velocity: F alpha -nu(gamma). The dynamics is described by a two-dimensional nonlinear area-contracting mapping obtained via the solution of Newton's second law of motion. We prove analytically that the decay of high energy is given by a continued fraction which recovers the following expressions: (i) linear for gamma = 1; (ii) exponential for gamma = 2; and (iii) second-degree polynomial type for gamma = 1.5. Our results are discussed for both the complete version and the simplified version. The procedure used in the present paper can be extended to many different kinds of system, including a class of billiards problems.

Uniform exponential dichotomy and continuity of attractors for singularly perturbed damped wave equations

For eta >= 0, we consider a family of damped wave equations u(u) + eta Lambda 1/2u(t) + au(t) + Lambda u = f(u), t > 0, x is an element of Omega subset of R-N, where -Lambda denotes the Laplacian with zero Dirichlet boundary condition in L-2(Omega). For a dissipative nonlinearity f satisfying a suitable growth restrictions these equations define on the phase space H-0(1)(Omega) x L-2(Omega) semigroups {T-eta(t) : t >= 0} which have global attractors A(eta) eta >= 0. We show that the family {A(eta)}(eta >= 0), behaves upper and lower semi-continuously as the parameter eta tends to 0(+).

Geometrically uniform hyperbolic codes

In this paper we generalize the concept of geometrically uniform codes, formerly employed in Euclidean spaces, to hyperbolic spaces. We also show a characterization of generalized coset codes through the concept of G-linear codes.

Uniform estimates of attracting sets of extended Lurie systems using LMIs

This research presents a systematic procedure to obtain estimates, via extended Lyapunov functions, of attracting sets of a class of nonlinear systems, as well as an estimate of their stability regions. The considered class of nonlinear systems, called in this note the extended Lurie system, consists of nonlinear systems like those of the Lurie problem where one of the nonlinear functions can violate the sector conditions of the Lurie problem around the origin. In case of nonautonomous systems the concept of absolute stability is extended and uniform estimates of the attracting set are obtained. Two classical nonlinear systems, the forced duffing equation and the Van der Pol system, are analyzed with the proposed procedure.

Embedded crack elements with non-uniform discontinuity modes

The consequences of the use of embedded crack finite elements with uniform discontinuity modes (opening and sliding) to simulate crack propagation in concrete are investigated. It is shown the circumstances in which the consideration of uniform discontinuity modes is not suitable to accurately model the kinematics induced by the crack and must be avoided. It is also proposed a technique to embed cracks with non-uniform discontinuity modes into standard displacement-based finite elements to overcome the shortcomings of the uniform discontinuity modes approach.