CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

airfoil boundary condition

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   February 4, 2007, 08:46
Default airfoil boundary condition
  #1
Bounecer
Guest
 
Posts: n/a
I have written a Fortran program to simulate flow around airfoil:

My problem is: Determination of separation on naca airfoil at low Reynolds numbers (1e3 to 1e5)

Flow caracteristics: * Steady, 2D * Incompressible, viscous and Laminar

Mathematical model: * Continuity equation * 2 momentum equations

Discretization method: * Structured C-grid generated by "mesh pilot" software * the farfield boundary is 15xcord * transformation of model from physical (x, y) to computational domain (ksi, eta) * Finte volumes discretization * Staggered grid for velocity and pressure * Simpler Algorithm * Second order Upwind scheem for convective terms and central differences for diffusive terms

Solution of Algebraic equations: * Sweeping * Thomas Algrithm (TDMA)

My program diverges after about 5000 itterations althought i've used low relaxation parameters: relaxation for velocity 4e-1 and for pressure 1e-1.

Is that because i've not well imposed boundary conditions or because of other problems?

  Reply With Quote

Old   February 5, 2007, 10:08
Default Re: airfoil boundary condition
  #2
ag
Guest
 
Posts: n/a
At those Reynolds numbers, the flow is going to separate very soon after the maximum thickness. The wake flow will be unsteady. I don't have a clue why your solver is diverging - it could be anything from a boundary condition to a coding error to the phase of the moon. But I would definitely consider using an unsteady solver for the problem you are seeking to simulate.
  Reply With Quote

Old   February 5, 2007, 23:11
Default Re: airfoil boundary condition
  #3
Ahmed
Guest
 
Posts: n/a
Discretization method: * Structured C-grid generated by "mesh pilot" software * the farfield boundary is 15xcord * transformation of model from physical (x, y) to computational domain (ksi, eta) * Finte volumes discretization *

The book by Blazek is an excellent reference about the Finite volume method, thereis no need to transform from physical to computational domains, the finite volume method works directly with the physical domain, I guess you have mixed up different concepts in your code, check the book, it has a CD with fortran code hope that will help
  Reply With Quote

Old   February 7, 2007, 08:30
Default Re: airfoil boundary condition
  #4
Bounecer
Guest
 
Posts: n/a
I have set the following boundary conditions: on airfoil (low side of comput domain): u=v=0 farfield (upper side of comput domain): u=1, v=0 inflow (left side):du/dksi=dv/dksi=0 outflow (right side): du/dksi=dv/dksi=0 u(ksi, eta) and v(ksi,eta) the compnents of velocity.

Is that the correct way to set boundary conditions?
  Reply With Quote

Old   February 7, 2007, 09:55
Default Re: airfoil boundary condition
  #5
ag
Guest
 
Posts: n/a
What are you doing with pressure on the boundary?
  Reply With Quote

Old   February 8, 2007, 11:06
Default Re: airfoil boundary condition
  #6
bounecer
Guest
 
Posts: n/a
I'm using SIMPLER Algorithm which allow to find pressure from velocity
  Reply With Quote

Old   June 28, 2010, 09:24
Default
  #7
Senior Member
 
morteza08's Avatar
 
Morteza
Join Date: May 2010
Location: Iran,Islamic Republic of
Posts: 159
Rep Power: 7
morteza08 is on a distinguished road
Quote:
Originally Posted by Bounecer
;47729
I have written a Fortran program to simulate flow around airfoil:

My problem is: Determination of separation on naca airfoil at low Reynolds numbers (1e3 to 1e5)

Flow caracteristics: * Steady, 2D * Incompressible, viscous and Laminar

Mathematical model: * Continuity equation * 2 momentum equations

Discretization method: * Structured C-grid generated by "mesh pilot" software * the farfield boundary is 15xcord * transformation of model from physical (x, y) to computational domain (ksi, eta) * Finte volumes discretization * Staggered grid for velocity and pressure * Simpler Algorithm * Second order Upwind scheem for convective terms and central differences for diffusive terms

Solution of Algebraic equations: * Sweeping * Thomas Algrithm (TDMA)

My program diverges after about 5000 itterations althought i've used low relaxation parameters: relaxation for velocity 4e-1 and for pressure 1e-1.

Is that because i've not well imposed boundary conditions or because of other problems?

hi i am an aerospace engineer....i have to simulate the flow over an airfoil, i have to write a code and i have not enough time to do that...
can i borrow your code?
i prefer spending my time on a pre-written code and try to modify it for my purpose..

thanks
morteza08 is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Wind turbine simulation Saturn CFX 34 October 16, 2014 05:27
problem with boundary condition??? smn CFX 5 November 24, 2009 07:37
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05
how to set up a wall boundary condition according to calculated wall shear stress? gameoverli OpenFOAM Pre-Processing 1 May 21, 2009 08:28
CFX Solver : Sudden crash Hervé CFX 2 June 16, 2008 06:40


All times are GMT -4. The time now is 16:03.