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

the internalField entry of U file is not valid in pisoFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 24, 2013, 05:32
Default the internalField entry of U file is not valid in pisoFoam
  #1
Member
 
Yao Lu
Join Date: May 2013
Posts: 32
Rep Power: 5
shuoxue is on a distinguished road
i was solving a 3-D cylinder flow case with pisoFoam.
the internalField entry of 0/U file was set to uniform (1 0 0).
however, the preview of 0s showed that the internalField entry failed to work. the U value of internalField remained to (0 0 0), see the figure below.
QQ??20131124164122.jpg

boundary
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

4
(
    inlet
    {
        type            patch;
        nFaces          1920;
        startFace       1747920;
    }
    outlet
    {
        type            patch;
        nFaces          1920;
        startFace       1746000;
    }
    cylinder
    {
        type            wall;
        nFaces          3840;
        startFace       1742160;
    }
    fb
    {
        type            wall;
        nFaces          36000;
        startFace       1706160;
    }
)

// ************************************************************************* //
0/U
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (1 0 0);

boundaryField
{
    inlet
    {
        type            freestream;
	freestreamValue $internalField;
    }
    outlet
    {
        type            freestream;
	freestreamValue $internalField;
    }
    cylinder
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    fb
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
}

// ************************************************************************* //
is there anything wrong with my code?
thanks in advance.
shuoxue is offline   Reply With Quote

Old   November 24, 2013, 06:23
Default
  #2
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,554
Blog Entries: 6
Rep Power: 27
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi,

did you check if your boundaries are correct.
It seems that inlet and outlet are not the patches you wanted to have.
Tobi is offline   Reply With Quote

Old   November 24, 2013, 07:56
Default
  #3
Member
 
Yao Lu
Join Date: May 2013
Posts: 32
Rep Power: 5
shuoxue is on a distinguished road
Quote:
Originally Posted by Tobi View Post
Hi,

did you check if your boundaries are correct.
It seems that inlet and outlet are not the patches you wanted to have.
Hi Tobias!

My computational domain is a circular cylinder. Inlet patch is on the left hand, and outlet patch is on the other side. Is there anything wrong about the boundary condition?

The similar geometry was used in a 2-D case with icoFoam. And internalField entry of U file worked.

shuoxue
shuoxue is offline   Reply With Quote

Old   November 24, 2013, 08:05
Default
  #4
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 1,164
Rep Power: 20
akidess will become famous soon enough
Do you have another file in 0/ which has the "object U" entry? E.g. a file U.org?

Note that since you are looking at the 0-folder, this has absolutely nothing to do with icoFoam or pisoFoam. The solver haven't touched your files yet.
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
*Join the OpenFOAM stackexchange Q&A site: http://area51.stackexchange.com/prop...oHPxcPqde7HtA2
akidess is offline   Reply With Quote

Old   November 24, 2013, 22:11
Default
  #5
Member
 
Yao Lu
Join Date: May 2013
Posts: 32
Rep Power: 5
shuoxue is on a distinguished road
Quote:
Originally Posted by akidess View Post
Do you have another file in 0/ which has the "object U" entry? E.g. a file U.org?

Note that since you are looking at the 0-folder, this has absolutely nothing to do with icoFoam or pisoFoam. The solver haven't touched your files yet.
Hi, I might found the primary issue. The fb patch is a wall type BC and its value of U is set to uniform (0 0 0). There is a transition layer between boundary and internalField. The real internalField part is surrounded by the layer, and it will be occurred when clipped.

U_Clip_0s.jpg
shuoxue is offline   Reply With Quote

Old   November 25, 2013, 05:05
Default
  #6
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 1,164
Rep Power: 20
akidess will become famous soon enough
The transition layer is probably only there because you chose to view point-interpolated field values in paraview.
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
*Join the OpenFOAM stackexchange Q&A site: http://area51.stackexchange.com/prop...oHPxcPqde7HtA2
akidess 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
OpenFOAM Installation for navalFoam sachinlb OpenFOAM Installation 21 June 23, 2014 08:07
Trouble compiling utilities using source-built OpenFOAM Artur OpenFOAM Programming & Development 14 October 29, 2013 11:59
swak4Foam-groovyBC build problem zxj160 OpenFOAM 18 July 30, 2013 13:14
build problem swak4Foam OF 2.2.0 mcathela OpenFOAM Installation 14 April 23, 2013 13:59
mesh airfoil NACA0012 anand_30 OpenFOAM Meshing & Mesh Conversion 12 December 12, 2011 05:16


All times are GMT -4. The time now is 09:22.