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)
(Example code for solving Smith-Hutton test)
Line 33: Line 33:
''Developer''
''Developer''
 +
 +
*[[ Example code for solving Smith-Hutton test - Fortran 90]]
== main.f90 ==
== main.f90 ==

Revision as of 01:21, 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

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

My wiki