A Smooth Discretization Bridging Finite Element and Mesh-free Methods Using Polynomial Reproducing Simplex Splines

This work sets forth a `hybrid' discretization scheme utilizing bivariate simplex splines as kernels in a polynomial reproducing scheme constructed over a conventional Finite Element Method (FEM)-like domain discretization based on Delaunay triangulation. Careful construction of the simplex spline knotset ensures the success of the polynomial reproduction procedure at all points in the domain of interest, a significant advancement over its precursor, the DMS-FEM. The shape functions in the proposed method inherit the global continuity (Cp-1) and local supports of the simplex splines of degree p. In the proposed scheme, the triangles comprising the domain discretization also serve as background cells for numerical integration which here are near-aligned to the supports of the shape functions (and their intersections), thus considerably ameliorating an oft-cited source of inaccuracy in the numerical integration of mesh-free (MF) schemes. Numerical experiments show the proposed method requires lower order quadrature rules for accurate evaluation of integrals in the Galerkin weak form. Numerical demonstrations of optimal convergence rates for a few test cases are given and the method is also implemented to compute crack-tip fields in a gradient-enhanced elasticity model.

Layered Exact-Repair Regenerating Codes via Embedded Error Correction and Block Designs

A new class of exact-repair regenerating codes is constructed by stitching together shorter erasure correction codes, where the stitching pattern can be viewed as block designs. The proposed codes have the help-by-transfer property where the helper nodes simply transfer part of the stored data directly, without performing any computation. This embedded error correction structure makes the decoding process straightforward, and in some cases the complexity is very low. We show that this construction is able to achieve performance better than space-sharing between the minimum storage regenerating codes and the minimum repair-bandwidth regenerating codes, and it is the first class of codes to achieve this performance. In fact, it is shown that the proposed construction can achieve a nontrivial point on the optimal functional-repair tradeoff, and it is asymptotically optimal at high rate, i.e., it asymptotically approaches the minimum storage and the minimum repair-bandwidth simultaneously.

Pt-Poisoning-Free Efficient CO Oxidation on Pt3Co Supported on MgO(100): An Ab Initio Study

Late-transition-metal-doped Pt clusters are prevalent in CO oxidation catalysis, as they exhibit better catalytic activity than pure Pt, while reducing the effective cost and poisoning However, completely eliminating the critical problem of Pt poisoning still poses a big challenge. Here, we report for the first time that, among the bimetallic clusters ((Pt3M where M = Co, Ni, and Cu)/MgO(100)), the CO adsorption site inverts for Pt3Co/MgO(100) from Pt to Co, due to the complete uptake of Pt d-states by lattice oxygen. While this resolves the problem of Pt poisoning, good reaction kinetics are predicted through low barriers for Langmuir-Hinshelwood and Mars van Krevelen (MvK) mechanisms of CO oxidation for Pt3Co/MgO(100) and Li-doped MgO(100), respectively. Li doping in MgO(100) compensates for the charge imbalance caused by a spontaneous oxygen vacancy formation. Pt-3 Co/Li-doped MgO(100) stands out as an exceptional CO oxidation catalyst, giving an MvK reaction barrier as low as 0.11 eV. We thereby propose a novel design strategy of d-band center inversion for CO oxidation catalysts with no Pt poisoning and excellent reaction kinetics.

A stabilization approach for mesh-free simulations of systems developing shocks or extreme strain localizations

A new stabilization scheme, based on a stochastic representation of the discretized field variables, is proposed with a view to reduce or even eliminate unphysical oscillations in the mesh-free numerical simulations of systems developing shocks or exhibiting localized bands of extreme deformation in the response. The origin of the stabilization scheme may be traced to nonlinear stochastic filtering and, consistent with a class of such filters, gain-based additive correction terms are applied to the simulated solution of the system, herein achieved through the element-free Galerkin method, in order to impose a set of constraints that help arresting the spurious oscillations. The method is numerically illustrated through its Applications to inviscid Burgers' equations, wherein shocks may develop as a result of intersections of the characteristics, and to a gradient plasticity model whose response is often characterized by a developing shear band as the external load is gradually increased. The potential of the method in stabilized yet accurate numerical simulations of such systems involving extreme gradient variations in the response is thus brought forth. (C) 2014 Elsevier Ltd. All rights reserved.

Structural and Magnetic Phase Transformations of Hydroxyapatite-Magnetite Composites under Inert and Ambient Sintering Atmospheres

The present work reports the impact of sintering conditions on the phase stability in hydroxyapatite (HA) magnetite (Fe3O4) bulk composites, which were densified using either pressureless sintering in air or by rapid densification via hot pressing in inert atmosphere. In particular, the phase abundances, structural and magnetic properties of the (1-x)HA-xFe(3)O(4) (x = 5, 10, 20, and 40 wt %) composites were quantified by corroborating results obtained from Rietveld refinement of the X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), and Mossbauer spectroscopy. Post heat treatment phase analysis revealed a major retention of Fe3O4 in argon atmosphere, while it was partially/completely oxidized to hematite (alpha-Fe2O3) in air. Mossbauer results suggest the high-temperature diffusion of Fe3+ into hydroxyapatite lattice, leading to the formation of Fe-doped HA. A preferential occupancy of Fe3+ at the Ca(1) and Ca(2) sites under hot-pressing and conventional sintering conditions, respectively, was observed. The lattice expansion in HA from Rietveld analysis correlated well with the amounts of Fe-doped HA determined from the Mossbauer spectra. Furthermore, hydroxyapatite in the monoliths and composites was delineated to exist in the monoclinic (P2(1)/b) structure as against the widely reported hexagonal (P6(3)/m) crystal lattice. The compositional similarity of iron doping in hydroxyapatite to that of tooth enamel and bone presents HA-Fe3O4 composites as potential orthopedic and dental implant materials.

Phase formation and stability of alloy phases in free nanoparticles: some insights

This paper explores phase formation and phase stability in free nanoparticles of binary alloys. A procedure for estimating the size and composition dependent free energies incorporating the contributions from the interfaces has been presented. Both single phase solid solution and two phase morphology containing interphase interfaces have been considered. A free energy scenario has been evaluated for two binary alloy systems Ag-Ni and Ag-Cu to predict the microstructure of the alloy nanoparticles at different size ranges as a function of composition. Both Ag-Cu and Ag-Ni systems exhibit wide bulk immiscibility. Ag-Ni nanoparticles were synthesized using the wet chemical synthesis technique whereas Ag-Cu nanoparticles were synthesized using laser ablation of a Ag-Cu target immersed in distilled water. Microstructural and compositional characterization of Ag-Ni and Ag-Cu nanoparticles on a single nanoparticle level was conducted using transmission electron microscopy. Nanoparticle microstructures observed from the microscopic investigation have been correlated with thermodynamic calculation results. It is shown that the observed two phase microstructure consisting of Ag-Ni solid solution in partial decomposed state coexisting with pure Ag phases in the case of Ag-Ni nanoparticles can be only be rationalized by invoking the tendency for phase separation of an initial solid solution with increase in nanoparticle size. Smaller sized Ag-Ni nanoparticles prefer a single phase solid solution microstructure. Due to an increase in particle size during the synthesis process the initial solid solution decomposes into an ultrafine scale phase separated microstructure. We have shown that it is necessary to invoke critical point phenomenon and wetting transition in systems showing a critical point that leads to phase separated Ag-Ni nanoparticles providing a catalytic substrate for the nucleation of equilibrium Ag over it. In the case of the Ag-Cu system, we report the experimental observation of a core shell structure at small sizes. This can be rationalized in terms of a metastable solid solution. It is argued that the nucleation barrier can prevent the formation of biphasic morphology with an internal interface. In such a situation, demixing of the solid solution can bring the system to a lower energy configuration. This has lead to the observed core-shell morphology in the Ag-Cu system during room temperature synthesis.

Heterochromatic paths in edge colored graphs without small cycles and heterochromatic-triangle-free graphs

Conditions for the existence of heterochromatic Hamiltonian paths and cycles in edge colored graphs are well investigated in literature. A related problem in this domain is to obtain good lower bounds for the length of a maximum heterochromatic path in an edge colored graph G. This problem is also well explored by now and the lower bounds are often specified as functions of the minimum color degree of G - the minimum number of distinct colors occurring at edges incident to any vertex of G - denoted by v(G). Initially, it was conjectured that the lower bound for the length of a maximum heterochromatic path for an edge colored graph G would be 2v(G)/3]. Chen and Li (2005) showed that the length of a maximum heterochromatic path in an edge colored graph G is at least v(G) - 1, if 1 <= v(G) <= 7, and at least 3v(G)/5] + 1 if v(G) >= 8. They conjectured that the tight lower bound would be v(G) - 1 and demonstrated some examples which achieve this bound. An unpublished manuscript from the same authors (Chen, Li) reported to show that if v(G) >= 8, then G contains a heterochromatic path of length at least 120 + 1. In this paper, we give lower bounds for the length of a maximum heterochromatic path in edge colored graphs without small cycles. We show that if G has no four cycles, then it contains a heterochromatic path of length at least v(G) - o(v(G)) and if the girth of G is at least 4 log(2)(v(G)) + 2, then it contains a heterochromatic path of length at least v(G) - 2, which is only one less than the bound conjectured by Chen and Li (2005). Other special cases considered include lower bounds for the length of a maximum heterochromatic path in edge colored bipartite graphs and triangle-free graphs: for triangle-free graphs we obtain a lower bound of 5v(G)/6] and for bipartite graphs we obtain a lower bound of 6v(G)-3/7]. In this paper, it is also shown that if the coloring is such that G has no heterochromatic triangles, then G contains a heterochromatic path of length at least 13v(G)/17)]. This improves the previously known 3v(G)/4] bound obtained by Chen and Li (2011). We also give a relatively shorter and simpler proof showing that any edge colored graph G contains a heterochromatic path of length at least (C) 2015 Elsevier Ltd. All rights reserved.

Chemical-free graphene by unzipping carbon nanotubes using cryo-milling

Current paper reports synthesis of chemical free graphene by unzipping of the carbon nanotubes (CNTs) using high strain rate deformation at 150K. A specially designed cryomill operating at 150 K was used for the experiments. The mechanism of unzipping was further explored using molecular dynamics (MD) simulations. Both experimental and simulation results reveal two modes of unzipping through radial and shear loading. (C) 2015 Elsevier Ltd. All rights reserved.

Structural crossover from nonmodulated to long-period modulated tetragonal phase and anomalous change in ferroelectric properties in the lead-free piezoelectric Na1/2Bi1/2TiO3-BaTiO3

The highly complex structure-property interrelationship in the lead-free piezoelectric (x) Na1/2Bi1/2TiO3 - (1 - x) BaTiO3 is a subject of considerable contemporary debate. Using comprehensive x-ray, neutron diffraction, dielectric, and ferroelectric studies, we have shown the existence of a new criticality in this system at x = 0.80, i.e., well within the conventional tetragonal phase field. This criticality manifests as a nonmonotonic variation of the tetragonality and coercivity and is shown to be associated with a crossover from a nonmodulated tetragonal phase (for x < 0.8) to a long-period modulated tetragonal phase (for x > 0.80). It is shown that the stabilization of long-period modulation introduces a characteristic depolarization temperature in the system. While differing qualitatively from the two-phase model often suggested for the critical compositions of this system, our results support the view with regard to the tendency in perovskites to stabilize long-period modulated structures as a result of complex interplay of antiferrodistortive modes Bellaiche and Iniguez, Phys. Rev. B 88, 014104 ( 2013); Prosandeev, Wang, Ren, Iniguez, ands Bellaiche, Adv. Funct. Mater. 23, 234 (2013)].

A Reliable Residual Based A Posteriori Error Estimator for a Quadratic Finite Element Method for the Elliptic Obstacle Problem

A residual based a posteriori error estimator is derived for a quadratic finite element method (FEM) for the elliptic obstacle problem. The error estimator involves various residuals consisting of the data of the problem, discrete solution and a Lagrange multiplier related to the obstacle constraint. The choice of the discrete Lagrange multiplier yields an error estimator that is comparable with the error estimator in the case of linear FEM. Further, an a priori error estimate is derived to show that the discrete Lagrange multiplier converges at the same rate as that of the discrete solution of the obstacle problem. The numerical experiments of adaptive FEM show optimal order convergence. This demonstrates that the quadratic FEM for obstacle problem exhibits optimal performance.

A Remark on the A Posteriori Error Analysis of Discontinuous Galerkin Methods for the Obstacle Problem

We revisit the a posteriori error analysis of discontinuous Galerkin methods for the obstacle problem derived in 25]. Under a mild assumption on the trace of obstacle, we derive a reliable a posteriori error estimator which does not involve min/max functions. A key in this approach is an auxiliary problem with discrete obstacle. Applications to various discontinuous Galerkin finite element methods are presented. Numerical experiments show that the new estimator obtained in this article performs better.

Effect of a mesh on boundary layer transitions induced by free-stream turbulence and an isolated roughness element

Streamwise streaks, their lift-up and streak instability are integral to the bypass transition process. An experimental study has been carried out to find the effect of a mesh placed normal to the flow and at different wall-normal locations in the late stages of two transitional flows induced by free-stream turbulence (FST) and an isolated roughness element. The mesh causes an approximately 30% reduction in the free-stream velocity, and mild acceleration, irrespective of its wall-normal location. Interestingly, when located near the wall, the mesh suppresses several transitional events leading to transition delay over a large downstream distance. The transition delay is found to be mainly caused by suppression of the lift-up of the high-shear layer and its distortion, along with modification of the spanwise streaky structure to an orderly one. However, with the mesh well away from the wall, the lifted-up shear layer remains largely unaffected, and the downstream boundary layer velocity profile develops an overshoot which is found to follow a plane mixing layer type profile up to the free stream. Reynolds stresses, and the size and strength of vortices increase in this mixing layer region. This high-intensity disturbance can possibly enhance transition of the accelerated flow far downstream, although a reduction in streamwise turbulence intensity occurs over a short distance downstream of the mesh. However, the shape of the large-scale streamwise structure in the wall-normal plane is found to be more or less the same as that without the mesh.

Electric field induced short range to long range structural ordering and its influence on the Eu+3 photoluminescence in the lead-free ferroelectric Na1/2Bi1/2TiO3

Eu+3 was incorporated into the lattice of a lead-free ferroelectric Na1/2Bi1/2TiO3 (NBT) as per the nominal formula Na0.5Bi0.5-xEuxTiO3. This system was investigated with regard to the Eu+3 photoluminescence (PL) and structural behaviour as a function of composition and electric field. Electric field was found to irreversibly change the features in the PL spectra and also in the x-ray diffraction patterns below the critical composition x = 0.025. Detailed analysis revealed that below the critical composition, electric field irreversibly suppresses the structural heterogeneity inherent of the host matrix NBT and brings about a long range ferroelectric state with rhombohedral (R3c) distortion. It is shown that the structural disorder on the nano-scale opens a new channel for radiative transition which manifests as a new emission line branching off from the main D-5(0)-> F-7(0) line along with a concomitant change in the relative intensity of the other crystal field induced Stark lines with different J values. The study suggests that Eu+3 luminescence can be used to probe the relative degree of field induced structural ordering in relaxor ferroelectrics and also in high performance piezoelectric alloys where electric field couples very strongly with the lattice and structural degrees of freedom. (C) 2015 AIP Publishing LLC.

Composition dependent multiple structural transformations of myoglobin in aqueous ethanol solution: A combined experimental and theoretical study

Experimental studies (circular dichroism and ultra-violet (UV) absorption spectra) and large scale atomistic molecular dynamics simulations (accompanied by order parameter analyses) are combined to establish a number of remarkable (and unforeseen) structural transformations of protein myoglobin in aqueous ethanol mixture at various ethanol concentrations. The following results are particularly striking. (1) Two well-defined structural regimes, one at x(EtOH) similar to 0.05 and the other at x(EtOH) similar to 0.25, characterized by formation of distinct partially folded conformations and separated by a unique partially unfolded intermediate state at x(EtOH) similar to 0.15, are identified. (2) Existence of non-monotonic composition dependence of (i) radius of gyration, (ii) long range contact order, (iii) residue specific solvent accessible surface area of tryptophan, and (iv) circular dichroism spectra and UV-absorption peaks are observed. Interestingly at x(EtOH) similar to 0.15, time averaged value of the contact order parameter of the protein reaches a minimum, implying that this conformational state can be identified as a molten globule state. Multiple structural transformations well known in water-ethanol binary mixture appear to have considerably stronger effects on conformation and dynamics of the protein. We compare the present results with studies in water-dimethyl sulfoxide mixture where also distinct structural transformations are observed along with variation of co-solvent composition. (C) 2015 AIP Publishing LLC.

Contrasting Effects of Graphene Oxide and Poly(ethylenimine) on the Polymorphism in Poly(vinylidene fluoride)

The nature of interaction between a heteronucleating agent (graphene oxide, GO) and a strongly polar macromolecule (poly(ethylenimine), PEI) with poly(vinylidene fluoride) (PVDF) influencing the crystalline structure and morphology has been systematically investigated in this work. PEI interacts with PVDF via ion-dipole interaction, which helps in lowering the free energy barrier for nucleation thereby promoting faster crystallization. In contrast, besides interacting with PVDF, GO also promotes heteronucleation in PVDF. We observed that both GO and PEI have very different effects on the overall crystalline morphology of PVDF. For instance, the neat PVDF showed a mixture of both alpha and beta phases when cooled from the melt. However, incorporation of 0.1 wt % GO resulted in phase transformation from the stable alpha-phase to polar beta-polymorph in PVDF. In contrast, PEI, which also resulted in faster crystallization in PVDF predominantly, resulted in the stable alpha- phase. Various techniques like Fourier transform infrared spectroscopy, X-ray diffraction, and differential scanning calorimetry were employed to confirm the phase transformations in PVDF. PEI was further grafted onto GO nanosheets to understand the combined effects of both GO and PEI on the polymorphism in PVDF. The PVDF/PEI-GO composite showed a mixture of phases, predominantly rich in a. These phenomenal effects were further analyzed and corroborated with the specific interaction between GO and PEI with PVDF using X-ray photon scattering (XPS) and NMR. In addition, the dielectric permittivity increased significantly in the presence of GO and PEI in the composites. For instance, PVDF/PEI-GO showed the highest permittivity of 39 at 100 Hz.