# Sample code for solving Smith-Hutton test - Fortran 90

(Difference between revisions)
 Revision as of 20:28, 20 September 2005 (view source)Michail (Talk | contribs)← Older edit Latest revision as of 18:02, 20 April 2012 (view source)Michail (Talk | contribs) (5 intermediate revisions not shown) Line 5: Line 5: Sample program for solving Smith-Hutton Test using different schemes of covective terms approximation Sample program for solving Smith-Hutton Test using different schemes of covective terms approximation - '''Copyright (C) 2005  Michail Kirichkov ,Lithuania, Marijampole''' + '''Copyright (C) 2005  Michail Kiričkov ,Lithuania, Marijampole''' This program is free software; you can redistribute it and/or This program is free software; you can redistribute it and/or Line 21: Line 21: Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA. Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA. + -------------------------------------------------------------------- + '''Short description''' - *[[Bound_1.f90 - Set boundary conditions]] + In this program general transport equation for scalar transport is solved. - *[[Check_koeff.f90 - Performing checking of koefficients]] + Because it was primarily developing program for learning, there was implemented an output for all variables in the txt format as well as for TECPLOT format. - *[[Coeff_11.f90 - calculate the coefficients]] + Structured orthogonal uniform grid is used (now). All variables are saved in array '''F(1:nx,1:ny,10)''' which mean that 2D array for 10 variables. - *[[Contr_gem1.f90]] + - *[[COnvrg_sub.f90 - calculationg ]] + Variables '''F(:,:,1)''' and '''F(:,:,2)''' represented U and V velocity components, and  '''F(:,:,5)'''represented scalar. - *[[Geom_1.f90]] + - *[[Grid_sub.f90]] + Arrays '''Xp''' and '''Yp''' represents the coordinates of the CV centres, and '''X''' and '''Y''' - coordinates of grids nodes - *[[HLPA.f90]] + - *[[ini.f90]] + ---------------------------------------------------------------------- - *[[input.f90]] + - *[[Main.f90 - Main modul]] + - *[[Out_array_1.f90]] + *[[Bound_1.f90     - Set boundary conditions]] - *[[Output_aa.f90]] + *[[Check_koeff.f90 - Performing checking of koefficients (output for visual control) ]] - *[[TDMA_11.f90 - Solution of system of linear equatrions by Thomas method]] + *[[Coeff_11.f90     - calculate the coefficients]] - *[[icomm_1.f90 - include module]] + *[[Contr_gem1.f90   - control and output of geometry]] + *[[COnvrg_sub.f90   - calculationg of covergence kriteria]] + *[[Geom_1.f90       - calculation of geometry properties]] + *[[Grid_sub.f90     - definition of the grid]] + *[[HLPA.f90         - HLPA scheme implementation modul]] + *[[QUICK.f90        - QUICK scheme implementation modul]] + *[[ini.f90         - setting "0" to all arrays]] + *[[input.f90       - determination of intial conditions]] + *[[Main.f90         - Main modul]] + *[[Out_array_1.f90 - output of separate array in .txt file for visual control ]] + *[[Output_aa.f90   - output of all necessary information ]] + *[[TDMA_11.f90     - Solution of system of linear equatrions by Thomas method]] + *[[icomm_1.f90     - include module]]

## Latest revision as of 18:02, 20 April 2012

Dear friends

It's just a scrap. Later I'll correct it, although it's a complete working code

Sample program for solving Smith-Hutton Test using different schemes of covective terms approximation Copyright (C) 2005 Michail Kiričkov ,Lithuania, Marijampole

This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.

Short description

In this program general transport equation for scalar transport is solved. Because it was primarily developing program for learning, there was implemented an output for all variables in the txt format as well as for TECPLOT format. Structured orthogonal uniform grid is used (now). All variables are saved in array F(1:nx,1:ny,10) which mean that 2D array for 10 variables.

Variables F(:,:,1) and F(:,:,2) represented U and V velocity components, and F(:,:,5)represented scalar.

Arrays Xp and Yp represents the coordinates of the CV centres, and X and Y - coordinates of grids nodes