Parameter sensitivity and predictive uncertainty in a new water quality model, Q(2)

A new dynamic model of water quality, Q(2), has recently been developed, capable of simulating large branched river systems. This paper describes the application of a generalized sensitivity analysis (GSA) to Q(2) for single reaches of the River Thames in southern England. Focusing on the simulation of dissolved oxygen (DO) (since this may be regarded as a proxy for the overall health of a river); the GSA is used to identify key parameters controlling model behavior and provide a probabilistic procedure for model calibration. It is shown that, in the River Thames at least, it is more important to obtain high quality forcing functions than to obtain improved parameter estimates once approximate values have been estimated. Furthermore, there is a need to ensure reasonable simulation of a range of water quality determinands, since a focus only on DO increases predictive uncertainty in the DO simulations. The Q(2) model has been applied here to the River Thames, but it has a broad utility for evaluating other systems in Europe and around the world.

Identifying the critical parameters of a cyanobacterial growth and movement model by using generalised sensitivity analysis

Bloom-forming and toxin-producing cyanobacteria remain a persistent nuisance across the world. Modelling of cyanobacteria in freshwaters is an important tool for understanding their population dynamics and predicting the location and timing of the bloom events in lakes and rivers. A new deterministic-mathematical model was developed, which simulates the growth and movement of cyanobacterial blooms in river systems. The model focuses on the mathematical description of the bloom formation, vertical migration and lateral transport of colonies within river environments by taking into account the major factors that affect the cyanobacterial bloom formation in rivers including, light, nutrients and temperature. A technique called generalised sensitivity analysis was applied to the model to identify the critical parameter uncertainties in the model and investigates the interaction between the chosen parameters of the model. The result of the analysis suggested that 8 out of 12 parameters were significant in obtaining the observed cyanobacterial behaviour in a simulation. It was found that there was a high degree of correlation between the half-saturation rate constants used in the model.

On the time for fold surfaces to order during the crystallization of polyethylene from the melt and its dependence on molecular parameters

New experiments underpin the interpretation of the basic division in crystallization behaviour of polyethylene in terms of whether or not there is time for the fold surface to order before the next molecular layer is added at the growth front. For typical growth rates, in Regime 11, polyethylene lamellae form with disordered {001} fold surfaces then transform, with lamellar thickening and twisting, towards the more-ordered condition found for slower crystallization in Regime 1, in which lamellae form with and retain {201} fold surfaces. Several linear and linear-low-density polyethylenes have been used to show that, for the same polymer crystallized alone or in a blend, the growth rate at which the change in initial lamellar condition occurs is reasonably constant thereby supporting the concept of a specific time for surfaces to attain the ordered {201}) state. This specific time, in the range from milliseconds to seconds, increases with molecular length, and in linear-low-density polymer, for higher branch contents. (c) 2006 Elsevier Ltd. All rights reserved.

A new Snow-SVAT to simulate the accumulation and ablation of seasonal snow cover beneath a forest canopy

A new snow-soil-vegetation-atmosphere transfer (Snow-SVAT) scheme, which simulates the accumulation and ablation of the snow cover beneath a forest canopy, is presented. The model was formulated by coupling a canopy optical and thermal radiation model to a physically-based multi-layer snow model. This canopy radiation model is physically-based yet requires few parameters, so can be used when extensive in-situ field measurements are not available. Other forest effects such as the reduction of wind speed, interception of snow on the canopy and the deposition of litter were incorporated within this combined model, SNOWCAN, which was tested with data taken as part of the Boreal Ecosystem-Atmosphere Study (BOREAS) international collaborative experiment. Snow depths beneath four different canopy types and at an open site were simulated. Agreement between observed and simulated snow depths was generally good, with correlation coefficients ranging between r^2=0.94 and r^2=0.98 for all sites where automatic measurements were available. However, the simulated date of total snowpack ablation generally occurred later than the observed date. A comparison between simulated solar radiation and limited measurements of sub-canopy radiation at one site indicates that the model simulates the sub-canopy downwelling solar radiation early in the season to within measurement uncertainty.

Syntheses of two new 1D and 3D networks of Cu(II) and Co(II) using malonate and urotropine as bridging ligands: Crystal structures and magnetic studies

Two new metal-organic based polymeric complexes, [Cu-4(O2CCH2CO2)(4)(L)].7H(2)O (1) and [CO2(O2CCH2CO2)(2)(L)].2H(2)O (2) [L = hexamethylenetetramine (urotropine)], have been synthesized and characterized by X-ray crystal structure determination and magnetic studies. Complex 1 is a 1D coordination polymer comprising a carboxylato, bridged Cu-4 moiety linked by a tetradentate bridging urotropine. Complex 2 is a 3D coordination polymer made of pseudo-two-dimensional layers of Co(II) ions linked by malonate anions in syn-anticonformation which are bridged by bidentate urotropine in trans fashion, Complex 1 crystallizes in the orthothombic system, space group Pmmn, with a = 14,80(2) Angstrom, b = 14.54(2) Angstrom, c = 7.325(10) Angstrom, beta = 90degrees, and Z = 4. Complex 2 crystallizes in the orthorhombic system, space group Imm2, a = 7.584(11) Angstrom, b = 15.80(2) Angstrom, c = 6.939(13) Angstrom, beta = 90.10degrees(1), and Z = 4. Variable temperature (300-2 K) magnetic behavior reveals the existence of ferro- and antiferromagnetic interactions in 1 and only antiferromagnetic interactions in 2. The best fitted parameters for complex 1 are J = 13.5 cm(-1), J = -18.1 cm(-1), and g = 2.14 considering only intra-Cu-4 interactions through carboxylate and urotropine pathways. In case of complex 2, the fit of the magnetic data considering intralayer interaction through carboxylate pathway as well as interlayer interaction via urotropine pathway gave no satisfactory result at this moment using any model known due to considerable orbital contribution of Co(II) ions to the magnetic moment and its complicated structure. Assuming isolated Co(II) ions (without any coupling, J = 0) the shape of the chi(M)T curve fits well with experimental data except at very low temperatures.

A new porous 2D coordination polymer built by copper(II) and trimesic acid

A 2D porous material, Cu-3(tmen)(3)(tma)(2)(H2O)(2)(.)6.5H(2)O [tmen = N,N,N',N'-tetramethylethane-1,2-diamine; tmaH(3) = 1,3,5-benzenetricarboxylic acid/trimesic acid], has been synthesized and characterized by X-ray single crystal analysis, variable temperature magnetic measurements, IR spectra and XRPD pattern. The complex consists of 2D layers built by three crystallographically independent Cu(tmen) moieties bridged by tma anions. Of the three copper ions, Cu(1) and Cu(2) present distorted square pyramidal coordination geometry, while the third exhibits a severely distorted octahedral environment. The Cu(1)(tmen) and Cu(2)(tmen) building blocks bridged by tma anions give rise to chains with a zig-zag motif, which are cross-connected by Cu(3)(tmen)-tma polymers sharing metal ions Cu(2) through pendant tma carboxylates. The resulting 2D architecture extends in the crystallographic ab-plane. The adjacent sheets are embedded through the Cu(3)(tmen) tma chains, leaving H2O-filled channels. There are 6.5 lattice water molecules per formula unit, some of which are disordered. Upon heating, the lattice water molecules get eliminated without destroying the crystal morphology and the compound rehydrated reversibly on exposure to humid atmosphere. Magnetic data of the complex have been fitted considering isolated irregular Cu-3 triangles (three different J parameters) by applying the CLUMAG program. The best fit indicates three close comparable J parameters and very weak antiferromagnetic interactions are operative between the metal centers. (C) 2004 Elsevier B.V. All rights reserved.

A new cyano-bridged one-dimensional (GdFeIII)-Fe-III coordination polymer with o-phenanthroline as the blocking ligand: synthesis, structure, and magnetic properties

Synthesis, structural characterization, and magnetic properties of a new cyano-bridged one-dimensional iron (III)-gadolinium (III) compound, trans-[Gd(o-phen)(2)(H2O)(2)(mu-CN)(2)Fe(CN)(4)], - 2no-phen (o-phen = 1,10-phenanthroline), have been described. The compound crystallizes in the triclinic P (1) over bar space group with the following unit cell parameters: a = 10.538(14) angstrom, b = 12.004(14) angstrom, c = 20.61(2) angstrom, alpha = 92.41(1)degrees, beta = 92.76(1)degrees, gamma = 11 2.72(1)degrees, and Z = 2. In this complex, each gadolinium (III) is coordinated to two nitrile nitrogens of the CN groups coming from two different ferricyanides, the mutually trans cyanides of each of which links another different Gd-III to create -NC-Fe(CN)(4)-CN-Gd-NC- type 1-D chain structure. The one-dimensional chains are self-assembled in two-dimensions via weak C-H center dot center dot center dot N hydrogen bonds. Both the variable-temperature (2-300 K, 0.01 T and 0.8 T) and variable-field (0-50 000 Gauss, 2 K) magnetic measurements reveal the existence of very weak interaction in this molecule. The temperature dependence of the susceptibilities has been analyzed using a model for a chain of alternating classic (7/2) and quantum (1/2) spins. (c) 2005 Elsevier B.V. All rights reserved.

Two new mu-(1,3-azido)-bridged polymers: alternating single and double bridges in a 1D nickel(II) complex and uniform bridge in a 2D copper(II) complex: Syntheses, single-crystal structures and magnetic studies

Two polymeric azido bridged complexes [Ni2L2(N-3)(3)](n)(ClO4). (1) and [Cu(bpdS)(2)(N-3)],(ClO4),(H2O)(2.5n) (2) [L = Schiff base, obtained from the condensation of pyridine-2-aldehyde with N,N,2,2-tetramethyl-1,3-propanediamine; bpds = 4,4'-bipyridyl disulfide] have been synthesized and their crystal structures have been determined. Complex 1, C26H42ClN15Ni2O4, crystallizes in a triclinic system, space group P1 with a 8.089(13), b = 9.392(14), c = 12.267(18) angstrom, a = 107.28(l), b 95.95(1), gamma = 96.92(1)degrees and Z = 2; complex 2, C20H21ClCuN7O6.5S4, crystallizes in an orthorhombic system, space group Pnna with a = 10.839(14), b = 13.208(17), c = 19.75(2) angstrom and Z = 4. The crystal structure of I consists of 1D polymers of nickel(L) units, alternatively connected by single and double bridging mu-(1,3-N-3) ligand with isolated perchlorate anions. Variable temperature magnetic susceptibility data of the complex have been measured and the fitting,of magnetic data was carried out applying the Borris-Almenar formula for such types of alternating one-dimensional S = 1 systems, based on the Hamiltonian H = -J Sigma(S2iS2i-1 + aS(2i)S(2i+1)). The best-fit parameters obtained are J = -106.7 +/- 2 cm(-1); a = 0.82 +/- 0.02; g = 2.21 +/- 0.02. Complex 2 is a 2D network of 4,4 topology with the nodes occupied by the Cu-II ions, and the edges formed by single azide and double bpds connectors. The perchlorate anions are located between pairs of bpds. The magnetic data have been fitted considering the complex as a pseudo-one-dimensional system, with all copper((II)) atoms linked by [mu(1,3-azido) bridging ligands at axial positions (long Cu...N-3 distances) since the coupling through long bpds is almost nil. The best-fit parameters obtained with this model are J = -1.21 +/- 0.2 cm(-1), g 2.14 +/- 0.02. (c) Wiley-VCH Verlag GmbH & Co. KGaA, 69451 Weinheim, Germany, 2005).

Developing a method to characterize indoor environmental parameters in naturally ventilated classrooms

Until recently, there has been little investigation concerning the poor indoor air quality (IAQ) in classrooms. Despite the evidence that the educational building systems in many of the UK institutions have significant defects that may degrade IAQ, systematic assessments of IAQ measurements has been rarely undertaken. When undertaking IAQ measurement, there is a difficult task of representing and characterizing the environment parameters. Although technologies exist to measure these parameters, direct measurements especially in a naturally ventilated spaces are often difficult. This paper presents a methodology for developing a method to characterize indoor environment flow parameters as well as the Carbon Dioxide (CO2) concentrations. Thus, CO2 concentration level can be influenced by the differences in the selection of sampling points and heights. However, because this research focuses on natural ventilation in classrooms, air exchange is provided mainly by air infiltration. It is hoped that the methodology developed and evaluated in this research can effectively simplify the process of estimating the parameters for a systematic assessment of IAQ measurements in a naturally ventilated classrooms.

New mathematical modeling approach for predicting microbial inactivation by high hydrostatic pressure

A new primary model based on a thermodynamically consistent first-order kinetic approach was constructed to describe non-log-linear inactivation kinetics of pressure-treated bacteria. The model assumes a first-order process in which the specific inactivation rate changes inversely with the square root of time. The model gave reasonable fits to experimental data over six to seven orders of magnitude. It was also tested on 138 published data sets and provided good fits in about 70% of cases in which the shape of the curve followed the typical convex upward form. In the remainder of published examples, curves contained additional shoulder regions or extended tail regions. Curves with shoulders could be accommodated by including an additional time delay parameter and curves with tails shoulders could be accommodated by omitting points in the tail beyond the point at which survival levels remained more or less constant. The model parameters varied regularly with pressure, which may reflect a genuine mechanistic basis for the model. This property also allowed the calculation of (a) parameters analogous to the decimal reduction time D and z, the temperature increase needed to change the D value by a factor of 10, in thermal processing, and hence the processing conditions needed to attain a desired level of inactivation; and (b) the apparent thermodynamic volumes of activation associated with the lethal events. The hypothesis that inactivation rates changed as a function of the square root of time would be consistent with a diffusion-limited process.

Robust identification for linear-in-the-parameters models

In this paper new robust nonlinear model construction algorithms for a large class of linear-in-the-parameters models are introduced to enhance model robustness, including three algorithms using combined A- or D-optimality or PRESS statistic (Predicted REsidual Sum of Squares) with regularised orthogonal least squares algorithm respectively. A common characteristic of these algorithms is that the inherent computation efficiency associated with the orthogonalisation scheme in orthogonal least squares or regularised orthogonal least squares has been extended such that the new algorithms are computationally efficient. A numerical example is included to demonstrate effectiveness of the algorithms. Copyright (C) 2003 IFAC.

A new solution of the rendering equation with stratified Monte Carlo approach

This paper is turned to the advanced Monte Carlo methods for realistic image creation. It offers a new stratified approach for solving the rendering equation. We consider the numerical solution of the rendering equation by separation of integration domain. The hemispherical integration domain is symmetrically separated into 16 parts. First 9 sub-domains are equal size of orthogonal spherical triangles. They are symmetric each to other and grouped with a common vertex around the normal vector to the surface. The hemispherical integration domain is completed with more 8 sub-domains of equal size spherical quadrangles, also symmetric each to other. All sub-domains have fixed vertices and computable parameters. The bijections of unit square into an orthogonal spherical triangle and into a spherical quadrangle are derived and used to generate sampling points. Then, the symmetric sampling scheme is applied to generate the sampling points distributed over the hemispherical integration domain. The necessary transformations are made and the stratified Monte Carlo estimator is presented. The rate of convergence is obtained and one can see that the algorithm is of super-convergent type.

A new RBF neural network with boundary value constraints

We present a novel topology of the radial basis function (RBF) neural network, referred to as the boundary value constraints (BVC)-RBF, which is able to automatically satisfy a set of BVC. Unlike most existing neural networks whereby the model is identified via learning from observational data only, the proposed BVC-RBF offers a generic framework by taking into account both the deterministic prior knowledge and the stochastic data in an intelligent manner. Like a conventional RBF, the proposed BVC-RBF has a linear-in-the-parameter structure, such that it is advantageous that many of the existing algorithms for linear-in-the-parameters models are directly applicable. The BVC satisfaction properties of the proposed BVC-RBF are discussed. Finally, numerical examples based on the combined D-optimality-based orthogonal least squares algorithm are utilized to illustrate the performance of the proposed BVC-RBF for completeness.

Using modelling tools to optimise the calculation and subsequent reduction of cable distortion parameters with application to oil well logging systems

This paper discusses how the use of computer-based modelling tools has aided the design of a telemetry unit for use with oil well logging. With the aid of modern computer-based simulation techniques, the new design is capable of operating at data rates of 2.5 times faster than previous designs.

A new approach to reduced-order modelling

An error polynomial is defined, the coefficients of which indicate the difference at any instant between a system and a model of lower order approximating the system. It is shown how Markov parameters and time series proportionals of the model can be matched with those of the system by setting error polynomial coefficients to zero. Also discussed is the way in which the error between system and model can be considered as being a filtered form of an error input function specified by means of model parameter selection.