Environmental Engineering Reference
In-Depth Information
p
φ
0
0.3
a
4
b
0.2
3
2
0.1
1
t
0.4
φ
1
5
10
0.1
0.2
0.3
φ
−
φ
0
obtained as an ensemble average of 10
6
real-
izations of model (
5.37
) with
a
Figure 5.18. (a) Behavior of
=−
1. The initial condition is
φ
=
0
.
1, and
s
gn
=
0
.
5
and
s
gn
=
5 for the lower and the upper curve, respectively. (b) Example of the
steady-state pdf of model (
5.38
)for
a
=−
1,
s
gn
=
5, and positive initial condition.
φ
=
In this case
f
eff
(
0
)
0 for any value of
s
gn
. The condition for marginal stability
is then
φ
0
=
d
f
eff
d
a
+
s
gn
=
0
,
(5.36)
φ
which provides the critical noise intensity
s
c
=−
a
. From condition (
B5.5-7
), it fol-
lows that
φ
0
is stable for
s
gn
<
−
a
and becomes unstable for
s
gn
>
−
a
. In the former
case, small disturbances around the steady state
φ
0
are dimmed whereas in the latter
they tend to grow exponentially for
t
0. However, this amplification is transient if
the spatial coupling is absent (see Fig.
5.17
): Fig.
5.18
(a) shows the temporal behav-
ior of the ensemble average of a number of numerically evaluated realizations of the
zero-dimensional stochastic model obtained when the spatial component is eliminated
from Eq. (
5.34
):
→
d
d
t
=
3
f
(
φ
)
+
g
(
φ
)
ξ
gn
(
t
)
=
a
φ
−
φ
+
φξ
gn
(
t
)
.
(5.37)
The results in Fig.
5.18
(a) show that the growth phase appears only when
s
gn
>
s
c
and
at the beginning of simulations (i.e., at short term), whereas the effect of the initial
perturbation disappears in the long term.
The steady-state pdf corresponding to model (
5.37
) can be found with the methods
described in Chapter 2:
2
2
s
gn
φ
a
s
gn
Ce
−
−
1
p
(
φ
)
=
φ
,
(5.38)
where
C
is a normalization constant and the domain of
p
(
φ
) is either the positive or
the negative semiaxis of
, depending on whether the initial condition is positive or
negative valued. In fact, the condition
g
(
φ
0 be a natural
boundary of the dynamics [an example of pdf is shown in Fig.
5.18
(b)]. Equation
(
5.38
) shows that the mode of the steady-state distribution is always at
φ
=
0)
=
0 entails that
φ
=
0forany
noise intensity. This confirms that in the long run perturbations of the deterministic
φ
=
Search WWH ::
Custom Search