Medical, Pharma, Engineering, Science, Technology and Business

Faculty of Science, Department of Mathematics, Mansoura University, Egypt

- *Corresponding Author:
- Sohalya MA

Faculty of Science, Department of Mathematics

Mansoura University, Egypt

**Tel:**+20 50 2383781

**E-mail:**m_stat2000@yahoo.com; mtyassen@yahoo.com; islamelbaz88@gmail.com

**Received date:** January 28, 2017; **Accepted date:** March 14, 2017; **Published date:** March 21, 2017

**Citation: **Sohalya MA, Yassena MT, Elbaza IM (2017) Stochastic Heat and Advection Cauchy. J Appl Computat Math 6:345. doi: 10.4172/2168-9679.1000345

**Copyright:** © 2017 Sohalya MA, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Applied & Computational Mathematics

In this paper, the solutions of Cauchy problems for the stochastic advection and stochastic diffusion equations are obtained using the finite difference method. In the case when the flow velocity is a function of stochastic flow velocity and also, the diffusion coefficient in the stochastic heat equation is a function of stochastic diffusion coefficient, the consistency and stability of the finite difference scheme we are used need to be performed under mean square calculus.

Stochastic differential equation; Mean square; Cauchy problem

Diffusion and advection problems involve randomness not only from the error measurements but may be also due to material defects or from impurities. The basic motivation of this paper focuses on this fundamental difficulty: the numerical consistency and stability of finite difference discretization of a Cauchy problems for stochastic advection and stochastic diffusion equations in the forms:

(1)

(2)

The independent variables are *t* (time) and *x* (space). The dependent variable is *U* which is a function of *t* and *x. p(x,ß)* is a discrete stochastic process (s.p.) defined on a common probability space (*Ω,F,P*) and depending of the stochastic variable *ß,(U0(x)),f(x)* are an initial deterministic data function.

In order to study the consistency and stability for the problems eqns (1) and (2) with the according difference schemes, it is necessary to use mean square calculus. Many papers have been studied the analytical solutions for stochastic differential models [1-7]. Also, in recent years numerical methods have also been applied to Cauchy problems such as heat equation with stochastic variable input [8]. Our paper is organized as follows, Section 2, deals with preliminaries of some points used in our discussing. Section 3, deals with constructing the stochastic finite difference scheme of eqns (1) and (2). Section 4, deals with the study of stability and consistency of (SFDS) of eqns (1) and (2).

Definition 1: A real random variable X on a probability space (Ω,F,P) and satisfying the property that is called 2^{nd} order random variable (*2.r.v*) where, E[ ] denotes the expectation value operator [9].

Remark 1:

If X∈L_{2}(Ω), Then the L_{2} norm is defined as

[9]

The approximations of the solution s.p. to eqns (1) and (2), [8], will be constructed in the sense of fixed station for the time, we will work in the following Banach space defined by

(3)

**The SFDS for stochastic cauchy diffusion eqn (1)**

The discretization of the problem eqn (1) will lead to stochastic difference scheme where at each node the unknown represents a r.v. Let us subdivide the domain (–∞,∞) × [0,∞] and define the grid cells for the space to be *Δx=(x _{k}–x_{k–1})* for

Using the first-order forward difference formula for approximating *U _{t}*

Then,

(4)

Using the first-order backward difference formula for approximating *p'(x,ß)*, *U _{x}* and the second-order centred difference formula for approximating

(5)

Substituting eqns (4) and (5) into eqn (1), we can motivate the scheme

Then,

(6)

**The SFDS for stochastic cauchy advection eqn (2)**

Now, we will construct the stochastic finite difference scheme by the same approach we use in the diffusion equation. Using eqn (4) for approximating *U _{t}* and first-order backward difference formula for approximating

(7)

Substituting eqns (4) and (7) into eqn (2), we get

Then,

(8)

**Consistency and stability of SFDS eqn (6)**

The consistency of numerical scheme is determined by studying the truncation error of the scheme. A finite difference scheme is consistent with the stochastic partial differential equation if the difference between the finite difference equation and the stochastic partial differential equation vanishes as the size of the grid spacing goes to zero independently. If this condition is met, then the numerical scheme is said to be consistent. Consistency of finite difference schemes is determined by using the technique of modified stochastic partial differential equation (MSPDE). This MSPDE is determined by expressing each term in the finite difference equation in a Taylor series at some base point. If the grid points approaches zero and the MSPDE changes back to original SPDE, then the scheme is consistent.

**Definition 2:** The stochastic finite difference scheme (SFDS)

is said to be mean square –consistent with the stochastic partial differential equation (RPDE) *Lu=F*, if the solution stochastic process (s.p.) of (RPDE), *U*, satisfies

(9)

**Theorem 1:** Let us consider the stochastic cauchy problem (1), and assume that

(1) *p(x)* and *p'(x)* are functions of the 2-r.v. ß such that *0<m<p(x,ω)<M*; m,M are positive constants.

(2) *U _{xxx}(x,t)(ω), U_{tt}(x,t)(ω) *are uniformly bounded for

**Proof:** Assuming that for each ω∈Ω, the solution U*(x,t)(ω) o*f (1) admits taylor expansion of the form

(10)

(11)

(12)

According to expression (9) with G=0, and using the expansions (10)-(12), we have

Since *p(x,ω)* is mean square differentiable, then we have

(13)

Assuming, *U _{xxx}(x,t)(ω), U_{tt}(x,t)(ω)* are uniformly bounded for

Then, from eqn (9) we get

Then,

*f _{ß}(ß)* is the probability density function of the stochastic variable. Hence, the SFDS (6) is mean square consistent.

A finite difference scheme is said to be stable if it produces a bounded solution for a stable stochastic partial differential equation and it is unstable if it produces an unbounded solution of the stochastic partial differential equations. There are several methods used to analyse stability of numerical schemes, for example the discrete perturbation method, the Von-Neumann method and the matrix method. In this work, we use the mean square calculus for stability analysis of numerical scheme.

**Definition 3:** A stochastic difference scheme is mean square stable, if there exist some positive constants ε,δ, non-negative constants *l,m* and *u ^{0}* is an initial data such that:

(14)

for all,

**Theorem 2:** The SFDS (6) that according to the problem (1) is mean square stable under the condition: , *0<m<p(x,ω)<M* where *m,M* are positive constants.

**Proof:** Since, our scheme is

Under the hypotheses , 0<Δ*t*≤δ and 0*<m<p(x,ω)<M, *we have Then,

(15)

Hence,

Then,

Finally, we have:

Therefore,

Hence, the SFDS (6) is mean square stable with l=1, m=0 and under the conditions eqn (15).

**Consistency and stability of SFDS (8)**

Taking into account definition eq (2), we can state that theorem:

**Theorem 3:** The SFDS (8) is mean square -consistent Let us consider the random cauchy problem (2), and assume that

(1) *p(x) and p'(x)* are functions of the 2-r.v. ß such that *0<m<p(x,ω)<M; m*,M are positive constants.

(2) *U _{xxx}(x,t)(ω), U_{tt}(x,t)(ω)* are uniformly bounded for

**Proof:** Assuming that for each ω∈Ω, the solution U*(x,t)(ω) o*f (2) admits taylor expansion of the form

(16)

(17)

According to expression (10) with G=0, and using the expansions (16) and (17), we have

Since *p(x,ω)* is mean square differentiable, then we have

(18)

Assuming, *U _{xxx}(x,t)(ω), U_{tt}(x,t)(ω)* are uniformly bounded for

Then, from eqn (9) we get

Then,

*f _{ß}(ß)* is the probability density function of the stochastic variable. Hence, the SFDS (8) is mean square consistent.

Now, taking into account definition (3), we can state that theorem

Theorem 4: The SFDS (8) that according to the problem (2) is mean square stable under the condition: , 0*<m<p(x,ω)<M* where *m,M* are positive constants.

**Proof:** Since, our scheme is

Under the hypotheses , 0<Δ*t*≤δand 0<*m<p(x,ω)<M, *we have for 0<*rM*≤1 Then,

(19)

Hence,

Then,

Finally, we have:

Therefore,

Hence, the SFDS (8) is mean square stable with l=1, m=0 and under the conditions eqn (19).

- Kloeden PE, Platen E (1995) Higher-order implicit strong numerical schemes for stochastic differential equations. Journal of Statistical Physics 66: 283-314.
- Yassen MT, Sohaly MA, Elbaz I (2016) Random Crank-Nicolson Scheme for Random Heat Equation in Mean Square Sense. American Journal of Computational Mathematics 6: 66.
- Cortés JC, Sevilla-Peris P, Jódar L (2005) Analytic-numerical approximating processes of diffusion equation with data uncertainty. Computers & Mathematics with Applications 49: 1255-1266.
- El-Tawil MA, Sohaly MA (2011) Mean square numericabl methods for initial value random differential equations. Open Journal of Discrete Mathematics 1: 66.
- Bishehniasar M, Soheili AR (2013) Approximation of stochastic advection-diffusion equation using compact finite difference technique. Iranian Journal of Science and Technology 37: 327-333.
- Wan X, Xiu D, Karniadakis GE (2004) Stochastic solutions for the two-dimensional advection-diffusion equation. SIAM Journal on Scientific Computing 26: 578-590.
- Villafuerte L, Cortés JC (2013) Solving Random Differential Equations by Means of Differential Transform Methods. In Advances in Dynamical Systems and Applications. Research India Publications 8: 413-425.
- Navarro-Quiles A, Romero JV, Roselló MD, Sohaly MA (2016) Approximating the solution stochastic process of the random Cauchy one-dimensional heat model. Abstract and Applied Analysis 2016: 5391368.
- Soong TT, Bogdanoff JL (1974) Random Differential Equations in Science and Engineering. Journal of Applied Mechanics 41: 1148.

Select your language of interest to view the total content in your interested language

- Adomian Decomposition Method
- Algebraic Geometry
- Analytical Geometry
- Applied Mathematics
- Axioms
- Balance Law
- Behaviometrics
- Big Data Analytics
- Binary and Non-normal Continuous Data
- Binomial Regression
- Biometrics
- Biostatistics methods
- Clinical Trail
- Complex Analysis
- Computational Model
- Convection Diffusion Equations
- Cross-Covariance and Cross-Correlation
- Differential Equations
- Differential Transform Method
- Fourier Analysis
- Fuzzy Boundary Value
- Fuzzy Environments
- Fuzzy Quasi-Metric Space
- Genetic Linkage
- Hamilton Mechanics
- Hypothesis Testing
- Integrated Analysis
- Integration
- Large-scale Survey Data
- Matrix
- Microarray Studies
- Mixed Initial-boundary Value
- Molecular Modelling
- Multivariate-Normal Model
- Noether's theorem
- Non rigid Image Registration
- Nonlinear Differential Equations
- Number Theory
- Numerical Solutions
- Physical Mathematics
- Quantum Mechanics
- Quantum electrodynamics
- Quasilinear Hyperbolic Systems
- Regressions
- Relativity
- Riemannian Geometry
- Robust Method
- Semi Analytical-Solution
- Sensitivity Analysis
- Smooth Complexities
- Soft biometrics
- Spatial Gaussian Markov Random Fields
- Statistical Methods
- Theoretical Physics
- Theory of Mathematical Modeling
- Three Dimensional Steady State
- Topology
- mirror symmetry
- vector bundle

- Total views:
**241** - [From(publication date):

March-2017 - Aug 17, 2017] - Breakdown by view type
- HTML page views :
**205** - PDF downloads :
**36**

Peer Reviewed Journals

International Conferences 2017-18