An extended finite element model coupled with level set method for stable pitting corrosion

Mass balance between metal and electrolytic solution, separated by a moving interface, in stable pit growth results in a set of governing equations which are solved for concentration field and interface position (pit boundary evolution), which requires only three inputs, namely the solid metal concentration, saturation concentration of the dissolved metal ions and diffusion coefficient. A combined eXtended Finite Element Model (XFEM) and level set method is developed in this paper. The extended finite element model handles the jump discontinuity in the metal concentrations at the interface, by using discontinuous-derivative enrichment formulation for concentration discontinuity at the interface. This eliminates the requirement of using front conforming mesh and re-meshing after each time step as in conventional finite element method. A numerical technique known as level set method tracks the position of the moving interface and updates it over time. Numerical analysis for pitting corrosion of stainless steel 304 is presented. The above proposed method is validated by comparing the numerical results with experimental results, exact solutions and some other approximate solutions.

Frequency Domain Based Solution for Certain Class of Wave Equations: An exhaustive study of Numerical Solutions

The paper discusses the frequency domain based solution for a certain class of wave equations such as: a second order partial differential equation in one variable with constant and varying coefficients (Cantilever beam) and a coupled second order partial differential equation in two variables with constant and varying coefficients (Timoshenko beam). The exact solution of the Cantilever beam with uniform and varying cross-section and the Timoshenko beam with uniform cross-section is available. However, the exact solution for Timoshenko beam with varying cross-section is not available. Laplace spectral methods are used to solve these problems exactly in frequency domain. The numerical solution in frequency domain is done by discretisation in space by approximating the unknown function using spectral functions like Chebyshev polynomials, Legendre polynomials and also Normal polynomials. Different numerical methods such as Galerkin Method, Petrov- Galerkin method, Method of moments and Collocation method or the Pseudo-spectral method in frequency domain are studied and compared with the available exact solution. An approximate solution is also obtained for the Timoshenko beam with varying cross-section using Laplace Spectral Element Method (LSEM). The group speeds are computed exactly for the Cantilever beam and Timoshenko beam with uniform cross-section and is compared with the group speeds obtained numerically. The shear mode and the bending modes of the Timoshenko beam with uniform cross-section are separated numerically by applying a modulated pulse as the shear force and the corresponding group speeds for varying taper parameter in are obtained numerically by varying the frequency of the input pulse. An approximate expression for calculating group speeds corresponding to the shear mode and the bending mode, and also the cut-off frequency is obtained. Finally, we show that the cut-off frequency disappears for large in, for epsilon > 0 and increases for large in, for epsilon < 0.

Microscopic description of the evolution of the local structure and an evaluation of the chemical pressure concept in a solid solution

Extended x-ray absorption fine-structure studies have been performed at the Zn K and Cd K edges for a series of solid solutions of wurtzite Zn1-xCdxS samples with x = 0.0, 0.1, 0.25, 0.5, 0.75, and 1.0, where the lattice parameter as a function of x evolves according to the well-known Vegard's law. In conjunction with extensive, large-scale first-principles electronic structure calculations with full geometry optimizations, these results establish that the percentage variation in the nearest-neighbor bond distances are lower by nearly an order of magnitude compared to what would be expected on the basis of lattice parameter variation, seriously undermining the chemical pressure concept. With experimental results that allow us to probe up to the third coordination shell distances, we provide a direct description of how the local structure, apparently inconsistent with the global structure, evolves very rapidly with interatomic distances to become consistent with it. We show that the basic features of this structural evolution with the composition can be visualized with nearly invariant Zn-S-4 and Cd-S-4 tetrahedral units retaining their structural integrity, while the tilts between these tetrahedral building blocks change with composition to conform to the changing lattice parameters according to the Vegard's law within a relatively short length scale. These results underline the limits of applicability of the chemical pressure concept that has been a favored tool of experimentalists to control physical properties of a large variety of condensed matter systems.

An extended finite-element model coupled with level set method for analysis of growth of corrosion pits in metallic structures

Mass balance between metal and electrolytic solution, separated by a moving interface, in stable pit growth results in a set of governing equations which are solved for concentration field and interface position (pit boundary evolution). The interface experiences a jump discontinuity in metal concentration. The extended finite-element model (XFEM) handles this jump discontinuity by using discontinuous-derivative enrichment formulation, eliminating the requirement of using front conforming mesh and re-meshing after each time step as in the conventional finite-element method. However, prior interface location is required so as to solve the governing equations for concentration field for which a numerical technique, the level set method, is used for tracking the interface explicitly and updating it over time. The level set method is chosen as it is independent of shape and location of the interface. Thus, a combined XFEM and level set method is developed in this paper. Numerical analysis for pitting corrosion of stainless steel 304 is presented. The above proposed model is validated by comparing the numerical results with experimental results, exact solutions and some other approximate solutions. An empirical model for pitting potential is also derived based on the finite-element results. Studies show that pitting profile depends on factors such as ion concentration, solution pH and temperature to a large extent. Studying the individual and combined effects of these factors on pitting potential is worth knowing, as pitting potential directly influences corrosion rate.

Current kinematics and dynamics of Africa and the East African Rift System

Although the East African Rift System (EARS) is an archetype continental rift, the forces driving its evolution remain debated. Some contend buoyancy forces arising from gravitational potential energy (GPE) gradients within the lithosphere drive rifting. Others argue for a major role of the diverging mantle flow associated with the African Superplume. Here we quantify the forces driving present-day continental rifting in East Africa by (1) solving the depth averaged 3-D force balance equations for 3-D deviatoric stress associated with GPE, (2) inverting for a stress field boundary condition that we interpret as originating from large-scale mantle tractions, (3) calculating dynamic velocities due to lithospheric buoyancy forces, lateral viscosity variations, and velocity boundary conditions, and (4) calculating dynamic velocities that result from the stress response of horizontal mantle tractions acting on a viscous lithosphere in Africa and surroundings. We find deviatoric stress associated with lithospheric GPE gradients are similar to 8-20 MPa in EARS, and the minimum deviatoric stress resulting from basal shear is similar to 1.6 MPa along the EARS. Our dynamic velocity calculations confirm that a force contribution from GPE gradients alone is sufficient to drive Nubia-Somalia divergence and that additional forcing from horizontal mantle tractions overestimates surface kinematics. Stresses from GPE gradients appear sufficient to sustain present-day rifting in East Africa; however, they are lower than the vertically integrated strength of the lithosphere along most of the EARS. This indicates additional processes are required to initiate rupture of continental lithosphere, but once it is initiated, lithospheric buoyancy forces are enough to maintain rifting.

Phase Evolution in the AuCu/Sn System by Solid-State Reactive Diffusion

The interfacial reactions between several Au(Cu) alloys and pure Sn were studied experimentally at 200A degrees C. Amounts of Cu in the AuSn4 and AuSn2 phases were as low as 1 at.%. On the basis of these experimental results there is no continuous solid solution between (Au,Cu)Sn and (Cu,Au)(6)Sn-5. The copper content of (Au,Cu)Sn was determined to be approximately 7-8 at.%. Substantial amounts of Au were present in the (Cu,Au)(6)Sn-5 and (Cu,Au)(3)Sn phases. Two ternary compounds were formed, one with stoichiometry varying from (Au40.5Cu39)Sn-20.5 to (Au20.2Cu59.3)Sn-20.5 (ternary ``B''), the other with the composition Au34Cu33Sn33 (ternary ``C''). The measured phase boundary compositions of the product phases are plotted on the available Au-Cu-Sn isotherm and the phase equilibria are discussed. The complexity and average thickness of the diffusion zone decreases with increasing Cu content except for the Au(40 at.%Cu) couple.

Effect of Rayleigh numbers on the evolution of double-diffusive salt fingers

This is a transient two-dimensional numerical study of double-diffusive salt fingers in a two-layer heat-salt system for a wide range of initial density stability ratio (R-rho 0) and thermal Rayleigh numbers (Ra-T similar to 10(3) - 10(11)). Salt fingers have been studied for several decades now, but several perplexing features of this rich and complex system remain unexplained. The work in question studies this problem and shows the morphological variation in fingers from low to high thermal Rayleigh numbers, which have been missed by the previous investigators. Considerable variations in convective structures and evolution pattern were observed in the range of Ra-T used in the simulation. Evolution of salt fingers was studied by monitoring the finger structures, kinetic energy, vertical profiles, velocity fields, and transient variation of R-rho(t). The results show that large scale convection that limits the finger length was observed only at high Rayleigh numbers. The transition from nonlinear to linear convection occurs at about Ra-T similar to 10(8). Contrary to the popular notion, R-rho(t) first decrease during diffusion before the onset time and then increase when convection begins at the interface. Decrease in R-rho(t) is substantial at low Ra-T and it decreases even below unity resulting in overturning of the system. Interestingly, all the finger system passes through the same state before the onset of convection irrespective of Rayleigh number and density stability ratio of the system. (C) 2014 AIP Publishing LLC.

Modified Adomian Decomposition Method for Van der Pol equations

In the paper, the well known Adomian Decomposition Method (ADM) is modified to solve the parabolic equations. The present method is quite different than the numerical method. The results are compared with the existing exact or analytical method. The already known existing Adomian Decomposition Method is modified to improve the accuracy and convergence. Thus, the modified method is named as Modified Adomian Decomposition Method (MADM). The Modified Adomian Decomposition Method results are found to converge very quickly and are more accurate compared to ADM and numerical methods. MADM is quite efficient and is practically well suited for use in these problems. Several examples are given to check the reliability of the present method. Modified Adomian Decomposition Method is a non-numerical method which can be adapted for solving parabolic equations. In the current paper, the principle of the decomposition method is described, and its advantages are shown in the form of parabolic equations. (C) 2014 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/3.0/).

Regimes of nonlinear depletion and regularity in the 3D Navier-Stokes equations

The periodic 3D Navier-Stokes equations are analyzed in terms of dimensionless, scaled, L-2m-norms of vorticity D-m (1 <= m <= infinity). The first in this hierarchy, D-1, is the global enstrophy. Three regimes naturally occur in the D-1-D-m plane. Solutions in the first regime, which lie between two concave curves, are shown to be regular, owing to strong nonlinear depletion. Moreover, numerical experiments have suggested, so far, that all dynamics lie in this heavily depleted regime 1]; new numerical evidence for this is presented. Estimates for the dimension of a global attractor and a corresponding inertial range are given for this regime. However, two more regimes can theoretically exist. In the second, which lies between the upper concave curve and a line, the depletion is insufficient to regularize solutions, so no more than Leray's weak solutions exist. In the third, which lies above this line, solutions are regular, but correspond to extreme initial conditions. The paper ends with a discussion on the possibility of transition between these regimes.

Angular momentum transport and evolution of lopsided galaxies

The surface brightness distribution in the majority of stellar galactic discs falls off exponentially. Often what lies beyond such a stellar disc is the neutral hydrogen gas whose distribution also follows a nearly exponential profile at least for a number of nearby disc galaxies. Both the stars and gas are commonly known to host lopsided asymmetry especially in the outer parts of a galaxy. The role of such asymmetry in the dynamical evolution of a galaxy has not been explored so far. Following Lindblad's original idea of kinematic density waves, we show that the outer part of an exponential disc is ideally suitable for hosting lopsided asymmetry. Further, we compute the transport of angular momentum in the combined stars and gas disc embedded in a dark matter halo. We show that in a pure star and gas disc, there is a transition point where the free precession frequency of a lopsided mode, Omega - kappa, changes from retrograde to prograde and this in turn reverses the direction of angular momentum flow in the disc leading to an unphysical behaviour. We show that this problem is overcome in the presence of a dark matter halo, which sets the angular momentum flow outwards as required for disc evolution, provided the lopsidedness is leading in nature. This, plus the well-known angular momentum transport in the inner parts due to spiral arms, can facilitate an inflow of gas from outside perhaps through the cosmic filaments.

Co-Evolution of Content Spread and Popularity in Mobile Opportunistic Networks

We consider a setting in which a single item of content is disseminated in a population of mobile nodes by opportunistic copying when pairs of nodes come in radio contact. The nodes in the population may either be interested in receiving the content (referred to as destinations) or not yet interested in receiving the content (referred to as relays). We consider a model for the evolution of popularity, the process by which relays get converted into destinations. A key contribution of our work is to model and study the joint evolution of content popularity and its spread in the population. Copying the content to relay nodes is beneficial since they can help spread the content to destinations, and could themselves be converted into destinations. We derive a fluid limit for the joint evolution model and obtain optimal policies for copying to relay nodes in order to deliver content to a desired fraction of destinations, while limiting the fraction of relay nodes that get the content but never turn into destinations. We prove that a time-threshold policy is optimal for controlling the copying to relays, i.e., there is an optimal time-threshold up to which all opportunities for copying to relays are exploited, and after which relays are not copied to. We then utilize simulations and numerical evaluations to provide insights into the effects of various system parameters on the optimally controlled co-evolution model.

Combinatorial effect of rolling and carbonaceous nanoparticles on the evolution of crystallographic texture and structural properties of ultra high molecular weight polyethylene

Ultra high molecular weight polyethylene (PE) is a structural polymer widely used in biomedical implants. The mechanical properties of PE can be improved either by controlled crystalline orientation (texture) or by the addition of reinforcing agents. However, the combinatorial effect has not received much attention. The objective of this study was to characterize the structure and mechanical properties of PE composites incorporating multiwall carbon nanotubes (MWCNT) and reduced graphene oxide (RGO) subjected to hot rolling. The wide angle X-ray diffraction studies revealed that mechanical deformation resulted in a mixture of orthorhombic and monoclinic crystals. Furthermore, the presence of nanoparticles resulted in lower crystallinity in PE with smaller crystallite size, more so in RGO than in MWCNT composites. Rolling strengthened the texture of both orthorhombic and the monoclinic phases in PE. Presence of RGO weakened the texture of both phases of PE after rolling whereas MWCNT only mildly weakened the texture. This resulted in a reduction in the elastic modulus of RGO composites whereas moduli of neat polymer and the MWCNT composite increased after rolling. This study provides new insight into the role of nanoparticles in texture evolution during polymer processing with implications for processing of structural polymer composites.

Transition from dissipative to conservative dynamics in equations of hydrodynamics

We show, by using direct numerical simulations and theory, how, by increasing the order of dissipativity (alpha) in equations of hydrodynamics, there is a transition from a dissipative to a conservative system. This remarkable result, already conjectured for the asymptotic case alpha -> infinity U. Frisch et al., Phys. Rev. Lett. 101, 144501 (2008)], is now shown to be true for any large, but finite, value of alpha greater than a crossover value alpha(crossover). We thus provide a self-consistent picture of how dissipative systems, under certain conditions, start behaving like conservative systems and hence elucidate the subtle connection between equilibrium statistical mechanics and out-of-equilibrium turbulent flows.

Active guests in the MoS2/MoSe2 host lattice: efficient hydrogen evolution using few-layer alloys of MoS2(1-x)Se2x

Few-layer transition metal dichalcogenide alloys based on molybdenum sulphoselenides MoS2(1-x)Se2x] possess higher hydrogen evolution (HER) activity compared to pristine few-layer MoS2 and MoSe2. Variation of the sulphur or selenium content in the parent dichalcogenides reveals a systematic structure-activity relationship for different compositions of alloys, and it is found that the composition MoS1.0Se1.0 shows the highest HER activity amongst the catalysts studied. The tunable electronic structure of MoS2/MoSe2 upon Se/S incorporation probably assists in the realization of high HER activity.

Adaptive evolution of a novel avian-origin influenza A/H7N9 virus

In China, the recent outbreak of novel influenza A/H7N9 virus has been assumed to be severe, and it may possibly turn brutal in the near future. In order to develop highly protective vaccines and drugs for the A/H7N9 virus, it is critical to find out the selection pressure of each amino acid site. In the present study, six different statistical methods consisting of four independent codon-based maximum likelihood (CML) methods, one hierarchical Bayesian (HB) method and one branch-site (BS) method, were employed to determine if each amino acid site of A/H7N9 virus is under natural selection pressure. Functions for both positively and negatively selected sites were inferred by annotating these sites with experimentally verified amino acid sites. Comprehensively, the single amino acid site 627 of PB2 protein was inferred as positively selected and it function was identified as a T-cell epitope (TCE). Among the 26 negatively selected amino acid sites of PB2, PB1, PA, HA, NP, NA, M1 and NS2 proteins, only 16 amino acid sites were identified to be involved in TCEs. In addition, 7 amino acid sites including, 608 and 609 of PA, 480 of NP, and 24, 25, 109 and 205 of M1, were identified to be involved in both B-cell epitopes (BCEs) and TCEs. Conversely, the function of positions 62 of PA, and, 43 and 113 of HA was unknown. In conclusion, the seven amino acid sites engaged in both BCEs and TCEs were identified as highly suitable targets, as these sites will be predicted to play a principal role in inducing strong humoral and cellular immune responses against A/H7N9 virus. (C) 2014 Elsevier Inc. All rights reserved.