Finite element analysis of steady-state natural convection problems in fluid-saturated porous media heated from below

In this paper, a progressive asymptotic approach procedure is presented for solving the steady-state Horton-Rogers-Lapwood problem in a fluid-saturated porous medium. The Horton-Rogers-Lapwood problem possesses a bifurcation and, therefore, makes the direct use of conventional finite element methods difficult. Even if the Rayleigh number is high enough to drive the occurrence of natural convection in a fluid-saturated porous medium, the conventional methods will often produce a trivial non-convective solution. This difficulty can be overcome using the progressive asymptotic approach procedure associated with the finite element method. The method considers a series of modified Horton-Rogers-Lapwood problems in which gravity is assumed to tilt a small angle away from vertical. The main idea behind the progressive asymptotic approach procedure is that through solving a sequence of such modified problems with decreasing tilt, an accurate non-zero velocity solution to the Horton-Rogers-Lapwood problem can be obtained. This solution provides a very good initial prediction for the solution to the original Horton-Rogers-Lapwood problem so that the non-zero velocity solution can be successfully obtained when the tilted angle is set to zero. Comparison of numerical solutions with analytical ones to a benchmark problem of any rectangular geometry has demonstrated the usefulness of the present progressive asymptotic approach procedure. Finally, the procedure has been used to investigate the effect of basin shapes on natural convection of pore-fluid in a porous medium. (C) 1997 by John Wiley & Sons, Ltd.

Trunculins G-I: New norsesterterpene cyclic peroxides from a southern Australian marine sponge, Latrunculia sp.

A Latrunculia sp, collected off Port Phillip Bay, Victoria, returned three new norsesterterpene cyclic peroxides. Trunculins G (9), H (10) and I (11) were isolated as their methyl esters (12), (13) and (14) respectively. Gross structures for these new trunculins were assigned on the basis of spectroscopic analysis, while the absolute stereochemistry about the cyclic peroxide terminus was established by application of the Horeau and Mosher procedures.

Nitric oxide inhibition of the rat olfactory cyclic nucleotide-gated cation channel

The effects of nitric oxide (NO) and other cysteine modifying agents were examined on cyclic nucleotide-gated (CNG) cation channels from rat olfactory receptor neurons. The NO compounds, S-nitroso-cysteine (SNC) and 3-morpholino-sydnonomine (SIN-1), did not activate the channels when applied for up to 10 min. The cysteine alkylating agent, N-ethylmaleimide (NEM), and the oxidising agent, dithionitrobensoate (DTNB), were also without agonist efficacy. Neither SNC nor DTNB altered the cAMP sensitivity of the channels. However, 2-min applications of SIN-1, SNC and DTNB inhibited the cAMP-gated current to approximately 50% of the control current level. This inhibition showed no spontaneous reversal for 5 min but was completely reversed by a 2-min exposure to DTT. The presence of cAMP protected the channels against NO-induced inhibition. These results indicate that inhibition is caused by S-nitrosylation of neighboring sulfhydryl groups leading to sulfhydryl bond formation. This reaction is favored in the closed channel state. Since recombinantly expressed rat olfactory alpha and beta CNG channel homomers and alpha/beta heteromers are activated and not inhibited by cysteine modification, the results of this study imply the existence of a novel subunit or tightly bound factor which dominates the effect of cysteine modification in the native channels. As CNG channels provide a pathway for calcum influx, the results may also have important implications for the physiological role of NO in mammalian olfactory receptor neurons.

Determination of pore size distributions by regularization and finite element collocation

The problem of extracting pore size distributions from characterization data is solved here with particular reference to adsorption. The technique developed is based on a finite element collocation discretization of the adsorption integral, with fitting of the isotherm data by least squares using regularization. A rapid and simple technique for ensuring non-negativity of the solutions is also developed which modifies the original solution having some negativity. The technique yields stable and converged solutions, and is implemented in a package RIDFEC. The package is demonstrated to be robust, yielding results which are less sensitive to experimental error than conventional methods, with fitting errors matching the known data error. It is shown that the choice of relative or absolute error norm in the least-squares analysis is best based on the kind of error in the data. (C) 1998 Elsevier Science Ltd. All rights reserved.

Finite element modelling of temperature gradient driven rock alteration and mineralization in porous rock masses

We present finite element simulations of temperature gradient driven rock alteration and mineralization in fluid saturated porous rock masses. In particular, we explore the significance of production/annihilation terms in the mass balance equations and the dependence of the spatial patterns of rock alteration upon the ratio of the roll over time of large scale convection cells to the relaxation time of the chemical reactions. Special concepts such as the gradient reaction criterion or rock alteration index (RAI) are discussed in light of the present, more general theory. In order to validate the finite element simulation, we derive an analytical solution for the rock alteration index of a benchmark problem on a two-dimensional rectangular domain. Since the geometry and boundary conditions of the benchmark problem can be easily and exactly modelled, the analytical solution is also useful for validating other numerical methods, such as the finite difference method and the boundary element method, when they are used to dear with this kind of problem. Finally, the potential of the theory is illustrated by means of finite element studies related to coupled flow problems in materially homogeneous and inhomogeneous porous rock masses. (C) 1998 Elsevier Science S.A. All rights reserved.

The solution structure of a copper(II) compound of a new cyclic octapeptide by EPR spectroscopy and force field calculations

A new cyclic octapeptide, cyclo(Ile-Ser-(Gly)Thz-Ile-Thr-(Gly)Thz) (PatN), related to patellamide A, has been synthesized and reacted with copper(II) and base to form mono- and dinuclear complexes. The coordination environments around copper(TI) have been characterized by EPR spectroscopy. The solution structure of the thermodynamically most stable product, a purple dicopper(TI) compound, has been examined by simulating weakly dipole-dipole coupled EPR spectra based upon structural parameters obtained from force field (MM and MD) calculations. The MM-EPR method produces a saddle-shaped structure for [Cu-2(PatN)(OH2)(6)] that is similar to the known solution structure of patellamide A and the known solid-state structure of [Cu-2(AscidH(2))CO3(OH2)(2)]. Compared with the latter, [Cu-2(PatN)] has no carbonate bridge and a significantly flatter topology. The MM-EPR approach to solution-structure determination for paramagnetic metallopeptides may find wide applications to other metallopeptides and metalloproteins.

Nuapapuin A and sigmosceptrellins D and E: New norterpene cyclic peroxides from a southern Australian marine sponge, Sigmosceptrella sp.

A Sigmosceptrella sp. from the Great Australian Eight, Australia, has yielded the new norditerpene cyclic peroxide, nuapapuin A (2a), and the norsesterterpene cyclic peroxide sigmosceptrellin D (3a), characterized as the corresponding methyl esters 2b and 3b. The crude methylated sponge extract also yielded the new norsesterterpene cyclic peroxide sigmosceptrellin E methyl ester (4). Relative stereochemistry about C2, C3, and C6 was assigned by established empirical rules and absolute stereochemistry by the advanced Mosher procedure. A plausible biosynthetic pathway has been proposed that rationalizes key transformations in the biosynthesis of all known norterpene cyclic peroxides and related norterpene ketones, dienes and sigmosceptrins.

Interaction of zinc(II) with the cyclic octapeptides, cyclo[Ile(Oxn)-D-Val(Thz)](2) and ascidiacyclamide, a cyclic peptide from Lissoclinum patella

The interactions between zinc salts and the naturally occurring cyclic octapeptide ascidiacyclamide in methanol, as well as a synthetic analogue cyclo[Ile(Oxn)-D-Val(Thz)](2), were monitored by H-1 NMR and CD spectroscopy. Three zinc complexes were identified, their relative amounts depending on the nature of the anion (perchlorate, triflate or chloride) and the presence or absence of base. Binding constants for two of the zinc species were calculated from CD or H-1 NMR spectra, [Zn(L - H)](+) (KZn(L-H) = [Zn(L - H)(+)]/[Zn2+][(L - H)(-)] = 10(7 +/- 2) M-1; 95% methanol/5% water, 298.0 K, NEt3/HClO4 buffer 0.04 M) and [ZnLCl](+) (K-ZnCIL = [ZnCIL+]/[Zn2+][Cl-][L] = 10(7.2) (+/-) (0.1) M-2; d(3)-methanol, 301 K).

Low-molecular-weight peptidic and cyclic antagonists of the receptor for the complement factor C5a

Activation of the human complement system of plasma proteins during immunological host defense can result in overproduction of potent proinflammatory peptides such as the anaphylatoxin C5a. Excessive levels of C5a are associated with numerous immunoinflammatory diseases, but there is as yet no clinically available antagonist to regulate the effects of C5a. We now describe a series of small molecules derived from the C-terminus of C5a, some of which are the most potent low-molecular-weight C5a receptor antagonists reported to date for the human polymorphonuclear leukocyte (PMN) C5a receptor. H-1 NMR spectroscopy was used to determine solution structures for two cyclic antagonists and to indicate that antagonism is related to a turn conformation, which can be stabilized in cyclic molecules that are preorganized for receptor binding. While several cyclic derivatives were of similar antagonistic potency, the most potent antagonist was a hexapeptide-derived macrocycle AcF[OPdChaWR] with an IC50 = 20 nM against a maximal concentration of C5a (100 nM) on intact human PMNs. Such potent C5a antagonists may be useful probes to investigate the role of C5a in host defenses and to develop therapeutic agents for the treatment of many currently intractable inflammatory conditions.

Binding of YY1 and Oct1 to a novel element that downregulates expression of IL-5 in human T cells

Background: IL-5 controls development of eosinophilia and has been shown to be involved in the pathogenesis of allergic diseases. In both atopic and nonatopic asthma, elevated IL-5 has been detected in peripheral blood and the airways. IL-5 is produced mainly by activated T cells, and its expression is regulated at the transcriptional level. Objective: This study focuses on the functional analysis of the human IL-5 (hIL-5) promoter and characterization of eis-regulatory elements and transcription factors involved in the suppression of IL-5 transcription in T cells. Methods: Methods used in this study include DNase I footprint assays, electrophoretic mobility shift assays, and functional analysis by mammalian cell transfection involving deletion analysis and site-directed mutagenesis. Results: We identified 5 protein binding regions (BRs) located within the proximal hIL-5 promoter. Functional analysis indicates that the BRs are involved in control of hIL-5 promoter activity. Two of these regions, BR3 and BR4 located at positions -102 to -73, have not previously been described as regulators of IL-5 expression in T cells. We show that the BR3 sequence contains a novel negative regulatory element located at positions -90 to -79 of the hIL-5 promoter, which binds Oct1, octamer-like, and YY1 nuclear factors. Substitution mutations, which abolished binding of these proteins to the BR3 sequence, significantly increased hIL-5 promoter activity in activated T cells. Conclusion: We suggest that Oct1, YY1, and octamer-like factors binding to the -90/-79 sequence within the proximal IL-5 promoter are involved in suppression of IL-5 transcription in T cells.

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

Finite element modelling of reactive mass transport problems in fluid-saturated porous media

We use the finite element method to solve reactive mass transport problems in fluid-saturated porous media. In particular, we discuss the mathematical expression of the chemical reaction terms involved in the mass transport equations for an isothermal, non-equilibrium chemical reaction. It has turned out that the Arrhenius law in chemistry is a good mathematical expression for such non-equilibrium chemical reactions especially from the computational point of view. Using the finite element method and the Arrhenius law, we investigate the distributions of PH (i.e. the concentration of H+) and the relevant reactive species in a groundwater system. Although the main focus of this study is on the contaminant transport problems in groundwater systems, the related numerical techniques and principles are equally applicable to the orebody formation problems in the geosciences. Copyright (C) 1999 John Wiley & Sons, Ltd.

Rectification of the olfactory cyclic nucleotide-gated channel by intracellular polyamines

Polyamine-induced inward rectification of cyclic nucleotide-gated channels was studied in inside-out patches from rat olfactory neurons. The polyamines, spermine, spermidine and putrescine, induced an 'instantaneous' voltage-dependent inhibition with K-d values at 0 mV of 39, 121 mu M and 2.7 mM, respectively. Hill coefficients for inhibition were significantly < 1, suggesting an allosteric inhibitory mechanism. The Woodhull model for voltage-dependent block predicted that all 3 polyamines bound to a site 1/3 of the electrical distance through the membrane from the internal side. Instantaneous inhibition was relieved at positive potentials, implying significant polyamine permeation. Spermine also induced exponential current relaxations to a 'steady-state' impermeant level. This inhibition was also mediated by a binding site 1/3 of the electrical distance through the pore, but with a K-d of 2.6 mM. Spermine inhibition was explained by postulating two spermine binding sites at a similar depth. Occupation of the first site occurs rapidly and with high affinity, but once a spermine molecule has bound, it inhibits spermine occupation of the second binding site via electrostatic repulsion. This repulsion is overcome at higher membrane potentials, but results in a lower apparent binding affinity for the second spermine molecule. The on-rate constant for the second spermine binding saturated at a low rate (similar to 200 sec(-1) at +120 mV), providing further evidence for an allosteric mechanism. Polyamine-induced inward rectification was significant at physiological concentrations.