Spelling suggestions: "subject:"(computacionais)""
31 |
GINGA - Graphical Interface for Comparative Genome Analysis: o desenvolvimento de um sistema computacional de visualização gráfica para a análise comparativa de genomas de bactérias / GINGA - Graphical Interface for comparative Genome Analysis: development of a computational system to visualize the comparative of bacterial genomes in a graphical viewPaschoal, Alexandre Rossi 23 March 2007 (has links)
Made available in DSpace on 2015-03-04T18:50:44Z (GMT). No. of bitstreams: 1
GINGA_Final_arp.pdf: 2617259 bytes, checksum: c988d46e191f2e44e3ce926beaa06fb0 (MD5)
Previous issue date: 2007-03-23 / This study aimed to develop a computational system applied to the comparative analysis of prokaryotic genomes in a graphical view. The system named GINGA Graphical Interface for comparative Genome Analysis was developed to analyse a draft genome sequence in comparison to a complete genome. The system shows the alignment between sequence of reads, contigs and scaffolds from partial sequenced genomes and the complete sequence of another genome and allows the identification shared and unique regions as well as rearrangements. GINGA is a web-based system developed using the PERL language to access a MySQL database where all the information regard to the comparative analysis is stored. The module of the interface to GD (Graphics Library) was used to help the construction of the graphical tool. The graphical view allows zoom in/out on the information on assembly, annotation and the organization of the sequences. Supplementary information can be accessed in the form of reports.
GINGA system was used to compare the genomes of Leifsonia xyli subsp. cynodontis (Lxc draft genome sequence) and Leifsonia xyli subsp. xyli (Lxx complete genome sequence). The mail goal was to identify genetic differences that may help to understand the pathogeniciy of Lxx towards sugarcane. A total of 9.754 reads assembled in 1.064 contigs and 317 scaffolds produced 1.470.731 of no redundant bases of Lxc genome and were used in the analysis. GINGA allowed the identification of 206.320 bp (~20%) of Lxc specific sequences organized in contigs and 56.884 bp organized in 19 scaffolds (5,9%), around 1 milion bp aligned to Lxx genome and at least 6 large scale genomic rearrangements. These results were presented in a graphical interface and allowed to guide the partial genome sequencing, helping to decide which regions should be further sequenced and at the same time allowing the formulation of hypothesis related to important biological aspects of these microorganisms / Esta dissertação resultou de um sistema computacional voltado para a visualização gráfica de análises comparativas entre genomas de procariotos. O sistema denominado de GINGA Graphical Interface for comparative Genome Analysis foi desenvolvido basicamente para analisar genomas parcialmente seqüenciados por meio da comparação com genomas completos. O sistema mostra a representação do alinhamento entre seqüências de reads, contigs e scaffolds do genoma parcial com a seqüência completa do outro genoma, permitindo a identificação de blocos comuns, regiões específicas e rearranjos. GINGA é um sistema web-based que foi desenvolvido em linguagem PERL para acessar um banco de dados MySQL, onde estão armazenadas as informações obtidas nas análises comparativas. O módulo de interface da biblioteca gráfica GD da linguagem PERL foi utilizado para a construção da ferramenta de visualização. A representação gráfica criada permite a navegação com opções de zoom in/out, disponibilizando as informações de montagem, anotação das seqüências codificadoras e da organização das seqüências entre os genomas. Relatórios são ainda disponibilizados como fonte complementar da apresentação dos resultados.
O sistema GINGA foi utilizado para analisar de maneira comparativa o genoma das bactérias Leifsonia xyli subsp. cynodontis (Lxc genoma parcialmente seqüenciado) e Leifsonia xyli subsp. xyli (Lxx genoma completamente seqüenciado). Lxx provoca o raquitismo da soqueria em cana-de-açúcar, enquanto Lxc é capaz de colonizar cana-de-açúcar sem provocar sintomas de doença. O objetivo foi revelar, ainda durante o processo de seqüenciamento do genoma de Lxc, diferenças genéticas existentes entre os genomas dessas duas bactérias. Fizeram parte das análises comparativas um total de 9.754 reads do genoma de Lxc que formaram 1.064 contigs e 317 scaffolds, totalizando 1.470.731 de bases não redundantes. GINGA permitiu a identificação de 206.320 bases (~19%) em seqüências de contigs específicos (contigs que não apresentaram alinhamento algum com o genoma completo de Lxx) e 19 scaffolds (5,9%) que totalizaram 56.884 bases específicas ao genoma de Lxc, além de aproximadamente 1 milhão de nucleotídeos alinhados ao genoma de Lxx e pelo menos 6 grandes rearranjos.
Estes resultados foram disponibilizados em uma interface gráfica e relatórios, permitindo orientar o andamento do projeto de seqüenciamento do genoma de Lxc quanto à seleção das regiões a serem seqüenciadas e, simultaneamente, oferecendo informações para a formalização de hipóteses relevantes à biologia destes microorganismos.
|
32 |
Uma filogenia mitocondrial de metazoários / A mitochondrial metazoan phylogenyGarcia, Marcelo 01 June 2007 (has links)
Made available in DSpace on 2015-03-04T18:50:50Z (GMT). No. of bitstreams: 1
Dissertacao Marcelo Garcia Jun-2007.pdf: 1765034 bytes, checksum: ccb95f7d696b4361f3d25db96fd2d70a (MD5)
Previous issue date: 2007-06-01 / Coordenacao de Aperfeicoamento de Pessoal de Nivel Superior / Inferring the evolutive relations between the animal phyla has been a formidable challenge to Science. The animal phyla represent quite distinct baupläne (body architectures) and are, therefore, difficult to compare. At the same time, their fossil record converges mostly to the same period on the geological scale. The recent availability of molecular data has, however, inaugurated a new front in animal phylogeny. The present work explores this opportunity by inferring a phylogeny with distance and maximum likelihood methods, employing all animal mitochondrial genomes ever sequenced. The results present only a few bigger groups with strong statistical support, like Diploblastica, Bilasteria, Protostomia and Deutorostomia, and many smaller groups of animals belonging to the same order or family. These results seems to confirm that the phyla radiated in such a short time interval that the phylogenetic signal did not hold out to produce a satisfactory resolution of the animal tree to date. Some limits may have yet to be tested, through models of evolution more fit to this scenario. For example was only recovered with the use of gamma distances for site-to-site substitution rate variability, at the expense of compressing the smaller branches throughout the tree. Nematodes and Platyhelminthes reveal a bias in GC and AT skew that cannot be adequately mapped by any reversible substitution pattern. Nevertheless, even if corrections are found for these issues, it is well possible that the hope of a better resolution in the animal tree will lie further on, by a better understanding of the evolutive process in a genomic scale. / Discernir as relações evolutivas entre os grandes grupos animais tem representado um formidável desafio para a Ciência. Os filos animais possuem arquiteturas corporais bastante distintas e por isso difíceis de serem comparadas. Ao mesmo tempo, seu registro fóssil converge aproximadamente para um mesmo intervalo na escala geológica, dificultando uma reconstrução filogenética com caracteres morfológicos. A disponibilidade de dados moleculares sobre os organismos abriu, contudo, novas possibilidades na filogenia animal. Esta dissertação buscou explorar essas possibilidades inferindo uma filogenia com métodos de distância e máxima verossimilhança, a partir de todos os genomas mitocondriais, completamente seqüenciados até o momento. No entanto, apenas alguns grandes agrupamentos como Diploblastica, Bilateria, Deuterostomia e Protostomia foram recuperados com forte suporte estatístico, além de pequenos agrupamentos de animais de mesma ordem ou família, indicando que os efeitos da rápida radiação no Cambriano se estenderam também ao registro molecular. Os resultados também indicam a necessidade de buscar modelos de evolução mais aderentes a este cenário. Deuterostomia, por exemplo, só foi recuperado monofileticamente assumindo-se a distribuição gama para variabilidade entre-sítios, ao custo, entretanto, da perda de definição nos ramos menores. Nematóides e Platelmintos, por sua vez, revelam um possível viés no skew (desvio) do conteúdo GC e AT de seus genes mitocondriais, que não é adequadamente mapeado por modelos de substituição reversíveis. Os indícios são de que a resolução da filogenia animal depende ainda de uma melhor compreensão da evolução molecular em escala genômica.
|
33 |
Modelagem Computacional de escoamento Bifásico em Meios Porosos Heterogêneos com Acoplamento Geomecânico / Computational Modelling of the Biphasic Flow in Heterogeneous Porous Media with Geomechanic CouplingMendes, Marcos Alcoforado 17 December 2007 (has links)
Made available in DSpace on 2015-03-04T18:50:55Z (GMT). No. of bitstreams: 1
TESE.pdf: 2353628 bytes, checksum: 5452725da7b8af8e0ba71da5e8c21c09 (MD5)
Previous issue date: 2007-12-17 / Conselho Nacional de Desenvolvimento Cientifico e Tecnologico / In this work we develop the computational modeling of the hydromechanical couplings which govern two-phase flow in a heterogeneous poroelastic media.
At the Darcy scale the governing equations are decomposed in two subsystems associated with the poromechanics and hydrodynamics. In this context new numerical methods are proposed for the computation of the Darcy velocity based on Petrov-Galerkin post processing techiniques in conjunction with locally conservative methods for the hyperbolic transport equation and for water saturation together with an operator splitting technique for the computation of the transient porosity effect upon the transport equation.
Numerical simulations allow us to identify different regimes of hydromechanics coupling during the secondary of petroleum withdrawal.In particular, we analyze the influence of the viscosity ratio and strength of heterogeneity upon the various coupling regimes. Among the many phenomena captured by the model we give particular emphasis on the delayed compaction of the reservoir due to the water flooding illustrating its effects upon the oil production curves.
Within the framework of the stochastic modelinf, we analyze the effects of the heterogeneity and the uncertainty in permeability and elastic coeddicientsupon the hydromechanical coupling regimes. The characterization of these regimes governed by yhe viscosity ratio and by the strength of heterogeneity is illustrated in numerical simulations of secondary oil withdrawal from a reservoir subject to the overburden due the weight of the overlaying formations. / Neste trabalho desenvolvemos a modelagem computacional dos fenômenos inerentes ao acoplamento hidromecânico que governam o escoamento de dois fluidos imiscíveis em uma matriz porosa heterogênea e deformável. As equações do modelo na escala de Darcy são decompostas em dois subsistemas associados a poromecânica e ao transporte dos fluidos. Neste contexto novos métodos numéricos são propostos para a computação da velocidade de Darcy baseados em técnicas de pós-processamento de Petrov-Galerkin em conjunção com métodos localmente conservativos para a equação hiperbólica de transporte da saturação aliados à técnica de decomposição de operadores para a computação dos efeitos da evolução temporal da porosidade induzidos pela deformação da matriz porosa sobre o transporte.
As simulações numéricas do modelo resultante nos permite identificar diferentes regimes do acoplamento hidromecânico durante o processo de extração secundária de petróleo. Em particular analisamos a influência da razão de viscosidade entre os fluidos e da heterogeneidade da matriz porosa sobre os diferentes regimes de acoplamento. Dentre os vários fenômenos capturados no modelo damos particular ênfase ao surgimento de um processo de compactação retardada do reservatório devido à inundação de água onde, em particular, ilustramos seu efeito sobre as curvas de produção de petróleo.
No contexto da modelagem estocástica analisamos também os efeitos das heterogeneidades e incertezas presentes nos coeficientes de permeabilidade e das constantes elásticas do meio poroso sobre os diferentes regimes do acoplamento geomecânico. A caracterização destes regimes governados pela razão de viscosidade e pelo coeficiente de variação das heterogeneidades é ilustrada por meio de simulações numéricas do processo de extração secundária de um reservatório sujeito ao peso das camadas superiores.
|
34 |
Método de Elementos Finitos Enriquecidos para uma Classe de Problemas Elípticos Não Lineares com Coeficientes Altamente OscilatóriosBarreda, Manuel Jesus Cruz 22 July 2009 (has links)
Made available in DSpace on 2015-03-04T18:51:13Z (GMT). No. of bitstreams: 1
thesis.pdf: 859231 bytes, checksum: 157dc07ef14060b96e2fca91dec344cb (MD5)
Previous issue date: 2010-07-22 / Fenômenos em materiais heterogêneos conduzem ao estudo de problemas em equações diferenciais parciais com coeficientes altamente oscilatórios. O tratamento numérico mediante o uso dos métodos tradicionais exige um alto custo computacional
ou é inviável. No presente trabalho pretendemos estender o método residual free bubbles com o intuito de gerar um procedimento de homegeneização numérica para o estudo de
uma classe de problemas elípticos não lineares com coeficientes que têm um comportamento altamente variável (problemas
multiescala).Mostramos que a formulação numérica decorrente da
metodologia residual free bubbles permite aproximar o problema
multiescala e, portanto, o problema efetivo. Para validar o procedimento proposto, apresentaremos estimativas de erro e resultados numéricos.
|
35 |
Evolução diferencial para problemas de otimização restrita / Differential evolution for constrained optimization problemsSilva, Eduardo Krempser da 04 March 2009 (has links)
Made available in DSpace on 2015-03-04T18:51:14Z (GMT). No. of bitstreams: 1
dissertacao_eduardo_krempser.pdf: 675502 bytes, checksum: 9909466baacd7f6cebcd6029c53ed6c4 (MD5)
Previous issue date: 2009-03-04 / Conselho Nacional de Desenvolvimento Cientifico e Tecnologico / Optimization is a large area of knowledge concerned with the need of a better use of resources and activities, becoming indispensable in the solution of several problems which arise from the study and formulation of real-world problems. Furthermore, the constraints that must be respected for each situation introduce in the methodologies of optimization an additional complication. Differential Evolution, which in its original formulation is applied only to unconstrained optimization problems in continuous space, also provides good results when applied to constrained optimization with discrete and continuous variables. This work presents the necessary improvements to Differential Evolution for its proper application to this class of problems, and proposes a new combination of techniques for this application, as well as a mechanism for dynamic selection of the appropriate variant of the technique. The initial proposal is a combination of Differential Evolution with a technique of adaptive penalty (APM) and the second proposal concerns the dynamic selection of variants during the search process. Several computational experiments are carried out confirming the competitiveness of the proposed algorithms. / A otimização é uma grande área de conhecimento voltada para a necessidade de um melhor aproveitamento de recursos e atividades, tornando-se indispensável na resolução de grande parte dos problemas oriundos de estudos e formulações de problemas reais. Além disso, as restrições que devem ser respeitadas para cada situação introduzem nas metodologias de otimização um complicador adicional. A Evolução Diferencial, que em sua formulação original é aplicada somente a problemas de otimização irrestrita e em espaços contínuos, apresenta também bons resultados quando aplicada à otimização restrita com variáveis contínuas e discretas. Este trabalho apresenta os aperfeiçoamentos necessários à Evolução Diferencial para sua adequada aplicação sobre essa classe de problemas, além de propor uma nova combinação de técnicas para essa aplicação, bem como um mecanismo de seleção dinâmica da variante adequada da técnica. A proposta inicial é a combinação da Evolução Diferencial com uma técnica adaptativa de penalização (APM) e a segunda proposta visa a seleção dinâmica de variantes durante o processo de busca. Vários experimentos computacionais são executados confirmando a competitividade dos algoritmos propostos.
|
36 |
Algoritmos genéticos assistidos por metamodelos baseados em similaridade / Genetic algorithms assisted by similarity-based metamodelsFonseca, Leonardo Goliatt da 24 June 2010 (has links)
Made available in DSpace on 2015-03-04T18:51:15Z (GMT). No. of bitstreams: 1
Goliatt tese.pdf: 3293381 bytes, checksum: e89a55efe1733bca4d491ccb2aca5af1 (MD5)
Previous issue date: 2010-06-24 / Coordenacao de Aperfeicoamento de Pessoal de Nivel Superior / Vários problemas de interesse em Ciência e Engenharia são formulados como problemas de otimização.
A complexidade dos problemas modernos tem levado ao desenvolvimento de modelos matemáticos de complexidade crescente, resultando em modelos de simulação computacionalmente custosos.
Algoritmos Genéticos (AG), inspirados na Teoria de Evolução por seleção natural, são ferramentas versáteis em problemas difíceis de busca e otimização.
Entretanto, eles usualmente requerem um elevado número de
avaliações até a obtenção de uma solução viável ou satisfatória. Em um cenário de simulações dispendiosas, o
uso de Algoritmos Genéticos pode tornar-se proibitivo.
Uma possível solução para este problema é o uso de um
metamodelo, para ser usado no processo de otimização no lugar do modelo de simulação.
Nesta tese desenvolveu-se uma metodologia para o uso combinado de AG e metamodelos para otimização mono- e
multi-objetivo de alto custo computacional, onde
metamodelos baseados em similaridade são
incorporados nos AG com o objetivo de melhorar o
seu desempenho.
A metodologia foi aplicada em problemas de otimização
coletados da literatura, e em problemas de Otimização Estrutural, demonstrando sua aplicabilidade e estabelecendo esta como uma alternativa para o melhoramento de soluções em um contexto de orçamento
fixo de simulações.
|
37 |
Um esquema de assimilação de dados oceanográficos para o modelo oceânico HYCOM ao largo da costa sudeste brasileira / A data assimilation scheme using the ocean model HYCOM for southeastern brazilian bightOliveira, Jean Felix de 22 December 2009 (has links)
Made available in DSpace on 2015-03-04T18:51:15Z (GMT). No. of bitstreams: 1
thesis.pdf: 19996358 bytes, checksum: b3a11077536e0bcd42efeade2b4a5bed (MD5)
Previous issue date: 2009-12-22 / The present work presents a data assimilation scheme customized to work with the Hybrid Coordinate Ocean Model (HYCOM) for the Southeastern Brazilian Bights. HYCOM uses hybrid vertical coordinates, i.e., it uses z coordinates in the mixed layer, isopycnal coordinates in the deep ocean and sigma-z coordinates in the continental shelf. However, since vertical profiles of the main ocean variables, like temperature, density and salinity, are observed in z -coordinates, the assimilation of these data into HYCOM is not trivial. For this reason, a technique to transform vertical profiles from isopycnal coordinates to z -coordinates is here proposed as an alternative to realize data assimilation in HYCOM. This technique uses Lagrangian multipliers with a optmization process that guarantees the conservation of the barotropic mass ux.
The technique of transformation is applied with the data assimilation method proposed by Ezer & Mellor (1997). The method uses statistical interpolation and correlations, a priori calculated with the model´s output, between the sea surface data - temperature (SST) and/or height (SSH) - and subsurface potential temperature and density structures.
Numerical experiments showed that the data assimilation scheme is able to reproduce eficiently the local ocean circulation. The best performance scheme included the correlation with both SST and SSH. / Neste trabalho é apresentado um esquema de assimilação de dados a ser realizado com o Modelo Oceânico de Coordenadas Híbridas HYCOM ao largo da costa sudeste brasileira. O HYCOM utiliza 3 diferentes coordenadas verticais, a saber: coordenada-z na camada de mistura, coordenada isopicnal no oceano profundo estratificado e coordenada sigma-z nas regiões mais rasas e costeiras. Entretanto, como os perfis verticais das principais variáveis oceânicas, como temperatura, salinidade e densidade, são observados e disponibilizados em coordenadas-z, a assimilação desses dados não é tão trivial. Por esse motivo, uma técnica de transformação de coordenadas verticais de isopicnal para z é aqui proposta como uma alternativa para a realização da assimilação de dados no HYCOM. Essa técnica utiliza multiplicadores de Lagrange juntamente com um processo de otimização que garante a conservação do fluxo de massa barotrópico. A técnica de transformação é aplicada juntamente com o método de assimilação de dados proposto por Ezer & Mellor (1997). Esse método utiliza interpolação estatística e correlações, calculadas a priori com resultados do modelo, entre dados de superfície - temperatura (TSM) e /ou altura (ASM) - e a estrutura de subsuperfície de temperatura e densidade potenciais.
Com base nos experimentos numéricos realizados, pode-se verificar que o esquema de assimilação de dados foi capaz de reproduzir eficientemente a circulação oceânica do domínio proposto e com os melhores resultados quando utilizando conjuntamente ASM e TSM nas correlações.
|
38 |
Métodos não lineares descontínuos submalha para a equação de convecção-difusão-reaçãoSilva, Natalia Cristina Braga Arruda Alves da 16 July 2010 (has links)
Made available in DSpace on 2015-03-04T18:51:19Z (GMT). No. of bitstreams: 1
TeseNatalia.pdf: 2041250 bytes, checksum: 574e8fdc58499aa0fe74fc2587a2a8d6 (MD5)
Previous issue date: 2010-07-16 / Coordenacao de Aperfeicoamento de Pessoal de Nivel Superior / Neste trabalho três novas formulações de Galekin descontínuo são desenvolvidas utilizando a estrutura submalha de duas escalas para resolver as equações de convecção-difusão-reação. A primeira, inspirada no método NSGS (Nonlinear Subgrid Stabilization) introduz um termo não-linear de viscosidade submalha agindo apenas nas escalas menores da malha de elementos finitos. A segunda, baseada nas arestas, introduz dois termos de estabilização submalha considerando o resíduo da solução aproximada na escala resolvida em cada elemento, e o salto da solução submalha nas arestas entre os elementos. A terceira, acrescenta viscosidade artificial que atua em todas as escalas, que também é controlada pela solução na escala resolvida a nível de elemento. Todas as formulações podem ser consideradas adaptativas, no sentido de que a estabilização atua somente onde é necessária. Para as duas primeiras formulações as estimativas de erro revelam taxas de convergência semelhantes `as de métodos estabilizados, embora se obtenha taxas ótimas para os três métodos em problemas regulares. Experimentos numéricos são realizados a fim de demonstrar o comportamento dos novos métodos em comparação com outros métodos encontrados na literatura
|
39 |
Uso de conceitos arquiteturais para aumento da expressividade na modelagem de workflows cientificos / Use of architectural concepts for improved expressiveness in the modeling of scientific workflowsMedeiros, Vivian 10 May 2013 (has links)
Made available in DSpace on 2015-03-04T18:57:51Z (GMT). No. of bitstreams: 1
thesis_vivian_final.pdf: 2459870 bytes, checksum: 86a3974bde419685b151ae0d526f328d (MD5)
Previous issue date: 2013-05-10 / In-silico experiments frequently take the form of workflows since they involve the execution of sets of interdependent tasks.The automation of these workflows through Scientific Workflow Management Systems (SWMSes) has been the subject of several papers in the area of e-Science.Traditionally, SWMSes provide fundamental automation mechanisms such as task dispatching and data movement, but recent studies have also shown interest in supporting non-functional attributes such as reliability and traceability in these systems.Nevertheless, the configuration of these attributes in current practice is limited in several aspects, particularly with respect to the interactions between tasks.Such characteristic leads to limited expressiveness in the workflow modeling languages provided by these systems. Aiming at improving such expressiveness, this work presents a proposal called Open Scientific Connectors (OSC), which adopts architectural concepts for the representation of interactions between tasks as first-class modeling elements. This approach provides better composition and reuse in workflow modeling, particularly in the configuration of mechanisms that deal with non-functional attributes. The proposal is evaluated through proofs of concept , in which the expressiveness, the reusability, the compositionality and the simplicity of OSC is compared with related work. A prototype of SWMS based on OSC also illustrates the feasibility of workflow models described in this language. / Experimentos in-silico não raramente tomam a forma de workflows por envolverem a execução de conjuntos de tarefas interdependentes. A automação desses workflows através de Sistemas Gerenciadores de Workflows Científicos (SGWfCs) tem sido alvo de vários trabalhos na área de e-ciência. Geralmente, os SGWfCs oferecem mecanismos fundamentais de automação como despacho de tarefas e movimentação de dados, mas trabalhos recentes têm demonstrado também o interesse no suporte a atributos não-funcionais, como confiabilidade e rastreabilidade, nesses sistemas.
A configuração desses atributos nos SGWfCs atuais é, contudo, limitada em vários aspectos, em particular no que se refere às interações entre tarefas. Essa limitação confere menor expressividade às linguagens de modelagem de workflows providas por esses SGWfCs. Objetivando o aprimoramento dessa expressividade, este trabalho apresenta a proposta Open Scientific Connectors (OSC), que usa conceitos arquiteturais para representar interações entre tarefas como elementos de primeira classe. Esta abordagem propicia uma maior capacidade de composição e de reuso na modelagem de workflows, particularmente nas configurações dos mecanismos que lidam com atributos não-funcionais. A proposta é avaliada por provas de conceito, onde aspectos como a expressividade, a capacidade de composição e de reuso e a simplicidade de OSC são comparados a trabalhos relacionados. Um protótipo de um SGWfC ilustra ainda a exequibilidade dos modelos descritos em OSC.
|
40 |
Integrative computational modeling & in-vivo characterization of residual deformations in hemodynamicsDamián Ares, Gonzalo 11 April 2016 (has links)
Submitted by Maria Cristina (library@lncc.br) on 2016-07-28T15:01:01Z
No. of bitstreams: 1
ThesisGAresFrente.pdf: 14194206 bytes, checksum: 2e35ae71aaffd676ba8015d68298aca6 (MD5) / Approved for entry into archive by Maria Cristina (library@lncc.br) on 2016-07-28T15:01:19Z (GMT) No. of bitstreams: 1
ThesisGAresFrente.pdf: 14194206 bytes, checksum: 2e35ae71aaffd676ba8015d68298aca6 (MD5) / Made available in DSpace on 2016-07-28T15:01:31Z (GMT). No. of bitstreams: 1
ThesisGAresFrente.pdf: 14194206 bytes, checksum: 2e35ae71aaffd676ba8015d68298aca6 (MD5)
Previous issue date: 2016-04-11 / Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ) / This thesis is concerned with two major problems arising in the modeling of the cardiovascular system. The first topic consists in a comprehensive approach for the simulation of arterial blood flow and its effect on the stress state of the arterial wall, and the second topic is concerned with the in-vivo characterization of residual deformations in arterial wall tissues, based on data provided by medical images.
Specifically, regarding the first topic, an original modeling framework is proposed for the treatment of hemodynamic problems with increased realism, featuring a combination of several modeling techniques in order to account for i) the fact that the initial (image-based) geometry corresponds to a configuration which is at equilibrium with an internal pressure acting over the lumen, and with tethering forces located at the artificial (axial) boundaries delimiting the arterial region of interest; ii) the fluid-structure interaction problem; iii) the complex constitutive behavior of the arterial wall; iv) the influence of surrounding tissues; v) the interaction of the vessel with the rest of the cardiovascular system; and iv) the influence of residual stresses.
In order to tackle the issues described above, the preload mechanical problem is solved in a first stage, finding the zero-load material configuration which is employed to define suitable constitutive equations. This is performed by finding the solution for the mechanical equilibrium of the given image configuration considering the vessel at this state to be loaded by an internal baseline pressure and an axial traction (caused by tethering forces) at the artificial boundaries. It is worthwhile to mention that this axial traction is such that a previously defined pre-stretch level is considered on the equilibrium image configuration.
Once the reference configuration is obtained, the complete 3D fluid-structure interaction simulation is carried out, coupled with a dimensionally reduced 1D model of the rest of the cardiovascular system. Strong coupling via fixed-point iterations is achieved for the fluid-structure interaction, while the dimensionally heterogeneous coupling is achieved through a Broyden method. Regarding the constitutive modeling, a fiber-reinforced hyperelastic constitutive law is considered. Furthermore, through the analysis of several numerical examples, the sensitivity with respect to the existence of the preload stresses is assessed to quantify the importance of this issue. These results indicate that the stress state of the arterial wall is strongly influenced by the existence of preload. Therefore, the consideration of such preload state is mandatory for the prediction of stresses in arterial tissue.
For the second topic, a conceptual framework is presented for the in-vivo estimation of residual deformations and stresses. As a given data, a set of known configurations for an arterial segment is considered, which can potentially be obtained from medical imaging techniques. The mechanical equilibrium equations corresponding to such configurations are introduced through a variational approach, highlighting the role of the residual deformations and associated stresses. In this context, a cost functional is proposed to measure the imbalance of the mechanical setting arising from the consideration of inconsistent residual deformations, based on the generalized residuals of the associated variational equations. Then, the characterization of residual deformations becomes an optimization problem, focused on the minimization of this cost functional. For this purpose, a simple gradient descent method and an interior-point algorithm for constrained optimization are explored in this work.
The proposed methodology is tested using three numerical examples based on manufactured solutions, a simple clamped bar, a thick-walled cylinder and a three-layered aorta artery. The obtained results are promising and suggest that the present method (or variants based on the present ideas), when coupled with adequate image acquisition techniques, could successfully lead to the in-vivo identification of residual deformations. / Esta tese aborda dois problemas de relevância na modelagem do sistema cardiovascular humano. O primeiro tema consiste no desenvolvimento de um enfoque abrangente para a simulação do escoamento sanguíneo e sua interação com a parede arterial, e o segundo tópico é a caracterização in-vivo de tensões e deformações residuais na parede arterial baseada em dados fornecidos por imagens médicas.
De maneira específica, em relação ao primeiro tópico, um marco de modelagem é proposto para o tratamento de problemas hemodinâmicos com um alto grau de realismo, apresentando uma combinação de diferentes técnicas de modelagem para levar em conta i) o fato que as geometrias iniciais obtidas a partir de imagens médicas são correspondentes a um sistema de carregamentos não nulos, definido pela existência da pressão interna no lumen e de tensões axiais localizadas nos contornos artificiais do segmento arterial; ii) o problema de interação fluido-estrutura; iii) o complexo comportamento constitutivo da parede arterial; iv) a interação do segmento de interesse com o resto do sistema cardiovascular; e v) a influência dos tecidos circundantes; e vi) a existência de tensões residuais.
Para a abordagem das questões descritas acima, o problema mecânico de precarregamento é resolvido em uma primeira etapa, encontrando a configuração material de carregamento nulo onde as equações constitutivas são usualmente definidas. Isto é realizado encontrando a solução do problema de equilíbrio mecânico da estrutura arterial dada, considerando que o vaso está submetido a um nível de pressão de base e uma tração axial nos contornos artificiais. Vale a pena ressaltar que esta tração axial é correspondente a um nível de pre-estiramento previamente definido.
Uma vez que a configuração de referência é obtida, a simulação fluido-estrutura 3D é realizada, acoplada com um modelo dimensionalmente reduzido do resto do sistema cardiovascular. Um acoplamento forte através de iterações de ponto fixo é empregado para representar a interação fluido-estrutura, equanto o acoplamento entre modelos dimensionalmente heterogêneos é conseguido usando um método tipo Broyden. Em relação à modelagem constitutiva, um modelo hyperelástico reforçado com fibras é considerado. Além disso, através da análise de vários exemplos numéricos, a sensibilidade com relação à existência de precarregamentos é quantificada para remarcar a relevância desta questão. Tais resultados indicam que o estado de tensão da parede arterial é fortemente influenciado pela existência de precarregamentos. Assim sendo, levar em consideração esse estado de precarga é fundamental para a predição de tensões no tecido arterial.
Em relação ao segundo tópico, um marco conceptual é apresentado para estimação de tensões e deformações residuais. Consideramos que os dados são um conjunto de configurações de um segmento arterial, as quais poderiam ser obtidas a partir do uso de técnicas de adquisição e , processamento e segmentação de imagens. Utilizando um enfoque variacional, são apresentadas as equações de equilíbrio mecânico para as configurações conhecidas, acentuando o papel desempenhado pelas deformações residuais. Neste contexto, apresenta-se um funcional custo que mede o desbalance mecânico que é originado se um campo de deformações residuais inconsistente é admitido. Este funcional custo está baseado no resíduo generalizado das equações variacionais previamente mencionadas. Como consequência, o problema de estimação de deformações residuais é transformado em um problema de otimização, no qual se procura minimizar o funcional custo proposto. Com este objetivo, neste trabalho de tese são considerados dois métodos, um método de gradiente e um algoritmo de ponto interior para problemas que apresentam restrições.
A metodologia proposta é testada em três exemplos numéricos baseados em soluções manufaturadas: um barra engastada, um cilindro de parede grossa, e uma artéria aorta composta por três camadas. Os resultados obtidos são promissores e sugerem que o método apresentado (ou variantes baseadas nas ideias aqui mostradas) junto com técnicas adequadas para a adquisição de imagens podem conduzir à identificação in-vivo de deformações residuais.
|
Page generated in 0.0558 seconds