Some experiments with high order compact methods using a computer algebra software - Part II (non-uniform grid)

We generalize a procedure proposed by Mancera and Hunt [P.F.A. Mancera, R. Hunt, Some experiments with high order compact methods using a computer algebra software-Part 1, Appl. Math. Comput., in press, doi: 10.1016/j.amc.2005.05.015] for obtaining a compact fourth-order method to the steady 2D Navier-Stokes equations in the streamfunction formulation-vorticity using the computer algebra system Maple, which includes conformal mappings and non-uniform grids. To analyse the procedure we have solved a constricted stepped channel problem, where a fine grid is placed near the re-entrant corner by transformation of the independent variables. (c) 2006 Elsevier 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.

Traveling waves in the Lethargic Crab Disease

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

Assessing the effects of vector control on dengue transmission

The dengue virus is an arbovirus transmitted to humans through mosquito Aedes aegypti. To describe the dynamics of dengue disease within these two populations we develop a compartment model taking into account chemical controls and mechanical control applied on the mosquitoes, which are the currently available controlling mechanisms to prevent dengue disease. To mimic seasonal variations, some parameters of the model are allowed to depend on time in order to divide a calendar year in favorable and unfavorable periods regarded to the development of vector population. (c) 2007 Elsevier B.V. All rights reserved.

Simulation results and applications of an advection bounded scheme to practical flows

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

Convexity of the extreme zeros of Gegenbauer and Laguerre polynomials

Let C-n(lambda)(x), n = 0, 1,..., lambda > -1/2, be the ultraspherical (Gegenbauer) polynomials, orthogonal. in (-1, 1) with respect to the weight function (1 - x(2))(lambda-1/2). Denote by X-nk(lambda), k = 1,....,n, the zeros of C-n(lambda)(x) enumerated in decreasing order. In this short note, we prove that, for any n is an element of N, the product (lambda + 1)(3/2)x(n1)(lambda) is a convex function of lambda if lambda greater than or equal to 0. The result is applied to obtain some inequalities for the largest zeros of C-n(lambda)(x). If X-nk(alpha), k = 1,...,n, are the zeros of Laguerre polynomial L-n(alpha)(x), also enumerated in decreasing order, we prove that x(n1)(lambda)/(alpha + 1) is a convex function of alpha for alpha > - 1. (C) 2002 Published by Elsevier B.V. B.V.

State constraints in optimal impulsive controls

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

Connection coefficients and zeros of orthogonal polynomials

We discuss an old theorem of Obrechkoff and some of its applications. Some curious historical facts around this theorem are presented. We make an attempt to look at some known results on connection coefficients, zeros and Wronskians of orthogonal polynomials from the perspective of Obrechkoff's theorem. Necessary conditions for the positivity of the connection coefficients of two families of orthogonal polynomials are provided. Inequalities between the kth zero of an orthogonal polynomial p(n)(x) and the largest (smallest) zero of another orthogonal polynomial q(n)(x) are given in terms of the signs of the connection coefficients of the families {p(n)(x)} and {q(n)(x)}, An inequality between the largest zeros of the Jacobi polynomials P-n((a,b)) (x) and P-n((alpha,beta)) (x) is also established. (C) 2001 Elsevier B.V. B.V. All rights reserved.

Zeros of Jacobi functions of second kind

The number of zeros in (- 1, 1) of the Jacobi function of second kind Q(n)((alpha, beta)) (x), alpha, beta > - 1, i.e. The second solution of the differential equation(1 - x(2))y (x) + (beta - alpha - (alpha + beta + 2)x)y' (x) + n(n + alpha + beta + 1)y(x) = 0,is determined for every n is an element of N and for all values of the parameters alpha > - 1 and beta > - 1. It turns out that this number depends essentially on alpha and beta as well as on the specific normalization of the function Q(n)((alpha, beta)) (x). Interlacing properties of the zeros are also obtained. As a consequence of the main result, we determine the number of zeros of Laguerre's and Hermite's functions of second kind. (c) 2005 Elsevier B.V. All rights reserved.

The Q-D algorithm for transforming series expansions into a corresponding continued fraction: an extension to cope with zero coefficients

An algorithm for deriving a continued fraction that corresponds to two series expansions simultaneously, when there are zero coefficients in one or both series, is given. It is based on using the Q-D algorithm to derive the corresponding fraction for two related series, and then transforming it into the required continued fraction. Two examples are given. (C) 2003 Elsevier B.V. All rights reserved.

On a symmetry in strong distributions

A strong Stieltjes distribution d psi(t) is called symmetric if it satisfies the propertyt(omega) d psi(beta(2)/t) = -(beta(2)/t)(omega) d psi(t), for t is an element of (a, b) subset of or equal to (0, infinity), 2 omega is an element of Z, and beta > 0.In this article some consequences of symmetry on the moments, the orthogonal L-polynomials and the quadrature formulae associated with the distribution are given. (C) 1999 Elsevier B.V. B.V. All rights reserved.

Edge detection and noise removal by use of a partial differential equation with automatic selection of parameters

This work deals with noise removal by the use of an edge preserving method whose parameters are automatically estimated, for any application, by simply providing information about the standard deviation noise level we wish to eliminate. The desired noiseless image u(x), in a Partial Differential Equation based model, can be viewed as the solution of an evolutionary differential equation u t(x) = F(u xx, u x, u, x, t) which means that the true solution will be reached when t ® ¥. In practical applications we should stop the time ''t'' at some moment during this evolutionary process. This work presents a sufficient condition, related to time t and to the standard deviation s of the noise we desire to remove, which gives a constant T such that u(x, T) is a good approximation of u(x). The approach here focused on edge preservation during the noise elimination process as its main characteristic. The balance between edge points and interior points is carried out by a function g which depends on the initial noisy image u(x, t0), the standard deviation of the noise we want to eliminate and a constant k. The k parameter estimation is also presented in this work therefore making, the proposed model automatic. The model's feasibility and the choice of the optimal time scale is evident through out the various experimental results.

Monotonicity of zeros of Laguerre polynomials

Denote by X(nk)(alpha), k = 1, ..., n, the zeros of the Laguerre polynomial L(n)((alpha))(X). We establish monotonicity with respect to the parameter at of certain functions involving X(nk)(alpha). As a consequence we obtain sharp upper bounds for the largest zero of L(n)((alpha))(X). (C) 2009 Elsevier B.V. All rights reserved.

Laguerre polynomials as Jensen polynomials of Laguerre-Polya entire functions

We prove that the only Jensen polynomials associated with an entire function in the Laguerre-Polya class that are orthogonal are the Laguerre polynomials. (C) 2009 Elsevier B.V. All rights reserved.