Uma das mais recentes abordagens propostas na literatura para tratar o problema do alto número de Weissenberg (We) é a Formulação Log-Conformação (FLC). Nesta formulação, a equação constitutiva viscoelástica utilizada é reescrita em termos de uma variável Ψ, que é o logaritmo do tensor conformação. Apesar do potencial de aplicação da FLC, pouca atenção tem sido dirigida para análise da acurácia da solução obtida para o campo de tensões quando se utiliza esta formulação. Assim, o objetivo do presente trabalho foi estudar a acurácia da solução obtida pela FLC na análise de escoamentos de fluidos viscoelásticos usando duas geometrias padrão de estudo: placas paralelas e cavidade quadrada com tampa móvel. Primeiramente, a FLC foi implementada no pacote de CFD OpenFOAM. Em seguida foram verificados os limites do número de Weissenberg na formulação numérica padrão (Welim,P), onde para a geometria de placas paralelas foi encontrado Welim,P = 0,3 e para a geometria da cavidade quadrada com tampa móvel foi encontrado Welim,P = 0,8. Depois o código implementado foi aplicado em ambas as geometrias, comparando-se a solução obtida pela FLC com aquela da formulação padrão na faixa de We < Welim,P. Os resultados obtidos na geometria de placas paralelas apresentaram boa concordância com a solução padrão e solução analítica. Para a geometria da cavidade quadrada com tampa móvel, que não possui solução analítica, boa concordância dos resultados também foi observada em comparação com a solução padrão. Posteriormente foram comparados os resultados obtidos pela FLC na faixa de We > Welim,P. Na geometria de placas paralelas, além da boa concordância com a solução analítica, obteve-se convergência em todos os casos estudados neste trabalho, com o maior número de Weissenberg utilizado sendo igual a 8 Os resultados da geometria da cavidade quadrada com tampa móvel também apresentaram boa concordância em comparação com dados da literatura, porém a convergência foi obtida até para We = 2. Com respeito à comparação das formulações numéricas com a solução analítica, feita apenas na geometria de placas paralelas, foi observado um erro máximo de 7,57% na solução padrão e de 12,33% na FLC. Em relação à análise da qualidade das tensões usando os resíduos da equação constitutiva viscoelástica como critério de acurácia, foi verificado nas duas geometrias que os valores de tensão obtidos usando a FLC são menos acurados que aqueles obtidos pela formulação explícita no tensor das tensões nos casos em que esta última converge. Também foi observado que a acurácia diminui com o aumento do We. Esse efeito pôde ser melhor notado na geometria de placas paralelas. Uma razão para a perda de acurácia da tensão provavelmente ocorre devido à natureza matemática da transformação algébrica inversa de Ψxx para τxx. O novo solver implementado neste trabalho apresentou convergência e soluções corretas para as duas geometrias, logo foi implementado corretamente. Ele também potencializa o solver de partida viscoelastiFluidFoam ao estender simulações para uma faixa maior do número de Weissenberg. / A recent approach proposed in the literature to deal with the High Weissenberg Number Problem is the Log-Conformation formulation (LCF). In this formulation the viscoelastic constitutive equation is rewritten in terms of the logarithm of the conformation tensor Ψ. Despite the great potential application of the LCF, little attention has been given in the literature to the accuracy of the obtained stress fields. The purpose of this work was to study the solution obtained by LCF in the analysis of viscoelastic flows using two benchmark geometries: parallel plates and lid driven cavity. Firstly, the LCF was implemented in the OpenFOAM CFD package. Then, the limits of Weissenberg number for the standard numerical formulation (Welim,P) were verified, obtaining Welim,P = 0.3 for the parallel plates and Welim,P = 0.8 for the lid driven cavity. When comparing the solution obtained by the LCF with that of the standard formulation in a range of We < Welim,P, the results obtained for the parallel plates geometry showed good agreement with the standard solution and the analytical solution. For the lid driven cavity geometry, for which there is not analytical solution, good agreement with the standard solution was also observed. For We > Welim,P in the parallel plates geometry, in addition to the good agreement with the analytical solution, it was possible to obtain convergence in all the cases studied in this work, with the largest number of Weissenberg used being equal to 8 The results of the lid driven cavity geometry also presented good agreement in comparison with literature data, but convergence was obtained up to We = 2. With respect to the comparison of the numerical formulations with the analytical solution for the parallel plates geometry, a maximum error of 7.57% was observed in the standard solution and of 12.33% in the LCF. When using the residues of the viscoelastic constitutive equation as a criterion of accuracy, it was verified that for the two geometries the stress values obtained using the LCF were less accurate than those obtained by the explicit formulation in the stress tensor. It has also been observed that accuracy decreases with increasing of We. One reason for the loss of stress accuracy probably occurs because of the mathematical nature of the inverse algebraic transformation from Ψxx to τxx. The new solver implemented in this work presented convergence and correct solutions for the two geometries, so it was implemented correctly. It also potentiates the viscoelastiFluidFoam starting solver by extending simulations to a larger range of Weissenberg number.

Orientador: Prof. Dr. Caetano Rodrigues Miranda / Dissertação (mestrado) - Universidade Federal do ABC, Programa de Pós-Graduação em Física, 2014.

Made available in DSpace on 2016-04-01T17:54:45Z (GMT). No. of bitstreams: 0 Previous issue date: 2008-06-20. Added 1 bitstream(s) on 2016-04-01T18:00:27Z : No. of bitstreams: 1 000556918.pdf: 1923814 bytes, checksum: 4705205551ed1cbedcc59e194357522a (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / O trabalho é baseado na elaboração de modelos matemáticos de suspensão que permitem avaliar parâmetros de desempenho de um veículo agrícola em função da excitação provocada pelas irregularidades do solo. Neste trabalho, estas irregularidades são representadas por entrada degrau, entrada seno e obstáculos de uma pista normalizada segundo a norma ISO 5008. As equações do movimento são obtidas a partir dos modelos matemáticos de um quarto de veículo e meio veículo. As simulações numéricas são executadas nos softwares Matlab®, Simulink® e Visual Nastran®. A partir da entrada conhecida, podem-se determinar as características dos elementos da suspensão para obter níveis desejáveis de conforto e segurança. Foram analisadas diferentes configurações do sistema, variando-se a relação de rigidez e a constante de amortecimento equivalente a partir de um modelo considerado padrão. Os resultados foram tratados estatisticamente e comparados entre si, verificando-se que o aumento da relação de rigidez resulta na redução da aceleração vertical e no aumento do curso da suspensão, melhorando o conforto e diminuindo a segurança. Por outro lado, a redução do amortecimento provoca a redução da aceleração vertical até um determinado ponto e, a partir deste, a aceleração volta a crescer. Já o curso da suspensão diminui à medida que o amortecimento é reduzido. / This paper is based on the elaboration of mathematical suspension models which can evaluate performance parameters of an agricultural vehicle based on the vibration caused by irregularities on the soil surface. In this paper, these irregularities are represented by step input, sine input and obstacles of leveled road according to ISO 5008 specifications. The motion equations are obtained from the mathematical models of a quarter car and half car. Numerical simulations are performed on Matlab®, Simulink® and Visual Nastran®. Based on known input, characteristics of suspension elements can be determined in order to obtain desired levels of comfort and safety. Different system configurations were analyzed changing the stiffness ratio and the damping coefficient based on a standard model. Results were treated statistically and compared among themselves. It was verified that the increase in stiffness ratio results in the reduction of the vertical acceleration and in the increase of the suspension travel, improving comfort but reducing safety. On the other hand, the reduction of damping coefficient causes the reduction of vertical acceleration up to a certain point, and from this point on the acceleration increases again. The suspension travel, however, reduces as damping coefficient is reduced.

Simulações Numéricas são executadas em um código numérico de alta precisão resolvendo as equações de Navier-Stokes e da continuidade para regimes de escoamento incompressíveis num contexto da turbulência bidimensional. Este código utiliza um esquema compacto de diferenças finitas de sexta ordem na aproximação das derivadas espaciais. As derivadas temporais são calculadas usando o esquema de Runge-Kuta de terceeira ordem com baixo armazenamento. Tal código numérico fornece uma representação melhorada para uma grande faixa de escalas de comprimento e de tempo. As técnicas dos contornos imersos acopladas ao método dos contornos virtuais permitem modelar escoamentos não-estacionários sobre geometrrias complexas, usando simplesmente uma malha Cartesiana uniforme. Por meio de procedimentos de aproximação/interpolação, as técnicas dos contornos imersos (aproximação Gaussiana, interpolação bilinear e redistribuição Gaussiana), permitem a representação do corpo sólido no interior do campo de escoamento, com a superfície não coincidindo com a malha computacional. O método dos contornos virtuais, proposto originalmente por Peskin, consiste, basicamente, na imposição na superfície e/ou no interior do corpo, de um termo de força temporal acrescentando às equações do momento. A aplicação deste campo de força local leva o fluido ao repouso na superfície do corpo, permitindo obter as condições de contorno de não-deslizamento e de não penetração de fluido na parede. A análise das oscilações induzidas no escoamento-contorno pelo processo de desprendimento de vórtices na esteira do cilindro circular e de geometria retangulares na incidência, para números de Reybolds variando de 40 a 400, confirma a eficiência computacional e a aplicabilidade das técncias implementadas.

Submitted by Milena Rubi (milenarubi@ufscar.br) on 2016-11-03T17:25:54Z No. of bitstreams: 1 SAMPAIO_Guilherme_2016.pdf: 12947275 bytes, checksum: e37468f50e4625a4fdc9875a5cf3f85f (MD5) / Approved for entry into archive by Milena Rubi (milenarubi@ufscar.br) on 2016-11-03T17:26:13Z (GMT) No. of bitstreams: 1 SAMPAIO_Guilherme_2016.pdf: 12947275 bytes, checksum: e37468f50e4625a4fdc9875a5cf3f85f (MD5) / Approved for entry into archive by Milena Rubi (milenarubi@ufscar.br) on 2016-11-03T17:26:21Z (GMT) No. of bitstreams: 1 SAMPAIO_Guilherme_2016.pdf: 12947275 bytes, checksum: e37468f50e4625a4fdc9875a5cf3f85f (MD5) / Made available in DSpace on 2016-11-03T17:26:28Z (GMT). No. of bitstreams: 1 SAMPAIO_Guilherme_2016.pdf: 12947275 bytes, checksum: e37468f50e4625a4fdc9875a5cf3f85f (MD5) Previous issue date: 2016-05-20 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / The comprehension of the phenomena related to magnetism allowed advances in various fields of knowledge and several modern applications that are now part of our daily lives. Specifically, the magnetic nanoparticles (NPs) has aroused a great interest among researchers in recent decades. This curiosity due to its unique properties is due to its small size and surface interactions. One of the main physical phenomena that appear when we reduce the size of a magnetic particle is called superparamagnetism. In this work, we have looked for a description and understanding of the mathematical formalism attached to magnetic phenomena at nanoscale through the application of the Langevin’s function to magnetic NPs in superparamagnetic regime and non-interacting model. Accordingly, we use theoretical simulations of magnetization curves as a function of applied field to study some parameters and verify its compatibility with magnetic characterization data obtained for samples of NPs magnetite and maghemite. In this way, we can get relevant information to the analysis of the magnetic behavior of nanoparticulate systems. Finally, we have lunched the potential of these magnetic NPs in biomedical applications both in diagnostic techniques and in disease treatments. / A compreensão dos fenômenos ligados ao magnetismo permitiu avanços nas mais variadas áreas do conhecimento e diversas aplicações modernas que hoje fazem parte do nosso cotidiano. Especificamente, as nanopartículas (NPs) magnéticas têm despertado um grande interesse por parte dos pesquisadores nas últimas décadas. Essa curiosidade devido às suas propriedades diferenciadas é decorrente de suas dimensões reduzidas e de interações de superfície. Um dos principais fenômenos físicos que surgem ao se reduzir o tamanho de partículas magnéticas é o chamado superparamagnetismo. Neste trabalho, procuramos por uma descrição e entendimento do formalismo matemático ligado ao fenômeno magnético em nanoescala por meio da aplicação do modelo proposto por Langevin a NPs magnéticas em regime superparamagnético não-interagente. Nesse sentido, utilizamos de simulações teóricas de curvas de magnetização em função do campo aplicado a fim de estudar alguns parâmetros e verificar sua compatibilidade com dados de caracterização magnética obtidos para amostras de NPs de magnetita e maghemita. Desta forma, podemos obter informações relevantes para a análise do comportamento magnético desses sistemas nanoparticulados. Por fim, evidenciamos as potencialidades dessas NPs magnéticas nas aplicações biomédicas tanto em técnicas de diagnóstico quanto em tratamentos de doenças.

Sistemas de ensaios não destrutivos por ondas guiadas despertam cada vez mais a atenção tanto da indústria, como da academia. Isso deve-se, principalmente, às possibilidades que as ondas guiadas permitem, como maior área de triagem que o ultrassom convencional. Porém a técnica tem suas limitações, já que esta apenas gera uma estimativa da localização de um defeito e não a sua dimensão. Nessa limitação, algoritmos tomográficos apresentam uma possibilidade de avanço da técnica, pois permitem determinar não só a localização de corrosões e defeitos, mas também seu dimensionamento. Este trabalho apresenta tomografia de difração como uma alternativa para avaliação de integridade estrutural. Primeiramente, utilizou-se análise por métodos numéricos para mostrar a validade dos algoritmos e posteriormente foi realizado um experimento em uma chapa real com o objetivo de reconstruir a imagem do defeito. / Guided waves nondestructive testing systems are increasingly attracting industrial and academic attention. The mainly reason for this attention is the possibility of screening a large area than conventional ultrasound technique. However, Guided Waves Testing has limitations, since it gives only an estimation of the location of a defect, but not the dimensions. Tomographic algorithms come up with an improvement of the technique, because it allows discovery not only the location of corrosions and defects, but the dimensions too. This work brings Diffraction Tomography as an alternative to structural health monitoring. First, a numerical analysis was implemented to demonstrate the validity of the algorithms, after that an experiment in a real plate was made with the objective to recover the defect image.

Uma das mais recentes abordagens propostas na literatura para tratar o problema do alto número de Weissenberg (We) é a Formulação Log-Conformação (FLC). Nesta formulação, a equação constitutiva viscoelástica utilizada é reescrita em termos de uma variável Ψ, que é o logaritmo do tensor conformação. Apesar do potencial de aplicação da FLC, pouca atenção tem sido dirigida para análise da acurácia da solução obtida para o campo de tensões quando se utiliza esta formulação. Assim, o objetivo do presente trabalho foi estudar a acurácia da solução obtida pela FLC na análise de escoamentos de fluidos viscoelásticos usando duas geometrias padrão de estudo: placas paralelas e cavidade quadrada com tampa móvel. Primeiramente, a FLC foi implementada no pacote de CFD OpenFOAM. Em seguida foram verificados os limites do número de Weissenberg na formulação numérica padrão (Welim,P), onde para a geometria de placas paralelas foi encontrado Welim,P = 0,3 e para a geometria da cavidade quadrada com tampa móvel foi encontrado Welim,P = 0,8. Depois o código implementado foi aplicado em ambas as geometrias, comparando-se a solução obtida pela FLC com aquela da formulação padrão na faixa de We < Welim,P. Os resultados obtidos na geometria de placas paralelas apresentaram boa concordância com a solução padrão e solução analítica. Para a geometria da cavidade quadrada com tampa móvel, que não possui solução analítica, boa concordância dos resultados também foi observada em comparação com a solução padrão. Posteriormente foram comparados os resultados obtidos pela FLC na faixa de We > Welim,P. Na geometria de placas paralelas, além da boa concordância com a solução analítica, obteve-se convergência em todos os casos estudados neste trabalho, com o maior número de Weissenberg utilizado sendo igual a 8 Os resultados da geometria da cavidade quadrada com tampa móvel também apresentaram boa concordância em comparação com dados da literatura, porém a convergência foi obtida até para We = 2. Com respeito à comparação das formulações numéricas com a solução analítica, feita apenas na geometria de placas paralelas, foi observado um erro máximo de 7,57% na solução padrão e de 12,33% na FLC. Em relação à análise da qualidade das tensões usando os resíduos da equação constitutiva viscoelástica como critério de acurácia, foi verificado nas duas geometrias que os valores de tensão obtidos usando a FLC são menos acurados que aqueles obtidos pela formulação explícita no tensor das tensões nos casos em que esta última converge. Também foi observado que a acurácia diminui com o aumento do We. Esse efeito pôde ser melhor notado na geometria de placas paralelas. Uma razão para a perda de acurácia da tensão provavelmente ocorre devido à natureza matemática da transformação algébrica inversa de Ψxx para τxx. O novo solver implementado neste trabalho apresentou convergência e soluções corretas para as duas geometrias, logo foi implementado corretamente. Ele também potencializa o solver de partida viscoelastiFluidFoam ao estender simulações para uma faixa maior do número de Weissenberg. / A recent approach proposed in the literature to deal with the High Weissenberg Number Problem is the Log-Conformation formulation (LCF). In this formulation the viscoelastic constitutive equation is rewritten in terms of the logarithm of the conformation tensor Ψ. Despite the great potential application of the LCF, little attention has been given in the literature to the accuracy of the obtained stress fields. The purpose of this work was to study the solution obtained by LCF in the analysis of viscoelastic flows using two benchmark geometries: parallel plates and lid driven cavity. Firstly, the LCF was implemented in the OpenFOAM CFD package. Then, the limits of Weissenberg number for the standard numerical formulation (Welim,P) were verified, obtaining Welim,P = 0.3 for the parallel plates and Welim,P = 0.8 for the lid driven cavity. When comparing the solution obtained by the LCF with that of the standard formulation in a range of We < Welim,P, the results obtained for the parallel plates geometry showed good agreement with the standard solution and the analytical solution. For the lid driven cavity geometry, for which there is not analytical solution, good agreement with the standard solution was also observed. For We > Welim,P in the parallel plates geometry, in addition to the good agreement with the analytical solution, it was possible to obtain convergence in all the cases studied in this work, with the largest number of Weissenberg used being equal to 8 The results of the lid driven cavity geometry also presented good agreement in comparison with literature data, but convergence was obtained up to We = 2. With respect to the comparison of the numerical formulations with the analytical solution for the parallel plates geometry, a maximum error of 7.57% was observed in the standard solution and of 12.33% in the LCF. When using the residues of the viscoelastic constitutive equation as a criterion of accuracy, it was verified that for the two geometries the stress values obtained using the LCF were less accurate than those obtained by the explicit formulation in the stress tensor. It has also been observed that accuracy decreases with increasing of We. One reason for the loss of stress accuracy probably occurs because of the mathematical nature of the inverse algebraic transformation from Ψxx to τxx. The new solver implemented in this work presented convergence and correct solutions for the two geometries, so it was implemented correctly. It also potentiates the viscoelastiFluidFoam starting solver by extending simulations to a larger range of Weissenberg number.

Neste trabalho modela-se simplificadamente o fluxo da Geleira de descarga Lange (Ilha Rei George, Antártica) para estimar seu estado de equilíbrio e a sua resposta a mudanças na taxa de acumulação. O modelo, bidimensional, representativo e numericamente estável, fornece resultados qualitativos para a velocidade e distribuição da temperatura do gelo, baseado no método das diferenças finitas. Para calibrar o modelo, medidas da velocidade da geleira foram obtidas durante a XVIII Operação Antártica Brasileira no verão austral de 1999/2000. Primeiramente, 21 estacas foram implantadas na superfície da Geleira Lange, a posição geográfica levantada por GPS diferencial duas vezes no início e final de um período de dois meses, permitindo a determinação da velocidade superficial. O fluxo de geometria complexa foi então simulado através de um sistema curvilíneo de coordenadas. O modelo computa a distribuição de velocidade ao longo do eixo longitudinal usando a lei do fluxo do gelo para simular a deformação interna da geleira, desconsiderando o desiizamento basal. O modelo indica que variações na taxa de acumulação líquida menores que 20% do valor atual não resultariam em mudanças significativas da morfologia superficial da geleira. Mesmo a redução em 50% na taxa de acumulação resultaria somente no rebaixamento da superfície em 26 m. Estes resultados reforçam a idéia de uma geleira perto do estado de equilíbrio ( steady-state). / A simplified numerical model simulates the flow of the outlet Lange Glacier (King George Island, Antarctica) to estimate its equilibrium state and morphologicai responses to accumulation rate changes. The bidimensional model uses the finite difference method to provide qualitative information about the surface velocity and the internai temperature distribution. Fieldwork, carried out in the austral summer of 1999/2000, during the XVIII Brazilian Antarctic Operation, provides ice velocity data for calibration. Firstly, 21 stakes were fixed on the surface of Lange Glacier, their geographical position surveyed twice within a period of two months, using differential GPS. The glacial flow was simulated using complex curvilinear coordinates. The model computes the velocity distribution along the main longitudinal axis using the ice flow la\v to calculate the internai deformation with no basal sliding and responses to changes in the net accumulation rate. Variations in the net accumulation rate of less than 20% of the present value do not change significantly the surface glacier morphology. Even a reduction of 50% of this rate willlower the surface by only 26 m in 100 years. These results suggest a present glacier condition near to the steady-state.

O objetivo deste trabalho é a obtenção dos modos normais e das freqüências naturais de vigas com condições de contorno não-clássicas, descritas pelo modelo estrutural de Euler-Bernoulli. Para a obtenção dos modos são consideradas duas bases, a base espectral clássica e a base dinâmica associada à resposta impulso do modelo. O deslocamento de uma viga sob ação de uma força periódica é obtido através da análise modal e apresentado de maneira gráfica para as diversas condições de contorno. / The objective of this work is the obtention of the normal modes and natural frequencies of beams with non-classical boundary conditions described by t he structural Euler-Bernoulli model. For this, two basis are considered, a classical spectral basis and the dynamical basis associated with the impulse response of the model. The displacement of a beam under the action of a periodic force is obtained through modal analysis and presented graphically for such various boundary conditions.

Nesta tese estudamos o espectro de absorção da Merocianina de Brooker (MB) e três derivados em solventes com diferentes polaridades. O interesse no estudo desses sistemas se dá pela presença de algumas propriedades particulares das merocianinas, por exemplo, apresentam um grande deslocamento solvatocrômico devido à mudança na polaridade do meio. Nós apresentamos os resultados teóricos obtidos para o espectro de absorção das moléculas considerando diferentes estruturas para o estado fundamental: formas trans, cis, zwiteriônica (zw), protonada (MBH+) e dímero(MB2). O efeito do solvente foi considerado utilizando diferentes modelos de solvatação: (i) modelo contínuo polarizável (PCM), (ii) através de uma configuração eletrostática média do solvente (ASEC), (iii) incluindo algumas moléculas explícitas do solvente juntamente com um ambiente eletrostático gerado pelas demais moléculas do solvente. Como uma investigação adicional, apresentamos nossos resultados de medidas experimentais para o espectro de absorção da MB em vários solventes variando acidez e concentração. Vimos, através de cálculos quânticos, que as formas cis/trans apresentam a excitação eletrônica na mesma região, e que a deformação estrutural gerada pela forma zw provoca um deslocamento para o vermelho na excitação eletrônica. Os nossos resultados teóricos e experimentais mostram que a forma MBH+ apresenta um solvatocromismo pequeno, com deslocamento de 20 nm provocado pela mudança água-clorofórmio. Realizando medidas experimentais através da titulação espectroscópica nós obtivemos o pKa associado ao processo de desprotonação/protonação da MB em água e metanol. Em água obtivemos um valor de 8.7, em boa concordância com valores da literatura. Apresentamos um valor inédito do pKap da MB em metanol, 9.9. Conseguimos uma boa descrição teórica para a excitação eletrônica da MB em solventes com alta polaridade, na região entre 430-500 nm, utilizando o método quântico TD-DFT com funcionais B3LYP e CAM-B3LYP e conjunto de funções base 6-311+G**, porém a excitação em solventes de baixa polaridade, que ocorre na região entre 550-650 nm, não é corretamente descrita considerando a MB nos solventes. Nós vimos, através de estudos experimentais, que a sonda MB pode agregar em solventes de baixa polaridade. Os cáculos teóricos para dímeros em solução mostraram a existência de uma excitação eletrônica de baixa intensidade nesta região. Adicionalmente, os espectros experimentais em solventes de baixa polaridade mostraram 2 bandas, onde a segunda se assemelha com a banda observada para a forma MBH+. Para explicar essas 2 bandas experimentais para a MB, apresentamos uma proposta teórica onde ocorre uma transferência de prótons (H+) entre os monômeros do dímero formado, gerando uma estrutura desprotonada (MBH-) e uma protonada (MBH+). Cálculos teóricos para a forma MBH- mostram que essa forma apresenta uma excitação eletrônica de intensidade moderada na região entre 550-650 nm. Com essa hipótese nós conseguimos descrever, através de cálculos teóricos, o solvatocromismo anômalo observado experimentalmente para o espectro eletrônico de absorção da MB nas duas regiões de polaridade dos solventes: alta polaridade, sendo descrita pela forma MB, e baixa polaridade, descrita pela forma MBH-. / In this thesis we studied the absorption spectrum of merocyanine Brooker (MB) and three derivatives in solvents with different polarities.The interest in this system is given by the presence of some particular properties of this molecule, for example, it presents a large solvatochromic shift due to the change in the polarity of the medium. We present the results for the absorption spectrum of the molecules considering different structures to the ground state: forms trans, cis, zwitterionic (zw), protonated MBH+ and dimer (MB2). The solvent effect was treaty by different ways: (i) continuous model using the PCM polarizable, (ii) by an average electrostatic configuration of solvent, ASEC, (iii) including some explicit solvent molecules with an electrostatic environment generated by other solvent molecules. As an additional investigation, we present results of experimental measurements in the thesis. We have seen, through quantum calculations, that the forms cis/trans have the electronic excitation in the same region and the structural deformation generated in the zw form causes a red shift. Our theoretical and experimental results show that the MBH+ form has a small solvatochromism, with displacement of 20 nm caused by water-chloroform change. Performing spectroscopic titration we got the pKa associated with the process of deprotonation/protonation of MB in water and methanol. In water we obtained a value of 8.7, in good agreement with the values reported in the literature. We present a unpublished pKap for the MB in methanol, 9.9. We got a good theoretical description for electronic excitation of MB in solvents with high polarity, in the region between 430-500 nm, using a method quantum TD-DFT B3LYP and CAM-B3LYP functional whit basic functions set 6-311+G** but the excitation in low polarity solvents, which occurs in the region between 550-650 nm, is not properly described considering the MB form in the solvents. We have seen, through experimental studies that the MB probe can aggregate in low polarity solvents. Theoretical calculations for dimer in solution showed the existence of a low intensity electron excitation in this region. Additionally, the experimental spectra in low polarity solvents showed 2 bands, where in the second band is similar to the observed to MBH+ form. To explain these two experimental bands for MB, we present a theoretical proposal where there is a proton transfer (H+) between the monomers of the dimer, generating a deprotonated structure (MBH-) and a protonated (MBH+). Theoretical calculations for the MBH- show that this form presents an electronic excitation of moderate intensity in the region between 550-650 nm. With this assumption we can describe, through theoretical calculations, the anomalous solvatochromism for the electronic absorption spectrum of MB in the two polarity regions of solvents: high polarity is described by the MB form, and low polarity, described by the MBH- form.