Sensibility analysis of the pearl model for pesticide leaching in the State of Mato Grosso do Sul, Brazil

For an accurate use of pesticide leaching models it is necessary to assess the sensitivity of input parameters. The aim of this work was to carry out sensitivity analysis of the pesticide leaching model PEARL for contrasting soil types of Dourados river watershed in the state of Mato Grosso do Sul, Brazil. Sensitivity analysis was done by carrying out many simulations with different input parameters and calculating their influence on the output values. The approach used was called one-at-a-time sensitivity analysis, which consists in varying independently input parameters one at a time and keeping all others constant with the standard scenario. Sensitivity analysis was automated using SESAN tool that was linked to the PEARL model. Results have shown that only soil characteristics influenced the simulated water flux resulting in none variation of this variable for scenarios with different pesticides and same soil. All input parameters that showed the greatest sensitivity with regard to leached pesticide are related to soil and pesticide properties. Sensitivity of all input parameters was scenario dependent, confirming the need of using more than one standard scenario for sensitivity analysis of pesticide leaching models.

Transmission loss allocation based on optimal power flow and sensitivity analysis

This paper presents a new approach to the transmission loss allocation problem in a deregulated system. This approach belongs to the set of incremental methods. It treats all the constraints of the network, i.e. control, state and functional constraints. The approach is based on the perturbation of optimum theorem. From a given optimal operating point obtained by the optimal power flow the loads are perturbed and a new optimal operating point that satisfies the constraints is determined by the sensibility analysis. This solution is used to obtain the allocation coefficients of the losses for the generators and loads of the network. Numerical results show the proposed approach in comparison to other methods obtained with well-known transmission networks, IEEE 14-bus. Other test emphasizes the importance of considering the operational constraints of the network. And finally the approach is applied to an actual Brazilian equivalent network composed of 787 buses, and it is compared with the technique used nowadays by the Brazilian Control Center. (c) 2007 Elsevier Ltd. All rights reserved.

A non-linear analysis of lying a floating pipeline on the seabed

In this article, a model for the determination of displacements, deformations and tensions of a submarine pipeline during the construction is presented. The process is carried out from an initial floating situation to the final laying position on the seabed. The existence of currents and small waves are also considered. Firstly, this technique, usually applied to polyethylene pipelines, is described in this paper as well as some real world examples, as well as the variables that can be modified to control the behavior of the structure. A detailed description of the actions in this process is considered, specially the ones related to marine environment, as Archimedes force, current and sea waves. The behavior of the pipeline is modeled with a non linear elasto dynamic model where geometric non linearities are taken into account. A 3-D beam model, without cross section deformation effects, is developed. Special care is taken in the numerical analysis, developed within an updated lagrangian formulation framework, with the sea bed contact, the follower forces due to the external water pressures and the dynamic actions. Finally, some subroutines are implemented into ANSYS to simulate the two dimensional case, where the whole construction process is achieved. With this software, a sensibility analysis of the bending moments, axial forces and stresses obtained with different values of the control variables in order to optimize the construction steps. These control variables are, the axial load in the pipe, the inundated inner length and the distance of the control barge from the coast.

Cost evaluation of cardiac magnetic resonance imaging versus coronary angiography for the work-up of coronary artery disease

Background: CMR has recently emerged as a robust and reliable technique to assess coronary artery disease (CAD). A negative perfusion CMR test predicts low event rates of 0.3-0.5%/year. Invasive coronary angiography (CA) remains the "gold standard" for the evaluation of CAD in many countries.Objective: Assessing the costs of the two strategies in the European CMR registry for the work-up of known or suspected CAD from a health care payer perspective. Strategy 1) a CA to all patients or 2) a CA only to patients who are diagnosed positive for ischemia in a prior CMR.Method and results: Using data of the European CMR registry (20 hospitals, 11'040 consecutive patients) we calculated the proportion of patients who were diagnosed positive (20.6%), uncertain (6.5%), and negative (72.9%) after the CMR test in patients with known or suspected CAD (n=2'717). No other medical test was performed to patients who were negative for ischemia. Positive diagnosed patients had a coronary angiography. Those with uncertain diagnosis had additional tests (84.7%: stress echocardiography, 13.1%: CCT, 2.3% SPECT), these costs were added to the CMR strategy costs. Information from costs for tests in Germany and Switzerland were used. A sensibility analysis was performed for inpatient CA. For costs see figure. Results - costs.Discussion: The CMR strategy costs less than the CA strategy for the health insurance systems both, in Germany and Switzerland. While lower in costs, the CMR strategy is a non-invasive one, does not expose to radiation, and yields additional information on cardiac function, viability, valves, and great vessels. Developing the use of CMR instead of CA might imply some reduction in costs together with superior patient safety and comfort, and a better utilization of resources at the hospital level. Document introduit le : 01.12.2011

Dynamics of experimental populations of native and introduced blowflies (Diptera: Calliphoridae): mathematical modelling and the transition from asymptotic equilibrium to bounded oscillations

The equilibrium dynamics of native and introduced blowflies is modelled using a density-dependent model of population growth that takes into account important features of the life-history in these flies. A theoretical analysis indicates that the product of maximum fecundity and survival is the primary determinant of the dynamics. Cochliomyia macellaria, a blowfly native to the Americas and the introduced Chrysomya megacephala and Chrysomya putoria, differ in their dynamics in that the first species shows a damping oscillatory behavior leading to a one-point equilibrium, whereas in the last two species population numbers show a two-point limit cycle. Simulations showed that variation in fecundity has a marked effect on the dynamics and indicates the possibility of transitions from one-point equilibrium to bounded oscillations and aperiodic behavior. Variation in survival has much less influence on the dynamics.

Modelo de programação linear para seleção de pulverizadores agrícolas de barras

A seleção de pulverizadores agrícolas que se adaptem às necessidades da propriedade, é um processo trabalhoso, sendo uma das etapas mais importantes dentro do processo produtivo. O objetivo do presente trabalho foi o de desenvolver e utilizar um modelo de programação linear para auxiliar na seleção de pulverizadores agrícolas de barras, baseado no menor custo horário do equipamento. Foram utilizadas as informações técnicas referentes a 20 modelos de pulverizadores disponíveis no mercado, sendo quatro autopropelidos, oito de arrasto e oito do tipo montado. A análise de sensibilidade dos componentes dos custos operacionais mostrou que as taxas de reparo e depreciação foram os fatores que mais interferiram na variação do custo horário do conjunto trator-pulverizador. O modelo matemático desenvolvido facilitou a realização da análise de sensibilidade que foi processada em um tempo muito pequeno.

Ajuste preventivo versus ajuste corretivo no controle on-line de processo do número de não-conformidades num item inspecionado

In production lines, the entire process is bound to unexpected happenings which may cost losing the production quality. Thus, it means losses to the manufacturer. Identify such causes and remove them is the task of the processing management. The on-line control system consists of periodic inspection of every month produced item. Once any of those items is quali ed as not t, it is admitted that a change in the fraction of the items occurred, and then the process is stopped for adjustments. This work is an extension of Quinino & Ho (2010) and has as objective main to make the monitoramento in a process through the control on-line of quality for the number of non-conformities about the inspected item. The strategy of decision to verify if the process is under control, is directly associated to the limits of the graphic control of non-conformities of the process. A policy of preventive adjustments is incorporated in order to enlarge the conforming fraction of the process. With the help of the R software, a sensibility analysis of the proposed model is done showing in which situations it is most interesting to execute the preventive adjustment

Dynamics of experimental populations of native and introduced blowflies (Diptera: Calliphoridae): mathematical modelling and the transition from asymptotic equilibrium to bounded oscillations

The equilibrium dynamics of native and introduced blowflies is modelled using a density-dependent model of population growth that takes into account important features of the life-history in these flies. A theoretical analysis indicates that the product of maximum fecundity and survival is the primary determinant of the dynamics. Cochliomyia macellaria, a blowfly native to the Americas and the introduced Chrysomya megacephala and Chrysomya putoria, differ in their dynamics in that the first species shows a damping oscillatory behavior leading to a one-point equilibrium, whereas in the last two species population numbers show a two-point limit cycle. Simulations showed that variation in fecundity has a marked effect on the dynamics and indicates the possibility of transitions from one-point equilibrium to bounded oscillations and aperiodic behavior. Variation in survival has much less influence on the dynamics.

A comparative study of the damping oscillation function of TCSC and UPFC

The main objective of this work is to analyze the ability of FACTS devices like TCSC and UPFC to damp low frequency oscillations and a POD controller is also included. A comparative study of damping effect of those devices IS carried out. The Power Sensitivity Model (PSM) is used to the representation of the electric power system. Sensibility analysis using the residue method shows the best place for the installation of FACTS and the procedure to determine POD parameters. ©2008 IEEE.

Análise de sensibilidade através de um modelo implicitamente acoplado para alívio de sobrecargas em redes de transmissão

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

A manutenção do equilíbrio econômico de contratos de obras públicas empreitadas: parâmetros para o estudo de riscos vinculados a atrasos no encaixe do preço

Dentre os diversos segmentos do Setor da Construção Civil, um dos principais é o de obras por empreitada global, que são aquelas contratadas por preço certo e total. As obras públicas, que têm como finalidade atender a uma utilidade pública e cujos contratos são regidos pela Lei 8.666/93, são geralmente contratadas sob esse regime de execução. O artigo 58 desta Lei determina que deva ser mantido o equilíbrio econômicofinanceiro dos contratos celebrados com o Poder Público. A fundamentação teórica trata dos conceitos de obras públicas, dos principais aspectos da Lei 8.666/93, de análise econômico-financeira e de análise de riscos. Este trabalho apresenta também uma sistemática e uma simulação numérica para formação de preços para obras empreitadas de Construção Civil. O objetivo geral deste trabalho é a criação de parâmetros para o estabelecimento de margens de proteção para a dilatação do prazo e respectivos encaixes de preços, visando, como determina a Lei 8.666/93, a manutenção do equilíbrio original. Com base na análise dos resultados de 63 obras, verifica-se a manutenção do equilíbrio econômico-financeiro, comparando-se a taxa de retorno restrita esperada pelas empresas no momento das licitações com a realmente alcançada na operação. A partir dos desvios detectados nas taxas de retorno nesses dois momentos, cria-se, pelo processo de simulação e análise de sensibilidade, parâmetros para cobertura dos referidos riscos. No momento da licitação, a taxa de retorno restrita média das obras analisadas é de 9,45% ao mês. No segundo momento, essa taxa cai para em média para 5,16% ao mês. Isto representa, em média, uma variação de 45,39%. A taxa média de atratividade das obras públicas analisadas está dentro do intervalo de 6% a 13% ao mês. Como parâmetro de cobertura de riscos quanto aos desvios de prazo e encaixe do preço, tem-se o intervalo de 1,5% a 11,5% do preço. Portanto, entre as obras analisadas, não existe o equilíbrio econômico-financeiro.

Determinação das velocidades intervalares usando a teoria paraxial do raio: aproximação de segunda ordem dos tempos de trânsito

Neste trabalho foi desenvolvido um método de solução ao problema inverso para modelos sísmicos compostos por camadas homogêneas e isotrópicas separadas por superfícies suaves, que determina as velocidades intervalares em profundidade e calcula a geometria das interfaces. O tempo de trânsito é expresso como uma função de parâmetros referidos a um sistema de coordenadas fixo no raio central, que é determinada numericamente na superfície superior do modelo. Essa função é posteriormente calculada na interface anterior que limita a camada não conhecida, através de um processo que determina a função característica em profundidade. A partir da função avaliada na interface anterior se calculam sua velocidade intervalar e a geometria da superfície posterior onde tem lugar a reflexão do raio. O procedimento se repete de uma forma recursiva nas camadas mais profundas obtendo assim a solução completa do modelo, não precisando em nenhum passo informação diferente à das camadas superiores. O método foi expresso num algoritmo e se desenvolveram programas de computador, os quais foram testados com dados sintéticos de modelos que representam feições estruturais comuns nas seções geológicas, fornecendo as velocidades em profundidade e permitindo a reconstrução das interfaces. Uma análise de sensibilidade sobre os programas mostrou que a determinação da função característica e a estimação das velocidades intervalares e geometria das interfaces são feitos por métodos considerados estáveis. O intervalo empírico de aplicabilidade das correções dinâmicas hiperbólicas foi tomado como uma estimativa da ordem de magnitude do intervalo válido para a aplicação do método.

Empilhamento pelo método superfície de reflexão comum 2-D com topografia e introdução ao caso 3-D

O método de empilhamento sísmico CRS simula seções sísmicas ZO a partir de dados de cobertura múltipla, independente do macro-modelo de velocidades. Para meios 2-D, a função tempo de trânsito de empilhamento depende de três parâmetros, a saber: do ângulo de emergência do raio de reflexão normal (em relação à normal da superfície) e das curvaturas das frentes de onda relacionadas às ondas hipotéticas, denominadas NIP e Normal. O empilhamento CRS consiste na soma das amplitudes dos traços sísmicos em dados de múltipla cobertura, ao longo da superfície definida pela função tempo de trânsito do empilhamento CRS, que melhor se ajusta aos dados. O resultado do empilhamento CRS é assinalado a pontos de uma malha pré-definida na seção ZO. Como resultado tem-se a simulação de uma seção sísmica ZO. Isto significa que para cada ponto da seção ZO deve-se estimar o trio de parâmetros ótimos que produz a máxima coerência entre os eventos de reflexão sísmica. Nesta Tese apresenta-se fórmulas para o método CRS 2-D e para a velocidade NMO, que consideram a topografia da superfície de medição. O algoritmo é baseado na estratégia de otimização dos parâmetros de fórmula CRS através de um processo em três etapas: 1) Busca dos parâmetros, o ângulo de emergência e a curvatura da onda NIP, aplicando uma otimização global, 2) busca de um parâmetro, a curvatura da onda N, aplicando uma otimização global, e 3) busca de três parâmetros aplicando uma otimização local para refinar os parâmetros estimados nas etapas anteriores. Na primeira e segunda etapas é usado o algoritmo Simulated Annealing (SA) e na terceira etapa é usado o algoritmo Variable Metric (VM). Para o caso de uma superfície de medição com variações topográficas suaves, foi considerada a curvatura desta superfície no algoritmo do método de empilhamento CRS 2-D, com aplicação a dados sintéticos. O resultado foi uma seção ZO simulada, de alta qualidade ao ser comparada com a seção ZO obtida por modelamento direto, com uma alta razão sinal-ruído, além da estimativa do trio de parâmetros da função tempo de trânsito. Foi realizada uma nálise de sensibilidade para a nova função de tempo de trânsito CRS em relação à curvatura da superfície de medição. Os resultados demonstraram que a função tempo de trânsito CRS é mais sensível nos pontos-médios afastados do ponto central e para grandes afastamentos. As expressões da velocidade NMO apresentadas foram aplicadas para estimar as velocidades e as profundidades dos refletores para um modelo 2-D com topografia suave. Para a inversão destas velocidades e profundidades dos refletores, foi considerado o algoritmo de inversão tipo Dix. A velocidade NMO para uma superfície de medição curva, permite estimar muito melhor estas velocidades e profundidades dos refletores, que as velocidades NMO referidas as superfícies planas. Também apresenta-se uma abordagem do empilhamento CRS no caso 3-D. neste caso a função tempo de trânsito depende de oito parâmetros. São abordadas cinco estratégias de busca destes parâmetros. A combinação de duas destas estratégias (estratégias das três aproximações dos tempos de trânsito e a estratégia das configurações e curvaturas arbitrárias) foi aplicada exitosamente no empilhamento CRS 3-D de dados sintéticos e reais.

Aproximações hiperbólicas dos tempos de trânsito de raios paraxiais a um raio central refletido e difratado

A simulação de uma seção sísmica de afastamento nulo (AN) a partir de dados de cobertura múltipla para um meio 2-D, através do empilhamento, é um método de imageamento de reflexão sísmica muito utilizado, que permite reduzir a quantidade de dados e melhor a razão sinal/ruído. Baseado na aproximação hiperbólica dos tempos de trânsito dependente de três parâmetros ou atributos cinemáticos de frentes de onda, recentemente, vem desenvolvendo-se um novo método para simular seções (AN) chamado método de empilhamento sísmico por Superfície de Reflexão Comum (ou empilhamento SRC). Também, seguindo este novo conceito de imageamento sísmico está surgindo um método para simular seções com afastamento comum (AC) a partir de dados de cobertura múltipla usando aproximações dos tempos de trânsito paraxiais na vizinhança de um raio central com afastamento finito. Esta nova aproximação dos tempos de trânsito depende de cinco atributos cinemáticos. Neste trabalho, a partir da aproximação dos tempos de trânsito paraxiais em relação a um raio central com afastamento finito, derivamos uma nova equação do tempo de trânsito usando a condição de um ponto difrator em profundidade, reduzindo a equação original para quatro parâmetros. Para ambas aproximações (reflexão e difração), mostramos a superfície de empilhamento SRC com afastamento finito (SRC-AF). Considerando um modelo sintético, realizamos um estudo comparativo das aproximações dos tempos de trânsito para as quatro configurações sísmicas (fonte comum (FC), receptor comum (RC), ponto-médio-comum (PMC) e afastamento comum (AC)). Para analisar o comportamento do operador SRC-AF, quando este é perturbado, discutimos sua sensibilidade em relação a cada um dos cinco parâmetros (K₁, K₂, K₃, β_S e β_G). Esta análise de sensibilidade é realizada em duas formas: Sensibilidade através da primeira derivada e Sensibilidade no Empilhamento SRC-AF. Após realizar a análise de sensibilidade utilizamos uma nova condição, K₂ = 0 e assim, obtemos uma nova aproximação, agora dependente de três parâmetros. Usando essas aproximações hiperbólicas (em função de cinco, quatro e três parâmetros), propomos um algoritmo para a simulação de seções AC a partir de dados de cobertura múltipla. Finalmente, é apresentado um estudo da zona de Fresnel, com o objetivo de determinar a delimitação da abertura da superfície de empilhamento SRC-AF.

Aproximações hiperbólicas dos tempos de trânsito com topografia

A simulação de uma seção sísmica de afastamento-nulo (AN) a partir de dados de cobertura múltipla pode ser realizada através do empilhamento sísmico, o qual é um método de imageamento de reflexão sísmica muito utilizado na indústria. O processo de empilhamento sísmico permite reduzir a quantidade de dados e melhorar a razão sinal/ruído. Baseado na aproximação hiperbólica dos tempos de trânsito dependente de três parâmetros ou atributos cinemáticos de frentes de onda. Recentemente, vem desenvolvendo-se um novo método para simular seções (AN) chamado método de empilhamento sísmico por Superfície de Reflexão Comum (ou empilhamento SRC). Este novo formalismo pode ser estendido para construir seções de afastamento-nulo (AN) a partir de dados de cobertura múltipla, usando aproximações dos tempos de trânsito paraxiais na vizinhança de um raio central com afastamento-nulo (AN), para o caso de uma linha de medição com topografia suave e rugosa. Essas duas aproximações de tempos de trânsito também dependem de três atributos cinemáticos de frentes de ondas. Nesta dissertação, apresenta-se uma revisão teórica da teoria paraxial do raio para a obtenção das aproximações dos tempos de trânsito paraxiais considerando uma linha de medição com topografia rugosa e suave. A partir das aproximações dos tempos de trânsito paraxiais em relação a um raio central com afastamento-nulo (AN), foram obtidas duas novas aproximações de tempos de trânsito usando a condição de um ponto difrator em profundidade, reduzindo as equações originais para dois parâmetros. Também foram obtidas as aproximações para o caso de raios paraxiais com afastamento-nulo (AN). Para as aproximações de reflexão e difração utilizando um mesmo modelo sintético, foram comparadas através da representação gráfica as superfícies de empilhamento das aproximações dos tempos de trânsito para a topografia suave e rugosa. Em seguida, para analisar o comportamento dos operadores associados a reflexão e a difração, quando estes são perturbados, discutimos suas sensibilidades em relação a cada um dos três parâmetros (β₀, K_PIN, K_N). Esta análise de sensibilidade foi realizada em duas formas: Sensibilidade através da perturbação de cada parâmetro visualizado nas superfícies de empilhamento SRC-TR e SDC-TR e da primeira derivada dos tempos de trânsito SRC-TR e SDC-TR. Finalmente, usando essas aproximações hiperbólicas em função de três e dois parâmetros e com base nos resultados da análise de sensibilidade, foi proposto um algoritmo para simular seções AN a partir de dados de cobertura múltipla.