A study of the stability of subcycling algorithms in structural dynamics

Algorithms for explicit integration of structural dynamics problems with multiple time steps (subcycling) are investigated. Only one such algorithm, due to Smolinski and Sleith has proved to be stable in a classical sense. A simplified version of this algorithm that retains its stability is presented. However, as with the original version, it can be shown to sacrifice accuracy to achieve stability. Another algorithm in use is shown to be only statistically stable, in that a probability of stability can be assigned if appropriate time step limits are observed. This probability improves rapidly with the number of degrees of freedom in a finite element model. The stability problems are shown to be a property of the central difference method itself, which is modified to give the subcycling algorithm. A related problem is shown to arise when a constraint equation in time is introduced into a time-continuous space-time finite element model. (C) 1998 Elsevier Science S.A.

Hydrothermal Stability of Modified Silica Membranes for Gas Separation

A new class of hybrid molecular sieve silica (MSS) membranes is developed and tested against standard and organic templated membranes. The hybrid membrane is synthesized by the standard sol-gel process, integrating a template (methyltriethoxysilane - MTES) and a C6 surfactant (triethylhexylammonium bromide) into the silica film matrix. After hydro treatment under a relative humidity of 96% for 50h, the hybrid membrane shows no changes in its gas separation capabilities or energy of mobility. The structural characteristics and integrity of the hybrid membrane are retained due to a high concentration of organophilic functional groups and alkoxides observed using 29 Si NMR. In contrast, the structural integrity of the membranes prepared with non-templated films deteriorated during the hydro treatment due to a large percentage of silanol groups (Si-OH) which react with water. The hybrid membranes underwent a decrease in the H2/CO2 selectivity of only 1% whereas for the non-templated membrane a 21% decrease was observed. The transport mechanism of the hybrid membranes is activated as permeation increased with temperature. The activation energy for the permeation of H2 is positive while negative for CO2. The H2 permeation obtained was 3x 10 -8 mol.m -2 .s -1 .Pa -1 and permselectivities for H2/CO2 and H2/N2 varied between 1-7 and 31-34, respectively.

A modified discrete random pore model allowing for different initial surface reactivity

We investigate here a modification of the discrete random pore model [Bhatia SK, Vartak BJ, Carbon 1996;34:1383], by including an additional rate constant which takes into account the different reactivity of the initial pore surface having attached functional groups and hydrogens, relative to the subsequently exposed surface. It is observed that the relative initial reactivity has a significant effect on the conversion and structural evolution, underscoring the importance of initial surface chemistry. The model is tested against experimental data on chemically controlled char oxidation and steam gasification at various temperatures. It is seen that the variations of the reaction rate and surface area with conversion are better represented by the present approach than earlier random pore models. The results clearly indicate the improvement of model predictions in the low conversion region, where the effect of the initially attached functional groups and hydrogens is more significant, particularly for char oxidation. It is also seen that, for the data examined, the initial surface chemistry is less important for steam gasification as compared to the oxidation reaction. Further development of the approach must also incorporate the dynamics of surface complexation, which is not considered here.

Algorithms for sequence analysis via mutagenesis

Despite many successes of conventional DNA sequencing methods, some DNAs remain difficult or impossible to sequence. Unsequenceable regions occur in the genomes of many biologically important organisms, including the human genome. Such regions range in length from tens to millions of bases, and may contain valuable information such as the sequences of important genes. The authors have recently developed a technique that renders a wide range of problematic DNAs amenable to sequencing. The technique is known as sequence analysis via mutagenesis (SAM). This paper presents a number of algorithms for analysing and interpreting data generated by this technique.

Comparison of BR and QR Eigenvalue Algorithms for Power System Small Signal Stability Analysis

The BR algorithm is a novel and efficient method to find all eigenvalues of upper Hessenberg matrices and has never been applied to eigenanalysis for power system small signal stability. This paper analyzes differences between the BR and the QR algorithms with performance comparison in terms of CPU time based on stopping criteria and storage requirement. The BR algorithm utilizes accelerating strategies to improve its performance when computing eigenvalues of narrowly banded, nearly tridiagonal upper Hessenberg matrices. These strategies significantly reduce the computation time at a reasonable level of precision. Compared with the QR algorithm, the BR algorithm requires fewer iteration steps and less storage space without depriving of appropriate precision in solving eigenvalue problems of large-scale power systems. Numerical examples demonstrate the efficiency of the BR algorithm in pursuing eigenanalysis tasks of 39-, 68-, 115-, 300-, and 600-bus systems. Experiment results suggest that the BR algorithm is a more efficient algorithm for large-scale power system small signal stability eigenanalysis.

Extended GCD and Hermite normal form algorithms via lattice basis reduction

Extended gcd calculation has a long history and plays an important role in computational number theory and linear algebra. Recent results have shown that finding optimal multipliers in extended gcd calculations is difficult. We present an algorithm which uses lattice basis reduction to produce small integer multipliers x(1), ..., x(m) for the equation s = gcd (s(1), ..., s(m)) = x(1)s(1) + ... + x(m)s(m), where s1, ... , s(m) are given integers. The method generalises to produce small unimodular transformation matrices for computing the Hermite normal form of an integer matrix.

The Modified Card Sorting Test: Test-retest stability and relationships with demographic variables in a healthy older adult sample

Objectives. To investigate the test-retest stability of a standardized version of Nelson's (1976) Modified Card Sorting Test (MCST) and its relationships with demographic variables in a sample of healthy older adults. Design. A standard card order and administration were devised for the MCST and administered to participants at an initial assessment, and again at a second session conducted a minimum of six months later in order to examine its test-retest stability. Participants were also administered the WAIS-R at initial assessment in order to provide a measure of psychometric intelligence. Methods. Thirty-six (24 female, 12 male) healthy older adults aged 52 to 77 years with mean education 12.42 years (SD = 3.53) completed the MCST on two occasions approximately 7.5 months (SD = 1.61) apart. Stability coefficients and test-retest differences were calculated for the range of scores. The effect of gender on MCST performance was examined. Correlations between MCST scores and age, education and WAIS-R IQs were also determined. Results. Stability coefficients ranged from .26 for the percent perseverative errors measure to .49 for the failure to maintain set measure. Several measures were significantly correlated with age, education and WAIS-R IQs, although no effect of gender on MCST performance was found. Conclusions. None of the stability coefficients reached the level required for clinical decision making. The results indicate that participants' age, education, and intelligence need to be considered when interpreting MCST performance. Normative studies of MCST performance as well as further studies with patients with executive dysfunction are needed.

Comparison of carbohydrate-deficient transferrin, immunoglobulin A antibodies reactive with acetaldehyde-modified protein and acetaldehyde-modified albumin with conventional markers of alcohol consumption

Carbohydrate-deficient transferrin (CDT) has emerged as the best new marker for alcohol abuse. Recently plasma immunoglobulin A (IgA) reactivity with acetaldehyde (AcH)-modified proteins, or the modified proteins per se, have been proposed as a markers for high levels of alcohol consumption. In this study, we have compared CDT, IgA reactivity with AcH adducts (IgA ASR), and AcH-modified albumin with conventional markers of high alcohol intake in groups with well-defined drinking histories, The plasma activity of ALT, AST, and gamma-glutamyltransferase increased steadily with increasing alcohol consumption, CDT and AcH-modified albumin showed a similar pattern, whereas IgA ASR appeared only to be elevated after a threshold level of consumption had been reached, Neither CDT IgA ASR or AcH-modified albumin correlated strongly with any of the conventional markers or each other. This study shows that CDT, IgA ASR, AcH-modified albumin, and the conventional markers are not related, but suggests that the concurrent use of CDT and IgA ASR may lead to better identification of high alcohol intake.

Effects of data characteristics on the results of reserve selection algorithms

We tested the effects of four data characteristics on the results of reserve selection algorithms. The data characteristics were nestedness of features (land types in this case), rarity of features, size variation of sites (potential reserves) and size of data sets (numbers of sites and features). We manipulated data sets to produce three levels, with replication, of each of these data characteristics while holding the other three characteristics constant. We then used an optimizing algorithm and three heuristic algorithms to select sites to solve several reservation problems. We measured efficiency as the number or total area of selected sites, indicating the relative cost of a reserve system. Higher nestedness increased the efficiency of all algorithms (reduced the total cost of new reserves). Higher rarity reduced the efficiency of all algorithms (increased the total cost of new reserves). More variation in site size increased the efficiency of all algorithms expressed in terms of total area of selected sites. We measured the suboptimality of heuristic algorithms as the percentage increase of their results over optimal (minimum possible) results. Suboptimality is a measure of the reliability of heuristics as indicative costing analyses. Higher rarity reduced the suboptimality of heuristics (increased their reliability) and there is some evidence that more size variation did the same for the total area of selected sites. We discuss the implications of these results for the use of reserve selection algorithms as indicative and real-world planning tools.

Protein sequence threading, the alignment problem, and a two-step strategy

Conventionally, protein structure prediction via threading relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. (C) 1999 John Wiley & Sons, Inc.

Characterization of the structural and surface properties of chemically modified MCM-41 material

Mesoporous Mobil catalytic materials of number 41 (MCM-41) silica was chemically modified using both inorganic and organic precursors and characterized using the techniques, XRD, XPS, MAS NMR, FTIR, W-Vis, and physical adsorption of nitrogen, hydrocarbons (hexane, benzene, acetone, and methanol) and water vapor. Modification using organic reagents was found to result in a significant loss in porosity and a shape change of surface properties (increased hydrophobicity and decreased acidity). With inorganic modifying reagents, the decrease in porosity was also observed while the surface properties were not significantly altered as reflected by the adsorption isotherms of organics and water vapors. Chemical modifications can greatly improve the hydrothermal stability of MCM-41 material because of the enhanced surface hydrophobicity (with organic modifiers) or increased pore wall thickness (with inorganic modifiers). (C) 2000 Elsevier Science B.V. All rights reserved.

Long-term quantification of the release of monomers from dental resin composites and a resin-modified glass ionomer cement

This study quantified the release of monomers from polymerized specimens of four commercially available resin composites and one glass ionomer cement immersed in water:ethanol solutions. Individual standard curves were prepared from five monomers: (1) triethylene glycol dimethacrylate (TEGDMA), (2) 2-hydroxy-ethyl methacrylate (HEMA), (3) urethane dimethacrylate (UDMA), (4) bisphenol A glycidyl dimethacrylate (BISGMA), and (5) bisphenol A. The concentration of the monomers was determined at Days 1, 7, 30, and 90 with the use of electrospray ionization/mass spectrometry. Data were expressed in mean mumol per mm(2) surface area of specimen and analyzed with Scheffe's test (P < 0.05). The following monomers were found in water: monomers (1) and (2) from Delton sealant, monomer (5) from ScotchBond Multipurpose Adhesive and Delton sealant, monomer (3) from Definite and monomer (4) from Fuji II LC, ScotchBond Multipurpose Adhesive, Synergy and Definite. All these monomers increased in concentration over time, with the exception of monomer (1) from Delton sealant. Monomers (3) and (5) were found in extracts of materials despite their absence from the manufacturer's published composition. All monomers were released in significantly higher concentrations in water:ethanol solutions than in water. The greatest release of monomers occurred in the first day. The effect of the measured concentrations of monomers (1-5) on human genes, cells, or tissues needs to be considered with the use of a biological model. (C) 2002 Wiley Periodicals, Inc.

Genetic algorithms for design of liquid retaining structure

In this paper, genetic algorithm (GA) is applied to the optimum design of reinforced concrete liquid retaining structures, which comprise three discrete design variables, including slab thickness, reinforcement diameter and reinforcement spacing. GA, being a search technique based on the mechanics of natural genetics, couples a Darwinian survival-of-the-fittest principle with a random yet structured information exchange amongst a population of artificial chromosomes. As a first step, a penalty-based strategy is entailed to transform the constrained design problem into an unconstrained problem, which is appropriate for GA application. A numerical example is then used to demonstrate strength and capability of the GA in this domain problem. It is shown that, only after the exploration of a minute portion of the search space, near-optimal solutions are obtained at an extremely converging speed. The method can be extended to application of even more complex optimization problems in other domains.