CFD Modelling of Direct Contact Condensation in Suppression Pools by Applying Condensation Models of Separated Flow

The condensation rate has to be high in the safety pressure suppression pool systems of Boiling Water Reactors (BWR) in order to fulfill their safety function. The phenomena due to such a high direct contact condensation (DCC) rate turn out to be very challenging to be analysed either with experiments or numerical simulations. In this thesis, the suppression pool experiments carried out in the POOLEX facility of Lappeenranta University of Technology were simulated. Two different condensation modes were modelled by using the 2-phase CFD codes NEPTUNE CFD and TransAT. The DCC models applied were the typical ones to be used for separated flows in channels, and their applicability to the rapidly condensing flow in the condensation pool context had not been tested earlier. A low Reynolds number case was the first to be simulated. The POOLEX experiment STB-31 was operated near the conditions between the ’quasi-steady oscillatory interface condensation’ mode and the ’condensation within the blowdown pipe’ mode. The condensation models of Lakehal et al. and Coste & Lavi´eville predicted the condensation rate quite accurately, while the other tested ones overestimated it. It was possible to get the direct phase change solution to settle near to the measured values, but a very high resolution of calculation grid was needed. Secondly, a high Reynolds number case corresponding to the ’chugging’ mode was simulated. The POOLEX experiment STB-28 was chosen, because various standard and highspeed video samples of bubbles were recorded during it. In order to extract numerical information from the video material, a pattern recognition procedure was programmed. The bubble size distributions and the frequencies of chugging were calculated with this procedure. With the statistical data of the bubble sizes and temporal data of the bubble/jet appearance, it was possible to compare the condensation rates between the experiment and the CFD simulations. In the chugging simulations, a spherically curvilinear calculation grid at the blowdown pipe exit improved the convergence and decreased the required cell count. The compressible flow solver with complete steam-tables was beneficial for the numerical success of the simulations. The Hughes-Duffey model and, to some extent, the Coste & Lavi´eville model produced realistic chugging behavior. The initial level of the steam/water interface was an important factor to determine the initiation of the chugging. If the interface was initialized with a water level high enough inside the blowdown pipe, the vigorous penetration of a water plug into the pool created a turbulent wake which invoked the chugging that was self-sustaining. A 3D simulation with a suitable DCC model produced qualitatively very realistic shapes of the chugging bubbles and jets. The comparative FFT analysis of the bubble size data and the pool bottom pressure data gave useful information to distinguish the eigenmodes of chugging, bubbling, and pool structure oscillations.

Management zones using fuzzy clustering based on spatial-temporal variability of soil and corn yield

Clustering soil and crop data can be used as a basis for the definition of management zones because the data are grouped into clusters based on the similar interaction of these variables. Therefore, the objective of this study was to identify management zones using fuzzy c-means clustering analysis based on the spatial and temporal variability of soil attributes and corn yield. The study site (18 by 250-m in size) was located in Jaboticabal, São Paulo/Brazil. Corn yield was measured in one hundred 4.5 by 10-m cells along four parallel transects (25 observations per transect) over five growing seasons between 2001 and 2010. Soil chemical and physical attributes were measured. SAS procedure MIXED was used to identify which variable(s) most influenced the spatial variability of corn yield over the five study years. Basis saturation (BS) was the variable that better related to corn yield, thus, semivariograms models were fitted for BS and corn yield and then, data values were krigged. Management Zone Analyst software was used to carry out the fuzzy c-means clustering algorithm. The optimum number of management zones can change over time, as well as the degree of agreement between the BS and corn yield management zone maps. Thus, it is very important take into account the temporal variability of crop yield and soil attributes to delineate management zones accurately.

Modeling temporal variations of Gracilaria Greville and Hypnea J.V. Lamouroux (Rhodophyta) assemblages on a midlittoral reef platform at Piedade Beach, Pernambuco State, Brazil

The diversity of algal banks composed of species out the genera Gracilaria Greville and Hypnea J.V. Lamouroux have been impacted by commercial exploitation and coastal eutrophication. The present study sought to construct dynamic models based on algal physiology to simulate seasonal variations in the biomasses of Gracilaria and Hypnea an intertidal reef at Piedade Beach in Jaboatão dos Guararapes, Pernambuco State, Brazil. Five 20 × 20 cm plots in a reef pool on a midlittoral reef platform were randomly sampled during April, June, August, October, and December/2009 and in January and March/2010. Water temperature, pH, irradiance, oxygen and salinity levels as well as the concentrations of ammonia, nitrate and phosphate were measured at the sampling site. Forcing functions were employed in the model to represent abiotic factors, and algal decay was simulated with a dispersal function. Algal growth was modeled using a logistic function and was found to be sensitive to temperature and salinity. Maximum absorption rates of ammonia and phosphate were higher in Hypnea than in Gracilaria, indicating that the former takes up nutrients more efficiently at higher concentrations. Gracilaria biomass peaked at approximately 120 g (dry weight m-2) in March/2010 and was significantly lower in August/2009; Hypnea biomasses, on the other hand, did not show any significant variations among the different months, indicating that resource competition may influence the productivity of these algae.

Hydraulic and fluvial geomorphological models for a bedrock channel reach of the Twenty Mile Creek /

Bedrock channels have been considered challenging geomorphic settings for the application of numerical models. Bedrock fluvial systems exhibit boundaries that are typically less mobile than alluvial systems, yet they are still dynamic systems with a high degree of spatial and temporal variability. To understand the variability of fluvial systems, numerical models have been developed to quantify flow magnitudes and patterns as the driving force for geomorphic change. Two types of numerical model were assessed for their efficacy in examining the bedrock channel system consisting of a high gradient portion of the Twenty Mile Creek in the Niagara Region of Ontario, Canada. A one-dimensional (1-D) flow model that utilizes energy equations, HEC RAS, was used to determine velocity distributions through the study reach for the mean annual flood (MAF), the 100-year return flood and the 1,000-year return flood. A two-dimensional (2-D) flow model that makes use of Navier-Stokes equations, RMA2, was created with the same objectives. The 2-D modeling effort was not successful due to the spatial complexity of the system (high slope and high variance). The successful 1 -D model runs were further extended using very high resolution geospatial interpolations inherent to the HEC RAS extension, HEC geoRAS. The modeled velocity data then formed the basis for the creation of a geomorphological analysis that focused upon large particles (boulders) and the forces needed to mobilize them. Several existing boulders were examined by collecting detailed measurements to derive three-dimensional physical models for the application of fluid and solid mechanics to predict movement in the study reach. An imaginary unit cuboid (1 metre by 1 metre by 1 metre) boulder was also envisioned to determine the general propensity for the movement of such a boulder through the bedrock system. The efforts and findings of this study provide a standardized means for the assessment of large particle movement in a bedrock fluvial system. Further efforts may expand upon this standardization by modeling differing boulder configurations (platy boulders, etc.) at a high level of resolution.

Inter-temporal Price Discrimination when Imports are Restricted by Quotas

A dominant firm holding import quota engages in inter-temporal price discrimination when facing a competitive fringe engaged in seasonal production. This causes a welfare loss that comes in addition the loss attributable to limitation of imports below the free trade level.

Bootstrap for panel data models with an application to the evaluation of public policies

Le but de cette thèse est d étendre la théorie du bootstrap aux modèles de données de panel. Les données de panel s obtiennent en observant plusieurs unités statistiques sur plusieurs périodes de temps. Leur double dimension individuelle et temporelle permet de contrôler l 'hétérogénéité non observable entre individus et entre les périodes de temps et donc de faire des études plus riches que les séries chronologiques ou les données en coupe instantanée. L 'avantage du bootstrap est de permettre d obtenir une inférence plus précise que celle avec la théorie asymptotique classique ou une inférence impossible en cas de paramètre de nuisance. La méthode consiste à tirer des échantillons aléatoires qui ressemblent le plus possible à l échantillon d analyse. L 'objet statitstique d intérêt est estimé sur chacun de ses échantillons aléatoires et on utilise l ensemble des valeurs estimées pour faire de l inférence. Il existe dans la littérature certaines application du bootstrap aux données de panels sans justi cation théorique rigoureuse ou sous de fortes hypothèses. Cette thèse propose une méthode de bootstrap plus appropriée aux données de panels. Les trois chapitres analysent sa validité et son application. Le premier chapitre postule un modèle simple avec un seul paramètre et s 'attaque aux propriétés théoriques de l estimateur de la moyenne. Nous montrons que le double rééchantillonnage que nous proposons et qui tient compte à la fois de la dimension individuelle et la dimension temporelle est valide avec ces modèles. Le rééchantillonnage seulement dans la dimension individuelle n est pas valide en présence d hétérogénéité temporelle. Le ré-échantillonnage dans la dimension temporelle n est pas valide en présence d'hétérogénéité individuelle. Le deuxième chapitre étend le précédent au modèle panel de régression. linéaire. Trois types de régresseurs sont considérés : les caractéristiques individuelles, les caractéristiques temporelles et les régresseurs qui évoluent dans le temps et par individu. En utilisant un modèle à erreurs composées doubles, l'estimateur des moindres carrés ordinaires et la méthode de bootstrap des résidus, on montre que le rééchantillonnage dans la seule dimension individuelle est valide pour l'inférence sur les coe¢ cients associés aux régresseurs qui changent uniquement par individu. Le rééchantillonnage dans la dimen- sion temporelle est valide seulement pour le sous vecteur des paramètres associés aux régresseurs qui évoluent uniquement dans le temps. Le double rééchantillonnage est quand à lui est valide pour faire de l inférence pour tout le vecteur des paramètres. Le troisième chapitre re-examine l exercice de l estimateur de différence en di¤érence de Bertrand, Duflo et Mullainathan (2004). Cet estimateur est couramment utilisé dans la littérature pour évaluer l impact de certaines poli- tiques publiques. L exercice empirique utilise des données de panel provenant du Current Population Survey sur le salaire des femmes dans les 50 états des Etats-Unis d Amérique de 1979 à 1999. Des variables de pseudo-interventions publiques au niveau des états sont générées et on s attend à ce que les tests arrivent à la conclusion qu il n y a pas d e¤et de ces politiques placebos sur le salaire des femmes. Bertrand, Du o et Mullainathan (2004) montre que la non-prise en compte de l hétérogénéité et de la dépendance temporelle entraîne d importantes distorsions de niveau de test lorsqu'on évalue l'impact de politiques publiques en utilisant des données de panel. Une des solutions préconisées est d utiliser la méthode de bootstrap. La méthode de double ré-échantillonnage développée dans cette thèse permet de corriger le problème de niveau de test et donc d'évaluer correctement l'impact des politiques publiques.

Spatio-temporal Variability of Atmospheric Ozone over Indian Subcontinent and its Relation to Meteorological Parameters

Ozone present in the atmosphere not only absorbs the biologically harmful ultraviolet radiation but also is an important ingredient of the climate system. The radiative absorption properties of ozone make it a determining factor in the structure of the atmosphere. Ozone in the troposphere has many negative impacts on humans and other living beings. Another significant aspect is the absorption of outgoing infrared radiation by ozone thus acting as a greenhouse gas. The variability of ozone in the atmosphere involves many interconnections with the incoming and outgoing radiation, temperature circulation etc. Hence ozone forms an important part of chemistry-climate as well as radiative transfer models. This aspect also makes the quantification of ozone more important. The discovery of Antarctic ozone hole and the role of anthropogenic activities in causing it made it possible to plan and implement necessary preventive measures. Continuous monitoring of ozone is also necessary to identify the effect of these preventive steps. The reactions involving the formation and destruction of ozone are influenced significantly by the temperature fluctuations of the atmosphere. On the other hand the variations in ozone can change the temperature structure of the atmosphere. Indian subcontinent is a region having large weather and climate variability which is evident from the large interannual variability of monsoon system over the region. Nearly half of Indian region comprises the tropical region. Most of ozone is formed in the tropical region and transported to higher latitudes. The formation and transport of ozone can be influenced by changes in solar radiation and various atmospheric circulation features. Besides industrial activities and vehicular traffic is more due to its large population. This may give rise to an increase in the production of tropospheric ozone which is greenhouse gas. Hence it becomes necessary to monitor the atmospheric ozone over this region. This study probes into the spatial distribution and temporal evolution of ozone over Indian subcontinent and discusses the contributing atmospheric parameters.

Generative Temporal Planning with Complex Processes

Autonomous vehicles are increasingly being used in mission-critical applications, and robust methods are needed for controlling these inherently unreliable and complex systems. This thesis advocates the use of model-based programming, which allows mission designers to program autonomous missions at the level of a coach or wing commander. To support such a system, this thesis presents the Spock generative planner. To generate plans, Spock must be able to piece together vehicle commands and team tactics that have a complex behavior represented by concurrent processes. This is in contrast to traditional planners, whose operators represent simple atomic or durative actions. Spock represents operators using the RMPL language, which describes behaviors using parallel and sequential compositions of state and activity episodes. RMPL is useful for controlling mobile autonomous missions because it allows mission designers to quickly encode expressive activity models using object-oriented design methods and an intuitive set of activity combinators. Spock also is significant in that it uniformly represents operators and plan-space processes in terms of Temporal Plan Networks, which support temporal flexibility for robust plan execution. Finally, Spock is implemented as a forward progression optimal planner that walks monotonically forward through plan processes, closing any open conditions and resolving any conflicts. This thesis describes the Spock algorithm in detail, along with example problems and test results.

Front spreading in population dynamic models. Theory and application to the Neolithic transition

This thesis presents population dynamics models that can be applied to predict the rate of spread of the Neolithic transition (change from hunter-gathering to farming economics) across the European continent, which took place about 9000 to 5000 years ago. The first models in this thesis provide predictions at a continental scale. We develop population dynamics models with explicit kernels and apply realistic data. We also derive a new time-delayed reaction-diffusion equation which yields speeds about a 10% slower than previous models. We also deal with a regional variability: the slowdown of the Neolithic front when reaching the North of Europe. We develop simple reaction-diffusion models that can predict the measured speeds in terms of the non-homogeneous distribution of pre-Neolithic (Mesolithic) population in Europe, which were present in higher densities at the North of the continent. Such models can explain the observed speeds.

How well do high resolution models reproduce tropical convection?

Cascade is a multi-institution project studying the temporal and spatial organization of tropical convective systems. While cloud resolving numerical models can reproduce the observed diurnal cycle of such systems they are sensitive to the chosen resolution. As part of this effort, we are comparing results from the Met. Office Unified Model to data from the Global Earth Radiation Budget satellite instrument over the African Monsoon Interdisciplinary Analyses region of North Africa. We use a variety of mathematical techniques to study the outgoing radiation and the evolution of properties such as the cloud size distribution. The effectiveness of various model resolutions is tested with a view to determining the optimum balance between resolution and the need to reproduce the observations.

Evaluation of two hybrid metric-conceptual models, for simulating phosphorus transfer from agricultural land in the river enborne, a lowland UK catchment

Across Europe, elevated phosphorus (P) concentrations in lowland rivers have made them particularly susceptible to eutrophication. This is compounded in southern and central UK by increasing pressures on water resources, which may be further enhanced by the potential effects of climate change. The EU Water Framework Directive requires an integrated approach to water resources management at the catchment scale and highlights the need for modelling tools that can distinguish relative contributions from multiple nutrient sources and are consistent with the information content of the available data. Two such models are introduced and evaluated within a stochastic framework using daily flow and total phosphorus concentrations recorded in a clay catchment typical of many areas of the lowland UK. Both models disaggregate empirical annual load estimates, derived from land use data, as a function of surface/near surface runoff, generated using a simple conceptual rainfall-runoff model. Estimates of the daily load from agricultural land, together with those from baseflow and point sources, feed into an in-stream routing algorithm. The first model assumes constant concentrations in runoff via surface/near surface pathways and incorporates an additional P store in the river-bed sediments, depleted above a critical discharge, to explicitly simulate resuspension. The second model, which is simpler, simulates P concentrations as a function of surface/near surface runoff, thus emphasising the influence of non-point source loads during flow peaks and mixing of baseflow and point sources during low flows. The temporal consistency of parameter estimates and thus the suitability of each approach is assessed dynamically following a new approach based on Monte-Carlo analysis. (c) 2004 Elsevier B.V. All rights reserved.

A temporal limit to the association between flag leaf life extension by fungicides and wheat yields

Three field experiments, each repeated over two or three seasons, on winter wheat investigated a possible limit to the association between grain yield and flag leaf life, as extended by fungicide application. The experiments involved up to six cultivars and different application rates, timings and frequencies of the strobilurin azoxystrobin and the triazole epoxiconazole. In the 2000/01 and 2001/02 seasons, the relationships between the thermal time to 37 % green flag leaf area (m) and yield deviated from linearity. 'Broken stick' models were fitted to cultivar x experiment combinations within each season and the limit to the benefit to yield associated with extending flag leaf life was 700 degrees C days (S.E. = 20.7) and 725 degrees C days (S.E. = 9.33) after anthesis in 2000/01 and 2001/02, respectively. In 2002/03, the relationship between yield and in did not deviate significantly (P > 0.05) from linearity, but in this latter year the fungicide application failed to increase In past 700 degrees C days. (c) 2004 Elsevier B.V. All rights reserved.

The effects of spatial and temporal heterogeneity on the population dynamics of four animal species in a Danish landscape

Background: Variation in carrying capacity and population return rates is generally ignored in traditional studies of population dynamics. Variation is hard to study in the field because of difficulties controlling the environment in order to obtain statistical replicates, and because of the scale and expense of experimenting on populations. There may also be ethical issues. To circumvent these problems we used detailed simulations of the simultaneous behaviours of interacting animals in an accurate facsimile of a real Danish landscape. The models incorporate as much as possible of the behaviour and ecology of skylarks Alauda arvensis, voles Microtus agrestis, a ground beetle Bembidion lampros and a linyphiid spider Erigone atra. This allows us to quantify and evaluate the importance of spatial and temporal heterogeneity on the population dynamics of the four species. Results: Both spatial and temporal heterogeneity affected the relationship between population growth rate and population density in all four species. Spatial heterogeneity accounted for 23–30% of the variance in population growth rate after accounting for the effects of density, reflecting big differences in local carrying capacity associated with the landscape features important to individual species. Temporal heterogeneity accounted for 3–13% of the variance in vole, skylark and spider, but 43% in beetles. The associated temporal variation in carrying capacity would be problematic in traditional analyses of density dependence. Return rates were less than one in all species and essentially invariant in skylarks, spiders and beetles. Return rates varied over the landscape in voles, being slower where there were larger fluctuations in local population sizes. Conclusion: Our analyses estimated the traditional parameters of carrying capacities and return rates, but these are now seen as varying continuously over the landscape depending on habitat quality and the mechanisms of density dependence. The importance of our results lies in our demonstration that the effects of spatial and temporal heterogeneity must be accounted for if we are to have accurate predictive models for use in management and conservation. This is an area which until now has lacked an adequate theoretical framework and methodology.

The temporal impulse response underlying saccadic decisions

Models of perceptual decision making often assume that sensory evidence is accumulated over time in favor of the various possible decisions, until the evidence in favor of one of them outweighs the evidence for the others. Saccadic eye movements are among the most frequent perceptual decisions that the human brain performs. We used stochastic visual stimuli to identify the temporal impulse response underlying saccadic eye movement decisions. Observers performed a contrast search task, with temporal variability in the visual signals. In experiment 1, we derived the temporal filter observers used to integrate the visual information. The integration window was restricted to the first similar to 100 ms after display onset. In experiment 2, we showed that observers cannot perform the task if there is no useful information to distinguish the target from the distractor within this time epoch. We conclude that (1) observers did not integrate sensory evidence up to a criterion level, (2) observers did not integrate visual information up to the start of the saccadic dead time, and (3) variability in saccade latency does not correspond to variability in the visual integration period. Instead, our results support a temporal filter model of saccadic decision making. The temporal impulse response identified by our methods corresponds well with estimates of integration times of V1 output neurons.

Bottom-up, not top-down, modulation of imitation by human and robotic models

Visual observation of human actions provokes more motor activation than observation of robotic actions. We investigated the extent to which this visuomotor priming effect is mediated by bottom-up or top-down processing. The bottom-up hypothesis suggests that robotic movements are less effective in activating the ‘mirror system’ via pathways from visual areas via the superior temporal sulcus to parietal and premotor cortices. The top-down hypothesis postulates that beliefs about the animacy of a movement stimulus modulate mirror system activity via descending pathways from areas such as the temporal pole and prefrontal cortex. In an automatic imitation task, subjects performed a prespecified movement (e.g. hand opening) on presentation of a human or robotic hand making a compatible (opening) or incompatible (closing) movement. The speed of responding on compatible trials, compared with incompatible trials, indexed visuomotor priming. In the first experiment, robotic stimuli were constructed by adding a metal and wire ‘wrist’ to a human hand. Questionnaire data indicated that subjects believed these movements to be less animate than those of the human stimuli but the visuomotor priming effects of the human and robotic stimuli did not differ. In the second experiment, when the robotic stimuli were more angular and symmetrical than the human stimuli, human movements elicited more visuomotor priming than the robotic movements. However, the subjects’ beliefs about the animacy of the stimuli did not affect their performance. These results suggest that bottom-up processing is primarily responsible for the visuomotor priming advantage of human stimuli.