Fourth-order method for solving the Navier-Stokes equations in a constricting channel

A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.

Structural and Biochemical Characterization of Peroxiredoxin Q beta from Xylella fastidiosa CATALYTIC MECHANISM and HIGH REACTIVITY

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

Avaliação de características morfométricas na relação solo-superfície da Bacia Hidrográfica do Córrego Rico, Jaboticabal (SP)

Com a análise das características morfométricas procura-se entender a relação solo-superfície, em decorrência dos processos erosivos sobre estruturas e litologias variadas. Neste trabalho, objetivou-se avaliar as características morfométricas do padrão de drenagem e do relevo em microbacias hidrográficas de 2ª ordem de magnitude, em Latossolo Vermelho-Escuro e Argissolo Vermelho-Amarelo, com vistas ao planejamento agroambiental da Bacia Hidrográfica do Córrego Rico, município de Jaboticabal, SP. O estudo foi baseado em fotografias aéreas, pancromáticas, verticais na escala 1:35.000 da região administrativa de Ribeirão Preto. Com técnicas de fotointerpretação, foram traçadas a rede de drenagem e respectivas microbacias hidrográficas, onde as variáveis selecionadas foram determinadas. As classes de risco de degradação, resultando no uso potencial do solo, foram determinadas de acordo com o coeficiente de rugosidade. As características do padrão de drenagem e do relevo repercutiram no comportamento hidrológico e litológico das microbacias hidrográficas, possibilitando a diferenciação entre as unidades de solos. A densidade de drenagem e a amplitude altimétrica foram as características que mais contribuíram para o poder discriminatório na relação solo-superfície. Na análise conjunta do potencial de uso do solo, verificou-se, principalmente, adequação para agricultura nos Latossolos e pecuária e, ou, reflorestamento nos Argissolos.

Morfometria de microbacias do Córrego Rico, afluente do Rio Mogi-Guaçu, Estado de São Paulo, Brasil

Este trabalho teve como objetivo avaliar as características morfométricas das microbacias (2ª, 3ª, 4ª e 5ª ordens de magnitude) da bacia hidrográfica do córrego Rico, sub-bacia do Rio Mogi-Guaçu, localizada na região administrativa de Ribeirão Preto, Estado de São Paulo, Brasil. Para tanto, foram determinados os parâmetros físicos e a configuração topográfica natural do sistema de drenagem. Os procedimentos para a obtenção dos dados foram fundamentados em técnicas de sensoriamento remoto e geoprocessamento. A partir da vetorização das cartas topográficas correspondentes à área de estudo, realizou-se a análise morfométrica quanto às características dimensionais, do padrão de drenagem e do relevo no sistema de informação geográfica ArcView. A microbacia é considerada de sexta ordem de magnitude, com área estimada de 542 km², com 85 microbacias de segunda ordem, 22 de terceira, sete de quarta ordem e duas de quinta. Utilizando o critério geométrico, na disposição fluvial das sub-bacias de cabeceiras observou-se a predominância dos modelos dendríticos e subdendríticos, enquanto a jusante predominava o modelo subparalelo, respectivamente, nas áreas de ocorrências dos arenitos Bauru e rochas efusivas básicas.

A simple and non-destructive model for individual leaf area estimation in citrus

Introduction. Leaf area is often related to plant growth, development, physiology and yield. Many non-destructive models have been proposed for leaf area estimation of several plant genotypes, demonstrating that leaf length, leaf width and leaf area are closely correlated. Thus, the objective of our study was to develop a reliable model for leaf area estimation from linear measurements of leaf dimensions for citrus genotypes. Materials and methods. Leaves of citrus genotypes were harvested, and their dimensions (length, width and area) were measured. Values of leaf area were regressed against length, width, the square of length, the square of width and the product (length x width). The most accurate equations, either linear or second-order polynomial, were regressed again with a new data set; then the most reliable equation was defined. Results and discussion. The first analysis showed that the variables length, width and the square of length gave better results in second-order polynomial equations, while the linear equations were more suitable and accurate when the width and the product (length x width) were used. When these equations were regressed with the new data set, the coefficient of determination (R(2)) and the agreement index 'd' were higher for the one that used the variable product (length x width), while the Mean Absolute Percentage Error was lower. Conclusion. The product of the simple leaf dimensions (length x width) can provide a reliable and simple non-destructive model for leaf area estimation across citrus genotypes.

Stochastic stability for Markovian jump linear systems associated with a finite number of jump times

This paper deals with a stochastic stability concept for discrete-time Markovian jump linear systems. The random jump parameter is associated to changes between the system operation modes due to failures or repairs, which can be well described by an underlying finite-state Markov chain. In the model studied, a fixed number of failures or repairs is allowed, after which, the system is brought to a halt for maintenance or for replacement. The usual concepts of stochastic stability are related to pure infinite horizon problems, and are not appropriate in this scenario. A new stability concept is introduced, named stochastic tau-stability that is tailored to the present setting. Necessary and sufficient conditions to ensure the stochastic tau-stability are provided, and the almost sure stability concept associated with this class of processes is also addressed. The paper also develops equivalences among second order concepts that parallels the results for infinite horizon problems. (C) 2003 Elsevier B.V. All rights reserved.

Optimal Conditions for Biomass and Recombinant Glycerol Kinase Production Using the Yeast Pichia pastoris

The extracellular glycerol kinase gene from Saccharomyces cerevisiae (GUT]) was cloned into the expression vector pPICZ alpha. A and integrated into the genome of the methylotrophic yeast Pichia pastoris X-33. The presence of the GUT1 insert was confirmed by PCR analysis. Four clones were selected and the functionality of the recombinant enzyme was assayed. Among the tested clones, one exhibited glycerol kinase activity of 0.32 U/mL, with specific activity of 0.025 U/mg of protein. A medium optimized for maximum biomass production by recombinant Pichia pastoris in shaker cultures was initially explored, using 2.31 % (by volume) glycerol as the carbon source. Optimization was carried out by response surface methodology (RSM). In preliminary experiments, following a Plackett-Burman design, glycerol volume fraction (phi(Gly)) and growth time (t) were selected as the most important factors in biomass production. Therefore, subsequent experiments, carried out to optimize biomass production, followed a central composite rotatable design as a function of phi(Gly) and time. Glycerol volume fraction proved to have a significant positive linear effect on biomass production. Also, time was a significant factor (at linear positive and quadratic levels) in biomass production. Experimental data were well fitted by a convex surface representing a second order polynomial model, in which biomass is a function of both factors (R(2)=0.946). Yield and specific activity of glycerol kinase were mainly affected by the additions of glycerol and methanol to the medium. The optimized medium composition for enzyme production was: 1 % yeast extract, 1 % peptone, 100 mM potassium phosphate buffer, pH=6.0, 1.34 % yeast nitrogen base (YNB), 4.10(-5) % biotin, 1 %, methanol and 1 %, glycerol, reaching 0.89 U/mL of glycerol kinase activity and 14.55 g/L of total protein in the medium after 48 h of growth.

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.

Determinação da solução ótima do problema de fluxo de potência ótimo via análise de sensibilidade

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

Black string corrections in variable tension braneworld scenarios

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

Third-body perturbation in the case of elliptic orbits for the disturbing body

This work presents a semi-analytical and numerical study of the perturbation caused in a spacecraft by a third-body using a double averaged analytical model with the disturbing function expanded in Legendre polynomials up to the second order. The important reason for this procedure is to eliminate terms due to the short periodic motion of the spacecraft and to show smooth curves for the evolution of the mean orbital elements for a long-time period. The aim of this study is to calculate the effect of lunar perturbations on the orbits of spacecrafts that are traveling around the Earth. An analysis of the stability of near-circular orbits is made, and a study to know under which conditions this orbit remains near circular completes this analysis. A study of the equatorial orbits is also performed. Copyright (C) 2008 R. C. Domingos et al.

Artificial neural networks for load flow and external equivalents studies

In this paper an artificial neural network (ANN) based methodology is proposed for (a) solving the basic load flow, (b) solving the load flow considering the reactive power limits of generation (PV) buses, (c) determining a good quality load flow starting point for ill-conditioned systems, and (d) computing static external equivalent circuits. An analysis of the input data required as well as the ANN architecture is presented. A multilayer perceptron trained with the Levenberg-Marquardt second order method is used. The proposed methodology was tested with the IEEE 30- and 57-bus, and an ill-conditioned 11-bus system. Normal operating conditions (base case) and several contingency situations including different load and generation scenarios have been considered. Simulation results show the excellent performance of the ANN for solving problems (a)-(d). (C) 2010 Elsevier B.V. All rights reserved.

Some experiments with high order compact methods using a computer algebra software - Part I

We consider a procedure for obtaining a compact fourth order method to the steady 2D Navier-Stokes equations in the streamfunction formulation using the computer algebra system Maple. The resulting code is short and from it we obtain the Fortran program for the method. To test the procedure we have solved many cavity-type problems which include one with an analytical solution and the results are compared with results obtained by second order central differences to moderate Reynolds numbers. (c) 2005 Elsevier B.V. All rights reserved.

An algorithm for arranging response surface designs in small blocks

We consider the problem of blocking response surface designs when the block sizes are prespecified to control variation efficiently and the treatment set is chosen independently of the block structure. We show how the loss of information due to blocking is related to scores defined by Mead and present an interchange algorithm based on scores to improve a given blocked design. Examples illustrating the performance of the algorithm are given and some comparisons with other designs are made. (C) 2000 Elsevier B.V. B.V. All rights reserved.

A study of a numerical solution of the steady two dimensions Navier-Stokes equations in a constricted channel problem by a compact fourth order method

We present a numerical solution for the steady 2D Navier-Stokes equations using a fourth order compact-type method. The geometry of the problem is a constricted symmetric channel, where the boundary can be varied, via a parameter, from a smooth constriction to one possessing a very sharp but smooth corner allowing us to analyse the behaviour of the errors when the solution is smooth or near singular. The set of non-linear equations is solved by the Newton method. Results have been obtained for Reynolds number up to 500. Estimates of the errors incurred have shown that the results are accurate and better than those of the corresponding second order method. (C) 2002 Elsevier B.V. All rights reserved.