Linear and weakly nonlinear stability analyses of cooperative car-following models

Stability analyses have been widely used to better understand the mechanism of traffic jam formation. In this paper, we consider the impact of cooperative systems (a.k.a. connected vehicles) on traffic dynamics and, more precisely, on flow stability. Cooperative systems are emerging technologies enabling communication between vehicles and/or with the infrastructure. In a distributed communication framework, equipped vehicles are able to send and receive information to/from other equipped vehicles. Here, the effects of cooperative traffic are modeled through a general bilateral multianticipative car-following law that improves cooperative drivers' perception of their surrounding traffic conditions within a given communication range. Linear stability analyses are performed for a broad class of car-following models. They point out different stability conditions in both multianticipative and nonmultianticipative situations. To better understand what happens in unstable conditions, information on the shock wave structure is studied in the weakly nonlinear regime by the mean of the reductive perturbation method. The shock wave equation is obtained for generic car-following models by deriving the Korteweg de Vries equations. We then derive traffic-state-dependent conditions for the sign of the solitary wave (soliton) amplitude. This analytical result is verified through simulations. Simulation results confirm the validity of the speed estimate. The variation of the soliton amplitude as a function of the communication range is provided. The performed linear and weakly nonlinear analyses help justify the potential benefits of vehicle-integrated communication systems and provide new insights supporting the future implementation of cooperative systems.

Studies on the meteriological aspects of air pollution over Cochin

Pollutants that once enter into the earth’s atmosphere become part of the atmosphere and hence their dispersion, dilution, direction of transportation etc. are governed by the meteorological conditions. The thesis deals with the study of the atmospheric dispersion capacity, wind climatology, atmospheric stability, pollutant distribution by means of a model and the suggestions for a comprehensive planning for the industrially developing city, Cochin. The definition, sources, types and effects of air pollution have been dealt with briefly. The influence of various meteorological parameters such as vector wind, temperature and its vertical structure and atmospheric stability in relation to pollutant dispersal have been studied. The importance of inversions, mixing heights, ventilation coefficients were brought out. The spatial variation of mixing heights studies for the first time on a microscale region, serves to delineate the regions of good and poor dispersal capacity. A study of wind direction fluctuation, σθ and its relation to stability and mixing heights were shown to be much useful. It was shown that there is a necessity to look into the method of σθ computation. The development of Gausssian Plume Model along with the application for multiple sources was presented. The pollutant chosen was sulphur dioxide and industrial sources alone were considered. The percentage frequency of occurrence of inversions and isothermals are found to be low in all months during the year. The spatial variation of mixing heights revealed that a single mixing height cannot be taken as a representative for the whole city have low mixing heights and monsoonal months showed lowest mixing heights. The study of ventilation co-efficients showed values less than the required optimum value 6000m2/5. However, the low values may be due to the consideration of surface wind alone instead of the vertically averaged wind. Relatively more calm conditions and light winds during night and strong winds during day time were observed. During the most of the year westerlies during day time and northeasterlies during night time are the dominant winds. Unstable conditions with high values of σθ during day time and stable conditions with lower values of σθ during night time are the prominent features. Monsoonal months showed neutral stability for most of the time. A study σθ of and Pasquill Stability category has revealed the difficulty in giving a unique value of for each stability category. For the first time regression equations have been developed relating mixing heights and σθ. A closer examination of σθ revealed that half of the range of wind direction fluctuations is to be taken, instead of one by sixth, to compute σθ. The spatial distribution of SO2 showed a more or less uniform distribution with a slight intrusion towards south. Winter months showed low concentrations contrary to the expectations. The variations of the concentration is found to be influenced more by the mixing height and the stack height rather than wind speed. In the densely populated areas the concentration is more than the threshold limit value. However, the values reported appear to be high, because no depletion of the material is assumed through dry or wet depositions and also because of the inclusion of calm conditions with a very light wind speed. A reduction of emission during night time with a consequent rise during day time would bring down the levels of pollution. The probable locations for the new industries could be the extreme southeast parts because the concentration towards the north falls off very quickly resulting low concentrations. In such a case pollutant spread would be towards south and west, thus keeping the city interior relatively free from pollution. A more detailed examination of the pollutant spread by means of models that would take the dry and wet depositions may be necessary. Nevertheless, the present model serves to give the trend of the distribution of pollutant concentration with which one can suggest the optimum locations for the new industries

Case-to-case variability of predictability of deep convection in a mesoscale model

Successful quantitative precipitation forecasts under convectively unstable conditions depend on the ability of the model to capture the location, timing and intensity of convection. Ensemble forecasts of two mesoscale convective outbreaks over the UK are examined with a view to understanding the nature and extent of their predictability. In addition to a control forecast, twelve ensemble members are run for each case with the same boundary conditions but with perturbations added to the boundary layer. The intention is to introduce perturbations of appropriate magnitude and scale so that the large-scale behaviour of the simulations is not changed. In one case, convection is in statistical equilibrium with the large-scale flow. This places a constraint on the total precipitation, but the location and intensity of individual storms varied. In contrast, the other case was characterised by a large-scale capping inversion. As a result, the location of individual storms was fixed, but their intensities and the total precipitation varied strongly. The ensemble shows case-to-case variability in the nature of predictability of convection in a mesoscale model, and provides additional useful information for quantitative precipitation forecasting.

Formulation of a tropical town energy budget (t-TEB) scheme

This work describes the tropical town energy budget (t-TEB) scheme addressed to simulate the diurnal occurrence of the urban heat island (UHI) as observed in the Metropolitan Area of Rio de Janeiro (MARJ; -22A degrees S; -44A degrees W) in Brazil. Reasoning about the tropical urban climate have guided the scheme implementation, starting from the original equations from Masson (Bound-Lay Meteorol 94:357-397, 2000). The modifications include (a) local scaling approaches for obtaining flux-gradient relationships in the roughness sub-layer, (b) the Monin-Obukhov similarity framework in the inertial sub-layer, (c) increasing aerodynamic conductance toward more unstable conditions, and (d) a modified urban subsurface drainage system to transfer the intercepted rainwater by roofs to the roads. Simulations along 2007 for the MARJ are obtained and compared with the climatology. The t-TEB simulation is consistent with the observations, suggesting that the timing and dynamics of the UHI in tropical cities could vary significantly from the familiar patterns observed in mid-latitude cities-with the peak heat island intensity occurring in the morning than at night. The simulations are suggesting that the thermal phase shift of this tropical diurnal UHI is a response of the surface energy budget to the large amount of solar radiation, intense evapotranspiration, and thermal response of the vegetated surfaces over a very humid soil layer.

Reproduction of a fish assemblage in the state of São Paulo, southeastern Brazil

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

Effects of cold fronts on MODIS-derived sensible and latent heat fluxes in Itumbiara reservoir (Central Brazil)

In this work we investigate the cold front passage effects on sensible and latent heat flux in a tropical hydroelectric reservoir. The study area, Itumbiara reservoir (Goiás State/Brazil) at the beginning of the austral winter, is characterized by the presence of a weak thermal stratification and the passage of several cold fronts from higher latitudes of South America. Sensible and latent heat fluxes were estimated considering the atmospheric boundary layer stability. In situ and MODIS water surface temperature data were used to adjust the coefficients for momentum and heat exchanges between water and atmosphere and spatialize the sensible and latent heat fluxes. The results showed that during a cold front event the sensible heat flux can be up to five times greater than the flux observed before. The latent heat flux tends to decrease during the cold front but increase again after the passage. The highest values of heat loss were observed at littoral zone and some Reservoir's embayment. The heat loss intensification can be separated in two moments: first, during the cold front passage, when the wind speed increases and the air temperature decreases; second, after the cold front passage, with air humidity decreasing. This can be considered a key process to understanding the heat loss in the Itumbiara reservoir. © 2013 COSPAR. Published by Elsevier Ltd. All rights reserved.

Avaliação de características dos regimes de umidade na flona de Caxiuanã-PA durante o experimento COBRA-PARÁ

Procura-se investigar a validade de um método de classificação de regimes de umidade, baseado na caracterização de diferentes "estados" da Camada Limite Atmosférica Tropical (CLAT), acima de uma área de floresta, de acordo com a metodologia proposta por Mahrt (1991). Para essas análises foram utilizados dados de radiossondagens e de uma torre micrometeorológica, coletados durante o período menos chuvoso da região, obtidos durante o experimento "COBRA-PARÁ" (realizado no período de 30/10 a 15/11 de 2006). A análise dos regimes de umidade consiste na representação em espaço de fase dos dados disponíveis da razão de Bowen (β), em função do parâmetro -h/L (onde h é a altura da camada de mistura turbulenta e L é o comprimento de Obukhov). Dependendo da localização dos dados nesse espaço foi possível caracterizar as seguintes classes: classe I - ar seco e instável; classe II - vento seco predominante; classe III - vento úmido; classe IV - condição úmida e instável; classe V - condensação de vapor d'água na superfície; classe VI - condição estável dominante; e classe VII -formação de orvalho induzido por radiação noturna resfriando a superfície. Das classes mencionadas, aquelas mais freqüentemente observadas em Caxiuanã, foram as III, IV e VI.

Calculation of the eigenfunctions of the two-group neutron diffusion equation and application to modal decomposition of BWR instabilities

In this thesis, numerical methods aiming at determining the eigenfunctions, their adjoint and the corresponding eigenvalues of the two-group neutron diffusion equations representing any heterogeneous system are investigated. First, the classical power iteration method is modified so that the calculation of modes higher than the fundamental mode is possible. Thereafter, the Explicitly-Restarted Arnoldi method, belonging to the class of Krylov subspace methods, is touched upon. Although the modified power iteration method is a computationally-expensive algorithm, its main advantage is its robustness, i.e. the method always converges to the desired eigenfunctions without any need from the user to set up any parameter in the algorithm. On the other hand, the Arnoldi method, which requires some parameters to be defined by the user, is a very efficient method for calculating eigenfunctions of large sparse system of equations with a minimum computational effort. These methods are thereafter used for off-line analysis of the stability of Boiling Water Reactors. Since several oscillation modes are usually excited (global and regional oscillations) when unstable conditions are encountered, the characterization of the stability of the reactor using for instance the Decay Ratio as a stability indicator might be difficult if the contribution from each of the modes are not separated from each other. Such a modal decomposition is applied to a stability test performed at the Swedish Ringhals-1 unit in September 2002, after the use of the Arnoldi method for pre-calculating the different eigenmodes of the neutron flux throughout the reactor. The modal decomposition clearly demonstrates the excitation of both the global and regional oscillations. Furthermore, such oscillations are found to be intermittent with a time-varying phase shift between the first and second azimuthal modes.

Age determination, planktic foraminifera and stable istopes of sediment core MSM05/5_712-2

A sediment core from the West Spitsbergen continental margin was studied to reconstruct climate and paleoceanographic variability during the last ~9 ka in the eastern Fram Strait. Our multiproxy evidence suggests that the establishment of the modern oceanographic configuration in the eastern Fram Strait occurred stepwise, in response to the postglacial sea-level rise and the related onset of modern sea-ice production on the shallow Siberian shelves. The late Early and Mid Holocene interval (9 to 5 ka) was generally characterized by relatively unstable conditions. High abundance of the subpolar planktic foraminifer species Turborotalita quinqueloba implies strong intensity of Atlantic Water (AW) inflow with high productivity and/or high AW temperatures, resulting in a strong heat flux to the Arctic. A series of short-lived cooling events (8.2, 6.9. and 6.1 ka) occurred superimposed on the warm late Early and Mid Holocene conditions. Our proxy data imply that simultaneous to the complete postglacial flooding of Arctic shallow shelves and the initiation of modern sea-ice production, strong advance of polar waters initiated modern oceanographic conditions in the eastern Fram Strait at ~5.2 ka. The Late Holocene was marked by the dominance of the polar planktic foraminifer species Neogloboquadrina pachyderma, a significant expansion of sea ice/icebergs, and strong stratification of the water column. Although planktic foraminiferal assemblages as well as sea surface and subsurface temperatures suggest a return of slightly strengthened advection of subsurface Atlantic Water after 3 ka, a relatively stable cold-water layer prevailed at the sea surface and the study site was probably located within the seasonally fluctuating marginal ice zone during the Neoglacial period.

Phenotypic evolution in a fossil gastropod species lineage: evidence for adaptive radiation?

Detecting speciation in the fossil record is a particularly challenging matter. Palaeontologists are usually confronted with poor preservation and limited knowledge on the palaeoenvironment. Even in the contrary case of adequate preservation and information, the linkage of pattern to process is often obscured by insufficient temporal resolution. Consequently, reliable documentations of speciation in fossils with discussions on underlying mechanisms are rare. Here we present a well-resolved pattern of morphological evolution in a fossil species lineage of the gastropod Melanopsis in the long-lived Lake Pannon. These developments are related to environmental changes, documented by isotope and stratigraphical data. After a long period of stasis, the ancestral species experiences a phenotypic change expressed as shift and expansion of the morphospace. The appearance of several new phenotypes along with changes in the environment is interpreted as adaptive radiation. Lake-level high stands affect distribution and availability of habitats and, as a result of varied functional demands on shell geometry, the distribution of phenotypes. The on-going divergence of the morphospace into two branches argues for increasing reproductive isolation, consistent with the model of ecological speciation. In the latest phase, however, progressively unstable conditions restrict availability of niches, allowing survival of one branch only.

Benthic fauna associated with mussel beds and shrimp swarms within hydrothermal fields on the Mid-Atlantic Ridge

Macrofaunal assemblages with prevalence of Bresiliidae shrimps and Mytilidae mussels are abundant in at hydrothermal vents along the Mid-Atlantic Ridge. Mussels inhabit zones of diffuse seeps of hydrothermal fluids with temperature abnormalities up to several degrees. Shrimps inhabit an extreme biotope in a mixed interface between seawater and hydrothermal fluids at temperature up to 20-30°C. We studied the mussel and shrimp assemblages in three hydrothermal vent fields: Rainbow, Broken Spur, and Snake Pit. Species richness of the mussel assemblages within at least two fields (Broken Spur and Snake Pit) is higher as compared with shrimps from the same hydrothermal vent fields. Fauna inhibiting shrimp swarms lack almost any taxa specific for particular assemblages: almost all the taxa are also present in the mussel beds. Structure of the shrimp assemblage is less homogeneous as compared with that of the mussel assemblage. Population prevalence of one taxon (Copepoda) in the shrimp assemblage is most likely connected with extreme and unstable conditions of the biotope occupied by the shrimps in a hydrothermal field. Taxonomic similarity between the mussel and shrimp assemblages within one hydrothermal vent field is higher as compared with similarity between the mussel (or shrimp) assemblages from different fields.

Late Cretaceous planktonic foraminiferal assemblages from the Exmouth Plateau, ODP Sites 122-762 and 122-763, eastern Indian Ocean

The evolution of planktonic foraminifera during the Late Cretaceous is marked in the Santonian by the disappearance of complex morphotypes (the marginotruncanids), and the contemporary increasing importance and diversification of another group of complex taxa, the globotruncanids. Upper Turonian to lower Campanian planktonic foraminiferal assemblages from Holes 762C and 763B (Ocean Drilling Program, Leg 122, Exmouth Plateau, 47°S palaeolatitude) were studied in detail to evaluate the compositional variations at the genus and species level based on the assumption that, in the Cretaceous oceans as in the modern, any faunal change was associated with changes in the characteristics and the degree of stability of the oceanic surface waters. Three major groups were recognised based on gross morphology, and following the assumption that Cretaceous planktonic foraminifera, although extinct, had life-history strategies comparable to those of modern planktonics: 1 - r-selected opportunists; 2 - k-selected specialists; 3 - r/k intermediate morphotypes which include all genera that display a range of trophic strategies in-between opportunist and specialist taxa. Although planktonic foraminiferal assemblages are characterised by a progressive appearance of complex taxa, this trend is discontinuous. Variation in number of species and specimens within genera has allowed recognition of five discrete intervals each of them reflecting different oceanic conditions based on fluctuations in diversity and abundance of the major morphotypes. Planktonic forms show cyclical fluctuations in diversity and abundance of cold (r-strategists) and warm taxa (k-strategists), perhaps representing alternating phases of unstable conditions (suggesting a weakly stratified upper water column in a mesotrophic environment), and well-stratified surface and near-surface waters (indicating a more oligotrophic environment). Interval 1, middle Turonian to early Coniacian in age, is dominated by the r/k intermediate morphotypes which alternate with r-strategists. These cyclical alternations are used to identify three additional subintervals. Interval 2, aged middle to late Coniacian, is characterised by the increasing number of species and relative abundance of k-strategists. After this maximum diversification the k-strategists show a progressive decrease reaching a minimum value in Interval 3 (early to late Santonian), which corresponds to the extinction of the genus Marginotruncana. In the Interval 4, latest Santonian in age, the k-strategists, represented mainly by the genera Globotruncana, increase again in diversity and abundance. The last Interval 5 (early Campanian) is dominated by juvenile globotruncanids and r-strategists which fluctuate in opposite phase. The positive peak (Interval 2) related to the maximum diversification of warm taxa (k-strategists) in the Coniacian seems to correspond to a warmer episode. It is followed by a marked decrease in the relative abundance of warm taxa (k-strategists crisis) with a minimum in the late Santonian (Interval 3), reflecting a decrease in temperature. Detailed analysis of faunal variations allows the Santonian faunal turnover to be ascribed to a cooling event strong enough to cause the extinction of the marginotruncanids.

Predicción de la velocidad crítica de flameo mediante el uso de modelos aeroelásticos extrapolados en tiempo real con datos de vuelo

La Aeroelasticidad fue definida por Arthur Collar en 1947 como "el estudio de la interacción mutua entre fuerzas inerciales, elásticas y aerodinámicas actuando sobre elementos estructurales expuestos a una corriente de aire". Actualmente, esta definición se ha extendido hasta abarcar la influencia del control („Aeroservoelasticidad‟) e, incluso, de la temperatura („Aerotermoelasticidad‟). En el ámbito de la Ingeniería Aeronáutica, los fenómenos aeroelásticos, tanto estáticos (divergencia, inversión de mando) como dinámicos (flameo, bataneo) son bien conocidos desde los inicios de la Aviación. Las lecciones aprendidas a lo largo de la Historia Aeronáutica han permitido establecer criterios de diseño destinados a mitigar la probabilidad de sufrir fenómenos aeroelásticos adversos durante la vida operativa de una aeronave. Adicionalmente, el gran avance experimentado durante esta última década en el campo de la Aerodinámica Computacional y en la modelización aeroelástica ha permitido mejorar la fiabilidad en el cálculo de las condiciones de flameo de una aeronave en su fase de diseño. Sin embargo, aún hoy, los ensayos en vuelo siguen siendo necesarios para validar modelos aeroelásticos, verificar que la aeronave está libre de inestabilidades aeroelásticas y certificar sus distintas envolventes. En particular, durante el proceso de expansión de la envolvente de una aeronave en altitud/velocidad, se requiere predecir en tiempo real las condiciones de flameo y, en consecuencia, evitarlas. A tal efecto, en el ámbito de los ensayos en vuelo, se han desarrollado diversas metodologías que predicen, en tiempo real, las condiciones de flameo en función de condiciones de vuelo ya verificadas como libres de inestabilidades aeroelásticas. De entre todas ellas, aquella que relaciona el amortiguamiento y la velocidad con un parámetro específico definido como „Margen de Flameo‟ (Flutter Margin), permanece como la técnica más común para proceder con la expansión de Envolventes en altitud/velocidad. No obstante, a pesar de su popularidad y facilidad de aplicación, dicha técnica no es adecuada cuando en la aeronave a ensayar se hallan presentes no-linealidades mecánicas como, por ejemplo, holguras. En particular, en vuelos de ensayo dedicados específicamente a expandir la envolvente en altitud/velocidad, las condiciones de „Oscilaciones de Ciclo Límite‟ (Limit Cycle Oscillations, LCOs) no pueden ser diferenciadas de manera precisa de las condiciones de flameo, llevando a una determinación excesivamente conservativa de la misma. La presente Tesis desarrolla una metodología novedosa, basada en el concepto de „Margen de Flameo‟, que permite predecir en tiempo real las condiciones de „Ciclo Límite‟, siempre que existan, distinguiéndolas de las de flameo. En una primera parte, se realiza una revisión bibliográfica de la literatura acerca de los diversos métodos de ensayo existentes para efectuar la expansión de la envolvente de una aeronave en altitud/velocidad, el efecto de las no-linealidades mecánicas en el comportamiento aeroelástico de dicha aeronave, así como una revisión de las Normas de Certificación civiles y militares respecto a este tema. En una segunda parte, se propone una metodología de expansión de envolvente en tiempo real, basada en el concepto de „Margen de Flameo‟, que tiene en cuenta la presencia de no-linealidades del tipo holgura en el sistema aeroelástico objeto de estudio. Adicionalmente, la metodología propuesta se valida contra un modelo aeroelástico bidimensional paramétrico e interactivo programado en Matlab. Para ello, se plantean las ecuaciones aeroelásticas no-estacionarias de un perfil bidimensional en la formulación espacio-estado y se incorpora la metodología anterior a través de un módulo de análisis de señal y otro módulo de predicción. En una tercera parte, se comparan las conclusiones obtenidas con las expuestas en la literatura actual y se aplica la metodología propuesta a resultados experimentales de ensayos en vuelo reales. En resumen, los principales resultados de esta Tesis son: 1. Resumen del estado del arte en los métodos de ensayo aplicados a la expansión de envolvente en altitud/velocidad y la influencia de no-linealidades mecánicas en la determinación de la misma. 2. Revisión de la normas de Certificación Civiles y las normas Militares en relación a la verificación aeroelástica de aeronaves y los límites permitidos en presencia de no-linealidades. 3. Desarrollo de una metodología de expansión de envolvente basada en el Margen de Flameo. 4. Validación de la metodología anterior contra un modelo aeroelástico bidimensional paramétrico e interactivo programado en Matlab/Simulink. 5. Análisis de los resultados obtenidos y comparación con resultados experimentales. ABSTRACT Aeroelasticity was defined by Arthur Collar in 1947 as “the study of the mutual interaction among inertia, elastic and aerodynamic forces when acting on structural elements surrounded by airflow”. Today, this definition has been updated to take into account the Controls („Aeroservoelasticity‟) and even the temperature („Aerothermoelasticity‟). Within the Aeronautical Engineering, aeroelastic phenomena, either static (divergence, aileron reversal) or dynamic (flutter, buzz), are well known since the early beginning of the Aviation. Lessons learned along the History of the Aeronautics have provided several design criteria in order to mitigate the probability of encountering adverse aeroelastic phenomena along the operational life of an aircraft. Additionally, last decade improvements experienced by the Computational Aerodynamics and aeroelastic modelization have refined the flutter onset speed calculations during the design phase of an aircraft. However, still today, flight test remains as a key tool to validate aeroelastic models, to verify flutter-free conditions and to certify the different envelopes of an aircraft. Specifically, during the envelope expansion in altitude/speed, real time prediction of flutter conditions is required in order to avoid them in flight. In that sense, within the flight test community, several methodologies have been developed to predict in real time flutter conditions based on free-flutter flight conditions. Among them, the damping versus velocity technique combined with a Flutter Margin implementation remains as the most common technique used to proceed with the envelope expansion in altitude/airspeed. However, although its popularity and „easy to implement‟ characteristics, several shortcomings can adversely affect to the identification of unstable conditions when mechanical non-linearties, as freeplay, are present. Specially, during test flights devoted to envelope expansion in altitude/airspeed, Limits Cycle Oscillations (LCOs) conditions can not be accurately distinguished from those of flutter and, in consequence, it leads to an excessively conservative envelope determination. The present Thesis develops a new methodology, based on the Flutter Margin concept, that enables in real time the prediction of the „Limit Cycle‟ conditions, whenever they exist, without degrading the capability of predicting the flutter onset speed. The first part of this Thesis presents a review of the state of the art regarding the test methods available to proceed with the envelope expansion of an aircraft in altitude/airspeed and the effect of mechanical non-linearities on the aeroelastic behavior. Also, both civil and military regulations are reviewed with respect aeroelastic investigation of air vehicles. The second part of this Thesis proposes a new methodology to perform envelope expansion in real time based on the Flutter Margin concept when non-linearities, as freeplay, are present. Additionally, this methodology is validated against a Matlab/Slimulink bidimensional aeroelastic model. This model, parametric and interactive, is formulated within the state-space field and it implements the proposed methodology through two main real time modules: A signal processing module and a prediction module. The third part of this Thesis compares the final conclusions derived from the proposed methodology with those stated by the flight test community and experimental results. In summary, the main results provided by this Thesis are: 1. State of the Art review of the test methods applied to envelope expansion in altitude/airspeed and the influence of mechanical non-linearities in its identification. 2. Review of the main civil and military regulations regarding the aeroelastic verification of air vehicles and the limits set when non-linearities are present. 3. Development of a methodology for envelope expansion based on the Flutter Margin concept. 4. A Matlab/Simulink 2D-[aeroelastic model], parametric and interactive, used as a tool to validate the proposed methodology. 5. Conclusions driven from the present Thesis and comparison with experimental results.

Shallow Landslide Risk Assessment of the Westerly Bluffs of the 7.5-minute Freeland Quadrangle, Whidbey Island, Washington

Landslides often occur on slopes rendered unstable by underlying geology, geomorphology, hydrology, weather-climate, slope modifications, or deforestation. Unfortunately, humans commonly exacerbate such unstable conditions through careless or imprudent development practices. Due to local geology, geography, and climatic conditions, Puget Sound of western Washington State is especially landslide-prone. Despite this known issue, detailed analyses of landslide risks for specific communities are few. This study aims to classify areas of high landslide risk on the westerly bluffs of the 7.5 minute Freeland quadrangle based on a combined approach: mapping using LiDAR imagery and the Landform Remote Identification Model (LRIM) to identify landslides, and implementation of the Shallow Slope Stability Model (SHALSTAB) to establish a landslide exceedance probability. The objective is to produce a risk assessment from two shallow landslide scenarios: (1) minimum bluff setback and runout and (2) maximum bluff setback and runout. A simple risk equation that takes into account the probability of hazard occurrence with physical and economic vulnerability (van Westen, 2004) was applied to both scenarios. Results indicate an possible total loss as much as $32.6b from shallow landslides, given a setback of 12 m and a runout of 235 m.

Online Access to the Encyclopaedia Slavica Sanctorum

Encyclopaedia Slavica Sanctorum project aims at building a repertoire of medieval and early modern Bulgarian texts for saints in combination with ethnological data and some visual sources. A basic project task is to produce an accessible on-line digital repository of this valuable cultural heritage treasure. The paper presents the Encyclopaedia Slavica Sanctorum environment, its architecture, functional specification, application modeling process and software implementation. The paper also discusses the specifics of the ―Encyclopaedia Slavica Sanctorum‖ project and its knowledge domain. The paper also presents the integration between the Encyclopaedia Slavica Sanctorum and the Bulgarian Iconographical Digital Library, a digital library keeping rare specimens, private collections of Orthodox icons, wall- paintings and other iconographical objects, selected from difficult-to-access storages, distant churches, chapels, and monasteries, objects in a risk environment or unstable conditions.