Finite element modelling of three-dimensional convection problems in fluid-saturated porous media heated from below

We use the finite element method to model three-dimensional convective pore-fluid flow in fluid-saturated porous media when they are heated from below. In particular, we employ the particle-tracking technique to mimic the trajectories of particles in three-dimensional fluid flow problems. The related numerical results demonstrated that: (1) The progressive asymptotic approach procedure, which was previously developed for the finite element modelling of two-dimensional convective pore-fluid flow problems, is equally applicable to the finite element modelling of three-dimensional convective pore-fluid flow in fluid-saturated porous media heated from below. (2) The perturbation of gravity at different planes has a significant effect on the pattern of three-dimensional convective pore-fluid flow and therefore, may influence the pattern of orebody formation and mineralization in three-dimensional hydrothermal systems. Copyright (C) 2001 John Wiley & Sons, Ltd.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Simultaneous steady state and dynamic design of process systems using structural indicators

The simultaneous design of the steady-state and dynamic performance of a process has the ability to satisfy much more demanding dynamic performance criteria than the design of dynamics only by the connection of a control system. A method for designing process dynamics based on the use of a linearised systems' eigenvalues has been developed. The eigenvalues are associated with system states using the unit perturbation spectral resolution (UPSR), characterising the dynamics of each state. The design method uses a homotopy approach to determine a final design which satisfies both steady-state and dynamic performance criteria. A highly interacting single stage forced circulation evaporator system, including control loops, was designed by this method with the goal of reducing the time taken for the liquid composition to reach steady-state. Initially the system was successfully redesigned to speed up the eigenvalue associated with the liquid composition state, but this did not result in an improved startup performance. Further analysis showed that the integral action of the composition controller was the source of the limiting eigenvalue. Design changes made to speed up this eigenvalue did result in an improved startup performance. The proposed approach provides a structured way to address the design-control interface, giving significant insight into the dynamic behaviour of the system such that a systematic design or redesign of an existing system can be undertaken with confidence.

Finite element modelling of rock alteration and metamorphic process in hydrothermal systems

We use the finite element method to simulate the rock alteration and metamorphic process in hydrothermal systems. In particular, we consider the fluid-rock interaction problems in pore-fluid saturated porous rocks. Since the fluid rock interaction takes place at the contact interface between the pore-fluid and solid minerals, it is governed by the chemical reaction which usually takes place very slowly at this contact interface, from the geochemical point of view. Due to the relative slowness of the rate of the chemical reaction to the velocity of the pore-fluid flow in the hydrothermal system to be considered, there exists a retardation zone, in which the conventional static theory in geochemistry does not hold true. Since this issue is often overlooked by some purely numerical modellers, it is emphasized in this paper. The related results from a typical rock alteration and metamorphic problem in a hydrothermal system have shown not only the detailed rock alteration and metamorphic process, but also the size of the retardation zone in the hydrothermal system. Copyright (C) 2001 John Wiley & Sons, Ltd.

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.

On entanglement and Lorentz transformations

We study the transformation of maximally entangled states under the action of Lorentz transformations in a fully relativistic setting. By explicit calculation of the Wigner rotation, we describe the relativistic analog of the Bell states as viewed from two inertial frames moving with constant velocity with respect to each other. Though the finite dimensional matrices describing the Lorentz transformations are non-unitary, each single particle state of the entangled pair undergoes an effective, momentum dependent, local unitary rotation, thereby preserving the entanglement fidelity of the bipartite state. The details of how these unitary transformations are manifested are explicitly worked out for the Bell states comprised of massive spin 1/2 particles and massless photon polarizations. The relevance of this work to non-inertial frames is briefly discussed.

Scope for further analytical solutions for constant flux infiltration into a semi-infinite soil profile or redistribution in a finite soil profile

[1] We attempt to generate new solutions for the moisture content form of the one-dimensional Richards' [1931] equation using the Lisle [1992] equivalence mapping. This mapping is used as no more general set of transformations exists for mapping the one-dimensional Richards' equation into itself. Starting from a given solution, the mapping has the potential to generate an infinite number of new solutions for a series of nonlinear diffusivity and hydraulic conductivity functions. We first seek new analytical solutions satisfying Richards' equation subject to a constant flux surface boundary condition for a semi-infinite dry soil, starting with the Burgers model. The first iteration produces an existing solution, while subsequent iterations are shown to endlessly reproduce this same solution. Next, we briefly consider the problem of redistribution in a finite-length soil. In this case, Lisle's equivalence mapping is generalized to account for arbitrary initial conditions. As was the case for infiltration, however, it is found that new analytical solutions are not generated using the equivalence mapping, although existing solutions are recovered.

John Waldridge, The Leaven of the Ancients: Suhrawardi and the Heritage of the Greeks, SUNY Series in Islam (Albany: State University of New York Press, 1999). Pp. 323

In the Leaven of the Ancients, John Walbridge studies the appropriation of non–Peripatetic philosophical ideas by an anti–Aristotelian Islamic philosopher, Shihab al-Din al-Suhrawardi (d. 1191). He proposes a comprehensive explanation of the origin of Suhrawardi's philosophical system, a revival of the “wisdom of the Ancients” and its philosophical affiliations “grounded” in Greek philosophy (p. xiii). Walbridge attempts to uncover the reasons for Suhrawardi's rejection of the prevailing neo–Aristotelian synthesis in Islamic philosophy, Suhrawardi's knowledge and understanding of non–Aristotelian Greek philosophy, the ancient philosophers Suhrawardi was attempting to follow, the relationship between Suhrawardi's specific philosophical teachings (logic, ontology, physics, and metaphysics), and his understanding of non–Aristotelian ancient philosophy and the relationship between Suhrawardi's system and the major Greek philosophers, schools, and traditions—in particular the Presocratics, Plato, and the Stoics (p. 8). Copyright © 2003 Cambridge University Press

Finite-difference time-domain-based studies of MRI pulsed field gradient-induced eddy currents inside the human body

In modern magnetic resonance imaging (MRI), patients are exposed to strong, rapidly switching magnetic gradient fields that, in extreme cases, may be able to elicit nerve stimulation. This paper presents theoretical investigations into the spatial distribution of induced current inside human tissues caused by pulsed z-gradient fields. A variety of gradient waveforms have been studied. The simulations are based on a new, high-definition, finite-difference time-domain method and a realistic inhomogeneous 10-mm resolution human body model with appropriate tissue parameters. it was found that the eddy current densities are affected not only by the pulse sequences but by many parameters such as the position of the body inside the gradient set, the local biological material properties and the geometry of the body. The discussion contains a comparison of these results with previous results found in the literature. This study and the new methods presented herein will help to further investigate the biological effects caused by the switched gradient fields in a MRI scan. (C) 2002 Wiley Periodicals, Inc.

Bringing the state to England: Andrew Tooke's translation of Samuel Pufendorf's 'De officio hominis et civis'

Andrew Tooke's 1691 English translation of Samuel Pufendorf's De officio hominis et civis, published as The Whole Duty of Man According to the Law of Nature, brought Pufendorf's manual fo statist natural law into English politics at a moment of temporary equilibrium in the unfinished contest between Crown and Parliament for the rights and powers of sovereignty. Drawing on the authors' re-edition of The Whole Duty of Man, this article describes and analyses a telling instance of how--by translation--the core political terms and concepts of the German natural jurist's 'absolutist' formulary were reshaped for reception in the different political culture of late seventeenth-century England.

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.

Solid-State Experimental and Theoretical Investigation of the Ammonium Salt of Croconate Violet, a Pseudo-Oxocarbon Ion

The present work describes the crystal structure, vibrational spectra, and theoretical calculations of ammonium salts of 3,5-bis-(dicyanomethylene)cyclopentane-1,2,4-trionate, (NH(4))(2)(C(11)N(4)O(3)) [(NH(4))(2)CV], also known as ammonium croconate violet. This compound crystallizes in triclinic P (1) over bar and contains two water molecules per unit formula. The crystal packing is stabilized by hydrogen bonds involving water molecules and ammonium cations, giving rise to a 3D polymeric arrangement. In this structure, a pi-stacking interaction is not observed, as the smaller centroid-centroid distance is 4.35 angstrom. Ab initio electronic structure calculations under periodic boundary conditions were performed to predict vibrational and electronic properties. The vibrational analysis was used to assist the assignments of the Raman and infrared bands. The solid structure was optimized and characterized as a minimum in the potential-energy surface. The stabilizing intermolecular hydrogen bonds in the crystal Structure were characterized by difference charge-density analysis. The analysis of the density of states of (NH(4))(2)CV gives an energy gap of 1.4 eV with a significant contribution of carbon and nitrogen 2p states for valence and conduction bands.

Production efficiency and financing of public health: an analysis of small municipalities in the state of Sao Paulo - Brazil

This paper measured the variations in performance of small municipalities in the State of Sao Paulo, Brazil, regarding the technical efficiency in the use of public funds in public primary health care actions concerning the funding profile, in a scenario of fiscal federalism. Technical efficiency is one of the parameters of evaluation of public sector performance and was measured by means of Data Envelopment Analysis (DEA). The correlation analysis of DEA score was used to verify possible associations between technical efficiency and the funding profile of expenses with health care. The results showed that 6.41% of the municipalities were considered efficient. They also showed that the level of municipality dependence to inter-governmental general purpose grants and the national health funding specific purpose grants have negative correlation with efficiency scores.

Systemic Model for Diagnosis of the Micro, Small and Medium Enterprises from Two Cities from the Countryside of the State of So Paulo in Brazil

This article presents a proposal of a systemic model composed for the micro and small companies (MSE) of the region of Ribeiro Preto and the agents which influenced their environment. The proposed model was based on Stafford Beer`s (Diagnosing the system for organizations. Chichester, Wiley, 1985) systemic methodologies VSM (Viable System Model) and on Werner Ulrich`s (1983) CSH (Critical Systems Heuristics). The VSM is a model for the diagnosis of the structure of an organization and of its flows of information through the application of the cybernetics concepts (Narvarte, In El Modelo del Sistema Viable-MSV: experiencias de su aplicacin en Chile. Proyecto Cerebro Colectivo del IAS, Santiago, 2001). On the other hand, CSH focus on the context of the social group applied to the systemic vision as a counterpoint to the organizational management view considered by the VSM. MSE of Ribeiro Preto and Sertozinho had been analyzed as organizations inserted in systems that relate and integrate with other systems concerning the public administration, entities of representation and promotion agencies. The research questions: which are the bonds of interaction among the subsystems in this process and who are the agents involved? The systemic approach not only diagnosed a social group, formed by MSE of Ribeiro Preto and Sertozinho, public authorities and support entities, but could also delineate answers that aimed the clarification of obscure questions generating financial assistance to the formularization of efficient actions for the development of this system.

Comparative analysis of population structure of endemic Western Atlantic hermit crab Loxopagurus loxochelis (Decapoda, Anomura) in two regions of state of Sao Paulo, Brazil

This work characterized the population structure of the hermit crab Loxopagurus loxochelis (Moreira, 1901) in terms of size frequency distribution and sex ratio. Specimens were collected monthly, over a period of one year (from July 2002 to June 2003), in seven transects (from 5 to 35 m of depth) using fishing boat equipped with two double-rig trawl nets, in Caraguatatuba and Ubatuba regions (state of Sao Paulo, Brazil). A total of 366 hermit crabs were collected in Caraguatatuba [222 males (60.65%), 114 non-ovigerous females (31.15%) and 30 ovigerous females (8.20%)] and 126 hermit crabs in Ubatuba [81 males (64.28%), 38 non-ovigerous females (30.16%) and seven ovigerous females (5.56%)]. In Caraguatatuba the highest incidence of ovigerous females occurred during winter (July 2002), whereas in Ubatuba, the number was incipient. The cephalothoracic shield length ranged from 2.0 to 7.9mm (5.29 +/- 0.96mm) in Caraguatatuba, and from 2.7 to 7.5mm (5.32 +/- 0.95mm) in Ubatuba. The mean size of males was significantly larger than the mean size of females in both regions. Overall sex ratio was in favor of males (1.54:1 in Caraguatatuba and 1.9:1 in Ubatuba). Sexual dimorphism was recorded to L. loxochelis by the presence of males in the largest size classes, following the standard pattern observed in Decapoda. There was an unimodal size distribution for both sexes, with normal distributions in both regions. The higher number of males in relation to females may indicate the existence of different growth and mortality rates between the sexes. Despite of the different geomorphologic characteristics between Caraguatatuba and Ubatuba regions, the dynamics of development was similar for both populations.