Codomain Rigidity of the Dirichlet to Neumann Operator for the Riemannian Wave Equation

We study the Dirichlet to Neumann operator for the Riemannian wave equation on a compact Riemannian manifold. If the Riemannian manifold is modelled as an elastic medium, this operator represents the data available to an observer on the boundary of the manifold when the manifold is set into motion through boundary vibrations. We study the Dirichlet to Neumann operator when vibrations are imposed and data recorded on disjoint sets, a useful setting for applications. We prove that this operator determines the Dirichlet to Neumann operator where sources and observations are on the same set, provided a spectral condition on the Laplace-Beltrami operator for the manifold is satisfied. We prove this by providing an implementable procedure for determining a portion of the Riemannian manifold near the area where sources are applied. Drawing on established results, an immediate corollary is that a compact Riemannian manifold can be reconstructed from the Dirichlet to Neumann operator where sources and observations are on disjoint sets.

Codomain Rigidity of the Dirichlet to Neumann Operator for the Riemannian Wave Equation

We study the Dirichlet to Neumann operator for the Riemannian wave equation on a compact Riemannian manifold. If the Riemannian manifold is modelled as an elastic medium, this operator represents the data available to an observer on the boundary of the manifold when the manifold is set into motion through boundary vibrations. We study the Dirichlet to Neumann operator when vibrations are imposed and data recorded on disjoint sets, a useful setting for applications. We prove that this operator determines the Dirichlet to Neumann operator where sources and observations are on the same set, provided a spectral condition on the Laplace-Beltrami operator for the manifold is satisfied. We prove this by providing an implementable procedure for determining a portion of the Riemannian manifold near the area where sources are applied. Drawing on established results, an immediate corollary is that a compact Riemannian manifold can be reconstructed from the Dirichlet to Neumann operator where sources and observations are on disjoint sets.

Variationelle Zeitdiskretisierungen höherer Ordnung der Konvektions-Diffusions-Reaktions-Gleichungen in zeitabhängigen Gebieten

Wir betrachten zeitabhängige Konvektions-Diffusions-Reaktions-Gleichungen in zeitabhängi- gen Gebieten, wobei die Bewegung des Gebietsrandes bekannt ist. Die zeitliche Entwicklung des Gebietes wird durch die ALE-Formulierung behandelt, die die Nachteile der klassischen Euler- und Lagrange-Betrachtungsweisen behebt. Die Position des Randes und seine Geschwindigkeit werden dabei so in das Gebietsinnere fortgesetzt, dass starke Gitterdeformationen verhindert werden. Als Zeitdiskretisierungen höherer Ordnung werden stetige Galerkin-Petrov-Verfahren (cGP) und unstetige Galerkin-Verfahren (dG) auf Probleme in zeitabhängigen Gebieten angewendet. Weiterhin werden das C 1 -stetige Galerkin-Petrov-Verfahren und das C 0 -stetige Galerkin- Verfahren vorgestellt. Deren Lösungen lassen sich auch in zeitabhängigen Gebieten durch ein einfaches einheitliches Postprocessing aus der Lösung des cGP-Problems bzw. dG-Problems erhalten. Für Problemstellungen in festen Gebieten und mit zeitlich konstanten Konvektions- und Reaktionstermen werden Stabilitätsresultate sowie optimale Fehlerabschätzungen für die nachbereiteten Lösungen der cGP-Verfahren und der dG-Verfahren angegeben. Für zeitabhängige Konvektions-Diffusions-Reaktions-Gleichungen in zeitabhängigen Gebieten präsentieren wir konservative und nicht-konservative Formulierungen, wobei eine besondere Aufmerksamkeit der Behandlung der Zeitableitung und der Gittergeschwindigkeit gilt. Stabilität und optimale Fehlerschätzungen für die in der Zeit semi-diskretisierten konservativen und nicht-konservativen Formulierungen werden vorgestellt. Abschließend wird das volldiskretisierte Problem betrachtet, wobei eine Finite-Elemente-Methode zur Ortsdiskretisierung der Konvektions-Diffusions-Reaktions-Gleichungen in zeitabhängigen Gebieten im ALE-Rahmen einbezogen wurde. Darüber hinaus wird eine lokale Projektionsstabilisierung (LPS) eingesetzt, um der Konvektionsdominanz Rechnung zu tragen. Weiterhin wird numerisch untersucht, wie sich die Approximation der Gebietsgeschwindigkeit auf die Genauigkeit der Zeitdiskretisierungsverfahren auswirkt.

Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico

LINS, Filipe C. A. et al. Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico. In: CONGRESSO BRASILEIRO DE PESQUISA E DESENVOLVIMENTO EM PETRÓLEO E GÁS, 5. 2009, Fortaleza, CE. Anais... Fortaleza: CBPDPetro, 2009.

Funzioni di Sobolev

In questa tesi cercherò di analizzare le funzioni di Sobolev su R}^{n}, seguendo le trattazioni Measure Theory and Fine Properties of Functions di L.C. Evans e R.F.Gariepy e l'elaborato Functional Analysis, Sobolev Spaces and Partial Differential Equations di H. Brezis. Le funzioni di Sobolev si caratterizzano per essere funzioni con le derivate prime deboli appartenenti a qualche spazio L^{p}. I vari spazi di Sobolev hanno buone proprietà di completezza e compattezza e conseguentemente sono spesso i giusti spazi per le applicazioni di analisi funzionale. Ora, come vedremo, per definizione, l'integrazione per parti è valida per le funzioni di Sobolev. È, invece, meno ovvio che altre regole di calcolo siano allo stesso modo valide. Così, ho inteso chiarire questa questione di carattere generale, con particolare attenzione alle proprietà puntuali delle funzioni di Sobolev. Abbiamo suddiviso il lavoro svolto in cinque capitoli. Il capitolo 1 contiene le definizioni di base necessarie per la trattazione svolta; nel secondo capitolo sono stati derivati vari modi di approssimazione delle funzioni di Sobolev con funzioni lisce e sono state fornite alcune regole di calcolo per tali funzioni. Il capitolo 3 darà un' interpretazione dei valori al bordo delle funzioni di Sobolev utilizzando l'operatore Traccia, mentre il capitolo 4 discute l' estensione su tutto R^{n} di tali funzioni. Proveremo infine le principali disuguaglianze di Sobolev nel Capitolo 5.

Nonlinear Dynamics Near the Cell Membrane of Chara Corallina

The phenomenon of patterned distribution of pH near the cell membrane of the algae Chara corallina upon illumination is well-known. In this paper, we develop a mathematical model, based on the detailed kinetic analysis of proton fluxes across the cell membrane, to explain this phenomenon. The model yields two coupled nonlinear partial differential equations which describe the spatial dynamics of proton concentration changes and transmembrane potential generation. The experimental observation of pH pattern formation, its period and amplitude of oscillation, and also its hysteresis in response to changing illumination, are all reproduced by our model. A comparison of experimental results and predictions of our theory is made. Finally, a mechanism for pattern formation in Chara corallina is proposed.

Some Inverse Problems in Analysis and Geometry

Thesis (Ph.D.)--University of Washington, 2016-08

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.

Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico

LINS, Filipe C. A. et al. Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico. In: CONGRESSO BRASILEIRO DE PESQUISA E DESENVOLVIMENTO EM PETRÓLEO E GÁS, 5. 2009, Fortaleza, CE. Anais... Fortaleza: CBPDPetro, 2009.

Um problema de dispersão de poluentes em rios e canais por meio do método de separação de variáveis

Neste trabalho obtém-se uma solução analítica para a equação de advecção-difusão aplicada a problemas de dispersão de poluentes em rios e canais. Para tanto, consideram-se os casos unidimensionais e bidimensionais em regime transiente com coeficientes de difusividade e velocidades constantes. A abordagem utilizada para a resolução deste problema é o método de Separação de Variáveis. Os modelos resolvidos foram simulados utilizando o MatLab. Apresentam-se os resultados das simulações numéricas em formato gráfico. Os resultados de algumas simulações numéricas existem na literatura e puderam ser comparados. O modelo proposto mostrou-se coerente em relação aos dados considerados. Para outras simulações não foram encontrados comparativos na literatura, todavia esses problemas governados por equações diferenciais parciais, mesmo lineares, não são de fácil solução analítica. Sendo que, muitas delas representam importantes problemas de matemática e física, com diversas aplicações na engenharia. Dessa forma, é de grande importância a disponibilidade de um maior número de problemas-teste para avaliação de desempenho de formulações numéricas, cada vez mais eficazes, já que soluções analíticas oferecem uma base mais segura para comparação de resultados.

Eigenfunctions and fundamental solutions of the Caputo fractional Laplace and Dirac operators

In this paper, by using the method of separation of variables, we obtain eigenfunctions and fundamental solutions for the three parameter fractional Laplace operator defined via fractional Caputo derivatives. The solutions are expressed using the Mittag-Leffler function and we show some graphical representations for some parameters. A family of fundamental solutions of the corresponding fractional Dirac operator is also obtained. Particular cases are considered in both cases.

A Posteriori Error Analysis of hp-Version Discontinuous Galerkin Finite Element Methods for Second-Order Quasilinear Elliptic Problems

We develop the a-posteriori error analysis of hp-version interior-penalty discontinuous Galerkin finite element methods for a class of second-order quasilinear elliptic partial differential equations. Computable upper and lower bounds on the error are derived in terms of a natural (mesh-dependent) energy norm. The bounds are explicit in the local mesh size and the local degree of the approximating polynomial. The performance of the proposed estimators within an automatic hp-adaptive refinement procedure is studied through numerical experiments.

Discontinuous Galerkin Methods for Advection-Diffusion-Reaction Problems on Anisotropically Refined Meshes

In this paper we consider the a posteriori and a priori error analysis of discontinuous Galerkin interior penalty methods for second-order partial differential equations with nonnegative characteristic form on anisotropically refined computational meshes. In particular, we discuss the question of error estimation for linear target functionals, such as the outflow flux and the local average of the solution. Based on our a posteriori error bound we design and implement the corresponding adaptive algorithm to ensure reliable and efficient control of the error in the prescribed functional to within a given tolerance. This involves exploiting both local isotropic and anisotropic mesh refinement. The theoretical results are illustrated by a series of numerical experiments.

Discontinuous Galerkin Methods on hp-Anisotropic Meshes I: A Priori Error Analysis

We consider the a priori error analysis of hp-version interior penalty discontinuous Galerkin methods for second-order partial differential equations with nonnegative characteristic form under weak assumptions on the mesh design and the local finite element spaces employed. In particular, we prove a priori hp-error bounds for linear target functionals of the solution, on (possibly) anisotropic computational meshes with anisotropic tensor-product polynomial basis functions. The theoretical results are illustrated by a numerical experiment.

A bistable reaction-diffusion system in a stretching flow

We examine the evolution of a bistable reaction in a one-dimensional stretching flow, as a model for chaotic advection. We derive two reduced systems of ordinary differential equations (ODEs) for the dynamics of the governing advection-reaction-diffusion partial differential equations (PDE), for pulse-like and for plateau-like solutions, based on a non-perturbative approach. This reduction allows us to study the dynamics in two cases: first, close to a saddle-node bifurcation at which a pair of nontrivial steady states are born as the dimensionless reaction rate (Damkoehler number) is increased, and, second, for large Damkoehler number, far away from the bifurcation. The main aim is to investigate the initial-value problem and to determine when an initial condition subject to chaotic stirring will decay to zero and when it will give rise to a nonzero final state. Comparisons with full PDE simulations show that the reduced pulse model accurately predicts the threshold amplitude for a pulse initial condition to give rise to a nontrivial final steady state, and that the reduced plateau model gives an accurate picture of the dynamics of the system at large Damkoehler number. Published in Physica D (2006)