Mathematical modelling of the oxidation of nickel titanium alloys

A one-dimensional bulk reaction model for the oxidation of nickeltitanium is formulated, with preferential oxidation of titaniumbeing included. The modelling is directed at the better understandingof the dominant mechanisms involved in the oxidation processand their significance for the biocompatibility of the alloy.Two different regimes for the relative diffusivities of oxygenand the metals are investigated. By assuming fast bulk reactions,different asymptotic structures emerge in different parameterregimes and the resulting models take the form of moving boundaryproblems. Different profiles of nickel concentration are obtained:in particular a nickel-rich layer (observed in practice) ispresent below the oxide/metal interface for the case when oxygenand the metals diffuse at comparable rates.     

Artificial time integration

Many recent algorithmic approaches involve the construction of a differential equation model for computational purposes, typically by introducing an artificial time variable. The actual computational model involves a discretization of the now time-dependent differential system, usually employing forward Euler. The resulting dynamics of such an algorithm is then a discrete dynamics, and it is expected to be “close enough” to the dynamics of the continuous system (which is typically easier to analyze) provided that small – hence many – time steps, or iterations, are taken. Indeed, recent papers in inverse problems and image processing routinely report results requiring thousands of iterations to converge. This makes one wonder if and how the computational modeling process can be improved to better reflect the actual properties sought. In this article we elaborate on several problem instances that illustrate the above observations. Algorithms may often lend themselves to a dual interpretation, in terms of a simply discretized differential equation with artificial time and in terms of a simple optimization algorithm; such a dual interpretation can be advantageous. We show how a broader computational modeling approach may possibly lead to algorithms with improved efficiency. AMS subject classification (2000)  65L05, 65M32, 65N21, 65N22, 65D18     

Sensitivity of depolarized lidar signals to cloud and aerosol particle properties

Measurements from depolarized lidars provide a promising method to retrieve both cloud and aerosol properties and a versatile complement to passive satellite-based sensors. For lidar observations of clouds and aerosols, multiple scattering plays an important role in the scattering process. Monte Carlo simulations are carried out to investigate the sensitivity of lidar backscattering depolarization to cloud and aerosol properties. Lidar parameters are chosen to be similar to those of the upcoming space-based CALIPSO lidar. Cases are considered that consist of a single cloud or aerosol layer, as well as a case in which cirrus clouds overlay different types of aerosols. It is demonstrated that besides thermodynamic cloud phase, the depolarized lidar signal may provide additional information on ice or aerosol particle shapes. However, our results show little sensitivity to ice or aerosol particle sizes. Additionally, for the case of multiple but overlapping layers involving both clouds and aerosols, the depolarized lidar contains information that can help identify the particle properties of each layer.     

Shape sensitivity of a plane crack front

The 3D‐elasticity model of a solid with a plane crack under the stress‐free boundary conditions at the crack is considered. We investigate variations of a solution and of energy functionals with respect to perturbations of the crack front in the plane. The corresponding expansions at least up to the second‐order terms are obtained. The strong derivatives of the solution are constructed as an iterative solution of the same elasticity problem with specified right‐hand sides. Using the expansion of the potential and surface energy, we consider an approximate quadratic form for local shape optimization of the crack front defined by the Griffith criterion. To specify its properties, a procedure of discrete optimization is proposed, which reduces to a matrix variational inequality. At least for a small load we prove its solvability and find a quasi‐static model of the crack growth depending on the loading parameter. Copyright © 2003 John Wiley & Sons, Ltd.     

Embedding of Dynamical Symmetry Groups U(1, 1) and U(2) of a Free Particle on AdS 2 and S 2 into Parasupersymmetry Algebra

Using two different types of the laddering equations realized simultaneously by the associated Gegenbauer functions, we show that all quantum states corresponding to the motion of a free particle on AdS ₂ and S ² are splitted into infinite direct sums of infinite-and finite-dimensional Hilbert subspaces which represent Lie algebras u(1, 1) and u(2) with infinite- and finite-fold degeneracies, respectively. In addition, it is shown that the representation bases of Lie algebras with rank 1, i.e., gl(2, C), realize the representation of nonunitary parasupersymmetry algebra of arbitrary order. The realization of the representation of parasupersymmetry algebra by the Hilbert subspaces which describe the motion of a free particle on AdS ₂ and S ² with the dynamical symmetry groups U(1, 1) and U(2) are concluded as well.     

Cold compaction of ceramic powder: Computational analysis of the effect of pressing method and die shape

This paper deals with a computational analysis of the influence of the pressing method and part geometry on the final density distribution in the cold compaction process of ceramic alumina powders. The analysis is based on the explicit finite-element model proposed and validated in a previous study. The mechanical behavior of the processing material is described using a multisurface elastoplastic model, the modified Drucker-Prager/Cap model Published in Prikladnaya Mekhanika, Vol. 43, No. 10, pp. 129–134, October 2007.     

Physicochemical studies on hydrophobic PEO‐PPO‐PEO triblock copolymer micelles in SDS micellar environment

The shape, size, aggregation, hydration, and correlation times of water insoluble PEO‐PPO‐PEO triblock copolymer micelles with sodium dodecylsulfate (SDS) micelles were investigated using transport studies and dynamic light scattering technique. From the conductance of micellar solutions of the polymer in 25 mM SDS and 5 mM NaCl, the hydration of polymer micelles were determined using the principle of obstruction of electrolyte migration by the polymer. The asymmetry of the micellar particles of polymer and polymer‐SDS mixed micellar systems in 5 mM NaCl and their average axial ratios were calculated using intrinsic viscosity and hydration data obeying Simha–Einstein equation. Hydration number and micellar sizes were variable with temperature. The shape of the polymer micelles has been ellipsoidal rather than spherical. The micellar volume, hydrodynamic radius, radius of gyration, diffusional coefficients as well as translational, rotational and effective correlation times have been calculated from the absolute values of the axes. The partial molal volume of polymer micelles has also been determined and its comparison with the molar volume of pure polymer suggested a volume contraction due to immobilization of the water phase by the hydrophilic head groups of the polymer. The thermodynamic activation parameters for viscous flow favor a more ordered water structure around polymer micelles at higher temperatures. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 2410–2420, 2007     

Identification of micro-scale calorimetric devices II

The miniaturized calorimetric devices furnish a reduced working flat surface and permits measurements with extremely low-mass quantities. The experimental sensitivity shows relevant position dependence with x-y surface coordinates and with z-distance. The device identification is realized via a 2-D model based in Fourier general equation. Using the Marquardt method the experimental flat surface device can be identified and the fitted parameters used to simulate the behavior of the experimental system. From the model, the effects of several dissipation configurations can be evaluated. Also, via the RC-analogy, a way to 3-D experimental devices is roughly described. This revised version was published online in July 2006 with corrections to the Cover Date.