| Path: | src/at_aq_galerkin_RRSS.f90 |
| Last Update: | Sat Mar 04 10:06:24 JST 2006 |
COPYRIGHT (c) 2006 Shin-ichi Takehiro. All rights reserved.
表題 at_aq_galerkin_RRSS1
チェビシェフ−ガラーキン法
非圧縮流体の流線関数・流れポテンシャル用
case RR: 両端粘着条件(両端で値と 1 階微分が 0)
f(:,i=0)=f(:,i=im)=f'(:,i=0)=f'(:,i=im)=0
[ f(:,x=xmax)=f(:,x=xmin)=f'(:,x=xmax)=0=f'(:,x=xmin)=0 ]
case SS: 両端自由すべり条件(両端で値と 2 階微分が 0)
f(:,i=0)=f(:,i=im)=f''(:,i=0)=f''(:,i=im)=0
[ f(:,x=xmax)=f(:,x=xmin)=f''(:,x=xmax)=0=f''(:,x=xmin)=0 ]
case RS: 片端粘着条件, 片端自由すべり条件
(両端で値が 0, 片側で 1 階微分, もう一方で 2 階微分が 0)
f(:,i=0)=f(:,i=im)=f'(:,i=0)=f''(:,i=im)=0
[ f(:,x=xmax)=f(:,x=xmin)=f''(:,x=xmax)=0=f'(:,x=xmin)=0 ]
case SR: 片端自由すべり条件, 片端粘着条件
(両端で値が 0, 片側で 2 階微分, もう一方で 1 階微分が 0)
f(:,i=0)=f(:,i=im)=f''(:,i=0)=f'(:,i=im)=0
[ f(:,x=xmax)=f(:,x=xmin)=f'(:,x=xmax)=0=f''(:,x=xmin)=0 ]
k 次のガラーキン基底が
\phi_k(x)=T_k(x) + C_{k-1}T_{k-1}(x) + C_{k-2}T_{k-2}(x)
+ C_{k-3}T_{k-3}(x) + C_{k-4}T_{k-4}(x)
の型式(Type1)
定式化については解説文書
「チェビシェフ関数展開を利用したガラーキン法」
(cheb_gal.pdf)を参照のこと.
履歴 2006/01/24 竹広真一 新規作成
2006/01/25 竹広真一 モジュール初期化フラッグ導入
2006/02/28 竹広真一 コメントを RDoc 用に変更