TY - JOUR
T1 - Partially reflected diffusion
AU - Singer, A.
AU - Schuss, Z.
AU - Osipov, A.
AU - Holcman, D.
PY - 2007
Y1 - 2007
N2 - The radiation (reactive or Robin) boundary condition for the diffusion equation is widely used in chemical and biological applications to express reactive boundaries. The underlying trajectories of the diffusing particles are believed to be partially absorbed and partially reflected at the reactive boundary; however, the relation between the reaction constant in the Robin boundary condition and the reflection probability is not well defined. In this paper we define the partially reflected process as a limit of the Markovian jump process generated by the Euler scheme for the underlying Itô dynamics with partial boundary reflection. Trajectories that cross the boundary are terminated with probability P√Δt and otherwise are reflected in a normal or oblique direction. We use boundary layer analysis of the corresponding master equation to resolve the nonuniform convergence of the probability density function of the numerical scheme to the solution of the Fokker-Planck equation in a half-space, with the Robin constant k. The boundary layer equation is of the Wiener-Hopf type. We show that the Robin boundary condition is recovered if and only if trajectories are reflected in the conormal direction σn, where σ is the (possibly anisotropic) constant diffusion matrix and n is the unit normal to the boundary. Otherwise, the density satisfies an oblique derivative boundary condition. The constant k is related to P by K = rP√σn, where r = 1/√π and σn = n Tσn. The reflection law and the relation are new for diffusion in higher dimensions.
AB - The radiation (reactive or Robin) boundary condition for the diffusion equation is widely used in chemical and biological applications to express reactive boundaries. The underlying trajectories of the diffusing particles are believed to be partially absorbed and partially reflected at the reactive boundary; however, the relation between the reaction constant in the Robin boundary condition and the reflection probability is not well defined. In this paper we define the partially reflected process as a limit of the Markovian jump process generated by the Euler scheme for the underlying Itô dynamics with partial boundary reflection. Trajectories that cross the boundary are terminated with probability P√Δt and otherwise are reflected in a normal or oblique direction. We use boundary layer analysis of the corresponding master equation to resolve the nonuniform convergence of the probability density function of the numerical scheme to the solution of the Fokker-Planck equation in a half-space, with the Robin constant k. The boundary layer equation is of the Wiener-Hopf type. We show that the Robin boundary condition is recovered if and only if trajectories are reflected in the conormal direction σn, where σ is the (possibly anisotropic) constant diffusion matrix and n is the unit normal to the boundary. Otherwise, the density satisfies an oblique derivative boundary condition. The constant k is related to P by K = rP√σn, where r = 1/√π and σn = n Tσn. The reflection law and the relation are new for diffusion in higher dimensions.
KW - Markovian jump process
KW - Reactive boundary condition
KW - Stochastic differential equations
KW - Wiener-Hopf boundary layer equation
UR - http://www.scopus.com/inward/record.url?scp=47149106223&partnerID=8YFLogxK
U2 - 10.1137/060663258
DO - 10.1137/060663258
M3 - ???researchoutput.researchoutputtypes.contributiontojournal.article???
AN - SCOPUS:47149106223
SN - 0036-1399
VL - 68
SP - 844
EP - 868
JO - SIAM Journal on Applied Mathematics
JF - SIAM Journal on Applied Mathematics
IS - 3
ER -