Direct Computation of Stochastic Flow in Reservoirs with Uncertain Parameters

A direct method is presented for determining the uncertainty in reservoir pressure, flow, and net present value (NPV) using the time-dependent, one phase, two- or three-dimensional equations of flow through a porous medium. The uncertainty in the solution is modelled as a probability distribution function and is computed from given statistical data for input parameters such as permeability. The method generates an expansion for the mean of the pressure about a deterministic solution to the system equations using a perturbation to the mean of the input parameters. Hierarchical equations that define approximations to the mean solution at each point and to the field covariance of the pressure are developed and solved numerically. The procedure is then used to find the statistics of the flow and the risked value of the field, defined by the NPV, for a given development scenario. This method involves only one (albeit complicated) solution of the equations and contrasts with the more usual Monte-Carlo approach where many such solutions are required. The procedure is applied easily to other physical systems modelled by linear or nonlinear partial differential equations with uncertain data.

Modelling the time-evolution of phytoplankton size spectra from satellite remote sensing

A dynamic size-structured model is developed for phytoplankton and nutrients in the oceanic mixed layer and applied to extract phytoplankton biomass at discrete size fractions from remotely sensed, ocean-colour data. General relationships between cell size and biophysical processes (such as sinking, grazing, and primary production) of phytoplankton were included in the model through a bottom–up approach. Time-dependent, mixed-layer depth was used as a forcing variable, and a sequential data-assimilation scheme was implemented to derive model trajectories. From a given time-series, the method produces estimates of size-structured biomass at every observation, so estimates seasonal succession of individual phytoplankton size, derived here from remote sensing for the first time. From these estimates, normalized phytoplankton biomass size spectra over a period of 9 years were calculated for one location in the North Atlantic. Further analysis demonstrated that strong relationships exist between the seasonal trends of the estimated size spectra and the mixed-layer depth, nutrient biomass, and total chlorophyll. The results contain useful information on the time-dependent biomass flux in the pelagic ecosystem.

Solvent-Dependent Formation of Os(0) Complexes by Electrochemical Reduction of [Os(CO)(2,2′-bipyridine)(L)Cl2]; L = Cl−,PrCN

Cyclic voltammetry and ultraviolet−visible/infrared (UV−vis/IR) spectroelectrochemistry were used to study the cathodic electrochemical behavior of the osmium complexes mer-[OsIII(CO) (bpy)Cl3] (bpy = 2,2′-bipyridine) and trans(Cl)-[OsII(CO) (PrCN)(bpy)Cl2] at variable temperature in different solvents (tetrahydrofuran (THF), butyronitrile (PrCN), acetonitrile (MeCN)) and electrolytes (Bu4NPF6, Bu4NCl). The precursors can be reduced to mer-[OsII(CO) (bpy•−)Cl3]2− and trans(Cl)-[OsII(CO)(PrCN) (bpy•−)Cl2]−, respectively, which react rapidly at room temperature, losing the chloride ligands and forming Os(0) species. mer-[OsIII(CO) (bpy)Cl3] is reduced in THF to give ultimately an Os−Os-bonded polymer, probably [Os0(CO) (THF)-(bpy)]n, whereas in PrCN the well-soluble, probably mononuclear [Os0(CO) (PrCN)(bpy)], species is formed. The same products were observed for the 2 electron reduction of trans(Cl)-[OsII(CO)(PrCN) (bpy)Cl2] in both solvents. In MeCN, similar to THF, the[Os0(CO) (MeCN)(bpy)]n polymer is produced. It is noteworthy that the bpy ligand in mononuclear [Os0(CO) (PrCN)(bpy)] is reduced to the corresponding radical anion at a significantly less negative potential than it is in polymeric [Os0(CO) (THF)(bpy)]n: ΔE1/2 = 0.67 V. Major differences also exist in the IR spectra of the Os(0) species: the polymer shows a broad ν(CO) band at much smaller wavenumbers compared to the soluble Os(0) monomer that exhibits a characteristic ν(Pr-CN) band below 2200 cm−1 in addition to the intense and narrow ν(CO) absorption band. For the first time, in this work the M0-bpy(M = Ru, Os) mono- and dicarbonyl species soluble in PrCN have been formulated as a mononuclear complex. Density functional theory (DFT) and time-dependent-DFT calculations confirm the Os(0) oxidation state and suggest that [Os0(CO)(PrCN)(bpy)] is a square planar moiety. The reversible bpy-based reduction of [Os0(CO) (PrCN)(bpy)] triggers catalytic reduction of CO2 to CO and HCOO−.

A numerical model of the ionospheric signatures of time-varying magneticreconnection: I. ionospheric convection

This paper presents a numerical model for predicting the evolution of the pattern of ionospheric convection in response to general time-dependent magnetic reconnection at the dayside magnetopause and in the cross-tail current sheet of the geomagnetic tail. The model quantifies the concepts of ionospheric flow excitation by Cowley and Lockwood (1992), assuming a uniform spatial distribution of ionospheric conductivity. The model is demonstrated using an example in which travelling reconnection pulses commence near noon and then move across the dayside magnetopause towards both dawn and dusk. Two such pulses, 8 min apart, are used and each causes the reconnection to be active for 1 min at every MLT that they pass over. This example demonstrates how the convection response to a given change in the interplanetary magnetic field (via the reconnection rate) depends on the previous reconnection history. The causes of this effect are explained. The inherent assumptions and the potential applications of the model are discussed.

Modelling the high-latitude ionosphere for time-varying plasma convection

The paper discusses how variations in the pattern of convective plasma flows should beincluded in self-consistent time-dependent models of the coupled ionosphere-thermosphere system. The author shows how these variations depend upon the mechanism by which the solar wind flow excites the convection. The modelling of these effects is not just of relevance to the polar ionosphere. This is because the influence of convection is not confined to high latitudes: the resultant heating and composition changes in the thermosphere are communicated to lower latitudes by the winds which are also greatly modified by the plasma convection. These thermospheric changes alter the global distribution of plasma by modulatingthe rates of the chemical reactions which areresponsible for the loss of plasma. Hence the modelling of these high-latitude processes is of relevanceto the design and operation of HF communication, radar and navigation systems worldwide.

Melatonin and the time window for the expression of the alpha 8 nicotinic acetylcholine receptor in the membrane of chick retinal cells in culture

We have previously shown that melatonin influences the development of alpha 8 nicotinic acetylcholine receptor (nAChR) by measurement of the acetylcholine-induced increase in the extracellular acidification rate (ECAR) in chick retinal cell cultures. Cellular differentiation that takes place between DIV (days in vitro) 4 and DIV 5 yields cells expressing alpha 8 nAChR and results in a significant increase in the ECAR acetylcholine-induced. Blocking melatonin receptors with luzindole for 48 h suppresses the development of functional alpha 8 nAChR. Here we investigated the time window for the effect of melatonin on retinal cell development in culture, and whether this effect was dependent on an increase in the expression of alpha 8 nAChR. First, we confirmed that luzindole was inhibiting the effects of endogenous melatonin, since it increases 2-[(125)I] iodomelatonin (23 pM) binding sites density in a time-dependent manner. Then we observed that acute (15, 60 min, or 12 h) luzindole treatment did not impair acetylcholine-induced increase in the ECAR mediated by activation of alpha 8 nAChR at DIV 5, while chronic treatment (from DIV 3 or DIV 4 till DIV 5, or DIV 3.5 till DIV 4.5) led to a time-dependent reduction of the increase in the acetylcholine-induced ECAR. The binding parameters for [(125)I]-alpha-bungarotoxin (10 nM) sites in membrane were unaffected by melatonin suppression that started at DIV 3. Thus, melatonin surges in the time window that occurs at the final stages of chick retinal cell differentiation in culture is essential for development of the cells expressing alpha 8 nAChR subtype in full functional form. (C) 2010 ISDN. Published by Elsevier Ltd. All rights reserved.

Hydrogen peroxide is a novel mediator of inflammatory hyperalgesia, acting via transient receptor potential vanilloid 1-dependent and independent mechanisms

Inflammatory diseases associated with pain are often difficult to treat in the clinic due to insufficient understanding of the nociceptive pathways involved. Recently, there has been considerable interest in the role of reactive oxygen species (ROS) in inflammatory disease, but little is known of the role of hydrogen peroxide (H(2)O(2)) in hyperalgesia. In the present study, intraplantar injection of H(2)O(2)-induced a significant dose- and time-dependent mechanical and thermal hyperalgesia in the mouse hind paw, with increased c-fos activity observed in the dorsal horn of the spinal cord. H(2)O(2) also induced significant nociceptive behavior Such as increased paw licking and decreased body liftings. H(2)O(2) levels were significantly raised in the carrageenan-induced hind paw inflammation model, showing that this ROS is produced endogenously in a model of inflammation. Moreover, superoxide dismutase and catalase significantly reduced carrageenan-induced mechanical and thermal hyperalgesia, providing evidence of a functionally significant endogenous role. Thermal, but not mechanical, hyperalgesia in response to H(2)O(2) (i.pl.) Was longer lasting in TRPV1 wild type mice compared to TRPV1 knockouts. It is unlikely that downstream lipid peroxidation was increased by H(2)O(2). In conclusion, we demonstrate a notable effect of H(2)O(2) in mediating inflammatory hyperalgesia, thus highlighting H(2)O(2) removal as a novel therapeutic target for anti-hyperalgesic drugs in the clinic. (C) 2008 International Association for the Study of Pain. Published by Elsevier B.V. All rights reserved.

A non-autonomous strongly damped wave equation: Existence and continuity of the pullback attractor

In this paper we consider the strongly damped wave equation with time-dependent terms u(tt) - Delta u - gamma(t)Delta u(t) + beta(epsilon)(t)u(t) = f(u), in a bounded domain Omega subset of R(n), under some restrictions on beta(epsilon)(t), gamma(t) and growth restrictions on the nonlinear term f. The function beta(epsilon)(t) depends on a parameter epsilon, beta(epsilon)(t) -> 0. We will prove, under suitable assumptions, local and global well-posedness (using the uniform sectorial operators theory), the existence and regularity of pullback attractors {A(epsilon)(t) : t is an element of R}, uniform bounds for these pullback attractors, characterization of these pullback attractors and their upper and lower semicontinuity at epsilon = 0. (C) 2010 Elsevier Ltd. All rights reserved.

Lower semicontinuity of attractors for non-autonomous dynamical systems

This paper is concerned with the lower semicontinuity of attractors for semilinear non-autonomous differential equations in Banach spaces. We require the unperturbed attractor to be given as the union of unstable manifolds of time-dependent hyperbolic solutions, generalizing previous results valid only for gradient-like systems in which the hyperbolic solutions are equilibria. The tools employed are a study of the continuity of the local unstable manifolds of the hyperbolic solutions and results on the continuity of the exponential dichotomy of the linearization around each of these solutions.

Dynamics of a class of ODEs more general than almost periodic

A temporally global solution, if it exists, of a nonautonomous ordinary differential equation need not be periodic, almost periodic or almost automorphic when the forcing term is periodic, almost periodic or almost automorphic, respectively. An alternative class of functions extending periodic and almost periodic functions which has the property that a bounded temporally global solution solution of a nonautonomous ordinary differential equation belongs to this class when the forcing term does is introduced here. Specifically, the class of functions consists of uniformly continuous functions, defined on the real line and taking values in a Banach space, which have pre-compact ranges. Besides periodic and almost periodic functions, this class also includes many nonrecurrent functions. Assuming a hyperbolic structure for the unperturbed linear equation and certain properties for the linear and nonlinear parts, the existence of a special bounded entire solution, as well the existence of stable and unstable manifolds of this solution are established. Moreover, it is shown that this solution and these manifolds inherit the temporal behaviour of the vector field equation. In the stable case it is shown that this special solution is the pullback attractor of the system. A class of infinite dimensional examples involving a linear operator consisting of a time independent part which generates a C(0)-semigroup plus a small time dependent part is presented and applied to systems of coupled heat and beam equations. (C) 2010 Elsevier Ltd. All rights reserved.

Susceptible-infected-recovered and susceptible-exposed-infected models

Two stochastic epidemic lattice models, the susceptible-infected-recovered and the susceptible-exposed-infected models, are studied on a Cayley tree of coordination number k. The spreading of the disease in the former is found to occur when the infection probability b is larger than b(c) = k/2(k - 1). In the latter, which is equivalent to a dynamic site percolation model, the spreading occurs when the infection probability p is greater than p(c) = 1/(k - 1). We set up and solve the time evolution equations for both models and determine the final and time-dependent properties, including the epidemic curve. We show that the two models are closely related by revealing that their relevant properties are exactly mapped into each other when p = b/[k - (k - 1) b]. These include the cluster size distribution and the density of individuals of each type, quantities that have been determined in closed forms.

A general Trotter-Kato formula for a class of evolution operators

In this article we prove new results concerning the existence and various properties of an evolution system U(A+B)(t, s)0 <= s <= t <= T generated by the sum -(A(t) + B(t)) of two linear, time-dependent, and generally unbounded operators defined on time-dependent domains in a complex and separable Banach space B. In particular, writing L(B) for the algebra of all linear bounded operators on B, we can express U(A+B)(t, s)0 <= s <= t <= T as the strong limit in C(8) of a product of the holomorphic contraction semigroups generated by -A (t) and - B(t), respectively, thereby proving a product formula of the Trotter-Kato type under very general conditions which allow the domain D(A(t) + B(t)) to evolve with time provided there exists a fixed set D subset of boolean AND(t is an element of)[0,T] D(A(t) + B(t)) everywhere dense in B. We obtain a special case of our formula when B(t) = 0, which, in effect, allows us to reconstruct U(A)(t, s)0 <=(s)<=(t)<=(T) very simply in terms of the semigroup generated by -A(t). We then illustrate our results by considering various examples of nonautonomous parabolic initial-boundary value problems, including one related to the theory of timedependent singular perturbations of self-adjoint operators. We finally mention what we think remains an open problem for the corresponding equations of Schrodinger type in quantum mechanics.

Velocity-dependent KAM islands

We revisit the non-dissipative time-dependent annular billiard and we consider the chaotic dynamics in two planes of conjugate variables in order to describe the behavior of the growth, or saturation, of the mean velocity of an ensemble of particles. We observed that the changes in the 4-d phase space occur without changing any parameter. They occur depending on where the initial conditions start. The emerging KAM islands interfere in the behavior of the particle dynamics especially in the Fermi acceleration mechanism. We show that Fermi acceleration can be suppressed, without dissipation, even considering the non-dissipative energy context. (C) 2011 Published by Elsevier Ltd.

Self-dual hopfions

We construct static and time-dependent exact soliton solutions with nontrivial Hopf topological charge for a field theory in 3 + 1 dimensions with the target space being the two dimensional sphere S(2). The model considered is a reduction of the so-called extended Skyrme-Faddeev theory by the removal of the quadratic term in derivatives of the fields. The solutions are constructed using an ansatz based on the conformal and target space symmetries. The solutions are said self-dual because they solve first order differential equations which together with some conditions on the coupling constants, imply the second order equations of motion. The solutions belong to a sub-sector of the theory with an infinite number of local conserved currents. The equation for the profile function of the ansatz corresponds to the Bogomolny equation for the sine-Gordon model.

Static Hopfions in the extended Skyrme-Faddeev model

We construct static soliton solutions with non-zero Hopf topological charges to a theory which is an extension of the Skyrme-Faddeev model by the addition of a further quartic term in derivatives. We use an axially symmetric ansatz based on toroidal coordinates, and solve the resulting two coupled non-linear partial differential equations in two variables by a successive over-relaxation (SOR) method. We construct numerical solutions with Hopf charge up to four, and calculate their analytical behavior in some limiting cases. The solutions present an interesting behavior under the changes of a special combination of the coupling constants of the quartic terms. Their energies and sizes tend to zero as that combination approaches a particular special value. We calculate the equivalent of the Vakulenko and Kapitanskii energy bound for the theory and find that it vanishes at that same special value of the coupling constants. In addition, the model presents an integrable sector with an in finite number of local conserved currents which apparently are not related to symmetries of the action. In the intersection of those two special sectors the theory possesses exact vortex solutions (static and time dependent) which were constructed in a previous paper by one of the authors. It is believed that such model describes some aspects of the low energy limit of the pure SU(2) Yang-Mills theory, and our results may be important in identifying important structures in that strong coupling regime.