Asymptotic and finite element analyses of mode III dynamic crack growth at a ductile-brittle interface

In this work, dynamic crack growth along a ductile-brittle interface under anti-plane strain conditions is studied. The ductile solid is taken to obey the J(2) flow theory of plasticity with linear isotropic strain hardening, while the substrate is assumed to exhibit linear elastic behavior. Firstly, the asymptotic near-tip stress and velocity fields are derived. These fields are assumed to be variable-separable with a power singularity in the radial coordinate centered at the crack tip. The effects of crack speed, strain hardening of the ductile phase and mismatch in elastic moduli of the two phases on the singularity exponent and the angular functions are studied. Secondly, full-field finite element analyses of the problem under small-scale yielding conditions are performed. The validity of the asymptotic fields and their range of dominance are determined by comparing them with the results of the full-field finite element analyses. Finally, theoretical predictions are made of the variations of the dynamic fracture toughness with crack velocity. The influence of the bi-material parameters on the above variation is investigated.

Dynamic stiffness matrix of a general cable element

A computational scheme for determining the dynamic stiffness coefficients of a linear, inclined, translating and viscously/hysteretically damped cable element is outlined. Also taken into account is the coupling between inplane transverse and longitudinal forms of cable vibration. The scheme is based on conversion of the governing set of quasistatic boundary value problems into a larger equivalent set of initial value problems, which are subsequently numerically integrated in a spatial domain using marching algorithms. Numerical results which bring out the nature of the dynamic stiffness coefficients are presented. A specific example of random vibration analysis of a long span cable subjected to earthquake support motions modeled as vector gaussian random processes is also discussed. The approach presented is versatile and capable of handling many complicating effects in cable dynamics in a unified manner.

Lectin binding to complex carbohydrate at the interface: a study by resonance energy transfer

The binding affinity of the oligosaccharide moiety of a neutral glycosphingolipid, asialoGM1, towards Ricinus communis agglutinin (RCAI) was determined for the first time by fluorescence resonance energy transfer (RET). The asialoGM1 was incorporated into a phospholipid (DMPC) vesicle doped with dansylated DPPE and then titrated with an increasing amount of the galactose specific RCAI. The efficiency of RET was determined by a saturable increase in the quenching of 'donor' fluorescence, i.e. the 'trp' residue of RCAI, due to the energy transfer from the 'acceptor' dansyl group on the surface of the vesicle. The apparent binding constant was found to be in the range of 10(5)-10(6) M-1 at 27 degrees C.

FIR system identification based on subspaces of a higher order cumulant matrix

A new linear algebraic approach for identification of a nonminimum phase FIR system of known order using only higher order (>2) cumulants of the output process is proposed. It is first shown that a matrix formed from a set of cumulants of arbitrary order can be expressed as a product of structured matrices. The subspaces of this matrix are then used to obtain the parameters of the FIR system using a set of linear equations. Theoretical analysis and numerical simulation studies are presented to characterize the performance of the proposed methods.

A transfer matrix approach for evaluation of the response of a multi-layer infinite plate to a two-dimensional pressure excitation

A 6 X 6 transfer matrix is presented to evaluate the response of a multi-layer infinite plate to a given two-dimensional pressure excitation on one of its faces or, alternatively, to evaluate the acoustic pressure distribution excited by the normal velocity components of the radiating surfaces. It is shown that the present transfer matrix is a general case embodying the transfer matrices of normal excitation and one-dimensional pressure excitation due to an oblique incident wave. It is also shown that the present transfer matrix obeys the necessary checks to categorize the physically symmetric multi-layer plate as dynamically symmetric. Expressions are derived to obtain the wave propagation parameters, such as the transmission, absorption and reflection coefficients, in terms of the elements of the transfer matrix presented. Numerical results for transmission loss and reflection coefficients of a two-layer configuration are presented to illustrate the effect of angles of incidence, layer characteristics and ambient media.

Calculation of elastic pinning strength of semicoherent Y2BaCuO5 precipitates in YBa2Cu3O7

The pinning energy due to the elastic interaction of a semicoherent Y2BaCuO5 precipitate with the YBa2Cu3O7 matrix is computed. This is achieved by setting up dislocation arrays at the interface. The elastic stresses generated by such arrays are integrated over a fluxoid volume to obtain the energy. It is seen that this elastic interaction energy makes an additive contribution to the total J(c) value.

Pipelined ring algorithm for matrix multiplication on transputer networks: performance analysis and estimation

A parallel matrix multiplication algorithm is presented, and studies of its performance and estimation are discussed. The algorithm is implemented on a network of transputers connected in a ring topology. An efficient scheme for partitioning the input matrices is introduced which enables overlapping computation with communication. This makes the algorithm achieve near-ideal speed-up for reasonably large matrices. Analytical expressions for the execution time of the algorithm have been derived by analysing its computation and communication characteristics. These expressions are validated by comparing the theoretical results of the performance with the experimental values obtained on a four-transputer network for both square and irregular matrices. The analytical model is also used to estimate the performance of the algorithm for a varying number of transputers and varying problem sizes. Although the algorithm is implemented on transputers, the methodology and the partitioning scheme presented in this paper are quite general and can be implemented on other processors which have the capability of overlapping computation with communication. The equations for performance prediction can also be extended to other multiprocessor systems.

Sequential Assembly of an Active RNA Polymerase Molecule at the Air-Water Interface

At the heart of understanding cellular processes lies our ability to explore the specific nature of communication between sequential information carrying biopolymers. However, the data extracted from conventional solution phase studies may not reflect the dynamics of communication between recognized partners as they occur in the crowded cellular milieu. We use the principle of immobilization of histidine-tagged biopolymers at a Ni(II)-encoded Langmuir monolayer to study sequence-specific protein-protein interactions in an artificially crowded environment The advantage of this technique lies in increasing the surface density of one of the interacting partners that allows us to study macromolecular interactions in a controlled crowded environment, but without compromising the speed of the reactions. We have taken advantage of this technique to follow the sequential assembly process of the multiprotein complex Escherichia coil RNA polymerase at the interface and also deciphered the role of one of the proteins, omega (omega), in the assembly pathway. Our reconstitution studies indicate that in the absence of molecular chaperones or other cofactors, omega (omega) plays a decisive role in refolding the largest protein beta prime (beta') and its recruitment into the multimeric assembly to reconstitute an active RNA polymerase. It was also observed that the monolayer had the ability to distinguish between sequence-specific and -nonspecific interactions despite the immobilization of one of the biomacromolecules. The technique provides a universal two-dimensional template for studying protein-ligand interactions while mimicking molecular crowding.

Density-matrix renormalization-group studies of the spin-1/2 Heisenberg system with dimerization and frustration

Using the density-matrix renormalization-group technique, we study the ground-state phase diagram and other low-energy properties of an isotropic antiferromagnetic spin-1/2 chain with both dimerization and frustration, i.e., an alternation delta of the nearest-neighbor exchanges and a next-nearest-neighbor exchange J(2). For delta = 0, the system is gapless for J(2) < J(2c) and has a gap for J(2) > J(2c) where J(2c) is about 0.241. For J(2) = J(2c) the gap above the ground state grows as delta to the power 0.667 +/- 0.001. In the J(2)-delta plane, there is a disorder line 2J(2) + delta = 1. To the left of this line, the peak in the static structure factor S(q) is at q(max) = pi (Neel phase), while to the right of the line, q(max) decreases from pi to pi/2 as J(2) is increased to large values (spiral phase). For delta = 1, the system is equivalent to two coupled chains as on a ladder and it is gapped for all values of the interchain coupling.

Engineered dimer interface mutants of triosephosphate isomerase: the role of inter-subunit interactions in enzyme function and stability

The role of inter-subunit interactions in maintaining optimal catalytic activity in triosephosphate isomerase (TIM) has been probed, using the Plasmodium falciparum enzyme as a model. Examination of subunit interface contacts in the crystal structures suggests that residue 75 (Thr, conserved) and residue 13 (Cys, variable) make the largest number of inter-subunit contacts. The mutants Cys13Asp (C13D) and Cys13Glu (C13E) have been constructed and display significant reduction in catalytic activity when compared with wild-type (WT) enzyme (similar to 7.4-fold decrease in k(cat) for the C13D and similar to 3.3-fold for the C13E mutants). Analytical gel filtration demonstrates that the C13D mutant dissociates at concentrations < 1.25 mu M, whereas the WT and the C13E enzymes retain the dimeric structure. The order of stability of the mutants in the presence of chemical denaturants, like urea and guanidium chloride, is WT > Cys13Glu > Cys13Asp. Irreversible thermal precipitation temperatures follow the same order as well. Modeling studies establish that the Cys13Asp mutation is likely to cause a significantly greater structural perturbation than Cys13Glu. Analysis of sequence and structural data for TIMs from diverse sources suggests that residues 13 and 82 form a pair of proximal sites, in which a limited number of residue pairs may be accommodated.

The Gaussian Toeplitz matrix

An analytical expression for the LL(T) decomposition for the Gaussian Toeplitz matrix with elements T(ij) = [1/(2-pi)1/2-sigma] exp[-(i - j)2/2-sigma-2] is derived. An exact expression for the determinant and bounds on the eigenvalues follows. An analytical expression for the inverse T-1 is also derived.

Tunable Brownian vortex at the interface

A general kind of Brownian vortices is demonstrated by applying an external nonconservative force field to a colloidal particle bound by a conservative optical trapping force at a liquid-air interface. As the liquid medium is translated at a constant velocity with the bead trapped at the interface, the drag force near the surface provides enough rotational component to bias the particle's thermal fluctuations in a circulatory motion. The interplay between the thermal fluctuations and the advection of the bead in constituting the vortex motions is studied, and we infer that the angular velocity of the circulatory motion offers a comparative measure of the interface fluctuations.

Symmetrizing a Hessenberg matrix: Designs for VLSI parallel processor arrays

A symmetrizer of a nonsymmetric matrix A is the symmetric matrix X that satisfies the equation XA = A(t)X, where t indicates the transpose. A symmetrizer is useful in converting a nonsymmetric eigenvalue problem into a symmetric one which is relatively easy to solve and finds applications in stability problems in control theory and in the study of general matrices. Three designs based on VLSI parallel processor arrays are presented to compute a symmetrizer of a lower Hessenberg matrix. Their scope is discussed. The first one is the Leiserson systolic design while the remaining two, viz., the double pipe design and the fitted diagonal design are the derived versions of the first design with improved performance.