Environmental Engineering Reference
In-Depth Information
Box 2.3: Direct derivation of the master equation for the processes
driven by WSN
In this box, we determine the master equation for process ( 2.57 ) driven by WSN, using
the approach by Cox and Miller ( 1965 ). We first transform generic multiplicative
Langevin equation ( 2.57 ) into an additive equation by using the transformation
z
= φ
φ )d
φ ,
1
/
g (
ρ
( z )
=
f [
φ
( z )]
/
g [
φ
( z )]. The Langevin equation for z is
d z
d t = ρ
( z )
+ ξ sn .
(B2.3-1)
In the infinitesimal interval d t the probability of having no jumps is (1
λ
d t )
+◦
(d t );
in this case, at time z
+
d t , the variable z will take the value
z ( t
+
d t )
=
z ( t )
z
,
(B2.3-2)
where
t + d t
z
=−
ρ
[ z (
τ
)]d
τ = ρ
[ z (
τ
)]d t
+◦
(d t )
,
(B2.3-3)
t
where
) is a high-order infinitesimal. The probability that a jump occurs in the same
interval d t is
(
·
λ
d t
+◦
(d t ). In this case,
z ( t
+
d t )
=
z ( t )
+
h
+
z
,
(B2.3-4)
where h is the size of the jump, with distribution p H ( h ). As a consequence, the
probability that the process takes a value in [ z
,
z
+
d z ] at time t
+
d t can be expressed as
p ( z
,
t
+
d t )d z
=
(1
λ
d t ) p ( z
+
z
,
t )d( z
+
z )
d t z
0
p ( z +
z ,
z )d( z +
z )d z
+ λ
t ) p H ( z
,
(B2.3-5)
where the second term on the right-hand side accounts for the case in which the process
reaches the z level because of a jump induced by the noise. Now, if expression ( B2.3-3 )
for
z is substituted into Eq. ( B2.3-5 ) and all terms of the order
(d t ) are neglected, we
obtain
,
+
=
λ
ρ
,
ρ
p ( z
t
d t )d z )
(1
d t ) p [ z
( z )d t
t ]d[ z
( z )d t ]
d t z
0
p [ z ρ
( z )d t
z )d[ z ρ
( z )d t ]d z
+ λ
,
t ] p H ( z
d t ) p ( z
t ) 1
( z )d t
( z )d t
=
(1
λ
,
t )
ρ
z p ( z
,
z ρ
d t d z z
0
p ( z ,
z )d z .
+ λ
t ) p H ( z
(B2.3-6)
Finally, dividing by d z , subtracting p ( z
,
t ) from both sides, dividing by d t , and taking
the limit as d t
0, we obtain (see Rodriguez-Iturbe et al. , 1999b ) the master equation:
z
p ( z
,
t )
=−
p ( z ,
z )d z .
z [ p ( z
,
t )
ρ
( z )]
λ
p ( z
,
t )
+ λ
t ) p H ( z
(B2.3-7)
t
0
Search WWH ::




Custom Search