WLS method for parameter estimation in water distribution networks

The weighted-least-squares method based on the Gauss-Newton minimization technique is used for parameter estimation in water distribution networks. The parameters considered are: element resistances (single and/or group resistances, Hazen-Williams coefficients, pump specifications) and consumptions (for single or multiple loading conditions). The measurements considered are: nodal pressure heads, pipe flows, head loss in pipes, and consumptions/inflows. An important feature of the study is a detailed consideration of the influence of different choice of weights on parameter estimation, for error-free data, noisy data, and noisy data which include bad data. The method is applied to three different networks including a real-life problem.

The determination of spatial pattern in Dictyostelium discoideum

Free-living amoebae of the cellular slime mould Dictyostelium discoideum aggregate when starved and give rise to a long and thin multicellular structure, the slug. The slug resembles a metazoan embryo, and as with other embryos it is possible to specify a fate map. In the case of Dictyostelium discoideum the map is especially simple: cells in the anterior fifth of the slug die and form a stalk while the majority of those in the posterior differentiate into spores. The genesis of this anterior-posterior distinction is the subject of our review. In particular, we ask: what are the relative roles of individual pre-aggregative predispositions and post-aggregative position in determining cell fate? We review the literature on the subject and conclude that both factors are important. Variations in nutritional status, or in cell cycle phase at starvation, can bias the probability that an amoeba differentiates into a stalk cell or a spore. On the other hand, isolates, or slug fragments, consisting of only prestalk cells or only prespore cells can regulate so as to result in a normal range of both cell types. We identify three levels of control, each being responsible for guiding patterning in normal development: (i) 'coin tossing', whereby a cell autonomously exhibits a preference for developing along either the stalk or the spore pathway with relative probabilities that can be influenced by the environment; (ii) 'chemical kinetics', whereby prestalk and prespore cells originate from undifferentiated amoebae on a probabilistic basis but, having originated, interact (e.g. via positive and negative feedbacks), and the interaction influences the possibility of conversion of one cell type into the other, and (iii) 'positional information', in which the spatial distribution of morphogens in the slug influences the pathway of differentiation. In the case of possibilities (i) and (ii), sorting out of like cell types leads to the final spatial pattern. In the case of possibility (iii), the pattern arises in situ.

Multi-installment load distribution in tree networks with delays

This paper presents a new strategy for load distribution in a single-level tree network equipped with or without front-ends. The load is distributed in more than one installment in an optimal manner to minimize the processing time. This is a deviation and an improvement over earlier studies in which the load distribution is done in only one installment. Recursive equations for the general case, and their closed form solutions for a special case in which the network has identical processors and identical links, are derived. An asymptotic analysis of the network performance with respect to the number of processors and the number of installments is carried out. Discussions of the results in terms of some practical issues like the tradeoff relationship between the number of processors and the number of installments are also presented.

Autocalibration of dynamic bacterial growth model for water distribution system using GA

The presence of residual chlorine and organic matter govern the bacterial regrowth within a water distribution system. The bacterial growth model is essential to predict the spatial and temporal variation of all these substances throughout the system. The parameters governing the bacterial growth and biodegradable dissolved organic carbon (BDOC) utilization are difficult to determine by experimentation. In the present study, the estimation of these parameters is addressed by using simulation-optimization procedure. The optimal solution by genetic algorithm (GA) has indicated that the proper combination of parameter values are significant rather than correct individual values. The applicability of the model is illustrated using synthetic data generated by introducing noise in to the error-free measurements. The GA was found to be a potential tool in estimating the parameters controlling the bacterial growth and BDOC utilization. Further, the GA was also used for evaluating the sensitivity issues relating parameter values and objective function. It was observed that mu and k(cl) are more significant and dominating compared to the other parameters. But the magnitude of the parameters is also an important issue in deciding the dominance of a particular parameter. GA is found to be a useful tool in autocalibration of bacterial growth model and a sensitivity study of parameters.

Velocity distribution function for a dilute granular material in shear flow

The velocity distribution function for the steady shear flow of disks (in two dimensions) and spheres (in three dimensions) in a channel is determined in the limit where the frequency of particle-wall collisions is large compared to particle-particle collisions. An asymptotic analysis is used in the small parameter epsilon, which is naL in two dimensions and na(2)L in three dimensions, where; n is the number density of particles (per unit area in two dimensions and per unit volume in three dimensions), L is the separation of the walls of the channel and a is the particle diameter. The particle-wall collisions are inelastic, and are described by simple relations which involve coefficients of restitution e(t) and e(n) in the tangential and normal directions, and both elastic and inelastic binary collisions between particles are considered. In the absence of binary collisions between particles, it is found that the particle velocities converge to two constant values (u(x), u(y)) = (+/-V, O) after repeated collisions with the wall, where u(x) and u(y) are the velocities tangential and normal to the wall, V = (1 - e(t))V-w/(1 + e(t)), and V-w and -V-w, are the tangential velocities of the walls of the channel. The effect of binary collisions is included using a self-consistent calculation, and the distribution function is determined using the condition that the net collisional flux of particles at any point in velocity space is zero at steady state. Certain approximations are made regarding the velocities of particles undergoing binary collisions :in order to obtain analytical results for the distribution function, and these approximations are justified analytically by showing that the error incurred decreases proportional to epsilon(1/2) in the limit epsilon --> 0. A numerical calculation of the mean square of the difference between the exact flux and the approximate flux confirms that the error decreases proportional to epsilon(1/2) in the limit epsilon --> 0. The moments of the velocity distribution function are evaluated, and it is found that [u(x)(2)] --> V-2, [u(y)(2)] similar to V-2 epsilon and -[u(x)u(y)] similar to V-2 epsilon log(epsilon(-1)) in the limit epsilon --> 0. It is found that the distribution function and the scaling laws for the velocity moments are similar for both two- and three-dimensional systems.

Residence time distribution for a class of Gaussian Markov processes

We study the distribution of residence time or equivalently that of "mean magnetization" for a family of Gaussian Markov processes indexed by a positive parameter alpha. The persistence exponent for these processes is simply given by theta=alpha but the residence time distribution is nontrivial. The shape of this distribution undergoes a qualitative change as theta increases, indicating a sharp change in the ergodic properties of the process. We develop two alternate methods to calculate exactly but recursively the moments of the distribution for arbitrary alpha. For some special values of alpha, we obtain closed form expressions of the distribution function. [S1063-651X(99)03306-1].

Stochastic dynamics modeling of the protein sequence length distribution in genomes: implications for microbial evolution

In this paper, we report an analysis of the protein sequence length distribution for 13 bacteria, four archaea and one eukaryote whose genomes have been completely sequenced, The frequency distribution of protein sequence length for all the 18 organisms are remarkably similar, independent of genome size and can be described in terms of a lognormal probability distribution function. A simple stochastic model based on multiplicative processes has been proposed to explain the sequence length distribution. The stochastic model supports the random-origin hypothesis of protein sequences in genomes. Distributions of large proteins deviate from the overall lognormal behavior. Their cumulative distribution follows a power-law analogous to Pareto's law used to describe the income distribution of the wealthy. The protein sequence length distribution in genomes of organisms has important implications for microbial evolution and applications. (C) 1999 Elsevier Science B.V. All rights reserved.

Velocity distribution for a dilute vibrated granular material

The velocity distribution for a vibrated granular material is determined in the dilute limit where the frequency of particle collisions with the vibrating surface is large compared to the frequency of binary collisions. The particle motion is driven by the source of energy due to particle collisions with the vibrating surface, and two dissipation mechanisms-inelastic collisions and air drag-are considered. In the latter case, a general form for the drag force is assumed. First, the distribution function for the vertical velocity for a single particle colliding with a vibrating surface is determined in the limit where the dissipation during a collision due to inelasticity or between successive collisions due to drag is small compared to the energy of a particle. In addition, two types of amplitude functions for the velocity of the surface, symmetric and asymmetric about zero velocity, are considered. In all cases, differential equations for the distribution of velocities at the vibrating surface are obtained using a flux balance condition in velocity space, and these are solved to determine the distribution function. It is found that the distribution function is a Gaussian distribution when the dissipation is due to inelastic collisions and the amplitude function is symmetric, and the mean square velocity scales as [[U-2](s)/(1 - e(2))], where [U-2](s) is the mean square velocity of the vibrating surface and e is the coefficient of restitution. The distribution function is very different from a Gaussian when the dissipation is due to air drag and the amplitude function is symmetric, and the mean square velocity scales as ([U-2](s)g/mu(m))(1/(m+2)) when the acceleration due to the fluid drag is -mu(m)u(y)\u(y)\(m-1), where g is the acceleration due to gravity. For an asymmetric amplitude function, the distribution function at the vibrating surface is found to be sharply peaked around [+/-2[U](s)/(1-e)] when the dissipation is due to inelastic collisions, and around +/-[(m +2)[U](s)g/mu(m)](1/(m+1)) when the dissipation is due to fluid drag, where [U](s) is the mean velocity of the surface. The distribution functions are compared with numerical simulations of a particle colliding with a vibrating surface, and excellent agreement is found with no adjustable parameters. The distribution function for a two-dimensional vibrated granular material that includes the first effect of binary collisions is determined for the system with dissipation due to inelastic collisions and the amplitude function for the velocity of the vibrating surface is symmetric in the limit delta(I)=(2nr)/(1 - e)much less than 1. Here, n is the number of particles per unit width and r is the particle radius. In this Limit, an asymptotic analysis is used about the Limit where there are no binary collisions. It is found that the distribution function has a power-law divergence proportional to \u(x)\((c delta l-1)) in the limit u(x)-->0, where u(x) is the horizontal velocity. The constant c and the moments of the distribution function are evaluated from the conservation equation in velocity space. It is found that the mean square velocity in the horizontal direction scales as O(delta(I)T), and the nontrivial third moments of the velocity distribution scale as O(delta(I)epsilon(I)T(3/2)) where epsilon(I) = (1 - e)(1/2). Here, T = [2[U2](s)/(1 - e)] is the mean square velocity of the particles.

Accelerated weathering of distribution class polymeric insulators

This paper deals with the long-term accelerated weathering of 11 kV polymeric insulators for 25000 h. Polymeric insulators were continuously subjected to accelerated weathering in a specially designed multistress-aging chamber under UV radiation, temperature and electric stress. Chemical, physical and electrical changes due to degradation have been assessed using various techniques. Some of the interesting results observed indicate that there is a significant reduction in the content of low molecular weight molecules, hydrophobicity was dynamic in nature and there is a significant increase in the surface roughness and oxidation levels with respect to the duration of the weathering.

Continuous distribution kinetics for the thermal degradation of LDPE in solution

Polymer degradation in solution has several advantages over melt pyrolysis, The degradation of low-density polyethylene (LDPE) occurs at much lower temperatures in solution (280-360degreesC) than in conventional melt pyrolysis (400-450degreesC). The thermal degradation kinetics of LDPE in solution was investigated in this work. LDPE was dissolved in liquid paraffin and degraded for 3 h at various temperatures (280-360degreesC). Samples were taken at specific times and analyzed with high-pressure liquid chromatography/gel permeation chromatography for the molecular weight distribution (MWD), The time evolution of the MWD was modeled with continuous distribution kinetics. Data indicated that LDPE followed random-chain-scission degradation. The rapid initial drop in molecular weight, observed up to 45 min, was attributed to the presence of weak links in the polymer. The rate coefficients for the breakage of weak and strong links were determined, and the corresponding average activation energies were calculated to be 88 and 24 kJ/mol, respectively. (C) 2002 John Wiley Sons, Inc.

Woody plant seedling distribution under invasive Lantana camara thickets in a dry-forest plot in Mudumalai, southern India

Lantana camara, a shrub of Central and South American origin, has become invasive across dry forests worldwide. The effect of the thicket-forming habit of L. camara as a dispersal and recruitment barrier in a community of native woody seedlings was examined in a 50-ha permanent plot located in the seasonally dry forest of Mudumalai, southern India. Sixty 100-m(2) plots were enumerated for native woody seedlings between 10-100 cm in height. Of these, 30 plots had no L. camara thickets, while the other 30 had dense thickets. The frequency of occurrence and abundance of seedlings were modelled as a function of dispersal mode (mammal, bird or mechanical) and affinities to forest habitats (dry forest, moist forest or ubiquitous) as well as presence or absence of dense L. camara thickets. Furthermore, frequency of occurrence and abundance of individual species were also compared between thickets and no L. camara. At the community level, L. camara density, dispersal mode and forest habitat affinities of species determined both frequency of occurrence and abundance of seedlings, with the abundance of dry-forest mammal-dispersed species and ubiquitous mechanically dispersed species being significantly lower under L. camara thickets. Phyllanthus emblica and Kydia calycina were found to be significantly less abundant under L. camara, whereas most other species were not affected by the presence of thickets. It was inferred that, by affecting the establishment of native tree seedlings, L. camara thickets could eventually alter the community composition of such forests.

Diperoxovanadate can substitute for H(2)O(2) at much lower concentration in inducing features of premature cellular senescence in mouse fibroblasts (NIH3T3)

Stress induced premature senescence (SIPS) in mammalian cells is an accelerated ageing response and experimentally obtained on treatment of cells with high concentrations of H(2)O(2), albeit at sub-lethal doses, because H(2)O(2) gets depleted by abundant cellular catalase. In the present study diperoxovanadate (DPV) was used as it is known to be stable at physiological pH, to be catalase-resistant and to substitute for H(2)O(2) in its activities at concentrations order of magnitudes lower. On treating NIH3T3 cells with DPV, SIPS-like morphology was observed along with an immediate response of rounding of the cells by disruption of actin cytoskeleton and transient G2/M arrest. DPV could bring about growth arrest and senescence associated features at 25 mu M dose, which were not seen with similar doses of either H(2)O(2) or vanadate. A minimal dose of 150 mu M of H(2)O(2) was required to induce similar affects as 25 mu M DPV. Increase in senescent associated markers such as p21, HMGA2 and PAI-1 was more prominent in DPV treated cells compared to similar dose of H(2)O(2). DPV-treated cells showed marked relocalization of Cyclin D1 from nucleus to cytoplasm. These results indicate that DPV, stable inorganic peroxide, is more efficient in inducing SIPS at lower concentrations compared to H(2)O(2). (C) 2011 Elsevier Ireland Ltd. All rights reserved.

Transformation in heterokaryons of Neurospora crassa is nuclear rather than cellular phenomenon

In Neurospora crassa, multinucleate macroconidia are used for genetic transformation. The barrier for such a transformation can be either at the cell membrane level or at the nuclear membrane level. For assessment of these possibilities, a forced heterokaryon (containing two genetically marked nuclei and auxotrophic for histidine) of Neurospora crassa was transformed with a plasmid containing his-3(+) gene. The transformants, which could grow without histidine supplementation, were then resolved into component homokaryons to determine into which nucleus or nuclei the plasmid had entered. Our results suggest that the barrier for transformation in Neurospora crassa is at the nuclear level, not at the cell membrane level. In a heterokaryon containing two genetically distinct nuclei, plasmid DNA integrated into only one of the nuclear types at any instance, but never into both nuclear types. Thus, in Neurospora crassa, the competent nucleus is essential for the transformation event to take place, and at a given time only one type of nucleus is competent to take up the exogenous DNA. Genomic Southern analysis showed that the transformants harbor both ectopic and homologous integrations of the plasmid DNA. The type and number of integrations were reflected at the post-translational level, since the specific activity of histidinol dehydrogenase (the translation product of his-3+ gene) was variable among several transformants and always less than the level of the wild type.