Extension And Testing Of A 2D Hydrodynamic Model For Direct Rainfall Runoff Simulation

While the simulation of flood risks originating from the overtopping of river banks is well covered within continuously evaluated programs to improve flood protection measures, flash flooding is not. Flash floods are triggered by short, local thunderstorm cells with high precipitation intensities. Small catchments have short response times and flow paths and convective thunder cells may result in potential flooding of endangered settlements. Assessing local flooding and pathways of flood requires a detailed hydraulic simulation of the surface runoff. Hydrological models usually do not incorporate surface runoff at this detailedness but rather empirical equations are applied for runoff detention. In return 2D hydrodynamic models usually do not allow distributed rainfall as input nor are any types of soil/surface interaction implemented as in hydrological models. Considering several cases of local flash flooding during the last years the issue emerged for practical reasons but as well as research topics to closing the model gap between distributed rainfall and distributed runoff formation. Therefore, a 2D hydrodynamic model, depth-averaged flow equations using the finite volume discretization, was extended to accept direct rainfall enabling to simulate the associated runoff formation. The model itself is used as numerical engine, rainfall is introduced via the modification of waterlevels at fixed time intervals. The paper not only deals with the general application of the software, but intends to test the numerical stability and reliability of simulation results. The performed tests are made using different artificial as well as measured rainfall series as input. Key parameters of the simulation such as losses, roughness or time intervals for water level manipulations are tested regarding their impact on the stability.

2D-drop model applied to the calculation of step formation energies on a (111) substrate

A model is presented for obtaining the step formation energy for metallic islands on (1 1 1) surfaces from Monte Carlo simulations. This model is applied to homo (Cu/Cu(1 1 1), Ag/Ag(1 1 1)) and heteroepitaxy (Ag/Pt(1 1 1)) systems. The embedded atom method is used to represent the interaction between the particles of the system, but any other type of potential could be used as well. The formulation can also be employed to consider the case of other single crystal surfaces, since the higher barriers for atom motion on other surfaces are not a hindrance for the simulation scheme proposed.

Coupling anchovy eggs and hydrodynamic model as an assessement tool for estuarine ecohydrological management

Nitrate and urban waste water directives have raised the need for a better understanding of coastal systems in European Union. The incorrect application of these directives can lead to important ecological or social penalties. In the paper this problem is addressed to Ria Formosa Coastal Lagoon. Ria Formosa hosts a Natural Park, important ports of the southern Portuguese coast and significant bivalve aquaculture activity. Four major urban waste water treatment plants discharging in the lagoon are considered in this study. Its treatment level must be selected, based on detailed information from a monitoring program and on a good knowledge of the processes determining the fate of the material discharged in the lagoon. In this paper the results of a monitoring program and simulations using a coupled hydrodynamic and water quality / ecological model, MOHID, are used to characterise the system and to understand the processes in Ria Formosa. It is shown that the water residence time in most of the lagoon is quite low, of the order of days, but it can be larger in the upper parts of the channels where land generated water is discharged. The main supply of nutrients to the lagoon comes from the open sea rather than from the urban discharges. For this reason the characteristics and behaviour of the general lagoon contrasts with the behaviour of the upper reaches of the channels where the influence of the waste water treatment plants are high. In this system the bottom mineralization was found to be an important mechanism, and the inclusion of that process in the model was essential to obtain good results.

A three-dimensional hydrodynamic model with generic vertical coordinate

A three-dimensional primitive equation model and its application to a tidal estuary is described. The model solves the primitive equations for incompressible fluids with Boussinesq and hydrostatic approximations. The discretization is based on the finite volume method and allows a general vertical coordinate. The computational code is implemented in such a way that different vertical coordinates can be used in different parts of the domain. The model was designed to be able to simulate the flow both in the open ocean and in coastal and estuarine zones and can be coupled in a simple way to ecological models. The model was implemented successfully in several estuarine and coastal areas. Results are show for the Sado estuary in Portugal to illustrate model accuracy and potential. Quantitative validation is based on field data (water levels and velocities) while qualitative verification is based on the analysis of secondary flows.

Hydrodynamic model intercomparison for the patos lagoon (Brazil)

An intercomparison study is carried out between two models with different formulations and spatial discretizations in order to overcome the limitations posed by the standard calibration and validation procedures and improve confidence in the hydrodynamic results for the Patos Lagoon. Numerical simulations were carried out applying the TELEMAC and MOHID models, based on the same boundary conditions and identical calibration coefficients so differences in calculated flow conditions result from the formulations and parameterizations of each model. Results from both models are compared with measurements from three stations inside the lagoon. Preliminary results indicate that both models compare well with the measurements and with each other. These results increase the confidence on hydrodynamic results for the Patos Lagoon and provide the first step towards water quality studies for the area.

Temperature Dynamics Investigation At Small And Shallow Lakes Using Hydrodynamic Model

A three-dimensional time-dependent hydrodynamic and heat transport model of Lake Binaba, a shallow and small dam reservoir in Ghana, emphasizing the simulation of dynamics and thermal structure has been developed. Most numerical studies of temperature dynamics in reservoirs are based on one- or two-dimensional models. These models are not applicable for reservoirs characterized with complex flow pattern and unsteady heat exchange between the atmosphere and water surface. Continuity, momentum and temperature transport equations have been solved. Proper assignment of boundary conditions, especially surface heat fluxes, has been found crucial in simulating the lake’s hydrothermal dynamics. This model is based on the Reynolds Average Navier-Stokes equations, using a Boussinesq approach, with a standard k − ε turbulence closure to solve the flow field. The thermal model includes a heat source term, which takes into account the short wave radiation and also heat convection at the free surface, which is function of air temperatures, wind velocity and stability conditions of atmospheric boundary layer over the water surface. The governing equations of the model have been solved by OpenFOAM; an open source, freely available CFD toolbox. As its core, OpenFOAM has a set of efficient C++ modules that are used to build solvers. It uses collocated, polyhedral numerics that can be applied on unstructured meshes and can be easily extended to run in parallel. A new solver has been developed to solve the hydrothermal model of lake. The simulated temperature was compared against a 15 days field data set. Simulated and measured temperature profiles in the probe locations show reasonable agreement. The model might be able to compute total heat storage of water bodies to estimate evaporation from water surface.

Dynamical collapse in a degenerate binary fermion mixture using a hydrodynamic model

We use a time-dependent dynamical hydrodynamic model to study a collapse in a degenerate fermion-fermion mixture ( DFFM) of different atoms. Due to a strong Pauli-blocking repulsion among identical spin-polarized fermions at short distances, there cannot be a collapse for repulsive interspecies fermion fermion interaction. However, there can be a collapse for a sufficiently attractive interspecies fermion-fermion interaction in a DFFM of different atoms. Using a variational analysis and numerical solution of the hydrodynamic model, we study different aspects of collapse in such a DFFM initiated by a jump in the interspecies fermion-fermion interaction ( scattering length) to a large negative ( attractive) value using a Feshbach resonance. Suggestion for experiments of collapse in a DFFM of distinct atoms is made.

Dynamic critical behavior of percolation observables in the 2d Ising model

We present preliminary results of our numerical study of the critical dynamics of percolation observables for the two-dimensional Ising model. We consider the (Monte-Carlo) short-time evolution of the system obtained with a local heat-bath method and with the global Swendsen-Wang algorithm. In both cases, we find qualitatively different dynamic behaviors for the magnetization and Omega, the order parameter of the percolation transition. This may have implications for the recent attempts to describe the dynamics of the QCD phase transition using cluster observables.

Dynamic behavior of cluster observables for the 2d Ising model

We present results of our numerical study of the critical dynamics of percolation observables for the two-dimensional Ising model. We consider the (Monte Carlo) short-time evolution of the system with small initial magnetization and heat-bath dynamics. We find qualitatively different dynamic behaviors for the magnetization M and for Ω, the so-called strength of the percolating cluster, which is the order parameter of the percolation transition. More precisely, we obtain a (leading) exponential form for Ω as a function of the Monte Carlo time t, to be compared with the power-law increase encountered for M at short times. Our results suggest that, although the descriptions in terms of magnetic or percolation order parameters may be equivalent in the equilibrium regime, greater care must be taken to interpret percolation observables at short times.

Topological lattice actions for the 2d XY model

We consider the 2d XY Model with topological lattice actions, which are invariant against small deformations of the field configuration. These actions constrain the angle between neighbouring spins by an upper bound, or they explicitly suppress vortices (and anti-vortices). Although topological actions do not have a classical limit, they still lead to the universal behaviour of the Berezinskii-Kosterlitz-Thouless (BKT) phase transition — at least up to moderate vortex suppression. In the massive phase, the analytically known Step Scaling Function (SSF) is reproduced in numerical simulations. However, deviations from the expected universal behaviour of the lattice artifacts are observed. In the massless phase, the BKT value of the critical exponent ηc is confirmed. Hence, even though for some topological actions vortices cost zero energy, they still drive the standard BKT transition. In addition we identify a vortex-free transition point, which deviates from the BKT behaviour.

Optimization of a label-free biosensor vertically characterized based on a periodic lattice of high aspect ratio SU-8 nano-pillars with a simplified 2D theoretical model

We have recently demonstrated a biosensor based on a lattice of SU8 pillars on a 1 μm SiO2/Si wafer by measuring vertically reflectivity as a function of wavelength. The biodetection has been proven with the combination of Bovine Serum Albumin (BSA) protein and its antibody (antiBSA). A BSA layer is attached to the pillars; the biorecognition of antiBSA involves a shift in the reflectivity curve, related with the concentration of antiBSA. A detection limit in the order of 2 ng/ml is achieved for a rhombic lattice of pillars with a lattice parameter (a) of 800 nm, a height (h) of 420 nm and a diameter(d) of 200 nm. These results correlate with calculations using 3D-finite difference time domain method. A 2D simplified model is proposed, consisting of a multilayer model where the pillars are turned into a 420 nm layer with an effective refractive index obtained by using Beam Propagation Method (BPM) algorithm. Results provided by this model are in good correlation with experimental data, reaching a reduction in time from one day to 15 minutes, giving a fast but accurate tool to optimize the design and maximizing sensitivity, and allows analyzing the influence of different variables (diameter, height and lattice parameter). Sensitivity is obtained for a variety of configurations, reaching a limit of detection under 1 ng/ml. Optimum design is not only chosen because of its sensitivity but also its feasibility, both from fabrication (limited by aspect ratio and proximity of the pillars) and fluidic point of view. (© 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

Parameters affecting the performance of wetting and drying in a two-dimensional finite element long wave hydrodynamic model

The marsh porosity method, a type of thin slot wetting and drying algorithm in a two-dimensional finite element long wave hydrodynamic model, is discussed and analyzed to assess model performance. Tests, including comparisons to simple examples and theoretical calculations, examine the effects of varying the marsh porosity parameters. The findings demonstrate that the wetting and drying concept of marsh porosity, often used in finite element hydrodynamic modeling, can behave in a more complex manner than initially expected.