Subcycling first- and second-order generalizations of the trapezoidal rule

Subcycling algorithms which employ multiple timesteps have been previously proposed for explicit direct integration of first- and second-order systems of equations arising in finite element analysis, as well as for integration using explicit/implicit partitions of a model. The author has recently extended this work to implicit/implicit multi-timestep partitions of both first- and second-order systems. In this paper, improved algorithms for multi-timestep implicit integration are introduced, that overcome some weaknesses of those proposed previously. In particular, in the second-order case, improved stability is obtained. Some of the energy conservation properties of the Newmark family of algorithms are shown to be preserved in the new multi-timestep extensions of the Newmark method. In the first-order case, the generalized trapezoidal rule is extended to multiple timesteps, in a simple way that permits an implicit/implicit partition. Explicit special cases of the present algorithms exist. These are compared to algorithms proposed previously. (C) 1998 John Wiley & Sons, Ltd.

The subcycled Newmark algorithm

The popular Newmark algorithm, used for implicit direct integration of structural dynamics, is extended by means of a nodal partition to permit use of different timesteps in different regions of a structural model. The algorithm developed has as a special case an explicit-explicit subcycling algorithm previously reported by Belytschko, Yen and Mullen. That algorithm has been shown, in the absence of damping or other energy dissipation, to exhibit instability over narrow timestep ranges that become narrower as the number of degrees of freedom increases, making them unlikely to be encountered in practice. The present algorithm avoids such instabilities in the case of a one to two timestep ratio (two subcycles), achieving unconditional stability in an exponential sense for a linear problem. However, with three or more subcycles, the trapezoidal rule exhibits stability that becomes conditional, falling towards that of the central difference method as the number of subcycles increases. Instabilities over narrow timestep ranges, that become narrower as the model size increases, also appear with three or more subcycles. However by moving the partition between timesteps one row of elements into the region suitable for integration with the larger timestep these the unstable timestep ranges become extremely narrow, even in simple systems with a few degrees of freedom. As well, accuracy is improved. Use of a version of the Newmark algorithm that dissipates high frequencies minimises or eliminates these narrow bands of instability. Viscous damping is also shown to remove these instabilities, at the expense of having more effect on the low frequency response.

Explicit/implicit partitioning and a new explicit form of the generalized alpha method

Most finite element packages use the Newmark algorithm for time integration of structural dynamics. Various algorithms have been proposed to better optimize the high frequency dissipation of this algorithm. Hulbert and Chung proposed both implicit and explicit forms of the generalized alpha method. The algorithms optimize high frequency dissipation effectively, and despite recent work on algorithms that possess momentum conserving/energy dissipative properties in a non-linear context, the generalized alpha method remains an efficient way to solve many problems, especially with adaptive timestep control. However, the implicit and explicit algorithms use incompatible parameter sets and cannot be used together in a spatial partition, whereas this can be done for the Newmark algorithm, as Hughes and Liu demonstrated, and for the HHT-alpha algorithm developed from it. The present paper shows that the explicit generalized alpha method can be rewritten so that it becomes compatible with the implicit form. All four algorithmic parameters can be matched between the explicit and implicit forms. An element interface between implicit and explicit partitions can then be used, analogous to that devised by Hughes and Liu to extend the Newmark method. The stability of the explicit/implicit algorithm is examined in a linear context and found to exceed that of the explicit partition. The element partition is significantly less dissipative of intermediate frequencies than one using the HHT-alpha method. The explicit algorithm can also be rewritten so that the discrete equation of motion evaluates forces from displacements and velocities found at the predicted mid-point of a cycle. Copyright (C) 2003 John Wiley Sons, Ltd.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Direct cadherin-activated cell signaling: a view from the plasma membrane

Classical cadherin adhesion molecules are key determinants of cell recognition and tissue morphogenesis, with diverse effects on cell behavior. Recent developments indicate that classical cadherins are adhesion-activated signaling receptors. In particular, early-immediate Rac signaling is emerging as a mechanism to coordinate cadherin-actin integration at the plasma membrane.

Role of Diadegma semiclausum (Hymenoptera: Ichneumonidae) in controlling Plutella xylostella (Lepidoptera: Plutellidae): Cage exclusion experiments and direct observation

We evaluated the role of the larval parasitoid, Diadegma semiclausum Hellen (Hymenoptera: Ichneumonidae), in controlling Plutella xylostella (L.) (Lepidoptera: Plutellidae) by cage exclusion experiments and direct field observation during the winter season in southern Queensland, Australia. The cage exclusion experiment involved uncaged, open cage and closed cage treatments. A higher percentage (54-83%) of P. xylostella larvae on sentinel plants were lost in the uncaged treatment than the closed (4-9%) or open cage treatments (11-29%). Of the larvae that remained in the uncaged treatment, 72-94% were parasitized by D. semiclausum , much higher than that in the open cage treatment (8-37% in first trial, and 38-63% in second trial). Direct observations showed a significant aggregation response of the field D. semiclausum populations to high host density plants in an experimental plot and to high host density plots that were artificially set-up near to the parasitoid source fields. The degree of aggregation varied in response to habitat quality of the parasitoid source field and scales of the manipulated host patches. As a result, density-dependence in the pattern of parasitism may depend on the relative degree of aggregation of the parasitoid population at a particular scale. A high degree of aggregation seems to be necessary to generate density-dependent parasitism by D. semiclausum . Integration of the cage exclusion experiment and direct observation demonstrated the active and dominant role of this parasitoid in controlling P. xylostella in the winter season. A biologically based IPM strategy, which incorporates the use of D. semiclausum with Bt, is suggested for the management of P. xylostella in seasons or regions with a mild temperature.

Nanocomposite Nafion-Silica membranes for direct methanol fuel cells

Commercially available proton exchange membranes such as Nafion do not meet the requirements for high power density direct methanol fuel cells, partly due to their high methanol permeability. The aim of this work is to develop a new class of high-proton conductivity membranes, with thermal and mechanical stability similar to Nafion and reduced methanol permeability. Nanocomposite membranes were produced by the in-situ sol-gel synthesis of silicon dioxide particles in preformed Nafion membranes. Microstructural modification of Nafion membranes with silica nanoparticles was shown in this work to reduce methanol crossover from 7.48x10-6 cm2s^-1 for pure Nafion® to 2.86 x10-6 cm2s^-1 for nanocomposite nafion membranes (Methanol 50% (v/v) solution, 75 degrees C). Best results were achieved with a silica composition of 2.6% (w/w). We propose that silica inhibits the conduction of methanol through Nafion by blocking sites necessary for methanol diffusion through the polymer electrolyte membrane. Effects of surface chemistry, nanoparticle formation and interactions with Nafion matrix are further addressed.

Chemical treatment of Escherichia coli. II. Direct extraction of recombinant protein from cytoplasmic inclusion bodies in intact cells

A method is presented for the direct extraction of the recombinant protein Long-R-3-IGF-I from inclusion bodies located in the cytoplasm of intact Escherichia coli cells. Chemical treatment with 6M urea, 3 mM EDTA, and 20 mM dithiothreitol (DTT) at pH 9.0 proved an effective combination for extracting recombinant protein from intact cells. Comparable levels of Long-R-3-IGF-I were recovered by direct extraction as achieved by in vitro dissolution following mechanical disruption. However, the purity of directly extracted recombinant protein was lower due to contamination by bacterial cell components. The kinetics of direct extraction are described using a first-order equation with the time constant of 3 min. Urea appears important for permeabilization of the cell and dissolution of the inclusion body. Conversely, EDTA is involved in permeabilization of the cell wall and DTT enhances protein release. pH proved to be important with lower levels of protein release achieved at low pH values (

The 1.1 angstrom resolution crystal structure of [Tyr(15)]EpI, a novel alpha-conotoxin from Conus episcopatus, solved by direct methods

Conotoxins are valuable probes of receptors and ion channels because of their small size and highly selective activity. alpha-Conotoxin EpI, a 16-residue peptide from the mollusk-hunting Conus episcopatus, has the amino acid sequence GCCSDPRCNMNNPDY(SO3H)C-NH2 and appears to be an extremely potent and selective inhibitor of the alpha 3 beta 2 and alpha 3 beta 4 neuronal subtypes of the nicotinic acetylcholine receptor (nAChR). The desulfated form of EpI ([Tyr(15)]EpI) has a potency and selectivity for the nAChR receptor similar to those of EpI. Here we describe the crystal structure of [Tyr(15)]EpI solved at a resolution of 1.1 Angstrom using SnB. The asymmetric unit has a total of 284 non-hydrogen atoms, making this one of the largest structures solved de novo try direct methods. The [Tyr(15)]EpI structure brings to six the number of alpha-conotoxin structures that have been determined to date. Four of these, [Tyr(15)]EpI, PnIA, PnIB, and MII, have an alpha 4/7 cysteine framework and are selective for the neuronal subtype of the nAChR. The structure of [Tyr(15)]EpI has the same backbone fold as the other alpha 4/7-conotoxin structures, supporting the notion that this conotoxin cysteine framework and spacing give rise to a conserved fold. The surface charge distribution of [Tyr(15)]EpI is similar to that of PnIA and PnIB but is likely to be different from that of MII, suggesting that [Tyr(15)]EpI and MII may have different binding modes for the same receptor subtype.

T-DNA tagging and characterisation of a novel meristem-specific promoter from tobacco

We have evaluated T-DNA mediated plant promoter tagging, with a left-border-linked promoterless firefly luciferase (luc) construct, as a strategy for the isolation of novel plant promoters. In a population of approximately 300 transformed tobacco plants, IO lines showed LUC activity, including novel tissue-specific and developmental patterns of expression. One line, showing LUC activity only in the shoot and root apical meristems, was further characterised. Inverse PCR was used to amplify a 1.5 kb fragment of plant DNA flanking the single-copy T-DNA insertion in this line. With the exception of a 249 bp highly repetitive element, this sequence is present as a single copy in the tobacco genome, and is not homologous to any previously characterised DNA sequences. Sequence analysis revealed the presence of several motifs that may be involved in transcriptional regulation. Transgenic tobacco plants transformed with a transcriptional fusion of this putative promoter sequence to the beta-glucuronidase (uidA) reporter gene, showed GUS activity confined to the shoot tip and mature pollen. This promoter may be useful to direct the expression of genes controlling the transition to flowering, or genes to reduce losses due to pests and stresses damaging plant apical meristems.

Critical sets in direct products of back circulant Latin squares

To date very Few families of critical sets for latin squares are known. The only previously known method for constructing critical sets involves taking a critical set which is known to satisfy certain strong initial conditions and using a doubling construction. This construction can be applied to the known critical sets in back circulant latin squares of even order. However, the doubling construction cannot be applied to critical sets in back circulant latin squares of odd order. In this paper a family of critical sets is identified for latin squares which are the product of the latin square of order 2 with a back circulant latin square of odd order. The proof that each element of the critical set is an essential part of the reconstruction process relies on the proof of the existence of a large number of latin interchanges.