Applying the extended Kalman filter to systems described by nonlinear differential-algebraic equations

This paper describes a method for the state estimation of nonlinear systems described by a class of differential-algebraic equation models using the extended Kalman filter. The method involves the use of a time-varying linearisation of a semi-explicit index one differential-algebraic equation. The estimation technique consists of a simplified extended Kalman filter that is integrated with the differential-algebraic equation model. The paper describes a simulation study using a model of a batch chemical reactor. It also reports a study based on experimental data obtained from a mixing process, where the model of the system is solved using the sequential modular method and the estimation involves a bank of extended Kalman filters.

Crossbred cow adoption and milk market participation in a multivariate count data framework

Cross-bred cow adoption is an important and potent policy variable precipitating subsistence household entry into emerging milk markets. This paper focuses on the problem of designing policies that encourage and sustain milkmarket expansion among a sample of subsistence households in the Ethiopian highlands. In this context it is desirable to measure households’ ‘proximity’ to market in terms of the level of deficiency of essential inputs. This problem is compounded by four factors. One is the existence of cross-bred cow numbers (count data) as an important, endogenous decision by the household; second is the lack of a multivariate generalization of the Poisson regression model; third is the censored nature of the milk sales data (sales from non-participating households are, essentially, censored at zero); and fourth is an important simultaneity that exists between the decision to adopt a cross-bred cow, the decision about how much milk to produce, the decision about how much milk to consume and the decision to market that milk which is produced but not consumed internally by the household. Routine application of Gibbs sampling and data augmentation overcome these problems in a relatively straightforward manner. We model the count data from two sites close to Addis Ababa in a latent, categorical-variable setting with known bin boundaries. The single-equation model is then extended to a multivariate system that accommodates the covariance between crossbred-cow adoption, milk-output, and milk-sales equations. The latent-variable procedure proves tractable in extension to the multivariate setting and provides important information for policy formation in emerging-market settings

A dynamical systems framework for intermittent data assimilation

We consider the problem of discrete time filtering (intermittent data assimilation) for differential equation models and discuss methods for its numerical approximation. The focus is on methods based on ensemble/particle techniques and on the ensemble Kalman filter technique in particular. We summarize as well as extend recent work on continuous ensemble Kalman filter formulations, which provide a concise dynamical systems formulation of the combined dynamics-assimilation problem. Possible extensions to fully nonlinear ensemble/particle based filters are also outlined using the framework of optimal transportation theory.

Intercomparison of methods of coupling between convection and large-scale circulation: 1. Comparison over uniform surface conditions

As part of an international intercomparison project, a set of single column models (SCMs) and cloud-resolving models (CRMs) are run under the weak temperature gradient (WTG) method and the damped gravity wave (DGW) method. For each model, the implementation of the WTG or DGW method involves a simulated column which is coupled to a reference state defined with profiles obtained from the same model in radiative-convective equilibrium. The simulated column has the same surface conditions as the reference state and is initialized with profiles from the reference state. We performed systematic comparison of the behavior of different models under a consistent implementation of the WTG method and the DGW method and systematic comparison of the WTG and DGW methods in models with different physics and numerics. CRMs and SCMs produce a variety of behaviors under both WTG and DGW methods. Some of the models reproduce the reference state while others sustain a large-scale circulation which results in either substantially lower or higher precipitation compared to the value of the reference state. CRMs show a fairly linear relationship between precipitation and circulation strength. SCMs display a wider range of behaviors than CRMs. Some SCMs under the WTG method produce zero precipitation. Within an individual SCM, a DGW simulation and a corresponding WTG simulation can produce different signed circulation. When initialized with a dry troposphere, DGW simulations always result in a precipitating equilibrium state. The greatest sensitivities to the initial moisture conditions occur for multiple stable equilibria in some WTG simulations, corresponding to either a dry equilibrium state when initialized as dry or a precipitating equilibrium state when initialized as moist. Multiple equilibria are seen in more WTG simulations for higher SST. In some models, the existence of multiple equilibria is sensitive to some parameters in the WTG calculations.

Regulation of early steps of GPVI signal transduction by phosphatases: a systems biology approach

We present a data-driven mathematical model of a key initiating step in platelet activation, a central process in the prevention of bleeding following Injury. In vascular disease, this process is activated inappropriately and causes thrombosis, heart attacks and stroke. The collagen receptor GPVI is the primary trigger for platelet activation at sites of injury. Understanding the complex molecular mechanisms initiated by this receptor is important for development of more effective antithrombotic medicines. In this work we developed a series of nonlinear ordinary differential equation models that are direct representations of biological hypotheses surrounding the initial steps in GPVI-stimulated signal transduction. At each stage model simulations were compared to our own quantitative, high-temporal experimental data that guides further experimental design, data collection and model refinement. Much is known about the linear forward reactions within platelet signalling pathways but knowledge of the roles of putative reverse reactions are poorly understood. An initial model, that includes a simple constitutively active phosphatase, was unable to explain experimental data. Model revisions, incorporating a complex pathway of interactions (and specifically the phosphatase TULA-2), provided a good description of the experimental data both based on observations of phosphorylation in samples from one donor and in those of a wider population. Our model was used to investigate the levels of proteins involved in regulating the pathway and the effect of low GPVI levels that have been associated with disease. Results indicate a clear separation in healthy and GPVI deficient states in respect of the signalling cascade dynamics associated with Syk tyrosine phosphorylation and activation. Our approach reveals the central importance of this negative feedback pathway that results in the temporal regulation of a specific class of protein tyrosine phosphatases in controlling the rate, and therefore extent, of GPVI-stimulated platelet activation.

Intercomparison of methods of coupling between convection and large-scale circulation. 2: Comparison over non-uniform surface conditions

As part of an international intercomparison project, the weak temperature gradient (WTG) and damped gravity wave (DGW) methods are used to parameterize large-scale dynamics in a set of cloud-resolving models (CRMs) and single column models (SCMs). The WTG or DGW method is implemented using a configuration that couples a model to a reference state defined with profiles obtained from the same model in radiative-convective equilibrium. We investigated the sensitivity of each model to changes in SST, given a fixed reference state. We performed a systematic comparison of the WTG and DGW methods in different models, and a systematic comparison of the behavior of those models using the WTG method and the DGW method. The sensitivity to the SST depends on both the large-scale parameterization method and the choice of the cloud model. In general, SCMs display a wider range of behaviors than CRMs. All CRMs using either the WTG or DGW method show an increase of precipitation with SST, while SCMs show sensitivities which are not always monotonic. CRMs using either the WTG or DGW method show a similar relationship between mean precipitation rate and column-relative humidity, while SCMs exhibit a much wider range of behaviors. DGW simulations produce large-scale velocity profiles which are smoother and less top-heavy compared to those produced by the WTG simulations. These large-scale parameterization methods provide a useful tool to identify the impact of parameterization differences on model behavior in the presence of two-way feedback between convection and the large-scale circulation.

Is bird incidence in Atlantic forest fragments influenced by landscape patterns at multiple scales?

The degree to which habitat fragmentation affects bird incidence is species specific and may depend on varying spatial scales. Selecting the correct scale of measurement is essential to appropriately assess the effects of habitat fragmentation on bird occurrence. Our objective was to determine which spatial scale of landscape measurement best describes the incidence of three bird species (Pyriglena leucoptera, Xiphorhynchus fuscus and Chiroxiphia caudata) in the fragmented Brazilian Atlantic forest and test if multi-scalar models perform better than single-scalar ones. Bird incidence was assessed in 80 forest fragments. The surrounding landscape structure was described with four indices measured at four spatial scales (400-, 600-, 800- and 1,000-m buffers around the sample points). The explanatory power of each scale in predicting bird incidence was assessed using logistic regression, bootstrapped with 1,000 repetitions. The best results varied between species (1,000-m radius for P. leucoptera; 800-m for X. fuscus and 600-m for C. caudata), probably due to their distinct feeding habits and foraging strategies. Multi-scale models always resulted in better predictions than single-scale models, suggesting that different aspects of the landscape structure are related to different ecological processes influencing bird incidence. In particular, our results suggest that local extinction and (re)colonisation processes might simultaneously act at different scales. Thus, single-scale models may not be good enough to properly describe complex pattern-process relationships. Selecting variables at multiple ecologically relevant scales is a reasonable procedure to optimise the accuracy of species incidence models.

Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows

In this paper we present a finite difference method for solving two-dimensional viscoelastic unsteady free surface flows governed by the single equation version of the eXtended Pom-Pom (XPP) model. The momentum equations are solved by a projection method which uncouples the velocity and pressure fields. We are interested in low Reynolds number flows and, to enhance the stability of the numerical method, an implicit technique for computing the pressure condition on the free surface is employed. This strategy is invoked to solve the governing equations within a Marker-and-Cell type approach while simultaneously calculating the correct normal stress condition on the free surface. The numerical code is validated by performing mesh refinement on a two-dimensional channel flow. Numerical results include an investigation of the influence of the parameters of the XPP equation on the extrudate swelling ratio and the simulation of the Barus effect for XPP fluids. (C) 2010 Elsevier B.V. All rights reserved.

Compartilhamento do conhecimento e confiança organizacional em ambientes fortemente hierarquizados

O compartilhamento do conhecimento e a confiança organizacional são fatores de grande interesse nas pesquisas sobre gestão do conhecimento. Esta tese tem por objetivo identificar como a confiança organizacional influencia a propensão ao compartilhamento do conhecimento em estruturas hierárquicas fortes, estudando os efeitos da satisfação com a carreira, do comprometimento organizacional e do tempo de exposição à hierarquia nesse processo. O método hipotético-dedutivo, aplicado com a técnica de modelagem de equações estruturais a uma amostra de 655 profissionais militares do Exército Brasileiro resultou na mediação do comprometimento organizacional afetivo no relacionamento entre a confiança organizacional e a propensão ao compartilhamento do conhecimento. Os resultados sugerem, ainda, que a percepção de utilidade do conhecimento recebido e o estado civil são variáveis significativas na explicação da variância da propensão ao compartilhamento do conhecimento. Por fim, o tempo de exposição à hierarquia impacta diretamente as variáveis estudadas sem, contudo, interferir no relacionamento entre o comprometimento organizacional e a propensão ao compartilhamento do conhecimento. Os resultados desta tese contribuem para o melhor entendimento do fenômeno de compartilhamento do conhecimento no ambiente organizacional.

Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows

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

Spontaneous symmetry breaking of Bose-Fermi mixtures in double-well potentials

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

Biometric relationship analysis of the wet weight by carapace width to the swimming crab Arenaeus-cribrarius (Lamarck, 1818) (Crustacea, Brachyura, Portunidae)

The biometric relationship between the weigth and the carapace width in the swimming crab A. cribrarius was compared with the results from other portunid crabs studied previously. During November/1988 to October/1989, a total of 403 specimens (189 males and 214 females) were collected with otter-trawl nets in north coast of the São Paulo State, Brazil. The animals were measured (carapace width excluding lateral spines = LC! and weighed (wet weight = PE). The empiric points of this relation were fit according to the power function (Y = a.X(b)) for each sex, maturation phases and total of individuals. The relation PE x LC indicates that the mole's growth changes during the ontogenesys from isometric (in juvenile phase) to allometric positive (in adult phase). For the females the growth is isometric in the two phases. The weight grows in a higher proportion than the carapace width variable (allometric positive growth). The data can be grouped in a single equation (PE = 7.85.10(-5).LC(3.14)) for the convertion between the variables there was a greater similarity between the equations obtained far each sex. In spite of this, the males present the fattening grade value (''a'') slightly higher than that of the females, possibly because of the greater size reached in its devellopment. The mean weight of the males is greater than the females one (p < 0.01). In the range 80 proves 90mm the males were more abundant, probably due to the females terminal ecdysis is near this size. The females only have the mean weight greater than the males in the 60 proves 70mm range (p < 0.01) when the puberty molt occurs and they present morphological changes in their reproductive system.

Turbulence modeling in engineering applications

A brief review is given of turbulence models in use today for engineering applications. The main categories covered are simple eddy-viscosity models, the k-ε two-equation model and Reynolds-stress-equation models as well as their algebraic stress derivatives. Calculation examples are presented for a variety of 2D and 3D flows.

Giant enhancement of phonon-assisted one-photon excited frequency upconversion in a Nd3+-doped tellurite glass

Changing the sample's temperature from 200 K to 535 K, we observed 670-fold enhancement of a phonon-assisted upconversion emission at ≈754 nm obtained from a Nd3+-doped tellurite glass excited by 5 ns laser pulses at 805 nm. A rate-equation model, including the relevant energy levels and temperature dependent transition rates, is proposed to describe the process. The results fit well with the data when one considers the nonradiative transitions contributing for the 754 nm luminescence are promoted by an effective phonon mode with energy of 700 cm-1. © 2013 American Institute of Physics.

The impact of climate change on the tourism sector in the Caribbean

The main aim of this study is to estimate the economic impact of climate change on nine countries in the Caribbean basin: Aruba, Barbados, Dominican Republic, Guyana, Jamaica, Montserrat, Netherlands Antilles, Saint Lucia and Trinidad and Tobago. A typical tourism demand function, with tourist arrivals as the dependent variable, is used in the analysis. To establish the baseline, the period under analysis is 1989-2007 and the independent variables are destination country GDP per capita and consumer price index, source country GDP, oil prices to proxy transportation costs between source and destination countries. At this preliminary stage the climate variables are used separately to augment the tourism demand function to establish a relationship, if any, among the variables. Various econometric models (single OLS models for each country, pooled regression, GMM estimation and random effects panel models) were considered in an attempt to find the best way to model the data. The best fit for the data (1989-2007) is the random effects panel data model augmented by both climate variables, i.e. temperature and precipitation. Projections of all variables in the model for the 2008-2100 period were done using forecasting techniques. Projections for the climate variables were undertaken by INSMET. The cost of climate change to the tourism sector was estimated under three scenarios: A2, B2 and BAU (the mid-point of the A2 and B2 scenarios). The estimated costs to tourism for the Caribbean subregion under the three scenarios are all very high and ranges from US$43.9 billion under the B2 scenario to US$46.3 billion under the BAU scenario.