Numerical prediction of gas concentrations and fluctuations above a triangular hill within a turbulent boundary layer

The objective of the present work is to propose a numerical and statistical approach, using computational fluid dynamics, for the study of the atmospheric pollutant dispersion. Modifications in the standard k-epsilon turbulence model and additional equations for the calculation of the variance of concentration are introduced to enhance the prediction of the flow field and scalar quantities. The flow field, the mean concentration and the variance of a flow over a two-dimensional triangular hill, with a finite-size point pollutant source, are calculated by a finite volume code and compared with published experimental results. A modified low Reynolds k-epsilon turbulence model was employed in this work, using the constant of the k-epsilon model C(mu)=0.03 to take into account the inactive atmospheric turbulence. The numerical results for the velocity profiles and the position of the reattachment point are in good agreement with the experimental results. The results for the mean and the variance of the concentration are also in good agreement with experimental results from the literature. (C) 2009 Elsevier Ltd. All rights reserved.

Dynamic Imaging in Electrical Impedance Tomography of the Human Chest With Online Transition Matrix Identification

One of the electrical impedance tomography objectives is to estimate the electrical resistivity distribution in a domain based only on electrical potential measurements at its boundary generated by an imposed electrical current distribution into the boundary. One of the methods used in dynamic estimation is the Kalman filter. In biomedical applications, the random walk model is frequently used as evolution model and, under this conditions, poor tracking ability of the extended Kalman filter (EKF) is achieved. An analytically developed evolution model is not feasible at this moment. The paper investigates the identification of the evolution model in parallel to the EKF and updating the evolution model with certain periodicity. The evolution model transition matrix is identified using the history of the estimated resistivity distribution obtained by a sensitivity matrix based algorithm and a Newton-Raphson algorithm. To numerically identify the linear evolution model, the Ibrahim time-domain method is used. The investigation is performed by numerical simulations of a domain with time-varying resistivity and by experimental data collected from the boundary of a human chest during normal breathing. The obtained dynamic resistivity values lie within the expected values for the tissues of a human chest. The EKF results suggest that the tracking ability is significantly improved with this approach.

Impedimetric immunosensor for electronegative low density lipoprotein (LDL(-)) based on monoclonal antibody adsorbed on (polyvinyl formal)-gold nanoparticles matrix

Monoclonal antibodies (MAb) have been commonly applied to measure LDL in vivo and to characterize modifications of the lipids and apoprotein of the LDL particles. The electronegative low density lipoprotein (LDL(-)) has an apolipoprotein B-100 modified at oxidized events in vivo. In this work, a novel LDL-electrochemical biosensor was developed by adsorption of anti-LDL(-) MAb on an (polyvinyl formal)-gold nanoparticles (PVF-AuNPs)-modified gold electrode. Electrochemical impedance spectroscopy (EIS) and cyclic voltammetry (CV) were used to characterize the recognition of LDL-. The interaction between MAb-LDL(-) leads to a blockage in the electron transfer of the [Fe(CN)(6)](4-)/K(4)[Fe(CN)(6)](3-) redox couple, which may could result in high change in the electron transfer resistance (R(CT)) and decrease in the amperometric responses in CV analysis. The compact antibody-antigen complex introduces the insulating layer on the assembled surface, which increases the diameter of the semicircle, resulting in a high R(CT), and the charge transferring rate constant k(0) decreases from 18.2 x 10(-6) m/s to 4.6 x 10(-6) m/s. Our results suggest that the interaction between MAb and lipoprotein can be quantitatively assessed by the modified electrode. The PVF-AuNPs-MAb system exhibited a sensitive response to LDL(-), which could be used as a biosensor to quantify plasmatic levels of LDL(-). (C) 2011 Elsevier B.V. All rights reserved.

Existence of Multiple Solutions for Second Order Boundary Value Problems

We give conditions on f involving pairs of lower and upper solutions which lead to the existence of at least three solutions of the two point boundary value problem y" + f(x, y, y') = 0, x epsilon [0, 1], y(0) = 0 = y(1). In the special case f(x, y, y') = f(y) greater than or equal to 0 we give growth conditions on f and apply our general result to show the existence of three positive solutions. We give an example showing this latter result is sharp. Our results extend those of Avery and of Lakshmikantham et al.

The interaction of metal ions and Marimastat with matrix metalloproteinase 9

The effect of a range of metal ions on the ability of Marimastat to inhibit matrix metalloproteinase 9 (MMP-9) was examined in a fluorescence based proteolytic assay. Whilst none of the metals examined significantly affected the inhibitory ability of Marimastat, several metal ions did have a significant effect on MMP-9 activity itself. In the absence of Marimastat, Zn(II) and Fe(II) significantly inhibited MMP-9 activity at metal ion concentrations of 10 and 100 muM, respectively. In both the absence and presence of Marimastat, Cd(II) significantly inhibited MMP-9 at 100 muM. In contrast, 1 mM Co(II) significantly upregulated MMP-9 proteolytic activity. (C) 2003 Elsevier Science Inc. All rights reserved.

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.

Mechanisms underlying the biphasic effect of vitamin K-1 (phylloquinone) on arterial blood pressure

Phylloquinone (vitamin K-1, VK1) is widely used therapeutically and intravenous administration of this quinone can induce hypotension. We aimed to investigate the mechanisms underlying the effects induced by VK1 on arterial blood pressure. With this purpose a catheter was inserted into the abdominal aorta of male Wistar rats for blood pressure and heart rate recording. Bolus intravenous injection of VK1 (0.5-20 mg kg(-1)) produced a transient increase in blood pressure followed by a fall. Both the pressor and depressor response induced by VK1 were dose-dependent. On the other hand, intravenous injection of VK1 did not alter heart rate. The nitric oxide synthase (NOS) inhibitor N-G-nitro-L-arginine methyl ester (L-NAME, 10 and 20 mg kg(-1)) reduced both the increase and decrease in blood pressure induced by VK1 (5 mgkg(-1)). On the other hand, indometacin (10 mg kg(-1)), a non-selective cyclooxygenase inhibitor, did not alter the increase in mean arterial pressure (MAP) induced by VK1. However, VK1-induced fall in MAP was significantly attenuated by indometacin. We concluded that VK1 induces a dose-dependent effect on blood pressure that consists of an acute increase followed by a more sustained decrease in MAP. The hypotension induced by VK1 involves the activation of the nitric oxide (NO) pathway and the release of vasodilator prostanoid(s).

Exact solution for the Bariev model with boundary fields

The Bariev model with open boundary conditions is introduced and analysed in detail in the framework of the Quantum Inverse Scattering Method. Two classes of independent boundary reflecting K-matrices leading to four different types of boundary fields are obtained by solving the reflection equations. The models are exactly solved by means of the algebraic nested Bethe ansatz method and the four sets or Bethe ansatz equations as well as their corresponding energy expressions are derived. (C) 2001 Elsevier Science B.V. All rights reserved.

Orbitally induced oscillations in the East Antarctic ice sheet at the Oligocene/Miocene boundary

Between 34 and 15 million years (Myr) ago, when planetary temperatures were 3-4 degreesC warmer than at present and atmospheric CO2 concentrations were twice as high as today(1), the Antarctic ice sheets may have been unstable(2-7). Oxygen isotope records from deep-sea sediment cores suggest that during this time fluctuations in global temperatures and high-latitude continental ice volumes were influenced by orbital cycles(8-10). But it has hitherto not been possible to calibrate the inferred changes in ice volume with direct evidence for oscillations of the Antarctic ice sheets(11). Here we present sediment data from shallow marine cores in the western Ross Sea that exhibit well dated cyclic variations, and which link the extent of the East Antarctic ice sheet directly to orbital cycles during the Oligocene/Miocene transition (24.1-23.7 Myr ago). Three rapidly deposited glaci-marine sequences are constrained to a period of less than 450 kyr by our age model, suggesting that orbital influences at the frequencies of obliquity (40 kyr) and eccentricity (125 kyr) controlled the oscillations of the ice margin at that time. An erosional hiatus covering 250 kyr provides direct evidence for a major episode of global cooling and ice-sheet expansion about 23.7 Myr ago, which had previously been inferred from oxygen isotope data (Mil event(5)).

Relationship of glycosaminoglycan and matrix changes to vascular smooth cell phenotype modulation in rabbit arteries after acute injury

Purpose: The phenotype of vascular smooth muscle cells (SMCs) is altered in several arterial pathologies, including the neointima formed after acute arterial injury. This study examined the time course of this phenotypic change in relation to changes in the amount and distribution of matrix glycosaminoglycans. Methods: The immunochemical staining of heparan sulphates (HS) and chondroitin sulphates (CS) in the extracellular matrix of the arterial wall was examined at early points after balloon catheter injury of the rabbit carotid artery. SMC phenotype was assessed by means of ultrastructural morphometry of the cytoplasmic volume fraction of myofilaments. The proportions of cell and matrix components in the media were analyzed with similar morphometric techniques. Results: HS and CS were shown in close association with SMCs of the uninjured arterial media as well as being more widespread within the matrix. Within 6 hours after arterial injury, there was loss of the regular pericellular distribution of both HS and CS, which was associated with a significant expansion in the extracellular space. This preceded the change in ultrastructural phenotype of the SMCs. The glycosaminoglycan loss was most exaggerated at 4 days, after which time the HS and CS reappeared around the medial SMCs. SMCs of the recovering media were able to rapidly replace their glycosaminoglycans, whereas SMCs of the developing neointima failed to produce HS as readily as they produced CS. Conclusions: These studies indicate that changes in glycosaminoglycans of the extracellular matrix precede changes in SMC phenotype after acute arterial injury. In the recovering arterial media, SMCs replace their matrix glycosaminoglycans rapidly, whereas the newly established neointima fails to produce similar amounts of heparan sulphates.

Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

Application of Petrov-Galerkin methods to transient boundary value problems in chemical engineering: adsorption with steep gradients in bidisperse solids

Petrov-Galerkin methods are known to be versatile techniques for the solution of a wide variety of convection-dispersion transport problems, including those involving steep gradients. but have hitherto received little attention by chemical engineers. We illustrate the technique by means of the well-known problem of simultaneous diffusion and adsorption in a spherical sorbent pellet comprised of spherical, non-overlapping microparticles of uniform size and investigate the uptake dynamics. Solutions to adsorption problems exhibit steep gradients when macropore diffusion controls or micropore diffusion controls, and the application of classical numerical methods to such problems can present difficulties. In this paper, a semi-discrete Petrov-Galerkin finite element method for numerically solving adsorption problems with steep gradients in bidisperse solids is presented. The numerical solution was found to match the analytical solution when the adsorption isotherm is linear and the diffusivities are constant. Computed results for the Langmuir isotherm and non-constant diffusivity in microparticle are numerically evaluated for comparison with results of a fitted-mesh collocation method, which was proposed by Liu and Bhatia (Comput. Chem. Engng. 23 (1999) 933-943). The new method is simple, highly efficient, and well-suited to a variety of adsorption and desorption problems involving steep gradients. (C) 2001 Elsevier Science Ltd. All rights reserved.