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.

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.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

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.

Determination of mechanical properties of PECVD silicon nitride thin films for tunable MEMS Fabry-Pérot optical filters

This paper reports an investigation on techniques for determining elastic modulus and intrinsic stress gradient in plasma-enhanced chemical vapor deposition (PECVD) silicon nitride thin films. The elastic property of the silicon nitride thin films was determined using the nanoindentation method on silicon nitride/silicon bilayer systems. A simple empirical formula was developed to deconvolute the film elastic modulus. The intrinsic stress gradient in the films was determined by using micrometric cantilever beams, cross-membrane structures and mechanical simulation. The deflections of the silicon nitride thin film cantilever beams and cross-membranes caused by in-thickness stress gradients were measured using optical interference microscopy. Finite-element beam models were built to compute the deflection induced by the stress gradient. Matching the deflection computed under a given gradient with that measured experimentally on fabricated samples allows the stress gradient of the PECVD silicon nitride thin films introduced from the fabrication process to be evaluated.

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.

Three-dimensional electrical impedance, tomography: A topology optimization approach

Electrical impedance tomography is a technique to estimate the impedance distribution within a domain, based on measurements on its boundary. In other words, given the mathematical model of the domain, its geometry and boundary conditions, a nonlinear inverse problem of estimating the electric impedance distribution can be solved. Several impedance estimation algorithms have been proposed to solve this problem. In this paper, we present a three-dimensional algorithm, based on the topology optimization method, as an alternative. A sequence of linear programming problems, allowing for constraints, is solved utilizing this method. In each iteration, the finite element method provides the electric potential field within the model of the domain. An electrode model is also proposed (thus, increasing the accuracy of the finite element results). The algorithm is tested using numerically simulated data and also experimental data, and absolute resistivity values are obtained. These results, corresponding to phantoms with two different conductive materials, exhibit relatively well-defined boundaries between them, and show that this is a practical and potentially useful technique to be applied to monitor lung aeration, including the possibility of imaging a pneumothorax.

A numerical study of effect of grain boundaries on elastic and plastic properties in nanocomposite materials

Nanocomposite materials have received considerable attention in recent years due to their novel properties. Grain boundaries are considered to play an important role in nanostructured materials. This work focuses on the finite element analysis of the effect of grain boundaries on the overall mechanical properties of aluminium/alumina composites. A grain boundary is incorporated into the commonly used unit cell model to investigate its effect on material properties. By combining the unit cell model with an indentation model, coupled with experimental indentation measurements, the ''effective'' plastic property of the grain boundary is estimated. In addition, the strengthening mechanism is also discussed based on the Estrin-Mecking model.

The subcycled Newmark algorithm

The popular Newmark algorithm, used for implicit direct integration of structural dynamics, is extended by means of a nodal partition to permit use of different timesteps in different regions of a structural model. The algorithm developed has as a special case an explicit-explicit subcycling algorithm previously reported by Belytschko, Yen and Mullen. That algorithm has been shown, in the absence of damping or other energy dissipation, to exhibit instability over narrow timestep ranges that become narrower as the number of degrees of freedom increases, making them unlikely to be encountered in practice. The present algorithm avoids such instabilities in the case of a one to two timestep ratio (two subcycles), achieving unconditional stability in an exponential sense for a linear problem. However, with three or more subcycles, the trapezoidal rule exhibits stability that becomes conditional, falling towards that of the central difference method as the number of subcycles increases. Instabilities over narrow timestep ranges, that become narrower as the model size increases, also appear with three or more subcycles. However by moving the partition between timesteps one row of elements into the region suitable for integration with the larger timestep these the unstable timestep ranges become extremely narrow, even in simple systems with a few degrees of freedom. As well, accuracy is improved. Use of a version of the Newmark algorithm that dissipates high frequencies minimises or eliminates these narrow bands of instability. Viscous damping is also shown to remove these instabilities, at the expense of having more effect on the low frequency response.

Mechanical behavior of prosthesis in Toucan beak (Ramphastos toco)

The purpose of this study is to characterize the structure of the beak of Toco Toucan (Ramphastos toco) and to investigate means for arresting fractures in the rhinotheca using acrylic resin. The structure of the rhamphastid bill has been described as a sandwich structured composite having a thin exterior comprised of keratin and a thick foam core constructed of mineralized collagenous rods (trabeculae). The keratinous rhamphotheca consists of superposed polygonal scales (approximately 50 pm in diameter and 1 mu m in thickness). In order to simulate the orientation of loading to which the beak is subjected during exertion of bite force, for example, we conducted flexure tests on the dorso-ventral axis of the maxilla. The initially intact (without induced fracture) beak fractured in the central portion when subjected to a force of 270 N, at a displacement of 23 mm. The location of this fracture served as a reference for the fractures induced in other beaks tested. The second beak was fractured and repaired by applying resin on both lateral surfaces. The repaired maxilla sustained a force of 70 N with 6.5 mm deflection. The third maxilla was repaired similarly except that it was conditioned in acid for 60s prior to fixation with resin. It resisted a force of up to 63 N at 6 mm of deflection. The experimental results were compared with finite element calculations for unfractured beak in bending configuration. The repaired specimens were found to have strength equal to only one third of the intact beak. Finite element simulations allow visualization of how the beak system (sandwich shell and cellular core) sustains high flexural strength. (C) 2010 Elsevier B.V. All rights reserved.

Low-fluence CO(2) laser irradiation decreases enamel solubility

This study investigated whether subablative-pulsed CO(2) laser (10.6 mu m) irradiation, using fluences lower than 1 J/cm(2), was capable of reducing enamel acid solubility. Fifty-one samples of bovine dental enamel were divided into three groups: control group, which was not irradiated (CG); group laser A (LA) irradiated with 0.3 J/cm ; and group laser B (LB) irradiated with 0.7 J/cm(2). After irradiation, the samples were subjected to demineralization in an acetate buffer solution and were then analyzed by SEM. A finite-element model was used to calculate the temperature increase. The calcium and phosphorous content in the demineralization solution were measured with an ICP-OES. ANOVA and the t-test pairwise comparison (p < 0.016) revealed that LB showed significantly lower mean Ca and P content values in the demineralization solution than other groups. A reduction in the enamel solubility can be obtained with pulsed CO(2) laser irradiation (0.7 J/cm(2), 135 mJ/pulse, 74 Hz, 100 mu s) without any surface photomodification and a less than 2 degrees C temperature increase at a 3-mm depth from the surface.

Understanding Contradictory Data in Contraction Stress Tests

The literature shows contradictory results regarding the role of composite shrinkage and elastic modulus as determinants of polymerization stress. The present study aimed at a better understanding of the test mechanics that could explain such divergences among studies. The hypothesis was that the effects of composite shrinkage and elastic modulus on stress depend upon the compliance of the testing system. A commonly used test apparatus was simulated by finite element analysis, with different compliance levels defined by the bonding substrate (steel, glass, composite, or acrylic). Composites with moduli between 1 and 12 GPa and shrinkage values between 0.5% and 6% were modeled. Shrinkage was simulated by thermal analogy. The hypothesis was confirmed. When shrinkage and modulus increased simultaneously, stress increased regardless of the substrate. However, if shrinkage and modulus were inversely related, their magnitudes and interaction with rod material determined the stress response.

Measuring Bond Strength Between Fiber Post and Root Dentin: A Comparison of Different Tests

Purpose The aim of this study was to evaluate the ability of bond strength tests to accurately measure the bond strength of fiber posts luted into root canals Materials and Methods The test methods studied were hourglass microtensile (HM), push-out (PS), modified push out (MP) and pull out (PL) The evaluated parameters were bond strength values, reliability (using Weibull analysis), failure mode (using confocal microscopy), and stress distribution (using finite element analysis) Forty human intact single rooted and endodontically treated teeth were divided into four groups Each group was assigned one of the test methods The samples in the HM and PS groups were 1 0 +/- 0 1 mm thick, the HM samples were hourglass shaped and the PS samples were disk shaped For the PL and MP groups, each 1 mm dentin slice was luted with a fiber post piece Three dimensional models of each group were made and stress was analyzed based on Von Mises criteria Results PL provided the highest values of bond strength followed by MP both of which also had greater amounts of adhesive failures PS showed the highest frequency of cohesive failures MP showed a more homogeneous stress distribution and a higher Weibull modulus Conclusion The specimen design directly influences the biomechanical behavior of bond strength tests

Can Fiber Posts Increase Root Stresses and Reduce Fracture?

The clinical success of fiber posts has been attributed to their lower elastic modulus. The tested hypothesis was that fiber posts could lead to lower risk of post debonding and lower risk of root fracture, despite an increase in root stresses. Stress analyses were carried out with a 3D finite element model of a premolar restored with a metallic or a fiber post. Bonded and non-bonded post/cement interface conditions were simulated. We calculated risk-of-fracture indices by determining the highest principal stress values divided by the tensile strength. Shear stresses along the post/cement interface were analyzed for the bonded models. Compared with the premolar restored with a metallic post, the fiber post generated lower stresses along the interface and higher stresses in the root. However, with the fiber post, fracture was less likely to occur in the root, since its core and post fracture indices were higher.