Stochastic perturbations to dynamical systems: a response theory approach

Using the formalism of the Ruelle response theory, we study how the invariant measure of an Axiom A dynamical system changes as a result of adding noise, and describe how the stochastic perturbation can be used to explore the properties of the underlying deterministic dynamics. We first find the expression for the change in the expectation value of a general observable when a white noise forcing is introduced in the system, both in the additive and in the multiplicative case. We also show that the difference between the expectation value of the power spectrum of an observable in the stochastically perturbed case and of the same observable in the unperturbed case is equal to the variance of the noise times the square of the modulus of the linear susceptibility describing the frequency-dependent response of the system to perturbations with the same spatial patterns as the considered stochastic forcing. This provides a conceptual bridge between the change in the fluctuation properties of the system due to the presence of noise and the response of the unperturbed system to deterministic forcings. Using Kramers-Kronig theory, it is then possible to derive the real and imaginary part of the susceptibility and thus deduce the Green function of the system for any desired observable. We then extend our results to rather general patterns of random forcing, from the case of several white noise forcings, to noise terms with memory, up to the case of a space-time random field. Explicit formulas are provided for each relevant case analysed. As a general result, we find, using an argument of positive-definiteness, that the power spectrum of the stochastically perturbed system is larger at all frequencies than the power spectrum of the unperturbed system. We provide an example of application of our results by considering the spatially extended chaotic Lorenz 96 model. These results clarify the property of stochastic stability of SRB measures in Axiom A flows, provide tools for analysing stochastic parameterisations and related closure ansatz to be implemented in modelling studies, and introduce new ways to study the response of a system to external perturbations. Taking into account the chaotic hypothesis, we expect that our results have practical relevance for a more general class of system than those belonging to Axiom A.

Bistable systems with stochastic noise: virtues and limits of effective one-dimensional Langevin equations

The understanding of the statistical properties and of the dynamics of multistable systems is gaining more and more importance in a vast variety of scientific fields. This is especially relevant for the investigation of the tipping points of complex systems. Sometimes, in order to understand the time series of given observables exhibiting bimodal distributions, simple one-dimensional Langevin models are fitted to reproduce the observed statistical properties, and used to investing-ate the projected dynamics of the observable. This is of great relevance for studying potential catastrophic changes in the properties of the underlying system or resonant behaviours like those related to stochastic resonance-like mechanisms. In this paper, we propose a framework for encasing this kind of studies, using simple box models of the oceanic circulation and choosing as observable the strength of the thermohaline circulation. We study the statistical properties of the transitions between the two modes of operation of the thermohaline circulation under symmetric boundary forcings and test their agreement with simplified one-dimensional phenomenological theories. We extend our analysis to include stochastic resonance-like amplification processes. We conclude that fitted one-dimensional Langevin models, when closely scrutinised, may result to be more ad-hoc than they seem, lacking robustness and/or well-posedness. They should be treated with care, more as an empiric descriptive tool than as methodology with predictive power.

Attention through self-synchronisation in the spiking neuron stochastic diffusion network

The paper discusses ensemble behaviour in the Spiking Neuron Stochastic Diffusion Network, SNSDN, a novel network exploring biologically plausible information processing based on higher order temporal coding. SNSDN was proposed as an alternative solution to the binding problem [1]. SNSDN operation resembles Stochastic Diffusin on Search, SDS, a non-deterministic search algorithm able to rapidly locate the best instantiation of a target pattern within a noisy search space ([3], [5]). In SNSDN, relevant information is encoded in the length of interspike intervals. Although every neuron operates in its own time, ‘attention’ to a pattern in the search space results in self-synchronised activity of a large population of neurons. When multiple patterns are present in the search space, ‘switching of at- tention’ results in a change of the synchronous activity. The qualitative effect of attention on the synchronicity of spiking behaviour in both time and frequency domain will be discussed.

Time complexity analysis of the stochastic diffusion search

The Stochastic Diffusion Search algorithm -an integral part of Stochastic Search Networks is investigated. Stochastic Diffusion Search is an alternative solution for invariant pattern recognition and focus of attention. It has been shown that the algorithm can be modelled as an ergodic, finite state Markov Chain under some non-restrictive assumptions. Sub-linear time complexity for some settings of parameters has been formulated and proved. Some properties of the algorithm are then characterised and numerical examples illustrating some features of the algorithm are presented.

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.

Cryopreservation of winter-dormant apple: III Bud water status and survival after cooling to -30°c and during recovery from cryopreservation

Abstract In a continuing study to improve the efficiency of dormant bud cryopreservation for tissues hardened in maritime climates, the water status of dormant buds was monitored between -4°C and recovery from liquid nitrogen (LN). Measurement of water content, simple thermal analysis and differential scanning calorimetry were employed. Buds did not lose water during cooling to, or holding at -30°C indicating that cryodehydration and/or other adaptive responses contributed during this essential step. A bud exotherm that was an artefact of warming was detected due to necessary handling at -4°C before cooling to -30°C. There were no significant differences between cultivars with respect to water status at -30°C or immediately upon rewarming from LN despite significant differences in post-LN survival. Buds rehydrated in 5 days, but up to 14 days may be needed for recovery for some cultivars. In some instances buds could be grafted without rehydration, taking up water across the early graft union.

Cryopreservation of winter-dormant apple buds: II - tissue water status after desiccation at -4°c and before further cooling

Abstract The established protocol for the cryopreservation of winter-dormant Malus buds requires that stem explants, containing a single, dormant bud are desiccated at -4°C, for up to 14 days, to reduce their water content to 25-30% of fresh weight. Using three apple cultivars, with known differences in response to cryopreservation, the pattern of evaporative water loss has been characterised, including early freezing events in the bud and cortical tissues that allow further desiccation by water migration to extracellular ice. There were no significant differences between cultivars in this respect or in the proportions of tissue water lost during the desiccation process. Differential Scanning Calorimetry (to -90°C) of intact buds indicated that bud tissues of the cultivar with the poorest response to cryopreservation had the highest residual water content at the end of the desiccation process and froze at the highest temperature Keywords: Malus, cryopreservation, dormant bud, dehydration

Alternatives to Sedum on green roofs: Can broad leaf perennial plants offer better ‘cooling service’?

Green roof plants alter the microclimate of building roofs and may improve roof insulation. They act by providing cooling by shading, but also through transpiration of water through their stomata. However, leaf surfaces can become warmer when plants close the stomata and decrease water loss in response to drying substrate (typically associated with green roofs during summers), also reducing transpirational cooling. By using a range of contrasting plant types (Sedum mix – an industry green roof ‘standard’, Stachys byzantina, Bergenia cordifolia and Hedera hibernica) we tested the hypothesis that plants differ in their ‘cooling potential’. We firstly examined how leaf morphology influenced leaf temperature and how drying substrate altered that response. Secondly, we investigated the relationship between leaf surface temperatures and the air temperatures immediately above the canopies (i.e. potential to provide aerial cooling). Finally we measured how the plant type influenced the substrate temperature below the canopy (i.e. potential for building cooling). In our experiments Stachys outperformed the other species in terms of leaf surface cooling (even in drying substrate, e.g. 5 oC cooler compared with Sedum), substrate cooling beneath its canopy (up to 12 oC) and even - during short intervals over hottest still periods - the air above the canopy (up to 1 oC, when soil moisture was not limited). We suggest that the choice of plant species on green roofs should not be entirely dictated by what survives on the shallow substrates of extensive systems, but consideration should be given to supporting those species providing the greatest eco-system service potential.

On the long time behavior of free stochastic Schrödinger evolutions

We discuss the time evolution of the wave function which is the solution of a stochastic Schrödinger equation describing the dynamics of a free quantum particle subject to spontaneous localizations in space. We prove global existence and uniqueness of solutions. We observe that there exist three time regimes: the collapse regime, the classical regime and the diffusive regime. Concerning the latter, we assert that the general solution converges almost surely to a diffusing Gaussian wave function having a finite spread both in position as well as in momentum. This paper corrects and completes earlier works on this issue.

On a stochastic Trotter formula with application to spontaneous localization models

We consider the relation between so called continuous localization models—i.e. non-linear stochastic Schrödinger evolutions—and the discrete GRW-model of wave function collapse. The former can be understood as scaling limit of the GRW process. The proof relies on a stochastic Trotter formula, which is of interest in its own right. Our Trotter formula also allows to complement results on existence theory of stochastic Schrödinger evolutions by Holevo and Mora/Rebolledo.

Analysis of Direct Outdoor Air Cooling Efficency for Combined Variable Air Volume Air-conditioning System in Stores in Cold Climates of China

Direct outdoor air cooling contributes a lot not only to the improvement of the indoor air quality but also to the energy saving. Its full use will reduce the water chiller’s running time especially in some stores where cooling load keeps much higher and longer than that in other buildings. A novel air-conditioning system named Combined Variable Air Volume system (CVAV), combining a normal AHU with a separate outdoor air supply system, was proposed firstly by the authors. The most attractive feature of the system is its full utilization of cooling capacity and freshness of outdoor air in the transition period of the year round. On the basis of the obtain of the dynamic cooling loads of the typical shopping malls in different four cities located in cold climates in China with the aid of DOE-2, the possibility of increasing the amount of outdoor air volume of CVAV system in the transition period instead of operating the water chillers was confirmed. Moreover, a new concept, Direct Outdoor Air Cooling Efficiency (DOACE), was defined as the ratio of cooling capacity of outdoor air to the water chiller, indicating the degree of outdoor air’s utilization. And the DOACE of the CVAV was calculated and compared with that of conventional all-air constant volume air-conditioning systems, the results showed that CVAV bear much more energy saving potential with the 10%~19% higher DOACE and it is a kind of energy efficient systems and can improve the indoor air quality as well.

Numerical study of a M-cycle cross-flow heat exchanger for indirect evaporative cooling

In this paper, numerical analyses of the thermal performance of an indirect evaporative air cooler incorporating a M-cycle cross-flow heat exchanger has been carried out. The numerical model was established from solving the coupled governing equations for heat and mass transfer between the product and working air, using the finite-element method. The model was developed using the EES (Engineering Equation Solver) environment and validated by published experimental data. Correlation between the cooling (wet-bulb) effectiveness, system COP and a number of air flow/exchanger parameters was developed. It is found that lower channel air velocity, lower inlet air relative humidity, and higher working-to-product air ratio yielded higher cooling effectiveness. The recommended average air velocities in dry and wet channels should not be greater than 1.77 m/s and 0.7 m/s, respectively. The optimum flow ratio of working-to-product air for this cooler is 50%. The channel geometric sizes, i.e. channel length and height, also impose significant impact to system performance. Longer channel length and smaller channel height contribute to increase of the system cooling effectiveness but lead to reduced system COP. The recommend channel height is 4 mm and the dimensionless channel length, i.e., ratio of the channel length to height, should be in the range 100 to 300. Numerical study results indicated that this new type of M-cycle heat and mass exchanger can achieve 16.7% higher cooling effectiveness compared with the conventional cross-flow heat and mass exchanger for the indirect evaporative cooler. The model of this kind is new and not yet reported in literatures. The results of the study help with design and performance analyses of such a new type of indirect evaporative air cooler, and in further, help increasing market rating of the technology within building air conditioning sector, which is currently dominated by the conventional compression refrigeration technology.