Aproximação para Problema de Controle Ótimo Impulsivo e Problema de Tempo Mínimo sobre Domínios Estratificados

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

An O(N log N) algorithm for shape modeling.

We present a shape-recovery technique in two dimensions and three dimensions with specific applications in modeling anatomical shapes from medical images. This algorithm models extremely corrugated structures like the brain, is topologically adaptable, and runs in O(N log N) time, where N is the total number of points in the domain. Our technique is based on a level set shape-recovery scheme recently introduced by the authors and the fast marching method for computing solutions to static Hamilton-Jacobi equations.

Gravidade de Lovelock e a correspondência AdS/CFT

A correspondência AdS/CFT é uma notável ferramenta no estudo de teorias de gauge fortemente acopladas que podem ser mapeadas em uma descrição gravitacional dual fracamente acoplada. A correspondência é melhor entendida no limite em que ambos $N$ e $\\lambda$, o rank do grupo de gauge e o acoplamento de \'t Hooft da teoria de gauge, respectivamente, são infinitos. Levar em consideração interações com termos de curvatura de ordem superior nos permite considerar correções de $\\lambda$ finito. Por exemplo, a primeira correção de acoplamento finito para supergravidade tipo IIB surge como um termo de curvatura com forma esquemática $\\alpha\'^3 R^4$. Neste trabalho investigamos correções de curvatura no contexto da gravidade de Lovelock, que é um cenário simples para investigar tais correções pois as suas equações de movimento ainda são de segunda ordem em derivadas. Esse cenário também é particularmente interessante do ponto de vista da correspondência AdS/CFT devido a sua grande classe de soluções de buracos negros assintoticamente AdS. Consideramos um sistema de gravidade AdS-axion-dilaton em cinco dimensões com um termo de Gauss-Bonnet e encontramos uma solução das equações de movimento, o que corresponde a uma black brane exibindo uma anisotropia espacial, onde a fonte da anisotropia é um campo escalar linear em uma das coordenadas espaciais. Estudamos suas propriedades termodinâmicas e realizamos a renormalização holográfica usando o método de Hamilton-Jacobi. Finalmente, usamos a solução obtida como dual gravitacional de um plasma anisotrópico fortemente acoplado com duas cargas centrais independentes, $a eq c$. Calculamos vários observáveis relevantes para o estudo do plasma, a saber, a viscosidade de cisalhamento sobre densidade de entropia, a força de arrasto, o parâmetro de jet quenching, o potencial entre um par quark-antiquark e a taxa de produção de fótons.

Sistemas integrables discretos, polinomios matriciales ortogonales y problemas de Riemann-Hilbert

El propósito de esta tesis doctoral es el estudio de la conexión, mediante el problema de Riemann-Hilbert, entre sistemas discretos y la teoría de polinomios matriciales ortogonales. La investigación de los modelos integrables se originó en la Mecánica Clásica, en relación a la resolución de las ecuaciones de Newton [2]. Los trabajos de Liouville, Hamilton, Jacobi y otros sentaron las bases de los sistemas integrables como prototipos modelos resolubles por cuadraturas, v.g., por integración directa [7]. Hay una cantidad importante de investigación dedicada a los aspectos geométricos de los sistemas clásicos integrables y superintegrables [66], [82], especialmente en relación a la separación de variables de la ecuación de Hamilton-Jacobi [75]. Fue la aplicación, en la segunda mitad del siglo pasado, de la transformada espectral inversa para la resolución del problema de Cauchy de la ecuación de Korteweg-de Vries [42, 43] la que marcó el inicio de una nueva etapa en este campo, el del estudio de sistemas integrables con un número infinito de grados de libertad, que generalmente se expresan en términos de jerarquías de ecuaciones no lineales en derivadas parciales. Particularmente reseñable, por su aplicación en la hidrodinámica y en la óptica cuántica, es la aparición de las soluciones a un número de solitones arbitrario. En las últimas tres décadas ha habido un importante interés por el estudio de modelos discretos, v.g., sistemas dinámicos de nidos en un retículo de puntos, y expresados en términos de ecuaciones no lineales en diferencia parciales. Muchas de las técnicas encontradas en el mundo continuo se extendieron a este nuevo contexto discreto. Hay dos razones fundamentales para este interés...

A Mean Value Theorem for non Differentiable Mappings in Banach Spaces

We prove that if f is a real valued lower semicontinuous function on a Banach space X and if there exists a C^1, real valued Lipschitz continuous function on X with bounded support and which is not identically equal to zero, then f is Lipschitz continuous of constant K provided all lower subgradients of f are bounded by K. As an application, we give a regularity result of viscosity supersolutions (or subsolutions) of Hamilton-Jacobi equations in infinite dimensions which satisfy a coercive condition. This last result slightly improves some earlier work by G. Barles and H. Ishii.

A Novel High Order Space-Time Spectral Method for the Time Fractional Fokker--Planck Equation

The fractional Fokker-Planck equation is an important physical model for simulating anomalous diffusions with external forces. Because of the non-local property of the fractional derivative an interesting problem is to explore high accuracy numerical methods for fractional differential equations. In this paper, a space-time spectral method is presented for the numerical solution of the time fractional Fokker-Planck initial-boundary value problem. The proposed method employs the Jacobi polynomials for the temporal discretization and Fourier-like basis functions for the spatial discretization. Due to the diagonalizable trait of the Fourier-like basis functions, this leads to a reduced representation of the inner product in the Galerkin analysis. We prove that the time fractional Fokker-Planck equation attains the same approximation order as the time fractional diffusion equation developed in [23] by using the present method. That indicates an exponential decay may be achieved if the exact solution is sufficiently smooth. Finally, some numerical results are given to demonstrate the high order accuracy and efficiency of the new numerical scheme. The results show that the errors of the numerical solutions obtained by the space-time spectral method decay exponentially.

New doubly periodic waves of the (2+1)-dimensional double sine-Gordon equation

New exact solutions of the (2 + 1)-dimensional double sine-Gordon equation are studied by introducing the modified mapping relations between the cubic nonlinear Klein-Gordon system and double sine-Gordon equation. Two arbitrary functions are included into the Jacobi elliptic function solutions. New doubly periodic wave solutions are obtained and displayed graphically by proper selections of the arbitrary functions.

Length-weight relationship in Nandus nandus (Hamilton) (Perciformes-Teleostei)

Length-weight relationship of the freshwater fish of Pamba River, Nandus nandus (Ham) has been worked out. The results showed that the slope values and elevations were not significant and hence a combined regression equation has been calculated for both the sexes (Log W=2.4130 Log L -0.3306). The’t’ test analyses were conducted and found that the growth departs significantly from the isometric growth. Thus the formula W=aL super(n) has to be applied in calculating the length-weight relationship of this species.

Sound Synthesis for Contact-Driven Musical Instruments via Discretisation of Hamilton's Equations

Physical modelling of musical instruments involves studying nonlinear interactions between parts of the instrument. These can pose several difficulties concerning the accuracy and stability of numerical algorithms. In particular, when the underlying forces are non-analytic functions of the phase-space variables, a stability proof can only be obtained in limited cases. An approach has been recently presented by the authors, leading to unconditionally stable simulations for lumped collision models. In that study, discretisation of Hamilton’s equations instead of the usual Newton’s equation of motion yields a numerical scheme that can be proven to be energy conserving. In this paper, the above approach is extended to collisions of distributed objects. Namely, the interaction of an ideal string with a flat barrier is considered. The problem is formulated within the Hamiltonian framework and subsequently discretised. The resulting nonlinearmatrix equation can be shown to possess a unique solution, that enables the update of the algorithm. Energy conservation and thus numerical stability follows in a way similar to the lumped collision model. The existence of an analytic description of this interaction allows the validation of the model’s accuracy. The proposed methodology can be used in sound synthesis applications involving musical instruments where collisions occur either in a confined (e.g. hammer-string interaction, mallet impact) or in a distributed region (e.g. string-bridge or reed-mouthpiece interaction).

Solution properties of the de Branges differential recurrence equation

In this 1984 proof of the Bieberbach and Milin conjectures de Branges used a positivity result of special functions which follows from an identity about Jacobi polynomial sums thas was published by Askey and Gasper in 1976. The de Branges functions Tn/k(t) are defined as the solutions of a system of differential recurrence equations with suitably given initial values. The essential fact used in the proof of the Bieberbach and Milin conjectures is the statement Tn/k(t)<=0. In 1991 Weinstein presented another proof of the Bieberbach and Milin conjectures, also using a special function system Λn/k(t) which (by Todorov and Wilf) was realized to be directly connected with de Branges', Tn/k(t)=-kΛn/k(t), and the positivity results in both proofs Tn/k(t)<=0 are essentially the same. In this paper we study differential recurrence equations equivalent to de Branges' original ones and show that many solutions of these differential recurrence equations don't change sign so that the above inequality is not as surprising as expected. Furthermore, we present a multiparameterized hypergeometric family of solutions of the de Branges differential recurrence equations showing that solutions are not rare at all.

Using the Cayley-Hamilton theorem with N-partitioned matrices

A definition is given for the characteristic equation of anN-partitioned matrix. It is then proved that this matrix satisfies its own characteristic equation. This can then be regarded as a version of the Cayley-Hamilton theorem, of use withN-dimensional systems.

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.

Fractional Extensions of Jacobi Polynomials and Gauss Hypergeometric Function

2000 Mathematics Subject Classification: 26A33, 33C45

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.