Simulation model for driving dynamics, energy use and power supply

Dissertação para obtenção do grau de Mestre em Engenharia Electrotécnica Ramo de Energia

Relativistic time effects in financial dynamics

Financial time series have a complex dynamic nature. Many techniques were adopted having in mind standard paradigms of time flow. This paper explores an alternative route involving relativistic effects. It is observed that the measuring perspective influences the results and that we can have different time textures.

Linear and sliding-mode control design for matrix converter-based unified power flow controllers

This paper presents the design and compares the performance of linear, decoupled and direct power controllers (DPC) for three-phase matrix converters operating as unified power flow controllers (UPFC). A simplified steady-state model of the matrix converter-based UPFC fitted with a modified Venturini high-frequency pulse width modulator is first used to design the linear controllers for the transmission line active (P) and reactive (Q) powers. In order to minimize the resulting cross coupling between P and Q power controllers, decoupled linear controllers (DLC) are synthesized using inverse dynamics linearization. DPC are then developed using sliding-mode control techniques, in order to guarantee both robustness and decoupled control. The designed P and Q power controllers are compared using simulations and experimental results. Linear controllers show acceptable steady-state behaviour but still exhibit coupling between P and Q powers in transient operation. DLC are free from cross coupling but are parameter sensitive. Results obtained by DPC show decoupled power control with zero error tracking and faster responses with no overshoot and no steady-state error. All the designed controllers were implemented using the same digital signal processing hardware.

Simulation and dynamics of freeway traffic

This paper presents the new package entitled Simulator of Intelligent Transportation Systems (SITS) and a computational oriented analysis of traffic dynamics. The SITS adopts a microscopic simulation approach to reproduce real traffic conditions considering different types of vehicles, drivers and roads. A set of experiments with the SITS reveal the dynamic phenomena exhibited by this kind of system. For this purpose a modelling formalism is developed that embeds the statistics and the Laplace transform. The results make possible the adoption of classical system theory tools and point out that it is possible to study traffic systems taking advantage of the knowledge gathered with automatic control algorithms. A complementary perspective for the analysis of the traffic flow is also quantified through the entropy measure.

A method to correct the flow distortion of offshore wind data using CFD simulation and experimental wind tunnel tests

The assessment of wind energy resource for the development of deep offshore wind plants requires the use of every possible source of data and, in many cases, includes data gathered at meteorological stations installed at islands, islets or even oil platforms—all structures that interfere with, and change, the flow characteristics. This work aims to contribute to the evaluation of such changes in the flow by developing a correction methodology and applying it to the case of Berlenga island, Portugal. The study is performed using computational fluid dynamic simulations (CFD) validated by wind tunnel tests. In order to simulate the incoming offshore flow with CFD models a wind profile, unknown a priori, was established using observations from two coastal wind stations and a power law wind profile was fitted to the existing data (a=0.165). The results show that the resulting horizontal wind speed at 80 m above sea level is 16% lower than the wind speed at 80 m above the island for the dominant wind direction sector.

Dinámica y cinética de los procesos físicos y químicos del radical OH con compuestos orgánicos volátiles (COVs) de relevancia atmosférica. Dynamics and kinetics of the physical and chemical processes of OH radicals with atmospheric relevant volatile organic compounds (VOCs)

El objetivo del presente proyecto es estudiar los procesos físicos y químicos del radical OH con compuestos orgánicos volátiles (COVs), con los cuales sea factible la formación de agregados de van der Waals (vdW) responsables de la curvatura en los gráficos de Arrhenius, empleando técnicas modernas, complementarias entre si y novedosas en el país. El problema será abordado desde tres perspectivas complementarias: 1) estudios cinéticos, 2) estudios mecanísticos y de distribución de productos y 3) estudios de la dinámica de los procesos físicos y químicos. La finalidad es alcanzar una mejor comprensión de los mecanismos que intervienen en el comportamiento químico de especies presentes en la atmósfera y obtener datos cinéticos de alta calidad que puedan alimentar modelos computacionales capaces de describir la composición de la atmósfera, presente y futura. Los objetivos son estudiar: 1) mediante fotólisis láser pulsada con detección por fluorescencia inducida por láser (PLP-LIF), en reactores de flujo, la cinética de reacción del radical OH(v”=0) con COVs que presentan gráficos de Arrhenius curvos con energías de activación negativas, tales como alcoholes insaturados, alquenos halogenados, éteres halogenados, ésteres alifáticos; 2) en una cámara de simulación de condiciones atmosféricas de gran volumen (4500 L), la identidad y el rendimiento de productos de las reacciones mencionadas, a fines de evaluar su impacto atmosférico y dilucidar los mecanismos de reacción; 3) mediante haces moleculares y espectroscopía láser, la estructura y reactividad de complejos de vdW entre alcoholes insaturados o aromáticos (cresoles) y el radical OH, como modelo de los aductos propuestos como responsables de la desviación al comportamiento de Arrhenius de las reacciones mencionadas; 4) mediante PLP-LIF y expansiones supersónicas, las constantes específicas estado a estado (ksts) de relajación/reacción del radical OH(v”=1-4) vibracionalmente excitado con los COVs mencionados. Los resultados experimentales obtenidos serán contrastados con cálculos ab-initio de estructura electrónica, los cuales apoyarán las interpretaciones, permitirán proponer estructuras de estados de transición y aductos colisionales, como así también calcular las frecuencias de vibración de los complejos de vdW para su posterior asignación en los espectros LIF y REMPI. Asimismo, los mecanismos de reacción propuestos y los parámetros cinéticos medidos experimentalmente serán comparados con aquellos obtenidos por cálculos teóricos. The aim of this project is to study the physical and chemical processes of OH radicals with volatile organic compounds (VOCs) with which the formation of van der Waals (vdW) clusters, responsible for the observed curvature in the Arrhenius plots, might be feasible. The problem will be addressed as follow : 1) kinetic studies; 2) products distribution and mechanistic studies and 3) dynamical studies of the physical and chemical processes. The purpose is to obtain a better understanding of the mechanisms that govern the chemical behavior of species present in the atmosphere and to obtain high quality kinetic data to be used as input to computational models. We will study: 1) the reaction kinetics of OH (v”=0) radicals with VOCs such as unsaturated alcohols, halogenated alkenes, halogenated ethers, aliphatic esters, which show curved Arrhenius plots and negative activation energies, by PLP-LIF, in flow systems; 2) in a large volume (4500 L) atmospheric simulation chamber, reaction products yields in order to evaluate their atmospheric impact and reaction mechanisms; 3) using molecular beams and laser spectroscopy, the structure and reactivity of the vdW complexes formed between the unsaturated or aromatic alcohols and the OH radicals as a model of the adducts proposed as responsible for the non-Arrhenius behavior; 4) the specific state-to-state relaxation/reaction rate constants (ksts) of the vibrationally excited OH (v”=1-4) radical with the VOCs by PLP-LIF and supersonic expansions. Ab-initio calculations will be carried out to support the interpretation of the experimental results, to obtain the transition state and collisional adducts structures, as well as to calculate the vibrational frequencies of the vdW complexes to assign to the LIF and REMPI spectra. Also, the proposed reaction mechanisms and the experimentally measured kinetic parameters will be compared with those obtained from theoretical calculations.

Optimized venous return with a self-expanding cannula: from computational fluid dynamics to clinical application.

The Smart canula concept allows for collapsed cannula insertion, and self-expansion within a vein of the body. (A) Computational fluid dynamics, and (B) bovine experiments (76+/-3.8 kg) were performed for comparative analyses, prior to (C) the first clinical application. For an 18F access, a given flow of 4 l/min (A) resulted in a pressure drop of 49 mmHg for smart cannula versus 140 mmHg for control. The corresponding Reynolds numbers are 680 versus 1170, respectively. (B) For an access of 28F, the maximal flow for smart cannula was 5.8+/-0.5 l/min versus 4.0+/-0.1 l/min for standard (P<0.0001), for 24F 5.5+/-0.6 l/min versus 3.2+/-0.4 l/min (P<0.0001), and for 20F 4.1+/-0.3 l/min versus 1.6+/-0.3 l/min (P<0.0001). The flow obtained with the smart cannula was 270+/-45% (20F), 172+/-26% (24F), and 134+/-13% (28F) of standard (one-way ANOVA, P=0.014). (C) First clinical application (1.42 m2) with a smart cannula showed 3.55 l/min (100% predicted) without additional fluids. All three assessment steps confirm the superior performance of the smart cannula design.

Time Variation in the Dynamics of Worker Flows: Evidence from the US and Canada

VAR methods have been used to model the inter-relationships between inflows and outfl ows into unemployment and vacancies using tools such as impulse response analysis. In order to investigate whether such impulse responses change over the course of the business cycle or or over time, this paper uses TVP-VARs for US and Canadian data. For the US, we find interesting differences between the most recent recession and earlier recessions and expansions. In particular, we find the immediate effect of a negative shock on both in ow and out flow hazards to be larger in 2008 than in earlier times. Furthermore, the effect of this shock takes longer to decay. For Canada, we fi nd less evidence of time-variation in impulse responses.

High-resolution numerical modelling of flow-vegetation interactions

In this paper, we present and apply a new three-dimensional model for the prediction of canopy-flow and turbulence dynamics in open-channel flow. The approach uses a dynamic immersed boundary technique that is coupled in a sequentially staggered manner to a large eddy simulation. Two different biomechanical models are developed depending on whether the vegetation is dominated by bending or tensile forces. For bending plants, a model structured on the Euler-Bernoulli beam equation has been developed, whilst for tensile plants, an N-pendula model has been developed. Validation against flume data shows good agreement and demonstrates that for a given stem density, the models are able to simulate the extraction of energy from the mean flow at the stem-scale which leads to the drag discontinuity and associated mixing layer.

The dynamics of dengue virus serotype 3 introduction and dispersion in the state of Bahia, Brazil

By 2002, dengue virus serotype 1 (DENV-1) and DENV-2 had circulated for more than a decade in Brazil. In 2002, the introduction of DENV-3 in the state of Bahia produced a massive epidemic and the first cases of dengue hemorrhagic fever. Based on the standardized frequency, timing and location of viral isolations by the state's Central Laboratory, DENV-3 probably entered Bahia through its capital, Salvador, and then rapidly disseminated to other cities, following the main roads. A linear regression model that included traffic flow, distance from the capital and DENV-1 circulation (r² = 0.24, p = 0.001) supported this hypothesis. This pattern was not seen for serotypes already in circulation and was not seen for DENV-3 in the following year. Human population density was another important factor in the intensity of viral circulation. Neither DENV-1 nor DENV-2 fit this model for 2001 or 2003. Since the vector has limited flight range and vector densities fail to correlate with intensity of viral circulation, this distribution represents the movement of infected people and to some extent mosquitoes. This pattern may mimic person-to-person spread of a new infection.

Population dynamics of Triatoma vitticeps (Stål, 1859) in Itanhomi, Minas Gerais, Brazil

Taxonomic markers (head structure morphometry, isoenzymes and randon amplified polymorphism of DNA - RAPD) were used to understand the population dynamics of Triatoma vitticeps, predominant triatomine species in Itanhomi district, using samples obtained from domestic, peridomiciliary and sylvatic habitats. Morphometric analysis revealed sexual dimorphism within the three samples although specimens could not be separated according to the habitat in which they were captured. Forty-two bands were analyzed from RAPD profiles generated using four primers. A dendrogram constructed from Dice's similarity coefficient values showed that migration of the insects between the habitats has occurred, without structuring of populations. Moreover, the dendrogram obtained from the genetic distance values showed an important gene flow between the sylvatic and domestic habitats. No polymorphism was found in the electrophoretic mobility of proteins for the ten enzymes studied. Our results revealed movement of triatomines between the three habitats, suggesting that the presence of T. vitticeps in houses should not be ignored. As invasion of houses by sylvatic insects is frequent and the natural infection indices of this species are among the highest known, epidemiological vigilance studies may reveal possible changes in T. vitticeps behaviour which could present future risks to public health.

Evolution of a debris flow channel monitored using a 3D terrestrial laser scanner

Like numerous torrents in mountainous regions, the Illgraben creek (canton of Wallis, SW Switzerland) produces almost every year several debris flows. The total area of the active catchment is only 4.7 km², but large events ranging from 50'000 to 400'000 m³ are common (Zimmermann 2000). Consequently, the pathway of the main channel often changes suddenly. One single event can for instance fill the whole river bed and dig new several-meters-deep channels somewhere else (Bardou et al. 2003). The quantification of both, the rhythm and the magnitude of these changes, is very important to assess the variability of the bed's cross section and long profile. These parameters are indispensable for numerical modelling, as they should be considered as initial conditions. To monitor the channel evolution an Optech ILRIS 3D terrestrial laser scanner (LIDAR) was used. LIDAR permits to make a complete high precision 3D model of the channel and its surroundings by scanning it from different view points. The 3D data are treated and interpreted with the software Polyworks from Innovmetric Software Inc. Sequential 3D models allow for the determination of the variation in the bed's cross section and long profile. These data will afterwards be used to quantify the erosion and the deposition in the torrent reaches. To complete the chronological evolution of the landforms, precise digital terrain models, obtained by high resolution photogrammetry based on old aerial photographs, will be used. A 500 m long section of the Illgraben channel was scanned on 18th of August 2005 and on 7th of April 2006. These two data sets permit identifying the changes of the channel that occurred during the winter season. An upcoming scanning campaign in September 2006 will allow for the determination of the changes during this summer. Preliminary results show huge variations in the pathway of the Illgraben channel, as well as important vertical and lateral erosion of the river bed. Here we present the results of a river bank on the left (north-western) flank of the channel (Figure 1). For the August 2005 model the scans from 3 viewpoints were superposed, whereas the April 2006 3D image was obtained by combining 5 separate scans. The bank was eroded. The bank got eroded essentially on its left part (up to 6.3 m), where it is hit by the river and the debris flows (Figures 2 and 3). A debris cone has also formed (Figure 3), which suggests that a part of the bank erosion is due to shallow landslides. They probably occur when the river erosion creates an undercut slope. These geometrical data allow for the monitoring of the alluvial dynamics (i.e. aggradation and degradation) on different time scales and the influence of debris flows occurrence on these changes. Finally, the resistance against erosion of the bed's cross section and long profile will be analysed to assess the variability of these two key parameters. This information may then be used in debris flow simulation.

Minimum price variations, time priority and quotes dynamics

We analyze the impact of a minimum price variation (tick) and timepriority on the dynamics of quotes and the trading costs when competitionfor the order flow is dynamic. We find that convergence to competitiveoutcomes can take time and that the speed of convergence is influencedby the tick size, the priority rule and the characteristics of the orderarrival process. We show also that a zero minimum price variation is neveroptimal when competition for the order flow is dynamic. We compare thetrading outcomes with and without time priority. Time priority is shownto guarantee that uncompetitive spreads cannot be sustained over time.However it can sometimes result in higher trading costs. Empiricalimplications are proposed. In particular, we relate the size of thetrading costs to the frequency of new offers and the dynamics of theinside spread to the state of the book.

Direct simulation of interface dynamics: linking capillary pressure, interfacial area and surface energy

We perform direct numerical simulations of drainage by solving Navier- Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and to model the transition from stable flow to viscous fingering, we focus on the definition of macroscopic capillary pressure. When the fluids are at rest, the difference between inlet and outlet pressures and the difference between the intrinsic phase average pressure coincide with the capillary pressure. However, when the fluids are in motion these quantities are dominated by viscous forces. In this case, only a definition based on the variation of the interfacial energy provides an accurate measure of the macroscopic capillary pressure and allows separating the viscous from the capillary pressure components.