Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

Computer simulations of coupled problems in geological and geochemical systems

The finite element method is used to simulate coupled problems, which describe the related physical and chemical processes of ore body formation and mineralization, in geological and geochemical systems. The main purpose of this paper is to illustrate some simulation results for different types of modelling problems in pore-fluid saturated rock masses. The aims of the simulation results presented in this paper are: (1) getting a better understanding of the processes and mechanisms of ore body formation and mineralization in the upper crust of the Earth; (2) demonstrating the usefulness and applicability of the finite element method in dealing with a wide range of coupled problems in geological and geochemical systems; (3) qualitatively establishing a set of showcase problems, against which any numerical method and computer package can be reasonably validated. (C) 2002 Published by Elsevier Science B.V.

Heroin-assisted treatment as a response to the public health problem of opiate dependence

Injection drug use (involving the injection of illicit opiates) poses serious public health problems in many countries. Research has indicated that injection drug users are at higher risk for morbidity in the form of HIV/AIDS and Hepatitis B and C, and drug-related mortality, as well as increased criminal activity. Methadone maintenance treatment is the most prominent form of pharmacotherapy treatment for illicit opiate dependence in several countries, and its application varies internationally with respect to treatment regulations and delivery modes. In order to effectively treat those patients who have previously been resistant to methadone maintenance treatment, several countries have been studying and/or considering heroin-assisted treatment as a complementary form of opiate pharmacotherapy treatment. This paper provides an overview of the prevalence of injection drug use and the opiate dependence problem internationally, the current opiate dependence treatment landscape in several countries, and the status of ongoing or planned heroin-assisted treatment trials in Australia, Canada and certain European countries.

Analysis of steady-state heat transfer through mid-crustal vertical cracks with upward throughflow in hydrothermal systems

We conduct a theoretical analysis of steady-state heat transfer problems through mid-crustal vertical cracks with upward throughflow in hydrothermal systems. In particular, we derive analytical solutions for both the far field and near field of the system. In order to investigate the contribution of the forced advection to the total temperature of the system, two concepts, namely the critical Peclet number and the critical permeability of the system, have been presented and discussed in this paper. The analytical solution for the far field of the system indicates that if the pore-fluid pressure gradient in the crust is lithostatic, the critical permeability of the system can be used to determine whether or not the contribution of the forced advection to the total temperature of the system is negligible. Otherwise, the critical Peclet number should be used. For a crust of moderate thickness, the critical permeability is of the order of magnitude of 10(-20) m(2), under which heat conduction is the overwhelming mechanism to transfer heat energy, even though the pore-fluid pressure gradient in the crust is lithostatic. Furthermore, the lower bound analytical solution for the near field of the system demonstrates that the permeable vertical cracks in the middle crust can efficiently transfer heat energy from the lower crust to the upper crust of the Earth. Copyright (C) 2002 John Wiley Sons, Ltd.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.

Theoretical and numerical analyses of convective instability in porous media with temperature-dependent viscosity

Exact analytical solutions of the critical Rayleigh numbers have been obtained for a hydrothermal system consisting of a horizontal porous layer with temperature-dependent viscosity. The boundary conditions considered are constant temperature and zero vertical Darcy velocity at both the top and bottom of the layer. Not only can the derived analytical solutions be readily used to examine the effect of the temperature-dependent viscosity on the temperature-gradient driven convective flow, but also they can be used to validate the numerical methods such as the finite-element method and finite-difference method for dealing with the same kind of problem. The related analytical and numerical results demonstrated that the temperature-dependent viscosity destabilizes the temperature-gradient driven convective flow and therefore, may affect the ore body formation and mineralization in the upper crust of the Earth. Copyright (C) 2003 John Wiley Sons, Ltd.

Electrooxidation of glyphosate herbicide at different DSA (R) compositions: pH, concentration and supporting electrolyte effect

This paper has investigated the electrochemical oxidation of glyphosate herbicide (GH) on RuO(2) and IrO(2) dimensionally stable anode (DSA (R)) electrodes. Electrolysis was achieved under galvanostatic control as a function of pH, GH concentration, supporting electrolyte, and current density. The influence of the oxide composition on GH degradation seems to be significant in the absence of chloride; Ti/Ir(0.30)Sn(0.70)O(2) is the best electrode material to oxidize GH. GH oxidation is favored at low pH values. The use of chloride medium increases the oxidizing power and the influence of the oxide composition is meaningless. At 30 mA cm(-2) and 4 h of electrolysis, complete GH removal from the electrolyzed solution has been obtained. In chloride medium, application of 50 mA cm(-2) leads to virtually total mineralization ( release of phosphate ions = 91%) for all the evaluated oxide materials. (C) 2008 Elsevier Ltd. All rights reserved.

Calcium Carbonate Particle Growth Depending on Coupling among Adjacent Layers in Hybrid LB/LbL Films

There are practical and academic situations that justify the study of calcium carbonate crystallization and especially of systems that are associated with organic matrices and a confined medium. Despite the fact that many different matrices have been studied, the use of well-behaved, thin organic films may provide new knowledge about this system. In this work, we have studied the growth of calcium carbonate particles on well-defined organic matrices that were formed by layer-by-layer (LbL) polyelectrolyte films deposited on phospholipid Langmuir-Blodgett films (LB). We were able to change the surface electrical charge density of the LB films by changing the proportions of a negatively charged lipid, the sodium salt of dimyristoyl-sn-glycero-phosphatidyl acid (DMPA), and a zwitterionic lipid. dimyristoyl-sn-glycero-phosphatidylethanolamine (DMPE). This affects the subsequent polyelectrolyte LbL film deposition, which also changes the the nature of the bonding (electrostatic interaction or hydrogen bonding). This approach allowed for the formation of calcium carbonate particles of different final shapes, roughnesses, and sizes. The masses of deposited lipids, polyelectrolytes, and calcium cabonate were quantified by the quartz crystal microbalance technique. The structures of obtained particles were analyzed by scanning electron microscopy.

Electrochemical oxidation of synthetic tannery wastewater in chloride-free aqueous media

The electrochemical treatment of a synthetic tannery wastewater, prepared with several compounds used by finishing tanneries, was studied in chloride-free media. Boron-doped diamond (Si/BDD), antimony-doped tin dioxide (Ti/SnO(2)-Sb), and iridium-antimony-doped tin dioxide (Ti/SnO(2)-Sb-Ir)were evaluated as anode. The influence of pH and current density on the treatment was assessed by means of the parameters used to measure the level of organic contaminants in the wastewater; i.e., total phenols, chemical oxygen demand (COD), total organic carbon (TOC), and absorbance. Results showed that faster decrease in these parameters occurred when the Si/BDD anode was used. Good results were obtained with the Ti/SnO(2)-Sb anode, but its complete deactivation was reached after 4h of electrolysis at 25 mA cm(-2), indicating that the service life of this electrode is short. The Ti/SnO(2)-Sb-Ir anode is chemically and electrochemically more stable than the Ti/SnO(2)-Sb anode, but it is not suitable for the electrochemical treatment under the studied conditions. No significant changes were observed for electrolyses performed at different pH conditions with Si/BDD, and this electrode led to almost complete mineralization after 4 h of electrolysis at 100mAcm(-2). The increase in current density resulted in faster wastewater oxidation, with lower current efficiency and higher energy consumption. Si/BBD proved to be the best electrodic material for the direct electrooxidation of tannery wastewaters. (C) 2010 Elsevier B.V. All rights reserved.

Electrochemical oxidation of acid black 210 dye on the boron-doped diamond electrode in the presence of phosphate ions: Effect of current density, pH, and chloride ions

The electrochemical oxidation of acid black 210 dye (AB-210) on the boron-doped diamond (BDD) was investigated under different pH conditions. The best performance for the AB-210 oxidation occurred in alkaline phosphate solution. This is probably due to oxidizing agents such as phosphate radicals and peroxodiphosphate ions, which can be electrochemically produced with good yields on the BDD anode, mainly in alkaline solution. Under this condition, the COD (chemical oxygen demand) removal was higher than that obtained from the model proposed by Comninellis. Electrolyses performed in phosphate buffer and in the presence of chloride ions resulted in faster COD and color removals in acid and neutral solutions, but in alkaline phosphate solution, a better performance in terms of TOC removal was obtained in the absence of chloride. Moreover, organochloride compounds were detected in all electrolyses performed in the presence of chloride. The AB-210 electrooxidation on BDD using phosphate as supporting electrolyte proved to be interesting since oxidizing species generated from phosphate ions were able to completely degrade the dye without producing organochloride compounds. (C) 2009 Elsevier Ltd. All rights reserved.

Electrochemical degradation of glyphosate formulations at DSA(A (R)) anodes in chloride medium: an AOX formation study

The electrochemical degradation of different glyphosate herbicide formulations on RuO(2) and IrO(2) DSA(A (R)) electrodes is investigated. Parameters that could influence the formation of organochloride compounds during electrolysis are studied. The effects of chloride concentration, electrodic composition, current density, and electrolysis time are reported. The influence of the oxide composition on herbicide degradation seems to be almost insignificant; however, there is a straight relationship between anode composition and organic halides formation. Commercial herbicide formulations have lower degradation rates and lead to the formation of a larger quantities of organochloride compounds. In high chloride concentrations, there is a significant increase in organic mineralization, and the relationship between chloride concentration and organic halides formation is direct. Only in low chloride medium investigated the organochloride concentration obtained was below the limit values allowed in Brazil. The determination of organic halides absorbable (AOX) during electrolysis increases significantly with the applied current. Even during long-term electrolysis, a large amount of organochloride compounds is formed.

MINERALOGY AND GENESIS OF SMECTITES IN AN ALKALINE-SALINE ENVIRONMENT OF PANTANAL WETLAND, BRAZIL

Smectite formation in alkaline-saline environments has been attributed to direct precipitation from solution and/or transformation from precursor minerals, but these mechanisms are not universally agreed upon in the literature. The objective of this work was to investigate the mineralogy of smectites in the soils surrounding a representative alkaline-saline lake of Nhecolandia, a sub-region of the Pantanal wetland, Brazil, and then to identify the mechanisms of their formation. Soils were sampled along a toposequence and analyzed by X-ray diffraction, transmission electron microscopy-energy dispersive X-ray analysis, and inductively coupled plasma-mass spectrometry. Water was collected along a transect involving the studied toposequence and equilibrium diagrams were calculated using the databases PHREEQC and AQUA. The fine-clay fraction is dominated by smectite, mica, and kaolinite. Smectites are concentrated at two places in the toposequence: an upper zone, which includes the soil horizons rarely reached by the lake-level variation; and a lower zone, which includes the surface horizon within the area of seasonal lake-level variation. Within the upper zone, the smectite is dioctahedral, rich in Al and Fe, and is classified as ferribeidellite. This phase is interstratified with mica and vermiculite and has an Fe content similar to that of the mica identified. These characteristics suggest that the ferribeidellite originates from transformation of micas and that vermiculite is an intermediate phase in this transformation. Within the lower zone, smectites are dominantly trioctahedral, Mg-rich, and are saponitic and stevensitic minerals. In addition, samples enriched in these minerals have much smaller rare-earth element (REE) contents than other soil samples. The water chemistry shows a geochemical control of Mg and saturation with respect to Mg-smectites in the more saline waters. The REE contents, water chemistry, and the presence of Mg-smectite where maximum evaporation is expected, suggest that saponitic and stevensitic minerals originate by chemical precipitation from the water column of the alkaline-saline lake.

Numerical modelling of chemical effects of magma solidification problems in porous rocks

The solidification of intruded magma in porous rocks can result in the following two consequences: (1) the heat release due to the solidification of the interface between the rock and intruded magma and (2) the mass release of the volatile fluids in the region where the intruded magma is solidified into the rock. Traditionally, the intruded magma solidification problem is treated as a moving interface (i.e. the solidification interface between the rock and intruded magma) problem to consider these consequences in conventional numerical methods. This paper presents an alternative new approach to simulate thermal and chemical consequences/effects of magma intrusion in geological systems, which are composed of porous rocks. In the proposed new approach and algorithm, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with the proposed mass source and physically equivalent heat source. The major advantage in using the proposed equivalent algorithm is that a fixed mesh of finite elements with a variable integration time-step can be employed to simulate the consequences and effects of the intruded magma solidification using the conventional finite element method. The correctness and usefulness of the proposed equivalent algorithm have been demonstrated by a benchmark magma solidification problem. Copyright (c) 2005 John Wiley & Sons, Ltd.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.