Unified Galerkin- and DAE-Based Approximation of Fractional Order Systems

We consider numerical solutions of nonlinear multiterm fractional integrodifferential equations, where the order of the highest derivative is fractional and positive but is otherwise arbitrary. Here, we extend and unify our previous work, where a Galerkin method was developed for efficiently approximating fractional order operators and where elements of the present differential algebraic equation (DAE) formulation were introduced. The DAE system developed here for arbitrary orders of the fractional derivative includes an added block of equations for each fractional order operator, as well as forcing terms arising from nonzero initial conditions. We motivate and explain the structure of the DAE in detail. We explain how nonzero initial conditions should be incorporated within the approximation. We point out that our approach approximates the system and not a specific solution. Consequently, some questions not easily accessible to solvers of initial value problems, such as stability analyses, can be tackled using our approach. Numerical examples show excellent accuracy. DOI: 10.1115/1.4002516]

Activation barriers for d.c. conductivity in ionic glasses: Classical calculations using the small-cluster approximation

A small-cluster approximation has been used to calculate the activation barriers for the d.c. conductivity in ionic glasses. The main emphasis of this approach is on the importance of the hitherto ignored polarization energy contribution to the total activation energy. For the first time it has been demonstrated that the d.c. conductivity activation energy can be calculated by considering ionic migration to a neighbouring vacancy in a smali cluster of ions consisting of face-sharing anion polyhedra. The activation energies from the model calculations have been compared with the experimental values in the case of highly modified lithium thioborate glasses.

Influence of non-geometrical factors on intracrystalline diffusion -- Role of sorbate-zeolite interactions

Molecular dynamics simulations on Xe in NaY and Ar in NaCaA zeolite are reported. Rates of cage-to-cage crossovers in the two zeolites exhibit trends which are contrary to that expected from geometrical considerations. The results suggest the important role of the sorbate-zeolite interactions in determining the molecular sieve properties of zeolites for small sized sorbates. The results are explained in terms of the barrier height for cage-to-cage crossover in the two zeolites.

Phase diagram of the two-dimensional disordered Hubbard model in the Hartree-Fock approximation

Ground-state properties of the two-dimensional Hubbard model with point-defect disorder are investigated numerically in the Hartree-Fock approximation. The phase diagram in the p(point defect concentration)-delta(deviation from half filling) plane exhibits antiferromagnetic, spin-density-wave, paramagnetic, and spin-glass-like phases. The disorder stabilizes the antiferromagnetic phase relative to the spin-density-wave phase. The presence of U strongly enhances the localization in the antiferromagnetic phase. The spin-density-wave and spin-glass-like phases are weakly localized.

Sorbate properties and cage-to-cage diffusion of argon in NaCaA: a molecular dynamics study

Detailed molecular dynamics simulations of argon in zeolite NaCaA are reported. Thermodynamic, structural, and dynamical properties of the sorbate as a function of temperature have been obtained. The properties calculated include various site-site radial distribution functions, different energy distribution functions, selfdiffusion coefficients, the power spectra, and properties relating to cage-to-cage diffusion. The results suggest that sorbate is delocalized above 300 K. Both modes of cage-to-cage diffusion-the surface-mediated and centralized diffusion-are associated with negative barrier heights. Surprisingly, rate of cage-to-cage diffusion is associated with negative and positive activation energies below and above 500 K. The observed differences in the behavior of the rate of cage-to-cage diffusion between Xe-NaY and Ar-NaCaA systems and the nature of the potential energy surface are discussed. Presence of sorbatezeolite interactions results in significant enhancement in the rate of cage-to-cage diffusion and rate of cage visits. It is shown that properties dependent on the long-time behavior such as the diffusion coefficient and the rate of cages visited exhibit the expected Arrhenius dependence on temperature.

New Approximation Algorithms for Minimum Cycle Bases of Graphs

We consider the problem of computing an approximate minimum cycle basis of an undirected non-negative edge-weighted graph G with m edges and n vertices; the extension to directed graphs is also discussed. In this problem, a {0,1} incidence vector is associated with each cycle and the vector space over F-2 generated by these vectors is the cycle space of G. A set of cycles is called a cycle basis of G if it forms a basis for its cycle space. A cycle basis where the sum of the weights of the cycles is minimum is called a minimum cycle basis of G. Cycle bases of low weight are useful in a number of contexts, e.g. the analysis of electrical networks, structural engineering, chemistry, and surface reconstruction. Although in most such applications any cycle basis can be used, a low weight cycle basis often translates to better performance and/or numerical stability. Despite the fact that the problem can be solved exactly in polynomial time, we design approximation algorithms since the performance of the exact algorithms may be too expensive for some practical applications. We present two new algorithms to compute an approximate minimum cycle basis. For any integer k >= 1, we give (2k - 1)-approximation algorithms with expected running time O(kmn(1+2/k) + mn((1+1/k)(omega-1))) and deterministic running time O(n(3+2/k) ), respectively. Here omega is the best exponent of matrix multiplication. It is presently known that omega < 2.376. Both algorithms are o(m(omega)) for dense graphs. This is the first time that any algorithm which computes sparse cycle bases with a guarantee drops below the Theta(m(omega) ) bound. We also present a 2-approximation algorithm with expected running time O(M-omega root n log n), a linear time 2-approximation algorithm for planar graphs and an O(n(3)) time 2.42-approximation algorithm for the complete Euclidean graph in the plane.

Growth mechanism of phases by interdiffusion and diffusion of species in the niobium-silicon system

The integrated diffusion coefficient of the phases and the tracer diffusion coefficients of the species are determined in the Nb-Si system by the diffusion couple technique. The diffusion rate of Si is found to be faster than that of Nb in both the NbSi2 and Nb5Si3 phases. The possible atomic mechanism of diffusion is discussed based on the crystal structure and on available details of the defect concentration data. The faster diffusion rate of Si in the Nb5Si3 phase is found to be unusual. The growth mechanism of the phases is also discussed on the basis of the data calculated in this study. (C) 2010 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Decoupled three-dimensional finite element computation of thermoelastic damping using Zener's approximation

We consider three dimensional finite element computations of thermoelastic damping ratios of arbitrary bodies using Zener's approach. In our small-damping formulation, unlike existing fully coupled formulations, the calculation is split into three smaller parts. Of these, the first sub-calculation involves routine undamped modal analysis using ANSYS. The second sub-calculation takes the mode shape, and solves on the same mesh a periodic heat conduction problem. Finally, the damping coefficient is a volume integral, evaluated elementwise. In the only other decoupled three dimensional computation of thermoelastic damping reported in the literature, the heat conduction problem is solved much less efficiently, using a modal expansion. We provide numerical examples using some beam-like geometries, for which Zener's and similar formulas are valid. Among these we examine tapered beams, including the limiting case of a sharp tip. The latter's higher-mode damping ratios dramatically exceed those of a comparable uniform beam.

Stochastic Approximation Algorithms for Constrained Optimization via Simulation

We develop four algorithms for simulation-based optimization under multiple inequality constraints. Both the cost and the constraint functions are considered to be long-run averages of certain state-dependent single-stage functions. We pose the problem in the simulation optimization framework by using the Lagrange multiplier method. Two of our algorithms estimate only the gradient of the Lagrangian, while the other two estimate both the gradient and the Hessian of it. In the process, we also develop various new estimators for the gradient and Hessian. All our algorithms use two simulations each. Two of these algorithms are based on the smoothed functional (SF) technique, while the other two are based on the simultaneous perturbation stochastic approximation (SPSA) method. We prove the convergence of our algorithms and show numerical experiments on a setting involving an open Jackson network. The Newton-based SF algorithm is seen to show the best overall performance.

Can one electrochemically measure the statistical morphology of a rough electrode

We have developed a theory for an electrochemical way of measuring the statistical properties of a nonfractally rough electrode. We obtained the expression for the current transient on a rough electrode which shows three times regions: short and long time limits and the transition region between them. The expressions for these time ranges are exploited to extract morphological information about the surface roughness. In the short and long time regimes, we extract information regarding various morphological features like the roughness factor, average roughness, curvature, correlation length, dimensionality of roughness, and polynomial approximation for the correlation function. The formulas for the surface structure factors (the measure of surface roughness) of rough surfaces in terms of measured reversible and diffusion-limited current transients are also obtained. Finally, we explore the feasibility of making such measurements.

An actor-critic algorithm with function approximation for discounted cost constrained Markov decision processes

We develop in this article the first actor-critic reinforcement learning algorithm with function approximation for a problem of control under multiple inequality constraints. We consider the infinite horizon discounted cost framework in which both the objective and the constraint functions are suitable expected policy-dependent discounted sums of certain sample path functions. We apply the Lagrange multiplier method to handle the inequality constraints. Our algorithm makes use of multi-timescale stochastic approximation and incorporates a temporal difference (TD) critic and an actor that makes a gradient search in the space of policy parameters using efficient simultaneous perturbation stochastic approximation (SPSA) gradient estimates. We prove the asymptotic almost sure convergence of our algorithm to a locally optimal policy. (C) 2010 Elsevier B.V. All rights reserved.

Effect of sodium diffusion on the structural and electrical properties of Cu2ZnSnS4 thin films

Cu2ZnSnS4 (CZTS) is a kesterite semiconductor consisting of abundantly available elements. It has a band gap of 1.5 eV and a large absorption coefficient. Hence, thin films made of this material can be used as absorber layers of a solar cell. CZTS films were deposited on soda lime and Na free borosilicate glass substrates through Ultrasonic Spray Pyrolysis. The diffusion of sodium from soda lime glass was found to have a profound effect on characteristics like grain size, crystal texture and conductivity of CZTS thin films. Copper ion concentration also varied during the deposition and it was observed that the carrier concentration was enhanced when there was a deficiency of copper in the films. The effect of sodium diffusion and copper deficiency in enhancing the structural and electrical properties of CZTS films are presented in this paper. (C) 2010 Elsevier B.V. All rights reserved.

Diffusion Parameters and Growth Mechanism of Phases in the Cu-Sn System

The tracer diffusion coefficients of the elements as well as the integrated interdiffusion coefficients are determined for the Cu3Sn and Cu6Sn5 intermetallic compounds using incremental diffusion couples and Kirkendall marker shift measurements. The activation energies are determined for the former between 498 K and 623 K (225 A degrees C and 350 A degrees C) and for the latter between 423 K and 473 K (150 A degrees C and 200 A degrees C). Sn is found to be a slightly faster diffuser in Cu6Sn5, and Cu is found to be the faster diffuser in Cu3Sn. The results from the incremental couples are used to predict the behavior of a Cu/Sn couple where simultaneous growth of both intermetallics occurs. The waviness at the Cu3Sn/Cu6Sn5 interface and possible reasons for not finding Kirkendall markers in both intermetallics in the Cu/Sn couple are discussed.