Applications of hydrodynamic and water quality models to the Sau and Boadella Reservoirs

Applicaciò d'un model de hydrodinàmica i de qualitat de l'aigua als embassaments de Sau i Boadella

Contribution of the double layer to transient faradaic processes: implications for hydrodynamic modulated voltammetry of nanostructures

The electrochemistry of Pt nanostructured electrodes is investigated using hydrodynamic modulated voltammetry (HMV). Here a liquid crystal templating process is used to produce platinum-modified electrodes with a range of surface areas (roughness factor 42.4-280.8). The electroreduction of molecular oxygen at these nanostructured platinum surfaces is used to demonstrate the ability of HMV to discriminate between faradaic and nonfaradaic electrode reactions. The HMV approach shows that the reduction of molecular oxygen experiences considerable signal loss within the high pseudocapacitive region of the voltammetry. Evidence for the contribution of the double layer to transient mass transfer events is presented. In addition, a model circuit and appropriate theoretical analysis are used to illustrate the transient responses of a time variant faradaic component. This in conjunction with the experimental evidence shows that, far from being a passive component in this system, the double layer can contribute to HMV faradaic reactions under certain conditions.

Integration of a 3D Variational data assimilation scheme with a coastal area morphodynamic model of Morecambe Bay

This paper describes the implementation of a 3D variational (3D-Var) data assimilation scheme for a morphodynamic model applied to Morecambe Bay, UK. A simple decoupled hydrodynamic and sediment transport model is combined with a data assimilation scheme to investigate the ability of such methods to improve the accuracy of the predicted bathymetry. The inverse forecast error covariance matrix is modelled using a Laplacian approximation which is calibrated for the length scale parameter required. Calibration is also performed for the Soulsby-van Rijn sediment transport equations. The data used for assimilation purposes comprises waterlines derived from SAR imagery covering the entire period of the model run, and swath bathymetry data collected by a ship-borne survey for one date towards the end of the model run. A LiDAR survey of the entire bay carried out in November 2005 is used for validation purposes. The comparison of the predictive ability of the model alone with the model-forecast-assimilation system demonstrates that using data assimilation significantly improves the forecast skill. An investigation of the assimilation of the swath bathymetry as well as the waterlines demonstrates that the overall improvement is initially large, but decreases over time as the bathymetry evolves away from that observed by the survey. The result of combining the calibration runs into a pseudo-ensemble provides a higher skill score than for a single optimized model run. A brief comparison of the Optimal Interpolation assimilation method with the 3D-Var method shows that the two schemes give similar results.

A depth-integrated 2D coastal and estuarine model with conformal boundary-fitted mesh generation

Details are given of the development and application of a 2D depth-integrated, conformal boundary-fitted, curvilinear model for predicting the depth-mean velocity field and the spatial concentration distribution in estuarine and coastal waters. A numerical method for conformal mesh generation, based on a boundary integral equation formulation, has been developed. By this method a general polygonal region with curved edges can be mapped onto a regular polygonal region with the same number of horizontal and vertical straight edges and a multiply connected region can be mapped onto a regular region with the same connectivity. A stretching transformation on the conformally generated mesh has also been used to provide greater detail where it is needed close to the coast, with larger mesh sizes further offshore, thereby minimizing the computing effort whilst maximizing accuracy. The curvilinear hydrodynamic and solute model has been developed based on a robust rectilinear model. The hydrodynamic equations are approximated using the ADI finite difference scheme with a staggered grid and the solute transport equation is approximated using a modified QUICK scheme. Three numerical examples have been chosen to test the curvilinear model, with an emphasis placed on complex practical applications

Cellular automaton model for molecular traffic jams

We consider the time evolution of an exactly solvable cellular automaton with random initial conditions both in the large-scale hydrodynamic limit and on the microscopic level. This model is a version of the totally asymmetric simple exclusion process with sublattice parallel update and thus may serve as a model for studying traffic jams in systems of self-driven particles. We study the emergence of shocks from the microscopic dynamics of the model. In particular, we introduce shock measures whose time evolution we can compute explicitly, both in the thermodynamic limit and for open boundaries where a boundary-induced phase transition driven by the motion of a shock occurs. The motion of the shock, which results from the collective dynamics of the exclusion particles, is a random walk with an internal degree of freedom that determines the jump direction. This type of hopping dynamics is reminiscent of some transport phenomena in biological systems.

Multi-scale rheological model for discontinuous phenomena in materials under deformation conditions

A study of possibilities given by the developed Cellular Automata–Finite Element (CAFE) multi-scale model for prediction of the initiation and propagation of micro-shear bands and shear bands in metallic materials subjected to plastic deformation is described in the paper. Particular emphasis in defining the criterion for initiation of micro-shear and shear bands, as well as in defining the transition rules for the cellular automata, is put on accounting for the physical aspects of those phenomena occurring in two different scales in the material. The proposed approach led to the creation of the real multi-scale model of strain localization. This model predicts material behavior in various thermo-mechanical processes. Selected examples of applications of the developed model to simulations of metal forming processes, which involve strain localization, are presented in the paper. An approach based on the Smoothed Particle Hydrodynamic, which allows to overcome difficulties with remeshing in the traditional CAFE method, is presented in the paper as well. In this approach remeshing becomes possible and mesh distortion, which limits application of the CAFE method to simple deformation processes, is eliminated.

Development of the multi-scale analysis model to simulate strain localization occurring during material processing

Abstract A detailed description of possibilities given by the developed Cellular Automata—Finite Element (CAFE) multi scale model for prediction of the initiation and propagation of micro shear bands and shear bands in metallic materials subjected to plastic deformation is presented in the work. Particular emphasis in defining the criterion for initiation of micro shear and shear bands, as well as in defining the transition rules for the cellular automata, is put on accounting for the physical aspects of these phenomena occurring in two different scales in the material. The proposed approach led to the creation of the real multi scale model of strain localization phenomena. This model predicts material behavior in various thermo-mechanical processes. Selected examples of applications of the developed model to simulations of metal forming processes, which involve strain localization, are presented in the work. An approach based on the Smoothed Particle Hydrodynamic, which allows to overcome difficulties with remeshing in the traditional CAFE method, is a subject of this work as well. In the developed model remeshing becomes possible and difficulties limiting application of the CAFE method to simple deformation processes are solved. Obtained results of numerical simulaA detailed description of possibilities given by the developed Cellular Automata—Finite Element (CAFE) multi scale model for prediction of the initiation and propagation of micro shear bands and shear bands in metallic materials subjected to plastic deformation is presented in the work. Particular emphasis in defining the criterion for initiation of micro shear and shear bands, as well as in defining the transition rules for the cellular automata, is put on accounting for the physical aspects of these phenomena occurring in two different scales in the material. The proposed approach led to the creation of the real multi scale model of strain localization phenomena. This model predicts material behavior in various thermo-mechanical processes. Selected examples of applications of the developed model to simulations of metal forming processes, which involve strain localization, are presented in the work. An approach based on the Smoothed Particle Hydrodynamic, which allows to overcome difficulties with remeshing in the traditional CAFE method, is a subject of this work as well. In the developed model remeshing becomes possible and difficulties limiting application of the CAFE method to simple deformation processes are solved. Obtained results of numerical simulations are compared with the experimental results of cold rolling process to show good predicative capabilities of the developed model.tions are compared with the experimental results of cold rolling process to show good predicative capabilities of the developed model.

Dynamics of interactions involving deformable drops : hydrodynamic dimpling under attractive and repulsive electrical double layer interactions

A model developed previously to analyze force measurements between two deformable droplets in the atomic force microscope [Langmuir 2005, 21, 2912-2922] is used to model the drainage of an aqueous film between a mica plate and a deformable mercury drop for both repulsive and attractive electrical double-layer interactions between the mica and the mercury. The predictions of the model are compared with previously published data [Faraday Discuss. 2003, 123, 193-206] on the evolution of the aqueous film whose thickness has been measured with subnanometer precision. Excellent agreement is found between theoretical results and experimental data. This supports the assumptions made in the model which include no-slip boundary conditions at both interfaces. Furthermore, the successful fit attests to the utility of the model as a tool to explore details of the drainage mechanisms of nanometer-thick films in which fluid flow, surface deformations, and colloidal forces are all involved. One interesting result is that the model can predict the time at which the aqueous film collapses when attractive mica-mercury forces are present without the need to invoke capillary waves or other local instabilities of the mercury/electrolyte interface.

The central-place territorial model does not apply to space-use by juvenile brook charr Salvelinus fontinalis in lakes

1. Territoriality is widely accepted as the mechanism responsible for density-dependent mortality, emigration, and 'self-thinning' of populations of juvenile salmonine fishes in streams. Numerous studies have focused on territoriality exclusively in stream (lotic) environments and thus have fostered a stereotyped view of juvenile salmonines as sedentary and territorial. We term this behavioural paradigm the central-place territorial model (CPTM).

2. We tested predictions characterizing the CPTM for young-of-the-year (YOY) brook charr (Salvelinus fontinalis) in two Canadian lakes to determine if territoriality may also potentially limit space and population size of brook charr in lakes.

3. Our findings were not consistent with the CPTM. Fish in both lakes were not central-place forages. Maximum displacement distance did not increase with body length as predicted by the general salmonine model of Grant and Kramer (1990). Net displacement distanced increased with the proportion of time spent moving. Aggressive frequency was greatest for fish which spent large proportions of time moving and did not defend from a central-place.

4. Fish in both lakes were rarely aggressive, highly active, and often moved back over the same areas. However, lake fish which migrated to a tributary stream had no net displacement (central-place foraging) illustrating the immediate effects of current on foraging tactics and space-use.

5. The effect of hydrodynamic environment (flowing vs. still water) on fish behaviour needs to be explicitly considered in future models of salmonine behaviour.

Hydrodynamic Simulations Of Vacuum-Driven Storage Tanks

As a result of urbanization, stormwater runoff flow rates and volumes are significantly increased due to increasing impervious land cover and the decreased availability of depression storage. Storage tanks are the basic devices to efficiently control the flow rate in drainage systems during wet weather. Presented in the paper conception of vacuum-driven detention tanks allows to increase the storage capacity by usage of space above the free surface water elevation at the inlet channel. Partial vacuum storage makes possible to gain cost savings by reduction of both the horizontal area of the detention tank and necessary depth of foundations. Simulation model of vacuum-driven storage tank has been developed to estimate potential profits of its application in urban drainage system. Although SWMM5 has no direct options for vacuum tanks an existing functions (i.e. control rules) have been used to reflect its operation phases. Rainfall data used in simulations were recorded at raingage in Czestochowa during years 2010÷2012 with time interval of 10minutes. Simulation results gives overview to practical operation and maintenance cost (energy demand) of vacuum driven storage tanks depending of the ratio: vacuum-driven volume to total storage capacity. The following conclusion can be drawn from this investigations: vacuum-driven storage tanks are characterized by uncomplicated construction and control systems, thus can be applied in newly developed as well as in the existing urban drainage systems. the application of vacuum in underground detention facilities makes possible to increase of the storage capacity of existing reservoirs by usage the space above the maximum depth. Possible increase of storage capacity can achieve even a few dozen percent at relatively low investment costs. vacuum driven storage tanks can be included in existing simulation software (i.e. SWMM) using options intended for pumping stations (including control and action rules ).

Hydropower Energy Simulation Using Mike 11 Model; A Case Study In South Germany's Small Run-Of-River Hydropower Plants

Renewable energy production is a basic supplement to stabilize rapidly increasing global energy demand and skyrocketing energy price as well as to balance the fluctuation of supply from non-renewable energy sources at electrical grid hubs. The European energy traders, government and private company energy providers and other stakeholders have been, since recently, a major beneficiary, customer and clients of Hydropower simulation solutions. The relationship between rainfall-runoff model outputs and energy productions of hydropower plants has not been clearly studied. In this research, association of rainfall, catchment characteristics, river network and runoff with energy production of a particular hydropower station is examined. The essence of this study is to justify the correspondence between runoff extracted from calibrated catchment and energy production of hydropower plant located at a catchment outlet; to employ a unique technique to convert runoff to energy based on statistical and graphical trend analysis of the two, and to provide environment for energy forecast. For rainfall-runoff model setup and calibration, MIKE 11 NAM model is applied, meanwhile MIKE 11 SO model is used to track, adopt and set a control strategy at hydropower location for runoff-energy correlation. The model is tested at two selected micro run-of-river hydropower plants located in South Germany. Two consecutive calibration is compromised to test the model; one for rainfall-runoff model and other for energy simulation. Calibration results and supporting verification plots of two case studies indicated that simulated discharge and energy production is comparable with the measured discharge and energy production respectively.

Software Tools For Large Scale Interactive Hydrodynamic Modeling

Running hydrodynamic models interactively allows both visual exploration and change of model state during simulation. One of the main characteristics of an interactive model is that it should provide immediate feedback to the user, for example respond to changes in model state or view settings. For this reason, such features are usually only available for models with a relatively small number of computational cells, which are used mainly for demonstration and educational purposes. It would be useful if interactive modeling would also work for models typically used in consultancy projects involving large scale simulations. This results in a number of technical challenges related to the combination of the model itself and the visualisation tools (scalability, implementation of an appropriate API for control and access to the internal state). While model parallelisation is increasingly addressed by the environmental modeling community, little effort has been spent on developing a high-performance interactive environment. What can we learn from other high-end visualisation domains such as 3D animation, gaming, virtual globes (Autodesk 3ds Max, Second Life, Google Earth) that also focus on efficient interaction with 3D environments? In these domains high efficiency is usually achieved by the use of computer graphics algorithms such as surface simplification depending on current view, distance to objects, and efficient caching of the aggregated representation of object meshes. We investigate how these algorithms can be re-used in the context of interactive hydrodynamic modeling without significant changes to the model code and allowing model operation on both multi-core CPU personal computers and high-performance computer clusters.

Numerical analysis of refrigerant flow along non-adiabatic capillary tubes using a two-fluid model

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

On negative flows of the AKNS hierarchy and a class of deformations of a bihamiltonian structure of hydrodynamic type

A deformation parameter of a bihamiltonian structure of hydrodynamic type is shown to parametrize different extensions of the AKNS hierarchy to include negative flows. This construction establishes a purely algebraic link between, on the one hand, two realizations of the first negative flow of the AKNS model and, on the other, two-component generalizations of Camassa-Holmand Dym-type equations. The two-component generalizations of Camassa-Holm- and Dym-type equations can be obtained from the negative-order Hamiltonians constructed from the Lenard relations recursively applied on the Casimir of the first Poisson bracket of hydrodynamic type. The positive-order Hamiltonians, which follow froth the Lenard scheme applied on the Casimir of the second Poisson bracket of hydrodynamic type, are shown to coincide with the Hamiltonians of the AKNS model. The AKNS Hamiltonians give rise to charges conserved with respect to equations of motion of two-component Camassa-Holm- and two-component Dym-type equations.

A two-fluid model for refrigerant flow through adiabatic capillary tubes

This work presents a numerical model to simulate refrigerant flow through capillary tubes, commonly used as expansion devices in refrigeration systems. The flow is divided in a single-phase region, where the refrigerant is in the subcooled liquid state, and a region of two-phase flow. The capillary tube is considered straight and horizontal. The flow is taken as one-dimensional and adiabatic. Steady-state condition is also assumed and the metastable flow phenomena are neglected. The two-fluid model, considering the hydrodynamic and thermal non-equilibrium between the liquid and vapor phases, is applied to the two-phase flow region. Comparisons are made with experimental measurements of the mass flow rate and pressure distribution along two capillary tubes working with refrigerant R-134a in different operating conditions. The results indicate that the present model provides a better estimation than the commonly employed homogeneous model. Some computational results referring to the quality, void fraction, velocities, and temperatures of each phase are presented and discussed.