Measures of observation impact in non-Gaussian data assimilation

ABSTRACT Non-Gaussian/non-linear data assimilation is becoming an increasingly important area of research in the Geosciences as the resolution and non-linearity of models are increased and more and more non-linear observation operators are being used. In this study, we look at the effect of relaxing the assumption of a Gaussian prior on the impact of observations within the data assimilation system. Three different measures of observation impact are studied: the sensitivity of the posterior mean to the observations, mutual information and relative entropy. The sensitivity of the posterior mean is derived analytically when the prior is modelled by a simplified Gaussian mixture and the observation errors are Gaussian. It is found that the sensitivity is a strong function of the value of the observation and proportional to the posterior variance. Similarly, relative entropy is found to be a strong function of the value of the observation. However, the errors in estimating these two measures using a Gaussian approximation to the prior can differ significantly. This hampers conclusions about the effect of the non-Gaussian prior on observation impact. Mutual information does not depend on the value of the observation and is seen to be close to its Gaussian approximation. These findings are illustrated with the particle filter applied to the Lorenz ’63 system. This article is concluded with a discussion of the appropriateness of these measures of observation impact for different situations.

Tests of non-linearity using LIFFE futures transactions price data

This paper presents and implements a number of tests for non-linear dependence and a test for chaos using transactions prices on three LIFFE futures contracts: the Short Sterling interest rate contract, the Long Gilt government bond contract, and the FTSE 100 stock index futures contract. While previous studies of high frequency futures market data use only those transactions which involve a price change, we use all of the transaction prices on these contracts whether they involve a price change or not. Our results indicate irrefutable evidence of non-linearity in two of the three contracts, although we find no evidence of a chaotic process in any of the series. We are also able to provide some indications of the effect of the duration of the trading day on the degree of non-linearity of the underlying contract. The trading day for the Long Gilt contract was extended in August 1994, and prior to this date there is no evidence of any structure in the return series. However, after the extension of the trading day we do find evidence of a non-linear return structure.

Testing for non-linearity in daily sterling exchange rates

A number of tests for non-linear dependence in time series are presented and implemented on a set of 10 daily sterling exchange rates covering the entire post Bretton-Woods era until the present day. Irrefutable evidence of non-linearity is shown in many of the series, but most of this dependence can apparently be explained by reference to the GARCH family of models. It is suggested that the literature in this area has reached an impasse, with the presence of ARCH effects clearly demonstrated in a large number of papers, but with the tests for non-linearity which are currently available being unable to classify any additional non-linear structure.

An evaluation of three upwinding approximations for numerical modeling the flow in tubular membrane of Newtonian and non-Newtonian fluids

In this paper, the laminar fluid flow of Newtonian and non-Newtonian of aqueous solutions in a tubular membrane is numerically studied. The mathematical formulation, with associated initial and boundary conditions for cylindrical coordinates, comprises the mass conservation, momentum conservation and mass transfer equations. These equations are discretized by using the finite-difference technique on a staggered grid system. Comparisons of the three upwinding schemes for discretization of the non-linear (convective) terms are presented. The effects of several physical parameters on the concentration profile are investigated. The numerical results compare favorably with experimental data and the analytical solutions. (C) 2011 Elsevier Inc. All rights reserved.

Utilização de modelos de programação não-linear para sistematização de terras para irrigação em áreas regulares

Objetivou-se com este trabalho, desenvolver modelos de programação não-linear para sistematização de terras, aplicáveis para áreas com formato regular e que minimizem a movimentação de terra, utilizando o software GAMS para o cálculo. Esses modelos foram comparados com o Método dos Quadrados Mínimos Generalizado, desenvolvido por Scaloppi & Willardson (1986), sendo o parâmetro de avaliação o volume de terra movimentado. Concluiu-se que, ambos os modelos de programação não-linear desenvolvidos nesta pesquisa mostraram-se adequados para aplicação em áreas regulares e forneceram menores valores de movimentação de terra quando comparados com o método dos quadrados mínimos.

Modelagem linear e identificação do modelo dinâmico de um robô móvel com acionamento diferencial

This work presents a modelling and identification method for a wheeled mobile robot, including the actuator dynamics. Instead of the classic modelling approach, where the robot position coordinates (x,y) are utilized as state variables (resulting in a non linear model), the proposed discrete model is based on the travelled distance increment Delta_l. Thus, the resulting model is linear and time invariant and it can be identified through classical methods such as Recursive Least Mean Squares. This approach has a problem: Delta_l can not be directly measured. In this paper, this problem is solved using an estimate of Delta_l based on a second order polynomial approximation. Experimental data were colected and the proposed method was used to identify the model of a real robot

Controlador preditivo não linear aplicado ao controle de golfadas em processos de produção de petróleo

Slugging is a well-known slugging phenomenon in multiphase flow, which may cause problems such as vibration in pipeline and high liquid level in the separator. It can be classified according to the place of its occurrence. The most severe, known as slugging in the riser, occurs in the vertical pipe which feeds the platform. Also known as severe slugging, it is capable of causing severe pressure fluctuations in the flow of the process, excessive vibration, flooding in separator tanks, limited production, nonscheduled stop of production, among other negative aspects that motivated the production of this work . A feasible solution to deal with this problem would be to design an effective method for the removal or reduction of the system, a controller. According to the literature, a conventional PID controller did not produce good results due to the high degree of nonlinearity of the process, fueling the development of advanced control techniques. Among these, the model predictive controller (MPC), where the control action results from the solution of an optimization problem, it is robust, can incorporate physical and /or security constraints. The objective of this work is to apply a non-conventional non-linear model predictive control technique to severe slugging, where the amount of liquid mass in the riser is controlled by the production valve and, indirectly, the oscillation of flow and pressure is suppressed, while looking for environmental and economic benefits. The proposed strategy is based on the use of the model linear approximations and repeatedly solving of a quadratic optimization problem, providing solutions that improve at each iteration. In the event where the convergence of this algorithm is satisfied, the predicted values of the process variables are the same as to those obtained by the original nonlinear model, ensuring that the constraints are satisfied for them along the prediction horizon. A mathematical model recently published in the literature, capable of representing characteristics of severe slugging in a real oil well, is used both for simulation and for the project of the proposed controller, whose performance is compared to a linear MPC

Conventional and microwave sintering of CaCu(3)Ti(4)O(12)/CaTiO(3) ceramic composites: non-ohmic and dielectric properties

The non-ohmic and dielectric properties as well as the dependence on the microstructural features of CaCu(3)Ti(4)O(12)/CaTiO(3) ceramic composites obtained by conventional and microwave sintering were investigated. It was demonstrated that the non-ohmic and dielectric properties depend strongly on the sintering conditions. It was found that the non-linear coefficient reaches values of 65 for microwave-sintered samples and 42 for samples sintered in a conventional furnace when a current density interval of 1-10 mA cm(-2) is considered. The non-linear coefficient value of 65 is equivalent to 1500 for samples sintered in the microwave if a current interval of 5-30 mA is considered as is shortly discussed by Chung et al (2004 Nature Mater. 3 774). Due to a high non-linear coefficient and a low leakage current (90 mu A) under both processing conditions, these samples are promising for varistor applications. The conventionally sintered samples exhibit a higher relative dielectric constant at 1 kHz (2960) compared with the samples sintered in the microwave furnace (2100). At high frequencies, the dielectric constant is also larger in the samples sintered in the conventional furnace. Depending on the application, one or another synthesis methodology is recommended, that is, for varistor applications sintered in a microwave furnace and for dielectric application sintered in a conventional furnace.

Comparative degradation of ZnO- and SnO(2)-based polycrystalline non-ohmic devices by current pulse stress

The degradation behaviour of SnO(2)-based varistors (SCNCr) due to current pulses (8/20 mu s) is reported here for the first time in comparison with the ZnO-based commercial varistors (ZnO). Puncturing and/or cracking failures were observed in ZnO-based varistors possessing inferior thermo-mechanical properties in comparison with that found in a SCNCr system free of failures. Both systems presented electric degradation related to the increase in the leakage current and decrease in the electric breakdown field, non-linear coefficient and average value of the potential barrier height. However, it was found that a more severe degradation occurred in the ZnO-based varistors concerning their non-ohmic behaviour, while in the SCNCr system, a strong non-ohmic behaviour remained after the degradation. These results indicate that the degradation in the metal oxide varistors is controlled by a defect diffusion process whose rate depends on the mobility, the concentration of meta-stable defects and the amount of electrically active interfaces. The improved behaviour of the SCNCr system is then inferred to be associated with the higher amount of electrically active interfaces (85%) and to a higher energy necessary to activate the diffusion of the specific defects.

Primary power distribution systems planning taking into account reliability, operation and expansion costs

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Signatures of primordial non-Gaussianities in the matter power-spectrum and bispectrum: the time-RG approach

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Short comments on self-synchronization of two non-ideal sources supported by a flexible portal frame structure

A practical problem of synchronization of a non-ideal (i.e. when the excitation is influenced by the response of the system) and non-linear vibrating system was posed and investigated by means of numerical simulations. Two rotating unbalanced motors compose the mathematical model considered here with limited power supply mounted on the horizontal beam of a simple portal frame. As a starting point, the problem is reduced to a four-degrees-of-freedom model and its equations of motion, derived elsewhere via a Lagrangian approach, are presented. The numerical results show the expected phenomena associated with the passage through resonance with limited power. Further, for a two-to-one relationship between the frequencies associated with the first symmetric mode and the sway mode, by using the variation of torque constants, the control of the self-synchronization and synchronization (in the system) are observed at certain levels of excitations.

A reducing of a chaotic movement to a periodic orbit, of a micro-electro-mechanical system, by using an optimal linear control design

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Study of the thermal behavior of the transition phase of Co(II)-diclofenac compound by non-isothermal method

The Co(II)-diclofenac complex was evaluated by simultaneous thermogravimetry-differential thermal analysis (TG-DTA) and differential scanning calorimetry (DSC). The DTA curve profile shows one exothermic peak because of the transition phase of the compound between 170 and 180 A degrees C, which was confirmed by X-ray powder diffractometry. The transition phase behavior was studied by DSC curves at several heating rates of a sample mass between 1 and 10 mg in nitrogen atmosphere and in a crucible with and without a lid. Thus, the kinetic parameters were evaluated using an isoconversional non-linear fitting proposed by Capela and Ribeiro. The results show that the activation energy and pre-exponential factor for the transition phase is dependant on the different experimental conditions. Nevertheless, these results indicate that the kinetic compensation effect shows a relationship between them.

Generalização de um modelo linear de impedância eletroquímica

Em geral, a função de um modelo de impedância para processos de eletrodo simples é deduzida de um modelo elétrico equivalente, denominado circuito de Randles. Neste trabalho estudou-se a generalização dessa função, mediante a introdução de um parâmetro não-elétrico, relacionado com a flexibilidade do ângulo de fase e da magnitude. A função foi ajustada às medidas experimentais de impedância obtidas de um sistema constituído de uma liga Ti-10%Al (m/m) em solução de cloreto de sódio 0,9%, variando-se a amplitude de perturbação. Verificou-se que a função generalizada foi adequada para descrever a impedância do sistema analisado, reduzindo as distorções entre a curva experimental e a curva teórica. Além disso, os melhores resultados foram obtidos para sinais de perturbação do sistema com amplitude igual a 10 mV.