Calibração de redes de distribuição de água pelo método inverso aplicado a escoamento transiente

Desenvolve-se um método para estimar os parâmetros de uma rede hidráulica a partir de dados observados de cargas hidráulicas transientes. Os parâmetros físicos da rede como fatores de atrito, rugosidades absolutas, diâmetros e a identificação e quantificação de vazamentos são as grandezas desconhecidas. O problema transiente inverso é resolvido utilizando uma abordagem indireta que compara os dados disponíveis de carga hidráulica transiente observados com os calculados através de um método matemático. O Método Transiente Inverso (MTI) com um Algoritmo Genético (AG) emprega o Método das Características (MOC) na solução das equações do movimento para escoamento transiente em redes de tubos. As condições de regime permanente são desconhecidas. Para avaliar a confiabilidade do MTI-AG desenvolvido aqui, uma rede-exemplo é usada para os vários problemas de calibração propostos. O comportamento transiente é imposto por duas manobras distintas de uma válvula de controle localizada em um dos nós da rede. Analisam-se, ainda, o desempenho do método proposto mediante a variabilidade do tamanho do registro transiente e de possíveis erros de leitura nas cargas hidráulicas. Ensaios numéricos realizados mostram que o método é viável e aplicável à solução de problema inverso em redes hidráulicas, sobretudo recorrendo-se a poucos dados observados e ao desconhecimento das condições iniciais de estado permanente. Nos diversos problemas de identificação, as informações transientes obtidas da manobra mais brusca produziu estimações mais eficientes.

Determinação da tenacidade à fratura em corte (GIIc) de adesivos estruturais pelo ensaio End-Notched Flexure (ENF)

Os adesivos têm sido alvo de estudo ao longo dos últimos anos para ligação de componentes a nível industrial. Devido à crescente utilização das juntas adesivas, torna-se necessária a existência de modelos de previsão de resistência que sejam fiáveis e robustos. Neste âmbito, a determinação das propriedades dos adesivos é fundamental para o projeto de ligações coladas. Uma abordagem recente consiste no uso de modelos de dano coesivo (MDC), que permitem simular o comportamento à fratura das juntas de forma bastante fiável. Esta técnica requer a definição das leis coesivas em tração e corte. Estas leis coesivas dependem essencialmente de 2 parâmetros: a tensão limite e a tenacidade no modo de solicitação respetivo. O ensaio End-Notched Flexure (ENF) é o mais utilizado para determinar a tenacidade em corte, porque é conhecido por ser o mais expedito e fiável para caraterizar este parâmetro. Neste ensaio, os provetes são sujeitos a flexão em 3 pontos, sendo apoiados nas extremidades e solicitados no ponto médio para promover a flexão entre substratos, o que se reflete numa solicitação de corte no adesivo. A partir deste ensaio, e após de definida a tenacidade em corte (GIIc), existem alguns métodos para estimativa da lei coesiva respetiva. Nesta dissertação são definidas as leis coesivas em corte de três adesivos estruturais através do ensaio ENF e um método inverso de ajuste dos dados experimentais. Para o efeito, foram realizados ensaios experimentais considerado um adesivo frágil, o Araldite® AV138, um adesivo moderadamente dúctil, o Araldite® 2015 e outro dúctil, o SikaForce® 7752. O trabalho experimental consistiu na realização dos ensaios ENF e respetivo tratamento dos dados para obtenção das curvas de resistência (curvas-R) através dos seguintes métodos: Compliance Calibration Method (CCM), Direct Beam Theory (DBT), Corrected Beam Theory (CBT) e Compliance-Based Beam Method (CBBM). Os ensaios foram simulados numericamente pelo código comercial ABAQUS®, recorrendo ao Métodos de Elementos Finitos (MEF) e um MDC triangular, com o intuito de estimar a lei coesiva de cada um dos adesivos em solicitação de corte. Após este estudo, foi feita uma análise de sensibilidade ao valor de GIIc e resistência coesiva ao corte (tS 0), para uma melhor compreensão do efeito destes parâmetros na curva P- do ensaio ENF. Com o objetivo de testar adequação dos 4 métodos de obtenção de GIIc usados neste trabalho, estes foram aplicados a curvas P- numéricas de cada um dos 3 adesivos, e os valores de GIIc previstos por estes métodos comparados com os respetivos valores introduzidos nos modelos numéricos. Como resultado do trabalho realizado, conseguiu-se obter uma lei coesiva única em corte para cada um dos 3 adesivos testados, que é capaz de reproduzir com precisão os resultados experimentais.

Avaliação de parâmetros hidráulicos para modelos de distribuição de água no solo sob gotejamento

As soluções analíticas de distribuição de água para ponto fonte e regime de fluxo não- permanente são dependentes de parâmetros de solos considerados constantes em suas deduções. Erros na determinação desses parâmetros implicam insucesso dessas soluções. Este trabalho teve por objetivo avaliar o comportamento dos parâmetros alfada equação de Gardner k(h) = Ks ealfah e o parâmetro de linearização k = dK(teta)/dteta em diferentes posições do bulbo molhado, considerando os processos de infiltração isolado e seguido de redistribuição de água, à semelhança de um ciclo de irrigação. Dados de potencial matricial foram coletados em diversas posições do bulbo molhado em duas situações:(i) início da irrigação até atingir regime permanente em todo o bulbo molhado (infiltração); e (ii) durante dois ciclos de irrigação envolvendo infiltração e redistribuição de água. Os resultados mostraram que os parâmetros alfa e k variaram nas posições do bulbo molhado em relação ao gotejador, de acordo com o regime de umidade a que tais posições estiveram sujeitas. A obtenção desses parâmetros pelo método inverso requer testes que considerem as fases de infiltração e redistribuição em pelo menos dois ciclos de irrigação.

Diseño de un sistema regional de alerta de tsunamis como parte integral de la operación portuaria

En este trabajo se presenta el desarrollo de una metodología para obtener un universo de funciones de Green y el algoritmo correspondiente, para estimar la altura de tsunamis a lo largo de la costa occidental de México en función del momento sísmico y de la extensión del área de ruptura de sismos interplaca localizados entre la costa y la Trinchera Mesoamericana. Tomando como caso de estudio el sismo ocurrido el 9 de octubre de 1995 en la costa de Jalisco-Colima, se estudiaron los efectos del tsunami originados en la hidrodinámica del Puerto de Manzanillo, México, con una propuesta metodológica que contempló lo siguiente: El primer paso de la metodología contempló la aplicación del método inverso de tsunamis para acotar los parámetros de la fuente sísmica mediante la confección de un universo de funciones de Green para la costa occidental de México. Tanto el momento sísmico como la localización y extensión del área de ruptura de sismos se prescribe en segmentos de planos de falla de 30 X 30 km. A cada uno de estos segmentos del plano de falla corresponde un conjunto de funciones de Green ubicadas en la isobata de 100 m, para 172 localidades a lo largo de la costa, separadas en promedio 12 km entre una y otra. El segundo paso de la metodología contempló el estudio de la hidrodinámica (velocidades de las corrientes y niveles del mar en el interior del puerto y el estudio del runup en la playa) originada por el tsunami, la cual se estudió en un modelo hidráulico de fondo fijo y en un modelo numérico, representando un tsunami sintético en la profundidad de 34 m como condición inicial, el cual se propagó a la costa con una señal de onda solitaria. Como resultado de la hidrodinámica del puerto de Manzanillo, se realizó un análisis de riesgo para la definición de las condiciones operativas del puerto en términos de las velocidades en el interior del mismo, y partiendo de las condiciones iniciales del terremoto de 1995, se definieron las condiciones límites de operación de los barcos en el interior y exterior del puerto. In this work is presented the development of a methodology in order to obtain a universe of Green's functions and the corresponding algorithm in order to estimate the tsunami wave height along the west coast of Mexico, in terms of seismic moment and the extent of the area of the rupture, in the interplate earthquakes located between the coast and the Middle America Trench. Taking as a case of study the earthquake occurred on October 9, 1995 on the coast of Jalisco-Colima, were studied the hydrodynamics effects of the tsunami caused in the Port of Manzanillo, Mexico, with a methodology that contemplated the following The first step of the methodology contemplated the implementation of the tsunami inverse method to narrow the parameters of the seismic source through the creation of a universe of Green's functions for the west coast of Mexico. Both the seismic moment as the location and extent of earthquake rupture area prescribed in segments fault planes of 30 X 30 km. Each of these segments of the fault plane corresponds a set of Green's functions located in the 100 m isobath, to 172 locations along the coast, separated on average 12 km from each other. The second step of the methodology contemplated the study of the hydrodynamics (speed and directions of currents and sea levels within the port and the study of the runup on the beach Las Brisas) caused by the tsunami, which was studied in a hydraulic model of fix bed and in a numerical model, representing a synthetic tsunami in the depth of 34 m as an initial condition which spread to the coast with a solitary wave signal. As a result of the hydrodynamics of the port of Manzanillo, a risk analysis to define the operating conditions of the port in terms of the velocities in the inner and outside of the port was made, taken in account the initial conditions of the earthquake and tsunami ocurred in Manzanillo port in 1995, were defined the limits conditions of operation of the ships inside and outside the port.

Influência da gestão do habitat na distribuição e abundância do coelho-bravo no Sítio Natura 2000 Monchique

O coelho-bravo, devido à sua importância ecológica e económica, tem sido alvo de diversos planos de gestão e vários esforços têm sido empreendidos no sentido de contrariar o decréscimo das suas populações. Este estudo foi realizado em três zonas de caça do Sítio Monchique e o principal objectivo é determinar se as medidas de gestão implementadas influenciam a distribuição e abundância da espécie na área de estudo. A abundância relativa foi interpolada com o método "Inverso do Peso da Distância" {IDW), e as relações entre presença de coelho e os descritores ambientais foram analisadas através de Modelos Lineares Generalizados (GLM). Os resultados da modelação estatística mostraram que as medidas de melhoria de habitat parecem ter sido determinantes para um aumento da área de distribuição do coelho-bravo nos locais intervencionados. São propostas novas medidas de gestão, cujo objectivo será promover a continuação do aumento da ocorrência e abundância da espécie neste local. /ABSTRACT: The wild rabbit, due to its ecological and economical role, has been the target of several management plans and considerable efforts have been made to enhance its populations. This study was held in three game estates located inside Monchique Natura 2000. Site and aims to determine if the habitat management actions implemented in the study area influence rabbit distribution and abundance. The relative abundance was interpolated to all study area with lnverse Distance Weight method {IDW} and the relationships between rabbit presence and the environmental descriptors were evaluated with Generalized Linear Models (GLM). The results of the statistical modelling showed that the management actions seem to have contributed significantly to an enhancement on the rabbit occurrence in the studied game estates. Several new management actions are proposed with the aim to continue to increase rabbit occurrence and abundance in this site.

Método de inversão matricial rápido GSGR para vectoring em sistemas xDSL, com implementação em chip DSP

Este trabalho apresenta um método rápido de inversão de matrizes densas, e uma possível aplicação com métodos de Vectoring, em pré-codificação e cancelamento de crosstalk de sistemas xDSL. A família de tecnologias xDSL utiliza os pares trançados de fios de cobre telefônicos como meio físico para transmitir dados digitais. O crosstalk é a principal causa de degradação de sinais na mais nova geração de sistemas xDSL, o G.fast, e para combatê-lo são utilizadas técnicas de pré-codificação e cancelamento, chamadas de Vectoring. O método proposto, chamado de GSGR, consiste em uma abordagem diferente para o método clássico de Squared Givens Rotations (SGR), adequado a implementações em plataformas embarcadas de processamento digital de sinais. Foram realizados testes comparativos do método GSGR com métodos diretos clássicos de inversão, utilizando uma plataforma digital multicore baseada no chip TI DSP TMS320C6670 e a plataforma de software Matlab. Os resultados dos testes de inversão de matrizes usando dados reais e dados simulados mostraram que o GSGR foi superior em velocidade de execução sem apresentar perdas significativas de acurácia para a aplicação em sistemas xDSL.

Estudo do método de inversão tomográfica 2D baseada nos atributos cinemáticos da onda NIP e construção de interface gráfica para o processo de picking

Neste trabalho, estuda-se um novo método de inversão tomográfica de reflexão para a determinação de um modelo isotrópico e suave de velocidade por meio da aplicação, em dados sintéticos e reais, do programa Niptomo que é uma implementação do método de inversão tomográfica dos atributos cinemáticos da onda hipotética do ponto de incidência normal (PIN). Os dados de entrada para a inversão tomográfica, isto é, o tempo de trânsito e os atributos da onda PIN (raio de curvatura da frente de onda emergente e ângulo de emergência), são retirados de uma série de pontos escolhidos na seção afastamento nulo (AN) simulada, obtida pelo método de empilhamento por superfícies de reflexão comum (SRC). Normalmente, a escolha destes pontos na seção AN é realizada utilizando-se programas de picking automático, que identificam eventos localmente coerentes na seção sísmica com base nos parâmetros fornecidos pelo usuário. O picking é um dos processos mais críticos dos métodos de inversão tomográfica, pois a inclusão de dados de eventos que não sejam de reflexões primárias podem ser incluídos neste processo, prejudicando assim o modelo de velocidades a ser obtido pela inversão tomográfica. Este trabalho tem por objetivo de construir um programa de picking interativo para fornecer ao usuário o controle da escolha dos pontos de reflexões sísmicas primárias, cujos dados serão utilizados na inversão tomográfica. Os processos de picking e inversão tomográfica são aplicados nos dados sintéticos Marmousi e nos dados da linha sísmica 50-RL-90 da Bacia do Tacutu. Os resultados obtidos mostraram que o picking interativo para a escolha de pontos sobre eventos de reflexões primárias favorece na obtenção de um modelo de velocidade mais preciso.

A primitivação pelo método de substutuição

Expõem-se neste texto os fundamentos do método de primitivação por substituição, estudado nas disciplinas de Análise Matemática de muitos dos cursos de licenciatura de natureza científica ou técnica. Trata-se de um método simples, embora gere mais dificuldades junto dos alunos do que o método de primitivação por partes. E, à semelhança do que se fez em trabalho similar sobre este último método de primitivação, ilustra-se o agora apresentado com a doutrina essencial e com uma coletânea de exemplos úteis.

Krigagem e inverso do quadrado da distância para interpolação dos parâmetros da equação de chuvas intensas

Em projetos de obras hidráulicas, é necessária a estimativa da chuva intensa adequada à realidade técnico-financeira dos projetos, para posterior cálculo de vazões máximas. O uso da equação de chuvas intensas é a forma mais usual para estimativa da chuva de projeto. No entanto, há dificuldade em obtê-la para locais desprovidos de dados pluviométricos, o que pode ser superado mediante técnicas computacionais para a espacialização dos seus parâmetros. Para isto, pode-se trabalhar com procedimentos estatísticos para interpolação dos parâmetros com base em estações meteorológicas vizinhas ao local. Objetivou-se, neste trabalho, comparar os dois procedimentos estatísticos mais exatos, conforme já constatado em outras aplicações físicas, que são a krigagem (interpolador geoestatístico) e o inverso do quadrado da distância, para interpolação dos parâmetros da equação de chuvas intensas. Utilizaram-se 140 estações meteorológicas do estado de São Paulo, as quais dispõem da respectiva equação de chuvas intensas estimada com base em dados pluviométricos, consideradas como ideais. Destas, 126 foram utilizadas para estudo da dependência espacial dos parâmetros e 14 para comparação de erros nas estimativas da chuva intensa. Constatou-se que ambos os métodos apresentaram boa precisão, mas a krigagem produziu menores erros médios para 11 estações, verificando-se melhorias consideráveis em especial para tempos de retorno e tempos de duração da precipitação usuais em projetos em pequenas bacias hidrográficas, podendo-se sugerir este método como o mais adequado.

Método para a espacialização dos elementos do balanço hídrico climatológico

O objetivo deste trabalho foi desenvolver e avaliar um método de espacialização dos elementos do balanço hídrico climatológico (BHC). O método utiliza modelos digitais do terreno que contêm a distribuição espacial dos elementos climáticos de entrada do BHC. A avaliação foi feita com dados reais de BHC referentes a 110 estações pluviométricas do Estado do Espírito Santo. Comparou-se o desempenho do método proposto com o dos seguintes interpoladores espaciais: krigagem e inverso de uma potência da distância. O método proposto apresentou melhores estimativas espaciais de todos os elementos do balanço hídrico - evapotranspiração real, excedente, deficiência e disponibilidade hídrica.

Solução analítica da equação de difusão-advenção pelo método GILTT aplicada à dispersão de poluentes atmosféricos

O objetivo deste trabalho é obter uma nova solução analítica para a equação de advecção-difusão. Para tanto, considera-se um problema bidimensional difusivo-advectivo estacionário com coeficiente de difusão turbulenta vertical variável que modela a dispersão de poluentes na atmosfera. São utilizados três coeficientes difusivos válidos na camada limite convectiva e que dependem da altura, da distância da fonte e do perfil de velocidade. A abordagem utilizada para a resolução do problema é a técnica da Transformada Integral Generalizada, na qual a equação transformada do problema difusivo-advectivo é resolvida pela técnica da Transformada de Laplace com inversão analítica. Nenhuma aproximação é feita durante a derivação da solução, sendo assim, esta é exata exceto pelo erro de truncamento. O modelo ´e avaliado em condições moderadamente instáveis usando o experimento de Copenhagen. Apresentam-se os resultados numéricos e estatísticos, comparando os resultados obtidos com dados experimentais e com os resultados da literatura. O modelo proposto mostrou-se satisfatório em relação aos dados dos experimentos difusivos considerados.

Solução analítica das equações difusivas da teoria geral de perturbação pelo método da transformada de Laplace

Neste trabalho, apresentamos uma solução analítica para as equações difusivas unidimensionais da Teoria Geral de Perturbação em uma placa heterogênea, isto é, apresentamos as soluções analíticas para os problemas de autovalor para o fluxo de nêutrons e para o fluxo adjunto de nêutrons, para o cálculo do fator de multiplicação efetivo (keff), para o problema de fonte fixa e para o problema de função auxiliar. Resolvemos todos os problemas mencionados aplicando a Transformada de Laplace em uma placa heterogênea considerando um modelo de dois grupos de energia e realizamos a inversão de Laplace do fluxo transformado analiticamente através da técnica da expansão de Heaviside. Conhecendo o fluxo de nêutrons, exceto pelas constantes de integração, aplicamos as condições de contorno e de interface e resolvemos as equações algébricas homogêneas para o fator de multiplicação efetivo pelo método da bissecção. Obtemos o fluxo de nêutrons através da avaliação das constantes de integração para uma potência prescrita. Exemplificamos a metodologia proposta para uma placa com duas regiões e comparamos os resultados obtidos com os existentes na literatura.

Solução da equação de difusão-advecção para uma CLP não-homogênea e não-estacionária pelo método GILTT

Neste trabalho é apresentada a solução da equação de difusão-advecção transiente para simular a dispersão de poluentes na Camada Limite Planetária. A solução é obtida através do método analítico GILTT (Generalized Integral Laplace Transform Technique) e da técnica de inversão numérica da quadratura de Gauss. A validação da solução é comprovada utilizando as concentraçãos obtidas a partir do modelo com as obtidas experimentalmente pelo Experimento de Copenhagen. Nesta comparação foram utilizados os perfis de vento potencial e logaritmo e os parâmetros de turbulência propostos por Degrazia et al (1997) [19] e (2002) [17]. Os melhores resultados foram obtidos utilizando o perfil de vento potencial e o coeficiente de difusão propostos por Degrazia et al (1997). A influência da velocidade vertical é mostrada através do comportamento das concentrações de poluentes na pluma. Além disso, as velocidades verticais e longitudinais geradas pelo Large Eddy Simulation (LES) foram colocadas no modelo para poder simular uma camada limite turbulenta mais realística, a qual apresentou resultados satisfatórios quando comparados com os disponíveis na literatura.

Um método para localização facial Invariante a rotação

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior

Aperfeiçoamento no projeto de um protótipo de misturador-decantador à inversão de fases visando aplicação na indústria do petróleo

In the oil industry the mixture oil/water occurs in the operations of production, transportation and refining, as well as during the use of its derivatives. The amount of water produced associated with the oil varies and can reach values of 90% in volume in the case of mature phase of the production fields. The present work deals with the development of new design of the Mixer Settler based on Phase Inversion (MDIF) in a laboratory scale. We envisage this application in industrial scale so the phases of project, construction and operation are considered. The modifications most significant, in comparison with the original prototype, include the materials of construction and the substitution of the equipment used in the mixing stage of the process. It was tested the viability of substitution of the original system of mechanical mixing by a static mixer. A statistical treatment by means of an experimental design of composed central type was used in order to evaluate the behavior of the main variables of the separation process as function of the efficiency of separation for the new device. This procedure is useful to delimit an optimal region of operation with the equipment. The variables of process considered on the experimental design were: oil concentration in the feeding water (mg/L); Total volumetric flow rate (L/h); Ratio organic/water on volumetric basis (O/A). The separation efficiency is calculated by comparison of the content of oil and greases in the inlet and outlet of the equipment. For determination of TOG (Total Oil and Grease), the method used was based in the absorption of radiation in the infra-red region. The equipment used for these determinations was InfraCal® TOG/TPH Model HATR-T2 of the Wilks Enterprise, Incorporation. It´s important to stand out that this method of measure has being used by PETROBRAS S.A. Results of global efficiency of separation oil/water varied from 75.3 to 97.7% for contaminated waters containing up to 1664,1 mg/L of oil. By means of tests carried out with a real sample of contaminated water supplied by PETROBRAS we have got an effluent specified in terms of the legal standards required for discharging. Thus, the new design of equipment constitutes a real alternative for the conventional systems of treatment of produced water in the oil industry