Theoretical and numerical analysis of large-scale heat transfer problems with temperature-dependent pore-fluid densities

Purpose - In many scientific and engineering fields, large-scale heat transfer problems with temperature-dependent pore-fluid densities are commonly encountered. For example, heat transfer from the mantle into the upper crust of the Earth is a typical problem of them. The main purpose of this paper is to develop and present a new combined methodology to solve large-scale heat transfer problems with temperature-dependent pore-fluid densities in the lithosphere and crust scales. Design/methodology/approach - The theoretical approach is used to determine the thickness and the related thermal boundary conditions of the continental crust on the lithospheric scale, so that some important information can be provided accurately for establishing a numerical model of the crustal scale. The numerical approach is then used to simulate the detailed structures and complicated geometries of the continental crust on the crustal scale. The main advantage in using the proposed combination method of the theoretical and numerical approaches is that if the thermal distribution in the crust is of the primary interest, the use of a reasonable numerical model on the crustal scale can result in a significant reduction in computer efforts. Findings - From the ore body formation and mineralization points of view, the present analytical and numerical solutions have demonstrated that the conductive-and-advective lithosphere with variable pore-fluid density is the most favorite lithosphere because it may result in the thinnest lithosphere so that the temperature at the near surface of the crust can be hot enough to generate the shallow ore deposits there. The upward throughflow (i.e. mantle mass flux) can have a significant effect on the thermal structure within the lithosphere. In addition, the emplacement of hot materials from the mantle may further reduce the thickness of the lithosphere. Originality/value - The present analytical solutions can be used to: validate numerical methods for solving large-scale heat transfer problems; provide correct thermal boundary conditions for numerically solving ore body formation and mineralization problems on the crustal scale; and investigate the fundamental issues related to thermal distributions within the lithosphere. The proposed finite element analysis can be effectively used to consider the geometrical and material complexities of large-scale heat transfer problems with temperature-dependent fluid densities.

Tailored conductivity in ion implanted polyetheretherketone

Ion implantation can be used to confer electrical conductivity upon conventional insulating polymers such as polyetheretherketone (PEEK). We have implanted PEEK films using three different types of ion implantation: conventional inert gas and metal ion implantation, and ion beam mixing. We have applied a number of analytical techniques to compare the chemical, structural and electrical properties of these films. The most effective means of increasing electrical conductivity appears to be via ion beam mixing of metals into the polymer, followed by metal ion implantation and finally, inert gas ion implantation. Our results suggest that in all cases, the conducting region corresponds to the implanted layer in the near surface to a depth of similar to750 Angstrom (ion beam mixed) to similar to5000 Angstrom (metal ion). This latter value is significantly higher than would be expected from a purely ballistic standpoint, and can only be attributed to thermal inter-diffusion. Our data also indicates that graphitic carbon is formed within the implant region by chain scission and subsequent cross-linking. All ion implanted samples retained their bulk mechanical properties, i.e. they remained flexible. The implant layers showed no signs of de-lamination. We believe this to be the first comparative study between different implantation techniques, and our results support the proposition that soft electronic circuitry and devices can be created by conductivity engineering with ion beams. (C) 2004 Elsevier B.V. All rights reserved.

Clay mineralogical and isotopic (K-Ar, delta O-18, delta D) constraints on the evolution of the North Anatolian Fault Zone, Turkey

This study presents the first attempt to constrain the evolution of the North Anatolian Fault Zone (NAFZ) by age dating and isotope tracing of clay minerals formed during near-surface faulting. Extensive illitic clay mineralisation occurred along the NAFZ related to hydrothermal alteration of the fault gouges and pseudotachylytes. Samples representing the pre-fault protoliths outside the fault zone do not contain authigenic illitic clay minerals indicating that hydrothermal processes were confined to the areas within the fault zone. K-Ar age data indicate that the hydrothermal system and the associated illite authigenesis initiated at similar to 57 Ma. This process is interpreted to reflect the onset of significant strike-slip or transtensional faulting immediately after the continental collision related to the closure of the Neotethys Ocean. Following the initiation of the fault movements in the latest Paleocene-Early Eocene, displacements along the NAFZ have continued, with probably intensified fault activities at similar to 26 Ma and later than similar to 8 Ma. Oxygen isotope compositions of the illitic clays from different locations along the NAFZ are similar, with narrow ranges in delta O-18 values indicating clay precipitation from fluids with similar oxygen isotope compositions and crystallisation temperatures. The delta O-18 and delta D values of the calculated fluid isotopic composition (delta O-18=5.9 parts per thousand to 11.2 parts per thousand, delta D=-59 parts per thousand to -73 parts per thousand) are consistent with metamorphic and magmatic origin of fluids mobilised during active tectonism. The interpretation of the fluid flow history of the NAFZ is in agreement with that reported previously for some well-known large-scale high-angle fault zones, which similarly developed along collisional-type orogenic belts and are commonly associated with significant mesothermal ore mineralisation. (c) 2005 Elsevier B.V. All rights reserved.

Stirring in 3-d spherical models of convection in the Earth's mantle

On a global scale basalts from mid-ocean ridges are strikingly more homogeneous than basalts from intraplate volcanism. The observed geochemical heterogeneity argues strongly for the existence of distinct reservoirs in the Earth's mantle. It is an unresolved problem of Geodynamics as to how these findings can be reconciled with large-scale convection. We review observational constraints, and investigate stirring properties of numerical models of mantle convection. Conditions in the early Earth may have supported layered convection with rapid stirring in the upper layers. Material that has been altered near the surface is transported downwards by small-scale convection. Thereby a layer of homogeneous depleted material develops above pristine mantle. As the mantle cools over Earth history, the effects leading to layering become reduced and models show the large-scale convection favoured for the Earth today. Laterally averaged, the upper mantle below the lithosphere is least affected by material that has experienced near-surface differentiation. The geochemical signature obtained during the previous episode of small-scale convection may be preserved there for the longest time. Additionally, stirring is less effective in the high viscosity layer of the central lower mantle [1, 2], supporting the survival of medium-scale heterogeneities there. These models are the first, using 3-d spherical geometry and mostly Earth-like parameters, to address the suggested change of convective style. Although the models are still far from reproducing our planet, we find that proposal might be helpful towards reconciling geochemical and geophysical constraints.