Extending lagrange interpolation to develop a 4-node quadrilateral element

Today finite element method is a well established tool in engineering analysis and design. Though there axe many two and three dimensional finite elements available, it is rare that a single element performs satisfactorily in majority of practical problems. The present work deals with the development of 4-node quadrilateral element using extended Lagrange interpolation functions. The classical univariate Lagrange interpolation is well developed for 1-D and is used for obtaining shape functions. We propose a new approach to extend the Lagrange interpolation to several variables. When variables axe more than one the method also gives the set of feasible bubble functions. We use the two to generate shape function for the 4-node arbitrary quadrilateral. It will require the incorporation of the condition of rigid body motion, constant strain and Navier equation by imposing necessary constraints. The procedure obviates the need for isoparametric transformation since interpolation functions are generated for arbitrary quadrilateral shapes. While generating the element stiffness matrix, integration can be carried out to the accuracy desired by dividing the quadrilateral into triangles. To validate the performance of the element which we call EXLQUAD4, we conduct several pathological tests available in the literature. EXLQUAD4 predicts both stresses and displacements accurately at every point in the element in all the constant stress fields. In tests involving higher order stress fields the element is assured to converge in the limit of discretisation. A method thus becomes available to generate shape functions directly for arbitrary quadrilateral. The method is applicable also for hexahedra. The approach should find use for development of finite elements for use with other field equations also.

Topology design of three-dimensional structures using hybrid finite elements

Conventional three-dimensional isoparametric elements are susceptible to problems of locking when used to model plate/shell geometries or when the meshes are distorted etc. Hybrid elements that are based on a two-field variational formulation are immune to most of these problems, and hence can be used to efficiently model both "chunky" three-dimensional and plate/shell type structures. Thus, only one type of element can be used to model "all" types of structures, and also allows us to use a standard dual algorithm for carrying out the topology optimization of the structure. We also address the issue of manufacturability of the designs.

Activities in liquid Ga-Te alloys at 1120 k

The activity of gallium in liquid Ga-Te alloys has been measured at 1120 K using a solid state galvanic cell incorporating yttria-stabilized thoria as the solid electrolyte. The cell can be schematically represented as (−) W,Re,Ga(1)+Ga2O3(s)|(Y2O3) ThO2|Ga-Te(1) + Ga2O3(s), Re, W (+) The activity of tellurium was derived by Gibbs-Duhem integration. The activity of gallium shows negative deviation from Raoult's law for XGa < 0.6 and positive deviation from ideality for XGa > 0.6. The activity of gallium was constant in the composition range 0.73 < XGa < 0.89, indicating liquid state immiscibility in this region. The Gibbs energy of mixing and the concentration-concentration structure factor at long wavelength limit show a minimum at XGa ≈ 0.4, suggesting strong interactions in the liquid phase with formation of ‘Ga2Te3‘-type complexes

One-dimensional fermions with incommensuration close to dimerization

We study a model of fermions hopping on a chain with a weak incommensuration close to dimerization; both q, the deviation of the wave number from pi, and delta, the strength of the incommensuration, are assumed to be small. For free fermions, we show that there are an infinite number of energy bands which meet at zero energy as q approaches zero. The number of states lying inside the q = 0 gap remains nonzero as q/delta --> 0. Thus the limit q --> 0 differs from q = 0, as can be seen clearly in the low-temperature specific heat. For interacting fermions or the XXZ spin-(1/2) chain, we use bosonization to argue that similar results hold. Finally, our results can be applied to the Azbel-Hofstadter problem of particles hopping on a two-dimensional lattice in the presence of a magnetic field.

Structure and dynamics of benzene in one-dimensional channels

Molecular dynamics investigation of benzene in one-dimensional channel systems A1PO(4)-5, VPI-5, and carbon nanotube is reported. The results suggest that, in all the three host systems, the plane of benzene is almost perpendicular to the channel axis when the molecule is near the center of the channel and the plane of benzene is parallel to the channel axis when the molecule is near the wall of the channel. The density distribution of benzene as a function of channel length, z and the radial distance, r, from the channel axis is also different in the three host structures. Anisotropy in translational diffusion coefficient, calculated in body-fixed frame of benzene, suggests that benzene prefers to move with its plane parallel to the direction of motion in A1PO(4)-5 and VPI-5 whereas in carbon nanotube the motion occurs predominantly with the plane of the benzene perpendicular to the direction of motion.;Anisotropy associated with the rotational motion is seen to alter significantly in confinement as compared to liquid benzene. In A1PO(4)-5, the rotational anisotropy is reversed as compared to liquid benzene thereby suggesting that anisotropy arising out of molecular geometry can be reduced. Reorientational correlation times for C-6 and C-2 axes Of benzene are reported. Apart from the inertial decay of reorientational correlation function due to free, rotation, two other distinct regimes of decay are observed in narrower channels (AIPO(4)-5 and carbon nanotube): (i) an initial fast decay (0.5-2 ps) and (ii) a slower decay (>2 ps) of reorientational correlation function where C-6 decays slower than C-2 Similar to what is observed in liquid benzene. In the initial fast decay, it is seen that the decay for C-6 is faster than C-2 which is in contrast to what is observed in liquid benzene or for benzene confined in VPI-5.

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.

Double diffusive unsteady free convection on two-dimensional and axisymmetric bodies in a porous medium

The unsteady laminar free convection flow of an incompressible electrically conducting fluid over two-dimensional and axisymmetric bodies embedded in a highly porous medium with an applied magnetic field has been studied. The unsteadiness in the flow field is caused by the variation of the wall temperature and concentration with time. The coupled nonlinear partial differential equations with three independent variables have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. It is observed that the skin friction, heat transfer and mass transfer increase with the permeability parameter but decrease with the magnetic parameter. The results are strongly dependent on the variation of wall temperature and concentration with time. The skin friction and heat transfer increase or decrease as the buoyancy forces from species diffusion assist or oppose the thermal buoyancy force. However, the mass transfer is found to be higher for small values of the ratio of the buoyancy parameters than for large values

Nonclassical rotational inertia in a two-dimensional bosonic solid containing grain boundaries

We study the occurrence of nonclassical rotational inertia (NCRI) arising from superfluidity along grain boundaries in a two-dimensionalbosonic system. We make use of a standard mapping between the zero-temperature properties of this system and the statistical mechanics of interacting vortex lines in the mixed phase of a type-II superconductor. In the mapping, the liquid phase of the vortex system corresponds to the superfluid bosonic phase. We consider numerically obtained polycrystalline configurations of the vortex lines in which the microcrystals are separated by liquidlike grain-boundary regions which widen as the vortex system temperature increases. The NCRI of the corresponding zero-temperature bosonic systems can then be numerically evaluated by solving the equations of superfluid hydrodynamics in the channels near the grain boundaries. We find that the NCRI increases very abruptly as the liquid regions in the vortex system (equivalently, superfluid regions in the bosonic system) form a connected, system-spanning structure with one or more closed loops. The implications of these results for experimentally observed supersolid phenomena are discussed.

Geometry changes induced by negative hyperconjugative interactions involving carbonyl and thiocarbonyl groups

MNDO geometry optimizations have been carried out on a series of acyclic and cyclic unsymmetrically disubstituted carbonyl and thiocarbonyl compounds. The C=X unit shows a consistent and often sizeable tilt towards one of the substituents, following the order O > Snot, vert, similarN > C > B. Reference ab initio calculations and available experimental results support the MNDO results. The effect, which is particularly dramatic in small rings, is attributed primarily to favorable negative hyperconjugative interaction between the lone pair on X and a low lying adjacent σ* orbital. Such an interaction can lead to highly distorted structures, including perhaps to a planar molecule with an inverted sp2 carbon center.

Three-Dimensional Unsteady Flow With Heat and Mass Transfer Over a Continuous Stretching Surface

A numerical solution of the unsteady boundary layer equations under similarity assumptions is obtained. The solution represents the three-dimensional unsteady fluid motion caused by the time-dependent stretching of a flat boundary. It has been shown that a self-similar solution exists when either the rate of stretching is decreasing with time or it is constant. Three different numerical techniques are applied and a comparison is made among them as well as with earlier results. Analysis is made for various situations like deceleration in stretching of the boundary, mass transfer at the surface, saddle and nodal point flows, and the effect of a magnetic field. Both the constant temperature and constant heat flux conditions at the wall have been studied.