A parameter uniform difference scheme for parabolic partial differential equation with a retarded argument

The paper deals with the analysis of a non-stationary parabolic partial differential equation with a time delay. The highest order derivative term is affected by the small parameter. This is precisely the case when the magnitude of the convective term becomes much larger compare to that of diffusion term. The solution of problem exhibits steep gradients in the narrow intervals of space and short interval of times. In these cases a dissipative loss turned out to be more complex. Even for the one spatial dimension and one temporal variable, not all difference scheme can capture these steep variation. Although the analysis is restricted to the model in one space dimension, the technique and comparison principles developed should prove useful in assessing the merits of numerical solution of other nonlinear model equations too.     

An optimal asymptotic-numerical method for convection dominated systems having exponential boundary layers

The convection dominated diffusion problems are studied. Higher order accurate numerical methods are presented for problems in one and two dimensions. The underlying technique utilizes a superposition of given problem into two independent problems. The first one is the reduced problem that refers to the outer or smooth solution. Stretching transformation is used to obtain the second problem for inner layer solution. The method considered for outer or degenerate problems are based on higher order Runge–Kutta methods and upwind finite differences. However, inner problem is solved analytically or asymptotically. The schemes presented are proved to be consistent and stable. Possible extensions to delay differential equations and to nonlinear problems are outlined. Numerical results for several test examples are illustrated and a comparative analysis is presented. It is observed that the method presented is highly accurate and easy to implement. Moreover, the numerical results obtained are not only comparable with the exact solution but also in agreement with the theoretical estimates.     

A Robust numerical approach for singularly perturbed time delayed parabolic partial differential equations

We consider the numerical solution of a initial boundary value problem with a time delay. The problem under consideration is singularly perturbed from the mathematical perspective. Assuming that the coefficients of the differential equation are smooth, we construct and analyze the finite difference method whose solutions converge pointwise at all points of the domain independently of the singular perturbation parameter. The method permits its extension to the case of adaptive meshes, which may be used to improve the solution. Numerical examples are presented to demonstrate the effectiveness of the method. The convergence obtained in practice satisfies the theoretical predictions.     

Redefining the transparency order

In this paper, we consider the multi-bit Differential Power Analysis (DPA) in the Hamming weight model. In this regard, we revisit the definition of Transparency Order (\(\mathsf {TO}\)) from the work of Prouff (FSE 2005) and find that the definition has certain limitations. Although this work has been quite well referred in the literature, surprisingly, these limitations remained unexplored for almost a decade. We analyse the definition from scratch, modify it and finally provide a definition with better insight that can theoretically capture DPA in Hamming weight model for hardware implementation with precharge logic. At the end, we confront the notion of (revised) transparency order with attack simulations in order to study to what extent the low transparency order of an s-box impacts the efficiency of a side channel attack against its processing. To the best of our knowledge, this is the first time that such a critical analysis is conducted (even considering the original notion of Prouff). It practically confirms that the transparency order is indeed related to the resistance of the s-box against side-channel attacks, but it also shows that it is not sufficient alone to directly achieve a satisfying level of security. Regarding this point, our conclusion is that the (revised) transparency order is a valuable criterion to consider when designing a cryptographic algorithm, and even if it does not preclude to also use classical countermeasures like masking or shuffling, it enables to improve their effectiveness.     

Uniform convergence analysis of hybrid numerical scheme for singularly perturbed problems of mixed type

In this article, we consider a class of singularly perturbed mixed parabolic‐elliptic problems whose solutions possess both boundary and interior layers. To solve these problems, a hybrid numerical scheme is proposed and it is constituted on a special rectangular mesh which consists of a layer resolving piecewise‐uniform Shishkin mesh in the spatial direction and a uniform mesh in the temporal direction. The domain under consideration is partitioned into two subdomains. For the spatial discretization, the proposed scheme is comprised of the classical central difference scheme in the first subdomain and a hybrid finite difference scheme in the second subdomain, whereas the time derivative in the given problem is discretized by the backward‐Euler method. We prove that the method converges uniformly with respect to the perturbation parameter with almost second‐order spatial accuracy in the discrete supremum norm. Numerical results are finally presented to validate the theoretical results.© 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1931–1960, 2014     

A Phase Field Approach for Martensitic Transformations and Crystal Plasticity

This work is motivated by cryogenic turning which allows end shape machining and simultaneously attaining a hardened surface due to deformation induced martensitic transformations. To study the process on the microscale, a multivariant phase field model for martensitic transformations in conjunction with a crystal plastic material model is introduced. The evolution of microstructure is assumed to follow a time-dependent Ginzburg-Landau equation. To solve the field equations the finite element method is used. Time integration is performed with Euler backward schemes, on the global level for the evolution equation of the phase field, and on the element level for the crystal plastic material law. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Empirical process approach to some two-sample problems based on ranked set samples

We study the asymptotic properties of both the horizontal and vertical shift functions based on independent ranked set samples drawn from continuous distributions. Several tests derived from these shift processes are developed. We show that by using balanced ranked set samples with bigger set sizes, one can decrease the width of the confidence band and hence increase the power of these tests. These theoretical findings are validated through small-scale simulation studies. An application of the proposed techniques to a cancer mortality data set is also provided. Part of the research was conducted while Kaushik Ghosh was visiting Statistical Research and Applications Branch of the National Cancer Institute on an Intergovernmental Personnel Assignment.     

Failed disk recovery in double erasure RAID arrays

Reliability is a major concern in the design of large disk arrays. In this paper, we examine the effect of encountering more failures than that for which the RAID array was initially designed. Erasure codes are incorporated to enable system recovery from a specified number of disk erasures, and strive beyond that threshold to recover the system as frequently, and as thoroughly, as is possible. Erasure codes for tolerating two disk failures are examined. For these double erasure codes, we establish a correspondence between system operation and acyclicity of its graph model. For the most compact double erasure code, the full 2-code, this underlies an efficient algorithm for the computation of system operation probability (all disks operating or recoverable).When the system has failed, some disks are nonetheless recoverable. We extend the graph model to determine the probability that d disks have failed, a of which are recoverable by solving one linear equation, b of which are further recoverable by solving systems of linear equations, and d−a−b of which cannot be recovered. These statistics are efficiently calculated for the full 2-code by developing a three variable ordinary generating function whose coefficients give the specified values. Finally, examples are given to illustrate the probability that an individual disk can be recovered, even when the system is in a failed state.