Models of thermodiffusion in 1D

The goal of the paper is to study the Cauchy problem for 1D models of thermodiffusion. We explain qualitative properties of solutions. In particular, we show which part of the model has a dominant influence on well‐posedness, propagation of singularities, L^p ? L^q decay estimates on the conjugate line, and on the diffusion phenomenon. Copyright © 2013 John Wiley & Sons, Ltd.     

Spike patterns in a reaction–diffusion ODE model with Turing instability

We explore a mechanism of pattern formation arising in processes described by a system of a single reaction–diffusion equation coupled with ordinary differential equations. Such systems of equations arise from the modeling of interactions between cellular processes and diffusing growth factors. We focus on the model of early carcinogenesis proposed by Marciniak‐Czochra and Kimmel, which is an example of a wider class of pattern formation models with an autocatalytic non‐diffusing component. We present a numerical study showing emergence of periodic and irregular spike patterns because of diffusion‐driven instability. To control the accuracy of simulations, we develop a numerical code on the basis of the finite‐element method and adaptive mesh grid. Simulations, supplemented by numerical analysis, indicate a novel pattern formation phenomenon on the basis of the emergence of nonstationary structures tending asymptotically to a sum of Dirac deltas. Copyright © 2013 John Wiley & Sons, Ltd.     

ON INTERSECTIONS OF INDEPENDENT NONDEGENERATE DIFFUSION PROCESSES

Let X(1)= {X(1)(s), s ∈ R+ } and X(2)= {X(2)(t), t ∈ R+ } be two independent nondegenerate diusion processes with values in Rd. The existence and fractal dimension of intersections of the sample paths of X(1)and X(2)are studied. More generally, let E1, E2 ■(0, ∞) and F  Rd be Borel sets. A necessary condition and a suffcient condition for P{X(1)(E1) ∩ X(2)(E2) ∩ F = φ} 0 are proved in terms of the Bessel-Riesz type capacity and Hausdor measure of E1 ×E2 ×F in the metric space(R+ ×R+ ×Rd, ρ), where ρ is an unsymmetric metric defined in R+ × R+ × Rd. Under reasonable conditions, results resembling those of Browian motion are obtained.     

A compact 9 point stencil based on integrated RBFs for the convection–diffusion equation

In this paper, a high-order compact stencil for solving the convection–diffusion equation in two dimensions is proposed. The convection and diffusion terms are both approximated by means of radial basis functions (RBFs) that are constructed over 3×3$3\times 3$ rectangular stencils. Salient features here are that (i) integration is employed to construct local RBF approximations; and (ii) through the constants of integration, values of the convection–diffusion equation at several selected nodes on the stencil are also enforced. Numerical results indicate that (i) the inclusion of the governing equation into the stencil leads to a significant improvement in accuracy; (ii) when the convection dominates, accurate solutions are obtained at a regime of the RBF width which makes the RBFs peaked; and (iii) high levels of accuracy are achieved using relatively coarse grids.     

A new fractional finite volume method for solving the fractional diffusion equation

The inherent heterogeneities of many geophysical systems often gives rise to fast and slow pathways to water and chemical movement. One approach to model solute transport through such media is by fractional diffusion equations with a space–time dependent variable coefficient. In this paper, a two-sided space fractional diffusion model with a space–time dependent variable coefficient and a nonlinear source term subject to zero Dirichlet boundary conditions is considered.Some finite volume methods to solve a fractional differential equation with a constant dispersion coefficient have been proposed. The spatial discretisation employs fractionally-shifted Grünwald formulas to discretise the Riemann–Liouville fractional derivatives at control volume faces in terms of function values at the nodes. However, these finite volume methods have not been extended to two-dimensional and three-dimensional problems in a natural manner. In this paper, a new weighted fractional finite volume method with a nonlocal operator (using nodal basis functions) for solving this two-sided space fractional diffusion equation is proposed. Some numerical results for the Crank–Nicholson fractional finite volume method are given to show the stability, consistency and convergence of our computational approach. This novel simulation technique provides excellent tools for practical problems even when a complex transition zone is involved. This technique can be extend to two-dimensional and three-dimensional problems with complex regions.     

The inverse solution of the atomic mixing equations by an operator-splitting method

The quantification problem of recovering the original material distribution from secondary ion mass spectrometry (SIMS) data is considered in this paper. It is an inverse problem, is ill-posed and hence it requires a special technique for its solution. The quantification problem is essentially an inverse diffusion or (classically) a backward heat conduction problem. In this paper an operator-splitting method (that is proposed in a previous paper by the first author for the solution of inverse diffusion problems) is developed for the solution of the problem of recovering the original structure from the SIMS data. A detailed development of the quantification method is given and it is applied to typical data to demonstrate its effectiveness.     

An accurate application of the integral method applied to the diffusion of oxygen in absorbing tissue

Accurate integral methods are applied to a one dimensional moving boundary problem describing the diffusion of oxygen in absorbing tissue. These methods have been well studied for classic Stefan problems but this situation is unusual because there is no condition which contains the velocity of the moving boundary explicitly. This paper begins by giving a short time solution and then discusses some of the previous integral methods found in the literature. The main drawbacks of these solutions are that they cannot be solved from $t=0$ and also cannot determine the end behaviour. This is due to the non-uniform initial profile which integral methods typically fail to capture. The use of a novel transformation removes this non-uniformity and, on applying optimal integral methods to the resulting system, leads to simple and yet very accurate approximate solutions that overcome the deficiencies of previous methods.     

The roles of diffusivity and curvature in patterns on surfaces of revolution

We address the question of finding sufficient conditions for existence as well as nonexistence of nonconstant stable stationary solution to the diffusion equation _ut=div(a∇u)+f(u)$u_t=\mathrm{div}(a\mathrm{\nabla }u)+f(u)$ on a surface of revolution with and without boundary. Conditions found relate the diffusivity function a and the geometry of the surface where diffusion takes place. In the case where f is a bistable function, necessary conditions for the development of inner transition layers are given.     

A sufficient condition for the existence of a principal eigenvalue for nonlocal diffusion equations with applications

Considerable work has gone into studying the properties of nonlocal diffusion equations. The existence of a principal eigenvalue has been a significant portion of this work. While there are good results for the existence of a principal eigenvalue equations on a bounded domain, few results exist for unbounded domains. On bounded domains, the Krein–Rutman theorem on Banach spaces is a common tool for showing existence. This article shows that generalized Krein–Rutman can be used on unbounded domains and that the theory of positive operators can serve as a powerful tool in the analysis of nonlocal diffusion equations. In particular, a useful sufficient condition for the existence of a principal eigenvalue is given.     

Random attractors for stochastic reaction‐diffusion equations with multiplicative noise in

In this paper, we study the random dynamical system generated by a stochastic reaction‐diffusion equation with multiplicative noise and prove the existence of an ‐random attractor for such a random dynamical system. The nonlinearity f is supposed to satisfy some growth of arbitrary order .