CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

RasInterFoam STRANGE RESULTS AT BOUNDARY

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 22, 2006, 20:32
Default Hello friends I have been u
  #1
Senior Member
 
kumar
Join Date: Mar 2009
Posts: 112
Rep Power: 8
kumar2 is on a distinguished road
Hello friends

I have been using rasInterFoam to study wave making of hydrofoils in a channel. the problem i am experiancing is that although i have defined an outlet the fluid seems to be flowing back into the computational domain . let me explain this in detail. figure 1 shows the boundary conditions i have imposed


The boundary conditions for the atmosphere are exactly same as that of damBreak case in rasInterFoam tutorial. The boundaries for hydrofoil is simple ( basically a noslip wall)






figure 2, 3 , 4 & 5 shows the results after only 1000 time steps . As you can see there seems to be a serious error at the outlet. The velocity seems to have significantly slowed down and the flow seems to change direction and flow back into the domain ( see for example figure 5). i wonder what is happening to the mass conservation equations?

Could anyone please comment on these results

Thanks in advance
kumar
kumar2 is offline   Reply With Quote

Old   March 22, 2006, 21:05
Default Can't see the images - could y
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Can't see the images - could you please try and upload them again?

Thanks,

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   March 22, 2006, 21:12
Default Sorry friends my images are
  #3
Senior Member
 
kumar
Join Date: Mar 2009
Posts: 112
Rep Power: 8
kumar2 is on a distinguished road
Sorry friends

my images are too big . so could not upload. but let me explain with 2 figures. figure 1 shows gamma. you can also notice the inlet . there are 2 patches for the inlet so that gamma is 0 in the top and 1 in the bottom. U is specified and equals 0.8 in top & bottom patches . pd has zero gradient. there are also small values for k & epsilon .



The top boundary is atmosphere and the boundary conditions are same as in the tutorial of rasinterFoam.

The outlet is specified as two outlet patches with zero gradient for most quantities ( including U ). But gamma is 1 for the bottom patch of outlet and 0 for top patch of outlet. pd is set to zero as shown below ( in file /0/pd)
patch atmosphere
{
type totalPressure;
p0 uniform 0;
value uniform 0;
}

the hydrofoil is a noslip boundary and bottom of channel is symmetryPlane.



figure 2 shows the velocity field. you can see that the flow is decelerating and turning back into the computational domain . something terribly wrong is happening at the outlet

all comments welcome

thanks in advance
kumar
kumar2 is offline   Reply With Quote

Old   March 22, 2006, 21:16
Default Dear Hrv I would love to he
  #4
Senior Member
 
kumar
Join Date: Mar 2009
Posts: 112
Rep Power: 8
kumar2 is on a distinguished road
Dear Hrv

I would love to hear your comments . i have just uploaded the images

regards

kumar
kumar2 is offline   Reply With Quote

Old   March 22, 2006, 21:42
Default You need to change the outlet
  #5
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
You need to change the outlet boundary condition. Definitely not totalPressure on pd - use either fixedValue zero or zeroGradient (my bet is on zeroGradient, but this kind of thing is worth trying out). On U, use zeroGradient and on gamma do inletOutlet, using some sensible distribution for the refValue.

Your problem is definitely caused by the pressure b.c.

Good luck,

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   March 22, 2006, 22:42
Default Dear Hrv thanks a lot for y
  #6
Senior Member
 
kumar
Join Date: Mar 2009
Posts: 112
Rep Power: 8
kumar2 is on a distinguished road
Dear Hrv

thanks a lot for your comments. now that i know there is a problem with the pressure.b.c i am going to try out your suggestions. i will let you know once i fix the problem.

thanks a lot once again

regards

kumar
kumar2 is offline   Reply With Quote

Old   May 23, 2006, 18:07
Default Hello Dear Friends The figu
  #7
Senior Member
 
kumar
Join Date: Mar 2009
Posts: 112
Rep Power: 8
kumar2 is on a distinguished road
Hello Dear Friends

The figures show waves generated by a hydrofoil with rasInterFoam + nuTilda turb.model. The Boundary Conditions are like in the prevoius posts , as recommended by HRV. i got good results with k epsilon model . But k epsilon model has problems in the low Re range ( the Cp distribution has problems which in turn is reflected in the free surface deflections) . so i switched over to spallart allmaras model and my results are blowing up ! The free surface profile , velocity magnitude and nuTilda ( just before the solution blows up) are shown in the 3 figures. What seems to be happening is that the nuTilda gets very large near the atmosphere and the solution blows up. Please also note the patches of zero velocity near the atmosphere.





Note the hydrofoil is at the bottom left corner.

The top wall is modeled as the atmosphere with U as ( 0 0 0) - similar to the damBreak case . or should this be (0.8 0 0) ? where 0.8 is the inlet velocity
///////////////////////////////////////////
The important BC's for nuTilda are

InternalField = 5.7*10^-6 ;

1. atmosphere
type inletOutlet;
inletValue uniform 5.7e-06;
value uniform 5.7e-06;

2. Both the inlets ( one for air and one for water have fixedValue nuTildas of 5.7e-06 )

//////////////////////////////////////////////
The nuTilda value that i have specified are different from those given in the rasInterFoam tutorial . I found that with those ( the tutorial case )there was severe fluctuations of the free surface profiles - because of concentrations of nuTilda near the free surface. Also the nuTildas i have specified above resemble more closely the 'k' boundary condition.


All comments are welcome.

Thanks a lot

Kumar
kumar2 is offline   Reply With Quote

Old   March 9, 2007, 17:22
Default Hi, Dear all, Can anybody tel
  #8
New Member
 
JZ
Join Date: Mar 2009
Location: PBC, Florida, USA
Posts: 14
Rep Power: 8
jack2000 is on a distinguished road
Hi, Dear all,
Can anybody tell me if it is possible to use pressure inlet boundary ( without giving velocities) in interFoam solver. I tested this kind of boundary in simpleFoam, it worked. When I tried to use it in interFoam Solver, if gamma=0, it worked; but when gamma=1, it would broken-up.

What I intend to use this kind of boundary is to simulate constant head of water head passing through a hump in a open channel, and study the flow passing capacities. Thus it seems I need gamma=1?

Thanks in advance!

Jack
jack2000 is offline   Reply With Quote

Old   March 24, 2008, 19:38
Default I am trying to solve a back st
  #9
New Member
 
Onur Dundar
Join Date: Mar 2009
Location: Davis, CA, US
Posts: 7
Rep Power: 8
onurdundar is on a distinguished road
I am trying to solve a back step in rasinterfoam with k epsilon turbulence model. I added boundary
conditions in FoamX but when I run the rasInterFoam
it gives error about boundary conditions. I run after making the given correction. The result is this
#0 Foam::error::printStack(Foam:stream&) in "/home/dundar/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"

I can not find a true boundary conditions.
fvSchemes and FvSolution is the same with damBreak
Thanks for your concern

Onur


The boundary conditions are
inleta atmosphere
inletw inlet
outlet atmospehre
lowerwall wall

############################################
U
boundaryField
{
inleta
{
type pressureInletOutletVelocity;
phi phi;
value uniform (0.4289 0 0);
}

inletw
{
type fixedValue;
value uniform (0.4289 0 0);
}

outlet
{
type pressureInletOutletVelocity;
phi phi;
value uniform (0 0 0);
}

lowerWall
{
type fixedValue;
value uniform (0 0 0);
}

atmosphere
{
type pressureInletOutletVelocity;
phi phi;
value nonuniform List<vector>
#############################################
pd
boundaryField
{
inleta
{
type totalPressure;
p0 uniform 0;
U U;
phi phi;
rho none;
psi none;
gamma 1;
value uniform 0;
}

inletw
{
type zeroGradient;
}

outlet
{
type totalPressure;
p0 uniform 0;
U U;
phi phi;
rho none;
psi none;
gamma 1;
value uniform 0;
}

lowerWall
{
type zeroGradient;
}

atmosphere
{
type totalPressure;
p0 uniform 0;
U U;
phi phi;
rho none;
psi none;
gamma 1;
value uniform 0;
}
##############################################
k
boundaryField
{
inleta
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

inletw
{
type fixedValue;
value uniform 0.00184;
}

outlet
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

lowerWall
{
type zeroGradient;
}

atmosphere
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
##############################################
gamma
boundaryField
{
inleta
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

inletw
{
type fixedValue;
value uniform 1;
}

outlet
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

lowerWall
{
type zeroGradient;
}

atmosphere
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

################################################
R
boundaryField
{
inleta
{
type inletOutlet;
inletValue uniform (0 0 0 0 0 0 0 0 0);
value uniform (0 0 0 0 0 0 0 0 0);
}

inletw
{
type fixedValue;
value uniform (0 0 0 0 0 0 0 0 0);
}

outlet
{
type inletOutlet;
inletValue uniform (0 0 0 0 0 0 0 0 0);
value uniform (0 0 0 0 0 0 0 0 0);
}

lowerWall
{
type zeroGradient;
}

atmosphere
{
type inletOutlet;
inletValue uniform (0 0 0 0 0 0 0 0 0);
value uniform (0 0 0 0 0 0 0 0 0);
}
#################################################
Epsilon
boundaryField
{
inleta
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

inletw
{
type fixedValue;
value uniform 0;
}

outlet
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

lowerWall
{
type zeroGradient;
}

atmosphere
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
#################################################
nutilda
boundaryField
{
inleta
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

inletw
{
type fixedValue;
value uniform 0;
}

outlet
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

lowerWall
{
type fixedValue;
value uniform 0;
}

atmosphere
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

fvSchemes and FvSolution is the same with damBreak
Thanks for your concern

Onur
onurdundar 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
Newton Raphson - strange results Bren Main CFD Forum 5 February 11, 2009 11:47
Strange Block Mesh results ivan_cozza OpenFOAM Native Meshers: blockMesh 1 December 4, 2008 09:24
Slip boundary condition giving strange results hemph OpenFOAM Running, Solving & CFD 3 January 18, 2007 06:15
UDS transport: strange results Ale FLUENT 1 November 6, 2003 07:28
Strange results ! Dimitri FLUENT 1 August 8, 2002 11:07


All times are GMT -4. The time now is 15:39.