Statistical modelling of extreme rainfall in Madeira Island

Extreme rainfall events have triggered a significant number of flash floods in Madeira Island along its past and recent history. Madeira is a volcanic island where the spatial rainfall distribution is strongly affected by its rugged topography. In this thesis, annual maximum of daily rainfall data from 25 rain gauge stations located in Madeira Island were modelled by the generalised extreme value distribution. Also, the hypothesis of a Gumbel distribution was tested by two methods and the existence of a linear trend in both distributions parameters was analysed. Estimates for the 50– and 100–year return levels were also obtained. Still in an univariate context, the assumption that a distribution function belongs to the domain of attraction of an extreme value distribution for monthly maximum rainfall data was tested for the rainy season. The available data was then analysed in order to find the most suitable domain of attraction for the sampled distribution. In a different approach, a search for thresholds was also performed for daily rainfall values through a graphical analysis. In a multivariate context, a study was made on the dependence between extreme rainfall values from the considered stations based on Kendall’s τ measure. This study suggests the influence of factors such as altitude, slope orientation, distance between stations and their proximity of the sea on the spatial distribution of extreme rainfall. Groups of three pairwise associated stations were also obtained and an adjustment was made to a family of extreme value copulas involving the Marshall–Olkin family, whose parameters can be written as a function of Kendall’s τ association measures of the obtained pairs.

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

Modelling the dynamics of dengue real epidemics

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

Structural modelling of Eu3+-based siloxane-poly( oxyethylene) nanohybrids

The modelling of the local structure of sol-gel derived Eu3+-based organic/inorganic hybrids is reported, based on Small-Angle X-ray Scattering (SAXS), photoluminescence and mid-infrared spectroscopy. The hybrid matrix of these organically modified silicates, classed as di-ureasils and termed U(2000) and U(600), is formed by poly( oxyethylene) (POE) chains of variable length grafted to siloxane domains by means of urea cross-linkages. Europium triflate, Eu(CF3SO3)(3), was incorporated in the two di-ureasil matrices with compositions 400 greater than or equal ton greater than or equal to 10, n is the molar ratio of ether oxygens per Eu3+. The SAXS data for undoped hybrids (n=infinity) show the presence of a well-defined peak attributed to the existence of a liquid-like spatial correlation of siloxane rich domains embedded in the polymer matrix and located at the ends of the organic segments. The obtained siloxane particle gyration radius Rg(1) is around 5 Angstrom (error within 10%), whereas the interparticle distance d is 25 +/-2 Angstrom and 40 +/-2 Angstrom, for U(600) and U(2000), respectively. For the Eu3+-based nanocomposites the formation of a two-level hierarchical local structure is discerned. The primary level is constituted by strongly spatially correlated siloxane particles of gyration radius Rg(1) (4-6 and 3-8 Angstrom, errors within 5%, for U(600())n Eu(CF3SO3)(3), 200 greater than or equal ton greater than or equal to 40, and U(2000)(n)Eu(CF3SO3)(3), 400 greater than or equal ton greater than or equal to 40, respectively) forming large clusters of gyration radius Rg(2) (approximate to 75 +/- 10 Angstrom). The local coordination of Eu3+ in both di-ureasil series is described combining the SAXS, photoluminescence and mid-infrared results. In the di-ureasils containing long polymer chains, U(2000)(n)Eu(CF3SO3)(3), the cations interact exclusively with the carbonyl oxygens atoms of the urea bridges at the siloxane-POE interface. In the hybrids containing shorter chains, U(600)(n)Eu(CF3SO3)(3) with n ranging from 200 to 60, the Eu3+ ions interact solely with the ether-type oxygens of the polymer chains. Nevertheless, in this latter family of hybrids a distinct Eu3+ local site environment involving the urea cross-linkages is detected when the europium content is increased up to n=40.

Modelling frequency distributions of 5 minute-averaged solar radiation indexes using Beta probability functions

Five minute-averaged values of sky clearness, direct and diffuse indices, were used to model the frequency distributions of these variables in terms of optical air mass. From more than four years of solar radiation observations it was found that variations in the frequency distributions of the three indices of optical air mass for Botucatu, Brazil, are similar to those in other places, as published in the literature. The proposed models were obtained by linear combination of normalized Beta probability functions, using the observed distributions derived from three years of data. The versatility of these functions allows modelling of all three irradiance indexes to similar levels of accuracy. A comparison with the observed distributions obtained from one year of observations indicate that the models are able to reproduce the observed frequency distributions of all three indices at the 95% confidence level.

Reductions in the Prevalence and Incidence of Geohelminth Infections following a City-wide Sanitation Program in a Brazilian Urban Centre

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

Anaerobic threshold estimation by statistical modelling.

Anaerobic threshold (AT) is usually estimated as a change point problem by visual analysis of the cardiorespiratory response to incremental dynamic exercise. In this study, two phase linear (TPL) models of the linear-linear and linear-quadratic type were used for the estimation of AT. The correlation coefficient between the classical and statistical approaches was 0.88, and 0.89 after outlier exclusion. The TPL models provide a simple method for estimating AT that can be easily implemented using a digital computer for the automatic pattern recognition of AT.

VAr planning using genetic algorithm and linear programming

A combined methodology consisting of successive linear programming (SLP) and a simple genetic algorithm (SGA) solves the reactive planning problem. The problem is divided into operating and planning subproblems; the operating subproblem, which is a nonlinear, ill-conditioned and nonconvex problem, consists of determining the voltage control and the adjustment of reactive sources. The planning subproblem consists of obtaining the optimal reactive source expansion considering operational, economical and physical characteristics of the system. SLP solves the optimal reactive dispatch problem related to real variables, while SGA is used to determine the necessary adjustments of both the binary and discrete variables existing in the modelling problem. Once the set of candidate busbars has been defined, the program implemented gives the location and size of the reactive sources needed, if any, to maintain the operating and security constraints.

On the analytical solution for the slewing flexible beam-like system: analysis of the linear part of the perturbed problem

The dynamical system investigated in this work is a nonlinear flexible beam-like structure in slewing motion. Non-dimensional and perturbed governing equations of motion are presented. The analytical solution for the linear part of these perturbed equations for ideal and for non-ideal cases are obtained. This solution is necessary for the investigation of the complete weak nonlinear problem where all nonlinearities are small perturbations around a linear known solution. This investigation shall help the analyst in the modelling of dynamical systems with structure- actuator interactions.

Bilevel approach for optimal location and contract pricing of distributed generation in radial distribution systems using mixed-integer linear programming

In this study, a novel approach for the optimal location and contract pricing of distributed generation (DG) is presented. Such an approach is designed for a market environment in which the distribution company (DisCo) can buy energy either from the wholesale energy market or from the DG units within its network. The location and contract pricing of DG is determined by the interaction between the DisCo and the owner of the distributed generators. The DisCo intends to minimise the payments incurred in meeting the expected demand, whereas the owner of the DG intends to maximise the profits obtained from the energy sold to the DisCo. This two-agent relationship is modelled in a bilevel scheme. The upper-level optimisation is for determining the allocation and contract prices of the DG units, whereas the lower-level optimisation is for modelling the reaction of the DisCo. The bilevel programming problem is turned into an equivalent single-level mixed-integer linear optimisation problem using duality properties, which is then solved using commercially available software. Results show the robustness and efficiency of the proposed model compared with other existing models. As regards to contract pricing, the proposed approach allowed to find better solutions than those reported in previous works. © The Institution of Engineering and Technology 2013.

Formulação não linear de ração para frangos de corte: modelagem e otimização

Pós-graduação em Ciência Animal - FMVA

Comportamento dinâmico de suspensão eletromecânica: análise da resposta devido ao amortecimento e à rigidez não linear

Pós-graduação em Engenharia Elétrica - FEB

Modelação de redes de distribuição de energia elétrica com carregamento não linear, incluindo elementos de regulação e compensação

Pós-graduação em Engenharia Elétrica - FEIS

Short-term expansion planning of radial electrical distribution systems using mixed-integer linear programming

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

Estimating the development of the fennel aphid, Hyadaphis foeniculi (Passerini) (Hemiptera: Aphiididae), using non-linear models

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