close

Вход

Забыли?

вход по аккаунту

?

Эффективный алгоритм расчета квадратичных нелинейностей в спектральной аппроксимации для дифференциальных уравнений второго порядка.

код для вставкиСкачать
?╩╚╩═©
??? 532.526, 532.529
?.?. ?????, ?.?. ?????????
??????????? ???????? ??????? ????????????
????????????? ? ???????????? ?????????????
??? ???????????????? ????????? ??????? ???????
D.I. Popov, A.M. Sagalakov
The Effective Algorithm of Quadratic Nonlinearity
Calculation in Spectral Approximation
for the Second Order Differential Equation
??????????? ??????????? ????? ???????????? ????????????? ??? ????????????? ?????????? ?????? ?
???? ???????????? ????? ? ??????????? ????????.
??????????? ?????????? ????????????.
The special form of spectral approximation is concerned to represent nonliner terms as quadratic form with
sparse matrix. The testing results are discussed.
???????? ?????: ???????????? ?????????????,
????????? ????????, ????? ?????????, ????????????? ???????.
Key words: spectral approximation, Burgers equation,
Galerkin method, orthogonal system.
????????. ??????????? ???????????? ?????
??????? ????????-??????? ?????? ??? ?????????
??????? ??????? ? ?????? ???????????? ???????.
???????? ???????? ? ?????? ????????? ????????????? ??????????? ??????????. ? ?????? ?????????????? ??????? ????? ???????????? ??????? ?????????? ??????????????? ??????????? ?? ????????????? ? ???????? ?????????? ?????????? ?? ????????????? ??????? ????? ?????, ????????, ? ??????? [1, 2]. ?????? ???????? ?????????? ??????????? ??? ???????, ????????? ??????????????? ???
?????????????. ? ?????? ????????????????? ????????? ??????? ??????????? ??????????? ????? ???????? ???????, ????? ???????? ??????????? ??? ??????????????? ??????? ?????????? ?????????? ???????????? ????????. ????? ???????, ??? ????????????? ?????? ????? ???????? ????????????????
????????? ????????????? ??????, ???????????? ???
???????????? ???????.
????? ?????????????. ???????? [3], ??? ?????? ????????-??????? ?????? ????????? ??????????????? ???? ????????? ?? ???????????? W21,0 (QT ) ,
x?I ), ???????????? ? ?????? [4]. ??? ???????????
???????????? ???????, ??????????? ??????????? ?
????? ? ( x) ? (1? x 2 ) ?1/ 2 , W2,1? ( I ) ? L2,? ( I ) , ?????????-
?????. ????????, ????????, ??? L2 ( I ) ?????????
??????? ? L2,? ( I ) ? ?.?., ????????? ????????
{Tk ( x)}?0 ????? ? ???????? ????????????? ???????
? L2,? ( I ) [5]. ?????? ?????????? ? L2,? ( I ) ?? ??????????? ?????????? ? L2 ( I ) , ??????? ??????? ???
???????? ??????? ?????????? ???????. ? ??????,
???????? ?????????? ?????? ????????? ??? ??????? {?k }0N [4], ????? ??????? ??????? ????????
????????? H o1,? ( I ) ?W2,j? ( I ) , ??? j? 2 , H o1,? ( I ) ?
????????? ????????? ??????
?????? ??????, ????? ???????????? ????????????? ????? ??????? u(x,t), ??????? ???????? ?????????? W21 ( I ) , ?.?. ????? ?????? ??????????
??????????? ?? x, ? ?????????? ??????? ?? ???????. ????? ????? ????????????? ????? {?k }0N , ???
?k ( x) ?Tk ( x) ?Tk ? 2 ( x) ( Tk ( x) ? ??????? ????????,
? ?????
1
2,?
W ( I ) . ?????????? ??????? ??????? ? ????
N
u N ( x, t ) ? ? zk (t )?k ( x) .
k ?0
????? ???????? ??????????? ?????????? ????????? ???
N
bh (u N ,u N ) ? (u N Dx u N , ? h ) ? ? ? zk z s (? k ? s? , ? h )? .
??? QT ? I ?(0, T ) , I ? (?1,1) , ? ????? W21,0 (QT ) ?????????? ??????????? ???????????? ?? ?????????
?????????????
(u, v) ? ? (uv ? u x vx )dxdt .
QT
Co? ( I )
k ,s ?0
?????? ????????? ???????????:
1
Gskh ? (?k ? s? , ? h )? ? ? ?k ( x)? s? ( x)? h ( x)? ( x)dx .
?1
??????? ? ?????????? ??????? ????? ????????
? ????
?n ? cos n? ? cos(n ? 2)? ? 2sin ? sin(n ? 1)? ,
??? x ? cos ? , ? ?[0, ? ] .
????? ??????????? ?????
d? n
2sin(n ?1)?
??
cos ? ? 2(n ?1) cos(n ?1)? .
sin ?
dx
184
?ыыбмзк?л?х ?н?с?кзо ??р?бз? м?????зк?л?? лбнклбхлсрзбх...
?????????? ????????
sin( s ?1)?
Gskh (1) ? 4 ?
cos ? sin ? sin(h ?1)? (cos k? ? cos(k ? 2)? )d? .
sin ?
0
?
???????? ??????????, ???
Gskh (1) ? M skh (? ks??1h ?? ks??3h ?? ks??1h ? 2 ?? ks??3h ? 2 ) ,
??? M skh ? (? /2)ci , c0 ? 2 , ci ?1, i ? 0 , ? ???????? i ???????????? ??????????? h, k, s.
???????? ? ?????????? ?????????:
?
Gskh (2) ? 2( s ?1) ? cos( s ?1)? (cos k? ? cos( k ? 2)? )(cosh ? ?cos(h ? 2)? )d? .
0
?????? ???????
Gskh (2) ? ( s ?1) M skh (? ks ? h ?1 ?? ks ? h ?1 ?? ks ? h ? 3 ?? ks ? h ?1 )
? ( s ?1) M skh (? ks??2h ? 3 ?? ks??2h ?1 ?? ks??2h ?1 ?? ks??2h ?1 ).
????? ???????, ?????????? ????????? ???????????? ? E N ?1 ?????????? ?? N ?1 ???????????? ????,
? ????????? ?????????? h-? ????????, ??? ????????
????????, ???????? ????????? ??????? O( N ) .
?????????? ????????????. ??? ????????????
????????????? ????????????? ?????????? ???????
????????? ????????
? u / ? t ? u? u / ? y ??? 2 u / ? y 2 ? 0
? ????????-???????? ?????????
u (??, t ) ?1 , u (?, t ) ? 0 ; u ( y, 0) ?1
( ?? ? y ? 0 ), u ( y, 0) ? 0 ( 0 ? y ? ? ),
? ????????? ????????????? ? ?????? ???????? ?
???????????? ?????? [6]. ??? ????????, ??? ??????
????????? ??????????????? ?????? ??????? ?????.
??? ??????? ???????? ????? ?????????? ??????
????????????, ? ??? ????? ????????? ???? ????????, ??????????????? ???????? ?????? ??????????
???????. ??????, ???????? ?? ??????????? ????????, ??????? ??????????????? ????????? ? ????????? ????????, ??????? ???????? ????? ????????
???????. ????????? ????? H ?????? ????????, ?????
??????????????? ????????? y ? Hx ????? ????????
???????? ?????? ? ?????????? ????:
? u / ? t ? (1/ H )u? u / ? x ? (? / H 2 )? 2 u / ? x 2 ? 0 ,
R ?1/? , u (?1, t ) ?1 , u (1, t ) ? 0 ; u ( x, 0) ?1
( ?1? x ? 0 ), u ( x, 0) ? 0 ( 0? x ?1 ).
??? ?????????????? ?? ??????? ??????????????
????? ??????-?????????, ??????? ? ?? ????????
?????,
??????? ????????? ????????? ?????? ?? ??????? ????? ????????????????? ?? ??????? ???????
??? R = 10, t = 0.5 ??? N = 16. ??????, ?????????
?????? ?????????? ?????????? ????? ? ???????
?????? ?????????? ? ????? ??????-??????? ????????
???????? ??????? N + 1 ? ??????????? ?????? ????????? ? ????????? ??????. ???????????????
???????? ??????????????????? ?????????? ?????????: u1a ?ue ? 0.0013 ? u2 a ?ue ? 0.0008 , ?????? ??????????? ????????? ??????????? ??????? 1/420.
?? ??????? 1 ????????? ????????????? ???????
??????? ??? ????????? N, ?????, ??? ???????????
???????????? ?? ???? ??????? ? ??????????? ? ?????? ???????? N. ???????? ?????????? ??????????? ?????? ?? ??????????????, ??? ?????????????
??????? ?????? ???????? ?? ?????????? ???????
???????, ???????????? ??????? ??????? ?????.
??????? 1 ??????????, ??? ??? ??? N = 16 ?????
???????????? ?????????? ????????. ?????? ?????????? ?????????? N ????? ????????????? ?????????? ?? ?????????? ???????????. ? ??????? 2
???????????? ?????????? ????????? ?????????????
?????? ? ????????????? ?????? ????????? (??. [6])
??? ????????? ????? ??????????.
???. 1. ??????????? ???????
????????? ???????? ??? R = 10
??????? 1
??????? ????????? ????????
???????????? ??????? R = 10, t = 0.5
(? t ) ?1 (uhm ?uhm ?1 , ?h ) ? (? /2)[uhm ? uhm ?1 , ? h ]?
?bh (u m ?1 , u m ?1 , ? h ) ? 0,
? ??????? ??????? ??????? ? ????
X
N
RMS
u N ( x, t ) ? (1? x)/2 ? ? zk (t )?k ( x) .
???????????? ???????
N=5
0.0183
N=7
0.0079
N=9
0.0040
N = 16
0.0008
N = 32
0.0007
k ?0
??????? 2
????????? ??????? ????????? ???????? ? ???????????? [6]
X
RMS
*
N = 9*
0.0083
R = 10, t = 0.92
N=7
0.0058
N = 16
0.0013
? ???????, ??????????? ???????????? ??????? ????????? ? [6].
185
N = 9*
0.1049
R = 100, t = 0.92
N=9
0.0974
N = 104
0.0006
?╩╚╩═©
?????? ?????????????????? ??????????? ?? ?????? ?????? ??????? 2 ??????????, ??? ??? ?????????? ?????????? ???????????? ????? ???? ????????? ?????? ??????????, ???? ? ?????? [6] ????????? ?? ??????? ???????????? ?????? ?????-?????
?????????? ???????. ??????? ??? N = 104 ???????????? ?????????? ?????? ? ??????????? ??????????
? ??????????? ????????????? ??????? ???????
?????? ??????? ??????? ?????.
? ??????? 3 ???????????? ?????????? ????????? ??
???????????? ??????? ?????? [6], ? ??????? ? ????????
?????? ??????????????? ???????? ????????, ??? ????????? ????? ??????????; ? ??????? 4 ? ?????????
????????????? ?????? ? ?????? ???????? ?????????
?????????. ???????, ??? ??????, ??????????? ???????? ??????????? ????????? ??????? ?? ??????? (??
????, ??????????? ????????-??????? ?????? ? ??????
????, ??? ???????????), ??????????? ? ?????????????
?????????? H. ? ????? ???????? ???????? H ??????? ??
??????????? ??????? ? ? ???????, ????? ?????????
??????? ???????????? ????????? ????????, ????????????? ?????????? ??????? ????????.
??????? 3
????????? ??????? ????????? ???????? ? ???????????? [6]
X
RMS
*
R = 1, t = 0.19
N = 9*
N=7
0.0029
0.0011
R = 10, t = 0.80
N = 9*
N=7
0.0086
0.0072
R = 100, t = 0.66
N = 9*
N = 104
0.1175
0.0005
? ???????, ??????????? ???????????? ??????? ????????? ? ?????????????? ????????? ???????? [6].
??????? 4
????????? ??????? ????????? ???????? ? ???????????? [6]
X
RMS
*
N = 21**
0.0011
R = 10, t = 0.5
N = 12
0.0008
N = 21*
0.0005
N = 16
0.0002
R = 100, t = 0.32
N = 21**
0.0119
N = 40
0.0107
? ???????, ??????????? ??????? ???????? ????????? ????????? ? ?????????????? ???????????? ?? [6];
? ???????, ??????????? ??????? ???????? ????????? ????????? ? ?????????????? ???????? ?? [6].
**
?? ???????? 2 ? 3 ????????????????? ????????
??????? ????? ??? ?????? ?????????? R = 10 ?
R = 100. ??????????? ????????????? ?????? ???????,
??? ???????????? ????? ????????????? ????????
??????????? ? ?????????????? ?????, ??????????
? ??????????, ??????? ????? ??????????? ??? ??????? ???????? ?????.
???. 2. ???????? ??????
??????? ????? ??? R = 10
???. 3. ???????? ??????
??????? ????? ??? R = 100
186
?ыыбмзк?л?х ?н?с?кзо ??р?бз? м?????зк?л?? лбнклбхлсрзбх...
????????????????? ??????
1. Shen J. and Wang L. Some recent advances on spectral methods for unbounded domains // Commun. Comput.
Phys. ? 2009. ? Vol. 5, No. 2?4,
2. Boyd J.P. Chebyshev and Fourier spectral methods. ?
2nd ed. ? Dover Publication Inc., Mineola. ? NY, 2001.
3. ??????????? ?.?., ?????????? ?.?., ????????? ?.?. ???????? ? ????????????? ????????? ??????????????? ????. ? ?., 1967.
4. Shen J. Efficient spectral-Galerkin methods II. Direct solvers of second and fourth order equations by using
Chebyshev polinomials // SIAM J. Sci. Comput. ? 1995. ?
Vol. 16, No. 1.
5. ???????? ?.?. ?????????????? ?????? ???????. ?
?., 1949.
6. ??????? ?. ????????? ?????? ?? ?????? ??????
?????????. ? ?., 1988.
187
1/--страниц
Пожаловаться на содержимое документа