Environmental Engineering Reference
In-Depth Information
γ
1
s gn
2
k
1
1.5
0.5
s gn
1
1
s gn
0.1
2
Figure 5.19. Dispersion relation for prototype model ( 5.34 ), with spatial coupling a
la Swift-Hohenberg for three values of noise intensity ( a
=−
1, D
=
2, and k 0 =
1).
stable state
s c .
Mathematically, we can understand the short-term instability by observing that when
φ
φ
0 tend to disappear in spite of their initial amplification when s gn
>
is close to zero the disturbance effect that is due to the (multiplicative) noise tends
to prevail on the restoring effect of f . Conversely, as
3 of
φ
grows, the leading term
φ
0 .
Once the presence of a short-term instability has been detected, the capability of
spatiotemporal stochastic model ( 5.34 ) to generate patterns can be investigated, for
example, by use of the stability analysis by normal modes (see Box 5.1). In the case
of prototype model ( 5.34 ) we obtain from Eq. ( B5.1-4 ) a dispersion relation
f (
φ
) prevails over g and the deterministic local kinetics f tends to restore the state
φ
D ( k 0
k 2 ) 2
γ
( k )
=
a
+
s gn
,
(5.39)
which provides the same threshold s c as Eq. ( 5.36 ) for neutral stability; the maximum
amplification is for the wave number k
=
k 0 (see Fig. 5.19 ). It follows that statistically
steady periodic patterns, with wavelength
λ =
2
π/
k 0 , emerge when the noise intensity
exceeds the threshold s c =−
a .
The analysis of the structure function can be used as a prognostic tool to confirm
these results. The linear evolution of the structure function from Eq. ( B5.2-6 )is
2 a
k 2 ) 2 S ( k
S ( k
,
t )
D ( k 0
=
+
s gn
,
t )
.
(5.40)
t
Thus neutral stability corresponds to a critical value of the noise intensity:
D ( k 0
k 2 ) 2
s c
=−
a
+
.
(5.41)
When s gn <
s c , the structure function tends to zero and no patterns occur. Conversely,
when threshold ( 5.41 ) is exceeded, the linear evolution of the structure function is
divergent and a range of wave numbers becomes unstable. Notice that the wave
number most prone to instability is k
a . Therefore
these results are consistent with those obtained from the stability analysis.
=
k 0 , which corresponds to s c
=−
 
Search WWH ::




Custom Search