at_aq_galerkin_RRSS.f90

Path: src/at_aq_galerkin_RRSS.f90
Last Update: Sat Mar 04 10:06:24 +0900 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 用に変更

This file provides following module.

Required files

[Validate]