A tomografia por sensoriamento térmico é muito utilizada em diferentes aplicações industriais, tais como a determinação de propriedades térmicas de novos materiais, o controle da produção de calor e a temperatura no processo de manufatura. Entretanto, o emprego de técnicas tomográficas em processos industriais envolvendo transferência de calor ainda carece de métodos robustos e computacionalmente eficientes. Nesse contexto, o principal objetivo deste trabalho é contribuir para o desenvolvimento de uma técnica não intrusiva de medição do coeficiente de convecção a partir de medidas externas de temperatura e fluxo de calor baseada na solução do problema térmico inverso. Para tanto é necessário resolver um problema de condução acoplado a um problema de convecção de calor. Este acoplamento ocorre através do coeficiente de convecção no interior do domínio do problema, cuja determinação pode ser feita através da aplicação de um fluxo de calor e medição das temperaturas resultantes na superfície externa. A tomografia térmica é tratada como um problema de minimização global, cuja função objetivo é um funcional de erro que quantifica a diferença entre as medidas externas não intrusivas (temperatura real) e as medidas calculadas no modelo numérico (temperatura aproximada). A natureza mal condicionada do problema assim formulado se manifesta na superfície de minimização por produzir topologias problemáticas tais como múltiplos mínimos locais, pontos de sela, vales ao redor da solução, platôs, etc. Desse modo, uma técnica bastante sofisticada, capaz de convergir para a solução correta mesmo na presença dessas patologias é necessária para obtenção da solução. Neste trabalho optou-se pelo método de Newton para a minimização deste funcional em que a inversa da matriz Hessiana é substituída por uma pseudo-inversa construída a partir da técnica de Decomposição em Valores Singulares Truncados. Os resultados mostram que a técnica proposta foi capaz de superar os problemas de convergência associados à natureza intrínseca mal condicionada do problema inverso e o coeficiente de convecção foi reconstruído com precisão razoável. / Tomography by thermal sensing is widely used in different industrial applications, such as the determination of thermal properties of new materials, the control of heat production and the temperature in manufacturing processes. However, the application of tomographic techniques in industrial processes involving heat transfer still lacks robust and computationally efficient methods. In this context, the main objective of this thesis is to contribute to the development of a non-intrusive technique for measuring of the convection coefficient from external temperature and heat flow measurements based on the solution of the inverse thermal problem. This requires solving a conduction problem coupled with a heat convection problem, which is coupled through an internal convection coefficient, determined by applying a heat flux and measuring the resulting temperatures on the external boundary. The thermal tomography is treated as a global minimization problem in which the fitness function is an error functional that quantifies the difference between non-intrusive external measurements (actual temperature) and measurements calculated in a numerical model (approximate temperature). The ill-conditioned nature of the problem manifests itself in the minimization problem for producing problematic topologies, such as multiple local minima, saddle points, valleys around the solution, plateaus, etc. Thus, a very sophisticated technique that can converge to the correct solution even in the presence of these pathologies is necessary to obtain the solution. In this thesis the Newton\'s method was used for the minimization of this functional in which the inverse Hessian matrix was replaced by a pseudo-inverse built from the truncated singular value decomposition technique. Results show that the proposed technique was capable of overcoming the convergence problems associated with the intrinsic ill-conditioned nature of the inverse problem and the convection coefficient was reconstructed within reasonable precision.

Submitted by Maria Cristina (library@lncc.br) on 2016-07-27T18:52:15Z No. of bitstreams: 1 tese_Suelen.pdf: 922374 bytes, checksum: f324427616027a422decc0eaf56c7ae2 (MD5) / Approved for entry into archive by Maria Cristina (library@lncc.br) on 2016-07-27T18:52:32Z (GMT) No. of bitstreams: 1 tese_Suelen.pdf: 922374 bytes, checksum: f324427616027a422decc0eaf56c7ae2 (MD5) / Made available in DSpace on 2016-07-27T18:52:42Z (GMT). No. of bitstreams: 1 tese_Suelen.pdf: 922374 bytes, checksum: f324427616027a422decc0eaf56c7ae2 (MD5) Previous issue date: 2016-04-15 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (Capes) / In this work a new method for obstacles reconstruction from partial boundary measurements is proposed. For a given boundary excitation, we want to determine the quantity, locations and sizes of a number of obstacles embedding whiting a geometrical domain, from partial boundary measurements related to such an excitation. This problem is written in the form of ill-posed and over-determinated partial differential equation. The idea therefore is to rewrite it as an optimization problem where a shape functional measuring the misfit between the boundary measurement and the solution to an auxiliary boundary value problem is minimized with respect to a set of ball-shaped holes. The topological derivative concept is used for solving the resulting topology optimization problem, leading to a second-order reconstruction algorithm free of initial guess. The resulting algorithm is non-iterative and thus very robust with respect to noisy data. Finally, some numerical results are presented in order to demonstrate the effectiveness of proposed reconstruction algorithm. / O objetivo deste trabalho é apresentar um novo método de reconstrução de obstáculos. Mais precisamente, dada uma excitação deseja-se obter a solução de um problema inverso de reconstrução consistindo na determinação da quantidade, localização e tamanho de obstáculos no interior de um dado domínio geométrico a partir de leituras parciais da resposta à referida excitação. Este problema é escrito na forma de uma equação diferencial parcial sobredeterminada. Essa dificuldade é contornada reescrevendo o problema inverso na forma de um problema de otimização. A ideia básica consiste em minimizar um funcional de forma que mede a diferença entre o dado lido e o calculado numericamente em relação ao próprio domínio geométrico. Em particular o conceito de derivada topológica é utilizado, o que conduz a um algoritmo de reconstrução de segunda ordem e independente de qualquer chute inicial. Como o algoritmo resultante é não-iterativo, o processo de reconstrução torna-se extremamente robusto à presença de ruído. Vários exemplos numéricos de reconstrução são apresentados donde se verifica a validade dos resultados obtidos.

Submitted by Maria Cristina (library@lncc.br) on 2017-08-14T18:34:35Z No. of bitstreams: 1 Welerson_Dissertação.pdf: 2360212 bytes, checksum: 58a44b5d4fff215888d80a0a417700de (MD5) / Approved for entry into archive by Maria Cristina (library@lncc.br) on 2017-08-14T18:34:46Z (GMT) No. of bitstreams: 1 Welerson_Dissertação.pdf: 2360212 bytes, checksum: 58a44b5d4fff215888d80a0a417700de (MD5) / Made available in DSpace on 2017-08-14T18:34:57Z (GMT). No. of bitstreams: 1 Welerson_Dissertação.pdf: 2360212 bytes, checksum: 58a44b5d4fff215888d80a0a417700de (MD5) Previous issue date: 2016-11-04 / In this work the optimal control problem with respect to a set of pointwise sources is studied. In particular, the control is given by a finite linear combination of Dirac mass and the state is solution to the associated elliptic boundary value problem. The basic idea consists in minimizing a functional which measures the distance between the state and a target function, with respect to the number, intensities and locations of pointwise loads. The sensitivity of the cost functional with respect to a number of pointwise sources in the set of admissible solutions is derived in its explicit form with help of auxiliaries boundary value problems. The obtained result is then used to devise a non-iterative second order reconstruction algorithm, independent of any initial guess and without introducing regularization techniques. Finally, the devised reconstruction algorithm is applied for numerically solving a set of control and inverse reconstruction problems. / Neste trabalho o problema de controle ótimo com respeito a um conjunto de fontes puntuais é estudado. Em particular, o controle é dado por uma combinação linear finita de massas de Dirac e o estado é solução de um problema de valor de contorno elíptico. Objetiva-se, portanto, minimizar um funcional, que mede a distância entre o estado e uma função alvo, com respeito ao número, intensidades e localizações das cargas puntuais. A sensibilidade do funcional de custo, em relação a um certo número de fontes puntuais no conjunto de soluções admissíveis, é analisada na sua forma explícita com o auxílio de problemas de valor de contorno auxiliares. O resultado obtido é então utilizado para conceber um algoritmo de reconstrução de segunda ordem não iterativo, independente de qualquer chute inicial e sem a introdução de técnicas de regularização. Finalmente, o algoritmo de reconstrução elaborado é aplicado para resolver numericamente um conjunto de problemas de controle e de problemas inversos de reconstrução de fontes.

A têmpera dos aços envolve a austenitização de uma peça seguida por um resfriamento rápido para promover a formação de microestrutura martensítica. É necessário avaliar os meios de têmpera para manter o processo de têmpera sob controle. Os parâmetros mais importantes no processo de resfriamento são o coeficiente de transferência de calor e/ou o fluxo de calor entre o meio de têmpera e a peça a ser resfriada. Um dos métodos de se avaliar os meios de têmpera (meios de resfriamento) e saber o que está acontecendo dentro da peça durante o resfriamento do ponto de vista térmico é o problema inverso de condução de calor. O problema inverso de condução de calor consiste na determinação de parâmetros como fluxo de calor, taxa de resfriamento e temperatura em qualquer posição através da peça, assim como o coeficiente de transferência de calor. Esses parâmetros são obtidos a partir de medições de temperatura em um ou mais pontos dentro da peça. O escopo deste trabalho foi desenvolver um software baseado no problema inverso condução de calor para avaliar meios de resfriamento para têmpera. A validação deste código foi feita usando água, óleo de soja, óleo mineral e solução aquosa de NaNO3. / Steels quenching involves part austenitization followed by a fast cooling to promote martensitic microstructure formation. It is necessary to evaluate quenchants in order to keep the quenching process under control. The most important cooling process parameters are the heat transfer coefficient and/or the heat flux between the quenchant and the part to be cooled. One of the methods to evaluate quenchants (cooling media) and to know what is happening inside the part during the cooling in the thermal point of view is the inverse heat conduction problem. The inverse heat conduction problem consists in the determination of parameters like heat flux, cooling rate and temperature in any position across the part, as well as the heat transfer coefficient. These parameters are obtained from temperature measurements in one or more points inside the part. The scope of this work was to develop a software based in the inverse heat conduction problem in order to evaluate quenchants for quenching. The validation of this code was made using water, soybean oil, mineral oil and NaNO3 aqueous solution.

[pt] A introdução de um objeto estranho ao corpo humano pode resultar de diferentes processos que abrangem desde eventos iatrogênicos oriundos de procedimentos cirúrgicos, até eventos traumáticos, ocasionados por acidentes ou violência, como disparos de arma de fogo e perfurações com objetos cortantes. A presença desses corpos estranhos dentro do organismo humano pode ocasionar problemas de saúde que vão desde dor e incômodos ao óbito. Para remoção cirúrgica, torna-se essencial a localização desses objetos, com elevada exatidão, para redução do tempo cirúrgico e garantia de sucesso do procedimento. O presente trabalho tem por objetivo ampliar o estudo do problema de localização de corpos estranhos magnéticos (intrinsicamente magnético ou por indução magnética externa) no organismo humano. Foram desenvolvidos algoritmos computacionais capazes de promover simulações dos padrões de densidade de fluxo magnético em um plano decorrentes de uma fonte extensa, representada por parâmetros modificáveis de comprimento, posição espacial e inclinação ao plano de medição. Essas simulações foram utilizadas como entrada para o treinamento de Redes Neurais Artificiais que, após treinamento, foram capazes de resolver o problema inverso, caracterizando, a partir do mapa de campo magnético, a posição espacial, tamanho (comprimento) e inclinação da fonte metálica. Os resultados obtidos indicaram melhor desempenho com o uso de sensores triaxiais, cujo erro quadrático médio, em 3640 testes, mostrou-se inferior a 2 mm na orientação espacial, a 8 mm no comprimento e a 17 graus Celsius para a inclinação da fonte magnética em relação ao plano de medição. / [en] The introduction of a foreign object into the human body may result from different processes ranging from iatrogenic events during surgical procedures to traumatic events caused by accidents or violence such as firing and piercing with sharp objects. The presence of these foreign bodies within the human body can cause health problems ranging from pain and discomfort to death. For surgical removal, it is essential to locate these objects, with high accuracy, to reduce surgical time and guarantee the success of the procedure. The objective of the present work is to contribute to the study of the localization problem of magnetic foreign bodies (intrinsically magnetic or by external magnetic induction) in the human body. The developed computational algorithms are capable of promoting simulations of magnetic flux density patterns in a plane from an extensive source, represented by modifiable parameters of length, spatial position and slope to the measurement plane. These simulations were used as input for the training of Artificial Neural Networks that, after training, were able to solve the inverse problem, characterizing, from the magnetic field map, the spatial position, size (length) and slope of the metallic source. The results indicated a better performance with the use of triaxial sensors, whose mean square error, in 3640 tests, was less than 2 mm in the spatial orientation, 8 mm in length and 17 Celsius degrees for the tilt of the magnetic source in relation to the measurement plane.

Neste projeto, propomos estudar a chamada \"Teoria de Corpos de Classe,\" que oferece uma descrição simples das extensões abelianas de corpos locais e globais, bem como algumas de suas aplicações, como os teoremas de Kronecker-Weber e Scholz-Reichardt / In this work, we study the so called \"Class Field Theory\", which give us a simple description of the abelian extension of local and global elds. We also study some applications, like the Kronecker-Weber and Scholz-Reichardt theorems

O objetivo deste trabalho é estudar métodos dinâmicos de detecção de dano em vigas, em especial os métodos baseados na variação da flexibilidade medida dinamicamente. Os métodos revisados formam parte das técnicas de Detecção de Dano Não Destrutivas (DDND). Nas técnicas DDND o dano é localizado por comparação entre o estado sadio e o danificado da estrutura. Neste trabalho, o problema de vibração inverso é apresentado e a matriz de flexibilidade estática da estrutura é determinada a partir de seus parâmetros modais. Com ajuda de um Modelo de Elementos Finitos (MEF) são mostrados os diferentes padrões de variação da matriz de flexibilidade produzidos pela presença do dano. Baseando-se nestes padrões é possível identificar a posição do dano dentro da estrutura, como indicado pelos diversos exemplos apresentados. / The purpose of this work is to study dynamic methods for damage detection in beam structures. The attention is devoted to the methods based on dynamically measured flexibility. The reviewed methods are part of Nondestructive Damage Detection techniques (NDD). In the NDD techniques the damage is determined through the comparison between the undamaged and damaged state of the structure. In this work the inverse vibration problem is presented and the structure\'s flexibility matrix calculated from his modal parameters. The Finite Elements Model (FEM) is employed to show that a clear pattern exist for the changes in the flexibility matrix produced due to the presence of damage. The flexibility matrix changes is used to identify and locate damage as indicated by the several examples presented.

Neste trabalho, foram utilizadas as vantagens da ferramenta matemática de análise temporal e espectral, a transformada wavelet discreta (DWT), além dos coeficientes de predição linear (LPC) e do algoritmo de inteligência artificial, Least Squares Support Vector Machines (LS-SVM), para aplicações em análise de sinais de voz e classificação de vozes patológicas. Inúmeros trabalhos na literatura têm demonstrado o grande interesse existente por ferramentas auxiliares ao diagnóstico de patologias da laringe. Os componentes da DWT forneceram parâmetros de medida para a análise e classificação das vozes patológicas, principalmente aquelas provenientes de pacientes com edema de Reinke e nódulo nas pregas vocais. O banco de dados com as vozes patológicas foi obtido do Departamento de Otorrinolaringologia e Cirurgia de Cabeça e Pescoço do Hospital das Clínicas da Faculdade de Medicina de Ribeirão Preto (FMRP-USP). Utilizando-se o algoritmo de reconhecimento de padrões, LS-SVM, mostrou-se que a combinação dos componentes da DWT de Daubechies com o filtro LP inverso levou a um classificador de bom desempenho alcançando mais de 90% de acerto na classificação das vozes patológicas. / The main objective of this work was to use the advantages of the time-frequency analysis mathematical tool, discrete wavelet transform (DWT), besides the linear prediction coefficients (LPC) and the artificial intelligence algorithm, Least Squares Support Vector Machines (LS-SVM), for applications in voice signal analysis and classification of pathological voices. A large number of works in the literature has been shown that there is a great interest for auxiliary tools to the diagnosis of laryngeal pathologies. DWT components gave measure parameters for the analysis and classification of pathological voices, mainly that ones from patients with Reinke\'s edema and nodule in the vocal folds. It was used a data bank with pathological voices from the Otolaryngology and the Head and Neck Surgery sector of the Clinical Hospital of the Faculty of Medicine at Ribeirão Preto, University of Sao Paulo (FMRP-USP), Brazil. Using the automatic learning algorithm applied in pattern recognition problems, LS-SVM, results have showed that the combination of Daubechies\' DWT components and inverse LP filter leads to a classifier with good performance reaching more than 90% of accuracy in the classification of the pathological voices.

A teoria de elasticidade linear clássica é utilizada no modelamento de problemas da física médica relacionados com a determinação de parâmetros elásticos de tecidos biológicos a partir da medição in vivo, ou, in vitro dos deslocamentos, ou, das deformações. Baseados em observações experimentais, as quais revelam que os tecidos biológicos anômalos têm comportamento mecânico diferente dos tecidos biológicos sadios, os pesquisadores têm modelado estes tecidos como sólidos elástico-lineares, isotrópicos, heterogêneos e incompressíveis. Neste trabalho, analisa-se uma classe de problemas planos relacionados à determinação do módulo de elasticidade ao cisalhamento µ de tecidos biológicos e propõe-se um procedimento numérico não-iterativo para obter soluções aproximadas para estes problemas a partir de campos de deslocamentos conhecidos de ensaios possíveis de serem realizados em laboratório. Os ensaios são estáticos e são simulados numericamente via método dos elementos finitos. Apresentam-se resultados obtidos das distribuições de µ em cilindros retos, longos e de seção retangular contendo inclusões cilíndricas circulares centradas, ou, excêntricas. Consideram-se inclusões mais, ou, menos rígidas do que o meio elástico circundante. Adicionalmente, os resultados obtidos no presente trabalho são comparados com resultados de outros pesquisadores que utilizam ensaios dinâmicos. Neste sentido, dois casos de inclusão circular centrada são resolvidos com as condições de contorno adaptadas do caso dinâmico para o caso estático. Resolve-se finalmente o caso de uma inclusão de forma geométrica complexa e seis vezes mais rígida do que o entorno. O cilindro contendo esta inclusão está submetido às condições de contorno propostas neste trabalho e também às condições de contorno adaptadas do caso dinâmico. Em todos os casos analisados os resultados são satisfatórios, apesar do emprego de uma quantidade reduzida de elementos finitos na reconstrução de µ. Deve-se ressaltar que nenhum método de regularização foi utilizado para tratar os deslocamentos obtidos dos ensaios simulados. Este trabalho é de grande interesse na detecção de tumores cancerígenos, tais como tumores nos seios e na próstata, e no diagnóstico diferenciado de tecidos biológicos. / The theory of classical linear elasticity is used to model of problems in medical physics that are related to the determination of elastic parameters of biological tissues from the measurement in vivo, or, in vitro of either displacements or strains. Based on experimental observations, which indicate that the abnormal biological tissues have different mechanical behavior from normal biological tissues, researchers have modeled these tissues as an incompressible, heterogeneous, and isotropic linear elastic solid. In this work a class of plane problems related to the determination of the shear elastic modulus µ of biological tissues is examined. A non-iterative numerical procedure to obtain an approximate solution to these problems from known displacement fields is proposed. The displacement fields are obtained from experiments that are possible to reproduce in laboratory. The experiments are quasi-static and are simulated numerically using the finite element method. Results for the distribution of µ in long, straight cylinders of rectangular cross-sections, containing either centered or eccentric circular inclusions that are more, or, less stiff than the surrounding elastic medium, are presented. Additionally, the results obtained in this study are compared with results of other researchers who use dynamical experiments. In this sense, two cases of centered circular inclusions are solved by using an adaptation of the dynamical case to the static case. Finally, the case of an inclusion with a complex geometry that is six times more rigid than the surrounding medium is solved. In all cases analyzed, the results are satisfactory, despite the fact that they were obtained with a reduced number of finite elements. It should be noted that no method of regularization has been used to treat the displacement data obtained from the simulated experiments. This work is of great interest in the detection of cancerous tumours, such as those in the breasts and in the prostate, and in the differential diagnosis of biological tissues.

Conselho Nacional de Desenvolvimento CientÃfico e TecnolÃgico / No sÃculo XXI a humanidade produziu mais novos dados (informaÃÃes) do que em toda sua histÃria. Entender a natureza dos diversos sistemas que geram essa abundÃncia de dados se tornou um dos grandes desafios desse sÃculo. Uma forma de analisar formalmente esses grandes bancos de dados à empregando a teoria da informaÃÃo desenvolvida por Claude Shannon. Essa teoria permite, usando o princÃpio da mÃxima entropia, encontrar as distribuiÃÃes de probabilidades que melhor descrevem os comportamentos coletivos desses sistemas. Nessa dissertaÃÃo, discutimos a possibilidade de usar modelos tipo Ising para descrever observaÃÃes de sistemas reais. Devido a suas limitaÃÃes, empregar o modelo de Ising implica em supor que os elementos que constituem o sistema real sà podem estar em dois estados, por exemplo ativo ou inativo. AlÃm disso, o modelo de Ising da conta apenas de interaÃÃes entre pares de elementos e desconsidera a possibilidade de interaÃÃes entre grupos maiores de elementos. Como discutiremos, mesmo com essas limitaÃÃes tal modelo pode descrever bem resultados observados em alguns sistemas naturais, como por exemplo redes de neurÃnios. Especificamente, discutiremos resultados de trabalhos anteriores que mostram que usando apenas as mÃdias de atividade de cada neurÃnio e a correlaÃÃo entre os mesmo, usando a teoria de Shannon, observa-se que os estados visitados pela rede seguem à distribuiÃÃo de Ising. Para testar a aplicabilidade desse mÃtodo em diversos sistemas geramos dados sintÃticos, obtidos de modelos tipo Ising em trÃs situaÃÃes: ferromagnÃtico, anti-ferro e vidro de spins (spin glass). NÃs chamamos o sistema que gera os dados sintÃticos de sistema subjacente. Usamos mÃtodos de maximizaÃÃo de entropia para tentar construir sistemas modelos que consigam reproduzir as mÃdia e correlaÃÃes observadas nos dados sintÃticos. Dessa forma, verificamos em que situaÃÃes nossos mÃtodos conseguem de fato gerar um sistema modelo que reproduza o sistema subjacente que gerou os dados. Esses resultados podem estabelecer um limite de aplicabilidade para a tÃcnica discutida.