Wavelets And Adaptive Grids For The Discontinuous Galerkin Method

In this paper, space adaptivity is introduced to control the error in the numerical solution of hyperbolic systems of conservation laws. The reference numerical scheme is a new version of the discontinuous Galerkin method, which uses an implicit diffusive term in the direction of the streamlines, for stability purposes. The decision whether to refine or to unrefine the grid in a certain location is taken according to the magnitude of wavelet coefficients, which are indicators of local smoothness of the numerical solution. Numerical solutions of the nonlinear Euler equations illustrate the efficiency of the method. © Springer 2005.

A systematic review of diapoma (teleostei: characiformes: characidae: stevardiinae: diapomini) with descriptions of two new species from southern Brazil

Diapoma is reviewed and four species are recognized: (1) Diapoma thauma, new species, from streams of the rio Jacuí basin, state of Rio Grande do Sul; (2) D. pyrrhopteryx, new species collected from the rio Canoas and streams flowing into this basin in the states of Rio Grande do Sul and Santa Catarina, Brazil; (3) Diapoma terofali, from streams flowing into rio Uruguay in Uruguay and Rio Grande do Sul, Brazil and streams flowing into rio de la Plata, Argentina; and (4) Diapoma speculiferum, from lowland coastal streams in Rio Grande do Sul, Brazil and Uruguay. Diapoma pyrrhopteryx possess the posteroventral opercular elongation typical of D. speculiferum, type species of the genus, but which is absent in D. thauma and D. terofali. Nonetheless, all the diapomin species have the caudal pouch organ about equally developed in both sexes and the dorsal portion of the pouch opening bordered by a series of 3 to 8 elongated scales, the two derived features that characterize the group. The two previously described species, D. speculiferum and D. terofali, are redescribed. Previous hypotheses of relationships among the diapomin genera Planaltina, Diapoma and Acrobrycon are discussed on the basis of preliminary morphological information. It is proposed that the Diapomini is a monophyletic group. An identification key, information on sexual dimorphism, gonad anatomy, reproductive mode and distribution of the species of Diapoma are provided.

OFFSPRING IN RESPONSE TO PARENTAL FEMALE DENSITIES IN THE FRUIT FLY PARASITOID DIACHASMIMORPHA LONGICAUDATA (HYMENOPTERA: BRACONIDAE: OPIINAE)

As part of an evaluation of the braconid parasitoid Diachasmimorpha longicaudata (Ashmead) as a biocontrol agent of Ceratitis capitata (Wiedemann) in Brazil, the aims in the current study were to find the best parental ratio of females to males in the rearing cages in order to get the highest female biased offspring in the parasitoid rearing process, and to verify the parasitism efficiency on C. capitata according to parental female densities. Three treatments were assessed: T1 (20 females: 20 males), T2 (60 females: 20 males) and T3 (100 females: 20 males). Ten late-third instars of C. capitata were offered daily to each female parasitoid from the 1st to the 12th d of age. The parental female productivity, fecundity, offspring sex ratio, percentage of parasitoid emergence, and daily mortality of parental females and males at different female/male densities were evaluated. The results indicated that numbers higher than 20 parental females did not affect offspring sex ratio, overall offspring production, nor the percent parasitism. Female biased offspring occurred in all three parental female/male ratios analyzed in this study, except that predominately males developed from parasitoid eggs laid in the age interval 1-2 d post emergence. Higher parasitoid female productivity and fecundity were found at the 1:1 female/male per cage density whereas lower productivity and fecundity were recorded at the 5:1 female/male ratio. Higher female/male ratio in the parental cages increased the mortality rate of females but did not influence the number of parental male deaths. The results may facilitate advancement of an optimum mass-rearing system to aid in control of C. capitata in Brazil.

Cosmic distance duality relation and the shape of galaxy clusters

Context. Observations in the cosmological domain are heavily dependent on the validity of the cosmic distance-duality (DD) relation, eta = D(L)(z)(1+ z)(2)/D(A)(z) = 1, an exact result required by the Etherington reciprocity theorem where D(L)(z) and D(A)(z) are, respectively, the luminosity and angular diameter distances. In the limit of very small redshifts D(A)(z) = D(L)(z) and this ratio is trivially satisfied. Measurements of Sunyaev-Zeldovich effect (SZE) and X-rays combined with the DD relation have been used to determine D(A)(z) from galaxy clusters. This combination offers the possibility of testing the validity of the DD relation, as well as determining which physical processes occur in galaxy clusters via their shapes. Aims. We use WMAP (7 years) results by fixing the conventional Lambda CDM model to verify the consistence between the validity of DD relation and different assumptions about galaxy cluster geometries usually adopted in the literature. Methods. We assume that. is a function of the redshift parametrized by two different relations: eta(z) = 1+eta(0)z, and eta(z) = 1+eta(0)z/(1+z), where eta(0) is a constant parameter quantifying the possible departure from the strict validity of the DD relation. In order to determine the probability density function (PDF) of eta(0), we consider the angular diameter distances from galaxy clusters recently studied by two different groups by assuming elliptical (isothermal) and spherical (non-isothermal) beta models. The strict validity of the DD relation will occur only if the maximum value of eta(0) PDF is centered on eta(0) = 0. Results. It was found that the elliptical beta model is in good agreement with the data, showing no violation of the DD relation (PDF peaked close to eta(0) = 0 at 1 sigma), while the spherical (non-isothermal) one is only marginally compatible at 3 sigma. Conclusions. The present results derived by combining the SZE and X-ray surface brightness data from galaxy clusters with the latest WMAP results (7-years) favors the elliptical geometry for galaxy clusters. It is remarkable that a local property like the geometry of galaxy clusters might be constrained by a global argument provided by the cosmic DD relation.

TESTING THE DISTANCE-DUALITY RELATION WITH GALAXY CLUSTERS AND TYPE Ia SUPERNOVAE

In this Letter, we propose a new and model-independent cosmological test for the distance-duality (DD) relation, eta = D(L)(z)(1 + z)(-2)/D(A)(z) = 1, where D(L) and D(A) are, respectively, the luminosity and angular diameter distances. For D(L) we consider two sub-samples of Type Ia supernovae (SNe Ia) taken from Constitution data whereas D(A) distances are provided by two samples of galaxy clusters compiled by De Filippis et al. and Bonamente et al. by combining Sunyaev-Zeldovich effect and X-ray surface brightness. The SNe Ia redshifts of each sub-sample were carefully chosen to coincide with the ones of the associated galaxy cluster sample (Delta z < 0.005), thereby allowing a direct test of the DD relation. Since for very low redshifts, D(A)(z) approximate to D(L)(z), we have tested the DD relation by assuming that. is a function of the redshift parameterized by two different expressions: eta(z) = 1 + eta(0)z and eta(z) = 1 +eta(0)z/(1 + z), where eta(0) is a constant parameter quantifying a possible departure from the strict validity of the reciprocity relation (eta(0) = 0). In the best scenario (linear parameterization), we obtain eta(0) = -0.28(-0.44)(+0.44) (2 sigma, statistical + systematic errors) for the De Filippis et al. sample (elliptical geometry), a result only marginally compatible with the DD relation. However, for the Bonamente et al. sample (spherical geometry) the constraint is eta(0) = -0.42(-0.34)(+0.34) (3 sigma, statistical + systematic errors), which is clearly incompatible with the duality-distance relation.

A novel aerobic-anoxic biological filter for nitrogen removal from UASB effluent using biogas compounds as electron donors for denitrification

The performance of a new trickling filter (TF) configuration composed of an upper compartment for nitrification and a lower compartment for denitrification of effluent from a UASB reactor treating domestic sewage was evaluated. The TF was packed with new plastic material characterized by its durability and high percentage of void spaces. The feasibility of using the reduced compounds present in the biogas produced by a UASB reactor as electron donor for denitrification was also evaluated. Efficient nitrification and denitrification was achieved for the mean hydraulic (5.6 m(3) m(-2) d(-1)) organic (0.26 kg COD m(-3) d(-1)) and ammonia-N (0.08 kg m(-3) d(-1)) loading rates applied, resulting in ammonia-N removal ranging from 60 to 74%. The final effluent presented ammonia-N lower than 13 mg L(-1). Despite the presence of dissolved oxygen (DO) in the denitrification compartment, its performance was considered quite satisfactory and final nitrate concentrations were lower than 10 mg L(-1). The results indicate that methane was the main electron donor used for denitrification. Additionally, denitrification can probably be improved by avoiding high DO concentration in the denitrification compartment and by enhancing biogas transfer in the anoxic zone.

Friction loss correction in flowing well discharge tests

Artesian confined aquifers do not need pumping energy, and water from the aquifer flows naturally at the wellhead. This study proposes correcting the method for analyzing flowing well tests presented by Jacob and Lohman (1952) by considering the head losses due to friction in the well casing. The application of the proposed correction allowed the determination of a transmissivity (T = 411 m(2)/d) and storage coefficient (S = 3 x 10(-4)) which appear to be representative for the confined Guarani Aquifer in the study area. Ignoring the correction due to head losses in the well casing, the error in transmissivity evaluation is about 18%. For the storage coefficient the error is of 5 orders of magnitude, resulting in physically unacceptable value. The effect of the proposed correction on the calculated radius of the cone of depression and corresponding well interference is also discussed.

Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method.

Flow and axial dispersion simulation for traveling axisymmetric Taylor vortices

Numerical experiments using a finite difference method were carried out to determine the motion of axisymmetric Taylor vortices for narrow-gap Taylor vortex flow. When a pressure gradient is imposed on the flow the vortices are observed to move with an axial speed of 1.16 +/- 0.005 times the mean axial flow velocity. The method of Brenner was used to calculate the long-time axial spread of material in the flow. For flows where there is no pressure gradient, the axial dispersion scales with the square root of the molecular diffusion, in agreement with the results of Rosen-bluth et al. for high Peclet number dispersion in spatially periodic flows with a roll structure. When a pressure gradient is imposed the dispersion increases by an amount approximately equal to 6.5 x 10(-4) (W) over bar(2)d(2)/D-m, where (W) over bar is the average axial velocity in the annulus, analogous to Taylor dispersion for laminar flow in an empty tube.

Two dimensional computational fluid dynamic models for waste stabilisation ponds

Traditional waste stabilisation pond (WSP) models encounter problems predicting pond performance because they cannot account for the influence of pond features, such as inlet structure or pond geometry, on fluid hydrodynamics. In this study, two dimensional (2-D) computational fluid dynamics (CFD) models were compared to experimental residence time distributions (RTD) from literature. In one of the-three geometries simulated, the 2-D CFD model successfully predicted the experimental RTD. However, flow patterns in the other two geometries were not well described due to the difficulty of representing the three dimensional (3-D) experimental inlet in the 2-D CFD model, and the sensitivity of the model results to the assumptions used to characterise the inlet. Neither a velocity similarity nor geometric similarity approach to inlet representation in 2-D gave results correlating with experimental data. However. it was shown that 2-D CFD models were not affected by changes in values of model parameters which are difficult to predict, particularly the turbulent inlet conditions. This work suggests that 2-D CFD models cannot be used a priori to give an adequate description of the hydrodynamic patterns in WSP. (C) 1998 Elsevier Science Ltd. All rights reserved.

Ester prodrugs of a potent analgesic, morphine-6-sulfate: syntheses, spectroscopic and physicochemical properties

The aim of this work is to develop 3-acyl prodrugs of the potent analgesic morphine-6-sulfate (M6S). These are expected to have higher potency and/or exhibit longer duration of analgesic action than the parent compound. M6S and the prodrugs were synthesized, then purified either by recrystallization or by semi-preparative HPLC and the structures confirmed by mass spectrometry, IR spectrophotometry and by detailed 1- and 2-D NMR studies. The lipophilicities of the compounds were assessed by a combination of shake-flask, group contribution and HPLC retention methods. The octanol-buffer partition coefficient could only be obtained directly for 3-heptanoylmorphine-6-sulfate, using the shake-flask method. The partition coefficients (P) for the remaining prodrugs were estimated from known methylene group contributions. A good linear relationship between log P and the HPLC log capacity factors was demonstrated. Hydrolysis of the 3-acetyl prodrug, as a representative of the group, was found to occur relatively slowly in buffers (pH range 6.15-8.01), with a small buffer catalysis contribution. The rates of enzymatic hydrolysis of the 3-acyl group in 10% rat blood and in 10% rat brain homogenate were investigated. The prodrugs followed apparent first order hydrolysis kinetics, with a significantly faster hydrolysis rate found in 10% rat brain homogenate than in 10% rat blood for all compounds. (C) 1998 Elsevier Science B.V. All rights reserved.

Brain region-specific studies of the excitatory behavioral effects of morphine-3-glucuronide

This study was designed to determine in rats whether morphine-3-glucuronide (M3G) produces its neuro-excitatory effects most potently in the ventral hippocampus (as has been reported previously for subanalgesic doses of opioid peptides). Guide cannulae were implanted into one of seven regions of the rat brain: lateral ventricle; ventral, CA1 and CA2-CA3 regions of the hippocampus; amygdala; striatum or cortex. After a 7 day recovery period, rats received intracerebral injections of (i) M3G (1.1 or 11 nmol) (ii) DADLE ([D-Ala(2),D-Leu(5)]enkephalin), (45 nmol, positive controls) or (iii) vehicle (deionised water), and behavioral excitation was quantified over 80 min. High-dose M3G (11 nmol) evoked behavioral excitation in all brain regions but the onset, severity and duration of these effects varied considerably among brain regions. By contrast, low-dose M3G (1.1 nmol) evoked excitatory behaviors only when administered into the ventral hippocampus and the amygdala, with the most potent effects being observed in the ventral hippocampus. Prior administration of the nonselective opioid antagonists, naloxone and beta-funaltrexamine into the ventral hippocampus, markedly attenuated low-dose M3G's excitatory effects but did not significantly alter levels of excitation evoked by high-dose M3G. Naloxone given 10 min after M3G (1.1 or 11 nmol) did not significantly attenuate behavioral excitation. Thus, M3G's excitatory behavioral effects occur most potently in the ventral hippocampus as reported previously for subanalgesic doses of opioid peptides, and appear to be mediated through at least two mechanisms, one possibly involving excitatory opioid receptors and the other, non-opioid receptors.

Induction of an antigen-specific CTL response by a conformationally biased agonist of human C5a anaphylatoxin as a molecular adjuvant

A conformationally biased decapeptide agonist of human C5a anaphylatoxin (YSPKPMPLaR) was used as a molecular adjuvant in stimulating an Ag-specific CTL response against murine P815S target cells expressing an Ld-restricted CTL epitope of the hepatitis B surface Ag (HBsAg), Groups of BALB/c mice (H-2(d)) were immunized with aqueous solutions of the HBsAg CTL epitopes (IPQSLDSWWTSL and IPQSLDSTaVTSLRR); the C5a agonist (YSFKPMPLaR); the C5a agonist and HBsAg CTL epitopes admired (IPQSLDSWWTSL and IPQSLDSWWTSLRR + YSFKPMPLaR); the C5a-active, HBsAg CTL epitope-C5a agonist constructs (IPQSLDSWWTSLYSFKPMPLaR, IPQSLDSWWTSLRRYSFKPMPLaR, and IPQSLDSWWTSLRVRRYSFPMPLaR); a C5a-inactive, reverse-moiety construct (YSFKPMPLaRRRIPQSLDSWWTSL); and a C5a-attenuated, carboxyl-terminal-blocked construct (IPQSLDSWWTSLRRYSFKPMPLaRG). Ag-specific CD8(+) CTL responses were observed after the secondary boost in the absence of any added adjuvant only in mice that were immunized with C5a-active contructs, IPQSLDSWWTSLRRYSFKPMPLaR and IPQSLDSWWTSLRVRRYSFKPMPLaR. These two C5a-active immunogens contained potential subtilisin-sensitive linker sequences between the HBsAg CTL epitope and the C5a agonist; i.e., a double-Arg (RR) and a furin protease sensitive sequence (RVRR), The introduction of these potentially cleavable sequences may be a method of increasing the likelihood of liberating the CTL epitope from the C5a agonist by intracellular proteases, thereby facilitating entry of the epitope into Ag-processing pathways via an exogenous route.

Papillomavirus virus-like particles for the delivery of multiple cytotoxic T cell epitopes

Chimeric papillomavirus (PV) virus-like particles (VLPs) based on the bovine papillomavirus type 1 (BPV-1) L1 protein were constructed by replacing the 23-carboxyl-terminal amino acids of the BPV1 major protein L1 with an artificial polytope minigene, containing known CTL epitopes of human PV16 E7 protein, HIV IIIB gp120 P18, Nef, and reverse transcriptase (RT) proteins, and an HPV16 E7 linear B epitope. The CTL epitopes were restricted by three different MHC class 1 alleles (H-2(b), H-2(d), HLA-A*0201). The chimeric L1 protein assembled into VLPs when expressed in SF-9 cells by recombinant baculovirus. After immunization of mice with polytope VLPs in the absence of adjuvant, serum antibodies were detected which reacted with both polytope VLPs and wild-type BPV1L1 VLPs, in addition to the HPV16E7 linear B cell epitope. CTL precursors specific for the HPV16 E7, HIV P18, and RT CTL epitopes were also detected in the spleen of immunized mice. Polytope VLPs can thus deliver multiple B and T epitopes as immunogens to the MHC class I and class II pathways, extending the utility of VLPs as self-adjuvanting immunogen delivery systems. (C) 2000 Academic Press.

On the theory of mixing-frequency electron spin-echo envelope modulation spectroscopy

It has recently been stated that the parametrization of the time variables in the one-dimensional (I-D) mixing-frequency electron spin-echo envelope modulation (MIF-ESEEM) experiment is incorrect and hence the wrong frequencies for correlated nuclear transitions are predicted. This paper is a direct response to such a claim, its purpose being to show that the parametrization in land 2-D MIF-ESEEM experiments possesses the same form as that used in other 4-pulse incrementation schemes and predicts the same correlation frequencies. We show that the parametrization represents a shearing transformation of the 2-D time-domain and relate the resulting frequency domain spectrum to the HYSCORE spectrum in terms of a skew-projection. It is emphasized that the parametrization of the time-domain variables may be chosen arbitrarily and affects neither the computation of the correct nuclear frequencies nor the resulting resolution. The usefulness or otherwise of the MIF parameters \gamma\ > 1 is addressed, together with the validity of the original claims of the authors with respect to resolution enhancement in cases of purely homogeneous and inhomogeneous broadening. Numerical simulations are provided to illustrate the main points.