The goal of this thesis is to study of the evolution of 3D electrically charged liquid droplets of fluid evolving under the influence of surface tension and electrostatic forces. In the first part of the thesis, an appropriate mathematical model of the problem is introduced and the linear stability analysis is developed by perturbing a sphere with spherical harmonics. In the second part, the numerical solution of the problem is described with the use of the boundary elements method (BEM) on an adaptive mesh of triangular elements. The numerical method is validated by comparison with exact solutions. Finally, various numerical results are presented. These include neck formation in droplets, the evolution of surfaces with holes, singularity formation on droplets with various symmetries and numerical evidence that oblate spheroids are unstable.

[pt] Apresentam-se as formulações, consolidando a nomenclatura e os principais conceitos dos métodos de elementos de contorno: convencional (MCCEC), híbrido de tensões (MHTEC), híbrido de deslocamentos (MHDEC) e híbrido simplificado de tensões (MHSTEC). proposto o método híbrido simplificado de deslocamentos (MHSDEC), em contrapartida ao MHSTEC, baseando-se nas mesmas hipóteses de aproximação de tensões e deslocamentos do MHDEC e supondo que a solução fundamental em termos de tensões seja válida no contorno. Como decorrência do MHSTEC e do MHSDEC, é apresentado também o método híbrido de malha reduzida dos elementos de contorno (MHMREC), com aplicação computacionalmente vantajosa a problemas no domínio da freqüência ou envolvendo materiais não-homogêneos. A partir da investigação das equações matriciais desses métodos, são identificadas quatro novas relações matriciais, das quais uma verifica-se como válida para a obtenção dos elementos das matrizes de flexibilidade e de deslocamento que não podem ser determinados por integração ou avaliação direta. Também é proposta a correta consideração, ainda não muito bem explicada na literatura, de que forças de superfície devem ser interpoladas em função de atributos de superfície e não de atributos nodais. São apresentadas aplicações numéricas para problemas de potencial para cada método mencionado, em que é verificada a validade das novas relações matriciais. / [en] A consolidated, unified formulation of the conventional (CCBEM), hybrid stress (HSBEM), hybrid displacement (HDBEM) and simplified hybrid stress (SHSBEM) boundary element methods is presented. As a counterpart of SHSBEM, the simplified hybrid displacement boundary element method (SHDBEM) is proposed on the basis of the same stress and displacement approximation hypotheses of the HDBEM and on the assumption that stress fundamental solutions are also valid on the boundary. A combination of the SHSBEM and the SHDBEM gives rise to a provisorily called mesh-reduced hybrid boundary element method (MRHBEM), which seems computationally advantageous when applied to frequency domain problems or non-homogeneous materials. Four new matrix relations are identified, one of which may be used to obtain the flexibility and displacement matrix coefficients that cannot be determined by integration or direct evaluation. It is also proposed the correct consideration, still not well explained in the technical literature, that traction forces should be interpolated as functions of surface and not of nodal attributes. Numerical examples of potential problems are presented for each method, in which the validity of the new matrix relations is verified.

[pt] O desenvolvimento consistente do método convencional dos elementos de contorno (CBEM), com a adição de conceitos da versão simplificada do método híbrido dos elementos de contorno (HBEM), proveniente do potencial variacional de Hellinger-Reissner, conduz-se a um processo computacionalmente mais econômico, sem a necessidade de ter sua precisão numérica reduzida para problemas de grande escala, podendo ser bidimensional ou tridimensional, de potencial ou elasticidade. Conseguiu-se mostrar que as matrizes de potencial duplo e simples do CBEM, H e G, respectivamente, cuja avaliação numérica requer a manipulação de integrais singulares e impróprias, podem ser obtidas de maneira expedita, eliminando-se quase toda a integração numérica, com exceção de algumas integrais regulares. Uma importante característica da formulação proposta, que advém da base variacional do HBEM, é a facilidade da obtenção de resultados em pontos internos, de maneira direta e sem a utilização de qualquer integral de contorno, já que a solução fundamental é a própria solução do problema. O presente trabalho pertence a um projeto cujo resultado final deve ser um código computacional para problemas de grande escala (milhões de graus de liberdade). Nesta fase, alguns exemplos numéricos foram testados para avaliar a aplicabilidade do método expedito, o seu esforço computacional e a convergência do resultado para as variáveis envolvidas no método. Para isso, foram implementados algoritmos para problemas bidimensionais de potencial e elasticidade - usando elementos lineares, quadráticos e cúbicos - e tridimensionais - usando elementos triangulares e quadrilaterais, lineares e quadráticos nos dois casos. Os códigos computacionais foram implementados focando na solução de problemas de grande escala. Espera-se que numa etapa final o projeto possa ser bem mais eficaz, com a incorporação de procedimentos do método fast multipole. / [en] The consistent development of the conventional boundary elements method (CBEM) by adding the concepts of the hybrid boundary element simplified method (HBEM) , from the Hellinger-Reissner variational potential leads to a computationally less intensive procedure, although not necessarily less accurate for large scale, two-dimensional or three-dimensional problems of potential and elasticity. It was shown that both single-layer and double-layer potential matrices, G and H, respectively, are obtained in an expeditious way that vanish almost any numerical integration, except for a few regular integrals, even G and H evaluation requires the handling of singular and improper integrals. The proposed formulation comes from the HBEM variational base and its evaluation at internal points is straightforward without the application of any boundary integral, since the fundamental solution is the analytical one. This work belongs to a project that aims a computer code for large-scale problems (millions of degrees of freedom). At this stage, some numerical examples were analyzed to evaluate the applicability of the method expeditious its computational effort and convergence of the results for the variables involved in the method. It was developed by the algorithms implementation for potential and elasticity problems. In the case of two-dimensional were employed linear, quadratic and cubic elements and to the three-dimensional case were employed triangular, quadrilateral, linear and quadratic elements in both cases. The computational codes were always implemented focused on solving largescale problems. It is expected that in a final stage of the project with the incorporation procedure of the method fast multipole, it can be more efficiently.

Neste trabalho desenvolveu-se uma formulação adaptativa do método de elementos de contorno (MEC) para a análise de problemas de fratura elástica linear. Foi utilizado o método da colocação para a formulação das equações integrais de deslocamento e de tensão. Para a discretização das equações integrais foram utilizados elementos lineares que possibilitaram a obtenção das expressões exatas das integrais (integração analítica) sobre elementos de contorno e fratura. Para a montagem do sistema de equações algébricas foram utilizadas apenas equações de deslocamento, apenas equações de forças de superfície, ou as duas escritas para nós opostos da fratura levando, portanto ao método dos elementos de contorno dual usualmente empregado na análise de fratura. Para o processo de crescimento da trinca foi desenvolvido um procedimento especial objetivando a correta determinação da direção de crescimento da trinca. Os fatores de intensidade de tensão são calculados por meio da conhecida técnica de correlação de deslocamentos a qual relaciona os deslocamentos atuantes nas faces da fissura. Após a determinação dos fatores de intensidade de tensão é utilizada a teoria da máxima tensão circunferencial para a determinação do ângulo de propagação. O modelo adaptativo empregado é do tipo h onde apenas a sub-divisão dos elementos é feita com base em erros estimados. O erro a ser considerado foi estimado a partir de normas onde se consideraram: a variação aproximada dos deslocamentos, a variação das forças de superfície e a variação da energia de deformação do sistema, calculada com a sua integração sobre o contorno. São apresentados exemplos numéricos para demonstrar a eficiência dos procedimentos propostos. / In this work, an adaptative formulation of the boundary element method is developed to analyze linear elastic fracture problems. The collocation point method was used to formulate the integral equations for the displacements and stresses (or tractions). To discretize the integral equations, linear elements were used to obtain the exact expressions of the integrals over boundary elements and fracture. To construct the linear system of equations were used only displacement equations, traction equations or both of them written for opposite nodes of the fracture, leading to the dual boundary element formulation usually employed in the fracture analyses. For the process of growth of the crack a special procedure was developed aiming at the correct determination of the direction of growth of the crack. The stress intensity factors, to calculate he crack growth angle, are calculated through of correlation displacements technique which relates the displacements actuants in the faces of the crack. The employed adaptative model is the h-type where only the sub-division of the elements is done based on error estimate. The error estimates considered in this work are based on the following norms: displacement, traction and strain energy variations, this last considered from the integration over the boundary. Numerical examples are presented to demonstrate the efficiency of the proposed procedures.

Neste trabalho, é desenvolvido um programa de computador para a análise estática e tridimensional de problemas de interação solo-estrutura. O programa permite considerar várias camadas de solo, cada qual com características físicas diferentes. Sobre este solo, o qual pode conter estacas, podem ser apoiados diversos tipos de estruturas, tais como placas e até um edifício. Todos os materiais considerados são homogêneos, isotrópicos, elásticos e lineares. O solo tridimensional é modelado com o método dos elementos de contorno (MEC), empregando as soluções fundamentais de Kelvin e uma técnica alternativa na consideração do maciço não-homogêneo. Esta técnica, que é uma contribuição original deste trabalho, é baseada no relacionamento das soluções fundamentais de deslocamento dos diferentes domínios, permitindo que sejam analisados como um único sólido sem a necessidade de equações de equilíbrio e compatibilidade. Isso reduz o sistema de equações final e melhora a precisão dos resultados, conforme comprovado nos exemplos apresentados. Para reduzir o custo computacional sem prejudicar a precisão dos resultados, é utilizada uma malha de elementos de contorno infinitos (ECI) nas bordas da malha de ECs para modelar o comportamento das variáveis de campo em longas distâncias. A formulação do ECI mapeado utilizado é outra contribuição original deste trabalho, sendo baseado em um EC triangular. É demonstrado por meio de exemplos que tal formulação é eficiente para a redução de malha, contribuindo de forma significativa na redução do custo computacional. Todas as estruturas que interagem com o solo, incluindo as de fundação, são simuladas empregando o método dos elementos finitos (MEF). Cada estaca é modelada como uma linha de carga empregando um único elemento finito com 14 parâmetros nodais, o qual utiliza funções de forma do quarto grau para aproximar os deslocamentos horizontais, do terceiro grau para as forças horizontais e deslocamentos verticais, do segundo grau para as forças cisalhantes verticais e constantes para as reações da base. Este elemento é empregado em outros trabalhos, no entanto os autores utilizam as soluções fundamentais de Mindlin na consideração da presença da estaca no solo. Desta forma, a formulação desenvolvida neste trabalho com as soluções fundamentais de Kelvin pode ser considerada mais uma contribuição original. No edifício, que pode incluir um radier como estrutura de fundação, são utilizados dois tipos de EFs. Os pilares e vigas são simulados com elementos de barra, os quais possuem dois nós e seis graus de liberdade por nó. As lajes e o radier são modelados empregando elementos planos, triangulares e com três nós. Nestes EFs triangulares são superpostos efeitos de membrana e flexão, totalizando também seis graus de liberdade por nó. O acoplamento MEC/MEF é feito transformando as cargas de superfície do MEC em carregamentos nodais reativos no MEF. Além de exemplos específicos nos Capítulos teóricos, um Capítulo inteiro é dedicado a demonstrar a abrangência e precisão da formulação desenvolvida, comparando-a com resultados de outros autores. / In this work, a computer code is developed for the static analysis of three-dimensional soil-structure interaction problems. The program allows considering a layered soil, which may contain piles. This soil may support several structures, such as shells or even an entire building. All materials are considered homogeneous, isotropic, elastic and linear. The three-dimensional soil is modeled with the boundary element method (BEM), employing Kelvin fundamental solutions and an alternative multi-region technique. This technique, which is an original contribution of this work, is based on relating the displacement fundamental solution of the different domains, allowing evaluating them as an unique solid and not requiring compatibility or equilibrium equations. In such a way, the final system of equations is reduced and more accurate results are obtained, as demonstrated in the presented examples. In order to reduce the computational cost maintaining the accuracy, an infinite boundary element (IBE) mesh is employed at the BE mesh limits to model the far field behavior. The mapped IBE utilized, based on a triangular EC, is another original contribution of this work. In the presented examples it is demonstrated that this IBE formulation is efficient for mesh reduction, implying on a significant computational cost reduction. All structures that interact with the soil, including the foundations, are simulated with de finite element method (FEM). The piles are modeled using a one-dimensional 14 parameter finite element, with forth degree shape functions for horizontal displacement approximation, third degree shape functions for horizontal forces and vertical displacement, second degree shape functions for vertical share force, and constant for the base reaction. This element is employed in other works, however the authors utilize Mindlin fundamental solutions for the pile presence consideration in the soil. In such a way, the formulation developed in this work with Kelvin fundamental solutions may be considered one more original contribution. The building, which may include a radier as a foundation structure, is modeled using two types os FEs. Piles and beams are simulated using bar FEs with two nodes and six degrees of freedom per node. The radier and pavements are modeled employing plane triangular three-node FEs. In these FEs plate and membrane effects are superposed, totalizing six degrees of freedom per node. FEM/BEM coupling is made by transforming the BEM tractions in nodal reactions in the FEM. Even though specific examples are presented in the theoretical Chapters, a role Chapter is dedicated for demonstrating the formulation accuracy and coverage. In most examples, the results are compared with the ones obtained by other authors.

Este trabalho apresenta uma nova metodologia para elastografia virtual em imagens simuladas de ultrassom utilizando métodos numéricos e métodos de visão computacional. O objetivo é estimar o módulo de elasticidade de diferentes tecidos tendo como entrada duas imagens da mesma seção transversal obtidas em instantes de tempo e pressões aplicadas diferentes. Esta metodologia consiste em calcular um campo de deslocamento das imagens com um método de fluxo óptico e aplicar um método iterativo para estimar os módulos de elasticidade (análise inversa) utilizando métodos numéricos. Para o cálculo dos deslocamentos, duas formulações são utilizadas para fluxo óptico: Lucas-Kanade e Brox. A análise inversa é realizada utilizando duas técnicas numéricas distintas: o Método dos Elementos Finitos (MEF) e o Método dos Elementos de Contorno (MEC), sendo ambos implementados em Unidades de Processamento Gráfico de uso geral, GpGPUs ( \"General Purpose Graphics Units\" ). Considerando uma quantidade qualquer de materiais a serem determinados, para a implementação do Método dos Elementos de Contorno é empregada a técnica de sub-regiões para acoplar as matrizes de diferentes estruturas identificadas na imagem. O processo de otimização utilizado para determinar as constantes elásticas é realizado de forma semi-analítica utilizando cálculo por variáveis complexas. A metodologia é testada em três etapas distintas, com simulações sem ruído, simulações com adição de ruído branco gaussiano e phantoms matemáticos utilizando rastreamento de ruído speckle. Os resultados das simulações apontam o uso do MEF como mais preciso, porém computacionalmente mais caro, enquanto o MEC apresenta erros toleráveis e maior velocidade no tempo de processamento. / This thesis presents a new methodology for computational elastography applied to simulated ultrasound images, using numerical methods and comptuter vision methods. The aim is to estimate the elastic moduli of diferent tissues using two diferent images of the same cross section acquired in diferent times and pressure conditions. The proposed methodology consists in evaluate the displacement field using optical flow techniques and then apply an inverse analysis using a numerical method. In order to evaluate the displacement field, two distinct formulations for optical flow are used: Lucas-Kanade and Brox. For the inverse analysis problem, the Finite Element Method and the Boundary Element Method are used, both implemented in general purpose graphic units, GpGPUs. Considering a number of materials that may be present in the images, the multiresgions boundary element method is used in order to couple diferent matrices for diferent materials. The optimization process is evaluated using complex variable method. The methodology is validated in three diferent steps: noiseless simulations; additive white gaussian noise simulations; and ultrasound mathematical phantom with speckle tracking. The results show that the Finite Element Method presents more accurate estimatives but a high computational cost, while the Boundary Element Method presents tolerable errors but a better processing time.

O presente trabalho apresenta uma formulação misto do MEC (Método dos Elementos de Contorno) e o MEF (Método dos Elementos Finitos). Nessa formulação, as estacas são modeladas através do MEF como elementos de barra e o solo através do MEC, como um meio contínuo, elástico linear, isótropo e homogêneo, utilizando as soluções fundamentais de MINDLIN (1936). Os sistemas de equações do solo e das estacas para elementos verticais são apresentados como uma combinação de ambos, originando um único sistema final de equações. Apresentam-se também as modificações necessárias para um sistema composto por estacas inclinadas. Após a resolução do sistema final, obtém-se os deslocamentos e as tensões de contato solo-estaca. A seguir, apresentam-se alguns exemplos numéricos obtidos a partir da formulação proposta e compara-se com modelos de outros autores. / This work presents a hybrid formulation of BEM (Boundary Elements Method) and FEM (Finite Elements Method). In that formulation, the piles are modeled through FEM as bar elements and the soil through BEM, as an isotropic, homogeneous, semi-infinite and linear-elastic continuum, using the fundamental solutions of MINDLIN (1936). The systems of equations of the soil and of the piles for vertical elements are presented as a combination of both, originating a single final system of equations. Some modifications are accomplished for the system of inclinated piles. After the resolution of the final system, the displacements and the contact tensions between soil and pile are obtained. Numeric examples are obtained starting from the proposed formulation and to proceed they are compared with other authors\' models.

Made available in DSpace on 2016-12-23T14:08:09Z (GMT). No. of bitstreams: 1 Parte Inicial.pdf: 580643 bytes, checksum: 1783483d80317ac5307ad55e7cbdb752 (MD5) Previous issue date: 2013-03-25 / O propósito da pesquisa aqui elaborada é mostrar a viabilidade da aplicação de Funções de Base Radial de Suporte Compacto (FBRSC) no processo de aproximação direta do núcleo da ação de domínio através de integração de contorno. Essa formulação utilizada no tratamento da integral de domínio é denominada como (Método dos Elementos de Contorno com Integração Direta de Contorno) MECIC. Com o intuito de se avaliar a efetividade das FBRSC como funções de interpolação, serão realizados diversos testes numéricos, onde se deseja calcular o volume de superfícies. Então, serão realizados testes bidimensionais de aproximação, variando-se o suporte das FBRSCs, a fim de analisar o comportamento dessas funções. Depois de verificar a efetividade e a precisão das FBRSCs no processo de interpolação, desenvolvem-se programas, no ambiente do Método dos Elementos de Contorno (MEC), para a solução de problemas governados pela Equação de Poisson com a Formulação MECIC associada ao conceito de interpolação com FBRSC com suporte devidamente otimizados. A aferição das soluções numéricas obtidas se dá a partir da comparação com as suas respectivas soluções analíticas, facilmente encontradas na literatura especializada. Assim, possibilita-se estimar o erro relativo e então a eficácia da Formulação MECIC com FBRSC. Uma vez comprovado a sua eficácia, a Formulação MECIC com FBRSC é testada também com o esquema de interpolação com ajuste de pontos. Durante todo o desenvolvimento, atenta-se para a importância do custo computacional da formulação, a partir da geração de tabelas com o tempo de processamento dos programas implementados no MEC. Dessa forma, avalia-se qualitativamente o desempenho das FBRSC na Formulação MECIC, visando futuras aplicações na área de propagação de ondas sísmicas / The purpose of this research is to show the viability of application of Compactly Supported Radial Basis Function (CSRBF) in the process of direct approximation of the core of the domain action through boundary integration. This formulation is termed as (Boundary Elements Method with Directs Boundary Integration) MECIC, and is used in the treatment of the domain integration. By evaluating the effectiveness of CSRBF as interpolation functions, it performed several numerical tests to calculate the volume of surfaces. Also; by varying the support of CSRBFs, it performed two-dimensional approximation tests to examine the behavior of these functions. After verifying the effectiveness and accuracy of CSRBFs in the interpolation process, it developed computational programs to solve physical problems using the MECIC formulation, which is governed by Poisson s Equation. That formulation is associated with the concept of CSRBF in which the support is properly optimized. The calibration of the numerical solutions is given by the confrontation with their respective analytical solutions, easily found in the specialized literature. In this way, it is possible to estimate the relative error and the effectiveness of the MECIC formulation in association with the CSRBF concept. It is tested also with the curve fitting interpolation scheme. Owing the importance of the computational cost of that formulation, it is generated several time tables showing the processing time of those Boundary Elements Method computational programs. Therefore; aiming future applications in the seismic propagation wave area, it was finally evaluated the qualitative performance of the CSRBF in MECIC s formulation

Neste trabalho desenvolveu-se uma formulação adaptativa do método de elementos de contorno (MEC) para a análise de problemas de fratura elástica linear. Foi utilizado o método da colocação para a formulação das equações integrais de deslocamento e de tensão. Para a discretização das equações integrais foram utilizados elementos lineares que possibilitaram a obtenção das expressões exatas das integrais (integração analítica) sobre elementos de contorno e fratura. Para a montagem do sistema de equações algébricas foram utilizadas apenas equações de deslocamento, apenas equações de forças de superfície, ou as duas escritas para nós opostos da fratura levando, portanto ao método dos elementos de contorno dual usualmente empregado na análise de fratura. Para o processo de crescimento da trinca foi desenvolvido um procedimento especial objetivando a correta determinação da direção de crescimento da trinca. Os fatores de intensidade de tensão são calculados por meio da conhecida técnica de correlação de deslocamentos a qual relaciona os deslocamentos atuantes nas faces da fissura. Após a determinação dos fatores de intensidade de tensão é utilizada a teoria da máxima tensão circunferencial para a determinação do ângulo de propagação. O modelo adaptativo empregado é do tipo h onde apenas a sub-divisão dos elementos é feita com base em erros estimados. O erro a ser considerado foi estimado a partir de normas onde se consideraram: a variação aproximada dos deslocamentos, a variação das forças de superfície e a variação da energia de deformação do sistema, calculada com a sua integração sobre o contorno. São apresentados exemplos numéricos para demonstrar a eficiência dos procedimentos propostos. / In this work, an adaptative formulation of the boundary element method is developed to analyze linear elastic fracture problems. The collocation point method was used to formulate the integral equations for the displacements and stresses (or tractions). To discretize the integral equations, linear elements were used to obtain the exact expressions of the integrals over boundary elements and fracture. To construct the linear system of equations were used only displacement equations, traction equations or both of them written for opposite nodes of the fracture, leading to the dual boundary element formulation usually employed in the fracture analyses. For the process of growth of the crack a special procedure was developed aiming at the correct determination of the direction of growth of the crack. The stress intensity factors, to calculate he crack growth angle, are calculated through of correlation displacements technique which relates the displacements actuants in the faces of the crack. The employed adaptative model is the h-type where only the sub-division of the elements is done based on error estimate. The error estimates considered in this work are based on the following norms: displacement, traction and strain energy variations, this last considered from the integration over the boundary. Numerical examples are presented to demonstrate the efficiency of the proposed procedures.

Este trabalho apresenta uma nova metodologia para elastografia virtual em imagens simuladas de ultrassom utilizando métodos numéricos e métodos de visão computacional. O objetivo é estimar o módulo de elasticidade de diferentes tecidos tendo como entrada duas imagens da mesma seção transversal obtidas em instantes de tempo e pressões aplicadas diferentes. Esta metodologia consiste em calcular um campo de deslocamento das imagens com um método de fluxo óptico e aplicar um método iterativo para estimar os módulos de elasticidade (análise inversa) utilizando métodos numéricos. Para o cálculo dos deslocamentos, duas formulações são utilizadas para fluxo óptico: Lucas-Kanade e Brox. A análise inversa é realizada utilizando duas técnicas numéricas distintas: o Método dos Elementos Finitos (MEF) e o Método dos Elementos de Contorno (MEC), sendo ambos implementados em Unidades de Processamento Gráfico de uso geral, GpGPUs ( \"General Purpose Graphics Units\" ). Considerando uma quantidade qualquer de materiais a serem determinados, para a implementação do Método dos Elementos de Contorno é empregada a técnica de sub-regiões para acoplar as matrizes de diferentes estruturas identificadas na imagem. O processo de otimização utilizado para determinar as constantes elásticas é realizado de forma semi-analítica utilizando cálculo por variáveis complexas. A metodologia é testada em três etapas distintas, com simulações sem ruído, simulações com adição de ruído branco gaussiano e phantoms matemáticos utilizando rastreamento de ruído speckle. Os resultados das simulações apontam o uso do MEF como mais preciso, porém computacionalmente mais caro, enquanto o MEC apresenta erros toleráveis e maior velocidade no tempo de processamento. / This thesis presents a new methodology for computational elastography applied to simulated ultrasound images, using numerical methods and comptuter vision methods. The aim is to estimate the elastic moduli of diferent tissues using two diferent images of the same cross section acquired in diferent times and pressure conditions. The proposed methodology consists in evaluate the displacement field using optical flow techniques and then apply an inverse analysis using a numerical method. In order to evaluate the displacement field, two distinct formulations for optical flow are used: Lucas-Kanade and Brox. For the inverse analysis problem, the Finite Element Method and the Boundary Element Method are used, both implemented in general purpose graphic units, GpGPUs. Considering a number of materials that may be present in the images, the multiresgions boundary element method is used in order to couple diferent matrices for diferent materials. The optimization process is evaluated using complex variable method. The methodology is validated in three diferent steps: noiseless simulations; additive white gaussian noise simulations; and ultrasound mathematical phantom with speckle tracking. The results show that the Finite Element Method presents more accurate estimatives but a high computational cost, while the Boundary Element Method presents tolerable errors but a better processing time.