Integrating strong and weak discontinuities without integration subcells and example applications in an XFEM/GFEM framework

Partition of unity methods, such as the extended finite element method, allows discontinuities to be simulated independently of the mesh (Int. J. Numer. Meth. Engng. 1999; 45:601-620). This eliminates the need for the mesh to be aligned with the discontinuity or cumbersome re-meshing, as the discontinuity evolves. However, to compute the stiffness matrix of the elements intersected by the discontinuity, a subdivision of the elements into quadrature subcells aligned with the discontinuity is commonly adopted. In this paper, we use a simple integration technique, proposed for polygonal domains (Int. J. Nuttier Meth. Engng 2009; 80(1):103-134. DOI: 10.1002/nme.2589) to suppress the need for element subdivision. Numerical results presented for a few benchmark problems in the context of linear elastic fracture mechanics and a multi-material problem show that the proposed method yields accurate results. Owing to its simplicity, the proposed integration technique can be easily integrated in any existing code. Copyright (C) 2010 John Wiley & Sons, Ltd.

Bearing capacity of interfering multiple strip footings by using lower bound finite elements limit analysis

The ultimate bearing capacity of a number of multiple strip footings, identically spaced and equally loaded to failure at the same time,is computed by using the lower bound limit analysis in combination with finite elements. The efficiency factor due to the component of soil unit weight, is computed with respect to changes in the clear spacing (xi(gamma)) between the footings. It is noted that the failure load for a footing in the group becomes always greater than that of a single isolated footing. The values of xi(gamma) for the smooth footings are found to be always lower than the rough footings. The values ofxi(gamma) are found to increase continuously with a decrease in the spacing between footings. As compared to the available theoretical and experimental results reported in literature, the present analysis provides generally a little lower values of xi(gamma). (C) 2010 Elsevier Ltd. All rights reserved.

A density matrix renormalization group study of low-lying excitations of polythiophene within a Pariser-Parr-Pople model

Symmetrized density-matrix-renormalization-group calculations have been carried out, within Pariser-Parr-Pople Hamiltonian, to explore the nature of the ground and low-lying excited states of long polythiophene oligomers. We have exploited C-2 symmetry and spin parity of the system to obtain excited states of experimental interest, and studied the lowest dipole allowed excited state and lowest dipole forbidden two photon state, for different oligomer sizes. In the long system limit, the dipole allowed excited state always lies below the lowest dipole forbidden two-photon state which implies, by Kasha rule, that polythiophene fluoresces strongly. The lowest triplet state lies below two-photon state as usual in conjugated polymers. We have doped the system with a hole and an electron and obtained the charge excitation gap and the binding energy of the 1(1)B(u)(-) exciton. We have calculated the charge density of the ground, one-photon and two-photon states for the longer system size of 10 thiophene rings to characterize these states. We have studied bond order in these states to get an idea about the equilibrium excited state geometry of the system. We have also studied the charge density distribution of the singly and doubly doped polarons for longer system size, and observe that polythiophenes do not support bipolarons.

A soluble random-matrix model for relaxation in quantum systems

We study the relaxation of a degenerate two-level system interacting with a heat bath, assuming a random-matrix model for the system-bath interaction. For times larger than the duration of a collision and smaller than the Poincaré recurrence time, the survival probability of still finding the system at timet in the same state in which it was prepared att=0 is exactly calculated.

On the hardness and elastic modulus of bulk metallic glass matrix composites

The influences of the amorphous matrix and crystalline dendrite phases on the hardness and elastic moduli of Zr/Ti-based bulk metallic glass matrix composites have been assessed. While the moduli of the composites correspond to those predicted by the rule of mixtures, the hardness of the composites is similar to that of the matrix, suggesting that the plastic flow in the composites under constrained conditions such as indentation is controlled by the flow resistance of the contiguous matrix. (C) 2010 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Computing a matrix symmetrizer exactly using modified multiple modulus residue arithmetic

A symmetric solution X satisfying the matrix equation XA = AtX is called a symmetrizer of the matrix A. A general algorithm to compute a matrix symmetrizer is obtained. A new multiple-modulus residue arithmetic called floating-point modular arithmetic is described and implemented on the algorithm to compute an error-free matrix symmetrizer.

PVA-SSA-HPA Mixed-Matrix-Membrane Electrolytes for DMFCs

Stabilized forms of heteropolyacids (HPAs), namely phosphomolybdic acid (PMA), phosphotungstic acid (PTA), and silicotungstic acid (STA), are incorporated into poly (vinyl alcohol) (PVA) cross-linked with sulfosuccinic acid (SSA) to form mixed-matrix membranes for application in direct methanol fuel cells (DMFCs). Bridging SSA between PVA molecules not only strengthens the network but also facilitates proton conduction in HPAs. The mixed-matrix membranes are characterized for their mechanical stability, sorption capability, ion-exchange capacity, and wetting in conjunction with their proton conductivity, methanol permeability, and DMFC performance. Methanol-release kinetics is studied ex situ by volume-localized NMR spectroscopy (employing point-resolved spectroscopy'') with the results clearly demonstrating that the incorporation of certain inorganic fillers in PVA-SSA viz., STA and PTA, retards the methanol-release kinetics under osmotic drag compared to Nafion, although PVA-SSA itself exhibits a still lower methanol permeability. The methanol crossover rate for PVA-SSA-HPA-bridged-mixed-matrix membranes decreases dramatically with increasing current density rendering higher DMFC performance in relation to a DMFC using a pristine PVA-SSA membrane. A peak power density of 150 mW/cm(2) at a load current density of 500 mA/cm(2) is achieved for the DMFC using a PVA-SSA-STA-bridged-mixed-matrix-membrane electrolyte. (C) 2010 The Electrochemical Society. [DOI: 10.1149/1.3465653] All rights reserved.

Elasto-plastic Cracking Analysis of Reinforced Concrete

A new elasto-plastic cracking constitutive model for reinforced concrete is presented. The nonlinear effects considered cover almost all the nonlinearities exhibited by reinforced concrete under short term monotonic loading. They include concrete cracking in tension, plasticity in compression, aggregate interlock, tension softening, elasto-plastic behavior of steel, bond-slip between concrete, and steel reinforcement and tension stiffening. A new procedure for incorporating bondslip in smeared steel elements is described. A modified Huber-Hencky-Mises failure criterion for plastic deformation of concrete, which fits the experimental results under biaxial stresses better, is proposed. Multiple cracking at Gauss points and their opening and closing are considered. Matrix expressions are developed and are incorporated in a nonlinear finite element program. After the objectivity of the model is demonstrated, the model is used to analyze two different types of problems: one, a set of four shear panels, and the other, a reinforced concrete beam without shear reinforcement. The results of the analysis agree favorably with the experimental results.

Partially conserved axial-vector current inS-matrix theory

Mandelstam�s argument that PCAC follows from assigning Lorentz quantum numberM=1 to the massless pion is examined in the context of multiparticle dual resonance model. We construct a factorisable dual model for pions which is formulated operatorially on the harmonic oscillator Fock space along the lines of Neveu-Schwarz model. The model has bothm ? andm ? as arbitrary parameters unconstrained by the duality requirement. Adler self-consistency condition is satisfied if and only if the conditionm?2?m?2=1/2 is imposed, in which case the model reduces to the chiral dual pion model of Neveu and Thorn, and Schwarz. The Lorentz quantum number of the pion in the dual model is shown to beM=0.

A rectangular laminated anisotropic shallow thin shell finite element

This report contains the details of the development of the stiffness matrix for a rectangular laminated anisotropic shallow thin shell finite element. The derivation is done under linear thin shell assumptions. Expressing the assumed displacement state over the middle surface of the shell as products of one-dimensional first-order Hermite interpolation polynomials, it is possible to insure that the displacement state for the assembled set of such elements, to be geometrically admissible. Monotonic convergence of the total potential energy is therefore possible as the modelling is successively refined. The element is systematically evaluated for its performance considering various examples for which analytical or other solutions are available

On computing an equivalent symmetric matrix for a nonsymmetric matrix

A real or a complex symmetric matrix is defined here as an equivalent symmetric matrix for a real nonsymmetric matrix if both have the same eigenvalues. An equivalent symmetric matrix is useful in computing the eigenvalues of a real nonsymmetric matrix. A procedure to compute equivalent symmetric matrices and its mathematical foundation are presented.

An integer arithmetic method to compute generalized matrix inverse and solve linear equations exactly

An algorithm that uses integer arithmetic is suggested. It transforms anm ×n matrix to a diagonal form (of the structure of Smith Normal Form). Then it computes a reflexive generalized inverse of the matrix exactly and hence solves a system of linear equations error-free.

Functional Identification of Regulatory Elements within the Promoter Region of Platelet-Derived Growth Factor 2

Human platelet-derived growth factor (PDGF) is composed of two polypeptide chains, PDGF-1 and PDGF-2,the human homolog of the v-sis oncogene. Deregulation of PDGF-2 expression can confer a growth advantage to cells possessing the cognate receptor and, thus, may contribute to the malignant phenotype. We investigated the regulation of PDGF-2 mRNA expression during megakaryocytic differentiation of K562 cells. Induction by 12-O-tetradecanoylphorbol-13-acetate (TPA) led to a greater than 200-fold increase in PDGF-2 transcript levels in these cells. Induction was dependent on protein synthesis and was not enhanced by cycloheximide exposure.In our initial investigation of the PDGF-2 promoter, a minimal promoter region, which included sequences extending only 42 base pairs upstream of the TATA signal, was found to be as efficient as 4 kilobase pairs upstream of the TATA signal in driving expression of a reporter gene in uninduced K562 cells. We also functionally identified different regulatory sequence elements of the PDGF-2 promoter in TPA-induced K562 cells. One region acted as a transcriptional silencer, while another region was necessary for maximal activity of the promoter in megakaryoblasts. This region was shown to bind nuclear factors and was the target for trans-activation in normal and tumor cells. In one tumor cell line, which expressed high PDGF-2 mRNA levels, the presence of the positive regulatory region resulted in a 30-fold increase in promoter activity. However, the ability of the minimal PDGF-2 promoter to drive reporter gene expression in uninduced K562 cells and normal fibroblasts, which contained no detectable PDGF-2 transcripts, implies the existence of other negative control mechanisms beyond the regulation of promoter activity.