Modelling the fate of persistent organic pollutants (POPs) in the North Sea system

The environmental fate of selected persistent organic pollutants (POPs) in the North Sea system is modelled with a high resolution Fate and Transport Ocean Model (FANTOM) that uses hydrodynamic model output from the Hamburg Shelf Ocean Model (HAMSOM). Large amounts of POPs enter the North Sea from the surrounding highly populated, industrialised and agricultural countries. Major pathways to the North Sea are atmospheric deposition and river inputs, with additional contributions coming from bottom sediments and adjacent seas. The model domain covers the entire North Sea region, extending northward as far as the Shetland Islands, and includes adjacent basins such as the Skagerrak, Kattegat, and the westernmost part of the Baltic Sea. Model resolution (for both models) is 1.5’ latitude x 2.5’ longitude (approximately 3 km horizontal resolution) with 30 vertical levels. The POP model also has 20 sediment layers. Important model processes controlling the fate of POPs in the North Sea system are discussed. Results focus on Lindane gamma- HCH or gamma-hexachlorocyclohexane) and PCB 153.

Numerical modelling of low temperature laser produced magnesium plasma

The two-dimensional laser-plasma-interaction hydrodynamic code POLLUX has been used to simulate the ablation of a magnesium target by a 30-ns, 248-nm KrF excimer laser at low laser fluences of ≤10 J cm2. This code, originally written for much higher laser intensities, has been recently extended to include a detailed description of the equation of state in order to treat changes of phase within the target material, and also includes a Thomas Fermi description of the electrons. The simulated temporal and spatial evolution of the plasma plume in the early phase of the expansion (≤100 ns) is compared with experimental interferometric measurements of electron density. The expansion dynamics are in good agreement, although the simulated electron number density is about 2.5 times higher than the experimental values.

A Review of Numerical Modelling of Wave Energy Converter Arrays

Large-scale commercial exploitation of wave energy is certain to require the deployment of wave energy converters (WECs) in arrays, creating ‘WEC farms’. An understanding of the hydrodynamic interactions in such arrays is essential for determining optimum layouts of WECs, as well as calculating the area of ocean that the farms will require. It is equally important to consider the potential impact of wave farms on the local and distal wave climates and coastal processes; a poor understanding of the resulting environmental impact may hamper progress, as it would make planning consents more difficult to obtain. It is therefore clear that an understanding the interactions between WECs within a farm is vital for the continued development of the wave energy industry.To support WEC farm design, a range of different numerical models have been developed, with both wave phase-resolving and wave phase-averaging models now available. Phase-resolving methods are primarily based on potential flow models and include semi-analytical techniques, boundary element methods and methods involving the mild-slope equations. Phase-averaging methods are all based around spectral wave models, with supra-grid and sub-grid wave farm models available as alternative implementations.The aims, underlying principles, strengths, weaknesses and obtained results of the main numerical methods currently used for modelling wave energy converter arrays are described in this paper, using a common framework. This allows a qualitative comparative analysis of the different methods to be performed at the end of the paper. This includes consideration of the conditions under which the models may be applied, the output of the models and the relationship between array size and computational effort. Guidance for developers is also presented on the most suitable numerical method to use for given aspects of WEC farm design. For instance, certain models are more suitable for studying near-field effects, whilst others are preferable for investigating far-field effects of the WEC farms. Furthermore, the analysis presented in this paper identifies areas in which the numerical modelling of WEC arrays is relatively weak and thus highlights those in which future developments are required.

Numerical modelling of radionuclide migration in near-surface systems

Em todo o mundo são usados, hoje em dia, modelos numéricos hidrogeoquímicos para simular fenómenos naturais e fenómenos decorrentes de actividades antrópicas. Estes modelos ajudam-nos a compreender o ambiente envolvente, a sua variabilidade espacial e evolução temporal. No presente trabalho apresenta-se o desenvolvimento de modelos numéricos hidrogeoquímicos aplicados no contexto do repositório geológico profundo para resíduos nucleares de elevada actividade. A avaliação da performance de um repositório geológico profundo inclui o estudo da evolução geoquímica do repositório, bem como a análise dos cenários de mau funcionamento do repositório, e respectivas consequências ambientais. Se se escaparem acidentalmente radionuclídeos de um repositório, estes poderão atravessar as barreiras de engenharia e barreiras naturais que constituem o repositório, atingindo eventualmente, os ecosistemas superficiais. Neste caso, os sedimentos subsuperficiais constituem a última barreira natural antes dos ecosistemas superficiais. No presente trabalho foram desenvolvidos modelos numéricos que integram processos biogeoquímicos, geoquímicos, hidrodinâmicos e de transporte de solutos, para entender e quantificar a influência destes processos na mobilidade de radionuclídeos em sistemas subsuperficiais. Os resultados alcançados reflectem a robustez dos instrumentos numéricos utilizados para desenvolver simulações descritivas e predictivas de processos hidrogeoquímicos que influenciam a mobilidade de radionuclídeos. A simulação (descritiva) de uma experiência laboratorial revela que a actividade microbiana induz a diminuição do potencial redox da água subterrânea que, por sua vez, favorece a retenção de radionuclídeos sensíveis ao potencial redox, como o urânio. As simulações predictivas indicam que processos de co-precipitação com minerais de elementos maioritários, precipitação de fases puras, intercâmbio catiónico e adsorção à superfície de minerais favorecem a retenção de U, Cs, Sr e Ra na fase sólida de uma argila glaciar e uma moreia rica em calcite. A etiquetagem dos radionuclídeos nas simulações numéricas permitiu concluir que a diluição isotópica joga um papel importante no potencial impacte dos radionuclídeos nos sistemas subsuperficiais. A partir dos resultados das simulações numéricas é possivel calcular coeficientes de distribuição efectivos. Esta metodologia proporciona a simulação de ensaios de traçadores de longa duração que não seriam exequíveis à escala da vida humana. A partir destas simulações podem ser obtidos coeficientes de retardamento que são úteis no contexto da avaliação da performance de repositórios geológicos profundos.

Calibration and validation of hydrodynamic, salt and heat transport models for Ria de Aveiro lagoon (Portugal)

Ria deAveiro is a very complex shallow water coastal lagoon located on the northwest of Portugal. Important issues would be left unanswered without a good understanding of hydrodynamic and transport processes occurring in the lagoon. Calibration and validation of hydrodynamic, salt and heat transport models for Ria de Aveiro lagoon are presented. The calibration of the hydrodynamic model was performed adjusting the bottom friction coefficient, through the comparison between measured and predicted time series of sea surface elevation for 22 stations. Harmonic analysis was performed in order to evaluate the model's accuracy. To validate the hydrodynamic model measured and predicted SSE values were compared for 11 stations, as well as main flow direction velocities for 10 stations. The salt and heat transport models were calibrated comparing measured and predicted time series of salinity and water temperature for 7 stations, and the RMS of the difference between the series was determined. These models were validated comparing the model results with an independent field data set. The hydrodynamic and the salt and heat transport models for Ria de Aveiro were successfully calibrated and validated. They reproduce accurately the barotropic flows and can therefore adequately represent the salt and heat transport and the heat transfer processes occurring in Ria deAveiro.

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.

Hydrological modelling of Maranhão and Pracana Basins, a first approach

Modelling the hydrology of hydrographic basins has shown itself as a useful tool in environment management. The hydrological models can be used for multiple purposes: estimate runoff from sequences of rainfall, access stream water quality, quantify the diffuse pollution that reaches water masses such as estuaries, rivers and lakes, etc. This study has as final objective to simulate and analyse the flow, sediment transport and water quality as a function of landuse and soil type in the basins of Maranhão and Pracana. The modelling system used is SWAT, Soil Water Assessment Tool. In this first phase of the study the hydrodynamic calibration of the model was performed using measurements of average daily flows in five stations. The model compares well with the measurements; the annual average flows are similar and the majority of the measured flow peaks coincide with the model peaks.

Hydroinformatics based decision.making integrating modelling, sampling and profiling: the Ria Formosa case

A Waste Water monitoring program aiming to help decision making is presented. The program includes traditional and inboard sensor sampling, hydrodynamic and water quality modeling and a GIS based database to help the decision making of manager authorities. The focus is in the quality of waters receiving discharges from Waste Water Treatment Plants. Data was used to feed model simulations and produce hydrodynamic, effluent dispersion and ecological results. The system was then used to run different scenarios of discharge flow, concentration and location. The results enable to access the current water quality state of the lagoon and are being used as a decision making tool by the waste water managers in the evaluation phase of the treatment plant project to decide the location and the level of treatment of the discharge.

Modelling antropogenic impact on primary production of sado estuary

The assessment of human impact on complex estuarine systems is a multidisciplinary task that is highly demanding in terms of measurements and fieldwork. Nowadays the use of inexpensive and reliably modeling tools can substantially reduce the amount of measurements needed to characterize a system. These tools are also a convenient way to forecast the future evolution of the system and to study the impact of different scenarios of human influence. In this communication a modeling system composed by hydrodynamic, transport and ecological models is used to assess the current trophic state of Sado Estuary (Portugal) and to predict the future trends of the system based on different scenarios of human intervention. Special care is taken to the impact of changing riverine nutrient loads. Sado estuary is a large European estuary that has been considered until now in good trophic conditions with eutrophication appearing only in some isolated spots. Nevertheless in recent years some studies point out that the situation is changing. Sado estuary is a system with strong environmental opposing interests. It hosts a major industrial and urban center around the city of Setúbal and the upper reaches are used to intensive cultures such as rice. On the other hand the estuary possess an important ecological value since it is used by several important species of fish as a spawning and nursery area and it’s wetlands are used by many species of birds as winter shelters. Due to it’s importance the majority of Sado Estuary is considered as Natural Reserve.

Hydrodynamic and sediment fluxes through the inlets of the Ria Formosa

Tese de dout., Ciências do Mar, da Terra e do Ambiente (Ciências do Mar-Oceanografia Física), Faculdade de Ciências e Tecnologia, Univ. do Algarve, 2011

Towards a hydrodynamic operational model of the Algarve coast and St. Vincent Cape

This study attempts to implement a hydrodynamic operational model which can ultimately be used for projecting oil spill dispersal patterns and also sewage, pollution and can also be used in wave forecasting. A two layer nested model was created using MOHID Water, which is powerful ocean modelling software. The first layer (father) is used to impose the boundary conditions for the second layer (son). This was repeated for two different wind dominant regimes, Easterly and Westerly winds respectively. A qualitative comparison was done between measured tidal data and the tidal output. Sea surface temperature was also qualitatively compared with the model’s results. The results from both simulations were analysed and compared to historical literature. The comparison was done at the surface layer, 100 metre depth and at 800m depth. In the surface layer the first simulation generated an upwelling event near Cape St. Vincent and within the Algarve. The second simulation generated a non-upwelling event within which the surface was flow reversed and the warm water mass was along the Algarve coastline and evening turning clockwise around Cape St. Vincent. At the 100 metre depth for both simulations, velocity vortexes were observed near Cape St. Vincent travelling northerly and southerly at various instances. At 800metre depth a strong oceanic flow was observed moving north westerly along the continental shelf.

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.

On the development of advanced techniques for mixed-elastohydrodynamic lubrification modelling of journal and sliding bearing systems.

The present thesis is focused on the development of a thorough mathematical modelling and computational solution framework aimed at the numerical simulation of journal and sliding bearing systems operating under a wide range of lubrication regimes (mixed, elastohydrodynamic and full film lubrication regimes) and working conditions (static, quasi-static and transient conditions). The fluid flow effects have been considered in terms of the Isothermal Generalized Equation of the Mechanics of the Viscous Thin Films (Reynolds equation), along with the massconserving p-Ø Elrod-Adams cavitation model that accordingly ensures the so-called JFO complementary boundary conditions for fluid film rupture. The variation of the lubricant rheological properties due to the viscous-pressure (Barus and Roelands equations), viscous-shear-thinning (Eyring and Carreau-Yasuda equations) and density-pressure (Dowson-Higginson equation) relationships have also been taken into account in the overall modelling. Generic models have been derived for the aforementioned bearing components in order to enable their applications in general multibody dynamic systems (MDS), and by including the effects of angular misalignments, superficial geometric defects (form/waviness deviations, EHL deformations, etc.) and axial motion. The bearing exibility (conformal EHL) has been incorporated by means of FEM model reduction (or condensation) techniques. The macroscopic in fluence of the mixedlubrication phenomena have been included into the modelling by the stochastic Patir and Cheng average ow model and the Greenwood-Williamson/Greenwood-Tripp formulations for rough contacts. Furthermore, a deterministic mixed-lubrication model with inter-asperity cavitation has also been proposed for full-scale simulations in the microscopic (roughness) level. According to the extensive mathematical modelling background established, three significant contributions have been accomplished. Firstly, a general numerical solution for the Reynolds lubrication equation with the mass-conserving p - Ø cavitation model has been developed based on the hybridtype Element-Based Finite Volume Method (EbFVM). This new solution scheme allows solving lubrication problems with complex geometries to be discretized by unstructured grids. The numerical method was validated in agreement with several example cases from the literature, and further used in numerical experiments to explore its exibility in coping with irregular meshes for reducing the number of nodes required in the solution of textured sliding bearings. Secondly, novel robust partitioned techniques, namely: Fixed Point Gauss-Seidel Method (PGMF), Point Gauss-Seidel Method with Aitken Acceleration (PGMA) and Interface Quasi-Newton Method with Inverse Jacobian from Least-Squares approximation (IQN-ILS), commonly adopted for solving uid-structure interaction problems have been introduced in the context of tribological simulations, particularly for the coupled calculation of dynamic conformal EHL contacts. The performance of such partitioned methods was evaluated according to simulations of dynamically loaded connecting-rod big-end bearings of both heavy-duty and high-speed engines. Finally, the proposed deterministic mixed-lubrication modelling was applied to investigate the in fluence of the cylinder liner wear after a 100h dynamometer engine test on the hydrodynamic pressure generation and friction of Twin-Land Oil Control Rings.

Mathematical modelling of the lower urinary tract

The lower urinary tract is one of the most complex biological systems of the human body as it involved hydrodynamic properties of urine and muscle. Moreover, its complexity is increased to be managed by voluntary and involuntary neural systems. In this paper, a mathematical model of the lower urinary tract it is proposed as a preliminary study to better understand its functioning. Furthermore, another goal of that mathematical model proposal is to provide a basis for developing artificial control systems. Lower urinary tract is comprised of two interacting systems: the mechanical system and the neural regulator. The latter has the function of controlling the mechanical system to perform the voiding process. The results of the tests reproduce experimental data with high degree of accuracy. Also, these results indicate that simulations not only with healthy patients but also of patients with dysfunctions with neurological etiology present urodynamic curves very similar to those obtained in clinical studies.