A Teoria das equações integrais, desde a segunda metade do século XX, tem assumido um papel cada vez maior no âmbito de problemas aplicados. Com isso, surge a necessidade do desenvolvimento de métodos numéricos cada vez mais eficazes para a resolução deste tipo de equação. Isso tem como consequência a possibilidade de resolução de uma gama cada vez maior de problemas. Nesse sentido, outros tipos de equações integrais estão sendo objeto de estudos, dentre elas as chamadas equações integro-diferenciais. O presente trabalho tem como objetivo o estudo das equações integro-diferenciais singulares lineares e não-lineares. Mais especificamente, no caso linear, apresentamos os principais resultados necessários para a obtenção de um método numérico e a formulação de suas propriedades de convergência. O caso não-linear é apresentado através de um modelo matemático para tubulações em um tipo específico de reator nuclear (LMFBR) no qual origina-se a equação integro-diferencial. A partir da equação integro-diferencial um modelo numérico é proposto com base nas condições físicas do problema / The theory of the integral equations, since the second half of the 20th century, has been assuming an ever more important role in the modelling of applied problems. Consequently, the development of new numerical methods for integral equations is called for and a larger range of problems has been possible to be solved by these new techniques. In this sense, many types of integral equations have been derived from applications and been the object of studies, among them the so called singular integro-differential equation. The present work has, as its main objective, the study of singular integrodifferential equations, both linear and non-linear. More specifically, in the linear case, we present our main results regarding the derivation of a numerical method and its uniform convergence properties. The non-linear case is introduced through the mathematical model of boiler tubes in a specific type of nuclear reactor (LMFBR) from which the integro-differential equation originates. For this integro-differential equation a numerical method is proposed based on the physical conditions of the problem

Submitted by Erika Demachki (erikademachki@gmail.com) on 2014-09-24T20:35:00Z No. of bitstreams: 2 Santos, Claudia Gomes de Oliveira - Dissertação - 2013.pdf: 1950939 bytes, checksum: 050c57553672656134c6b1264cb562a6 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Approved for entry into archive by Jaqueline Silva (jtas29@gmail.com) on 2014-09-24T20:42:50Z (GMT) No. of bitstreams: 2 Santos, Claudia Gomes de Oliveira - Dissertação - 2013.pdf: 1950939 bytes, checksum: 050c57553672656134c6b1264cb562a6 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Made available in DSpace on 2014-09-24T20:42:50Z (GMT). No. of bitstreams: 2 Santos, Claudia Gomes de Oliveira - Dissertação - 2013.pdf: 1950939 bytes, checksum: 050c57553672656134c6b1264cb562a6 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Previous issue date: 2013-07-31 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / The Boundary Element Method (BEM) has been successfully employed in the analysis of various engineering problems. The BEM consists in a mathematical modeling, for a numerical solution of a system of integral equations, and in their cores may appear singularities. This paper presents the Classical and Hypersingular formulation of the Boundary Element Method for dimensional elastostatic problems with smooth boundary geometry. The improper integrals arising from the singularities of the core in the hypersingular formulation are treated by Hadamard finite parts. In the discretization process two types of interpolation are used, one traditional and the other special. Traditional interpolation is used in all bondary elements that have no point , special interpolation ensures the continuity of the tangential derivative of displacements on the element that contains the point . To accomplish this, a theoretical mathematics study of related topics was performed. The hypersingular formulation developed in this work was implemented through the Intel Visual Fortran compiler. Some problems were analyzed and the obtained results were compared with those of analytical solution or through the Finite Element Method. The results achieved were satisfactory validating the proposed formulation / O Método dos Elementos de Contorno (MEC) vem sendo empregado com sucesso na análise de diversos problemas de engenharia. O MEC consisti em uma modelagem matemática, para resolução numérica de um sistema de equações integrais, e que em seus núcleos podem aparecer singularidades. Nesse trabalho apresenta a formulação Clássica e Hipersingular do Método dos Elementos de Contorno para problemas de elastostática bidimensional com geometria de contornos não suaves. As integrais impróprias que surgem da singularidade do núcleo na formulação hipersingular são tratados por partes finitas de Hadamard. No processo de discretização utiliza-se de dois tipos de interpolação, uma tradicional e outra especial. A interpolação tradicional é utilizada em todos os elementos de contorno que não tem o ponto , a interpolação especial garante a continuidade da derivada tangencial dos deslocamentos no elemento que contém o ponto . Para a realização deste, foi realizado um estudo teórico-matemático dos tópicos afins. Implementou-se a formulação hipersingular desenvolvidas no trabalho através do compilador Intel Visual FORTRAN. Foram analisados alguns problemas e os resultados obtidos comparados àqueles de solução analítica ou através do Método dos Elementos Finitos. Os resultados alcançados mostraram-se satisfatórios validando a formulação proposta.

In this paper, the complete double layer boundary integral equation formulation for Stokes flows is extended to viscoelastic fluids to solve the mobility problem for a system of particles, where the non-linearity is handled by particular solutions of the Stokes inhomogeneous equation. Some techniques of the meshless method are employed and a point-wise solver is used to solve the viscoelastic constitutive equation. Hence volume meshing is avoided. The method is tested against the numerical solution for a sphere settling in the Odroyd-B fluid and some results on a prolate motion in shear flow of the Oldroyd-B fluid are reported and compared with some theoretical and experimental results. / Singapore-MIT Alliance (SMA)

With increasing complexity of microwave integrated circuits and tendency towards building integrated modules, real estate in printed circuit boards becomes more at premium. On the other hand, building MIC's on a single semiconductor substrate such as GaAs has other drawbacks as substrate requirements for different components are sometimes contradictory. This has motivated researchers to consider multi-layer and stacked designs. Multi-layer planar circuits offer advantages that cannot be equaled by traditional single layer designs. In this respect, a new class of planar structures, based upon a multi-layered stack of dual-mode stripline or microstrip patches is becoming increasingly popular. In the new stacked design coupling between planar circuits separated by a ground plane is accomplished through coupling apertures in the common ground plane. <br ><br /> This thesis is about developing a new approximate multiport network model for fast analysis of multi-layered planar structures with ground plane slots. To extend applicability of multiport network model (MNM) to the class of planar structures containing ground plane slots, a generalized network formulation for aperture problems is combined with traditional MNM to account for the presence of the slot. To this end, the slot is replaced by an unknown equivalent surface magnetic current. Slot ports are defined in terms of electric and magnetic fields over the slot in accordance with the generalized network formulation for aperture problems. While traditional MNM for planar circuits is based on generalized impedance matrices, we adopt a hybrid matrix approach for multi-layer structures. The hybrid matrix consists of four sub-matrices that relate terminal voltages and currents of edge and slot ports. The same generalized impedance matrix in the absence of the slot can be used to relate terminal voltages and currents of edge ports when the slot ports are short-circuited. Open circuit voltage at edge ports due to terminal voltages at slot ports and terminal currents at slot ports due to input currents at edge ports are represented by two transfer matrices. Both these transfer matrices can be calculated from 2D analysis which only considers <em>TM^z</em> modes. <br ><br /> Interaction among slot ports, represented by a generalized admittance matrix, however, requires considering both <em>TM^z</em> and <em>TE^z</em> modes. This generalized admittance matrix is obtained from tangential component of the magnetic field over the slot due to the equivalent surface magnetic current and relates terminal voltages and currents of slot ports. Full modal expansion consisting of both <em>TM^z</em> and <em>TE^z</em> modes is used to compute the generalized admittance matrix of a slot in a regularly shaped planar cavity. For irregularly shaped patches, modal expansion is not available. Instead, a new contour integral equation for magnetic field, derived for the first time in this thesis, is combined with complex images method for calculation of generalized admittance matrix of a slot radiating in a planar cavity of arbitrary shape. <br ><br /> Once the hybrid matrix representation of a planar circuit on a ground plane containing a slot is derived, it can be connected to the hybrid matrix of any other planar circuit on the other side of the ground plane. This can be done by enforcing network equivalent of continuity of tangential fields across the slot. This leads to a generalized impedance matrix for the multi-layer structure relating terminal voltages and currents of edge ports of both planar circuits. <br ><br /> To show the accuracy of the proposed method of analysis, several proof-of-concept structures have been analyzed by both this method and ANSOFT HFSS full-wave simulator as a reference. In most cases excellent agreement is achieved in predicting the return loss and radiation patterns of these multi-layer structures which proves the validity of the proposed approach for fast analysis and design of multi-layer planar structures.

With increasing complexity of microwave integrated circuits and tendency towards building integrated modules, real estate in printed circuit boards becomes more at premium. On the other hand, building MIC's on a single semiconductor substrate such as GaAs has other drawbacks as substrate requirements for different components are sometimes contradictory. This has motivated researchers to consider multi-layer and stacked designs. Multi-layer planar circuits offer advantages that cannot be equaled by traditional single layer designs. In this respect, a new class of planar structures, based upon a multi-layered stack of dual-mode stripline or microstrip patches is becoming increasingly popular. In the new stacked design coupling between planar circuits separated by a ground plane is accomplished through coupling apertures in the common ground plane. <br ><br /> This thesis is about developing a new approximate multiport network model for fast analysis of multi-layered planar structures with ground plane slots. To extend applicability of multiport network model (MNM) to the class of planar structures containing ground plane slots, a generalized network formulation for aperture problems is combined with traditional MNM to account for the presence of the slot. To this end, the slot is replaced by an unknown equivalent surface magnetic current. Slot ports are defined in terms of electric and magnetic fields over the slot in accordance with the generalized network formulation for aperture problems. While traditional MNM for planar circuits is based on generalized impedance matrices, we adopt a hybrid matrix approach for multi-layer structures. The hybrid matrix consists of four sub-matrices that relate terminal voltages and currents of edge and slot ports. The same generalized impedance matrix in the absence of the slot can be used to relate terminal voltages and currents of edge ports when the slot ports are short-circuited. Open circuit voltage at edge ports due to terminal voltages at slot ports and terminal currents at slot ports due to input currents at edge ports are represented by two transfer matrices. Both these transfer matrices can be calculated from 2D analysis which only considers <em>TM^z</em> modes. <br ><br /> Interaction among slot ports, represented by a generalized admittance matrix, however, requires considering both <em>TM^z</em> and <em>TE^z</em> modes. This generalized admittance matrix is obtained from tangential component of the magnetic field over the slot due to the equivalent surface magnetic current and relates terminal voltages and currents of slot ports. Full modal expansion consisting of both <em>TM^z</em> and <em>TE^z</em> modes is used to compute the generalized admittance matrix of a slot in a regularly shaped planar cavity. For irregularly shaped patches, modal expansion is not available. Instead, a new contour integral equation for magnetic field, derived for the first time in this thesis, is combined with complex images method for calculation of generalized admittance matrix of a slot radiating in a planar cavity of arbitrary shape. <br ><br /> Once the hybrid matrix representation of a planar circuit on a ground plane containing a slot is derived, it can be connected to the hybrid matrix of any other planar circuit on the other side of the ground plane. This can be done by enforcing network equivalent of continuity of tangential fields across the slot. This leads to a generalized impedance matrix for the multi-layer structure relating terminal voltages and currents of edge ports of both planar circuits. <br ><br /> To show the accuracy of the proposed method of analysis, several proof-of-concept structures have been analyzed by both this method and ANSOFT HFSS full-wave simulator as a reference. In most cases excellent agreement is achieved in predicting the return loss and radiation patterns of these multi-layer structures which proves the validity of the proposed approach for fast analysis and design of multi-layer planar structures.

Computer simulation of wave propagation is an active research area as wave phenomena are prevalent in many applications. Examples include wireless communication, radar cross section, underwater acoustics, and seismology. For high frequency waves, this is a challenging multiscale problem, where the small scale is given by the wavelength while the large scale corresponds to the overall size of the computational domain. Research into wave equation modeling can be divided into two regimes: time domain and frequency domain. In each regime, there are two further popular research directions for the numerical simulation of the scattered wave. One relies on direct discretization of the wave equation as a hyperbolic partial differential equation in the full physical domain. The other direction aims at solving an equivalent integral equation on the surface of the scatterer. In this dissertation, we present three new techniques for the frequency domain, boundary integral equations. / text

<p>This dissertation focuses on a FFT based forward EM solver and its application in inverse problems. The main contributions of this work are two folded. On the one hand, it presents the first scaled lab experiment system in the oil and gas industry for through casing hydraulic fracture evaluation. This system is established to validate the feasibility of contrasts enhanced fractures evaluation. On the other hand, this work proposes a FFT based VIE solver for hydraulic fracture evaluation. This efficient solver is needed for numerical analysis of such problem. The solver is then generalized to accommodate scattering simulations for anisotropic inhomogeneous magnetodielectric objects. The inverse problem on anisotropic objects are also studied.</p><p>Before going into details of specific applications, some background knowledge is presented. This dissertation starts with an introduction to inverse problems. Then algorithms for forward and inverse problems are discussed. The discussion on forward problem focuses on the VIE formulation and a frequency domain solver. Discussion on inverse problems focuses on iterative methods.</p><p>The rest of the dissertation is organized by the two categories of inverse problems, namely the inverse source problem and the inverse scattering problem. </p><p>The inverse source problem is studied via an application in microelectronics. In this application, a FFT based inverse source solver is applied to process near field data obtained by near field scanners. Examples show that, with the help of this inverse source solver, the resolution of unknown current source images on a device under test is greatly improved. Due to the improvement in resolution, more flexibility is given to the near field scan system.</p><p>Both the forward and inverse solver for inverse scattering problems are studied in detail. As a forward solver for inverse scattering problems, a fast FFT based method for solving VIE of magnetodielectric objects with large electromagnetic contrasts are presented due to the increasing interest in contrasts enhanced full wave EM imaging. This newly developed VIE solver assigns different basis functions of different orders to expand flux densities and vector potentials. Thus, it is called the mixed ordered BCGS-FFT method. The mixed order BCGS-FFT method maintains benefits of high order basis functions for VIE while keeping correct boundary conditions for flux densities and vector potentials. Examples show that this method has an excellent performance on both isotropic and anisotropic objects with high contrasts. Examples also verify that this method is valid in both high and low frequencies. Based on the mixed order BCGS-FFT method, an inverse scattering solver for anisotropic objects is studied. The inverse solver is formulated and solved by the variational born iterative method. An example given in this section shows a successful inversion on an anisotropic magnetodielectric object. </p><p>Finally, a lab scale hydraulic fractures evaluation system for oil/gas reservoir based on previous discussed inverse solver is presented. This system has been setup to verify the numerical results obtained from previously described inverse solvers. These scaled experiments verify the accuracy of the forward solver as well as the performance of the inverse solver. Examples show that the inverse scattering model is able to evaluate contrasts enhanced hydraulic fractures in a shale formation. Furthermore, this system, for the first time in the oil and gas industry, verifies that hydraulic fractures can be imaged through a metallic casing.</p> / Dissertation

Neste trabalho investigamos, dos pontos de vistas analítico e numérico, o comportamento assintótico da solução da equação do aerofólio, com uma singularidade do tipo Cauchy, de nida sobre um intervalo com uma pequena abertura. Exibimos um modelo matemático com uma solução f" para o intervalo disjunto G" = (−1,−ε) ∪ (ε, 1) e uma solução f0 que corresponde ao limite de f" quando (ε → 0), relacionando esta última com a solução da equação do aerofólio f no intervalo (−1, 1). Além do mais, demonstramos casos particulares de funções ψ = Tm e ψ = Un(onde Tm e Un são os polinômios de Tchebychev do primeiro e segundo tipo respectivamente) em que temos a igualdade f = f0 e conseqüentemente f" ≈ f. Apresentamos e comparamos numericamente as soluções f", f0 e f para diferentes funções ψ e valores de ε no intervalo G". Mostramos ainda soluções quase polinomiais analíticas da equação do aerofólio, e propomos um método espectral para a equação do aerofólio generalizada. Por m, obtemos soluções analíticas das equações do aerofólio para os intervalos G", (−1, 1)\ {0} e (−1, 1) para diferentes funções ψ(t) através da expansão em série da densidade da integral singular com núcleo Cauchy. / In this work we investigate, of the analytical and numerical points of views, the asymptotic behavior of the airfoil equation solution with a singularity of the Cauchy type, de ned over a interval with a small opening. We display a mathematical model with a f" solution to the disjoint interval G" = (−1,−ε)∪(ε, 1) and a f0 solution corresponding to limit of f" when (ε → 0), linking the latter with the solution of the airfoil equation f in the interval (−1, 1). Furthermore, we demonstrate particular cases of functions ψ = Tm and ψ = Un (where Tm and Un are the Chebyshev polynomials of the rst and second type respectively) where we have equality f = f0 and then f" ≈ f. We present and compare numerically the solutions f", f0 and f for di erent functions ψ and values of ε in G". We also show almost polynomial analytical solutions for the airfoil equation, and we propose a spectral method for the generalized airfoil equation. Finally, we obtain analytical solutions of the airfoil equations to the interval G", (−1, 1)\ {0} and (−1, 1) for various functions ψ(t) by expanding in series the density of the Cauchy singular integral.

Made available in DSpace on 2015-05-08T14:59:51Z (GMT). No. of bitstreams: 1 arquivototal.pdf: 2509612 bytes, checksum: f63968d7b069e4d6a16ca4bafd49aa7b (MD5) Previous issue date: 2013-07-26 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / The laminated composite beams are structural components which have aroused great interest in mechanical industry due to providing components of high strength-to-weight ratio, which make them particularly suitable to the automotive and aerospace industry. In recent decades the solutions of mathematical models for laminated beam have been established in a analytical or numerical forms. The latter have been built based on finite element method (FEM) philosophy. In this work a new solution based on Boundary Element Method (BEM) for both classical and for first order laminated beam theory for static loading is established. In addition, the BEM solution is extended to deal with harmonic loading under classic beams theory hypothesis. Discussions on mathematical steps to write down both integral equations and fundamental solutions for laminated beam problem are properly made. From the formulations here proposed, codes are implemented in C++, providing displacements and efforts at the boundary and domain of the beam. In addition, numerical results for typical cases are presented as well. / As vigas de compósitos laminados são componentes estruturais que têm despertado grande interesse na indústria mecânica, principalmente por prover componentes de alta razão resistência-peso, o que as tornam particularmente aplicável à indústria automobilística e aeronáutica. Nas últimas décadas as soluções dos modelos matemáticos de vigas laminadas têm sido estabelecidas na forma analítica ou numérica. Para o último caso, as soluções têm sido construídas na filosofia do método dos elementos finitos (MEF). Nesta dissertação uma nova solução baseada no Método dos Elementos de Contorno (MEC) para as teorias de vigas laminada clássica e de primeira ordem são estabelecidas para carregamentos estáticos. Além disso, a solução para o MEC é estendida para lidar com carregamento dinâmico harmônico sob a hipótese da Teoria Clássica de viga. Nas discussões sobre as etapas matemáticas são descritas as equações integrais e as soluções fundamentais para os problemas de vigas laminadas. A partir das formulações aqui propostas, códigos são implementados na linguagem C++, capaz de fornecer deslocamentos e esforços no contorno e no domínio da viga. Assim como, as frequências naturais. Além disso, os resultados numéricos, comparados com as soluções analíticas disponíveis na literatura, se mostraram de boa qualidade.

Neste trabalho investigamos, dos pontos de vistas analítico e numérico, o comportamento assintótico da solução da equação do aerofólio, com uma singularidade do tipo Cauchy, de nida sobre um intervalo com uma pequena abertura. Exibimos um modelo matemático com uma solução f" para o intervalo disjunto G" = (−1,−ε) ∪ (ε, 1) e uma solução f0 que corresponde ao limite de f" quando (ε → 0), relacionando esta última com a solução da equação do aerofólio f no intervalo (−1, 1). Além do mais, demonstramos casos particulares de funções ψ = Tm e ψ = Un(onde Tm e Un são os polinômios de Tchebychev do primeiro e segundo tipo respectivamente) em que temos a igualdade f = f0 e conseqüentemente f" ≈ f. Apresentamos e comparamos numericamente as soluções f", f0 e f para diferentes funções ψ e valores de ε no intervalo G". Mostramos ainda soluções quase polinomiais analíticas da equação do aerofólio, e propomos um método espectral para a equação do aerofólio generalizada. Por m, obtemos soluções analíticas das equações do aerofólio para os intervalos G", (−1, 1)\ {0} e (−1, 1) para diferentes funções ψ(t) através da expansão em série da densidade da integral singular com núcleo Cauchy. / In this work we investigate, of the analytical and numerical points of views, the asymptotic behavior of the airfoil equation solution with a singularity of the Cauchy type, de ned over a interval with a small opening. We display a mathematical model with a f" solution to the disjoint interval G" = (−1,−ε)∪(ε, 1) and a f0 solution corresponding to limit of f" when (ε → 0), linking the latter with the solution of the airfoil equation f in the interval (−1, 1). Furthermore, we demonstrate particular cases of functions ψ = Tm and ψ = Un (where Tm and Un are the Chebyshev polynomials of the rst and second type respectively) where we have equality f = f0 and then f" ≈ f. We present and compare numerically the solutions f", f0 and f for di erent functions ψ and values of ε in G". We also show almost polynomial analytical solutions for the airfoil equation, and we propose a spectral method for the generalized airfoil equation. Finally, we obtain analytical solutions of the airfoil equations to the interval G", (−1, 1)\ {0} and (−1, 1) for various functions ψ(t) by expanding in series the density of the Cauchy singular integral.