Strictly Stable High Order Di erence Approximations for ... - NO

mustardarchaeologistMechanics

Feb 22, 2014 (3 years and 8 months ago)

76 views

￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿
￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿
￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿
￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿
￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿ ￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿
￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿￿￿￿
￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿
￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿
￿￿￿￿￿￿￿ ￿￿￿ ￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿
￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿
￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿
￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿
￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿ ￿￿
￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿
￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿
￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿
￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿
￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿ ￿￿ ￿￿￿￿￿
￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿
￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿
￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿
￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿
￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿
￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿
U
0
t
+ F
0
x
+G
0
y
= S;￿￿￿
￿￿￿￿￿
U
0
=
0
B
B
@
½
0
(½u)
0
(½v)
0
(½E)
0
1
C
C
A
;
F
0
=
0
B
B
@
(½u)
0
(½u)
0
u
0
+p
0
(½v)
0
u
0
(½H)
0
u
0
+(½H)
0
u
0
1
C
C
A
;
G
0
=
0
B
B
@
(½v)
0
(½u)
0
v
0
(½v)
0
v
0
+p
0
(½H)
0
v
0
+(½H)
0
v
0
1
C
C
A
;
S = ¡
1
y
0
B
B
@
(½v)
0
(½u)
0
v
0
(½v)
0
v
0
(½H)
0
v
0
+(½H)
0
v
0
1
C
C
A
:
￿￿￿￿
½
0
= ½ ¡½
0
;(½u)
0
= ½u;u
0
=
(½u)
0
½
0

0
;
(½E)
0
= ½E ¡(½E)
0
;(½H)
0
= (½E)
0
+p
0
;
p
0
= (° ¡1)[(½E)
0
¡
1
2
(½u)
0
¢ u
0
]:
t ￿￿ ￿￿￿￿￿ ￿￿￿ x ￿￿￿ y ￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿
￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ½ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ u ￿￿￿ v ￿￿￿
x￿ ￿￿￿ y￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ E ￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿
￿￿￿￿￿ ￿￿￿ p ￿￿￿ ￿￿￿￿￿￿￿￿￿ ½
0
￿ (½E)
0
￿￿￿ (½H)
0
￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿
￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿
x = x(»;´)￿ y = y(»;´)￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿
￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿
￿￿ ￿￿￿￿
￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ S ´ 0 ￿￿
￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿ ￿￿￿
￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ U
0
= 0￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿
￿
￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿ ￿￿￿￿ ￿￿￿￿
￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿
￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿ ￿￿￿ ￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿
￿￿￿￿￿￿￿￿ ￿￿ ½
0
￿ (½u)
0
￿ (½E)
0
￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿
(½v)
0
￿ ￿￿ ￿ ￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿
￿￿ ￿￿￿￿￿
￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ [0;1] ￿￿ ￿￿￿￿￿￿￿￿￿￿￿
￿￿ N +1 ￿￿￿￿ ￿￿￿￿￿￿ x
m
= mh￿ m = 0;1;:::;N￿ ￿￿￿￿
h =
1
N
￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ Q ￿￿￿ ￿d=dx￿
￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿
(u;Qv)
h
= u
T
N
v
N
¡u
T
0
v
0
¡ (Qu;v)
h
;￿￿￿
￿￿￿￿￿ u;v 2 R
N+1
￿ ￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿
(u;v)
h
= hu
T
Hv;jjujj
2
h
= (u;u)
h
;
￿￿￿￿￿ H ￿￿ ￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ (N + 1) £
(N+1) ￿￿￿￿￿￿￿ ￿￿￿￿￿ H ￿￿ ￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿
￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿ ￿ ￿￿￿
￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ Q ￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿ ￿￿ ￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿
￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿
￿￿ H ￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿
￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿
￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿
￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿
￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿
￿￿ ￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿
￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿ ￿￿￿￿ ￿￿￿￿￿
￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿
￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿ ￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿ ￿ ￿￿￿￿￿
￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿
￿￿￿￿￿
1
h
l
X
j=1
®
j
(u
m+j
¡u
m¡j
)
￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿ du(x
m
)=dx ￿￿￿ ￿￿￿￿￿￿￿￿ ®
j
￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ 2(l ¡ 1)￿ ￿￿￿ ￿￿￿￿￿
￿￿￿￿￿￿ ®
j
￿￿￿￿ ￿￿￿￿￿￿￿ ￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿
˜
k ￿￿
h
˜
k = 2
l
X
j=1
®
j
sin(jhk)
￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿
Z
¼=2
¡¼=2
jhk ¡h
˜
kj
2
d(hk)
￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿
￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ h
˜
k = h
˜
k(hk) ￿￿ ￿￿￿￿￿￿￿￿
0
0.5
1
1.5
2
2.5
3
3.5
0
0.5
1
1.5
2
2.5
3
3.5
kh
approximate wave number
a
b
c
d
e
f
g
￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿
￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿
￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿ ￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿￿ ￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ 2¿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿
®
j
;j = 1:::¿ + 1￿ ￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ Q ￿￿￿￿￿ ￿￿
￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ¿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿ ￿￿￿￿￿￿
￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿ H ￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿
￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿
d
dx
￿￿￿ x = [0;1)￿
￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ Q ￿￿￿ ￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿
1;:::;2¿￿
Q ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿
hQ = H
¡1
µ
B C
¡C
T
D

;
￿￿￿￿￿ B ￿￿ ￿ ￿￿￿￿ 2¿ ￿￿ 2¿ ￿￿￿￿￿￿￿ C ￿￿ ￿ ￿￿￿￿￿ ￿￿￿￿￿￿
￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿ D ￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿
￿￿ ￿￿￿￿￿￿￿ HQ ￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿
h
00
q
00
= ¡1=2￿ ￿￿￿ Q ￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿ HQ ￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿
￿￿￿￿￿￿￿￿￿ ¿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ Q ￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿
￿￿￿￿ 0;:::;¿ ￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿
￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ H
￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ HQ￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿
￿￿￿￿ ￿￿ ￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿
￿￿ HQ￿ ￿￿￿￿￿￿￿ HQ ￿￿￿ H￿ Q ￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿
￿
￿￿￿￿￿￿￿
￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿
￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿
u
t
+
µ
1 0
0 ¡1

u
x
= 0;u =
µ
u
I
u
II

;
0 · x · 1;t ¸ 0;
￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿
u
I
(x;0) = sin2¼x;u
II
(x;0) = ¡sin2¼x;
￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿
u
I
(0;t) = u
II
(0;t);u
II
(1;t) = u
I
(1;t);t ¸ 0;
￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿
u
I
= sin2¼(x ¡t);u
II
= ¡sin2¼(x +t):
￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿
￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ u
I
￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿
￿￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿ t = 1:5￿
￿ ￿￿ ￿￿￿￿ ￿￿￿￿￿￿ ￿
￿￿￿ ￿￿￿
￿￿￿ ￿￿￿￿￿￿
￿￿￿
￿￿￿
￿￿￿￿￿￿
￿￿￿￿￿￿
￿￿￿
￿￿￿￿￿￿
￿￿￿￿￿￿
￿￿￿
￿￿￿￿￿￿
￿￿￿￿￿￿
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
-1
-0.5
0
0.5
1
1.5
x 10
5
x
Error
Error for SBP24: Time: 15 N=201
￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ N = 200 ￿￿ ￿￿￿￿
= 0:3￿
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
-12
-10
-8
-6
-4
-2
0
2
4
6
8
x 10
6
x
Error
Error for SBP24(6): Time: 15 N=201
￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ N = 200 ￿￿ ￿￿￿￿
= 0:3￿
￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿
￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿
￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿
￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿
￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿ ￿￿￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿
￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿
￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿ ￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿
￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿
￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿
￿￿￿ ￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿ ￿ ￿￿￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿ ￿￿￿￿￿
￿￿￿ ￿￿￿￿ ￿￿￿￿￿
￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿ ￿
￿￿ ￿￿￿￿￿￿￿￿ ￿ ￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿
￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿
￿￿￿￿￿ ￿￿ ￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿
￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿￿
￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿
￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿
￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿ ￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿
Axisymmetric boundary
With 3 rows of ghost points
Rocket
Duct
Farfield boundary
Solid wall boundary
Inlet/outlet boundary
Solid wall boundary
y
x
￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿
￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ x￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿ y = 0￿
￿
￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿ ￿￿￿￿ ￿￿ ￿￿￿￿ ￿￿ ￿￿￿
￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ 50m
￿￿￿ ￿￿￿ ￿￿￿￿￿￿ 5m￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿
￿￿ ￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿
￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿ 10m￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿
￿￿ ￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ u
0
="sin(!t)￿ ￿￿￿￿￿
"= 10
¡4
￿ ￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ · ￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿￿
￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿ u
0
￿￿
x = 0;y = 10m ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿ ￿￿ · ￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿
￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿
￿￿ ￿￿￿￿ ￿￿￿￿￿￿ 33Hz￿
￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿
￿￿￿￿￿￿ 240£240 ￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿ 38% ￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿
￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿￿
￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿ ￿￿ ￿￿￿ ￿￿￿￿
￿￿ ￿￿￿ ￿￿￿￿￿￿￿
10
20
30
40
50
60
70
80
90
100
110
-60
-40
-20
0
20
40
60
0
20
40
60
80
100
y
x
Acoustic pressure Around Ariane V
￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ p
0
￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿ ￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿
￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ 0￿
￿￿￿￿￿￿￿￿￿￿￿
￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿
￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿
￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿
￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿
￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿
￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿
￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿
￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿
￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿
￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿
￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿
￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿
￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿
￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿
￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿
￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿
￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿
￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿
￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿
￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿
￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿
￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿￿
￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿
￿￿￿ ￿￿￿￿￿
￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿
￿￿￿ ￿￿ ￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿
￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿
￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿
￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿ ￿￿￿
￿￿￿￿￿￿￿￿ ￿￿￿￿￿
￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿￿
￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿
￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿
￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿
￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿ ￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿
￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿
￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿ d=dx￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿
￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿
￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿
￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿
￿￿￿￿ ￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿
￿￿￿￿ ￿￿￿￿￿ ￿￿￿ ￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿￿
￿￿￿￿￿￿￿￿ ￿￿￿￿￿
￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿
￿￿￿￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿
￿￿￿￿￿￿￿ ￿￿￿￿￿
￿￿￿￿ ￿￿￿￿ ￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿￿ ￿￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿
￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿￿￿￿￿ ￿￿ ￿￿￿￿￿￿￿ ￿￿￿￿￿￿￿￿
￿￿￿￿ ￿￿￿￿ ￿￿￿ ￿￿￿￿￿￿ ￿￿￿￿￿
￿