Mixed finite elements for electromagnetic analysis

The occurrence of spurious solutions is a well-known limitation of the standard nodal finite element method when applied to electromagnetic problems. The two commonly used remedies that are used to address this problem are (i) The addition of a penalty term with the penalty factor based on the local dielectric constant, and which reduces to a Helmholtz form on homogeneous domains (regularized formulation); (ii) A formulation based on a vector and a scalar potential. Both these strategies have some shortcomings. The penalty method does not completely get rid of the spurious modes, and both methods are incapable of predicting singular eigenvalues in non-convex domains. Some non-zero spurious eigenvalues are also predicted by these methods on non-convex domains. In this work, we develop mixed finite element formulations which predict the eigenfrequencies (including their multiplicities) accurately, even for nonconvex domains. The main feature of the proposed mixed finite element formulation is that no ad-hoc terms are added to the formulation as in the penalty formulation, and the improvement is achieved purely by an appropriate choice of finite element spaces for the different variables. We show that the formulation works even for inhomogeneous domains where `double noding' is used to enforce the appropriate continuity requirements at an interface. For two-dimensional problems, the shape of the domain can be arbitrary, while for the three-dimensional ones, with our current formulation, only regular domains (which can be nonconvex) can be modeled. Since eigenfrequencies are modeled accurately, these elements also yield accurate results for driven problems. (C) 2014 Elsevier Ltd. All rights reserved.

An unusual demixing behavior in PS-PVME blends in the presence of nanoparticles

The effect of silver nanoparticles (sNP) on the demixing and the evolution of morphology in off-critical blends of 90/10 and 10/90 (wt/wt) PS/PVME polystyrene/poly(vinyl methyl ether)] was probed here using shear rheology and optical microscopy. The faster component (PVME) has a higher molecular weight (80 kDa) than the slower component (PS, 35 kDa), which makes this system quite interesting to study with respect to the evolving morphology, as the blends transit through the binodal and the spinodal envelopes. An unusual demixing behavior was observed in both PVME rich and PS rich blends. Temperature modulated differential scanning calorimetry measurements showed that the T-g value for the blends with sNP was slightly lower than that of the neat blends. A decreased volume of cooperativity at T-g suggests confined segmental dynamics in the presence of sNP. Although, the addition of sNP had no influence on the thermodynamic demixing temperature, it significantly altered the elasticity of the minor component during the transition of the blend from the homogeneous to the heterogeneous state. This is manifested from energetically driven localization of the sNP in the PVME phase during demixing. As a direct consequence of this, the formation of the microstructures upon demixing was observed to be delayed in the presence of sNP. Interestingly, in the intermediate quench depth, the higher viscoelastic phase evolved as an interconnected network, which subsequently coarsened into discrete droplets in the late stages for the 90/10 PS/PVME blends. Similar observations were made for 10/90 PS/PVME blends where threads of PVME appeared at deeper quench depths in the presence of sNP. The interconnected network formation of the minor phase (here PVME), which is also the faster component in the blend, was different from the usual demixing behavior.

Optical and structural properties of highly porous shell structured Fe doped TiO2 thin films

TiO2 thin films with 0.2 wt%, 0.4 wt%, 0.6 wt%, and 0.8 wt% Fe were prepared on glass and silicon substrates using sol-gel spin coating technique. The optical cut-off points are increasingly red-shifted and the absorption edge is shifted over the higher wavelength region with Fe content increasing. As Fe content increases, the optical band gap decreases from 3.03 to 2.48 eV whereas the tail width increases from 0.26 to 1.43 eV. The X-ray diffraction (XRD) patterns for doped films at 0.2 wt% and 0.8 wt% Fe reveal no characteristic peaks, indicating that the film is amorphous whereas undoped TiO2 exhibits (101) orientation with anatase phase. Thin films of higher Fe content exhibit a homogeneous, uniform, and nano-structured highly porous shell morphology.

Influence of Spatial Variability on Pavement Responses Using Latin Hypercube Sampling on Two-Dimensional Random Fields

Although uncertainties in material properties have been addressed in the design of flexible pavements, most current modeling techniques assume that pavement layers are homogeneous. The paper addresses the influence of the spatial variability of the resilient moduli of pavement layers by evaluating the effect of the variance and correlation length on the pavement responses to loading. The integration of the spatially varying log-normal random field with the finite-difference method has been achieved through an exponential autocorrelation function. The variation in the correlation length was found to have a marginal effect on the mean values of the critical strains and a noticeable effect on the standard deviation which decreases with decreases in correlation length. This reduction in the variance arises because of the spatial averaging phenomenon over the softer and stiffer zones generated because of spatial variability. The increase in the mean value of critical strains with decreasing correlation length, although minor, illustrates that pavement performance is adversely affected by the presence of spatially varying layers. The study also confirmed that the higher the variability in the pavement layer moduli, introduced through a higher value of coefficient of variation (COV), the higher the variability in the pavement response. The study concludes that ignoring spatial variability by modeling the pavement layers as homogeneous that have very short correlation lengths can result in the underestimation of the critical strains and thus an inaccurate assessment of the pavement performance. (C) 2014 American Society of Civil Engineers.

CLAP: A web-server for automatic classification of proteins with special reference to multi-domain proteins

Background: The function of a protein can be deciphered with higher accuracy from its structure than from its amino acid sequence. Due to the huge gap in the available protein sequence and structural space, tools that can generate functionally homogeneous clusters using only the sequence information, hold great importance. For this, traditional alignment-based tools work well in most cases and clustering is performed on the basis of sequence similarity. But, in the case of multi-domain proteins, the alignment quality might be poor due to varied lengths of the proteins, domain shuffling or circular permutations. Multi-domain proteins are ubiquitous in nature, hence alignment-free tools, which overcome the shortcomings of alignment-based protein comparison methods, are required. Further, existing tools classify proteins using only domain-level information and hence miss out on the information encoded in the tethered regions or accessory domains. Our method, on the other hand, takes into account the full-length sequence of a protein, consolidating the complete sequence information to understand a given protein better. Results: Our web-server, CLAP (Classification of Proteins), is one such alignment-free software for automatic classification of protein sequences. It utilizes a pattern-matching algorithm that assigns local matching scores (LMS) to residues that are a part of the matched patterns between two sequences being compared. CLAP works on full-length sequences and does not require prior domain definitions. Pilot studies undertaken previously on protein kinases and immunoglobulins have shown that CLAP yields clusters, which have high functional and domain architectural similarity. Moreover, parsing at a statistically determined cut-off resulted in clusters that corroborated with the sub-family level classification of that particular domain family. Conclusions: CLAP is a useful protein-clustering tool, independent of domain assignment, domain order, sequence length and domain diversity. Our method can be used for any set of protein sequences, yielding functionally relevant clusters with high domain architectural homogeneity. The CLAP web server is freely available for academic use at http://nslab.mbu.iisc.ernet.in/clap/.

A meshless adaptive multiscale method for fracture

The paper presents a multiscale method for crack propagation. The coarse region is modelled by the differential reproducing kernel particle method. Fracture in the coarse scale region is modelled with the Phantom node method. A molecular statics approach is employed in the fine scale where crack propagation is modelled naturally by breaking of bonds. The triangular lattice corresponds to the lattice structure of the (111) plane of an FCC crystal in the fine scale region. The Lennard-Jones potential is used to model the atom-atom interactions. The coupling between the coarse scale and fine scale is realized through ghost atoms. The ghost atom positions are interpolated from the coarse scale solution and enforced as boundary conditions on the fine scale. The fine scale region is adaptively refined and coarsened as the crack propagates. The centro symmetry parameter is used to detect the crack tip location. The method is implemented in two dimensions. The results are compared to pure atomistic simulations and show excellent agreement. (C) 2014 Elsevier B. V. All rights reserved.

Cavitation in brittle metallic glasses - Effects of stress state and distributed weak zones

Experimental studies and atomistic simulations have shown that brittle metallic glasses fail by a cavitation mechanism whose origin has been traced to the presence of intrinsic atomic density fluctuations which give rise to weak zones with reduced yield strength. It has been shown recently through continuum analysis that the presence of these zones can lower the cavitation stress considerably under equibiaxial loading. The objective of the present work is to study the effect of the applied stress state on the cavitation behavior of such a heterogeneous plastic solid with distributed weak zones. To this end, 2D plane strain finite element simulations are performed by subjecting a unit cell containing a weak zone to different (biaxiality) stress ratios. The volume fraction and yield strength of the weak zone are varied over a wide range. The results show that unlike in a homogeneous plastic solid, the cavitation stress of the heterogeneous aggregate does not reduce appreciably as the stress ratio decreases from unity when the yield strength of the weak zone is low. It is found that a non-dimensional parameter characterizing the stress state prevailing in the weak zone and its yield properties uniquely control the cavitation stress. The nature of cavitation bifurcation may change from unstable bifurcation to the left at sufficiently low stress ratio to one involving snap cavitation at high stress ratio. (C) 2014 Elsevier Ltd. All rights reserved.

Relaxor behaviour in BaBi4Ti4O15 ceramics fabricated using the powders obtained by mechanochemically assisted synthesis route

Mechanochemically activated reactants were found to facilitate the synthesis of fine powders comprising 200-400 nm range crystallites of BaBi4Ti4O15 at a significantly lower temperature (700 A degrees C) than that of solid-state reaction route. Reactants (CaCO3, Bi2O3 and TiO2) in stoichiometric ratio were ball milled for 48 h to obtain homogeneous mixture. The evolution of the BaBi4Ti4O15 phase was systematically followed using X-ray powder diffraction (XRD) technique. Scanning electron microscopy (SEM) and transmission electron microscopy (TEM) were employed to probe its structural and microstructural details. The electron diffraction studies established the presence of correlated octahedral rotations and associated long-range polar ordering. High-resolution TEM imaging nevertheless revealed structural inhomogeneities leading to intergrowth defects. Dense BaBi4Ti4O15 ceramics with an average grain size of 0.9 mu m were fabricated using mechanochemically assisted synthesized powders at relatively low temperature (1000 A degrees C). The effect of grain size on the dielectric and relaxor behaviour of BaBi4Ti4O15 ceramics was investigated. Fine-grained ceramics (average grain size similar to 0.9 mu m) showed higher diffusion in phase transition, lower temperature of phase transition, lower Vogel-Fulcher freezing temperature and higher activation energy for the polarization reversal than those for coarse-grained ceramics (average grain size similar to 7 mu m) fabricated via the conventional solid-state reaction route.

Texture discriminability in monkey inferotemporal cortex predicts human texture perception

Shape and texture are both important properties of visual objects, but texture is relatively less understood. Here, we characterized neuronal responses to discrete textures in monkey inferotemporal (IT) cortex and asked whether they can explain classic findings in human texture perception. We focused on three classic findings on texture discrimination: 1) it can be easy or hard depending on the constituent elements; 2) it can have asymmetries, and 3) it is reduced for textures with randomly oriented elements. We recorded neuronal activity from monkey inferotemporal (IT) cortex and measured texture perception in humans for a variety of textures. Our main findings are as follows: 1) IT neurons show congruent selectivity for textures across array size; 2) textures that were easy for humans to discriminate also elicited distinct patterns of neuronal activity in monkey IT; 3) texture pairs with asymmetries in humans also exhibited asymmetric variation in firing rate across monkey IT; and 4) neuronal responses to randomly oriented textures were explained by an average of responses to homogeneous textures, which rendered them less discriminable. The reduction in discriminability of monkey IT neurons predicted the reduced discriminability in humans during texture discrimination. Taken together, our results suggest that texture perception in humans is likely based on neuronal representations similar to those in monkey IT.

Matrix Crack Detection in Composite Plate with Spatially Random Material Properties using Fractal Dimension

Fractal dimension based damage detection method is investigated for a composite plate with random material properties. Composite material shows spatially varying random material properties because of complex manufacturing processes. Matrix cracks are considered as damage in the composite plate. Such cracks are often seen as the initial damage mechanism in composites under fatigue loading and also occur due to low velocity impact. Static deflection of the cantilevered composite plate with uniform loading is calculated using the finite element method. Damage detection is carried out based on sliding window fractal dimension operator using the static deflection. Two dimensional homogeneous Gaussian random field is generated using Karhunen-Loeve (KL) expansion to represent the spatial variation of composite material property. The robustness of fractal dimension based damage detection method is demonstrated considering the composite material properties as a two dimensional random field.

Properties of the chalcogenide-carbon nano tubes and graphene composite materials

Composite can deliver more than the individual elemental property of the material. Specifically chalcogenide- multi walled carbon nano tubes and chalcogenide- bilayer graphene composite materials could be interesting for the investigation, which have been less covered by the investigators. We describe micro structural properties of Se55Te25Ge20, Se55Te25Ge20 + 0.025% multi walled carbon nano tubes and Se55Te25Ge20 + 0.025% bilayer graphene materials. This gives realization of the alloying constituents inclusion/or diffusion inside the multi walled carbon nano tubes and bilayer graphene under the homogeneous parent alloy configuration. Raman spectroscopy, X-ray photoelectron spectroscopy, UV/Visible spectroscopy and Fourier transmission infrared spectroscopy have also been carried out under the discussion. A considerable core energy levels peak shifts have been noticed for the composite materials by the X-ray photoelectron spectroscopy. The optical energy band gaps are measured to be varied in between 1.2 and 1.3 eV. In comparison to parent (Se55Te25Ge20) alloy a higher infrared transmission has been observed for the composite materials. Subsequently, variation in physical properties has been explained on the basis of bond formation in solids. (C) 2014 Elsevier B. V. All rights reserved.

Study of Microstructure Evolution during Semi-Solid Processing of an In-Situ Al Alloy Composite

In the present work, effect of pouring temperature (650 degrees C, 655 degrees C, and 660 degrees C) on semi-solid microstructure evolution of in-situ magnesium silicide (Mg2Si) reinforced aluminum (Al) alloy composite has been studied. The shear force exerted by the cooling slope during gravity driven flow of the melt facilitates the formation of near spherical primary Mg2Si and primary Al grains. Shear driven melt flow along the cooling slope and grain fragmentation have been identified as the responsible mechanisms for refinement of primary Mg2Si and Al grains with improved sphericity. Results show that, while flowing down the cooling slope, morphology of primary Mg2Si and primary Al transformed gradually from coarse dendritic to mixture of near spherical particles, rosettes, and degenerated dendrites. In terms of minimum grain size and maximum sphericity, 650 degrees C has been identified as the ideal pouring temperature for the cooling slope semi-solid processing of present Al alloy composite. Formation of spheroidal grains with homogeneous distribution of reinforcing phase (Mg2Si) improves the isotropic property of the said composite, which is desirable in most of the engineering applications.

Matrix crack detection in spatially random composite structures using fractal dimension

Fractal dimension based damage detection method is studied for a composite structure with random material properties. A composite plate with localized matrix crack is considered. Matrix cracks are often seen as the initial damage mechanism in composites. Fractal dimension based method is applied to the static deformation curve of the structure to detect localized damage. Static deflection of a cantilevered composite plate under uniform loading is calculated using the finite element method. Composite material shows spatially varying random material properties because of complex manufacturing processes. Spatial variation of material property is represented as a two dimensional homogeneous Gaussian random field. Karhunen-Loeve (KL) expansion is used to generate a random field. The robustness of fractal dimension based damage detection methods is studied considering the composite plate with spatial variation in material properties.

Faster computation of the Karhunen-Loeve expansion using its domain independence property

The goal of this work is to reduce the cost of computing the coefficients in the Karhunen-Loeve (KL) expansion. The KL expansion serves as a useful and efficient tool for discretizing second-order stochastic processes with known covariance function. Its applications in engineering mechanics include discretizing random field models for elastic moduli, fluid properties, and structural response. The main computational cost of finding the coefficients of this expansion arises from numerically solving an integral eigenvalue problem with the covariance function as the integration kernel. Mathematically this is a homogeneous Fredholm equation of second type. One widely used method for solving this integral eigenvalue problem is to use finite element (FE) bases for discretizing the eigenfunctions, followed by a Galerkin projection. This method is computationally expensive. In the current work it is first shown that the shape of the physical domain in a random field does not affect the realizations of the field estimated using KL expansion, although the individual KL terms are affected. Based on this domain independence property, a numerical integration based scheme accompanied by a modification of the domain, is proposed. In addition to presenting mathematical arguments to establish the domain independence, numerical studies are also conducted to demonstrate and test the proposed method. Numerically it is demonstrated that compared to the Galerkin method the computational speed gain in the proposed method is of three to four orders of magnitude for a two dimensional example, and of one to two orders of magnitude for a three dimensional example, while retaining the same level of accuracy. It is also shown that for separable covariance kernels a further cost reduction of three to four orders of magnitude can be achieved. Both normal and lognormal fields are considered in the numerical studies. (c) 2014 Elsevier B.V. All rights reserved.