17.03.2015 Views

DIFFERENtIAl & DIFFERENCE EqUAtIONS ANd APPlICAtIONS

DIFFERENtIAl & DIFFERENCE EqUAtIONS ANd APPlICAtIONS

DIFFERENtIAl & DIFFERENCE EqUAtIONS ANd APPlICAtIONS

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

A DISCRETE EIGENFUNCTIONS METHOD FOR NUMERICAL<br />

SOLUTION OF RANDOM DIFFUSION MODELS<br />

L. JÓDAR, J. C. CORTÉS, AND L. VILLAFUERTE<br />

This paper deals with the construction of numerical solution of random diffusion models<br />

whose coefficients functions and the source term are stochastic processes depending<br />

on a common random variable and an initial condition which depends on a different<br />

one. After discretization, the random difference scheme is solved using a random discrete<br />

eigenfunctions method. Mean-square stability of the numerical solution is studied, and<br />

a procedure for computing the expectation and the variance of the discrete approximate<br />

stochastic process is given.<br />

Copyright © 2006 L. Jódar et al. This is an open access article distributed under the Creative<br />

Commons Attribution License, which permits unrestricted use, distribution, and<br />

reproduction in any medium, provided the original work is properly cited.<br />

1. Introduction<br />

Mathematical models regarding spatial uncertainty are frequent in geostatistic description<br />

of natural variables [4] and modeling hydrology problems [5, 6]. Wave propagation<br />

in random media has been treated in [2] and fishering problems are modeled in [3] using<br />

stochastic processes. In this paper, we consider random diffusion models of the form<br />

u t = [ p(x,β)u x<br />

]x − q(x,β)u + F(x,t,β), 0 0, ∣∣ a 1<br />

∣ ∣ + ∣∣ a 2<br />

∣ ∣ > 0,<br />

b 1 u(1,t)+b 2 u x (1,t) = 0, t>0, ∣∣ a 1<br />

∣ ∣ + ∣∣ a 2<br />

∣ ∣ > 0,<br />

u(x,0)= f (x,γ), 0 ≤ x ≤ 1,<br />

(1.1)<br />

where the unknown u(x,t)aswellascoefficient p(x,β), the initial condition f (x,γ), and<br />

the source term F(x,t,β) are second-order stochastic processes depending on mutually<br />

independent second-order random variables β, andγ defined on a common probability<br />

space (Ω,,P).<br />

Model (1.1) assumes that random variations of the internal influences of the system<br />

undergoing diffusion are stochastic processes depending on the random variable β and<br />

that random external sources to the medium in which the diffusion takes places is also a<br />

Hindawi Publishing Corporation<br />

Proceedings of the Conference on Differential & Difference Equations and Applications, pp. 457–466

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!