Apresentamos cálculos de estrutura eletrônica de oligotiofenos e oligoanilinas conectados a contatos metálicos. Diversos aglomerados de ouro e acoplamentos com as moléculas orgânicas foram testados. Expressões analíticas para a função de transmissão, de acordo com a teoria de Landauer, foram obtidas por meio de um modelo tight-binding, juntamente com resultados numéricos da teoria do funcional da densidade (DFT). O pacote computacional Gaussian 03 foi empregado para realizar os c_alculos DFT com o funcional híbrido B3LYP. Um conjunto de base misto foi escolhido: LANL2DZ para os átomo de ouro e 6-31G* para os átomos leves S, C, N e H. A geometria das cadeias orgânicas conectadas a eletrodos metálicos depende fortemente do tipo de conexão e fracamente do número de átomos representando o eletrodo. Em nossos modelos o átomo de enxofre conecta-se ao eletrodo nas geometrias bridge e atop. A principal diferença entre estes acoplamentos é uma significativa transferência de elétrons nos oligômeros ligados na forma atop. A redistribuição de carga resulta numa cadeia orgânica carregada positivamente acompanhada por uma mudança estrutural. Oligotiofenos apresentam uma distorção quinoide, enquanto que nas pernigranilinas(PE) sua alternância entre anéis quinoides e aromáticos é invertida. Um campo elétrico externo foi aplicado nos sistemas. A resposta da transmissão eletrônica _e correlacionada com o estado de oxidação das cadeias e acoplamento com os eletrodos. Nos sistemas de leucoesmeraldina (LE) e esmeraldina(EM) conectados de forma atop, o campo elétrico induz a localização dos orbitais de fronteira levando à diminuição da transmitância. Cadeias de LE e EM conectadas de forma bridge e as PE apresentam um aumento na transmissão em função do campo elétrico. Nas cadeias de PE as propriedades são menos dependentes da conexão com o eletrodo. Aplicamos uma tensão nas cadeias orgânicas forçando um alongamento, similarmente ao que acontece nos experimentos de quebra de junção. Cadeias de LE e EM conectadas de forma bridge bem como PE em ambas conexões com os eletrodos, apresentam a mesma dependência do transporte em função do estiramento. O aumento inicial da transmissão é associado à planarização da cadeia, que induz a deslocalização dos orbitais. Entretanto, o subsequente alongamento localiza os orbitais por aumentar os comprimentos de ligação, o que leva a uma baixa transmitância. Cadeias de LE e EM conectadas de forma atop são fracamente afetadas pelo alongamento. Em oligotiofenos nos dois tipos de conexão com os eletrodos, o estiramento provoca uma significativa redução da transmissão. Por este mecanismo, oligômeros mais longos apresentam valores de condutância similares aos mais curtos. / We present electronic structure and transport calculations of oligothiophenes and oligoanilines bonded to metallic contacts. Several gold clusters and couplings to the organic molecule were tested. Analytical expressions for the transmission function, within Landauer theory, were obtained by means of a tight-binding model, along with numerical results from density functional theory (DFT). The software Gaussian 03 was employed to perform DFT calculations with the hybrid functional B3LYP. A mixed basis set was chosen: LANL2DZ for gold atoms and 6-31G* for light atoms S, C, N and H. The geometry of organic chains bonded to metallic electrodes depend strongly on the type of connection with the metallic electrode and weakly on the number of atoms representing the electrode. In our models the sulfur atom connects to the electrodes in either atop or bridge geometries. The main diference between these two couplings is a significant electronic transfer on the atop bonded oligomers. The charge reorganization results in a positively charged organic part accompained by a structural change. Oligothiophenes present a quinoidal distortion, while in pernigranilines the alternation between aromatic and quinoid ring is reversed. An external electric field was applied to the systems. The electronic transmission response is correlated to the oxidation state and coupling to electrodes. On the atop bonded systems leucoemeraldine(LE) and emeraldine (EM), the electric field induces a localization of the frontier orbitals leading to a decreased transmission. Pernigranilines(PE) and bridge-bonded LE and EM present an increase of transmission as a function of the electric field. In the PE chains the properties are less dependent on the connection to the electrodes. We applied a tension on the organic chains by forcing an elongation, similarly to what happens in break-junction experiments. Bridge-bonded LE and EM, as well as PE in both types of connection to electrodes, present the same dependence of transport as function of stress. The initial increasing transmission is associated to a planarization of the chain that induces an orbital delocalization. However, the subsequent elongation localize orbitals by increasing the bond lengths and leads to a low transmission. Atop bonded LE and EM are weakly sensitive to elongation. On oligothiophenes in both types of connection to the electrodes, the stretching induces a significant decrease of transmission. By this mechanism longer oligomers present conductance values similar to small ones.

Made available in DSpace on 2015-03-05T13:56:58Z (GMT). No. of bitstreams: 0 Previous issue date: 9 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Este trabalho propõe uma metodologia para simular o transporte de poluentes leves sob a ação de campos de ventos em ambientes complexos. O método consiste na divisão da clássica equação de difusão advectiva em duas componentes, uma laminar e outra turbulenta. A primeira representa o campo de ventos, determinado a partir de poucas amostras devidamente coletadas na fronteira do ambiente. A estimativa dos campos de ventos resulta da interpolação destas amostras, para a determinação das condições de contorno do ambiente de simulação, e da solução da equação da continuidade para escoamentos potenciais. A segunda componente, que representa a parte turbulenta do escoamento, consiste em uma grade de valores nãohomogêneos para os coeficientes de difusão sobre o ambiente, determinada a partir de um modelo típico da difusão turbulenta, com base na fenomenologia da mecânica dos fluidos. Esta metodologia evita a necessidade da resolução das equações de Navier-Stokes, e permite maior eficiência na simulação. O método p / This work proposes a method for simulating the transport of pollutant, driven by wind fields in complex environments. The method consists in dividing the advective-difusion equation into two components, laminar and turbulent. The former represents the wind fields, obtained from a few samples properly collected at the environment boundary. The wind field estimation comes from the interpolation of those samples, for determining the boundary conditions of the simulation environment, and from numerical solution of the continuity equation for potential flows. The second component, representing the turbulent part of the flow, consists of a non-homogeneous scalar field of difusion coefcients defined in the ambient, in accordance with a typical model for the turbulent difusion, based on the fluid mechanics phenomenology. The proposed methodology avoids the need of solving the Navier-Stokes equations and gives the simulation more performance. The model, besides being computationally cheap, gives good precision on th

Submitted by Cátia Araújo (catia.araujo@unipampa.edu.br) on 2017-01-25T11:32:18Z No. of bitstreams: 2 license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) Um modelo para dispersão de poluentes na camada limite planetária com coeficientes de difusão dependentes da distância da fonte.pdf: 1813716 bytes, checksum: 2b7788e4583766741e3b5c4e7ae8da30 (MD5) / Approved for entry into archive by Cátia Araújo (catia.araujo@unipampa.edu.br) on 2017-01-25T11:33:16Z (GMT) No. of bitstreams: 2 license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) Um modelo para dispersão de poluentes na camada limite planetária com coeficientes de difusão dependentes da distância da fonte.pdf: 1813716 bytes, checksum: 2b7788e4583766741e3b5c4e7ae8da30 (MD5) / Made available in DSpace on 2017-01-25T11:33:16Z (GMT). No. of bitstreams: 2 license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) Um modelo para dispersão de poluentes na camada limite planetária com coeficientes de difusão dependentes da distância da fonte.pdf: 1813716 bytes, checksum: 2b7788e4583766741e3b5c4e7ae8da30 (MD5) Previous issue date: 2013-07-15 / Este trabalho apresenta a solução da equação da difusão-advecção bidimensional estacionária para simular a dispersão de poluentes na Camada Limite Planetária. A solução é obtida através do método ADMM (Analytical Dispersion Multilayer Model) e da técnica de inversão numérica utilizando o algoritmo de Fixed Talbot. A validação da solução é comprovada, mediante os parâmetros estatísticos, através do confrontamento das concentrações calculadas a partir do modelo com as obtidas experimentalmente pelo experimento de Prairie Grass. Para a determinação das concentrações utiliza-se o perfil do vento segundo o modelo de similaridade de Monin-Obukhov e os parâmetros de turbulência com dependência da distância longitudinal da fonte e da altura vertical, considerando a componente vertical do espectro Euleriano e de acordo com o modelo sugerido por Hɸjstrup que divide os espectros em alta e baixa frequência. Para efeito comparativo utiliza-se um coeficiente de difusão para grandes tempos de difusão. Os melhores resultados foram alcançados com a utilização dos coeficientes de difusão considerando a distância longitudinal da fonte e a altura vertical. / This work presents the solution of two-dimensional advection-diffusion equation stationary to simulate the dispersion of pollutants in the Planetary Boundary Layer. The solution is obtained through the ADMM method (Analytical Multilayer Dispersion Model) and the numerical inversion technique using the algorithm Fixed Talbot. Validation of the solution is proven, statistical parameters, through the confrontation of the concentrations calculated from the model with those obtained experimentally by the Experiment of Prairie Grass. For the determination of the concentration profile of the wind the form of Monin-Obukhov similarity and turbulence parameters with longitudinal distance dependence of source and of vertical height, considering the vertical component of the Eulerian spectrum and according to the model proposed by Hɸjstrup that divides the high and low frequency spectra. To use the comparative effect diffusion coefficient for large diffusion times according. The best results were achieved with the use of diffusion coefficients considering the longitudinal distance from the source and the vertical height.

Three methods were developed to better understand and characterize the near-field dynamic processes of rotary bell atomization. The methods were developed with the goal of possible integration into industry to identify equipment changes through changes in the primary atomization of the bell. The first technique utilized high-speed imaging to capture qualitative ligament breakup and, in combination with a developed image processing technique and PIV software, was able to gain statistical size and velocity information about both ligaments and droplets in the image data. A second technique, using an Nd:YAG laser with an optical filter, was used to capture size statistics at even higher rotational speeds than the first technique, and was utilized to find differences between serrated and unserrated bell ligament and droplet data. The final technique was incorporating proper orthogonal decomposition (POD) into image data of a side-profile view of a damaged and undamaged bell during operation. This was done to capture differences between the data sets to come up with a characterization for identifying if a bell is damaged or not for future industrial integration.

Solid oxide cells are best known in the energy sector as novel power generation devices through solid oxide fuel cells (SOFCs), which enable the direct conversion of chemical energy to electrical energy and result in high efficiency power generation. However, solid oxide electrolysis cells (SOECs) are receiving increased attention as a hydrogen production technology through high temperature electrolysis applications. The development of higher fidelity methods for modeling transport phenomena within solid oxide cells is necessary for the advancement of these key technologies. The proposed thesis analyzes the increased transport path lengths caused by constriction resistance effects in prevalent solid oxide cell designs. Such effects are so named because they arise from reductions in active transport area. Constriction resistance effects of SOFC geometry on continuum level mass and electronic transport through SOFC anodes are simulated. These effects are explored via analytic solutions of the Laplace equation with model verification achieved by computational methods such as finite element analysis (FEA). Parametric studies of cell geometry and fuel stream composition are performed based upon the models developed. These studies reveal a competition of losses present between mass and electronic transport losses and demonstrate the benefits of smaller SOFC unit cell geometry. Furthermore, the models developed for SOFC transport phenomena are applied toward the analysis of SOECs. The resulting parametric studies demonstrate that geometric configurations that demonstrate enhanced performance within SOFC operation also demonstrate enhanced performance within SOEC operation. Secondarily, the electrochemical degradation of SOFCs is explored with respect to delamination cracking phenomena about and within the critical electrolyte-anode interface. For thin electrolytes, constriction resistance effects may lead to the loss of electro-active area at both anode-electrolyte and cathode-electrolyte interfaces. This effect (referred to as masking) results in regions of unutilized electrolyte cross-sectional area, which can be a critical performance hindrance. Again analytic and computational means are employed in analyzing such degradation issues.

Diffusive mass transport within porous SOFC components is explored using two modeling approaches that can better inform the SOFC electrode design process. These approaches include performance metrics for electrode cross-sectional design and a fractal approach for modeling mass transport within the pore structure of the electrode reaction zone. The performance metrics presented are based on existing analytical models for transport within SOFC electrodes. These metrics include a correction factor for button-cell partial pressure predictions and two forms of dimensionless reactant depletion current density. The performance impacts of multi-dimensional transport phenomena are addressed through the development of design maps that capture the trade-offs inherent in the reduction of mass transport losses within SOFC electrode cross-sections. As a complement to these bulk electrode models, a fractal model is presented for modeling diffusion within the electrochemically active region of an SOFC electrode. The porous electrode is separated into bulk and reaction zone regions, with the bulk electrode modeled in one-dimension based on the dusty-gas formalism. The reaction zone is modeled in detail with a two-dimensional finite element model using a regular Koch pore cross-section as a fractal template for the pore structure. Drawing on concepts from the analysis of porous catalysts, this model leads to a straightforward means of assessing the performance impacts of reaction zone microstructure. Together, the modeling approaches presented provide key insights into the impacts of bulk and microstructural geometry on the performance of porous SOFC components.

Past nuclear weapon production activities have left a significant legacy of uranium (U) contamination in the vadose zone (VZ) of the Department of Energy (DOE) Hanford Site. This U is a source of groundwater (GW) contamination. There is a concern that elevated U concentration would slowly infiltrate through the VZ, reach the GW water table, and then end up in nearby rivers and lakes. Remediation of U-contaminated low moisture content soil is a challenging task considering the VZ depth, where contamination is found between 70 and 100 m below the ground surface, and the formation of highly soluble and stable CaUO2CO3 complexes is influenced by Hanford’s soil rich in carbonate. Injection of reactive gasses (e.g., NH3) is a promising technology to decrease U migration in through the VZ. The NH3 injection creates alkaline conditions that would alter the pore water chemistry (e.g., dissolving some aluminosilicates). Over time as the pH neutralizes, U(VI) could precipitate as uranyl mineral (e.g., Na-boltwoodite). Also, the dissolved U(VI) could be incorporated into the structure of some mineral phases or be coated by non-U minerals. These chemical reactions could control the U(VI) mobility to the GW. However, there is a lack of knowledge on how the VZ pore water constituents (e.g., Si, Al3+, HCO3-, and Ca2+) would affect U(VI) removal/precipitation in alkaline conditions. This study quantified the role of the major pore water constituents on the U(VI) removal and evaluated the uranyl minerals that could precipitate from a variety of SPW solutions. Results showed that the percentage of U(VI) removal was controlled by Si/Al ratios and Ca2+ concentration regardless of HCO3- concentration tested. XRD revealed the presence of uranyl minerals by analyzing precipitates formed from SPW solutions, but none of them were identified as uranyl silicates as expected from speciation modeling. The SEM images displayed dense amorphous regions high in silica content, where EDS elemental analysis unveiled higher U atomic percentage in some samples. U(VI) silicate and carbonate minerals were predicted by the speciation modeling.

Apresentamos cálculos de estrutura eletrônica de oligotiofenos e oligoanilinas conectados a contatos metálicos. Diversos aglomerados de ouro e acoplamentos com as moléculas orgânicas foram testados. Expressões analíticas para a função de transmissão, de acordo com a teoria de Landauer, foram obtidas por meio de um modelo tight-binding, juntamente com resultados numéricos da teoria do funcional da densidade (DFT). O pacote computacional Gaussian 03 foi empregado para realizar os c_alculos DFT com o funcional híbrido B3LYP. Um conjunto de base misto foi escolhido: LANL2DZ para os átomo de ouro e 6-31G* para os átomos leves S, C, N e H. A geometria das cadeias orgânicas conectadas a eletrodos metálicos depende fortemente do tipo de conexão e fracamente do número de átomos representando o eletrodo. Em nossos modelos o átomo de enxofre conecta-se ao eletrodo nas geometrias bridge e atop. A principal diferença entre estes acoplamentos é uma significativa transferência de elétrons nos oligômeros ligados na forma atop. A redistribuição de carga resulta numa cadeia orgânica carregada positivamente acompanhada por uma mudança estrutural. Oligotiofenos apresentam uma distorção quinoide, enquanto que nas pernigranilinas(PE) sua alternância entre anéis quinoides e aromáticos é invertida. Um campo elétrico externo foi aplicado nos sistemas. A resposta da transmissão eletrônica _e correlacionada com o estado de oxidação das cadeias e acoplamento com os eletrodos. Nos sistemas de leucoesmeraldina (LE) e esmeraldina(EM) conectados de forma atop, o campo elétrico induz a localização dos orbitais de fronteira levando à diminuição da transmitância. Cadeias de LE e EM conectadas de forma bridge e as PE apresentam um aumento na transmissão em função do campo elétrico. Nas cadeias de PE as propriedades são menos dependentes da conexão com o eletrodo. Aplicamos uma tensão nas cadeias orgânicas forçando um alongamento, similarmente ao que acontece nos experimentos de quebra de junção. Cadeias de LE e EM conectadas de forma bridge bem como PE em ambas conexões com os eletrodos, apresentam a mesma dependência do transporte em função do estiramento. O aumento inicial da transmissão é associado à planarização da cadeia, que induz a deslocalização dos orbitais. Entretanto, o subsequente alongamento localiza os orbitais por aumentar os comprimentos de ligação, o que leva a uma baixa transmitância. Cadeias de LE e EM conectadas de forma atop são fracamente afetadas pelo alongamento. Em oligotiofenos nos dois tipos de conexão com os eletrodos, o estiramento provoca uma significativa redução da transmissão. Por este mecanismo, oligômeros mais longos apresentam valores de condutância similares aos mais curtos. / We present electronic structure and transport calculations of oligothiophenes and oligoanilines bonded to metallic contacts. Several gold clusters and couplings to the organic molecule were tested. Analytical expressions for the transmission function, within Landauer theory, were obtained by means of a tight-binding model, along with numerical results from density functional theory (DFT). The software Gaussian 03 was employed to perform DFT calculations with the hybrid functional B3LYP. A mixed basis set was chosen: LANL2DZ for gold atoms and 6-31G* for light atoms S, C, N and H. The geometry of organic chains bonded to metallic electrodes depend strongly on the type of connection with the metallic electrode and weakly on the number of atoms representing the electrode. In our models the sulfur atom connects to the electrodes in either atop or bridge geometries. The main diference between these two couplings is a significant electronic transfer on the atop bonded oligomers. The charge reorganization results in a positively charged organic part accompained by a structural change. Oligothiophenes present a quinoidal distortion, while in pernigranilines the alternation between aromatic and quinoid ring is reversed. An external electric field was applied to the systems. The electronic transmission response is correlated to the oxidation state and coupling to electrodes. On the atop bonded systems leucoemeraldine(LE) and emeraldine (EM), the electric field induces a localization of the frontier orbitals leading to a decreased transmission. Pernigranilines(PE) and bridge-bonded LE and EM present an increase of transmission as a function of the electric field. In the PE chains the properties are less dependent on the connection to the electrodes. We applied a tension on the organic chains by forcing an elongation, similarly to what happens in break-junction experiments. Bridge-bonded LE and EM, as well as PE in both types of connection to electrodes, present the same dependence of transport as function of stress. The initial increasing transmission is associated to a planarization of the chain that induces an orbital delocalization. However, the subsequent elongation localize orbitals by increasing the bond lengths and leads to a low transmission. Atop bonded LE and EM are weakly sensitive to elongation. On oligothiophenes in both types of connection to the electrodes, the stretching induces a significant decrease of transmission. By this mechanism longer oligomers present conductance values similar to small ones.

Conselho Nacional de Desenvolvimento CientÃfico e TecnolÃgico / Neste trabalho, foram estudados vÃrios tipos de escoamentos laminares de fluidos incompressÃveis Newtonianos e nÃo-Newtonianos. Isto foi feito atravÃs do estudo da interaÃÃo destes escoamentos com geometrias complexas atravÃs de modelagem computacional e da soluÃÃo numÃrica das equaÃÃes de conservaÃÃo do momento e continuidade de massa. Numa primeira etapa, a modelagem computacional de uma rede de poros foi usada para a gerar padrÃes de agregados granulares resultando de mecanismos de erosÃo-deposiÃÃo de grÃos leves. A geometria da rede de poros foi alterada dinamicamente de acordo com a transferÃncia de momento do escoamento para as partÃculas localizadas em cada vÃrtice da rede de poros. Os resultados mostraram que, para esse processo irreversÃvel, o modelo foi capaz de reproduzir padrÃes tÃpicos de processos de erosÃo bem-conhecidos. Numa segunda etapa, um separador de partÃculas semelhante à estrutura pulmonar foi proposto com base nas propriedades de escoamento em uma estrutura ramificada e nas propriedades de transporte inercial das partÃculas, quantificadas atravÃs do nÃmero de Stokes. Os resultados indicaram que a variaÃÃo dos parÃmetros de construÃÃo da estrutura ramificada leva a um regime eficiente do processo de separaÃÃo em um amplo espectro de valores do nÃmero de Stokes. Por Ãltimo, o escoamento de vÃrios fluidos nÃo-Newtonianos atravÃs de meios porosos desordenados em trÃs-dimensÃes foi estudado. Os resultados mostraram, para fluidos do tipo lei-de-potÃncia, que o escoamento pode ser descrito como uma curva universal se o nÃmero de Reynolds e a permeabilidade hidrÃulica forem redefinidos de maneira apropriada. Fluidos de Bingham tambÃm foram estudados atravÃs do modelo de Herschel-Bulkley. Neste caso, as simulaÃÃes revelaram que as interaÃÃes entre a geometria complexa do espaÃo poroso, as propriedades reolÃgicas do fluido e os efeitos inerciais do escoamento sÃo responsÃveis por uma melhora substancial da permeabilidade hidrÃulica do sistema em valores intermediÃrios do nÃmero de Reynolds. / In this work many types of incompressible laminar Newtonian and Non-Newtonian flows are studied. The interplay of these flows with complex geometries was investigated using computational modeling and numerical solution of the conservation of momentum and mass continuity equations. As a first step, the computational modeling of a network of pores was adopted to reveal the formation patterns caused by the mechanism of erosion-deposition of light grains. The geometry of the pore network was changed dynamically according to the flow momentum transfer for particles located on each vertex of the pore network. The results showed that, for this irreversible processes, the model is capable of reproducing patterns of formation of well-known erosion processes. In a second step, a particle separator inspired on the lung structure was proposed based on the flow properties in a branched structure and transport of inertial particles, quantified in terms of the Stokes number. The results indicated that the variation of construction parameters of the branched structure leads to an efficient design of the separation process in a wide range of values of the Stokes number. Finally, the flow of non-Newtonian fluids through three-dimensional disordered porous media has been studied. The results showed, for power-law fluids that the flow can de described as a universal curve if the Reynolds number and the hydraulic permeability are redefined properly. The flow of Bingham fluids was also studied using the model of Herschel-Bulkley. In this case, the simulations showed that the interaction between the complex geometry of the pore space, the rheological properties of the fluid and the inertial effects of the flow is responsible for a substantial improvement of the hydraulic permeability of the system at intermediate values of the Reynolds number.

Submitted by Renata Lopes (renatasil82@gmail.com) on 2017-03-06T18:36:40Z No. of bitstreams: 1 rafaeladeoliveiraamaral.pdf: 3948829 bytes, checksum: 00e0ed24c70d94662bfaf016bfc787cc (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2017-03-06T20:27:03Z (GMT) No. of bitstreams: 1 rafaeladeoliveiraamaral.pdf: 3948829 bytes, checksum: 00e0ed24c70d94662bfaf016bfc787cc (MD5) / Made available in DSpace on 2017-03-06T20:27:03Z (GMT). No. of bitstreams: 1 rafaeladeoliveiraamaral.pdf: 3948829 bytes, checksum: 00e0ed24c70d94662bfaf016bfc787cc (MD5) Previous issue date: 2014-02-21 / FAPEMIG - Fundação de Amparo à Pesquisa do Estado de Minas Gerais / O concreto é um material heterogêneo e poroso, sujeito a alterações físicas e químicas quando exposto a condições ambientais extremas, como as temperaturas elevadas. Como a reação de hidratação do cimento é reversível e termo-ativada, a exposição do concreto a temperaturas elevadas pode ter efeitos deletérios, com a ocorrência de desidratação da matriz a base de cimento, fissuração devido a pressões internas geradas pela evaporação da água de amassamento remanescente da mistura e ao desplacamento superficial (spalling). A exposição a temperaturas elevadas pode ocorrer, basicamente em duas situações bastante distintas: _ causas acidentais, como em incêndios, onde os gradientes de temperatura são muito altos e a exposição ao fogo se dá por um período de tempo relativamente curto; _ em condições normais de serviço, como os repositórios de rejeitos radioativos em que a estrutura fica sujeita a temperaturas muito elevadas, por longos períodos de tempo (décadas) . Em um trabalho anterior - FERREIRA (2011, [1]) -, desenvolvido no PGMC em colaboração com a Universidade de Cergy-Pontoise, foi realizada a modelagem termo-hídrico do comportamento de estruturas compostas por camadas de rocha e concreto, semelhantes ao que se encontra em túneis, elementos de instalações hidrelétricas e repositórios destinados a rejeitos nucleares. Foi, então, empregado o código computacional livre CAST3M 2000, desenvolvido no CEA (Commissariat à l'Energie Atomique, França), que é baseado nos conceitos de operadores e objetos e, para validação, foram usados resultados experimentais e numéricos. Propõe-se aqui a incorporação naquele modelo do acoplamento mecânico, através do modelo de dano de Mazars. Tal modelo é válido para situações de carregamento continuamente crescente tem como hipóteses básicas: _ O dano local resulta de deformações de alongamento evidenciadas por sinais positivos de, ao menos, um dos componentes principais de deformação; _ O dano é representado por uma variável escalar D (entre 0 e 1) cuja evolução ocorre quando o valor de alongamento equivalente excede o valor da deformação de referência. A análise dos resultados se baseia na comparação com os resultados obtidos por FERREIRA (2011, [1]) e com os dados experimentais disponíveis para os corpos-de-prova em bicamada rocha-concreto. Nota-se uma clara correspondência numérico/experimental. Discrepâncias observadas indicam a necessidade de estudos mais aprofundados acerca do efeito da degradação mecânica sobre as propriedades térmicas do concreto. / Concrete is a heterogeneous porous material subjected to physical and chemical alterations when exposed to extreme environmental conditions, such as elevated temperatures. Cement paste dehydration may lead to cracking - and even to the ocorrence of spalling - due to internal pressures generated by water vapour formation. The exposure to high temperatures may occur in two different situations: _ accidental causes, as fires, when the temperature gradients are rather elevated and fire exposure takes a relatively short period of time; _ ordinary service conditions, like repositories for radioactive wastes. A previous M.Sc. work developed in the PGMC by FERREIRA (2011, [1]), as part of a scientific cooperation with the University of Cergy Pontoise, dealt with the thermo-hydric modelling of structures composed of concrete-rock bilayers - so as to represent elements present in tunnels, hydroeletric powerplant instalations, for instance. That work employed the free code Cast3M 2000 (developed in the CEA, Commissariat à l'Energie Atomique, França), based on the concepts of operators and objects. Experimental and numerical data were employed for validation. The present work proposes the incorporation of the mechanical coupling to that TH model, through the Mazars Damage Model - with the following basic assumptions: _ Local damage results from elongation deformations; _ Damage is represented by a scalar quantity D. The numerical results are compared with data from FERREIRA (2011, [1]) and experimental information available for concrete-rock bilayer samples. It is observed a clear numerical/experimental correspondence. Resulting discrepancies indicate the need for further investigation concerning the effects of mechanical degradation on the thermical and mechanical properties of concrete.