Finite element modelling of reactive mass transport problems in fluid-saturated porous media

We use the finite element method to solve reactive mass transport problems in fluid-saturated porous media. In particular, we discuss the mathematical expression of the chemical reaction terms involved in the mass transport equations for an isothermal, non-equilibrium chemical reaction. It has turned out that the Arrhenius law in chemistry is a good mathematical expression for such non-equilibrium chemical reactions especially from the computational point of view. Using the finite element method and the Arrhenius law, we investigate the distributions of PH (i.e. the concentration of H+) and the relevant reactive species in a groundwater system. Although the main focus of this study is on the contaminant transport problems in groundwater systems, the related numerical techniques and principles are equally applicable to the orebody formation problems in the geosciences. Copyright (C) 1999 John Wiley & Sons, Ltd.

A numerical study of pore-fluid, thermal and mass flow iu fluid-saturated porous rock basins

We present a numerical methodology for the study of convective pore-fluid, thermal and mass flow in fluid-saturated porous rock basins. lit particular, we investigate the occurrence and distribution pattern of temperature gradient driven convective pore-fluid flow and hydrocarbon transport in the Australian North West Shelf basin. The related numerical results have demonstrated that: (1) The finite element method combined with the progressive asymptotic approach procedure is a useful tool for dealing with temperature gradient driven pore-fluid flow and mass transport in fluid-saturated hydrothermal basins; (2) Convective pore-fluid flow generally becomes focused in more permeable layers, especially when the layers are thick enough to accommodate the appropriate convective cells; (3) Large dislocation of strata has a significant influence off the distribution patterns of convective pore;fluid flow, thermal flow and hydrocarbon transport in the North West Shelf basin; (4) As a direct consequence of the formation of convective pore-fluid cells, the hydrocarbon concentration is highly localized in the range bounded by two major faults in the basin.

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

Finite element analysis of heat transfer and mineralization in layered hydrothermal systems with upward throughflow

We use the finite element method to solve the coupled problem between convective pore-fluid flow, heat transfer and mineralization in layered hydrothermal systems with upward throughflow. In particular, we present the improved rock alteration index (IRAI) concept for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in the systems. To validate the numerical method used in the computation, analytical solutions to a benchmark problem have been derived. After the numerical method is validated, it is used to investigate the pattern of pore-fluid Aom, the distribution of temperature and the mineralization pattern of gold minerals in a layered hydrothermal system with upward throughflow. The related numerical results have demonstrated that the present concept of IRAI is useful and applicable for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in hydrothermal systems. (C) 2000 Elsevier Science S.A. All rights reserved.

Risk Assessment of Modern Landfill Structures in Finland

The purpose of this thesis was to investigate environmental permits of landfills with respect to the appropriateness of risk assessments focusing on contaminant migration, structures capable to protect the environment, waste and leachate management and existing environmental impacts of landfills. According to the requirements, a risk assessment is always required to demonstrate compliance with environmental protection requirements if the environmental permit decision deviates from the set requirements. However, there is a reason to doubt that all relevant risk factors are identified in current risk assessment practices in order to protect people end environment. In this dissertation, risk factors were recognized in 12 randomly selected landfills. Based on this analysis, a structural risk assessment method was created. The method was verified with two case examples. Several development needs were found in the risk assessments of the environmental permit decisions. The risk analysis equations used in the decisions did not adequately take into account all the determining factors like waste prospects, total risk quantification or human delineated factors. Instead of focusing on crucial factors, the landfill environmental protection capability is simply expressed via technical factors like hydraulic conductivity. In this thesis, it could be shown, that using adequate risk assessment approaches the most essential environmental impacts can be taken into account by consideration of contaminant transport mechanisms, leachate effects, and artificial landfill structures. The developed structural risk analysing (SRA) method shows, that landfills structures could be designed in a more cost-efficient way taking advantage of recycled or by-products. Additionally, the research results demonstrate that the environmental protection requirements of landfills should be updated to correspond to the capability to protect the environment instead of the current simplified requirements related to advective transport only.

Fluxo da água subterrânea em sistema de encosta-rio, município de Paulínia (SP): caracterização hidrogeológica e simulação numérica

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Mineralogical and isotopic record of diagenesis from the Opalinus Clay formation at Benken, Switzerland: Implications for the modeling of pore-water chemistry in a clay formation

Argillaceous rocks are considered to be a suitable geological barrier for the long-term containment of wastes. Their efficiency at retarding contaminant migration is assessed using reactive-transport experiments and modeling, the latter requiring a sound understanding of pore-water chemistry. The building of a pore-water model, which is mandatory for laboratory experiments mimicking in situ conditions, requires a detailed knowledge of the rock mineralogy and of minerals at equilibrium with present-day pore waters. Using a combination of petrological, mineralogical, and isotopic studies, the present study focused on the reduced Opalinus Clay formation (Fm) of the Benken borehole (30 km north of Zurich) which is intended for nuclear-waste disposal in Switzerland. A diagenetic sequence is proposed, which serves as a basis for determining the minerals stable in the formation and their textural relationships. Early cementation of dominant calcite, rare dolomite, and pyrite formed by bacterial sulfate reduction, was followed by formation of iron-rich calcite, ankerite, siderite, glauconite, (Ba, Sr) sulfates, and traces of sphalerite and galena. The distribution and abundance of siderite depends heavily on the depositional environment (and consequently on the water column). Benken sediment deposition during Aalenian times corresponds to an offshore environment with the early formation of siderite concretions at the water/sediment interface at the fluctuating boundary between the suboxic iron reduction and the sulfate reduction zones. Diagenetic minerals (carbonates except dolomite, sulfates, silicates) remained stable from their formation to the present. Based on these mineralogical and geochemical data, the mineral assemblage previously used for the geochemical model of the pore waters at Mont Terri may be applied to Benken without significant changes. These further investigations demonstrate the need for detailed mineralogical and geochemical study to refine the model of pore-water chemistry in a clay formation.

(Table 1) Grain size parameters of intertidal channel sand bars near Sylt, North-Frisian Wadden Sea

The aim of the present study is to investigate directional asymmetric properties and internal structures of the bedforms on the intertidal sand bars in comparison with the migration problems of the sand bodies developed in the channel systems of the tidal basin off the west coast of Schleswig-Holstein. The tidal channel sand bodies studied have 'V'-shaped outlines and are asymmetric in cross-section. Based on such knowledge it was hoped to understand and find possible factors for application to recent and ancient tidal depositional environments. The V-shaped intertidal channel sand bodies developed in the tidal environments between Sylt and Föhr Island are constantly migrating sand bars. The migration directions are in good agreement with the resultant vector mean directions of internal cross-stratification structures of asymmetric sedimentary bedforms. Finally, it is shown that the orientation of the apex of V-shaped sand bar as an equilibrium form alone can not indicate the migration direction, but that the orientation of the resultant vector mean of internal structures of sedimentary bedforms does indicate the migration direction. Based on the analyses of textural parameters of the migrating intertidal bar sands, it seems that sands of typical intertidal sand bars are negatively skewed and well sorted. The high rounding of quarz sand grains of these tidal channel sand bars seems to be an additional characteristical criterion for tidal depositional environments, as also indicated by Balazs and Klein (1972).

Problemas económico-sociales sobre migración y urbanización en la ciudad de Esmeraldas, República del Ecuador = Socio-economic problems of migration and urbanization in the town of Esmeraldas, Ecuador

Los documentos del Seminario fueron publicados por UNESCO en 1961 con el título: La urbanización en América Latina/Urbanization in Latin America

Influence of seaward boundary condition on contaminant transport in unconfined coastal aquifers

Contaminant transport in coastal aquifers is complicated partly due to the conditions at the seaward boundary including seawater intrusion and tidal variations of sea level. Their inclusion in modelling this system will be computationally expensive. Therefore, it will be instructive to investigate the consequence of simplifying the seaward boundary condition by neglecting the seawater density and tidal variations in numerical predictions of contaminant transport in this zone. This paper presents a comparison of numerical predictions for a simplified seaward boundary condition with experimental results for a corresponding realistic one including a saltwater interface and tidal variations. Different densities for contaminants are considered. The comparison suggests that the neglect of the seawater intrusion and tidal variations does not affect noticeably the overall migration rate of the plume before it reaches the saltwater interface. However, numerical prediction shows that a more dense contaminant travels further seaward and part of the solute mass exits under the sea if the seawater density is not included. This is not consistent with the experimental result, which shows that the contaminant travels upwards to the shoreline along the saltwater interface. Neglect of seawater density, therefore, will result in an underestimation of the exit rate of solute mass around the coastline and fictitious migration paths under the seabed. For a less dense contaminant, neglect of seawater density has little effect on numerical prediction of migration paths. (C) 2001 Elsevier Science B.V. All rights reserved.

Numerical modelling of contaminant transport in coastal aquifers

This paper employs a two-dimensional variable density flow and transport model to investigate the transport of a dense contaminant plume in an unconfined coastal aquifer. Experimental results are also presented to show the contaminant plume in a freshwater-seawater flow system. Both the numerical and experimental results suggest that the neglect of the seawater interface does not noticeably affect the horizontal migration rate of the plume before it reaches the interface. However, the contaminant will travel further seaward and part of the solute mass will exit under the sea if the higher seawater density is not included. If the seawater density is included, the contaminant will travel upwards towards the beach along the freshwater-saltwater interface as shown experimentally. Neglect of seawater density, therefore, will result in an underestimate of solute mass rate exiting around the coastline. (C) 2002 IMACS. Published by Elsevier Science B.V. All rights reserved.

Structural change in the ruhr region: problems, potentials and developments

Due to the decline of the heavy industries in the Ruhr region, the area has to reinvent itself. The orientation towards service industries proves to be a difficult task for the district and its population. This paper examines the challenges, problems and potentials of the Ruhr region against the backdrop of its economical history out of a sociological perspective. Thereby the economical situation and its outcome towards the population stand in the foreground of the paper.

Productivity gains from migration: An analysis of Cape Verdean return migrants

Are return migrants more productive than non-migrants? If so, is it a causal effect or simply self-selection? Existing literature has not reached a consensus on the role of return migration for origin countries. To answer these research questions, an empirical analysis was performed based on household data collected in Cape Verde. One of the most common identification problems in the migration literature is the presence of migrant self-selection. In order to disentangle potential selection bias, we use instrumental variable estimation using variation provided by unemployment rates in migrant destination countries, which is compared with OLS and Nearest Neighbor Matching (NNM) methods. The results using the instrumental variable approach provide evidence of labour income gains due to return migration, while OLS underestimates the coefficient of interest. This bias points towards negative self-selection of return migrants on unobserved characteristics, although the different estimates cannot be distinguished statistically. Interestingly, migration duration and occupational changes after migration do not seem to influence post-migration income. There is weak evidence that return migrants from the United States have higher income gains caused by migration than the ones who returned from Portugal.

Epidemiology of Chagas disease in non endemic countries: the role of international migration

Human infection with the protozoa Trypanosoma cruzi extends through North, Central, and South America, affecting 21 countries. Most human infections in the Western Hemisphere occur through contact with infected bloodsucking insects of the triatomine species. As T. cruzi can be detected in the blood of untreated infected individuals, decades after infection took place; the infection can be also transmitted through blood transfusion and organ transplant, which is considered the second most common mode of transmission for T. cruzi. The third mode of transmission is congenital infection. Economic hardship, political problems, or both, have spurred migration from Chagas endemic countries to developed countries. The main destination of this immigration is Australia, Canada, Spain, and the United States. In fact, human infection through blood or organ transplantation, as well as confirmed or potential cases of congenital infections has been described in Spain and in the United States. Estimates reported here indicates that in Australia in 2005-2006, 1067 of the 65,255 Latin American immigrants (16 per 1000) may be infected with T. cruzi, and in Canada, in 2001, 1218 of the 131,135 immigrants (9 per 1000) whose country of origin was identified may have been also infected. In Spain, a magnet for Latin American immigrants since the 2000, 5125 of 241,866 legal immigrants in 2003 (25 per 1000), could be infected. In the United States, 56,028 to 357,205 of the 7,20 million, legal immigrants (8 to 50 per 1000), depending on the scenario, from the period 1981-2005 may be infected with T. cruzi. On the other hand, 33,193 to 336,097 of the estimated 5,6 million undocumented immigrants in 2000 (6 to 59 per 1000) could be infected. Non endemic countries receiving immigrants from the endemic ones should develop policies to protect organ recipients from T. cruzi infection, prevent tainting the blood supply with T. cruzi, and implement secondary prevention of congenital Chagas disease.

Macrophage migration inhibitory factor is critically involved in basal and fluoxetine-stimulated adult hippocampal cell proliferation and in anxiety, depression, and memory-related behaviors.

Intensive research is devoted to unravel the neurobiological mechanisms mediating adult hippocampal neurogenesis, its regulation by antidepressants, and its behavioral consequences. Macrophage migration inhibitory factor (MIF) is a pro-inflammatory cytokine that is expressed in the CNS, where its function is unknown. Here, we show, for the first time, the relevance of MIF expression for adult hippocampal neurogenesis. We identify MIF expression in neurogenic cells (in stem cells, cells undergoing proliferation, and in newly proliferated cells undergoing maturation) in the subgranular zone of the rodent dentate gyrus. A causal function for MIF in cell proliferation was shown using genetic (MIF gene deletion) and pharmacological (treatment with the MIF antagonist Iso-1) approaches. Behaviorally, genetic deletion of MIF resulted in increased anxiety- and depression-like behaviors, as well as of impaired hippocampus-dependent memory. Together, our studies provide evidence supporting a pivotal function for MIF in both basal and antidepressant-stimulated adult hippocampal cell proliferation. Moreover, loss of MIF results in a behavioral phenotype that, to a large extent, corresponds with alterations predicted to arise from reduced hippocampal neurogenesis. These findings underscore MIF as a potentially relevant molecular target for the development of treatments linked to deficits in neurogenesis, as well as to problems related to anxiety, depression, and cognition.