CNMAC 2012 - 34º Congresso Nacional de Matemática Aplicada e Computacional

Conferências

Data e Horário Nome e Instituição área
Terça-feira, 18/09
16:30-17:30 h
C7 - Nick Kingsbury
University of Cambridge
Processamento de Sinais
Email: ngk@eng.cam.ac.uk
Segunda-feira, 17/09
16:30-17:30 h
C2 - Edriss Saleh Titi
U.C.Irvine and Weizmann Institute, Israel
PDE, Dynamical Systems, Numerical Analysis, Turbulence
Email: etiti@math.uci.edu
Quarta-feira, 19/09
16:30-17:30 h
C11 - I-Shih Liu
Instituto de Matemática - UFRJ
Mecânica do Contínuo/Matemática Aplicada
Email: liu@im.ufrj.br
Segunda-feira, 17/09
16:30-17:30 h
C3 - José Raimundo Coelho
Agencia Espacial Brasileira
Análise Numérica, Diretrizes do Programa Espacial Brasileiro e a Matemática Aplicada no Brasil
Email: presidente@aeb.gov.br
Terça-feira, 18/09
11:30-12:30 h
C6 - Jose Eduardo Souza de Cursi
INSA de Rouen
Modelagem Estocastica; Otimização Robusta
Email: souza@insa-rouen.fr
Quarta-feira, 19/09
16:30-17:30 h
C13 - Leslei Richard Foulds (Les Foulds)
Universidade federal de Goiás
Otimização Combinatória
Email: lesfoulds@gmail.com
Quinta-feira, 20/09
16:30-17:30 h
C15 - Yurico Y. Baldin
UFSCAR
O que aprendemos com as Oficinas para professores do projeto Klein de Matemática em português: olhares distintos sobre os artigos Klein.
Email: yuriko@dm.ufscar.br
Sexta-feira, 21/09
11:30-12:30 h
C16 - Pedro Dias
LNCC
Matemática Aplicada à Meteorologia
Email: pldsdias@lncc.br
Terça-feira, 18/09
11:30-12:30 h
C5 - Antonio C. G. Leitão
UFSC
Algoritimos
Email: acgleitao@gmail.com
Quinta-feira, 20/09
16:30-17:30 h
C14 - Marie Farge
ENS- Paris
Wavelets
Email: farge@lmd.ens.fr
Quarta-feira, 19/09
11:30-12:30 h
C10 -Luis G. Nonato
ICMC
Visualização e Processamento Geométrico
Email: gnonato@icmc.usp.br
Quarta-feira, 19/09
11:30-12:30 h
C9 - Weldon Lodwick
Univ. Colorado
Otimização Fuzzy
Email: weldon.lodwick@ucdenver.edu
Terça-feira, 18/09
16:30-17:30 h
C8 - Alexei A. Mailybaev
IMPA
Matemática Aplicada; Dinâmica dos Fluidos; Leis de Conservação
Email: a.mailybaev@gmail.com
Quarta-feira, 19/09
16:30-17:30 h
C12 - Osvaldo Rosso
UFAL
Análise de séries temporais, dinâmica não-linear, teoria wavelet, Teoria da Informação e suas aplicações às áreas da Biologia e Medicina
Email: oarosso@gmail.com
Terça-feira, 18/09
11:30-12:30 h
C4 - Antonio F. B. de A. Prado Engenharia Aeroespacial, com ênfase em Trajetórias Espaciais
Email: prado@dem.inpe.br
Sexta-feira, 21/09
11:30-12:30 h
C17 - Teresa E. Pérez
Universidade de Granada
Approximation Theory, Orthogonal Polinomials
Email: tperez@ugr.es
Segunda-feira, 17/09
11:00-12:00 h
C1 - José Mário Martinez
IMECC- Unicamp
Otimização
Email: martinez@ime.unicamp.br

 

Nick Kingsbury

Conferências

Dual-Tree Complex Wavelets - their key properties and a range of image-processing applications

Resumo: We will describe the Dual-Tree Complex Wavelet Transform (DT CWT), a form of discrete wavelet transform which generates complex coefficients by using two trees of wavelet filters in parallel to obtain their real and imaginary parts. This introduces limited redundancy (2^m : 1 for m-dimensional signals) and the Hilbert-pair property is designed into the two wavelet bases such that the transform provides approximate shift invariance and directionally selective filters (properties lacking in the traditional wavelet transform), while preserving the usual properties of perfect reconstruction and computational efficiency with good well-balanced frequency responses. Energy is preserved in the Q-shift wavelet filters, resulting in tight-frame properties too. In the decade since the dual tree was first proposed, it has been applied to images and 3D datasets and has been shown to outperform traditional wavelets in a range of applications, which will be summarized.

Complex-valued wavelets, the dual tree, and Hilbert pairs: why these lead to shift invariance and directional m-D wavelets?

Resumo: This mini-course will look in detail at why the dual-tree produces Hilbert pairs of wavelets and how to design good Q-shift filters. It will also discuss applications such as displacement estimation and 3D data registration, as well as deconvolution with sparse regularization.

We will discuss some of the interesting mathematical properties of the dual tree wavelet construction, when it is used to generate analytic complex-valued wavelets. We shall show how the requirement for analyticity (zero energy at negative frequencies) implies that the wavelet bases from the two parts of the dual tree should form a Hilbert transform pair, and how this leads to the half-sample delay condition between the lowpass (scaling function) filters of each tree. It will then be shown how this results in approximate shift invariance of the complex wavelet transform, and the elegant Q-shift filter solution to this condition will be given. We will then discuss how the analytic nature of the wavelets in 1-D allows separable wavelet filters in higher dimensional spaces to be highly directionally selective.

Finally we will demonstrate how the analytic wavelets can be used for a range of applications. In particular, we shall discuss how the approximately linear phase vs displacement properties of the complex wavelets allow for very efficient displacement estimation in 2D and 3D spaces. We shall also show how their tight-frame, sparsifying and decorrelation properties produce effective iterative deconvolution algorithms when combined with sparsity-seeking regularization methods.

Biografia

Nick Kingsbury is Professor of Signal Processing and head of the Signal Processing and Communications research group at the University of Cambridge, Department of Engineering. He has worked in the areas of digital communications, audio analysis and coding, and image processing. He has developed the dual-tree complex wavelet transform and is especially interested in the application of complex wavelets and related multiscale and multiresolution methods to the analysis of images and 3-D datasets.

Edriss Saleh Titi

Conferências

The Navier-Stokes, Euler and Other Related Equations

Resumo: In this talk I will survey the status of, and the most recent advances concerning, the questions of global regularity of solutions to the three-dimensional Navier-Stokes and Euler equations of incompressible uids. Furthermore, I will also present recent global regularity results concerning certain three-dimensional geophysical ows, including the three-dimensional viscous \primitive equations" of oceanic and atmospheric dynamics.

Jose Eduardo Souza de Cursi

Conferências

Quantificação da incerteza de soluções de problemas numéricos com parâmetros aleatórios

Resumo: Os problemas práticos envolvem parâmetros que devem ser estimados a partir de testes, observações ou experimentos et que são, portanto, incertos ou, no mínimo afetados por erros e variabilidade. Por exemplo, os parâmetros físicos de uma estrutura ou de um material, assim como os parâmetros de modelos econômicos sofrem de diferentes variabilidades causadas por condições operacionais e do meio ambiente. Tais variabilidades introduzem incertezas no funcionamento do sistema, de modo que é importante quantificá-las por razões de segurança. De maneira análoga, a variabilidade das condições futuras faz com que as predições da resposta do sistema contenham incertezas cuja quantificação é indispensável para uma avaliação correta dos riscos.

Na prática, a quantificação dos riscos resulta de uma análise da distribuição de probabilidades das soluções das equações de um modelo ou da otimização de um critério. Assim, é necessário determinar – de uma forma cuja utilização prática seja simples – a distribuição de probabilidades das variáveis e interesse.

Esta conferencia apresenta, de forma os principais métodos de determinação numérica de tais distribuições de probabilidade, com ilustrações em situações simples, mas significativas e exemplos de aplicações práticas.

I-Shih Liu

Conferências

Mathematical modeling of large deformation in salt tectonics

Resumo: We consider instability of a two-layered solid body of a denser material on top of a lighter one. This problem is widely known to geoscientist in sediment-salt migration as salt diapirism. In the literature, this problem has often been treated as Rayleigh-Taylor instability in viscous fluids instead of solid bodies.

To treat problems of large deformation, we propose a new algorithm for large deformation of solid bodies in general. It is a successive incremental linear approximation based on the well-known idea of small deformation superposed on large defornation at every numerical step. Numerical example for bending of a rectangular block into a circular section is presented for comparision with the exact solution.

Modeling sediment-salt as elastic-viscoelastic solid body, the proposed method is used for numerical simulation of salt tectonics showing the formation of salt diapirs. Experimental varification of material models is also shown.

Yurico Y. Baldin

Conferências

O que aprendemos com as Oficinas para professores do projeto Klein de Matemática em português: olhares distintos sobre os artigos Klein.

Resumo: O objetivo desta palestra é apresentar o Projeto Klein de Matemática em Português, um projeto da SBM, com apoio da SBMAC, SBEM, SBHMat e OBMEP, como uma contribuição organizada para o projeto internacional The Klein Project for 21st century da ICMI-IMU. A apresentação conterá os princípios do projeto para que todos conheçam seus objetivos e as metodologias, assim como fazer uma retrospectiva dos resultados alcançados pelas atividades realizadas desde seu início em 2009. A palestra apresentará também um olhar crítico sobre os impactos das suas diferentes atividades sobre a formação de licenciandos e docentes do Ensino Médio por meio dos pequenos textos Klein que constituem eixo norteador do projeto.

Les Foulds

Conferências

THE STEINER TREE PROBLEM FOR THE N-CUBE AND EXTENSIONS TO MATHEMATICAL BIOLOGY

Resumo: For a metric space (S, d), define a weighted graph G= G(S, d) with vertex set S so that each edge {s, t} has weight d(s, t). For a finite subset X of S, a minimum spanning tree T(X) for X is a tree (that is, a connected, acyclic subgraph) with vertex set X such that the sum of the weights of T(X) is a minimum. A Steiner minimal tree S(X) for X is a tree having the minimum possible length over all trees in G that contain X in their vertex sets.

For a fixed alphabet A = {0, 1} and a fixed positive integer N, let d denote the Hamming distance on AN, that is, d((a1,...,aN), (a1',...,aN')) is equal to the number of indices i such that ai ≠ ai'. Denote AN by QN.. The graph G = G(QN, d) is the 1-skeleton of the N-cube. We establish that the complexity of finding a Steiner minimal tree for QN is NP-complete.

We also discuss the extension of the above problem to the case where A has a given, finite number of symbols. For this more general problem, we discuss approaches for tree construction, the question of tree minimality and a tree comparison metric.

Finally, we mention how the results may be useful in the task of constructing evolutionary trees from molecular databases.

Biografia

Tenho trabalhado na área de pesquisa de operações, com ênfase em engenharia de transito. Tenho experiência de ensino da optimização, pesquisa de operações, ciência de gerência, estatísticas, ciência de decisão, modelos de gerência quantitativos e fabricação. Até junho 2007, fui professor na Universidade de Waikato, Hamilton, Nova Zelândia. Recebi um BSc em matemática, 1970 e um MSc com honras em matemática, 1972, ambos na Universidade de Auckland, Nova Zelândia. Fiz doutorado em engenharia industrial na Virginia Tech, USA, 1974. Minha área de pesquisa principal é o desenvolvimento dos modelos e das técnicas de optimização para modelagem e simulação de transito. Publiquei mais de 200 artigos de pesquisa nesses tópicos, em jornais academicos. Por exemplo: Nature, Operations Research, Management Science, the European Journal of Operational Research, the Journal of the Operational Research Society, the Annals of Operations Research e the Asia-Pacific Journal of Operational Research. Desde setembro 2007, estou vivendo permanentemente em GO, Brasil e continuando minhas pesquisas em modelagem e simulação de transito, sendo Professor Visitante - Instituto de Informática - UFG - Goiânia - GO.

Pedro Leite da Silva Dias

Conferências

Processamento de Alto Desempenho: papel, serviços e futuro

Resumo:A modelagem computacional é usada para compreender, analisar e resolver problemas científicos e tecnológicos e constitui-se numa grande área de pesquisa interdisciplinar. Ë baseada no desenvolvimento de modelos matemáticos, algorítimos e técnicas de processamento computacional para simular, testar hipóteses e prever a evolução de fenômenos e processos. Tem aplicações nas áreas de saúde, engenharias, ciências da terra, física, química, biologia e nas ciências humanas. A dimensão dos problemas tipicamente abordados na modelagem computacional cresceu substancialmente nas últimas décadas em função os seguintes aspectos: (a) capacidade de processamento (aumento da ordem de mil vezes entre 1997 e 2008, (b) capacidade de armazenamento (um milhão de vezes no mesmo período) e (c) avanço das tecnologias de informação e computação (redes mutiGigabit). O Brasil esteve razoavelmente bem colocado no ranking internacional no final dos anos 90 mas hoje está bem longe do desempenho computacional disponível nos demais países do BRIC. Mas existem "ilhas" no Brasil, como no caso da Previsão de Tempo e Clima (CPTEC) e na Petrobras. Serão discutidos alguns exemplos de aplicações de computação, de alto desempenho no Brasil, as necessidades atuais e perspectivas futuras do SINAPAD (Sistema Nacional de Processamento de Alto Desempenho).

Antonio C. G. Leitão

Conferências

Algoritmos iterativos para solução de problemas matemáticos desafiadores na sociedade contemporânea

Resumo:Na primeira parte da palestra será apresentada, por meio de exemplos, uma família de problemas matemáticos que se caracterizam por objetivar a determinação de um parâmetro em um sistema de equações (algébricas, diferenciais, integrais, etc). Tais problemas, denominados problemas inversos, oferecem os seguintes desafios:

i) as medições disponíveis do parâmetro desconhecido são indiretas e inexatas;
ii) o operador matemático que mapeia o parâmetro nos dados observados é instável (isto é, mal-posto).

Uma vez apresentados os problemas de interesse, a palestra será focada na discussão de 2 famílias de algoritmos iterativos que visam obter soluções aproximadas dos mesmos, a saber:

i) métodos tipo gradiente;
ii) métodos tipo Newton.

A análise de convergência desses métodos será apresentada de forma didática: Que perguntas devem ser respondidas para que se possa determinar e/ou quantificar a eficiência de um método iterativo para resolver problemas mal-condicionados?

Marie Farge

Conferências

How wavelets are used to study turbulent transport in fluid and plasma flows

Resumo:We will review wavelet methods to study fluid and plasma flows in the fully-developed turbulent regime for bounded and unbounded geometries. We will present new results concerning the emergence of dissipative structures in fluid and plasma turbulence. We will present a wavelet denoising method that we will apply to videos taken with a fast camera implemented in the tokamak Tore-Supra to study coherent structures observed in edge plasma. We will discuss perspectives concerning the use of wavelets to study turbulent transport in confined fluid and plasma flows and conjecture the existence of non viscous and non resistive dissipative structures in the fully-developed turbulent regime.

This work is done in collaboration with Romain Nguyen van yen (Freie Universität, Berlin) and Kai Schneider, Université de Provence, Marseille. It is financed by the Fédération de Recherche Fusion par Confinement Magnétique-ITER. It was published in Phys. Rev. Lett. (106, 184502, 2011) and in Nucl. Fusion (52, 013005, 2012). These papers can be downloaded from http://wavelets.ens.fr in 'Publications'.

Luis G. Nonato

Conferências

Multidimensional Projection for Visualizing and Interacting with High-Dimensional Data

Resumo: Interactive visualization has become a central component of most massive data analysis and exploration systems. However, enabling simple, intuitive and effective interactive visual resources is not a straightforward task, mainly when dealing with complex high-dimensional data. In this talk we will present the recent mathematical and computational tools we have developed towards manipulating large high-dimensional data sets. Much of our research relies on multidimensional projection methods designed to map high-dimensional instances of data onto low-dimensional visual spaces so as to preserve neighborhood structures as much as possible. The flexibility and usefulness of our framework will be demonstrated through a variety of applications such as interactive summarization of document collection, image colorization, vector field visualization, music playlist construction, visualization of web search results, among others. Challenges and future directions of this stimulating research field will also be discussed.

Weldon A. Lodwick

Conferências

Constrained Intervals, Interval Spaces, and Optimization Under Generalized Uncertainty

Resumo: Constrained intervals, intervals as a mapping from [0,1] to linear functions with non-negative slopes, and arithmetic on constrained intervals [1], generate a space that turns out to be a cancellative abelian monoid albeit with a richer set of properties than the usual (standard) interval arithmetic. This means that not only do we have the classical embedding as developed by H. Radström [3] and S. Markov [2] but the properties of these polynomials. We study a little of the geometry of the embedding of intervals into a quasi-linear vector space and some of the properties of the mapping of constrained intervals into a space of polynomials. Thus, there are three parts to this talk. (1) The representation of intervals as linear polynomials with non-negative slopes. (2) The algebraic structure of this new representation. (3) The link between (constrained) intervals and generalized uncertainty as they are used in optimization. The theoretical reason for considering a new representation of intervals is to have a formalization in (a subset of) polynomial space with the view to evaluate expressions (functions) of intervals. The theoretical reason for considering the algebraic structure of the embedding into a space with inverses is to solve equations and optimization problems.

Alexei Mailybaev

Conferências

Resonance in low-temperature oxidation waves for porous media

Resumo: We analyze traveling wave profiles possessing an internal resonance point for a class of systems of partial differential equations describing oxidation and vaporization of liquid fuel in a porous medium when an oxidizer (air) is injected. It is shown that the resonance is characterized by a saddle point of an associated vector field defined on a folded surface in state space.

We prove existence and uniqueness of this singularity for an open set of parameters. The singularity yields an extra restriction on wave parameters. This restriction is found explicitly in the physically relevant case of a small ratio of reaction/vaporization rates. We find the general large time asymptotic solutions of the problem as a sequence of waves, and show that resonance waves play a key role in determining these solutions. Efficiency of the resonant regime for oil recovery process is discussed.

Osvaldo A. Rosso

Conferências

GENERALIZED STATISTICAL COMPLEXITY MEASURE: A NEW TOOL FOR DYNAMICAL SYSTEMS.

Resumo: López-Ruiz, Mancini and Calbet have proposed a statistical complexity measure, based on the notion of "disequilibrium", as a quantifier of the degree of physical structure in a time series Given a probability distribution P and its associate information measure I[P], an amount of "disorder" associated to the state of a system, the LMC-statistical complexity, is the product of a normalized entropy H (normalized Shannon-entropy) times the disequilibrium Q, given by the Euclidean "distance" from P to the uniform distribution Pe. The statistical complexity vanishes both for a totally random process and for a purely periodic one. Martín, Plastino and Rosso improved on this measure by suitably modifying the distance-component (in the concomitant probability space). The obtained MPR-statistical complexity is (i) able to grasp essential details of the dynamics, (ii) an intensive quantity, and (iii) capable of discerning among different degrees of periodicity and chaos. The MPR-statistical complexity can be viewed as a functional that characterizes the probability distribution P associated to the time series generated by the dynamical system under study. It quantifies not only randomness but also the presence of correlational structures. In this seminar, selection of the information measure I and generalized disorder, as well as, election of distance D and generalized disequilibrium are reviewed. Evaluation of the probability distribution P associated to a dynamical system or time series under study is a physical problem. Additional improvements can be expected if the underlying probability distribution is "extracted" by more appropriate consideration regarding causal effects in the system's dynamics. Several well-known modelgenerated time series, usually regarded as being of either stochastic or chaotic nature, are analyzed. The main achievement of this approach is the possibility of clearly distinguishing between them in the H x MPR-Complexity representation space, something that is rather difficult otherwise. In addition, recent applications to time series from biological (epileptic EEG records, deformation of red blood cells, cancer progression, etc.) and physics systems (stochastic resonance, econophysics, literature, ENSO/El Niño, etc.) will be reviewer.

Antonio Fernando Bertachini de Almeida Prado

Conferências

TRANSFER ORBITS IN THE PLANAR RESTRICTED THREE-BODY PROBLEM

Resumo: In this research, the goal is to find transfer orbits for a spacecraft from one celestial body back to the same celestial body or to the Lagrangian points of that system. The mathematical model for the system is the planar restricted three-body problem. To avoid numerical problems during close approaches, the version of the model that has the global Lamaître regularization is used. First of all, attention is given to the family of transfer orbits involving no more than one revolution of the spacecraft. The systems under study are the ones with more important practical applications: the Earth-Sun and the Earth-Moon systems.

Five families of transfer orbits are found in the region studied and the results are plotted in terms of the true anomaly. The same plots also show the evolution of the Jacobian constant. A special effort is made to find transfer orbits with small *V.

After that, the problem of sending a spacecraft from the Earth to the Lagrangian points L4 and L5 is treated as a natural extension of the problem of sending a spacecraft from one body back to the same body. Two transfer orbits from the Earth to L4 and two transfer orbits from the Earth to L5 are found. Next, the numerical integration is extended beyond the desired Lagrangian point and it is found that, for all four orbits, the spacecraft passes near the Lagrangian points L3, L4 and L5 and comes back to the neighborhood of the Earth. Integration after this first close approach with the Earth has different results for each of the orbits studied, and it is necessary to describe them separately. In general, the orbits found here can be applied to:

i) Transfer a spacecraft between any two points in the group formed by the Earth and the Lagrangian points L3, L4, L5 (in the Earth-Sun system) with near-zero *V;

ii) Make a tour to the Lagrangian points for reconnaissance purposes with near-zero *V for the entire tour. The small relative velocities during the close approaches are ideal for the data acquisition phase or for a rendezvous with another spacecraft. There is also a possibility to recover the spacecraft after the tour (for re utilization and/or study), since it returns to the Earth's neighborhood;

iii) Build a cycler transportation system linking all the points involved or only two of them (as the two options for a cycler transportation system between the Earth and L5 or L4, that is shown in the present paper). In a system like that, a heavy spacecraft can stay in one of the orbits showed here and a small spacecraft can make a "taxi service" and rendezvous with the heavy vehicle to transport persons and/or materials to/from it. The "swing-by" with the Earth plays a very important role in these systems and will be explained in detail in the present paper.

Teresa E. Pérez

Conferências

On Multivariate Orthogonal Polynomials

Resumo: Classical orthogonal polynomials in one variable can be defined as orthogonal polynomials associated to a moment functional satisfying a Pearson differential equation. We extend this concept to several variables defining classical multivariate orthogonal polynomials as those associated to a moment functional satisfying a matrix analogue of the Pearson differential equation ([2]). In this talk, we present old and new results on classical multivariate orthogonal polynomials. In particular, classical multivariate orthogonal polynomials on the unit ball and some useful modifications will be studied. We will present the so–called Krall modification ([3]) obtained by adding a finite set of equally spaced mass points on the border of the unit ball, and some Sobolev modifications ([1], [4]), obtained by adding to the ball functional the evaluation of the gradients in one or several points. Sobolev orthogonal polynomials in several variables appears for the first time in [5] in a problem related to dwell time for polishing tools in fabricating optical surfaces.

José Mário Martinez

Conferências

Otimização em Estruturas Eletrônicas

Resumo: O Cálculo de Estruturas eletrônicas é um tema desafiante da matemática aplicada que envolve a resolução de problemas com enorme número de variáveis. Diversas técnicas de otimização são usadas para a resolução deste problema. Nesta palestra faremos uma revisão das técnicas mais bem sucedidas e relataremos nossas experiências na área.