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.

Wave propagation analysis in carbon nanotube embedded composite using wavelet based spectral finite elements

In this paper, elastic wave propagation is studied in a nanocomposite reinforced with multiwall carbon nanotubes (CNTs). Analysis is performed on a representative volume element of square cross section. The frequency content of the exciting signal is at the terahertz level. Here, the composite is modeled as a higher order shear deformable beam using layerwise theory, to account for partial shear stress transfer between the CNTs and the matrix. The walls of the multiwall CNTs are considered to be connected throughout their length by distributed springs, whose stiffness is governed by the van der Waals force acting between the walls of nanotubes. The analyses in both the frequency and time domains are done using the wavelet-based spectral finite element method (WSFEM). The method uses the Daubechies wavelet basis approximation in time to reduce the governing PDE to a set of ODEs. These transformed ODEs are solved using a finite element (FE) technique by deriving an exact interpolating function in the transformed domain to obtain the exact dynamic stiffness matrix. Numerical analyses are performed to study the spectrum and dispersion relations for different matrix materials and also for different beam models. The effects of partial shear stress transfer between CNTs and matrix on the frequency response function (FRF) and the time response due to broadband impulse loading are investigated for different matrix materials. The simultaneous existence of four coupled propagating modes in a double-walled CNT-composite is also captured using modulated sinusoidal excitation.

Groove Binding Ligands for the Interaction with Parallel-Stranded ps-Duplex DNA and Triplex DNA

DNA adopts different conformations not only based on novel base pairs, but also with different chain polarities. Besides several duplex structures (A, B, Z, parallel stranded (ps)-DNA, etc.), DNA also forms higher-order structures like triplex, tetraplex, and i-motif. Each of these structures has its own biological significance. The ps-duplexes have been found to be resistant to certain nucleases and endonucleases. Molecules that promote triple-helix formation have significant potential. These investigations have many therapeutic advantages which may be useful in the regulation of the expression of genes responsible for certain diseases by locking either their transcription (antigene) or translation (antisense). Each DNA minor groove binding ligand (MGBL) interacts with DNA through helical minor groove recognition in a sequence-specific manner, and this interferes with several DNA-associated processes. Incidentally, these ligands interact with some non-B-DNA and with higher-order DNA structures including ps-DNA and triplexes. While the design and recognition of minor grooves of duplex DNA by specific MGBLs have been a topic of many reports, limited information is available on the binding behavior of MGBLs with nonduplex DNA. In this review, we summarize various attempts of the interaction of MGBLs with ps-DNA and DNA triplexes.

A systolic evaluator for linear, quadratic, and cubic expressions

High-speed evaluation of a large number of linear, quadratic, and cubic expressions is very important for the modeling and real-time display of objects in computer graphics. Using VLSI techniques, chips called pixel planes have actually been built by H. Fuchs and his group to evaluate linear expressions. In this paper, we describe a topological variant of Fuchs' pixel planes which can evaluate linear, quadratic, cubic, and higher-order polynomials. In our design, we make use of local interconnections only, i.e., interconnections between neighboring processing cells. This leads to the concept of tiling the processing cells for VLSI implementation.

Noise minimization in RC-active filters using generalized impedance converters

A procedure has been given for minimizing the total output noise of a Generalized Impedance Converter (GIC), subject to constraints dictated by signal handling capability of the Operational Amplifiers and ease of microcircuit fabrication. The noise reduction is achieved only by the adjustment of RC elements of the GIC, and the total output noise after optimization in the example cited is close to the theoretical lower limit. The output noise of a higher-order filter can be reduced by RC-optimizing the individual GIC's of the active realization. Experimental results on a 20–24 kHz channel bank band-pass filter demonstrate the effectiveness of the above procedure.

Thermodynamics of Cobalt (II, III) Oxide (Co3O4): Evidence of Phase Transition

The Gibbs' energy change for the reaction, 3CoO (r.s.)+1/2O2(g)→Co3O4(sp), has been measured between 730 and 1250 K using a solid state galvanic cell: Pt, CuO+Cu2O|(CaO)ZrO2|CoO+Co3O4,Pt. The emf of this cell varies nonlinearly with temperature between 1075 and 1150 K, indicating a second or higher order phase transition in Co3O4around 1120 (±20) K, associated with an entropy change of ∼43 Jmol-1K-1. The phase transition is accompanied by an anomalous increase in lattice parameter and electrical conductivity. The cubic spinel structure is retained during the transition, which is caused by the change in CO+3 ions from low spin to high spin state. The octahedral site preference energy of CO+3 ion in the high spin state has been evaluated as -24.8 kJ mol-1. This is more positive than the value for CO+2 ion (-32.9 kJ mol-1). The cation distribution therefore changes from normal to inverse side during the phase transition. The transformation is unique, coupling spin unpairing in CO+3 ion with cation rearrangement on the spinel lattice, DTA in pure oxygen revealed a small peak corresponding to the transition, which could be differentiated from the large peak due to decomposition. TGA showed that the stoichiometry of oxide is not significantly altered during the transition. The Gibbs' energy of formation of Co3O4 from CoO and O2 below and above phase transition can be represented by the equations:ΔG0=-205,685+170.79T(±200) J mol-1(730-1080 K) and ΔG0=-157,235+127.53T(±200) J mol-1(1150-1250 K).

Structure and representation of correlation functions and the density matrix for a statistical wave field in optics

A systematic structure analysis of the correlation functions of statistical quantum optics is carried out. From a suitably defined auxiliary two‐point function we are able to identify the excited modes in the wave field. The relative simplicity of the higher order correlation functions emerge as a byproduct and the conditions under which these are made pure are derived. These results depend in a crucial manner on the notion of coherence indices and of unimodular coherence indices. A new class of approximate expressions for the density operator of a statistical wave field is worked out based on discrete characteristic sets. These are even more economical than the diagonal coherent state representations. An appreciation of the subtleties of quantum theory obtains. Certain implications for the physics of light beams are cited.

Filter Optimization for Grid Interactive Voltage Source Inverters

Higher order LCL filters are essential in meeting the interconnection standard requirement for grid-connected voltage source converters. LCL filters offer better harmonic attenuation and better efficiency at a smaller size when compared to the traditional L filters. The focus of this paper is to analyze the LCL filter design procedure from the point of view of power loss and efficiency. The IEEE 1547-2008 specifications for high-frequency current ripple are used as a major constraint early in the design to ensure that all subsequent optimizations are still compliant with the standards. Power loss in each individual filter component is calculated on a per-phase basis. The total inductance per unit of the LCL filter is varied, and LCL parameter values which give the highest efficiency while simultaneously meeting the stringent standard requirements are identified. The power loss and harmonic output spectrum of the grid-connected LCL filter is experimentally verified, and measurements confirm the predicted trends.

A Smooth Finite Element Method Based on Reproducing Kernel DMS-Splines

The element-based piecewise smooth functional approximation in the conventional finite element method (FEM) results in discontinuous first and higher order derivatives across element boundaries Despite the significant advantages of the FEM in modelling complicated geometries, a motivation in developing mesh-free methods has been the ease with which higher order globally smooth shape functions can be derived via the reproduction of polynomials There is thus a case for combining these advantages in a so-called hybrid scheme or a `smooth FEM' that, whilst retaining the popular mesh-based discretization, obtains shape functions with uniform C-p (p >= 1) continuity One such recent attempt, a NURBS based parametric bridging method (Shaw et al 2008b), uses polynomial reproducing, tensor-product non-uniform rational B-splines (NURBS) over a typical FE mesh and relies upon a (possibly piecewise) bijective geometric map between the physical domain and a rectangular (cuboidal) parametric domain The present work aims at a significant extension and improvement of this concept by replacing NURBS with DMS-splines (say, of degree n > 0) that are defined over triangles and provide Cn-1 continuity across the triangle edges This relieves the need for a geometric map that could precipitate ill-conditioning of the discretized equations Delaunay triangulation is used to discretize the physical domain and shape functions are constructed via the polynomial reproduction condition, which quite remarkably relieves the solution of its sensitive dependence on the selected knotsets Derivatives of shape functions are also constructed based on the principle of reproduction of derivatives of polynomials (Shaw and Roy 2008a) Within the present scheme, the triangles also serve as background integration cells in weak formulations thereby overcoming non-conformability issues Numerical examples involving the evaluation of derivatives of targeted functions up to the fourth order and applications of the method to a few boundary value problems of general interest in solid mechanics over (non-simply connected) bounded domains in 2D are presented towards the end of the paper

Effect of wall thickness on the end corrections of the extended inlet and outlet of a double-tuned expansion chamber

Simple expansion chambers, the simplest of the muffler configurations, have very limited practical application due to the presence of periodic troughs in the transmission loss spectrum which drastically lower the overall transmission loss of the muffler. Tuned extended inlet and outlet can be designed to nullify three-fourths of these troughs, making use of the plane wave theory. These cancellations would not occur unless one altered the geometric lengths for the extended tube in order to incorporate the effect of evanescent higher-order modes (multidimensional effect) through end corrections or lumped inertance approximation at the area discontinuities or junctions. End corrections of the extended inlet and outlet have been studied by several researchers. However the effect of wall thickness of the inlet/outlet duct on end correction has not been studied explicitly. This has significant effect on the tuning of an extended inlet/outlet expansion chamber. It is investigated here experimentally as well as numerically (through use of 3-D FEM software) for stationary medium. Crown Copyright (C) 2010 Published by Elsevier Ltd. All rights reserved.

Q-Band beam wave-guide

Consideration is given to a 25-foot long Q-band (8 mm) confocal, zoned dielectric lens beam waveguide. Numerical expressions for the axial and radial fields are presented. The experimental set-up consisted of uniformly spaced zoned dielectric lenses, a transmitting horn and a receiving horn. It was found that: (1) the wave beam is reiterated when confocal, zoned dielectric lenses act as phase transformers in place of smooth surfaced transformers in beam waveguides; (2) the axial field is oscillatory near the source and the oscillation persists for about 25 cm from the source; (3) the oscillation disappears after one lens is used; (4) higher order modes with higher attenuation rates die out faster than fundamental modes; (5) phase transformers do not alter beam modes; (6) without any lens the beam cross-section broadens significantly in the Z-direction; (7) with one lens the beam exhibits the reiteration phenomenon; and (8) inserting a second lens on the axial and cross-sectional field distribution shows further the reiteration principle.

Random-Restart Reactive Tabu Search Algorithm for Detection in Large-MIMO Systems

We present a low-complexity algorithm based on reactive tabu search (RTS) for near maximum likelihood (ML) detection in large-MIMO systems. The conventional RTS algorithm achieves near-ML performance for 4-QAM in large-MIMO systems. But its performance for higher-order QAM is far from ML performance. Here, we propose a random-restart RTS (R3TS) algorithm which achieves significantly better bit error rate (BER) performance compared to that of the conventional RTS algorithm in higher-order QAM. The key idea is to run multiple tabu searches, each search starting with a random initial vector and choosing the best among the resulting solution vectors. A criterion to limit the number of searches is also proposed. Computer simulations show that the R3TS algorithm achieves almost the ML performance in 16 x 16 V-BLAST MIMO system with 16-QAM and 64-QAM at significantly less complexities than the sphere decoder. Also, in a 32 x 32 V-BLAST MIMO system, the R3TS performs close to ML lower bound within 1.6 dB for 16-QAM (128 bps/Hz), and within 2.4 dB for 64-QAM (192 bps/Hz) at 10(-3) BER.

A simple method for potential flow simulation of cascades

A simple method using a combination of conformal mapping and vortex panel method to simulate potential flow in cascades is presented. The cascade is first transformed to a single body using a conformal mapping, and the potential flow over this body is solved using a simple higher order vortex panel method. The advantage of this method over existing methodologies is that it enables the use of higher order panel methods, as are used to solve flow past an isolated airfoil, to solve the cascade problem without the need for any numerical integrations or iterations. The fluid loading on the blades, such as the normal force and pitching moment, may be easily calculated from the resultant velocity field. The coefficient of pressure on cascade blades calculated with this methodology shows good agreement with previous numerical and experimental results.

Unveiling the Role of Dps in the Organization of Mycobacterial Nucleoid

In order to preserve genetic information in stress conditions, bacterial DNA is organized into higher order nucleoid structure. In this paper, with the help of Atomic Force Microscopy, we show the different structural changes in mycobacterial nucleoid at different points of growth in the presence of different concentrations of glucose in the medium. We also observe that in Mycobacterium smegmatis, two different Dps proteins (Dps1 and Dps2) promote two types of nucleoid organizations. At the late stationary phase, under low glucose availability, Dps1 binds to DNA to form a very stable toroid structure. On the other hand, under the same condition, Dps2-DNA complex forms an incompletely condensed toroid and finally forms a further stable coral reef structure in the presence of RNA. This coral reef structure is stable in high concentration of bivalent ion like Mg2+.

Extended gravitational action and novel consequences

It is found that the inclusion of higher derivative terms in the gravitational action along with concepts of phase transition and spontaneous symmetry breaking leads to some novel consequence. The Ricci scalar plays the dual role, like a physical field as well as a geometrical field. One gets Klein-Gordon equation for the emerging field and the corresponding quanta of geometry are called Riccions. For the early universe the model removes singularity along with inflation. In higher dimensional gravity the Riccions can break into spin half particle and antiparticle along with breaking of left-right symmetry. Most tantalizing consequences is the emergence of the physical universe from the geometry in the extreme past. Riccions can Bose condense and may account for the dark matter.