An approach to verifying and debugging simulation models governed by ordinary differential equations: Part 2. Residuals analysis and a case study

In Part 1 of this paper a methodology for back-to-back testing of simulation software was described. Residuals with error-dependent geometric properties were generated. A set of potential coding errors was enumerated, along with a corresponding set of feature matrices, which describe the geometric properties imposed on the residuals by each of the errors. In this part of the paper, an algorithm is developed to isolate the coding errors present by analysing the residuals. A set of errors is isolated when the subspace spanned by their combined feature matrices corresponds to that of the residuals. Individual feature matrices are compared to the residuals and classified as 'definite', 'possible' or 'impossible'. The status of 'possible' errors is resolved using a dynamic subset testing algorithm. To demonstrate and validate the testing methodology presented in Part 1 and the isolation algorithm presented in Part 2, a case study is presented using a model for biological wastewater treatment. Both single and simultaneous errors that are deliberately introduced into the simulation code are correctly detected and isolated. Copyright (C) 2003 John Wiley Sons, Ltd.

A comprehensive insight into floc characteristics and their impact on compressibility and settleability of activated sludge

This paper presents a comprehensive study of sludge floc characteristics and their impact on compressibility and settleability of activated sludge in full scale wastewater treatment processes. The sludge flocs were characterised by morphological (floc size distribution, fractal dimension, filament index), physical (flocculating ability, viscosity, hydrophobicity and surface charge) and chemical (polymeric constituents and metal content) parameters. Compressibility and settleability were defined in terms of the sludge volume index (SVI) and zone settling velocity (ZSV). The floc morphological and physical properties had important influence on the sludge compressibility and settleability. Sludges containing large flocs and high quantities of filaments, corresponding to lower values of fractal dimension (D-f), demonstrated poor compressibility and settleability. Sludge flocs with high flocculating ability had lower SVI and higher ZSV, whereas high values of hydrophobicity, negative surface charge and viscosity of the sludge flocs correlated to high SVI and low ZSV. The quantity of the polymeric compounds protein. humic substances and carbohydrate in the sludge and the extracted extracellular polymeric substances (EPS) had significant positive correlations with SVI. The ZSV was quantitatively independent of the polymeric constituents. High concentrations of the extracted EPS were related to poor compressibility and settleability. The cationic ions Ca, Mg, Al and Fe in the sludge improved significantly the sludge compressibility and settleability. (C) 2003 Elsevier Science B.V. All rights reserved.

Transdermal penetration of vasoconstrictors - Present understanding and assessment of the human epidermal flux and retention of free bases and ion-pairs

Purpose. As reductions in dermal clearance increase the residence time of solutes in the skin and underlying tissues we compared the topical penetration of potentially useful vasoconstrictors (VCs) through human epidermis as both free bases and ion-pairs with salicylic acid (SA). Methods. We determined the in vitro epidermal flux of ephedrine, naphazoline, oxymetazoline, phenylephrine, and xylometazoline applied as saturated solutions in propylene glycol: water (1: 1) and of ephedrine, naphazoline and tetrahydrozoline as 10% solutions of 1: 1 molar ratio ion-pairs with SA in liquid paraffin. Results. As free bases, ephedrine had the highest maximal flux, Jmax = 77.4 +/- 11.7 mug/cm(2)/h, being 4-fold higher than tetrahydrozoline and xylometazoline, 6-fold higher than phenylephrine, 10-fold higher than naphazoline and 100-fold higher than oxymetazoline. Stepwise regression of solute physicochemical properties identified melting point as the most significant predictor of flux. As ion-pairs with SA, ephedrine and naphazoline had similar fluxes (11.5 +/- 2.3 and 12.0 +/- 1.6 mug/cm(2)/h respectively), whereas tetrahydrozoline was approximately 3-fold slower. Corresponding fluxes of SA from the ion-pairs were 18.6 +/- 0.6, 7.8 +/- 0.8 and 1.1 +/- 0.1 respectively. Transdermal transport of VC's is discussed. Conclusions. Epidermal retention of VCs and SA did not correspond to their molar ratio on application and confirmed that following partitioning into the stratum corneum, ion-pairs separate and further penetration is governed by individual solute characteristics.

Structure of the HERG K+ channel S5P extracellular linker - Role of an amphipathic alpha-helix in C-type inactivation

The HERG K+ channel has very unusual kinetic behavior that includes slow activation but rapid inactivation. These features are critical for normal cardiac repolarization as well as in preventing lethal ventricular arrhythmias. Mutagenesis studies have shown that the extracellular peptide linker joining the fifth transmembrane domain to the pore helix is critical for rapid inactivation of the HERG K+ channel. This peptide linker is also considerably longer in HERG K+ channels, 40 amino acids, than in most other voltage-gated K+ channels. In this study we show that a synthetic 42-residue peptide corresponding to this linker region of the HERG K+ channel does not have defined structural elements in aqueous solution; however, it displays two well defined helical regions when in the presence of SDS micelles. The helices correspond to Trp(585)-Ile(593) and Gly(604)-Tyr(611) of the channel. The Trp(585)-Ile(593) helix has distinct hydrophilic and hydrophobic surfaces. The Gly(604)-Tyr(611) helix corresponds to an N-terminal extension of the pore helix. Electrophysiological studies of HERG currents following application of exogenous S5P peptides show that the amphipathic helix in the S5P linker interacts with the pore region of the channel in a voltage-dependent manner.

Impact of bronchopulmonary dysplasia, brain injury, and severe retinopathy on the outcome of extremely low-birth-weight infants at 18 months

CONTEXT: Despite more than 2 decades of outcomes research after very preterm birth, clinicians remain uncertain about the extent to which neonatal morbidities predict poor long-term outcomes of extremely low-birth-weight (ELBW) infants. OBJECTIVE: To determine the individual and combined prognostic effects of bronchopulmonary dysplasia (BPD), ultrasonographic signs of brain injury, and severe retinopathy of prematurity (ROP) on 18-month outcomes of ELBW infants. DESIGN: Inception cohort assembled for the Trial of Indomethacin Prophylaxis in Preterms (TIPP). SETTING AND PARTICIPANTS: A total of 910 infants with birth weights of 500 to 999 g who were admitted to 1 of 32 neonatal intensive care units in Canada, the United States, Australia, New Zealand, and Hong Kong between 1996 and 1998 and who survived to a postmenstrual age of 36 weeks. MAIN OUTCOME MEASURES: Combined end point of death or survival to 18 months with 1 or more of cerebral palsy, cognitive delay, severe hearing loss, and bilateral blindness. RESULTS: Each of the neonatal morbidities was similarly and independently correlated with a poor 18-month outcome. Odds ratios were 2.4 (95% confidence interval [CI], 1.8-3.2) for BPD, 3.7 (95% CI, 2.6-5.3) for brain injury, and 3.1 (95% CI, 1.9-5.0) for severe ROP. In children who were free of BPD, brain injury, and severe ROP the rate of poor long-term outcomes was 18% (95% CI, 14%-22%). Corresponding rates with any 1, any 2, and all 3 neonatal morbidities were 42% (95% CI, 37%-47%), 62% (95% CI, 53%-70%), and 88% (64%-99%), respectively. CONCLUSION: In ELBW infants who survive to a postmenstrual age of 36 weeks, a simple count of 3 common neonatal morbidities strongly predicts the risk of later death or neurosensory impairment.

Water table waves in an unconfined aquifer: Experiments and modeling

[1] Comprehensive measurements are presented of the piezometric head in an unconfined aquifer during steady, simple harmonic oscillations driven by a hydrostatic clear water reservoir through a vertical interface. The results are analyzed and used to test existing hydrostatic and nonhydrostatic, small-amplitude theories along with capillary fringe effects. As expected, the amplitude of the water table wave decays exponentially. However, the decay rates and phase lags indicate the influence of both vertical flow and capillary effects. The capillary effects are reconciled with observations of water table oscillations in a sand column with the same sand. The effects of vertical flows and the corresponding nonhydrostatic pressure are reasonably well described by small-amplitude theory for water table waves in finite depth aquifers. That includes the oscillation amplitudes being greater at the bottom than at the top and the phase lead of the bottom compared with the top. The main problems with respect to interpreting the measurements through existing theory relate to the complicated boundary condition at the interface between the driving head reservoir and the aquifer. That is, the small-amplitude, finite depth expansion solution, which matches a hydrostatic boundary condition between the bottom and the mean driving head level, is unrealistic with respect to the pressure variation above this level. Hence it cannot describe the finer details of the multiple mode behavior close to the driving head boundary. The mean water table height initially increases with distance from the forcing boundary but then decreases again, and its asymptotic value is considerably smaller than that previously predicted for finite depth aquifers without capillary effects. Just as the mean water table over-height is smaller than predicted by capillarity-free shallow aquifer models, so is the amplitude of the second harmonic. In fact, there is no indication of extra second harmonics ( in addition to that contained in the driving head) being generated at the interface or in the interior.

Asymmetric contribution of α and β subunits to the activation of αβ heteromeric glycine receptors

This study investigated the role of beta subunits in the activation of alphabeta heteromeric glycine receptor (GlyR) chloride channels recombinantly expressed in HEK293 cells. The approach involved incorporating mutations into corresponding positions in alpha and beta subunits and comparing their effects on receptor function. Although cysteine-substitution mutations to residues in the N-terminal half of the alpha subunit M2-M3 loop dramatically impaired the gating efficacy, the same mutations exerted little effect when incorporated into corresponding positions of the beta subunit. Furthermore, although the alpha subunit M2-M3 loop cysteines were modified by a cysteine-specific reagent, the corresponding beta subunit cysteines showed no evidence of reactivity. These observations suggest structural or functional differences between alpha and beta subunit M2-M3 loops. In addition, a threonine-->leucine mutation at the 9' position in the beta subunit M2 pore-lining domain dramatically increased the glycine sensitivity. By analogy with the effects of the same mutation in other ligand-gated ion channels, it was concluded that the mutation affected the GlyR activation mechanism. This supports the idea that the GlyR beta subunit is involved in receptor gating. In conclusion, this study demonstrates that beta subunits contribute to the activation of the GlyR, but that their involvement in this process is significantly different to that of the alpha subunit.

A class B switch-mode assisted linear amplifier

A switch-mode assisted linear amplifier (SMALA) combining a linear (Class B) and a switch-mode (Class D) amplifier is presented. The usual single hysteretic controlled half-bridge current dumping stage is replaced by two parallel buck converter stages, in a parallel voltage controlled topology. These operate independently: one buck converter sources current to assist the upper Class B output device, and a complementary converter sinks current to assist the lower device. This topology lends itself to a novel control approach of a dead-band at low power levels where neither class D amplifier assists, allowing the class B amplifier to supply the load without interference, ensuring high fidelity. A 20 W implementation demonstrates 85% efficiency, with distortion below 0.08% measured across the full audio bandwidth at 15 W. The class D amplifier begins assisting at 2 W, and below this value, the distortion was below 0.03%. Complete circuitry is given, showing the simplicity of the additional class D amplifier and its corresponding control circuitry.

Dissolution of sucrose crystals in the anhydrous sorbitol melt

The dissolution of a sugar (sucrose as a model) with higher melting point was studied in a molten food polyol (sorbitol as a model) with lower melting point, both in anhydrous state. A DSC and optical examination revealed the dissolution of anhydrous sucrose crystals (mp 192 degreesC) in anhydrous sorbitol (mp 99 degreesC) liquid melt. The sucrose-sorbitol crystal mixtures at the proportions of 10, 30, 60, 100 and 150 g of sucrose per 100 g of sorbitol were heat scanned in a DSC to above melting endotherm of sorbitol but well below the onset temperature of melting of sucrose at three different temperatures 110, 130 and 150 degreesC. The heat scanning modes used were with or without isothermal holding. The dissolution of sucrose in the sorbitol liquid melt was manifested by an increase in the glass transition temperature of the melt and corresponding decrease in endothermic melting enthalpy of sucrose. At given experimental conditions, as high as 25 and 85% of sucrose dissolved in the sorbitol melt during 1 h of isothermal holding at 110 and 150 degreesC, respectively. Optical microscopic observation also clearly showed the reduction in the size of sucrose crystals in sorbitol melt during the isothermal holding at those temperatures. (C) 2003 Elsevier Science Ltd. All rights reserved.

G-values for scission and crosslinking on γ-radiolysis of Ultem at 303 K

The radiation chemical yields G(S) and G(X) for H-linking and Y-linking models for Ultem have been calculated from molecular weight analysis by gel permeation chromatography. These G-values have been compared with the G-values obtained from analysis of soluble fractions above the gel dose, Which have been reported in previous works. An analysis of the molecular weight data in terms of H-linking and Y-linking mechanisms yielded values of G(S-H) = 1.0 x 10(-3) and G(H) = 6.0 x 10(-3) and G(S-Y) = 1.3 x 10(-2) and G(Y) = 1.8 x 10(-2). The corresponding values obtained from the solubility data were G(SH) = 0.53 x 10(-2), G(H) = 1.39 x 10(-2), G(S-Y) = 4.2 x 10(-2) and G(Y) - 4.6 x 10(-2). The origin of the disagreement between the molecular weight and solubility values is not clear, but it could arise as a result of observed microgel formation below the reported gel dose of 0.13 MGy. Whether the crosslink mechanism proceeds by an H-linking or Y-linking process is also unclear and will require direct observation of the crosslinking structures.

Effect of pore-network connectivity on multicomponent adsorption of large molecules

The effect of pore-network connectivity on binary liquid-phase adsorption equilibria using the ideal adsorbed solution theory (LAST) was studied. The liquid-phase binary adsorption experiments used ethyl propionate, ethyl butyrate, and ethyl isovalerate as the adsorbates and commercial activated carbons Filtrasorb-400 and Norit ROW 0.8 as adsorbents. As the single-component isotherm, a modified Dubinin-Radushkevich equation was used. A comparison with experimental data shows that incorporating the connectivity of the pore network and considering percolation processes associated with different molecular sizes of the adsorptives in the mixture, as well as their different corresponding accessibility, can improve the prediction of binary adsorption equilibria using the LAST Selectivity of adsorption for the larger molecule in binary systems increases with an increase in the pore-network coordination number, as well with an increase in the mean pore width and in the spread of the pore-size distribution.

Activation energy distribution of thermal annealing of a bituminous coal

A bituminous coal was pyrolyzed in a nitrogen stream in an entrained flow reactor at various temperatures from 700 to 1475 degreesC. Char samples were collected at different positions along the reactor. Each collected sample was oxidized nonisothermally in a TGA for reactivity determination. The reactivity of the coal char was found to decrease rapidly with residence time until 0.5 s, after which it decreased only slightly. On the bases of the reactivity data at various temperatures, a new approach was utilized to obtaining the true activation energy distribution function for thermal annealing without the assumption of any distribution function form or a constant preexponential factor. It appears that the true activation energy distribution function consists of two separate parts corresponding to different temperature ranges, suggesting different mechanisms in different temperature ranges. Partially burnt coal chars were also collected along the reactor when the coal was oxidized in air at various temperatures from 700 to 1475 degreesC. The collected samples were analyzed for the residual carbon content and the specific reaction rate was estimated. The characteristic time of thermal deactivation was compared with that of oxidation under realistic conditions. The characteristic times were found to be close to each other, indicating the importance of thermal deactivation during combustion of the coal studied.

Prediction of multilayer adsorption and capillary condensation phenomena in cylindrical mesopores

MCM-41 periodic mesoporous silicates with a high degree of structural ordering are synthesized and used as model adsorbents to study the isotherm prediction of nitrogen adsorption. The nitrogen adsorption isotherm at 77 K for a macroporous silica is measured and used in high-resolution alpha(s)-plot comparative analysis to determine the external surface area, total surface area and primary mesopore volume of the MCM-41 materials. Adsorption equilibrium data of nitrogen on the different pore size MCM-41 samples (pore diameters from 2.40 to 4.92 nm) are also obtained. Based on the Broekhoff and de Boer' thermodynamic analysis, the nitrogen adsorption isotherms for the different pore size MCM-41 samples are interpreted using a novel strategy, in which the parameters of an empirical expression, used to represent the potential of interaction between the adsorbate and adsorbent, are obtained by fitting only the multilayer region prior to capillary condensation for C-16 MCM-41. Subsequently the entire isotherm, including the phase transition, is predicted for all the different pore size MCM-41 samples without any fitting. The results show that the prediction of multilayer adsorption and total adsorbed amount are in good agreement with the experimental isotherms. The predictions of the relative pressure corresponding to capillary equilibrium (coexistence) transition agree remarkably with experimental data on the adsorption branch even for hysteretic isotherms, confirming that this is the branch appropriate for pore size distribution analysis. The impact of pore radius on the adsorption film thickness and capillary coexistence pressure is also investigated, and found to agree with the experimental data. (C) 2003 Elsevier Inc. All rights reserved.

Variação radial da retratibilidade e densidade básica da madeira de Eucalyptus saligna Sm.

A ultra-estrutura e a composição química da madeira, bem como suas propriedades físicas e mecânicas, variam significativamente entre espécies, entre árvores de uma mesma espécie e, mesmo, entre diferentes partes de uma mesma árvore. Com este trabalho objetivou-se o estudo dos parâmetros de retratibilidade e de densidade básica da madeira Eucalyptus saligna, com idade de 16 anos, proveniente de talhões experimentais da EMBRAPA Florestas, de Colombo, Paraná. As amostras foram retiradas à altura do DAP de quatro posições eqüidistantes a partir da medula em direção à periferia, correspondendo a 0, 33, 66 e 100%, com dimensões nominais de 1,0 x 2,0 x 3,0 cm, sendo a última dimensão no sentido longitudinal. Elas foram mantidas em câmara fechada com ventilação, próximo de soluções salinas supersaturadas, com o objetivo de proporcionar diferentes condições de umidade relativa. Uma vez atingidas as distintas condições de umidade de equilíbrio, as amostras foram secas em estufa a 105 ºC e obtidos os dados de retratibilidade e densidade básica da madeira nas posições mencionadas. Constataram-se valores de contração volumétrica mais baixos na região medular, apresentando um acréscimo para as demais posições. Comportamento semelhante foi observado para os coeficientes das contrações lineares nas direções tangencial e radial. O fator anisotrópico foi consideravelmente mais elevado na região medular, decrescendo substancialmente em direção ao alburno. A densidade básica não mostrou sinais efetivos de estabilidade, apesar de mostrar tendência de aumento em direção à periferia do tronco.

Basis set convergence in Hartree-Fock calculations of some diatomic molecules containing first and second-row atoms

Investiga-se a convergência de conjuntos de bases em direção ao limite numérico da energia Hartree-Fock (HF) total para as seqüências hierárquicas dos conjuntos de bases XZP e ccpVXZ. Para as duas hierarquias, melhoramentos significativos são obtidos com cada incremento em X. Para estimar o limite do conjunto de base completo, uma forma exponencial foi usada. Entre as várias aproximações consideradas aqui, uma extrapolação exponencial de três parâmetros aplicada aos resultados TZP, QZP e 5ZP deu os limites do conjunto de bases mais precisos. Em adição, energias HF dos orbitais moleculares ocupados mais altos de algumas moléculas diatômicas foram calculadas com o conjunto 5ZP e comparadas com as correspondentes obtidas com o conjunto cc-pV5Z e com um método numérico HF.