New twelve node serendipity quadrilateral plate bending element based on Mindlin-Reissner theory using Integrated Force Method

The Integrated Force Method (IFM) is a novel matrix formulation developed for analyzing the civil, mechanical and aerospace engineering structures. In this method all independent/internal forces are treated as unknown variables which are calculated by simultaneously imposing equations of equilibrium and compatibility conditions. This paper presents a new 12-node serendipity quadrilateral plate bending element MQP12 for the analysis of thin and thick plate problems using IFM. The Mindlin-Reissner plate theory has been employed in the formulation which accounts the effect of shear deformation. The performance of this new element with respect to accuracy and convergence is studied by analyzing many standard benchmark plate bending problems. The results of the new element MQP12 are compared with those of displacement-based 12-node plate bending elements available in the literature. The results are also compared with exact solutions. The new element MQP12 is free from shear locking and performs excellent for both thin and moderately thick plate bending situations.

Exact solution for the unsteady motion of a viscous-fluid in a porous annulus

An exact solution to the problem of time-dependent motion of a viscous fluid in an annulus with porous walls is obtained under the assumption that the rate of suction at one wall is equal to the rate of injection at the other. Finite Hankel transform is used to obtain a closed-form solution for the axial velocity. The average axial velocity profiles are depicted graphically.

Magnetic, electric, and dielectric properties of FeCo alloy nanoparticles dispersed in amorphous matrix

Nanocrystalline Fe53Co47 alloy was synthesized by a single-step transmetallation chemical method at room temperature. The Fe53Co47 alloy nanoparticles of 77 and 47 wt% were dispersed in silica matrix by the sol-gel process using tetraethyl orthosilcate. Structural studies reveal that the as-prepared alloy powders are in bcc phase and silica is in an amorphous state. The phase-transition temperature and Mossbauer spectra analysis of the Fe-Co alloy establishes the homogeneous alloy formation. A saturation magnetization of 218 emu/g was obtained for pure FeCo alloy at room temperature. Scanning electron microscopic analysis demonstrates the hollow-sphere morphology for FeCo alloy particles. Magnetic nanocomposite consisting of 47 wt% FeCo-silica shows enhanced thermal stability over the native FeCo alloy. Electrical and dielectric properties of 47 wt% FeCo-silica nanocomposites were investigated as a function of frequency and temperature. It was found that the dielectric constants and dielectric loss were stable throughout the measured temperature (310-373 K). Our results indicate that FeCo-silica nanocomposite is a promising candidate for high-frequency applications. (C) 2010 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim

An optimal parallel algorithm for sparse matrix bandwidth reduction on a linear array

A simple and efficient algorithm for the bandwidth reduction of sparse symmetric matrices is proposed. It involves column-row permutations and is well-suited to map onto the linear array topology of the SIMD architectures. The efficiency of the algorithm is compared with the other existing algorithms. The interconnectivity and the memory requirement of the linear array are discussed and the complexity of its layout area is derived. The parallel version of the algorithm mapped onto the linear array is then introduced and is explained with the help of an example. The optimality of the parallel algorithm is proved by deriving the time complexities of the algorithm on a single processor and the linear array.

Exact traveling- wave solution for model geophysical systems

Exact traveling-wave solutions of time-dependent nonlinear inhomogeneous PDEs, describing several model systems in geophysical fluid dynamics, are found. The reduced nonlinear ODEs are treated as systems of linear algebraic equations in the derivatives. A variety of solutions are found, depending on the rank of the algebraic systems. The geophysical systems include acoustic gravity waves, inertial waves, and Rossby waves. The solutions describe waves which are, in general, either periodic or monoclinic. The present approach is compared with the earlier one due to Grundland (1974) for finding exact solutions of inhomogeneous systems of nonlinear PDEs.

Exact nonlinear travelling hydromagnetic wave solutions

Exact travelling wave solutions for hydromagnetic waves in an exponentially stratified incompressible medium are obtained. With the help of two integrals it becomes possible to reduce the system of seven nonlinear PDE's to a second order nonlinear ODE which describes an one dimensional harmonic oscillator with a nonlinear friction term. This equation is studied in detail in the phase plane. The travelling waves are periodic only when they propagate either horizontally or vertically. The reduced second order nonlinear differential equation describing the travelling waves in inhomogeneous conducting media has rather ubiquitous nature in that it also appears in other geophysical systems such as internal waves, Rossby waves and topographic Rossby waves in the ocean.

Exact analysis of Burgers equation on semiline with flux condition at the origin

The initial boundary value problem for the Burgers equation in the domain x greater-or-equal, slanted 0, t > 0 with flux boundary condition at x = 0 has been solved exactly. The behaviour of the solution as t tends to infinity is studied and the “asymptotic profile at infinity” is obtained. In addition, the uniqueness of the solution of the initial boundary value problem is proved and its inviscid limit as var epsilon → 0 is obtained.

Exact eigenstates of the Pariser-Parr-Pople model for anthracene

The pi-electronic structure of anthracene is discussed by combining exact solutions of the Pariser-Parr-Pople (PPP) model and semiempirical PM3 calculations. Symmetry adaptation of the 2.8 million singlet valence-bond (VB) diagrams is explicitly demonstrated for D2h and electron-hole symmetry. Standard PPP parameters provide a comprehensive fit to one- and two-photon anthracene spectra and intensities up to the strong 1 B-1(3u)-absorption at 5.24 eV, the 10th excited state in the dense correlated spectrum, and indicate a reassignment of two-photon absorptions. The singlet-triplet gap and fine-structure constants also agree with experiment. Fully-relaxed PM3 geometries are obtained for the anthracene ground state and for singlet, triplet, and charged bipolarons. The PM3 bond lengths correlate well with PPP bond orders for the idealized structure. Single-determinantal PM3 excitation and relaxation energies for bipolarons are consistent with exact PPP results and contrast all-valence electron with pi-electron calculations. Several correlation effects are noted in the rich pi-spectra of anthracene in connection with improved PPP modeling of conjugated molecules and polymers.

A monolithic strategy for fluid-structure interaction problems

In this work, we present a new monolithic strategy for solving fluid-structure interaction problems involving incompressible fluids, within the context of the finite element method. This strategy, similar to the continuum dynamics, conserves certain properties, and thus provides a rational basis for the design of the time-stepping strategy; detailed proofs of the conservation of these properties are provided. The proposed algorithm works with displacement and velocity variables for the structure and fluid, respectively, and introduces no new variables to enforce velocity or traction continuity. Any existing structural dynamics algorithm can be used without change in the proposed method. Use of the exact tangent stiffness matrix ensures that the algorithm converges quadratically within each time step. An analytical solution is presented for one of the benchmark problems used in the literature, namely, the piston problem. A number of benchmark problems including problems involving free surfaces such as sloshing and the breaking dam problem are used to demonstrate the good performance of the proposed method. Copyright (C) 2010 John Wiley & Sons, Ltd.

Processing maps for hot-working of powder metallurgy 1100 Al-10 vol % SiC-particulate metal-matrix composite

The hot-working characteristics of the metal-matrix composite (MMC) Al-10 vol % SiC-particulate (SiCp) powder metallurgy compacts in as-sintered and in hot-extruded conditions were studied using hot compression testing. On the basis of the stress-strain data as a function of temperature and strain rate, processing maps depicting the variation in the efficiency of power dissipation, given by eegr = 2m/(m+1), where m is the strain rate sensitivity of flow stress, have been established and are interpreted on the basis of the dynamic materials model. The as-sintered MMC exhibited a domain of dynamic recrystallization (DRX) with a peak efficiency of about 30% at a temperature of about 500°C and a strain rate of 0.01 s�1. At temperatures below 350°C and in the strain rate range 0.001�0.01 s�1 the MMC exhibited dynamic recovery. The as-sintered MMC was extruded at 500°C using a ram speed of 3 mm s�1 and an extrusion ratio of 10ratio1. A processing map was established on the extruded product, and this map showed that the DRX domain had shifted to lower temperature (450°C) and higher strain rate (1 s�1). The optimum temperature and strain rate combination for powder metallurgy billet conditioning are 500°C and 0.01 s�1, and the secondary metal-working on the extruded product may be done at a higher strain rate of 1 s�1 and a lower temperature of 425°C.

Relationship Between The Impedance Matrix And The Transfer Matrix With Specific Reference To Symmetrical, Reciprocal And Conservative Systems

Impedance matrix and transfer matrix methods are often used in the analysis of linear dynamical systems. In this paper, general relationships between these matrices are derived. The properties of the impedance matrix and the transfer matrix of symmetrical systems, reciprocal systems and conservative systems are investigated. In the process, the following observations are made: (a) symmetrical systems are not a subset of reciprocal systems, as is often misunderstood; (b) the cascading of reciprocal systems again results in a reciprocal system, whereas cascading of symmetrical systems does not necessarily result in a symmetrical system; (c) the determinant of the transfer matrix, being ±1, is a property of both symmetrical systems and reciprocal systems, but this condition, however, is not sufficient to establish either the reciprocity or the symmetry of the system; (d) the impedance matrix of a conservative system is skew-Hermitian.

Some new results regarding spikes and a heuristic for spike construction

his paper addresses the problem of minimizing the number of columns with superdiagonal nonzeroes (viz., spiked columns) in a square, nonsingular linear system of equations which is to be solved by Gaussian elimination. The exact focus is on a class of min-spike heuristics in which the rows and columns of the coefficient matrix are first permuted to block lower-triangular form. Subsequently, the number of spiked columns in each irreducible block and their heights above the diagonal are minimized heuristically. We show that ifevery column in an irreducible block has exactly two nonzeroes, i.e., is a doubleton, then there is exactly one spiked column. Further, if there is at least one non-doubleton column, there isalways an optimal permutation of rows and columns under whichnone of the doubleton columns are spiked. An analysis of a few benchmark linear programs suggests that singleton and doubleton columns can abound in practice. Hence, it appears that the results of this paper can be practically useful. In the rest of the paper, we develop a polynomial-time min-spike heuristic based on the above results and on a graph-theoretic interpretation of doubleton columns.

The superheating of Pb embedded in a Zn matrix the role of interface melting

Evidence of a shape-dependent superheating of entrained nanosized Pb particles in a Zn matrix has been presented. It is shown that size dependence and pressure effects cannot explain the observed differences in melting points. The importance of crystallography and morphology at the microlevel at the interphase interface in controlling interfacial melting has been emphasized in order to explain the melting of entrained particles.

New Approximation Algorithms for Minimum Cycle Bases of Graphs

We consider the problem of computing an approximate minimum cycle basis of an undirected non-negative edge-weighted graph G with m edges and n vertices; the extension to directed graphs is also discussed. In this problem, a {0,1} incidence vector is associated with each cycle and the vector space over F-2 generated by these vectors is the cycle space of G. A set of cycles is called a cycle basis of G if it forms a basis for its cycle space. A cycle basis where the sum of the weights of the cycles is minimum is called a minimum cycle basis of G. Cycle bases of low weight are useful in a number of contexts, e.g. the analysis of electrical networks, structural engineering, chemistry, and surface reconstruction. Although in most such applications any cycle basis can be used, a low weight cycle basis often translates to better performance and/or numerical stability. Despite the fact that the problem can be solved exactly in polynomial time, we design approximation algorithms since the performance of the exact algorithms may be too expensive for some practical applications. We present two new algorithms to compute an approximate minimum cycle basis. For any integer k >= 1, we give (2k - 1)-approximation algorithms with expected running time O(kmn(1+2/k) + mn((1+1/k)(omega-1))) and deterministic running time O(n(3+2/k) ), respectively. Here omega is the best exponent of matrix multiplication. It is presently known that omega < 2.376. Both algorithms are o(m(omega)) for dense graphs. This is the first time that any algorithm which computes sparse cycle bases with a guarantee drops below the Theta(m(omega) ) bound. We also present a 2-approximation algorithm with expected running time O(M-omega root n log n), a linear time 2-approximation algorithm for planar graphs and an O(n(3)) time 2.42-approximation algorithm for the complete Euclidean graph in the plane.