Application of hpDGFEM to mechanisms at channel microband electrodes

We extend our earlier work (Harriman et al., Oxford University Computing Laboratory Technical Report NA04/19) on hp-DGFEM for disc electrodes to the case of reaction mechanisms to the increasingly popular channel microband electrode configuration. We present results for the simple E reaction mechanism (convection-diffusion equation), for the ECE and EC2E reaction mechanisms (linear and nonlinear systems of reaction-convection- diffusion equations, respectively) and for the DISP1 and DISP2 reaction mechanisms (linear and nonlinear coupled systems of reaction-convection-diffusion equations, respectively). In all cases we demonstrate excellent agreement with previous results using relatively coarse meshes and without the need for streamline-diffusion stabilisation, even at high flow rates

Finite-time singularities in the dynamical evolution of contact lines

We study finite-time singularities in the linear advection-diffusion equation with a variable speed on a semi-infinite line. The variable speed is determined by an additional condition at the boundary, which models the dynamics of a contact line of a hydrodynamic flow at a 180 contact angle. Using apriori energy estimates, we derive conditions on variable speed that guarantee that a sufficiently smooth solution of the linear advection--diffusion equation blows up in a finite time. Using the class of self-similar solutions to the linear advection-diffusion equation, we find the blow-up rate of singularity formation. This blow-up rate does not agree with previous numerical simulations of the model problem.Comment: 9 pages, 2 figure

Modelling diffusion in crystals under high internal stress gradients

Diffusion of vacancies and impurities in metals is important in many processes occurring in structural materials. This diffusion often takes place in the presence of spatially rapidly varying stresses. Diffusion under stress is frequently modelled by local approximations to the vacancy formation and diffusion activation enthalpies which are linear in the stress, in order to account for its dependence on the local stress state and its gradient. Here, more accurate local approximations to the vacancy formation and diffusion activation enthalpies, and the simulation methods needed to implement them, are introduced. The accuracy of both these approximations and the linear approximations are assessed via comparison to full atomistic studies for the problem of vacancies around a Lomer dislocation in Aluminium. Results show that the local and linear approximations for the vacancy formation enthalpy and diffusion activation enthalpy are accurate to within 0.05 eV outside a radius of about 13 Å (local) and 17 Å (linear) from the centre of the dislocation core or, more generally, for a strain gradient of roughly up to 6 × 10^6 m^-1 and 3 × 10^6 m^-1, respectively. These results provide a basis for the development of multiscale models of diffusion under highly non-uniform stress

Bistable reaction equations with doubly nonlinear diffusion

Reaction-diffusion equations appear in biology and chemistry, and combine linear diffusion with different kind of reaction terms. Some of them are remarkable from the mathematical point of view, since they admit families of travelling waves that describe the asymptotic behaviour of a larger class of solutions $0\leq u(x,t)\leq 1$ of the problem posed in the real line. We investigate here the existence of waves with constant propagation speed, when the linear diffusion is replaced by the "slow" doubly nonlinear diffusion. In the present setting we consider bistable reaction terms, which present interesting differences w.r.t. the Fisher-KPP framework recently studied in \cite{AA-JLV:art}. We find different families of travelling waves that are employed to describe the wave propagation of more general solutions and to study the stability/instability of the steady states, even when we extend the study to several space dimensions. A similar study is performed in the critical case that we call "pseudo-linear", i.e., when the operator is still nonlinear but has homogeneity one. With respect to the classical model and the "pseudo-linear" case, the travelling waves of the "slow" diffusion setting exhibit free boundaries. \\ Finally, as a complement of \cite{AA-JLV:art}, we study the asymptotic behaviour of more general solutions in the presence of a "heterozygote superior" reaction function and doubly nonlinear diffusion ("slow" and "pseudo-linear").Comment: 42 pages, 11 figures. Accepted version on Discrete Contin. Dyn. Sys

On the linear fractional self-attracting diffusion

In this paper, we introduce the linear fractional self-attracting diffusion driven by a fractional Brownian motion with Hurst index 1/2<H<1, which is analogous to the linear self-attracting diffusion. For 1-dimensional process we study its convergence and the corresponding weighted local time. For 2-dimensional process, as a related problem, we show that the renormalized self-intersection local time exists in L^2 if $\frac12<H<\frac3{4}$.Comment: 14 Pages. To appear in Journal of Theoretical Probabilit