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.

Interfacial microrheology as a tool to study viscoelastic transitions in nanoconfined soft matter

We present a method to perform in situ microrheological measurements on monolayers of soft materials undergoing viscoelastic transitions under compression. Using the combination of a Langmuir trough mounted on the inverted microscope stage of a laser scanning confocal microscope we track the motion of individual fluorescent quantum dots partly dispersed in monolayers spread at the air-water interface. From the calculated mean square displacement of the probe particles and extending a well established scheme of the generalized Stokes-Einstein relation in bulk to the interface we arrive at the viscoelastic modulus for the respective monolayers as a function of surface density. Measurements on monolayers of glassy as well as nonglassy polymers and a standard fatty acid clearly show sensitivity of our technique to subtle variations, in the viscoelastic properties of the highly confined materials under compression. Evidence for possible spatial variations of such viscoelastic properties at a given surface density for the fatty acid monolayer is also provided.

Molecular theory of solvation and solvation dynamics of a classical ion in a dipolar liquid

A microscopic theory of equilibrium solvation and solvation dynamics of a classical, polar, solute molecule in dipolar solvent is presented. Density functional theory is used to explicitly calculate the polarization structure around a solvated ion. The calculated solvent polarization structure is different from the continuum model prediction in several respects. The value of the polarization at the surface of the ion is less than the continuum value. The solvent polarization also exhibits small oscillations in space near the ion. We show that, under certain approximations, our linear equilibrium theory reduces to the nonlocal electrostatic theory, with the dielectric function (c(k)) of the liquid now wave vector (k) dependent. It is further shown that the nonlocal electrostatic estimate of solvation energy, with a microscopic c(k), is close to the estimate of linearized equilibrium theories of polar liquids. The study of solvation dynamics is based on a generalized Smoluchowski equation with a mean-field force term to take into account the effects of intermolecular interactions. This study incorporates the local distortion of the solvent structure near the ion and also the effects of the translational modes of the solvent molecules.The latter contribution, if significant, can considerably accelerate the relaxation of solvent polarization and can even give rise to a long time decay that agrees with the continuum model prediction. The significance of these results is discussed.

A molecular theory of collective orientational relaxation in pure and binary dipolar liquids

A molecular theory of collective orientational relaxation of dipolar molecules in a dense liquid is presented. Our work is based on a generalized, nonlinear, Smoluchowski equation (GSE) that includes the effects of intermolecular interactions through a mean‐field force term. The effects of translational motion of the liquid molecules on the orientational relaxation is also included self‐consistently in the GSE. Analytic expressions for the wave‐vector‐dependent orientational correlation functions are obtained for one component, pure liquid and also for binary mixtures. We find that for a dipolar liquid of spherical molecules, the correlation function ϕ(k,t) for l=1, where l is the rank of the spherical harmonics, is biexponential. At zero wave‐vector, one time constant becomes identical with the dielectric relaxation time of the polar liquid. The second time constant is the longitudinal relaxation time, but the contribution of this second component is small. We find that polar forces do not affect the higher order correlation functions (l>1) of spherical dipolar molecules in a linearized theory. The expression of ϕ(k,t) for a binary liquid is a sum of four exponential terms. We also find that the wave‐vector‐dependent relaxation times depend strongly on the microscopic structure of the dense liquid. At intermediate wave vectors, the translational diffusion greatly accelerates the rate of orientational relaxation. The present study indicates that one must pay proper attention to the microscopic structure of the liquid while treating the translational effects. An analysis of the nonlinear terms of the GSE is also presented. An interesting coupling between the number density fluctuation and the orientational fluctuation is uncovered.

Polarization Caging in Diffusion-Controlled Electron Transfer Reactions in Solution

In some bimolecular diffusion-controlled electron transfer (ET) reactions such as ion recombination (IR), both solvent polarization relaxation and the mutual diffusion of the reacting ion pair may determine the rate and even the yield of the reaction. However, a full treatment with these two reaction coordinates is a challenging task and has been left mostly unsolved. In this work, we address this problem by developing a dynamic theory by combining the ideas from ET reaction literature and barrierless chemical reactions. Two-dimensional coupled Smoluchowski equations are employed to compute the time evolution of joint probability distribution for the reactant (P-(1)(X,R,t)) and the product (p((2))(X,R,t)), where X, as is usual in ET reactions, describes the solvent polarization coordinate and R is the distance between the reacting ion pair. The reaction is described by a reaction line (sink) which is a function of X and R obtained by imposing a condition of equal energy on the initial and final states of a reacting ion pair. The resulting two-dimensional coupled equations of motion have been solved numerically using an alternate direction implicit (ADI) scheme (Peaceman and Rachford, J. Soc. Ind. Appl. Math. 1955, 3, 28). The results reveal interesting interplay between polarization relaxation and translational dynamics. The following new results have been obtained. (i) For solvents with slow longitudinal polarization relaxation, the escape probability decreases drastically as the polarization relaxation time increases. We attribute this to caging by polarization of the surrounding solvent, As expected, for the solvents having fast polarization relaxation, the escape probability is independent of the polarization relaxation time. (ii) In the slow relaxation limit, there is a significant dependence of escape probability and average rate on the initial solvent polarization, again displaying the effects of polarization caging. Escape probability increases, and the average rate decreases on increasing the initial polarization. Again, in the fast polarization relaxation limit, there is no effect of initial polarization on the escape probability and the average rate of IR. (iii) For normal and barrierless regions the dependence of escape probability and the rate of IR on initial polarization is stronger than in the inverted region. (iv) Because of the involvement of dynamics along R coordinate, the asymmetrical parabolic (that is, non-Marcus) energy gap dependence of the rate is observed.

A mathematical analysis of nonlinear waves in a fluid filled visco-elastic tube

Our investigations in this paper are centred around the mathematical analysis of a ldquomodal waverdquo problem. We have considered the axisymmetric flow of an inviscid liquid in a thinwalled viscoelastic tube under certain simplifying assumptions. We have first derived the propagation space equations in the long wave limit and also given a general procedure to derive these equations for arbitrary wave length, when the flow is irrotational. We have used the method of operators of multiple scales to derive the nonlinear Schrödinger equation governing the modulation of periodic waves and we have elaborated on the ldquolong modulated wavesrdquo and the ldquomodulated long wavesrdquo. We have also examined the existence and stability of Stokes waves in this system. This is followed by a discussion of the progressive wave solutions of the long wave equations. One of the most important results of our paper is that the propagation space equations are no longer partial differential equations but they are in terms of pseudo-differential operators.Die vorliegenden Untersuchungen beziehen sich auf die mathematische Behandlung des ldquorModalwellenrdquo-Problems. Die achsensymmetrische Strömung einer nichtviskosen Flüssigkeit in einem dünnwandigen viskoelastischen Rohr, unter bestimmten vereinfachenden Annahmen, wird betrachtet. Zuerst werden die Gleichungen des Ausbreitungsraumes im Langwellenbereich abgeleitet und eine allgemeine Methode zur Herleitung dieser Gleichungen für beliebige Wellenlängen bei nichtrotierender Strömung angegeben. Eine Operatorenmethode mit multiplem Maßstab wird verwendet zur Herleitung der nichtlinearen Schrödinger-Gleichung für die Modulation der periodischen Wellen, und die ldquorlangmodulierten Wellenrdquo sowie die ldquormodulierten Langwellenrdquo werden aufgezeigt. Weiters wird die Existenz und die Stabilität der Stokes-Wellen im System untersucht. Anschließend werden die progressiven Wellenlösungen der Langwellengleichungen diskutiert. Eines der wichtigsten Ergebnisse dieser Arbeit ist, daß die Gleichungen des Ausbreitungsraumes keine partiellen Differentialgleichungen mehr sind, sondern Ausdrücke von Pseudo-Differentialoperatoren.

Solvent effect on the spectroscopic properties of 6MAMC and 7MAMC

The absorption and emission spectra of two dyes namely 6MAMC and 7MAMC have been recorded at room temperature in solvents of different polarities. The ground-state dipole moments (mu(g)) of these two were determined experimentally by Guggenheim method and were compared with theoretical values obtained using quantum chemical method. The exited state (mu(e))dipole moments were estimated from Lippert's, Bakhshiev's and Chamma-Viallet's equations by using the variation of the Stokes shift with the solvent dielectric constant and refractive index. The ground and excited-state dipole moments were calculated by means of the solvatochromic shift method and also the excited-state dipole moments are determined in combination with ground-state dipole moments. It was observed that dipole moments of excited state were higher than those of the ground state, indicating a substantial redistribution of the pi-electron densities in a more polar excited state for these two dyes. (C) 2010 Elsevier B.V. All rights reserved.

Solid finite elements through three decades

conventionally, solid finite elements have been looked upon as just generalizations of two-dimensional finite elements. In this article we trace their development starting from the days of their inception. Keeping in tune with our perceptions on developing finite elements, without taking recourse to any extra variational techniques, we discuss a few of the techniques which have been applied to solid finite elements. Finally we critically examine our own work on formulating solid finite elements based on the solutions to the Navier equations.

Single-molecule thermodynamics: the heat distribution function of a charged particle in a static magnetic field

Interest in the applicability of fluctuation theorems to the thermodynamics of single molecules in external potentials has recently led to calculations of the work and total entropy distributions of Brownian oscillators in static and time-dependent electromagnetic fields. These calculations, which are based on solutions to a Smoluchowski equation, are not easily extended to a consideration of the other thermodynamic quantity of interest in such systems-the heat exchanges of the particle alone-because of the nonlinear dependence of the heat on a particle's stochastic trajectory. In this paper, we show that a path integral approach provides an exact expression for the distribution of the heat fluctuations of a charged Brownian oscillator in a static magnetic field. This approach is an extension of a similar path integral approach applied earlier by our group to the calculation of the heat distribution function of a trapped Brownian particle, which was found, in the limit of long times, to be consistent with experimental data on the thermal interactions of single micron-sized colloids in a viscous solvent.

Shear-induced enhancement of self-diffusion in interacting colloidal suspensions

We set up the generalized Langevin equations describing coupled single-particle and collective motion in a suspension of interacting colloidal particles in a shear how and use these to show that the measured self-diffusion coefficients in these systems should be strongly dependent on shear rate epsilon. Three regimes are found: (i) an initial const+epsilon(.2), followed by (ii) a large regime of epsilon(.1/2) behavior, crossing over to an asymptotic power-law approach (iii) D-o - const x epsilon(.-1/2) to the Stokes-Einstein value D-o. The shear dependence is isotropic up to very large shear rates and increases with the interparticle interaction strength. Our results provide a straightforward explanation of recent experiments and simulations on sheared colloids.

Generation of an eight node brick element using Papcovitch-Neuber functions

Some conventional finite elements suffer from drawbacks, such as shear locking, membrane locking, etc. To overcome them researchers have developed various techniques, termed as tricks by some and variational crimes by others. Many attempts have been made, but satisfactory explanations for why some of these techniques work have not been obtained, especially in the case of solid elements. This paper attempts a simple non-conforming solid element using assumed displacement fields which satisfy the Navier equation exactly. Its behaviour under simple loadings like bending, torsion and tension is examined and comparisons are made with existing elements.

PL due to Discrete Gap Levels in Some Chalcogenide Glasses- A Configurational Coordinate Diagram Illustration

Photoluminescence (PL) studies were carried out on a-Se and a few Ge20Se80−xBix and Ge20Se70−xBixTe10 bulk glassy semiconductors at 4.2 K with Ar+ laser as excitation source. While a-Se and samples with lesser at% of Bi show fine structured PL with a large Stokes shift, samples with higher at% of Bi did not show any detectable PL. The investigations show at least three radiative recombination transitions. Features extracted by deconvoluting the experimental spectra show that the discrete gap levels associated with the inherent coordination defects are involved in the PL transitions. Absence of PL in samples with higher Bi at% are explained on the basis of nonradiative transition mechanisms. Overall PL mechanism involving gap levels in chalcogenide glasses is illustrated with the help of a configurational coordinate diagram.

String-like propagation of the 5-coordinated defect state in supercooled water: molecular origin of dynamic and thermodynamic anomalies

We find that at low temperature water, large amplitude (similar to 60 degrees) rotational jumps propagate like a string, with the length of propagation increasing with lowering temperature. The strings are formed by mobile 5-coordinated water molecules which move like a Glarum defect (J. Chem. Phys., 1960, 33, 1371), causing water molecules on the path to change from 4-coordinated to 5-coordinated and again back to 4-coordinated water, and in the process cause the tagged water molecule to jump, by following essentially the Laage-Hynes mechanism (Science, 2006, 311, 832-835). The effects on relaxation of the propagating defect causing large amplitude jumps are manifested most dramatically in the mean square displacement (MSD) and also in the rotational time correlation function of the O-H bond of the molecule that is visited by the defect (transient transition to the 5-coordinated state). The MSD and the decay of rotational time correlation function, both remain quenched in the absence of any visit by the defect, as postulated by Glarum long time ago. We establish a direct connection between these propagating events and the known thermodynamic and dynamic anomalies in supercooled water. These strings are found largely in the regions that surround the relatively rigid domains of 4-coordinated water molecules. The propagating strings give rise to a noticeable dynamical heterogeneity, quantified here by a sharp rise in the peak of the four-point density response function, chi(4)(t). This dynamics heterogeneity is also responsible for the breakdown of the Stokes-Einstein relation.

Scale effects on buckling analysis of orthotropic nanoplates based on nonlocal two-variable refined plate theory

This article presents the buckling analysis of orthotropic nanoplates such as graphene using the two-variable refined plate theory and nonlocal small-scale effects. The two-variable refined plate theory takes account of transverse shear effects and parabolic distribution of the transverse shear strains through the thickness of the plate, hence it is unnecessary to use shear correction factors. Nonlocal governing equations of motion for the monolayer graphene are derived from the principle of virtual displacements. The closed-form solution for buckling load of a simply supported rectangular orthotropic nanoplate subjected to in-plane loading has been obtained by using the Navier's method. Numerical results obtained by the present theory are compared with first-order shear deformation theory for various shear correction factors. It has been proven that the nondimensional buckling load of the orthotropic nanoplate is always smaller than that of the isotropic nanoplate. It is also shown that small-scale effects contribute significantly to the mechanical behavior of orthotropic graphene sheets and cannot be neglected. Further, buckling load decreases with the increase of the nonlocal scale parameter value. The effects of the mode number, compression ratio and aspect ratio on the buckling load of the orthotropic nanoplate are also captured and discussed in detail. The results presented in this work may provide useful guidance for design and development of orthotropic graphene based nanodevices that make use of the buckling properties of orthotropic nanoplates.

Solvent effect on the fluorescence quenching of biologically active carboxamide by aniline and carbon tetrachloride in different solvents using S-V plots

The fluorescence quenching studies of carboxamide namely (E)-N-(3-Chlorophenyl)-2-(3,4,5-trimethoxybenzylideneamino)-4,5,6,7 tetrahydrobenzob]thiophene-3-carboxamide ENCTTTC] by aniline and carbon tetrachloride in six different solvents namely toluene, cyclohexane, n-hexane, n-heptane, n-decane and n-pentane have been carried out at room temperature with a view to understand the quenching mechanisms. The Stern-Volmer (S-V) plots have been found to be nonlinear with a positive deviation for all the solvents studied. In order to interpret these results we have invoked the ground state complex formation and sphere of action static quenching models. Using these models various quenching rate parameters have been determined. The magnitudes of these parameters suggest that sphere of action static quenching model agrees well with the experimental results. Hence the positive deviation is attributed to the static and dynamic quenching. Further, with the use of Finite Sink approximation model, it was possible to check these bimolecular reactions as diffusion-limited and to estimate independently distance parameter R' and mutual diffusion coefficient D. Finally an effort has been made to correlate the values of R' and D with the values of the encounter distance R and the mutual coefficient D determined using the Edward's empirical relation and Stokes Einstein relation. (C) 2011 Elsevier B.V. All rights reserved.