Fuel supplier selection for large scale UK bioenergy schemes

This article presents a potential method to assist developers of future bioenergy schemes when selecting from available suppliers of biomass materials. The method aims to allow tacit requirements made on biomass suppliers to be considered at the design stage of new developments. The method used is a combination of the Analytical Hierarchy Process and the Quality Function Deployment methods (AHP-QFD). The output of the method is a ranking and relative weighting of the available suppliers which could be used to improve optimization algorithms such as linear and goal programming. The paper is at a conceptual stage and no results have been obtained. The aim is to use the AHP-QFD method to bridge the gap between treatment of explicit and tacit requirements of bioenergy schemes; allowing decision makers to identify the most successful supply strategy available.

Bit-error probability for direct detection of optical RZ signal degraded by ASE noise and timing jitter

We analyze theoretically the interplay between optical return-to-zero signal degradation due to timing jitter and additive amplified-spontaneous-emission noise. The impact of these two factors on the performance of a square-law direct detection receiver is also investigated. We derive an analytical expression for the bit-error probability and quantitatively determine the conditions when the contributions of the effects of timing jitter and additive noise to the bit error rate can be treated separately. The analysis of patterning effects is also presented. © 2007 IEEE.

Some Comments on Q-Irresolute and Quasi-Irresolute Functions

The aim of this paper is to continue the study of θ-irresolute and quasi-irresolute functions as well as to give an example of a function which is θ-irresolute but neither quasi-irresolute nor an R-map and thus give an answer to a question posed by Ganster, Noiri and Reilly. We prove that RS-compactness is preserved under open, quasi-irresolute surjections.

Relação entre função sexual, sintomatologia depressiva e qualidade de vida em mulheres grávidas

INTRODUCTION: Human sexuality is recognized as one of the pillars of quality of life. In women, sexual function is influenced throughout life by many factors that can lead to the appearance of changes in the cycle of sexual response, and hence the quality of life (QOL). Pregnancy is a period of change, leaving them physically and mentally vulnerable, which may affect sexual function and quality of life during pregnancy. OBJECTIVE: To investigate the relationship between sexual function, presence of depressive symptoms and quality of life in pregnant women. METHODS: The study included 207 pregnant women attending prenatal examination of the Maternity Divine Love, Parnamirim / RN and the participants of the Course for Pregnant Women of the Department of Physical Therapy at UFRN (central campus). Initially it was applied, a questionnaire containing questions about sociodemographic, gynecological and obstetric data, as well as body and sexual self-knowledge. Sexual function was assessed using the Sexual Function Index Female (Female Sexual Function Index - FSFI). To assess the quality of life, we used the Quality Index Ferrans Life & Powers mom. The presence of depressive symptoms was verified by applying the Beck Depression Inventory. The Shapiro-Wilk test for normality was carried variables, Mann-Whitney test for carrying out the comparisons and the Wilcoxon test for comparing the monthly sexual frequency before and during pregnancy. Multiple linear regression was used to verify the relationship between sexual function, depressive symptoms and quality of life. We used the Spearman correlation to check correlation between the variables. Ap value <0.05 was adopted. RESULTS: Sexual function and depressive symptoms were related quality of life (R2 = 0.30, p <0.001). Depression had a moderate negative correlation with quality of life (0.53; p <0.001), whereas sexual function showed a positive correlation with low quality of life (0.22; p = 0.001). The planning of pregnancy, education and income shown to influence depression scores. With respect to sexual function, it was seen that during pregnancy, a reduction in the monthly frequency of sexual partner (Z = -10.56; p <0.001). Among the sexual domain, just the pain, showed a statistically significant difference compared between the second and third quarter (Z = -1.91, p <0.05). The score of the quality of life of women with sexual dysfunction was xvii significantly lower than that pregnant women without dysfunction (Z = -2.87, p = 0.004). Conclusion: Sexual function and the presence of depressive symptoms are related to the quality of life of pregnant women.

Occurrence and Function of Hoogsteen Base Pairs in Nucleic Acids

Nucleic acids (DNA and RNA) play essential roles in the central dogma of biology for the storage and transfer of genetic information. The unique chemical and conformational structures of nucleic acids – the double helix composed of complementary Watson-Crick base pairs, provide the structural basis to carry out their biological functions. DNA double helix can dynamically accommodate Watson-Crick and Hoogsteen base-pairing, in which the purine base is flipped by ~180° degrees to adopt syn rather than anti conformation as in Watson-Crick base pairs. There is growing evidence that Hoogsteen base pairs play important roles in DNA replication, recognition, damage or mispair accommodation and repair. Here, we constructed a database for existing Hoogsteen base pairs in DNA duplexes by a structure-based survey from the Protein Data Bank, and structural analyses based on the resulted Hoogsteen structures revealed that Hoogsteen base pairs occur in a wide variety of biological contexts and can induce DNA kinking towards the major groove. As there were documented difficulties in modeling Hoogsteen or Watson-Crick by crystallography, we collaborated with the Richardsons’ lab and identified potential Hoogsteen base pairs that were mis-modeled as Watson-Crick base pairs which suggested that Hoogsteen can be more prevalent than it was thought to be. We developed solution NMR method combined with the site-specific isotope labeling to characterize the formation of, or conformational exchange with Hoogsteen base pairs in large DNA-protein complexes under solution conditions, in the absence of the crystal packing force. We showed that there are enhanced chemical exchange, potentially between Watson-Crick and Hoogsteen, at a sharp kink site in the complex formed by DNA and the Integration Host Factor protein. In stark contrast to B-form DNA, we found that Hoogsteen base pairs are strongly disfavored in A-form RNA duplex. Chemical modifications N1-methyl adenosine and N1-methyl guanosine that block Watson-Crick base-pairing, can be absorbed as Hoogsteen base pairs in DNA, but rather potently destabilized A-form RNA and caused helix melting. The intrinsic instability of Hoogsteen base pairs in A-form RNA endows the N1-methylation as a functioning post-transcriptional modification that was known to facilitate RNA folding, translation and potentially play roles in the epitranscriptome. On the other hand, the dynamic property of DNA that can accommodate Hoogsteen base pairs could be critical to maintaining the genome stability.

Crack propagation in non-homogenous materials: Evaluation of mixed-mode SIFs, T-stress and kinking angle using a variant EFG method

A new variant of the Element-Free Galerkin (EFG) method, that combines the diffraction method, to characterize the crack tip solution, and the Heaviside enrichment function for representing discontinuity due to a crack, has been used to model crack propagation through non-homogenous materials. In the case of interface crack propagation, the kink angle is predicted by applying the maximum tangential principal stress (MTPS) criterion in conjunction with consideration of the energy release rate (ERR). The MTPS criterion is applied to the crack tip stress field described by both the stress intensity factor (SIF) and the T-stress, which are extracted using the interaction integral method. The proposed EFG method has been developed and applied for 2D case studies involving a crack in an orthotropic material, crack along an interface and a crack terminating at a bi-material interface, under mechanical or thermal loading; this is done to demonstrate the advantages and efficiency of the proposed methodology. The computed SIFs, T-stress and the predicted interface crack kink angles are compared with existing results in the literature and are found to be in good agreement. An example of crack growth through a particle-reinforced composite materials, which may involve crack meandering around the particle, is reported.

Eigenfunctions and fundamental solutions of the fractional Laplace and Dirac operators: the Riemann-Liouville case

In this paper we study eigenfunctions and fundamental solutions for the three parameter fractional Laplace operator $\Delta_+^{(\alpha,\beta,\gamma)}:= D_{x_0^+}^{1+\alpha} +D_{y_0^+}^{1+\beta} +D_{z_0^+}^{1+\gamma},$ where $(\alpha, \beta, \gamma) \in \,]0,1]^3$, and the fractional derivatives $D_{x_0^+}^{1+\alpha}$, $D_{y_0^+}^{1+\beta}$, $D_{z_0^+}^{1+\gamma}$ are in the Riemann-Liouville sense. Applying operational techniques via two-dimensional Laplace transform we describe a complete family of eigenfunctions and fundamental solutions of the operator $\Delta_+^{(\alpha,\beta,\gamma)}$ in classes of functions admitting a summable fractional derivative. Making use of the Mittag-Leffler function, a symbolic operational form of the solutions is presented. From the obtained family of fundamental solutions we deduce a family of fundamental solutions of the fractional Dirac operator, which factorizes the fractional Laplace operator. We apply also the method of separation of variables to obtain eigenfunctions and fundamental solutions.

Estudo clínico-laboratorial e histopatológico de cães naturalmente infectados por Leishmania chagasi com diferentes graus de manifestação física

Canine Visceral Leishmania (CVL) is an important zoonotic disease that has a world wide distribution and has a large impact on public health on the American Continent, especially in Brazil, where the nature of endemic diseases in humans affects a large part of the nation. The influence of the prevalence of CVL in the increased rate of human cases in endemic areas and in the unleashing of epidemic outbreaks shows the need for a more profound understanding, that would generate significant advances in the current measures used to control the reservoirs of sickness that are practiced by the Programa Nacional de Vigilância e Controle da Leishmaniose Visceral. The present work describes and compares the clinical-laboratorial and histopathological findings of twenty-three dogs that were naturally infected by Leishmania chagasi, from endemic areas in metropolitan Natal, Rio Grande do Norte, Brazil. These animals, that were selected and given physical and serological exams (IFI and ELISA rK-39), were classified according to the degree of clinical severity and had blood samples drawn (whole blood and serum) for a complete hemogram and a coagulogram to be done as well as biochemical tests for kidney and liver function. The confirmation of infection by L. chagasi was done after the euthanasia of the animals, through the direct demonstration of the parasite in the impression of the spleen and liver crowned with GIEMSA and through a cultivation by means of NNN/Schneider. According to the clinical evaluation, the animals were classified as asymptomatic (7), oligosymptomatic (7) and polysymptomatic (9). Among the animals that were chosen to be autopsied, there were 2 asymptomatic, 3 oligosymptomatic and 3 polysymptomatic, for the purpose of studying their histopathology, having collected fragments of the spleen, liver, kidneys and skin and were fixed in 10% tamponed formol. The comparison between the average parameters of the clinical-laboratory tested animals in the groups was done through the Student t test (a<0.05). The main clinical signals observed were lymphadenomegaly, alopecy, dermatitis, exfoliation, cutaneous ulcers, onicogriphosis and emaciation. The main clinical-laboratorial alterations established, mainly in the polysymptomatic group, were anemia, hyperproteinemia, hyperglobulinemia, alterations in the albumin/globulin ratio and increased ALT activity. Renal alterations were not verified (urea and creatinine levels were normal). Thrombocytopenia was observed in three clinical groups. However, the other indicators of coagulation function (TAP and TTPA) did not have abnormal variations. There were inflammatory infiltrations and leishmania amastigotes in the skin of polysymptomatic dogs, however, they were not found in the skin of asymptomatic animals. Hypertrophy and hyperplasia of the phagocyte mononuclear system, leishmania amastigote parasites were found in the macrophages, extramedullary hematopoiesis and degenerative alterations were detected in the spleen and liver of 8 of the animals submitted to histopathological exams. In accord with these results, it was demonstrated that the expected alterations in the hematological and biochemical parameters in function of their viscerotropic nature of CVL are mainly observed in the more advanced stages of the disease. The absence of inflammatory infiltration and parasite load in the skin suggest that infected animals without symptoms may have an importance irrelevant to the infectiousness of the vector

Desenvolvimento de sistema de análise e descrição de funções para uma PME

Mestrado em Gestão de Recursos Humanos

Comparison of Photovoltaic panel’s standard and simplified models

Shockley diode equation is basic for single diode model equation, which is overly used for characterizing the photovoltaic cell output and behavior. In the standard equation, it includes series resistance (Rs) and shunt resistance (Rsh) with different types of parameters. Maximum simulation and modeling work done previously, related to single diode photovoltaic cell used this equation. However, there is another form of the standard equation which has not included Series Resistance (Rs) and Shunt Resistance (Rsh) yet, as the Shunt Resistance is much bigger than the load resistance and the load resistance is much bigger than the Series Resistance. For this phenomena, very small power loss occurs within a photovoltaic cell. This research focuses on the comparison of two forms of basic Shockley diode equation. This analysis describes a deep understanding of the photovoltaic cell, as well as gives understanding about Series Resistance (Rs) and Shunt Resistance (Rsh) behavior in the Photovoltaic cell. For making estimation of a real time photovoltaic system, faster calculation is needed. The equation without Series Resistance and Shunt Resistance is appropriate for the real time environment. Error function for both Series resistance (Rs) and Shunt resistances (Rsh) have been analyzed which shows that the total system is not affected by this two parameters' behavior.

Sensitivity analysis of parameters of a photovoltaic cell under different condition

A photovoltaic cell is a component which converts light energy into electrical energy. Different environmental parameters and internal parameters have a great impact on the output of the photovoltaic cell. To identify its characteristics and estimate the output, the well known Shockley diode equation is used. This equation contains all the parameters, as one environmental and different internal. The properties of these parameters were studied and their sensitivity have been analyzed through the use of an error function; this error function allows the study of the behaviour of the parameters and their characteristics against the output of the photovoltaic cell through the analysis of its curves giving the sensitivity of the different parameters to the output of the photovoltaic cell. Using these results the impact of the parameters of the photovoltaic cell has been clearly identified. White noise is included both with the ideal values and the simulation and the ideal value is imposed to get the real time environment flavor. This work analyses both systems with and without white noise.

Similarity solution of axisymmetric flow in porous media

Applications of the axisymmetric Boussinesq equation to groundwater hydrology and reservoir engineering have long been recognised. An archetypal example is invasion by drilling fluid into a permeable bed where there is initially no such fluid present, a circumstance of some importance in the oil industry. It is well known that the governing Boussinesq model can be reduced to a nonlinear ordinary differential equation using a similarity variable, a transformation that is valid for a certain time-dependent flux at the origin. Here, a new analytical approximation is obtained for this case. The new solution,, which has a simple form, is demonstrated to be highly accurate. (c) 2005 Elsevier Ltd. All rights reserved.

Study of the voltammetric behaviour of metam and its application to an amperometric flow system

The electrochemical behaviour of the pesticide metam (MT) at a glassy carbon working electrode (GCE) and at a hanging mercury drop electrode (HMDE) was investigated. Different voltammetric techniques, including cyclic voltammetry (CV) and square wave voltammetry (SWV), were used. An anodic peak (independent of pH) at +1.46 V vs AgCl/Ag was observed in MTaqueous solution using the GCE. SWV calibration curves were plotted under optimized conditions (pH 2.5 and frequency 50 Hz), which showed a linear response for 17–29 mg L−1. Electrochemical reduction was also explored, using the HMDE. A well defined cathodic peak was recorded at −0.72 V vs AgCl/ Ag, dependent on pH. After optimizing the operating conditions (pH 10.1, frequency 150 Hz, potential deposition −0.20 V for 10 s), calibration curves was measured in the concentration range 2.5×10−1 to 1.0 mg L−1 using SWV. The electrochemical behaviour of this compound facilitated the development of a flow injection analysis (FIA) system with amperometric detection for the quantification of MT in commercial formulations and spiked water samples. An assessment of the optimal FIA conditions indicated that the best analytical results were obtained at a potential of +1.30 V, an injection volume of 207 μL and an overall flow rate of 2.4 ml min−1. Real samples were analysed via calibration curves over the concentration range 1.3×10−2 to 1.3 mg L−1. Recoveries from the real samples (spiked waters and commercial formulations) were between 97.4 and 105.5%. The precision of the proposed method was evaluated by assessing the relative standard deviation (RSD %) of ten consecutive determinations of one sample (1.0 mg L−1), and the value obtained was 1.5%.

Multiscale descriptions of density-driven flow instabilities in porous media

Les instabilités engendrées par des gradients de densité interviennent dans une variété d'écoulements. Un exemple est celui de la séquestration géologique du dioxyde de carbone en milieux poreux. Ce gaz est injecté à haute pression dans des aquifères salines et profondes. La différence de densité entre la saumure saturée en CO2 dissous et la saumure environnante induit des courants favorables qui le transportent vers les couches géologiques profondes. Les gradients de densité peuvent aussi être la cause du transport indésirable de matières toxiques, ce qui peut éventuellement conduire à la pollution des sols et des eaux. La gamme d'échelles intervenant dans ce type de phénomènes est très large. Elle s'étend de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères à laquelle interviennent les phénomènes à temps long. Une reproduction fiable de la physique par la simulation numérique demeure donc un défi en raison du caractère multi-échelles aussi bien au niveau spatial et temporel de ces phénomènes. Il requiert donc le développement d'algorithmes performants et l'utilisation d'outils de calculs modernes. En conjugaison avec les méthodes de résolution itératives, les méthodes multi-échelles permettent de résoudre les grands systèmes d'équations algébriques de manière efficace. Ces méthodes ont été introduites comme méthodes d'upscaling et de downscaling pour la simulation d'écoulements en milieux poreux afin de traiter de fortes hétérogénéités du champ de perméabilité. Le principe repose sur l'utilisation parallèle de deux maillages, le premier est choisi en fonction de la résolution du champ de perméabilité (grille fine), alors que le second (grille grossière) est utilisé pour approximer le problème fin à moindre coût. La qualité de la solution multi-échelles peut être améliorée de manière itérative pour empêcher des erreurs trop importantes si le champ de perméabilité est complexe. Les méthodes adaptatives qui restreignent les procédures de mise à jour aux régions à forts gradients permettent de limiter les coûts de calculs additionnels. Dans le cas d'instabilités induites par des gradients de densité, l'échelle des phénomènes varie au cours du temps. En conséquence, des méthodes multi-échelles adaptatives sont requises pour tenir compte de cette dynamique. L'objectif de cette thèse est de développer des algorithmes multi-échelles adaptatifs et efficaces pour la simulation des instabilités induites par des gradients de densité. Pour cela, nous nous basons sur la méthode des volumes finis multi-échelles (MsFV) qui offre l'avantage de résoudre les phénomènes de transport tout en conservant la masse de manière exacte. Dans la première partie, nous pouvons démontrer que les approximations de la méthode MsFV engendrent des phénomènes de digitation non-physiques dont la suppression requiert des opérations de correction itératives. Les coûts de calculs additionnels de ces opérations peuvent toutefois être compensés par des méthodes adaptatives. Nous proposons aussi l'utilisation de la méthode MsFV comme méthode de downscaling: la grille grossière étant utilisée dans les zones où l'écoulement est relativement homogène alors que la grille plus fine est utilisée pour résoudre les forts gradients. Dans la seconde partie, la méthode multi-échelle est étendue à un nombre arbitraire de niveaux. Nous prouvons que la méthode généralisée est performante pour la résolution de grands systèmes d'équations algébriques. Dans la dernière partie, nous focalisons notre étude sur les échelles qui déterminent l'évolution des instabilités engendrées par des gradients de densité. L'identification de la structure locale ainsi que globale de l'écoulement permet de procéder à un upscaling des instabilités à temps long alors que les structures à petite échelle sont conservées lors du déclenchement de l'instabilité. Les résultats présentés dans ce travail permettent d'étendre les connaissances des méthodes MsFV et offrent des formulations multi-échelles efficaces pour la simulation des instabilités engendrées par des gradients de densité. - Density-driven instabilities in porous media are of interest for a wide range of applications, for instance, for geological sequestration of CO2, during which CO2 is injected at high pressure into deep saline aquifers. Due to the density difference between the C02-saturated brine and the surrounding brine, a downward migration of CO2 into deeper regions, where the risk of leakage is reduced, takes place. Similarly, undesired spontaneous mobilization of potentially hazardous substances that might endanger groundwater quality can be triggered by density differences. Over the last years, these effects have been investigated with the help of numerical groundwater models. Major challenges in simulating density-driven instabilities arise from the different scales of interest involved, i.e., the scale at which instabilities are triggered and the aquifer scale over which long-term processes take place. An accurate numerical reproduction is possible, only if the finest scale is captured. For large aquifers, this leads to problems with a large number of unknowns. Advanced numerical methods are required to efficiently solve these problems with today's available computational resources. Beside efficient iterative solvers, multiscale methods are available to solve large numerical systems. Originally, multiscale methods have been developed as upscaling-downscaling techniques to resolve strong permeability contrasts. In this case, two static grids are used: one is chosen with respect to the resolution of the permeability field (fine grid); the other (coarse grid) is used to approximate the fine-scale problem at low computational costs. The quality of the multiscale solution can be iteratively improved to avoid large errors in case of complex permeability structures. Adaptive formulations, which restrict the iterative update to domains with large gradients, enable limiting the additional computational costs of the iterations. In case of density-driven instabilities, additional spatial scales appear which change with time. Flexible adaptive methods are required to account for these emerging dynamic scales. The objective of this work is to develop an adaptive multiscale formulation for the efficient and accurate simulation of density-driven instabilities. We consider the Multiscale Finite-Volume (MsFV) method, which is well suited for simulations including the solution of transport problems as it guarantees a conservative velocity field. In the first part of this thesis, we investigate the applicability of the standard MsFV method to density- driven flow problems. We demonstrate that approximations in MsFV may trigger unphysical fingers and iterative corrections are necessary. Adaptive formulations (e.g., limiting a refined solution to domains with large concentration gradients where fingers form) can be used to balance the extra costs. We also propose to use the MsFV method as downscaling technique: the coarse discretization is used in areas without significant change in the flow field whereas the problem is refined in the zones of interest. This enables accounting for the dynamic change in scales of density-driven instabilities. In the second part of the thesis the MsFV algorithm, which originally employs one coarse level, is extended to an arbitrary number of coarse levels. We prove that this keeps the MsFV method efficient for problems with a large number of unknowns. In the last part of this thesis, we focus on the scales that control the evolution of density fingers. The identification of local and global flow patterns allows a coarse description at late times while conserving fine-scale details during onset stage. Results presented in this work advance the understanding of the Multiscale Finite-Volume method and offer efficient dynamic multiscale formulations to simulate density-driven instabilities. - Les nappes phréatiques caractérisées par des structures poreuses et des fractures très perméables représentent un intérêt particulier pour les hydrogéologues et ingénieurs environnementaux. Dans ces milieux, une large variété d'écoulements peut être observée. Les plus communs sont le transport de contaminants par les eaux souterraines, le transport réactif ou l'écoulement simultané de plusieurs phases non miscibles, comme le pétrole et l'eau. L'échelle qui caractérise ces écoulements est définie par l'interaction de l'hétérogénéité géologique et des processus physiques. Un fluide au repos dans l'espace interstitiel d'un milieu poreux peut être déstabilisé par des gradients de densité. Ils peuvent être induits par des changements locaux de température ou par dissolution d'un composé chimique. Les instabilités engendrées par des gradients de densité revêtent un intérêt particulier puisque qu'elles peuvent éventuellement compromettre la qualité des eaux. Un exemple frappant est la salinisation de l'eau douce dans les nappes phréatiques par pénétration d'eau salée plus dense dans les régions profondes. Dans le cas des écoulements gouvernés par les gradients de densité, les échelles caractéristiques de l'écoulement s'étendent de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères sur laquelle interviennent les phénomènes à temps long. Etant donné que les investigations in-situ sont pratiquement impossibles, les modèles numériques sont utilisés pour prédire et évaluer les risques liés aux instabilités engendrées par les gradients de densité. Une description correcte de ces phénomènes repose sur la description de toutes les échelles de l'écoulement dont la gamme peut s'étendre sur huit à dix ordres de grandeur dans le cas de grands aquifères. Il en résulte des problèmes numériques de grande taille qui sont très couteux à résoudre. Des schémas numériques sophistiqués sont donc nécessaires pour effectuer des simulations précises d'instabilités hydro-dynamiques à grande échelle. Dans ce travail, nous présentons différentes méthodes numériques qui permettent de simuler efficacement et avec précision les instabilités dues aux gradients de densité. Ces nouvelles méthodes sont basées sur les volumes finis multi-échelles. L'idée est de projeter le problème original à une échelle plus grande où il est moins coûteux à résoudre puis de relever la solution grossière vers l'échelle de départ. Cette technique est particulièrement adaptée pour résoudre des problèmes où une large gamme d'échelle intervient et évolue de manière spatio-temporelle. Ceci permet de réduire les coûts de calculs en limitant la description détaillée du problème aux régions qui contiennent un front de concentration mobile. Les aboutissements sont illustrés par la simulation de phénomènes tels que l'intrusion d'eau salée ou la séquestration de dioxyde de carbone.

Determination of Minoxidil by Bleaching the Permanganate Carrier Solution in a Flow-Based Spectrophotometric System

The determination of minoxidil (MX) with potassium permanganate as a carrier in a flow injection method is described. The detection at 550nm was linear from 1.0x10-5 to 5.0x10-4mol L-1. The limit of detection (3 sigma/slope) was 8.92x10-6mol L-1, with an analytical frequency of 32h-1. The proposed method was applied to commercial samples, with recoveries from 104.7 to 106.4%. Comparison with the HPLC procedure reveled relative errors from 0.48 to 1.4%, and the results agreed within a 95% confidence level.