# Inviscid 2D Airfoil Case

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

 May 15, 2007, 12:37 I am running an inviscid case #1 Member   Doug Hunsaker Join Date: Mar 2009 Location: Logan, UT Posts: 63 Rep Power: 9 I am running an inviscid case for a 2D airfoil and am running into some snags. I'm using the potentialFoam solver and have tried to mimic the potentialFoam/cylinder case in the tutorials directory. However, my boundary conditions must be a little off because the potential flow field I get as the solution is not correct. I'm attaching a picture of the resultant flowfield. If anyone has any suggestions, I would greatly appreciate them. Thanks. -Doug inviscid.jpg

 May 17, 2007, 17:35 Hi Doug, Reading your threa #2 Senior Member   Vincent RIVOLA Join Date: Mar 2009 Location: France Posts: 281 Rep Power: 10 Hi Doug, Reading your thread I'm not really surprised you get velocity vectors perpendicular to the exit plane. Actualy, it comes straight forward from your boundary condition which is zeroGradient at the exit. The zeroGradient boundary condition trough a plane will always give you vectors perpendicular to this plane. Indeed if vectors were not perpendicular, you would end up with a difference in velocity between two points side to side the exit boundary plane. And so, you would have a gradient. If you really want to have vectors with 45°angle at the exit (which I'm not sure it's correct), you have to prescribe them as you do at the entrance. Hope that helped, Vincent

 May 17, 2007, 17:44 Vincent- Doesn't the zeroGr #3 Member   Doug Hunsaker Join Date: Mar 2009 Location: Logan, UT Posts: 63 Rep Power: 9 Vincent- Doesn't the zeroGradient velocity boundary condition constrain the velocity vector on one side of the boundary to match the velocity vector on the other side of the boundary? This shouldn't constrain the velocity vector to be perpendicular to the boundary. It should just constrain the difference in velocity vectors from one side of the boundary to the next to be zero. Also, specifying the velocity field at the exit doesn't allow for mass conservation across the domain. Therefore, I cannot specify the exit velocity profile. I realize the velocity profile at the exit won't be exactly uniform, but it should be very close to uniform because the flowfield is based off of potential flow. Right? Thanks for your help. -Doug

 May 17, 2007, 19:31 Hi, What about the lateral bo #4 Super Moderator   Takuya OSHIMA Join Date: Mar 2009 Location: Niigata City, Japan Posts: 518 Blog Entries: 1 Rep Power: 12 Hi, What about the lateral boundaries? If they are of fixed wall types the outflow naturally gets perpendicular... Takuya

 May 18, 2007, 11:16 Takuya- Do you mean the upp #5 Member   Doug Hunsaker Join Date: Mar 2009 Location: Logan, UT Posts: 63 Rep Power: 9 Takuya- Do you mean the upper and lower sides of the parabola? I'm using a C-grid, so there is a parabolic shape around the front of the airfoil. Yes, the entire length of that parabola is set to a constant velocity and zeroGradient pressure. Is that correct for inlet/outlet boundary conditions of this type? I can't imagine that this is causing the exit velocity to align perpendicular to the exit boundary. The parabola acts as a velocity inlet and is not set as a wall. Therefore, it should not cause the flow to be forced out perpendicular to the exit boundary. If by lateral boundaries you mean the walls in the z direction (this case is 2D), I have those set to empty as I believe they should be. -doug

 April 15, 2008, 03:51 Hi all, I am wondering we d #6 New Member   Richard Jones Join Date: Mar 2009 Location: Adelaide, South Australia, Australia Posts: 22 Rep Power: 9 Hi all, I am wondering we deal with the kutta-joukowski condition at the trailing edge in OF potential solutions? Any help would be greatly appreciated, Richard.

 September 9, 2010, 07:12 #7 New Member   Join Date: Jul 2010 Posts: 20 Rep Power: 8 Hello all, I'm having the same issue with outlet velocity vectors that seem to be forced in a perpendicular direction to the outlet patch face. I'm dealing with a stator blade which is part of a cascade. Thus, I use periodic boundaries with translational periodicity in the y-direction. Due to this periodicity, inlet and outlet boundaries are planar in y-direction with normal facing in x-direction. The flow enters the domain with a certain angle to the x direction (achieved with pressureDirectedInletOutletVelocity boundary condition for the velocity at the inlet). Due to the deflection of the flow caused by the blade, the flow should exit the domain with a certain angle, which differs form that at the inlet face (but is unknown a priori, thus it's not possible to give this direction as condition...). All simulations I've done so far show that the velocity vectors at the outlet are forced to be in positive x-direction, perpendicular to the face. This is not what one would physically expect and I'm very confused about that behaviour. Does anybody else have some experience with these kinds of problems? Do I have to use different outlet boundary conditions? Many thanks in advance!

 September 9, 2010, 08:49 #8 New Member     Darío Montes Join Date: Aug 2009 Location: Córdoba, Argentina Posts: 12 Rep Power: 9 Hi Doug! Did you fix your problem with the boundaries? Im having similar problems with a 0015 in the trailing edge with high velocities (60% and 70% over Uinf). The rest of the flowfield appears to be ok except the t.e. Thanks, regards!!

 October 13, 2010, 07:06 #9 New Member   Join Date: Jul 2010 Posts: 20 Rep Power: 8 Hi all, just to let you know: I fixed the problem of forced perpendicular velocity vectors at the outlet. The problem was something completely different: I had to recompile my gcc with another set of optimization flags. Apparently they have a significant influence on the outcome of the simulation results. I might be wrong on that, maybe it was not only the optimization but also the compiler version. Anyway, the problem has been resolved for me now. Maybe you want to have a look on your compiler and compiler flags if you face similar problems.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post doug OpenFOAM Running, Solving & CFD 1 September 9, 2010 17:24 R ARUN BALAJI FLUENT 0 April 15, 2008 23:33 S.T Main CFD Forum 6 December 10, 2007 09:07 Akbar Main CFD Forum 1 July 22, 2005 11:49 Akbar FLUENT 0 July 15, 2005 06:49

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