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.

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.

A model of seller holdout

We model a buyer who wishes to combine objects owned by two separate sellers in order to realize higher value. Sellers are able to avoid entering into negotiations with the buyer, so that the order in which they negotiate is endogenous. Holdout occurs if at least one of the sellers is not present in the first round of negotiations. We demonstrate that complementarity of the buyer's technology is a necessary condition for equilibrium holdout. Moreover, a rise in complementarity leads to an increased likelihood of holdout, and an increased efficiency loss. Applications include patents, the land assembly problem, and mergers.

Demand elasticities for food products in Brazil: a two-stage budgeting system

The object of this article is to estimate demand elasticities for a basket of staple food important for providing the caloric needs of Brazilian households. These elasticities are useful in the measurement of the impact of structural reforms on poverty. A two-stage demand system was constructed, based on data from Household Expenditure Surveys (POF) produced by IBGE (The Brazilian Bureau of Statistics) in 1987/88 and 1995/96. We have used panel data to estimate the model, and have calculated income, own-price, and cross-price elasticities for eight groups of goods and services and, in the second stage, for 11 sub groups of staple food products. We estimated those elasticities for the whole sample of consumers and for two income groups.

Thermal Characteristics of Polybia scutellaris Nests (Hymenoptera: Vespidae) Using Computational Fluid Dynamics: A Possible Adaptation to Tropical Climates

Polybia scutellaris constructs huge nests characterized by numerous spinal projections on the surface. We investigated the thermal characteristics of P scutellaris nests in order to determine whether their nest temperature is homeothermically maintained and whether the spines play a role in the thermoregulation of the nests. In order to examine these hypotheses, we measured the nest temperature in a active nest and in an abandoned nest. The temperature in the active nest was almost stable at 27 degrees C, whereas that of the abandoned nest varied with changes in the ambient temperature, suggesting that nest temperature was maintained by the thermogenesis of colony individuals. In order to predict the thermal properties of the spines, a numerical simulation was employed. To construct a 3D-model of a P scutellaris nest, the nest architecture was simplified into an outer envelope and the surface spines, for both of which the initial temperature was set at 27 degrees C. The physical properties of the simulated nest were regarded to be those of wood since the nest of this species is constructed from plant materials. When the model was exposed to cool air (12 degrees C), the temperature was lower in the models with more spines. On the other hand, when the nest was heated (42 degrees C), the temperature increase was smaller in models with more spines. It is suggested that the spines act as a heat radiator, not as an insulator, against the changes in ambient temperature.

Matching occupation and self: Does matching theory adequately model children's thinking?

The present exploratory-descriptive cross-national study focused on the career development of 11- to 14-yr.-old children, in particular whether they can match their personal characteristics with their occupational aspirations. Further, the study explored whether their matching may be explained in terms of a fit between person and environment using Holland's theory as an example. Participants included 511 South African and 372 Australian children. Findings relate to two items of the Revised Career Awareness Survey that require children to relate personal-social knowledge to their favorite occupation. Data were analyzed in three stages using descriptive statistics, i.e., mean scores, frequencies, and percentage agreement. The study indicated that children perceived their personal characteristics to be related to their occupational aspirations. However, how this matching takes place is not adequately accounted for in terms of a career theory such as that of Holland.

Chitosan as a Bioadhesive Agent Between Porphyrins and Phospholipids in a Biomembrane Model

Porphyrins are currently used in photodynamic therapy as photosensitizers. In this paper we studied the interaction of two charged porphyrins, 5, 10, 15, 20-mesotetrakis(N-metyl-4-pyridyl) porphyrin, (TMPyP/chloride salt) cationic, and 5, 10, 15, 20-meso-tetrakis(sulfonatophenyl) porphyrin, (TPPS(4)/sodium salt) anionic, nanoassembled in phospholipid Langmuir monolayers and Langmuir-Blodgett films. Furthermore, we used chitosan to mediate the interaction between the porphyrins and the model membrane, aiming to understand the role of the polysaccharide in a molecular level. The effect of the interaction of the photosensitizers on the fluidity of the lipid monolayer was investigated by using dilatational surface elasticity. We also used photoluminescence (PL) spectroscopy to identify the porphyrins adsorbed in the phospholipid films. We observed an expansion of the monolayer promoted by the adsorption of the porphyrins into the lipid-air interface which was more pronounced in the case of TMPyP, as a consequence of a strong electrostatic interaction with the anionic monolayer. The chitosan promoted a higher adsorption of the porphyrins on the phospholipid monolayers and enabled the porphyrin to stay in its monomeric form (as confirmed by PL spectroscopy), thus demonstrating that chitosan can be pointed out as a potential photosensitizer delivery system in photodynamic therapy.

Non-diagonal solutions of the reflection equation for the trigonometric A((1)) (n-1) vertex model

We obtain a class of non-diagonal solutions of the reflection equation for the trigonometric A(n-1)((1)) vertex model. The solutions can be expressed in terms of intertwinner matrix and its inverse, which intertwine two trigonometric R-matrices. In addition to a discrete (positive integer) parameter l, 1 less than or equal to l less than or equal to n, the solution contains n + 2 continuous boundary parameters.

A model of migrants' remittances with human capital investment and intrafamilial transfers

This study analyzes data on migrants' remittances using a two-period theory of intergenerational transfers based on an informal, intrafamilial loan arrangement using weak altruism, a behavior between strong altruism and pure self-interest. The model provides an integrated theory of migrants' remittances, human capital investment decisions, and intrafamilial transfers applicable to low-income countries with no official pension schemes and imperfect capital markets. Propositions, derived from the theory, are tested, re-analyzing original survey data on remittances of Pacific island migrants in Sydney. When weak altruism and strong altruism yield opposite predictions, the econometric results tend to confirm the former hypothesis and invalidate the latter.

Amphetamine modulates cellular recruitment and airway reactivity in a rat model of allergic lung inflammation

Asthma is characterized by pulmonary cellular infiltration, vascular exudation and airway hyperresponsiveness. Several drugs that modify central nervous system (CNS) activity can modulate the course of asthma. Amphetamine (AMPH) is a highly abused drug that presents potent stimulating effects on the CNS and has been shown to induce behavioral, biochemical and immunological effects. The purpose of this study was to investigate the effects of AMPH on pulmonary cellular influx, vascular permeability and airway reactivity. AMPH effects on adhesion molecule expression, IL-10 and IL-4 release and mast cell degranulation were also studied. Male Wistar rats were sensitized with ovalbumin (OVA) plus alum via subcutaneous injection. One week later, the rats received another injection of OVA-alum (booster). Two weeks after this booster, the rats were subjected to AMPH treatment 12 h prior to the OVA airway challenge. In rats treated with AMPH, the OVA challenge reduced cell recruitment into the lung, the vascular permeability and the cellular expression of ICAM-1 and Mac-1. Additionally, elevated levels of IL-10 and IL-4 were found in samples of lung explants from allergic rats. AMPH treatment, in comparison, increased IL-10 levels but reduced those of IL-4 in the lung explants. Moreover, the tracheal responsiveness to methacholine (MCh), as well as to an in vitro OVA challenge, was reduced by AMPH treatment, and levels of PCA titers were not modified by the drug. Our findings suggest that single AMPH treatment down-regulates several parameters of lung inflammation, such as cellular migration, vascular permeability and tracheal responsiveness. These results also indicate that AMPH actions on allergic lung inflammation include endothelium-leukocyte interaction mechanisms, cytokine release and mast cell degranulation. (C) 2010 Elsevier Ireland Ltd. All rights reserved.

Developmental vitamin D deficiency alters adult behaviour: An informative animal model of schizophrenia

Background: There is growing evidence that vitamin D is active in the brain but until recently there was a lack of evidence about its role during brain development. Guided by certain features of the epidemiology of schizophrenia, we have explored the role of vitamin D in the developing brain and behaviour using whole animal models. Methods: Sprague-Dawley rats were fed a vitamin D deficient diet (DVD) or control diet 6 weeks prior to mating and housed under UVB-free lighting conditions. On the day of birth all rats were fed a control diet for the remainder of the study. We observed behaviour at two timepoints; on the day of birth to study maternal behaviour, and at 10 weeks of age to study offspring behaviour in adulthood, under baseline and drug induced conditions (MK-801, haloperidol, amphetamine). Results: Prenatal vitamin D deficiency results in subtle alterations in maternal behaviour as well as long lasting effects on the adult offspring, despite a return to normal vitamin D levels during postnatal life. These affects were specific to transient prenatal vitamin D depletion as adult vitamin D depletion, combined prenatal and chronic postnatal vitamin D depletion, or ablation of the vitamin D receptor in mice led to markedly different outcomes. Conclusions: The developmental vitamin D (DVD) model now draws strength from epidemiological evidence of schizophrenia and animal experiments. Although the DVD model does not replicate every aspect of schizophrenia, it has several attractive features: (1) the exposure is based on clues from epidemiology; (2) it reproduces the increase in lateral ventricles; (3) it reproduces well-regarded behavioural phenotypes associated with schizophrenia (e.g. MK- 801 induced hyperlocomotion); and (4) it implicates a disturbance in dopamine signaling. In summary, low prenatal levels of vitamin D can influence critical components of orderly brain development and that this has a long lasting effect on behaviour.

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.

Value of real time three-dimensional echocardiography in patients with hypertrophic cardiomyopathy: Comparison with two-dimensional echocardiography and magnetic resonance imaging

Real time three-dimensional echocardiography (RT3DE) has been demonstrated to be an accurate technique to quantify left ventricular (LV) volumes and function in different patient populations. We sought to determine the value of RT3DE for evaluating patients with hypertrophic cardiomyopathy (HCM), in comparison with cardiac magnetic resonance imaging (MRI). Methods: We studied 20 consecutive patients with HCM who underwent two-dimensional echocardiography (2DE), RT3DE, and MRI. Parameters analyzed by echocardiography and MRI included: wall thickness, LV volumes, ejection fraction (LVEF), mass, geometric index, and dyssynchrony index. Statistical analysis was performed by Lin agreement coefficient, Pearson linear correlation and Bland-Altman model. Results: There was excellent agreement between 2DE and RT3DE (Rc = 0.92), 2DE and MRI (Rc = 0.85), and RT3DE and MRI (Rc = 0.90) for linear measurements. Agreement indexes for LV end-diastolic and end-systolic volumes were Rc = 0.91 and Rc = 0.91 between 2DE and RT3DE, Rc = 0.94 and Rc = 0.95 between RT3DE and MRI, and Rc = 0.89 and Rc = 0.88 between 2DE and MRI, respectively. Satisfactory agreement was observed between 2DE and RT3DE (Rc = 0.75), RT3DE and MRI (Rc = 0.83), and 2DE and MRI (Rc = 0.73) for determining LVEF, with a mild underestimation of LVEF by 2DE, and smaller variability between RT3DE and MRI. Regarding LV mass, excellent agreement was observed between RT3DE and MRI (Rc = 0.96), with bias of -6.3 g (limits of concordance = 42.22 to -54.73 g). Conclusion: In patients with HCM, RT3DE demonstrated superior performance than 2DE for the evaluation of myocardial hypertrophy, LV volumes, LVEF, and LV mass.

AT(1) blockade during lactation as a model of chronic nephropathy: mechanisms of renal injury

Suppression of the renin-angiotensin system during lactation causes irreversible renal structural changes. In this study we investigated 1) the time course and the mechanisms underlying the chronic kidney disease caused by administration of the AT(1) receptor blocker losartan during lactation, and 2) whether this untoward effect can be used to engender a new model of chronic kidney disease. Male Munich-Wistar pups were divided into two groups: C, whose mothers were untreated, and L(Lact), whose mothers received oral losartan (250 mg.kg(-1).day(-1)) during the first 20 days after delivery. At 3 mo of life, both nephron number and the glomerular filtration rate were reduced in L(Lact) rats, whereas glomerular pressure was elevated. Unselective proteinuria and decreased expression of the zonula occludens-1 protein were also observed, along with modest glomerulosclerosis, significant interstitial expansion and inflammation, and wide glomerular volume variation, with a stable subpopulation of exceedingly small glomeruli. In addition, the urine osmolality was persistently lower in L(Lact) rats. At 10 mo of age, L(Lact) rats exhibited systemic hypertension, heavy albuminuria, substantial glomerulosclerosis, severe renal interstitial expansion and inflammation, and creatinine retention. Conclusions are that 1) oral losartan during lactation can be used as a simple and easily reproducible model of chronic kidney disease in adult life, associated with low mortality and no arterial hypertension until advanced stages; and 2) the mechanisms involved in the progression of renal injury in this model include glomerular hypertension, glomerular hypertrophy, podocyte injury, and interstitial inflammation.

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier Ltd. All rights reserved.