CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

trouble with `pressureDirectedInletVelocity'

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

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

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?


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

Old   August 10, 2011, 11:42
New Member
Ali Q Raeini
Join Date: Feb 2010
Posts: 21
Rep Power: 9
aliqasemi is on a distinguished road
Originally Posted by yohey View Post

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

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, 11:15
Smile Inlet Angle
New Member
Akash Sharma
Join Date: May 2013
Location: Paris
Posts: 15
Rep Power: 6
akashjangid is on a distinguished road

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


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 02:20
Trouble in Gambit Tom FLUENT 8 July 23, 2002 09:17
equivalencing trouble Steve Collie CFX 6 August 22, 2001 18:26
trouble with fluxsum Jordi Siemens 2 September 29, 2000 03:57

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