CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

trouble with `pressureDirectedInletVelocity'

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 6, 2009, 04:38
Question trouble with `pressureDirectedInletVelocity'
  #1
New Member
 
Yohey Ishizu
Join Date: May 2009
Location: Tokyo, Japan
Posts: 6
Rep Power: 7
yohey is on a distinguished road
Hi.

I'm using OpenFOAM-1.6.x on Ubuntu 8.04 and in trouble with `pressureDirectedInletVelocity' boundary.


I'm working on the attached case. This case contains direction-fixed boundary as
inletDirection uniform (1 1 0);
But, with this condition, the velocity in the mesh next to boundary is dobuled.
This is because of `operator=' in `pressureDirectedInletVelocityFvPatchVectorField.C '.
fvPatchField<vector>::operator=(inletDir_*(inletDi r_ & pvf));
So, I modified the `inletDirection' as (1/sqrt(2) 1/sqrt(2) 0).
inletDirection uniform (0.707107 0.707107 0);
This worked correctly.


(1) Next, I tried with more precise value.
inletDirection uniform (0.70710678 0.70710678 0);
But, calculation fails with this boundary.


(2) Actually this problem is derived from `pressureDirectedInletVelocityFvPatchVectorField.C '. So, I modified such sentence in this file as
fvPatchField<vector>::operator=(inletDir_*(inletDi r_ & pvf)/magSqr(inletDir_));
This modification also crashes calculation.


I can't understand these behavior of (1) and (2).

Would anyone help me to find the cause of these problems?


Sincerely.

Yohey
Attached Files
File Type: gz inletTest.tar.gz (2.0 KB, 45 views)
yohey is offline   Reply With Quote

Old   August 10, 2011, 12:42
Default
  #2
New Member
 
Ali Q Raeini
Join Date: Feb 2010
Posts: 21
Rep Power: 6
aliqasemi is on a distinguished road
Yohey[/QUOTE]
Quote:
Originally Posted by yohey View Post
Hi.

(1) Next, I tried with more precise value.
inletDirection uniform (0.70710678 0.70710678 0);
But, calculation fails with this boundary.

Yohey
This is because some of the faces in the patch are normal to the vector you specified for inletDirection, which causes a division by zero in the file pressureDirectedInletOutletVelocityFvPatchVectorFi eld.C. Changing the code or the value for inletDirection should fix the problem. I had the same problem using using OF.1.7.1, nearly two years later.
aliqasemi is offline   Reply With Quote

Old   October 4, 2013, 12:15
Smile Inlet Angle
  #3
New Member
 
Akash Sharma
Join Date: May 2013
Location: Paris
Posts: 15
Rep Power: 3
akashjangid is on a distinguished road
Hello

I am doing simulation on a rotor blade.
Does anyone know how to give the inlet angle like this (which depends on z and y)
dir_cos1 = 4.04
dir_cos2 = EXPR{57.2958*(acos(0.07045*(z/sqrt(z*z+y*y))))}
dir_cos3 = EXPR{57.2958*(acos(-0.07045*(y/sqrt(z*z+y*y))))}

Thanks !!
akashjangid is offline   Reply With Quote

Reply

Tags
boundary condition, inlet

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
[GAMBIT] Trouble meshing complex VOF geometry RPJones ANSYS Meshing & Geometry 2 February 14, 2011 19:54
Finned Pipe Trouble fabioxyz CFX 8 August 31, 2009 03:20
Trouble in Gambit Tom FLUENT 8 July 23, 2002 10:17
equivalencing trouble Steve Collie CFX 6 August 22, 2001 19:26
trouble with fluxsum Jordi CD-adapco 2 September 29, 2000 04:57


All times are GMT -4. The time now is 21:08.