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.     

Optical Tunnelling from One-Dimensional Square-Well Potentials

The authors study a one-dimensional model for optical tunnellingwith a refractive index in the shape of a square well. The relevanceof the model and its limiting cases are discussed. The mainresult is the leading behaviour of the exponentially small imaginarypart of the eigenvalue which determines the radiation loss.The leading behaviour of the imaginary part is calculated usingBerry's formula which controls the asymptotic expansion of theAiry function Bi(z) to better-than-exponential accuracy.     

Fully-developed pipe and planar flows of multimode viscoelastic fluids

Two solutions are presented for fully-developed pipe and planar flows of multimode viscoelastic models. The fluids have a Newtonian solvent contribution and the polymer modes are described by the Phan-Thien—Tanner (PTT), the FENE-P or the Giesekus equation. The first solution is exact and can handle any number of modes, but is only semi-analytical. The second solution, which is presented only for the PTT model with a linear stress coefficient and the FENE-P model, can also handle any number of modes. It is based on a truncated series expansion and is completely analytical, but provides only an approximated solution. The complexity of the multimode solutions is investigated first with the exact semi-analytical method and it is shown that at high Deborah number flows the high-order stresses can become as important as the stress of the first mode. It is also under these conditions that the approximated analytical solution deviates from the exact semi-analytical solution. A criterion for the accurate use of the approximated solution is presented. Fortran codes are provided to obtain these solutions at the internet address at the end.     

Die-swell,splashing drop and a numerical technique for solving the Oldroyd B model for axisymmetric free surface flows

This work deals with the development of a numerical method for simulating viscoelastic axisymmetric free surface flow of an Oldroyd B fluid. A novel formulation is developed for the computation of the non-Newtonian extra-stress components on rigid boundaries and on the symmetry axis. The full free surface stress conditions are employed. The resulting governing equations are solved by finite differences on a Marker-and-cell (MAC) type grid. Validation is provided by simulating a pipe flow problem. The classical die-swell problem is solved and swelling ratios are provided. The height of the splash caused by a falling liquid drop for various Reynolds and Weissenberg numbers is then studied, and the height of the splash is shown to diminish with increasing viscoelasticity.     

One-dimensional global attractor for strongly damped wave equations

We consider the dynamical behavior of the strongly damped wave equations under homogeneous Neumann boundary condition. By the property of limit set of asymptotic autonomous differential equations, we prove that in certain parameter region, the system has a one-dimensional global attractor, which is a periodic horizontal curve.     

A class of subgrid-scale models preserving the symmetry group of Navier–Stokes equations

Navier–Stokes equations (NS) admit transformations which transform a solution to another solution (galilean transformation, scaling transformation, …). They also admit viscosity dependent transformations which transform a solution to a solution of another NS with different viscosity. These particular transformations are called symmetries of NS. Each of them has a physical role (such as conservation laws, …). A consistent turbulence model should then remain invariant under these symmetry transformations. Unfortunately, this is not the case of several models.In this article, a class of subgrid-scale models preserving the symmetries of NS is built. This class is then refined such that the models respect the second law of thermodynamics. One of the simplest models of the class is tested to the flow in a ventilated room. Better results than those provided by Smagorinsky and dynamic models are obtained.     

A difference characteristic for one-dimensional deterministic systems

A numerical characteristic for one-dimensional deterministic systems reflecting its higher order difference structure is introduced. The comparison with Lyapunov exponent is given. A difference analogy for Eggleston theorem as well as an estimate for Hausdorff dimension of the difference attractor, formulated in terms of the new characteristic is proved.