1
$\begingroup$

I want to solve Lid driven cavity flow with the help of projection method i want to make own fortran program to solve two dimensional unsteady flow.

I propose to solve the equations in the following way: apply a second order Adam Bashforth method for convection terms and Crank nicloson scheme for viscous term by neglecting pressure. I also want to calculate intermediate velocity first. i need help to develop a code for this in fortran 90.

$\endgroup$
3
  • 1
    $\begingroup$ You may want to look at this code, accompanied by this paper. $\endgroup$ – Paul Dec 19 '13 at 14:49
  • $\begingroup$ hi your paper is good u sent me but i want to work on some specific instructions which are given to me.need some help in programming related to fortran $\endgroup$ – sajjad ur Dec 19 '13 at 16:01
  • 2
    $\begingroup$ An "anonymous user", possibly the OP not logged in, proposed the following edit (which I rejected as a radical change): " I want to work in the following way 1.apply a second order Adam Bashforth method for convection terms and Crank nicloson scheme for viscous term by neglecting pressure. and want to calculate intermediate velocity first. i need help in fortran 90 code". I'm of the opinion it doesn't suffice to reopen the Q, even if the restatement is that of the OP. $\endgroup$ – hardmath Dec 26 '13 at 13:08
0
$\begingroup$

In your case, since you want to code yourself, I would start with a paper that comprehensively describes discretization.

For your needs D. Kim, S. Choi - A Second-Order Time-Accurate Finite Volume Method for Unsteady Incompressible Flow on Hybrid Unstructured Grids, JCP, 162 (2000) would probably be the right thing.

There you will find time-integration description, interpolations are shown with nice illustrations, a few words on matrix solver, and finally lid driven cavity flow is solved, so you can compare your results with theirs and verify your code.

$\endgroup$
5
  • $\begingroup$ thanx for ur concern do u have any fortran code for this type of problem actually i want to work on specific way $\endgroup$ – sajjad ur Dec 19 '13 at 15:56
  • $\begingroup$ In that case the alternative may be to look into a code of Bendiks Jan Boersma at TU Delft for simple temporal shear layer flow. It is a finite difference code, with staggered grid and using a projection method of Van Kan. You can modify it for lid driven flow. PDF document, describing the algorithm and implementation is included. It's here: ahd.tudelft.nl/~bendiks/demo/intro.html $\endgroup$ – Johntra Volta Dec 19 '13 at 21:38
  • $\begingroup$ hi thanx a lot for ur concern i want to specify my work with you hopefully u provide a help me $\endgroup$ – sajjad ur Dec 20 '13 at 5:55
  • $\begingroup$ I want to work in the following way 1.apply a second order Adam Bashforth method for convection terms and Crank nicloson scheme for viscous term by neglecting pressure. and want to calculate intermediate velocity first. i need help in fortran 90 code $\endgroup$ – sajjad ur Dec 20 '13 at 8:56
  • $\begingroup$ Follow the link given above (Boersma DNS code for shear layers). The code uses CN-AB combination in a Van Kan projection algorithm. Just look into the given fortran code and pdf with explanation. $\endgroup$ – Johntra Volta Dec 20 '13 at 11:37

Not the answer you're looking for? Browse other questions tagged or ask your own question.