Stabilized Petrov-Galerkin methods for the convection-diffusion-reaction and the Helmholtz equations

We present two new stabilized high-resolution numerical methods for the convection–diffusion–reaction (CDR) and the Helmholtz equations respectively. The work embarks upon a priori analysis of some consistency recovery procedures for some stabilization methods belonging to the Petrov–Galerkin framework. It was found that the use of some standard practices (e.g. M-Matrices theory) for the design of essentially non-oscillatory numerical methods is not feasible when consistency recovery methods are employed. Hence, with respect to convective stabilization, such recovery methods are not preferred. Next, we present the design of a high-resolution Petrov–Galerkin (HRPG) method for the 1D CDR problem. The problem is studied from a fresh point of view, including practical implications on the formulation of the maximum principle, M-Matrices theory, monotonicity and total variation diminishing (TVD) finite volume schemes. The current method is next in line to earlier methods that may be viewed as an upwinding plus a discontinuity-capturing operator. Finally, some remarks are made on the extension of the HRPG method to multidimensions. Next, we present a new numerical scheme for the Helmholtz equation resulting in quasi-exact solutions. The focus is on the approximation of the solution to the Helmholtz equation in the interior of the domain using compact stencils. Piecewise linear/bilinear polynomial interpolation are considered on a structured mesh/grid. The only a priori requirement is to provide a mesh/grid resolution of at least eight elements per wavelength. No stabilization parameters are involved in the definition of the scheme. The scheme consists of taking the average of the equation stencils obtained by the standard Galerkin finite element method and the classical finite difference method. Dispersion analysis in 1D and 2D illustrate the quasi-exact properties of this scheme. Finally, some remarks are made on the extension of the scheme to unstructured meshes by designing a method within the Petrov–Galerkin framework.

Diffusion of ICT related problems among students : the teachers??? experience

Resumen tomado de la publicaci??n

Continuity of attractors for parabolic problems with localized large diffusion

In this paper we study the continuity of asymptotics of semilinear parabolic problems of the form u(t) - div(p(x)del u) + lambda u =f(u) in a bounded smooth domain ohm subset of R `` with Dirichlet boundary conditions when the diffusion coefficient p becomes large in a subregion ohm(0) which is interior to the physical domain ohm. We prove, under suitable assumptions, that the family of attractors behave upper and lower semicontinuously as the diffusion blows up in ohm(0). (c) 2006 Elsevier Ltd. All rights reserved.

Detection of Mycoplasma hyopneumoniae by polymerase chain reaction in swine presenting respiratory problems

A Nested-PCR (N-PCR) tem como objetivo melhorar a sensibilidade do diagnóstico direto da Pneumonia Enzoótica Suína, pois o isolamento do Mycoplasma hyopneumoniae é trabalhoso tornando-se inviável na rotina. Neste trabalho, foi realizado um projeto piloto para a otimização da técnica de N-PCR, utilizando três variáveis: tipo de amostra biológica, meio de transporte da amostra e método de extração do DNA, utilizando oito animais. Os resultados obtidos foram empregados no segundo experimento para a validação do teste utilizando 40 animais. Os resultados obtidos, pela otimização da N-PCR, neste trabalho, permite sugerir esta prova como método de diagnóstico de rotina no monitoramento das infecções por Mycoplasma hyopneumoniae em granjas de suínos.

Positive solutions for parametric Dirichlet problems with indefinite potential and superdiffusive reaction

We consider a parametric semilinear Dirichlet problem driven by the Laplacian plus an indefinite unbounded potential and with a reaction of superdifissive type. Using variational and truncation techniques, we show that there exists a critical parameter value λ_{∗}>0 such that for all λ> λ_{∗} the problem has least two positive solutions, for λ= λ_{∗} the problem has at least one positive solutions, and no positive solutions exist when λ∈(0,λ_{∗}). Also, we show that for λ≥ λ_{∗} the problem has a smallest positive solution.

Nonlinear periodic problems with a jumping reaction

We consider a periodic problem driven by the scalar $p-$Laplacian and with a jumping (asymmetric) reaction. We prove two multiplicity theorems. The first concerns the nonlinear problem ($1

Asymptotic stability of a coupled Advection-Diffusion-Reaction system arising in bioreactor processes.

In this work, we perform an asymptotic analysis of a coupled system of two Advection-Diffusion-Reaction equations with Danckwerts boundary conditions, which models the interaction between a microbial population (e.g., bacterias), called biomass, and a diluted organic contaminant (e.g., nitrates), called substrate, in a continuous flow bioreactor. This system exhibits, under suitable conditions, two stable equilibrium states: one steady state in which the biomass becomes extinct and no reaction is produced, called washout, and another steady state, which corresponds to the partial elimination of the substrate. We use the method of linearization to give sufficient conditions for the asymptotic stability of the two stable equilibrium configurations. Finally, we compare our asymptotic analysis with the usual asymptotic analysis associated to the continuous bioreactor when it is modeled with ordinary differential equations.

Reduced models of chemical reaction in chaotic flows

We describe and evaluate two reduced models for nonlinear chemical reactions in a chaotic laminar flow. Each model involves two separate steps to compute the chemical composition at a given location and time. The “manifold tracking model” first tracks backwards in time a segment of the stable manifold of the requisite point. This then provides a sample of the initial conditions appropriate for the second step, which requires solving one-dimensional problems for the reaction in Lagrangian coordinates. By contrast, the first step of the “branching trajectories model” simulates both the advection and diffusion of fluid particles that terminate at the appropriate point; the chemical reaction equations are then solved along each of the branched trajectories in a second step. Results from each model are compared with full numerical simulations of the reaction processes in a chaotic laminar flow.

Reaction time assessment in children with ADHD

Attention deficit, impulsivity and hyperactivity are the cardinal features of attention deficit hyperactivity disorder (ADHD) but executive function (EF) disorders, as problems with inhibitory control, working memory and reaction time, besides others EFs, may underlie many of the disturbs associated with the disorder. OBJECTIVE: To examine the reaction time in a computerized test in children with ADHD and normal controls. METHOD: Twenty-three boys (aged 9 to 12) with ADHD diagnosis according to Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition, 2000 (DSM-IV) criteria clinical, without comorbidities, Intelligence Quotient (IQ) >89, never treated with stimulant and fifteen normal controls, age matched were investigated during performance on a voluntary attention psychophysical test. RESULTS: Children with ADHD showed reaction time higher than normal controls. CONCLUSION: A slower reaction time occurred in our patients with ADHD. This findings may be related to problems with the attentional system, that could not maintain an adequate capacity of perceptual input processes and/or in motor output processes, to respond consistently during continuous or repetitive activity.

Mechanistic studies of the 'blue' Cu enzyme, bilirubin oxidase, as a highly efficient electrocatalyst for the oxygen reduction reaction

The 'blue copper' enzyme bilirubin oxidase from Myrothecium verrucaria shows significantly enhanced adsorption on a pyrolytic graphite 'edge' (PGE) electrode that has been covalently modified with naphthyl-2-carboxylate functionalities by diazonium coupling. Modified electrodes coated with bilirubin oxidase show electrocatalytic voltammograms for the direct, four-electron reduction of O(2) by bilirubin oxidase with up to four times the current density of an unmodified PGE electrode. Electrocatalytic voltammograms measured with a rapidly rotating electrode (to remove effects of O(2) diffusion limitation) have a complex shape (an almost linear dependence of current on potential below pH 6) that is similar regardless of how PGE is chemically modified. Importantly, the same waveform is observed if bilirubin oxidase is adsorbed on Au(111) or Pt(111) single-crystal electrodes (at which activity is short-lived). The electrocatalytic behavior of bilirubin oxidase, including its enhanced response on chemically-modified PGE, therefore reflects inherent properties that do not depend on the electrode material. The variation of voltammetric waveshapes and potential-dependent (O(2)) Michaelis constants with pH and analysis in terms of the dispersion model are consistent with a change in rate-determining step over the pH range 5-8: at pH 5, the high activity is limited by the rate of interfacial redox cycling of the Type 1 copper whereas at pH 8 activity is much lower and a sigmoidal shape is approached, showing that interfacial electron transfer is no longer a limiting factor. The electrocatalytic activity of bilirubin oxidase on Pt(111) appears as a prominent pre-wave to electrocatalysis by Pt surface atoms, thus substantiating in a single, direct experiment that the minimum overpotential required for O(2) reduction by the enzyme is substantially smaller than required at Pt. At pH 8, the onset of O(2) reduction lies within 0.14 V of the four-electron O(2)/2H(2)O potential.

A wavelet-based adaptive technique for adsorption problems involving steep gradients

A general, fast wavelet-based adaptive collocation method is formulated for heat and mass transfer problems involving a steep moving profile of the dependent variable. The technique of grid adaptation is based on sparse point representation (SPR). The method is applied and tested for the case of a gas–solid non-catalytic reaction in a porous solid at high Thiele modulus. Accurate and convergent steep profiles are obtained for Thiele modulus as large as 100 for the case of slab and found to match the analytical solution.

Computationally efficient solution techniques for adsorption problems involving steep gradients in bidisperse particles

A piecewise uniform fitted mesh method turns out to be sufficient for the solution of a surprisingly wide variety of singularly perturbed problems involving steep gradients. The technique is applied to a model of adsorption in bidisperse solids for which two fitted mesh techniques, a fitted-mesh finite difference method (FMFDM) and fitted mesh collocation method (FMCM) are presented. A combination (FMCMD) of FMCM and the DASSL integration package is found to be most effective in solving the problems. Numerical solutions (FMFDM and FMCMD) were found to match the analytical solution when the adsorption isotherm is linear, even under conditions involving steep gradients for which global collocation fails. In particular, FMCMD is highly efficient for macropore diffusion control or micropore diffusion control. These techniques are simple and there is no limit on the range of the parameters. The techniques can be applied to a variety of adsorption and desorption problems in bidisperse solids with non-linear isotherm and for arbitrary particle geometry.