Novos modelos para cálculo de energia livre de solvatação em simulações de dinâmica molecular

Combinando métodos computacionais da eletrostática e dinâmica molecular, este trabalho teve como objetivo descrever processos de solvatação em sistemas quimicamente importantes. Foram obtidas propriedades termodinâmicas necessárias para o entendimento do processo de solvatação. Estão descritos e testados modelos para descrever a interação soluto-solvente, possibilitando, assim, aprimorar a descrição físico-química dos processos de solvatação. Utilizaram-se programas desenvolvidos em nosso grupo e programas comerciais que permitem os cálculos de dinâmica molecular e química quântica. Uma nova abordagem para o cálculo de energia livre de solvatação foi desenvolvida proporcionando a obtenção acurada e eficiente dessa propriedade, dentro do enfoque da dinâmica molecular. Nessa nova abordagem, novas metodologias para a geração de cavidades moleculares foram propostas e avaliadas. As energias livres de solvatação obtidas estão em boa concordância com os valores experimentais.

Análise estática e dinâmica, linear e não-linear geométrica, através de elementos hexaédricos de oito nós com um ponto de integração

Para a análise estática e dinâmica, linear e não-linear de placas, cascas e vigas, implementa-se neste trabalho o elemento hexaédrico com integração reduzida, livre de travamento volumétrico e travamento de cisalhamento e que não apresenta modos espúrios. Na formulação do elemento, utiliza-se apenas um ponto de integração. Desta forma, a matriz de rigidez é dada de forma explícita e o tempo computacional é significativamente reduzido, especialmente em análise não-linear. Os modos espúrios são suprimidos através de um procedimento de estabilização que não exige parâmetros especificados pelo usuário. Para evitar o travamento de cisalhamento, desenvolve-se o vetor de deformações num sistema co-rotacional e remove-se certos termos não constantes nas componentes de deformações de cisalhamento. O travamento volumétrico é resolvido fazendo-se com que a parte dilatacional (esférica) da matriz gradiente seja avaliada apenas no ponto central do elemento. Como a eliminação do travamento de cisalhamento depende de uma abordagem no sistema local, emprega-se um procedimento co-rotacional para obter o incremento de deformação no sistema local e atualizar os vetores de tensões e forças internas na análise não-linear Para a solução das equações de equilíbrio na análise estática, utilizam-se métodos diretos baseados na eliminação de Gauss ou métodos iterativos de Gradientes Conjugados Precondicionado elemento-por-elemento (EBE). Para a análise dinâmica, as equações de equilíbrio são integradas através do método explícito de Taylor-Galerkin ou do método implícito de Newmark. Através de exemplos numéricos demonstra-se a eficiência e o potencial do elemento tridimensional na análise de casca, placas e vigas submetidas a grandes deslocamentos e grande rotações. Os resultados são comparados com trabalhos que utilizam elementos clássicos de placa e casca.

Parâmetros bioecológicos de Gryon gallardoi (Brethes) (Hymenoptera: Scelionidae) e modelagem da dinâmica espaço-temporal da sua interação com Spartocera dentiventris (Berg) (Hemiptera: Coreidae) através da simulação de múltiplos agentes

Estudos referentes à dinâmica espaço-temporal da interação hospedeiro-parasitóide são de fundamental importância para o desenvolvimento de técnicas de manejo de insetos em agroecossistemas. Uma recente abordagem têm sido o estudo de populações através da simulação computacional dos indivíduos que a compõem. Este método tem possibilitado a incorporação de importantes características como individualidade e estrutura espacial nos modelos teóricos. Neste sentido, o presente estudo visou a elaboração de um modelo baseado no comportamento individual para simular a dinâmica espaço-temporal da interação Spartocera dentiventris (Berg) (Hemiptera: Coreidae), inseto associado à cultura do fumo, e de seu inimigo natural, o parasitóide de ovos Gryon gallardoi (Brethes) (Hymenoptera: Scelionidae). Para a elaboração do modelo, estudos sobre a bioecologia de G. gallardoi foram desenvolvidos, de modo a investigar (i) o efeito da temperatura no seu desenvolvimento e viabilidade, (ii) a mortalidade de imaturos em campo, (iii) os parâmetros reprodutivos e a longevidade dos adultos, (iv) os padrões de dispersão em campo e (v) a resposta funcional e a interferência mútua dos parasitóides em diferentes densidades de ovos de S. dentiventris. No intuito de avaliar a capacidade de predição do modelo elaborado, a dinâmica espaço-temporal de ovos sadios e parasitados de S. dentiventris foi acompanhada em um cultivo de fumo, estabelecido em área experimental do Departamento de Fitossanidade da Faculdade de Agronomia da UFRGS, Porto Alegre (30o01’S e 51o13’O), RS, Brasil A viabilidade do desenvolvimento dos parasitóides na faixa de 20 a 30oC não diferiu significativamente, alcançando 98,8%. O tempo de desenvolvimento ovo-adulto de machos e fêmeas foi inversamente proporcional ao aumento da temperatura. Os valores estimados para o limite térmico inferior de desenvolvimento e para a constante térmica foram 15,5oC e 185,19GD para machos e 15,6oC e 192,31GD para fêmeas. Diversos fatores afetam o sucesso de imaturos de G. gallardoi em campo, sendo que o malogro e a predação por sugadores são os principais responsáveis pela emergência de apenas 37,87% dos adultos. Foi observado um período médio de oviposição de 10,1 ± 1,74 dias, com o pico de oviposição no segundo dia, sendo depositados ao longo dos mesmos uma média de 67,5 ± 11,29 ovos por fêmea. Fêmeas de G. gallardoi foram significativamente mais longevas que os machos, vivendo, respectivamente, 13,7 ± 1,94 e 10,6 ± 1,78 dias. A razão sexual total observada foi de 0,79. A dispersão de G. gallardoi em cultivo de fumo não foi sensivelmente influenciada pelo vento, sendo estimada uma capacidade de locomoção diária média para fêmeas de no mínimo 7,62 m. O padrão de parasitismo de G. gallardoi ajustou-se perfeitamente à resposta funcional do tipo II (pseudo-r2= 0,99), sendo obtidos os valores de 0,0557 e 0,9989 h para os componentes a’ e Tm , respectivamente O aumento da densidade de parasitóides acarretou, de maneira geral, uma diminuição no número de ovos parasitados por parasitóide. Foi obtido para o índice de interferência mútua “m” o valor 0,626. A dinâmica temporal de S. dentiventris – G. gallardoi simulada ajustou-se muito bem aos dados obtidos em campo (r2= 0,82 para ovos sadios e r2= 0,72 para ovos parasitados). Da mesma forma, os arranjos espaciais dos ovos e da taxa de parasitismo observados em campo foram satisfatoriamente reproduzidos através das simulações computacionais. Os resultados indicam que G. gallardoi apresenta um bom potencial para controlar populações de S. dentiventris, seja pela sua viabilidade de desenvolvimento em uma larga faixa de temperatura ou pela sua capacidade reprodutiva e padrões de dispersão similares àqueles observados para outras espécies já utilizadas no controle biológico de insetos. Por outro lado, a sensibilidade às baixas temperaturas, a alta mortalidade de imaturos observada em campo e uma resposta funcional do tipo II podem ser fatores que se contrapõe ao sucesso do parasitóide no controle de populações de S. dentiventris. Considerando a independência entre os dados observados e os dados obtidos através de simulações, o modelo elaborado obteve pleno sucesso em simular a dinâmica da interação S. dentiventris – G. gallardoi a partir do comportamento individual dos diferentes agentes Duas importantes características associadas à estabilidade de sistemas hospedeiro-parasitóide foram observadas no estudo de campo e reproduzidas no modelo – um padrão de parasitismo inversamente dependente da densidade dos hospedeiros e uma resposta agregativa dos parasitóides em algumas regiões do sistema. Desta forma, futuros estudos podem avaliar, a partir de um enfoque baseado no indivíduo, quais mecanismos são responsáveis por tais comportamentos e qual é o seu efeito na dinâmica do sistema. Conclui-se, através deste trabalho, que o estudo da dinâmica de populações através de uma abordagem baseada no indivíduo, com base em modelos computacionais, pode ser uma alternativa viável para superar as limitações e a complexidade inerente ao uso de modelos populacionais analíticos. Desta forma, o maior realismo implícito nestes modelos os torna ferramentas muito úteis para diversas áreas da ecologia aplicada, particularmente para a entomologia agrícola.

Ferramenta computacional para a medição de campos de velocidade utilizando processamento digital de imagens

A proposta deste trabalho, consiste na elaboração de uma ferramenta computacional para a medição de campos de velocidades em escoamentos com baixas velocidades (< 0,5 m/s) utilizando o processamento digital de imagens. Ao longo dos anos, inúmeras técnicas foram desenvolvidas com este objetivo. Para cada tipo de aplicação, uma técnica se aplica com maior ou menor eficiência do que outras. Para o caso de estudos em fluídos transparentes, onde o escoamento pode ser visualizado, técnicas que utilizam processamento digital de imagens vêm ganhando um grande impulso tecnológico nos últimos anos. Este impulso, é devido a fatores como: câmaras vídeo filmadoras de última geração, dispositivos de aquisição de imagens e componentes de processamento e armazenamento de dados cada vez mais poderosos. Neste contexto, está a velocimetria por processamento de imagens de partículas cuja sigla é PIV (particle image velocimetry). Existem várias formas de se implementar um sistema do tipo PIV. As variantes dependem, basicamente, do equipamento utilizado. Para sua implementação é necessário, inicialmente, um sistema de iluminação que incide em partículas traçadoras adicionadas ao fluido em estudo. Após, as partículas em movimento são filmadas ou fotografadas e suas imagens adquiridas por um computador através de dispositivos de captura de imagens. As imagens das partículas são então processadas, para a obtenção dos vetores velocidade. Existem diferentes formas de processamento para a obtenção das velocidades. Para o trabalho em questão, devido às características dos equipamentos disponíveis, optou-se por uma metodologia de determinação da trajetória de partículas individuais, que, apesar de limitada em termos de módulo de velocidade, pode ser aplicada a muitos escoamentos reais sob condições controladas Para validar a ferramenta computacional desenvolvida, imagens ideais de partículas foram simuladas como se estivessem em escoamento, através do deslocamento conhecido de vários pixels. Seguindo o objetivo de validação, foi utilizada ainda uma imagem real de partículas, obtida com o auxílio de um plano de iluminação de luz coerente (LASER) e câmaras de vídeo tipo CCD. O programa desenvolvido foi aplicado em situações de escoamento real e os resultados obtidos foram satisfatórios dentro da escala de velocidades inicialmente presumida.

Estudo computacional de líquidos iônicos do tipo dialquilimidazólio

Neste trabalho foram desenvolvidos parâmetros de campo de força dos ânions tetrafenilborato (BPh4 -) e hexafluorfosfato (PF6 -) dentro da metodologia AMBER para a simulação computacional por dinâmica molecular de Líquidos Iônicos formados por estes ânions e cátions do tipo dialquilimidazólio, o 1-n-butil-3-metilimidazólio (BMI+), o 1-etil-3-metilimidazólio (EMI+) e o 1,3-dimetilimidazólio (MMI+). A validação destes parâmetros foi realizada por comparação entre as freqüências dos modos normais obtido por cálculo ab initio com aquelas obtidas por mecânica molecular, juntamente com uma comparação entre as estruturas moleculares e momentos multipolares obtidos pelas duas metodologias. Seguiu-se então a validação por comparação dos resultados dos cálculos de dinâmica molecular com dados experimentais, como densidades, entalpias de vaporização, condutividade elétrica, estrutura radial e espacial e também dados de difração de nêutrons. Foi atingida uma concordância bastante grande entre dados experimentais e cálculo teórico principalmente no que diz respeito à estrutura dos Líquidos Iônicos e foi possível racionalizar em termos de tamanho do grupo alquila do cátion (n-butil, metil ou etil) e do tamanho do ânion tanto propriedades estruturais deste líquidos quanto características dinâmicas dos mesmos.

Modelagem estratigráfica de clinoformas deposicionais : construção e aplicação de um modelo computacional baseado em mecânica estatística

Sob a premissa de que a geometria do perfil deposicional das clinoformas pode conter informações sobre a composição textural dos sedimentos que as compõem, bem como sobre a energia da coluna d’água, foi desenvolvido um modelo computacional para a análise teórica dos fatores controladores do tipo de curvatura dos foresets das clinoformas, e, por extensão dos taludes submarinos. Um modelo análogo de suspensão dinâmica foi implementado experimentalmente com um programa em código Matlab, denominado MAXWELL, e é classificado como um modelo estratigráfico, bidimensional, analítico, determinístico, dinâmico, parcialmente dedutivo e parcialmente baseado em regras. Contém um modelo de escoamento de fluido em linhas de corrente, e trata indiretamente a tensão de cisalhamento no domínio de um sistema fechado bifásico líquido-vapor análogo, a partir a proposta de analogias entre os processos de evaporação e de suspensão de sedimentos. É uma abordagem baseada na competência de transporte do fluxo d’água, pois considera temperatura e velocidade de fluxo combinado onda-corrente como variáveis controladoras análogas da intensidade de transferência de energia. Os processos deposicionais marinhos são reduzidos a um fenômeno da superfície deposicional, que é considerada análoga a uma interface líquidovapor dentro de um sistema fechado. A equação de distribuição estatística de velocidades moleculares de Maxwell é usada para calcular a proporção de moléculas na fase líquida, considerada análoga à proporção de sedimentos depositados, bem como a proporção na fase vapor, tomada como análoga à proporção de sedimentos em suspensão. O estado energético do sistema é parametrizado por três potenciais: energia interna, como função do tamanho de grão (areia, silte ou argila); energia térmica do meio, como função da energia hidrodinâmica, e energia gravitacional, como função do declive topográfico As simulações indicam que os principais fatores controladores do perfil deposicional das clinoformas, e por extensão, dos taludes submarinos em geral, são a energia hidrodinâmica da coluna d’água e a granulometria (ou coesão) do sedimento, que podem ser consideradas dois parâmetros comutáveis, isto é, grãos grossos ou coesos podem produzir sobre a geometria das clinoformas o mesmo efeito da baixa energia hidrodinâmica, e vice-versa. Com base no fator de decaimento da energia hidrodinâmica com o aumento da batimetria, foram definidas três ordens de grandeza da intensidade da energia da coluna d’água: baixa energia (10–1), alta x energia (10-2) e muito alta energia (10-3). Com base nesse critério, foram caracterizados quatro tipos de perfis deposicionais de clinoformas: dois tipos sigmoidais e dois tipos exponenciais. Os tipos sigmoidais podem ser de alta energia ou de muito alta energia, e distinguem-se pela granulometria do sedimento e pela distribuição de declividades em relação à dimensão na direção horizontal. Os perfis de tipo exponencial podem ser tanto de baixa energia, quanto de alta energia. A subida relativa do nível do mar afeta o tipo geométrico do perfil deposicional apenas de forma indireta, a partir da retenção da fração grossa do influxo sedimentar na plataforma ou topset. Os principais fatores controladores do ângulo de declividade do perfil das clinoformas são a granulometria do sedimento, a energia hidrodinâmica, a profundidade d’água da bacia e o desvio lateral da corrente de fundo. A dedução da litofácies a partir da geometria das clinoformas parece promissora apenas para os tipos teóricos sigmoidais, que são distintos na forma e no conteúdo sedimentar.

Dinâmica monetária eficiente sob encontros aleatórios: uma classe de métodos numéricos que exploram concavidade

A dificuldade em se caracterizar alocações ou equilíbrios não estacionários é uma das principais explicações para a utilização de conceitos e hipóteses que trivializam a dinâmica da economia. Tal dificuldade é especialmente crítica em Teoria Monetária, em que a dimensionalidade do problema é alta mesmo para modelos muito simples. Neste contexto, o presente trabalho relata a estratégia computacional de implementação do método recursivo proposto por Monteiro e Cavalcanti (2006), o qual permite calcular a sequência ótima (possivelmente não estacionária) de distribuições de moeda em uma extensão do modelo proposto por Kiyotaki e Wright (1989). Três aspectos deste cálculo são enfatizados: (i) a implementação computacional do problema do planejador envolve a escolha de variáveis contínuas e discretas que maximizem uma função não linear e satisfaçam restrições não lineares; (ii) a função objetivo deste problema não é côncava e as restrições não são convexas; e (iii) o conjunto de escolhas admissíveis não é conhecido a priori. O objetivo é documentar as dificuldades envolvidas, as soluções propostas e os métodos e recursos disponíveis para a implementação numérica da caracterização da dinâmica monetária eficiente sob a hipótese de encontros aleatórios.

Ferramentas para análise dinâmica e estrutural de um reboque de linha leve

Neste trabalho é desenvolvida uma metodologia de projeto para identificar as regiões críticas da estrutura de um reboque de linha leve sendo tracionado em pavimentos do tipo rodovia de baixa qualidade e estrada secundária de muito baixa qualidade. Para tanto, são levantados alguns dados experimentais da estrutura, necessários para a aproximação e simulação dinâmica de um modelo simplificado. A excitação da base é realizada por atuadores que simulam as oscilações verticais de um perfil de estrada, a qual é definida de acordo com os estudos realizados por Dodds e Robson (1973). Isto permite a determinação de um histórico de carregamentos das regiões da estrutura do chassi sob a ação das molas da suspensão. Em seguida, é gerado um modelo estrutural simplificado do reboque em elementos finitos, chamado de global, no qual são determinadas as regiões sob ação das maiores tensões. Tendo identificada a região mais crítica da estrutura, é criado um modelo local desta parte, onde se pode observar a distribuição de tensões com mais detalhe, permitindo a identificação dos pontos de concentração de tensões. Desta forma, com a aplicação do método de análise global-local é possível a obtenção de resultados detalhados quanto aos esforços da estrutura com um menor custo computacional.

Determinação da resposta dinâmica de feixes de condutores em linhas aéreas de transmissão

Até os dias de hoje, a abordagem estática tem sido a usualmente empregada para a avaliação da resposta de cabos de Linhas de Transmissão (LTs), apesar da resposta dinâmica ser, em muitos casos, reconhecidamente importante na avaliação do desempenho dos cabos. Para uma análise completa, não se pode desconhecer a natureza dinâmica da maioria dos fenômenos a que as LTs estão submetidas, sendo exemplos típicos as excitações mecânicas causadas pela ação do vento (carregamentos de baixa e alta freqüência) e a ruptura de cabos. Cabe salientar que mesmo a análise estática empregada tem sido simplificada, considerando normalmente apenas casos contemplados com soluções analíticas. A justificativa para emprego da análise estática baseia-se no menor esforço numérico exigido. Entretanto, essa justificativa já não se sustenta dado os grandes avanços na área computacional, possibilitando o estudo do desempenho de LTs na ocorrência de fenômenos que provocam carregamentos dinâmicos. Adicionalmente, há uma crescente demanda para que se obtenha um melhor entendimento de muitas questões relativas ao comportamento dinâmico das LTs, que possam explicar desempenhos observados e/ou sustentar o desenvolvimento de novas alternativas mais arrojadas de projeto e construção Entre os diversos métodos utilizados na engenharia estrutural, a integração direta das equações do movimento, através de métodos numéricos como as diferenças finitas centrais, se constitui numa poderosa ferramenta de cálculo. Esta ferramenta possibilita o tratamento de problemas envolvendo não linearidades geométricas e do material, como são os casos onde se avalia a resposta de cabos suspensos submetidos a carregamentos variáveis no tempo. Esta pesquisa objetiva a aplicação do método da integração direta das equações do movimento na análise de feixes de cabos de LT, quando submetidos à ação de carregamentos mecânicos variáveis no tempo, principalmente à excitação de ventos oriundos de fenômenos com natureza complexa (tormentas elétricas, por exemplo). É apresentado um método para determinar a resposta dinâmica de feixes que considera a interação entre o vento incidente e o movimento do condutor. A solução é obtida por integração numérica, no domínio do tempo, das equações de movimento de um modelo tridimensional não-linear discreto de feixe, as quais definem as forças nos espaçadores e cabos através de coeficientes aerodinâmicos obtidos experimentalmente. Também são apresentados modelos dos fenômenos meteorológicos mais comuns em nosso país (Brasil): tormenta extratropical (EPS) e tormenta elétrica (TS) Como ilustração, são apresentados exemplos de modelagem de vãos de LTs com condutor singelo e com feixes de condutores. Os exemplos demonstram a capacidade de avaliação de Estados Limites relacionados à distância relativa entre subcondutores, à estabilidade do feixe, à representação das suas propriedades e ao comportamento dinâmico, bem como aos carregamentos transmitidos às estruturas. A análise emprega conhecidas relações constitutivas para representar o comportamento tensão-deformação dos cabos. O enfoque utilizado possibilita a avaliação mais precisa de casos reais que ainda não podiam ser convenientemente tratados, além de permitir a extensão para estudos bem mais complexos, tais como feixes com disposições assimétricas.

Movimento de partículas carregadas em fluidos ionizados : fundamentos matemáticos da teoria de eletroforese capilar

Neste trabalho, discutimos o movimento de uma macromolécula carregada em um fluido ionizado. A interação do campo elétrico é descrita pela equação de Poisson-Boltzmann acoplada às equações governantes para a dinâmica do fluido e às equações dinâmicas da partícula. Uma formulação fraca é introduzida no caso em que o domínio ocupado pelo fluido é finito e um teorema de existência de soluções fracas, local em tempo, é estabelecido. Dois modelos são considerados: fluxos não-estacionários e estacionários. No primeiro caso, a hidrodinâmica do sistema é governada pelas equações de Navier-Stokes, considerando-se um termo forçante relacionado ao potencial elétrico; no segundo caso, uma velocidade de deslizamento, a qual depende não linearmente sobre os potenciais, é introduzida como uma condição de contorno para um problema estacionário de Stokes. O caso de um fluido ocupando uma região infinita é também discutido supondo-se uma hipótese de aproximação sobre o campo elétrico.

Influência de parâmetros moleculares em funções de correlação temporal na dinâmica de solvatação mecânica

No presente trabalho descrevemos nossos resultados relativos à investigação da dinâmica de solvatação mecânica por meio de simulações por dinâmica molecular, respeitando o regime da resposta linear, em sistemas-modelo de argônio líquido com um soluto monoatômico ou diatômico dissolvido. Estudamos sistematicamente a influência dos parâmetros moleculares dos solutos (tamanho, polarizabilidade) e da densidade frente a vários modelos de solvatação. Funções de Correlação Temporal da Energia de Solvatação foram calculadas com relação à correlações de n-corpos (n = 2; 3) distinguindo interações repulsivas e atrativas para ambos os sistemas líquidos. Também obtivemos segundas derivadas temporais dessas funções referindo-se à parcelas translacionais, rotacionais e roto-translacionais na solução do diatômico. Encontramos que funções de correlação temporal coletivas podem ser razoavelmente bem aproximadas por correlações binárias a densidades baixas e, a densidades altas, correlações ternárias tornam-se mais importantes produzindo um descorrelacionamento mais rápido das funções coletivas devido a efeitos de cancelamento parciais. As funções de correlação para interações repulsivas e atrativas exibem comportamentos dinâmicos independentes do modelo de solvatação devido a fatores de escalonamento linear que afetam apenas as amplitudes das dessas funções de correlação temporal. Em geral, os sistemas com grau de liberdade rotacional apresentam tempos de correlação mais curtos para a dinâmica coletiva e tempos de correlação mais longos para as funções binárias e ternárias. Finalmente, esse estudo mostra que os sistemas contendo o diatômico relaxam-se predominantemente por mecanismos translacionais binários em modelos de solvatação envolvendo alterações apenas na polarizabilidade do soluto, e por mecanismos rotacionais atrativos binários em modelos envolvendo alterações no comprimento de ligação.

Análise estática e dinâmica de placas e cascas de materiais compósitos laminados usando elementos finitos hexaédricos de oito nós com integração reduzida

Neste trabalho implementou-se o elemento hexaédrico com um ponto de integração para análise estática e dinâmica de placas e cascas de materiais compósitos laminados com ou sem enrijecedores. O elemento está livre de travamento volumétrico e travamento de cisalhamento, não apresentando modos espúrios. São também incluídos problemas com não-linearidade geométrica A matriz de rigidez e de massa são dadas de forma explícita, reduzindo o tempo computacional, especialmente na análise não-linear. Para evitar o travamento de cisalhamento as componentes de deformações são referidas a um sistema co-rotacional. O travamento volumétrico é também eliminado, já que a parte dilatacional (esférica) da matriz gradiente é avaliada apenas no ponto central do elemento. Para a solução das equações de equilíbrio na análise estática, utilizam-se um método direto baseado na eliminação de Gauss ou um método iterativo de gradientes conjugados com precondicionamento executado através da eliminação incompleta de Choleski. Para a análise dinâmica, as equações de equilíbrio são integradas através do método explícito ou implícito de Taylor-Galerkin ou do método implícito de Newmark. Para análise não-linear utiliza-se o Método Generalizado de Controle dos Deslocamentos. Através de exemplos numéricos demonstra-se a eficiência e o potencial do elemento tridimensional na análise linear e não-linear de placas e cascas de materiais laminados. Os resultados são comparados com trabalhos que utilizam diferentes elementos de placas e cascas.

Ventilação e prescrições urbanísticas: uma aplicação simulada para a Orla da Praia do Meio em Natal/RN

This work analyses a study on natural ventilation and its relation to the urban legislation versus the building types in an urban fraction of coastal area of Praia do Meio in the city of Natal/RN, approaching the type or types of land use most appropriate to this limited urban fraction. The objective of this study is to analyse the effects of the present legislation as well as the types of buildings in this area on the natural ventilation. This urban fraction was selected because it is one of the sites from where the wind flows into the city of Natal. This research is based on the hypothesis stating that the reduction on the porosity of the urban soil (decrease in the set back/boundary clearance), and an increase in the form (height of the buildings) rise the level of the ventilation gradient, consequently causing a reduction on the wind speed at the lowest part of the buildings. Three-dimensional computational models were used to produce the modes of occupation allowed in the urban fraction within the area under study. A Computational Fluid Dynamics (CFD) software was also used to analyse the modes of land occupation. Following simulation, a statistical assessment was carried out for validation of the hypothesis. It was concluded that the reduction in the soil porosity as a consequence of the rates that defined the minimum boundary clearance between the building and the boundary of the plot (and consequently the set back), as well as the increase in the building form (height of the buildings) caused a reduction in the wind speed, thus creating heat islands

Análise de componentes arquitetônicos para potencialização da ventilação natural com ênfase em captadores de vento

The building envelope is the principal mean of interaction between indoors and environment, with direct influence on thermal and energy performance of the building. By intervening in the envelope, with the proposal of specific architectural elements, it is possible to promote the use of passive strategies of conditioning, such as natural ventilation. The cross ventilation is recommended by the NBR 15220-3 as the bioclimatic main strategy for the hot and humid climate of Natal/RN, offering among other benefits, the thermal comfort of occupants. The analysis tools of natural ventilation, on the other hand, cover a variety of techniques, from the simplified calculation methods to computer fluid dynamics, whose limitations are discussed in several papers, but without detailing the problems encountered. In this sense, the present study aims to evaluate the potential of wind catchers, envelope elements used to increase natural ventilation in the building, through CFD simplified simulation. Moreover, it seeks to quantify the limitations encountered during the analysis. For this, the procedure adopted to evaluate the elements implementation and efficiency was the CFD simulation, abbreviation for Computer Fluid Dynamics, with the software DesignBuilder CFD. It was defined a base case, where wind catchers were added with various settings, to compare them with each other and appreciate the differences in flows and air speeds encountered. Initially there has been done sensitivity tests for familiarization with the software and observe simulation patterns, mapping the settings used and simulation time for each case simulated. The results show the limitations encountered during the simulation process, as well as an overview of the efficiency and potential of wind catchers, with the increase of ventilation with the use of catchers, differences in air flow patterns and significant increase in air speeds indoors, besides changes found due to different element geometries. It is considered that the software used can help designers during preliminary analysis in the early stages of design