CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Wiki > Approximation Schemes for convective term - structured grids -...

Approximation Schemes for convective term - structured grids - Summary of Discretizations Schemes and examples

From CFD-Wiki

(Difference between revisions)
Jump to: navigation, search
(Example code for solving Smith-Hutton test)
(main.f90)
Line 35: Line 35:
*[[ Example code for solving Smith-Hutton test - Fortran 90]]
*[[ Example code for solving Smith-Hutton test - Fortran 90]]
-
 
-
== main.f90 ==
 
-
 
-
Program Main
 
-
 
-
include 'icomm_1.f90'
 
-
 
-
Dimension F_out(nx,ny)
 
-
Character  Filename*10
 
-
 
-
Call Ini
 
-
 
-
Call Grid_rectangular
 
-
Call input_all
 
-
Call Geom
 
-
Call Init_all
 
-
    Call Control_gem
 
-
 
-
Do 100 Niter=1,100
 
-
 
-
write(*,*) '------------------------',Niter,'------------------------',Niter
 
-
 
-
Call Coef_1(5)
 
-
Call Check_koeff
 
-
 
-
!  Alfa = 0.95   
 
-
!  Call Under_Relaxation(Alfa,5)
 
-
 
-
    Call TDMA_1(5)
 
-
 
 
-
Call Bound_condition
 
-
 
-
  Call Convergence_Criteria(5)
 
-
 
-
100 continue
 
-
 
-
 
-
!----------------------------------------------------------------
 
-
NImax = NXmaxp
 
-
NJmax = NYmaxp
 
-
 
-
F_out    = F(:,:,1)
 
-
 
-
!         1234567890
 
-
Filename  ='1_U_s.txt'
 
-
 
-
!    Call  Out_array(F_out,NImax,NJmax,Filename)
 
-
!-------------------------------------------------------------------
 
-
!----------------------------------------------------------------
 
-
NImax = NXmaxp
 
-
NJmax = NYmaxp
 
-
 
-
F_out    = F(:,:,2)
 
-
 
-
!         1234567890
 
-
Filename  ='1_V_s.txt'
 
-
 
-
!    Call  Out_array(F_out,NImax,NJmax,Filename)
 
-
!-------------------------------------------------------------------
 
-
!----------------------------------------------------------------
 
-
NImax = NXmaxp
 
-
NJmax = NYmaxp
 
-
 
-
F_out    = F(:,:,5)
 
-
 
-
!         1234567890
 
-
Filename  ='1_T_s.txt'
 
-
 
-
!    Call  Out_array(F_out,NImax,NJmax,Filename)
 
-
!-------------------------------------------------------------------
 
-
 
-
 
-
 
-
 
-
Call output_all
 
-
 
-
WRITE(*,*) 'PRIVET'
 
-
STOP
 
-
END
 

Revision as of 01:22, 19 September 2005

When we shall fill this page, I offer to make common identifications, because in different issues was used different notation.

Also we beg everybody to help me with original works. Later I shall write, what is necessary. If anyone have literature connected with convective schemes, please drop me a line.

We shall be very glad and grateful to hear any critical suggestion (please drop a few lines at Wiki Forum)

It is just a skeleton, but we hope that it will be developed into the good thing


Contents

Discretizations Schemes Estimation of order

Discretizations Schemes Estimation of error

Selection advice

Comparison of Discretizations Schemes

Numerical examples

Pure convection of a scalar step by a rotating velocity field (Smith-Hutton test)

R.M.Smith and A.G.Hutton (1982), "The numerical treatment of advection: A performance comparison of current methods", Numerical Heat Transfer, Vol. 5, p439.


Square Lid-driven cavity flow

Example code for solving Smith-Hutton test

Dear friends

It's jasy a scetch. Later I'll correct it, although it's a working code

Developer

My wiki