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.

Evaluation of a bounded high order upwind scheme for 3D incompressible free surface flow computations

fit the context of normalized variable formulation (NVF) of Leonard and total variation diminishing (TVD) constraints of Harten. this paper presents an extension of it previous work by the authors for solving unsteady incompressible flow problems. The main contributions of the paper are threefold. First, it presents the results of the development and implementation of a bounded high order upwind adaptative QUICKEST scheme in the 3D robust code (Freeflow), for the numerical solution of the full incompressible Navier-Stokes equations. Second, it reports numerical simulation results for 1D hock tube problem, 2D impinging jet and 2D/3D broken clam flows. Furthermore, these results are compared with existing analytical and experimental data. and third, it presents the application of the numerical method for solving 3D free surface flow problems. (C) 2007 IMACS. Published by Elsevier B.V. All rights reserved,

Computations of the flow past a still sphere at moderate reynolds numbers using an immersed boundary method

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

3D relativistic hydrodynamic computations using lattice-QCD-inspired equations of state

In this communication, we report results of three-dimensional hydrodynamic computations, by using equations of state with a critical end Point as suggested by the lattice QCD. Some of the results are an increase of the multiplicity in the mid-rapidity region and a larger elliptic-flow parameter nu(2). We discuss also the effcts of the initial-condition fluctuations and the continuous emission.

Some superstring amplitude computations with the non-minimal pure spinor formalism

We use the non-minimal pure spinor formalism to compute in a super-Poincare covariant manner the four-point massless one and two-loop open superstring amplitudes, and the gauge anomaly of the six-point one-loop amplitude. All of these amplitudes are expressed as integrals of ten-dimensional superfields in a pure spinor superspace which involves five theta coordinates covariantly contracted with three pure spinors. The bosonic contribution to these amplitudes agrees with the standard results, and we demonstrate identities which show how the t(8) and epsilon(10) tensors naturally emerge from integrals over pure spinor superspace.

Wavelet-galerkin method for computations of electromagnetic fields-computation of connection coefficients

One of the key issues which makes the waveletGalerkin method unsuitable for solving general electromagnetic problems is a lack of exact representations of the connection coefficients. This paper presents the mathematical formulae and computer procedures for computing some common connection coefficients. The characteristic of the present formulae and procedures is that the arbitrary point values of the connection coefficients, rather than the dyadic point values, can be determined. A numerical example is also given to demonstrate the feasibility of using the wavelet-Galerkin method to solve engineering field problems. © 2000 IEEE.

Regional differences in the spatial and temporal heterogeneity of oceanographic habitat used by Steller sea lions

Over the past three decades, the decline and altered spatial distribution of the western stock of Steller sea lions (Eumetopias jubatus) in Alaska have been attributed to changes in the distribution or abundance of their prey due to the cumulative effects of fisheries and environmental perturbations. During this period, dietary prey occurrence and diet diversity were related to population decline within metapopulation regions of the western stock of Steller sea lions, suggesting that environmental conditions may be variable among regions. The objective of this study, therefore, was to examine regional differences in the spatial and temporal heterogeneity of oceanographic habitat used by Steller sea lions within the context of recent measures of diet diversity and population trajectories. Habitat use was assessed by deploying satellite-depth recorders and satellite relay data loggers on juvenile Steller sea lions (n = 45) over a five-year period (2000–2004) within four regions of the western stock, including the western, central, and eastern Aleutian Islands, and central Gulf of Alaska. Areas used by sea lions during summer months (June, July, and August) were demarcated using satellite telemetry data and characterized by environmental variables (sea surface temperature [SST] and chlorophyll a [chl a]), which possibly serve as proxies for environmental processes or prey. Spatial patterns of SST diversity and Steller sea lion population trends among regions were fairly consistent with trends reported for diet studies, possibly indicating a link between environmental diversity, prey diversity, and distribution or abundance of Steller sea lions. Overall, maximum spatial heterogeneity coupled with minimal temporal variability of SST appeared to be beneficial for Steller sea lions. In contrast, these patterns were not consistent for chl a, and there appeared to be an ecological threshold. Understanding how Steller sea lions respond to measures of environmental heterogeneity will ultimately be useful for implementing ecosystem management approaches and developing additional conservation strategies.

A DESCRIPTIVE ANALYSIS OF THE SEASONAL VARIATION OF PHYSICAL OCEANOGRAPHIC CHARACTERISTICS IN THE NORTHERN REGION OF THE TODOS OS SANTOS BAY (BAHIA, BRAZIL)

The results obtained in the August and December 2003, August 2004 and January 2005 oceanographic campaigns in the northern region of the Todos os Santos Bay (lat. 12 degrees 44.5`S; long. 038 degrees 35.00`W) between the Madre de Deus and Mare islands are analyzed. Instruments of continuous and discrete samplings were used to measure hydrographic properties currents and tides. The water mass of the northern region of the bay is forced by semidiurnal and mesotides of form number 0.08 and the lunar component M(2) height was estimated at 91cm. The time series of the surface currents indicated movements in the N/S direction, forced by the tide with maximum magnitudes of 0.73 m.s(-1) on the December 2003 campaign. However, in August 2004 the currents were dominated by the wind stress forcing, with a maximum speed of 1.85 m.s(-1) and SE direction. Near the bottom, the influence of the tide is not as evident, with a decrease in intensity due to internal and bottom friction, with a maximum velocity of 0.17 m.s(-1). The thermal and haline structures were weakly horizontally, as well as vertically stratified, with extreme values varying in the intervals 23 degrees C (August, 2004) to 28 degrees C (December, 2003) and 31.0 psu (August, 2003) to 36.0 psu (December, 2003), respectively. Some conclusions may be drawn from these results: i) The signs of the dilution of the fresh water discharges of the Caipe, Mataripe and Sao Paulo rivers in the region under the influence of the RLAM were observed only during the winter periods, but in the summer the region was flooded by waters of oceanic origin and the salinities above 36.0 indicated TW mass intrusion; ii) The N-S circulation near the RLAM is strongly dominated by the tide, and the importance of the M(2) component was unequivocal, however, the E-W component presented some tidal modulation away from abrupt bottom topographical changes, and iii) The residual series, calculated as the difference between the original and modeled, is about 1/4 of the original and confirmed its semidiurnal character.

Interannual variation in commercial oyster (Crassostrea gigas) farming in the sea (Florianopolis, Brazil, 27 degrees 44 ' S; 48 degrees 33 ' W) in relation to temperature, chlorophyll a and associated oceanographic conditions

Aquaculture of filter-feeding bivalve mollusks involves the fruitful conversion of marine particulate organic matter into premium protein of high nutritive value. Culture performance of bivalves is largely dependent on hydrological conditions and directly affected by e. g. temperature and chlorophyll levels. Accordingly, these parameters may be related with seasonality but also with oceanographic features combined with climate events. Yields of Pacific cupped oyster (Crassostrea gigas) reared at commercial procedures in suspended structures (long-lines) in a sheltered bay in Southern Brazil (Santa Catarina State, 27S 43'; 48 W 30') were evaluated in relation to local environmental conditions: sea surface temperature, chlorophyll a concentration, and associate effects of cold fronts events and El Nino and La Nina periods. Outputs from four consecutive commercial crop years were analyzed (2005/06, 2006/07, 2007/08, 2008/09) in terms of oyster survival and development time during the following grow-out phases of the culture cycle: seed to juvenile, juvenile to adult, adult to marketable. Since culture management and genetics were standardized significant differences verified among crop performance could be mostly related to environmental effects. Time series of temperature and chlorophyll a (remote sensing data) from crop periods displayed significant seasonal and interannual variation. As expected, performance during initial grow-out stages (seed to juvenile) was critical for final crop yield. Temperature was the main factor affecting survival in these initial stages with a trend of negative correlation, though not statistically significant. On the other hand, oyster development rate was significantly and positively affected by chlorophyll a concentration. Chlorophyll a values could be increased by upwelled cold nutrient-rich South Atlantic Central Water (SACW, related to predominant Northern winds) though further dependent on occurrence of Southern winds (cold fronts) to assist seawater penetration into the sheltered farming area. Lower salinity nutrient-rich northward drifted waters from La Plata River discharge may also result in chlorophyll a rise in the farming area. The El Nino period (July 2006 to February 2007) coincided with lower chlorophyll a levels in the farming site that may be related to both decreased number of cold fronts as well as predominance of Northern winds that retain northward spreading of La Plata River discharge waters. In contrast, the La Nina period (August 2007 to June 2008) corresponded to higher chlorophyll a values in the farming area by both upwelling of SACW and penetration of La Plata River discharge water assisted by increased occurrence of Southern winds and cold fronts. The recognition of the potentially changing climate and effects upon the environment will be an important step in planning future development of bivalve aquaculture.

Identification, Abundance and Rearing of Mysids of Gran Canaria: aplication to laboratory and oceanographic respiratory metabolism studies

Programa de doctorado en Oceanografía