Neste trabalho, é feita a análise dos modos de plasmon que se propagam em um filme metálico que cobre uma fibra óptica generalizada. Os modos de plasmon estudados são: Fuga pela Cobertura (lcv), Ligado Simétrico (Sb), Fuga pelo Núcleo (lcr) e Ligado Assimétrico (ab). Os filmes metálicos, para efeito de comparação, utilizados neste trabalho, são: a prata, o ouro e o paládio. Desenvolveu-se um modelo matemático do fenômeno eletromagnético e um software, que gerou um banco de dados que facilitasse a análise de estruturas, com diversas combinações de parâmetros. Com o banco de dados, foram obtidos diversos gráficos, que permitiram: analisar os modos de plasmon, verificar a atenuação das ondas e o comportamento do campo eletromagnético em cada região da estrutura. As confrontações entre as estruturas com filmes de: prata, ouro e paládio, permitiram concluir que aquelas elaboradas com os filmes de prata e de ouro são as que apresentam menores perdas, portanto, as recomendadas na confecção de sensores. Como a prata é mais acessível que o ouro, aconselha-se a sua utilização. A análise e os resultados deste trabalho são originais na literatura especializada. / In this work it is done the analysis of the ways of plasmon that developed themselves in a mettalic film that covers an optic and generalized fiber. The studied ways of plasmon are: Cover-Leaky (lcv), Symmetric Bounded (Sb), Core-Leaky (lcr) and Asymmetrical Bounded (ab). The metallic films, for reason of comparison, used in this work, are: the silver, the gold and the paladium. It was developed a mathematical model of the electricmagnetic phenomenon and a software, that created a database which could make it easy the analysis of structures with several combinations of parameters. With the database have been got several graphs that let: analyse the types of plasmon, see the reduction of waves and the behaviour of the electricmagnetic field in each area of the structure. The confrontations between the structures with films of silver, gold and paladium made it possible to conclude that those elaborated with films of silver and gold are those that present smaller losses, so, are recommended in the confection of sensors. As silver is more accecible than gold it is advisable its utilization. The analysis and the results of this work are original in the specialized literature.

Neste trabalho são analisados os quatro modos de plasmon, ligados simétrico (Sb) e assimétrico (ab), fuga pelo núcleo (ln) e fuga pela cobertura (lc), que se propagam em uma fibra óptica fracamente guiada envolta por um filme metálico. No filme metálico é depositada uma camada dielétrica extra e acima desta, uma outra denominada cobertura. A análise será desenvolvida para filmes metálicos de prata, paládio e ouro. Esta estrutura é muito útil na confecção de sensores ópticos. / In this work the four Plasmon modes are analyzed, the symmetrical (Sb) and asymmetrical bounded (ab); the core (ln) and covering leaky modes (lc), that propagate in weakly guided optical fibers with a metallic film around that. In the metallic film a layer extra dielectric is deposited and above this, another layer denominated covering. The analysis will be developed for metallic films of the Silver, Palladium and Gold. This structure is very useful to making optical sensors.

Neste trabalho, é feita a análise dos modos de plasmon que se propagam em um filme metálico que cobre uma fibra óptica generalizada. Os modos de plasmon estudados são: Fuga pela Cobertura (lcv), Ligado Simétrico (Sb), Fuga pelo Núcleo (lcr) e Ligado Assimétrico (ab). Os filmes metálicos, para efeito de comparação, utilizados neste trabalho, são: a prata, o ouro e o paládio. Desenvolveu-se um modelo matemático do fenômeno eletromagnético e um software, que gerou um banco de dados que facilitasse a análise de estruturas, com diversas combinações de parâmetros. Com o banco de dados, foram obtidos diversos gráficos, que permitiram: analisar os modos de plasmon, verificar a atenuação das ondas e o comportamento do campo eletromagnético em cada região da estrutura. As confrontações entre as estruturas com filmes de: prata, ouro e paládio, permitiram concluir que aquelas elaboradas com os filmes de prata e de ouro são as que apresentam menores perdas, portanto, as recomendadas na confecção de sensores. Como a prata é mais acessível que o ouro, aconselha-se a sua utilização. A análise e os resultados deste trabalho são originais na literatura especializada. / In this work it is done the analysis of the ways of plasmon that developed themselves in a mettalic film that covers an optic and generalized fiber. The studied ways of plasmon are: Cover-Leaky (lcv), Symmetric Bounded (Sb), Core-Leaky (lcr) and Asymmetrical Bounded (ab). The metallic films, for reason of comparison, used in this work, are: the silver, the gold and the paladium. It was developed a mathematical model of the electricmagnetic phenomenon and a software, that created a database which could make it easy the analysis of structures with several combinations of parameters. With the database have been got several graphs that let: analyse the types of plasmon, see the reduction of waves and the behaviour of the electricmagnetic field in each area of the structure. The confrontations between the structures with films of silver, gold and paladium made it possible to conclude that those elaborated with films of silver and gold are those that present smaller losses, so, are recommended in the confection of sensors. As silver is more accecible than gold it is advisable its utilization. The analysis and the results of this work are original in the specialized literature.

A evolução nos sistemas digitais de comunicação está intrinsicamente relacionada ao desenvolvimento da tecnologia de fibras ópticas. Desde a sua criação, na década de 60, inúmeras pesquisas vem sendo realizadas com o intuito de aumentar a capacidade de informação transmitida, por meio da redução da atenuação, controle da dispersão cromática e eliminação das não-linearidades. Neste contexto, as Fibras de Bragg surgem como uma estrutura de grande potencialidade para se minimizar tais inconvenientes. As fibras de Bragg possuem um mecanismo de operação diferente em relação às fibras tradicionais de suportar os modos confinados. Nelas, o núcleo possui um baixo índice de refração, e a casca é constituída por anéis dielétricos de diferentes índices de refração, alocados alternadamente. Para uma fibra de Bragg com núcleo oco, como a considerada neste trabalho, há perdas decorrentes dos modos de fuga. Portanto, a análise da dispersão destas estruturas se situa no plano complexo, tornando-a muito difícil. Esta dissertação será fundamentada em uma estratégia imprescindível à análise dos modos transversais TE0m, TM0m e dos híbridos. Os resultados encontrados são validados confrontando-os com os obtidos na literatura. O trabalho discutirá as perdas e dispersões dos modos citados, e os resultados obtidos poderão nortear as pesquisas das fibras de Bragg. / The evolution of digital communication systems is intrinsically related to the development of optical fiber technology. Since its creation in the 1960s, many studies have been conducted in order to increase the system capacity, such as the attenuation reduction, chromatic dispersion control and elimination of nonlinearities. In this context, Bragg fibers appear as a structure with great potential to mitigate these drawbacks. Bragg fibers have a different operational mechanism with respect to traditional fibers to support the confined modes. Their core has a low refractive index, and the cladding consists of dielectric rings of different refractive indices, allocated alternately. For a Bragg fiber with hollow core, as considered in this paper, there are losses due to the occurrence of leaky modes. Therefore, the dispersion analysis of these structures falls in the complex plane, making it even harder. This dissertation will be based on a strategy essential to the analysis of transverse modes: TE0m, TM0m and hybrids. The found results have been validated by comparing them with those obtained in the literature. The paper discusses the losses and dispersions of the mentioned modes, and the results obtained will serve to guide the research on Bragg fibers.

The primary goal of this work is to arrive at direction tomography (DT) algorithms freed from the severe linearization in the formulation, and as-assumptions on variation of the refractive index distribution (RID), involved in the earlier approaches based on Born and Royton approximations and the Fourier di reaction theorem (FDT). To start with, a direct single-step re-covery of RID from intensity measurements is demonstrated, replacing the common two-step procedure involving, rest the recovery of phase from in-density followed by the inversion of scattered led for the RID. The information loss, unavoidable in a two-step procedure is thus successfully addressed. Secondly, an iterative method which works with a forward model obtained directly from the Helmholtz equation is developed. This forward model, though has simplifying assumptions, is more general and can accommodate larger variations in RID than that allowed in the previous linear models. The iterative procedure has an update step which uses a linearization of the forward model and a re-linearization step at the updated RID. The procedure which directly employs the measured intensities is used as part of a deterministic Gauss-Newton algorithm and a stochastic optimization algorithm which uses the ensemble Kalman lter to arrive at the recursive update. The stochastic method is found to be more noise-tolerant and efficient to take care of process model inaccuracies. The proof is seen in better reconstructions from experimental data for two example objects, namely, a graded-index optical bre and a photonic-crystal bre. It is further ob-served that the reconstructions from photonic crystal bre are blurred, noisy and less accurate. Identifying the inaccurate implementation of the modemed Helmholtz equation for large k values employing the current sampling rate as the shortcoming, a new procedure, which splits the bandwidth into smaller components using short-time Fourier Transform is developed. The set of equations arrived at, each t for a narrow frequency band, is solved and the solutions are reassembled to obtain the scattered led for the original problem. The simulated di rated intensities so obtained are better matched to their measured experimental counterparts. However, the impel-mentation of the mode end procedure is computation-intensive, for which a parallel-processing machine can be a good solution. The recovery of RID with this mode cation is not attempted in this work and is left for future implementation.

A evolução nos sistemas digitais de comunicação está intrinsicamente relacionada ao desenvolvimento da tecnologia de fibras ópticas. Desde a sua criação, na década de 60, inúmeras pesquisas vem sendo realizadas com o intuito de aumentar a capacidade de informação transmitida, por meio da redução da atenuação, controle da dispersão cromática e eliminação das não-linearidades. Neste contexto, as Fibras de Bragg surgem como uma estrutura de grande potencialidade para se minimizar tais inconvenientes. As fibras de Bragg possuem um mecanismo de operação diferente em relação às fibras tradicionais de suportar os modos confinados. Nelas, o núcleo possui um baixo índice de refração, e a casca é constituída por anéis dielétricos de diferentes índices de refração, alocados alternadamente. Para uma fibra de Bragg com núcleo oco, como a considerada neste trabalho, há perdas decorrentes dos modos de fuga. Portanto, a análise da dispersão destas estruturas se situa no plano complexo, tornando-a muito difícil. Esta dissertação será fundamentada em uma estratégia imprescindível à análise dos modos transversais TE0m, TM0m e dos híbridos. Os resultados encontrados são validados confrontando-os com os obtidos na literatura. O trabalho discutirá as perdas e dispersões dos modos citados, e os resultados obtidos poderão nortear as pesquisas das fibras de Bragg. / The evolution of digital communication systems is intrinsically related to the development of optical fiber technology. Since its creation in the 1960s, many studies have been conducted in order to increase the system capacity, such as the attenuation reduction, chromatic dispersion control and elimination of nonlinearities. In this context, Bragg fibers appear as a structure with great potential to mitigate these drawbacks. Bragg fibers have a different operational mechanism with respect to traditional fibers to support the confined modes. Their core has a low refractive index, and the cladding consists of dielectric rings of different refractive indices, allocated alternately. For a Bragg fiber with hollow core, as considered in this paper, there are losses due to the occurrence of leaky modes. Therefore, the dispersion analysis of these structures falls in the complex plane, making it even harder. This dissertation will be based on a strategy essential to the analysis of transverse modes: TE0m, TM0m and hybrids. The found results have been validated by comparing them with those obtained in the literature. The paper discusses the losses and dispersions of the mentioned modes, and the results obtained will serve to guide the research on Bragg fibers.

Dans cette thèse on étudie la régularité de problèmes elliptiques (Laplace, Helmholtz) ou paraboliques (équation de la chaleur) avec donnée mesure dans divers cadres géométriques. Ainsi, on considère pour les seconds membres des masses de Dirac en un point, sur une ligne infinie, semi-infinie ou finie, et également sur une courbe régulière. Les solutions de ces problèmes étant singulières sur la fracture (modélisée par la masse de Dirac dans le second membre), on étudie la régularité dans des espaces de Sobolev avec poids. Dans le cas d'une fracture droite, on utilise une technique classique qui consiste à appliquer une transformée de Fourier ou de Mellin à l'équation de Laplace. Ceci nous amène à étudier l'équation de Helmholtz en 2D. Pour ce dernier, on montre des estimations uniformes qui permettent ensuite de prendre la transformée inverse et d'obtenir le résultat de régularité attendu. De même, la transformée de Laplace transforme l'équation de la chaleur dans la même équation de Helmholtz en 2D. Dans le cas d'une fracture courbe régulière, grâce aux résultats de [D'angelo:2012], en utilisant un argument de localisation et un recouvrement dyadique, on obtient une régularité améliorée de la solution toujours dans les espaces de Sobolev avec poids. / In this thesis, we study the regularity of elliptic problems (Laplace, Helmholtz) or parabolic problems (heat equation) with measure data in different geometric frames. Thus, we consider for the second members, Dirac masses at a point, on a line, on a half-line, or on a bounded segment, and also on a regular curve.  As the solutions of these problems are singular on the fracture (modeled by Dirac mass in the second member), we study their regularity in weighted Sobolev spaces.   In the case of a straight fracture, using Fourier or Mellin technique reduces the problem in dimension three to a Helmholtz problem in dimension two. For the latter, we prove uniform estimates, which are then used to apply the inverse transform and to obtain the expected regularity result. Similarly, the Laplace transformation transforms the heat equation into the same Helmholtz equation in 2D.  In the case of a smooth curve fracture, thanks to the results of [D'angelo:2012], using a localization argument and a dyadic recovery we get an improved smoothness of the solution always in weighted Sobolev spaces.

La présente thèse est consacrée à l'étude de quelques problèmes inverses pour l'équation de Helmholtz jauge-covariante, dont des cas particuliers comprennent l'équation de Schrödinger pour une particule élémentaire chargée dans un champ magnétique et l'équation d'onde harmonique en temps qui décrive des ondes acoustiques dans un fluide en écoulement. Ces problèmes ont comme motivation des applications dans des tomographies différentes, qui comprennent la tomographie acoustique, la tomographie qui utilise des particules élémentaires et la tomographie d'impédance électrique. En particulier, nous étudions des problèmes inverses motivés par des applications en tomographie acoustique de fluide en écoulement. Nous proposons des formules et équations qui permettent de réduire le problème de tomographie acoustique à un problème de diffusion inverse approprié. En suivant, nous développons un algorithme fonctionnel-analytique pour la résolution de ce problème de diffusion inverse. Cependant, en général, la solution de ce problème n'est unique qu'à une transformation de jauge appropriée près. À cet égard, nous établissons des formules qui permettent de se débarrasser de cette non-unicité de jauge et retrouver des paramètres du fluide, en mesurant des ondes acoustiques à des plusieurs fréquences. Nous présentons également des exemples des fluides qui ne sont pas distinguable dans le cadre de tomographie acoustique considérée. En suivant, nous considérons le problème de diffusion inverse sans information de phase. Ce problème est motivé par des applications en tomographie qui utilise des particules élémentaires, où seulement le module de l'amplitude de diffusion peut être mesuré facilement. Nous établissons des estimations dans l'espace de configuration pour les reconstructions sans phase de type Borne, qui sont requises pour le développement des méthodes de diffusion inverse précises. Finalement, nous considérons le problème de détermination d'une surface de Riemann dans le plan projectif à partir de son bord. Ce problème survient comme une partie du problème de Dirichlet-Neumann inverse pour l'équation de Laplace sur une surface inconnue, qui est motivé par des applications en tomographie d'impédance électrique. / This work is devoted to study of some inverse problems for the gauge-covariant Helmholtz equation, whose particular cases include the Schrödinger equation for a charged elementary particle in a magnetic field and the time-harmonic wave equation describing sound waves in a moving fluid. These problems are mainly motivated by applications in different tomographies, including acoustic tomography, tomography using elementary particles and electrical impedance tomography. In particular, we study inverse problems motivated by applications in acoustic tomography of moving fluid. We present formulas and equations which allow to reduce the acoustic tomography problem to an appropriate inverse scattering problem. Next, we develop a functional-analytic algorithm for solving this inverse scattering problem. However, in general, the solution to the latter problem is unique only up to an appropriate gauge transformation. In this connection, we give formulas and equations which allow to get rid of this gauge non-uniqueness and recover the fluid parameters, by measuring acoustic fields at several frequencies. We also present examples of fluids which are not distinguishable in this acoustic tomography setting. Next, we consider the inverse scattering problem without phase information. This problem is motivated by applications in tomography using elementary particles, where only the absolute value of the scattering amplitude can be measured relatively easily. We give estimates in the configuration space for the phaseless Born-type reconstructions, which are needed for the further development of precise inverse scattering algorithms. Finally, we consider the problem of determination of a Riemann surface in the complex projective plane from its boundary. This problem arises as a part of the inverse Dirichlet-to-Neumann problem for the Laplace equation on an unknown 2-dimensional surface, and is motivated by applications in electrical impedance tomography.

In this thesis we look at acoustic wave propagation under water. We look in particular at waves generated by a point source and what happens with the propagation when we model the bottom as flat or as curvilinear. We assume the source to be working at a certain frequency and therefore we model this problem by solving the Helmholtz equation. Since Helmholtz equation has some unwanted numerical properties we are interested in finding new numerical methods that could accelerate the solver. In this thesis we use the Waveholtz iteration, which solves Helmholtz equation by connecting it to the time-dependent wave equation. We use finite differences and the SBP-SAT method to approximate the spatial problem numerically and for modelling the sea bottom we use curvilinear coordinates.  To compare the Waveholtz iteration we also solve Helmholtz equation with a naive solver. The naive solver consists of approximating the equation with finite differences and then solving the linear system of equation by some iterative solver, which for our tests will be GMRES. The results show that the Waveholtz iteration converges in less iterations than our naive solver. It also shows that the number of iterations stays unchanged when changing our discretization, which otherwise is a big problem for our naive solver. This allows us to increase the accuracy of our numerical solution without changing the computation time too much.  We show that the number of iterations increases according to theory for an increasing frequency, and that for open problems we even see a smaller increase. For certain resonant frequencies in Helmholtz equation we do not expect the Waveholtz iteration to converge. In the neighbourhood of these frequencies the convergence becomes slow and we need many iterations for a solution of a certain accuracy. By reformulating the Waveholtz iteration as a Krylov solution we can see that resonances in Helmholtz equation have a smaller impact of the convergence. / I detta examensarbete undersöker vi akustisk vågutbredning i vatten. Vi kollar specifikt på vågor som genereras av en punktkälla och vad som sker när vi modellerar botten som plan eller som kurvlinjär. Då vi antar att punktkällan arbetar vid en bestämd frekvens, kommer vi modellera det fysikaliska problemet genom att lösa Helmholtz ekvation. Helmholtz ekvation har dock några numeriska egenskaper som är oönskade, och därför finns ett intresse av att hitta nya numeriska metoder som löser ekvationen. I detta examensarbete undersöker vi Waveholtz iteration, som löser Helmholtz ekvation genom att koppla den till den tidsberoende vågekvationen. Vi använder finita differenser och SBP-SAT metoden för att approximera det rumsliga problemet numeriskt. För att ge en detaljerad beskrivning av botten använder vi kurvlinjära koordinater. För att jämföra Waveholtz iterationen med något löser vi även Helmholtz med hjälp av en naiv lösare. Den naiva lösaren består av att approximera problemet med finita differenser och sedan lösa det linjära systemet rakt av med en iterativ lösare (vilket för våra fall kommer vara GMRES). Resultatet visar att Waveholtz iteration konvergerar på ett lägre antal iterationer än vår naiva lösare. Det visar även att antalet iterationer inte förändras när vi ändrar diskretisering, vilket annars är ett problem för vår naiva lösare. Detta innebär att vi kan få en högre noggrannhet utan att förlänga beräkningstiden alltför mycket.  Vi visar även att antalet iterationer ökar som förväntat med en ökad frekvens, samt att för öppna problem så ökar antalet iteration mindre än enligt teorin. Vid vissa resonanta frekvenser i Helmholtz ekvation förväntar vi oss att Waveholtz iteration inte kommer konvergerar. I närheten av dessa frekvenser blir konvergensen långsam och vi behöver många iterationer för att lösa problemet. Genom att formulera Waveholtz iteration som en Krylov lösning kommer resonanser i Helmholtz ekvation ge en mindre negativ effekt på konvergensen än om den är formulerad som en fixpunkts iteration.

We present three new stabilized finite element (FE) based Petrov-Galerkin methods for the convection-diffusionreaction (CDR), the Helmholtz and the Stokes problems, respectively. The work embarks upon a priori analysis of a consistency recovery procedure for some stabilization methods belonging to the Petrov- Galerkin framework. It was ound that the use of some standard practices (e.g. M-Matrices theory) for the design of essentially non-oscillatory numerical methods is not appropriate when consistency recovery methods are employed. Hence, with respect to convective stabilization, such recovery methods are not preferred. Next, we present the design of a high-resolution Petrov-Galerkin (HRPG) method for the CDR problem. The structure of the method in 1 D is identical to the consistent approximate upwind (CAU) Petrov-Galerkin method [doi: 10.1016/0045-7825(88)90108-9] except for the definitions of he stabilization parameters. Such a structure may also be attained via the Finite Calculus (FIC) procedure [doi: 10.1 016/S0045-7825(97)00119-9] by an appropriate definition of the characteristic length. The prefix high-resolution is used here in the sense popularized by Harten, i.e. second order accuracy for smooth/regular regimes and good shock-capturing in non-regular re9jmes. The design procedure in 1 D embarks on the problem of circumventing the Gibbs phenomenon observed in L projections. Next, we study the conditions on the stabilization parameters to ircumvent the global oscillations due to the convective term. A conjuncture of the two results is made to deal with the problem at hand that is usually plagued by Gibbs, global and dispersive oscillations in the numerical solution. A multi dimensional extension of the HRPG method using multi-linear block finite elements is also presented. Next, we propose a higher-order compact scheme (involving two parameters) on structured meshes for the Helmholtz equation. Making the parameters equal, we recover the alpha-interpolation of the Galerkin finite element method (FEM) and the classical central finite difference method. In 1 D this scheme is identical to the alpha-interpolation method [doi: 10.1 016/0771 -050X(82)90002-X] and in 2D choosing the value 0.5 for both the parameters, we recover he generalized fourth-order compact Pade approximation [doi: 10.1 006/jcph.1995.1134, doi: 10.1016/S0045- 7825(98)00023-1] (therein using the parameter V = 2). We follow [doi: 10.1 016/0045-7825(95)00890-X] for the analysis of this scheme and its performance on square meshes is compared with that of the quasi-stabilized FEM [doi: 10.1016/0045-7825(95)00890-X]. Generic expressions for the parameters are given that guarantees a dispersion accuracy of sixth-order should the parameters be distinct and fourth-order should they be equal. In the later case, an expression for the parameter is given that minimizes the maximum relative phase error in 2D. A Petrov-Galerkin ormulation that yields the aforesaid scheme on structured meshes is also presented. Convergence studies of the error in the L2 norm, the H1 semi-norm and the I ~ Euclidean norm is done and the pollution effect is found to be small. / Presentamos tres nuevos metodos estabilizados de tipo Petrov- Galerkin basado en elementos finitos (FE) para los problemas de convecci6n-difusi6n- reacci6n (CDR), de Helmholtz y de Stokes, respectivamente. El trabajo comienza con un analisis a priori de un metodo de recuperaci6n de la consistencia de algunos metodos de estabilizaci6n que pertenecen al marco de Petrov-Galerkin. Hallamos que el uso de algunas de las practicas estandar (por ejemplo, la eoria de Matriz-M) para el diserio de metodos numericos esencialmente no oscilatorios no es apropiado cuando utilizamos los metodos de recu eraci6n de la consistencia. Por 10 tanto, con res ecto a la estabilizaci6n de conveccion, no preferimos tales metodos de recuperacion . A continuacion, presentamos el diser'io de un metodo de Petrov-Galerkin de alta-resolucion (HRPG) para el problema CDR. La estructura del metodo en 10 es identico al metodo CAU [doi: 10.1016/0045-7825(88)90108-9] excepto en la definicion de los parametros de estabilizacion. Esta estructura tambien se puede obtener a traves de la formulacion del calculo finito (FIC) [doi: 10.1 016/S0045- 7825(97)00119-9] usando una definicion adecuada de la longitud caracteristica. El prefijo de "alta-resolucion" se utiliza aqui en el sentido popularizado por Harten, es decir, tener una solucion con una precision de segundo orden en los regimenes suaves y ser esencialmente no oscilatoria en los regimenes no regulares. El diser'io en 10 se embarca en el problema de eludir el fenomeno de Gibbs observado en las proyecciones de tipo L2. A continuacion, estudiamos las condiciones de los parametros de estabilizacion para evitar las oscilaciones globales debido al ermino convectivo. Combinamos los dos resultados (una conjetura) para tratar el problema COR, cuya solucion numerica sufre de oscilaciones numericas del tipo global, Gibbs y dispersiva. Tambien presentamos una extension multidimensional del metodo HRPG utilizando los elementos finitos multi-lineales. fa. continuacion, proponemos un esquema compacto de orden superior (que incluye dos parametros) en mallas estructuradas para la ecuacion de Helmholtz. Haciendo igual ambos parametros, se recupera la interpolacion lineal del metodo de elementos finitos (FEM) de tipo Galerkin y el clasico metodo de diferencias finitas centradas. En 10 este esquema es identico al metodo AIM [doi: 10.1 016/0771 -050X(82)90002-X] y en 20 eligiendo el valor de 0,5 para ambos parametros, se recupera el esquema compacto de cuarto orden de Pade generalizada en [doi: 10.1 006/jcph.1 995.1134, doi: 10.1 016/S0045-7825(98)00023-1] (con el parametro V = 2). Seguimos [doi: 10.1 016/0045-7825(95)00890-X] para el analisis de este esquema y comparamos su rendimiento en las mallas uniformes con el de "FEM cuasi-estabilizado" (QSFEM) [doi: 10.1016/0045-7825 (95) 00890-X]. Presentamos expresiones genericas de los para metros que garantiza una precision dispersiva de sexto orden si ambos parametros son distintos y de cuarto orden en caso de ser iguales. En este ultimo caso, presentamos la expresion del parametro que minimiza el error maxima de fase relativa en 20. Tambien proponemos una formulacion de tipo Petrov-Galerkin ~ue recupera los esquemas antes mencionados en mallas estructuradas. Presentamos estudios de convergencia del error en la norma de tipo L2, la semi-norma de tipo H1 y la norma Euclidiana tipo I~ y mostramos que la perdida de estabilidad del operador de Helmholtz ("pollution effect") es incluso pequer'ia para grandes numeros de onda. Por ultimo, presentamos una coleccion de metodos FE estabilizado para el problema de Stokes desarrollados a raves del metodo FIC de primer orden y de segundo orden. Mostramos que varios metodos FE de estabilizacion existentes y conocidos como el metodo de penalizacion, el metodo de Galerkin de minimos cuadrados (GLS) [doi: 10.1016/0045-7825(86)90025-3], el metodo PGP (estabilizado a traves de la proyeccion del gradiente de presion) [doi: 10.1 016/S0045-7825(96)01154-1] Y el metodo OSS (estabilizado a traves de las sub-escalas ortogonales) [doi: 10.1016/S0045-7825(00)00254-1] se recuperan del marco general de FIC. Oesarrollamos una nueva familia de metodos FE, en adelante denominado como PLS (estabilizado a traves del Laplaciano de presion) con las formas no lineales y consistentes de los parametros de estabilizacion. Una caracteristica distintiva de la familia de los metodos PLS es que son no lineales y basados en el residuo, es decir, los terminos de estabilizacion dependera de los residuos discretos del momento y/o las ecuaciones de incompresibilidad. Oiscutimos las ventajas y desventajas de estas tecnicas de estabilizaci6n y presentamos varios ejemplos de aplicacion