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

InterFoam error for two phase flow simulation in channel

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 7, 2014, 15:34
Default InterFoam error for two phase flow simulation in channel
  #1
Senior Member
 
Ali reza
Join Date: Mar 2014
Posts: 110
Rep Power: 12
1988 is on a distinguished road
hello
This is a two phase flow and the solver is interFoam in openFoam software.
As you can see here,I have made this geometry in gambit software and because of meshing problems I had to define some faces between volumes for example face number 1 and 2 and I used "specify continuum types"option in gambit to let the flow to pass.(there is flow in this geometry that it needs to pass these faces).
I checked the mesh and it was right but as soon as I want to see the geometry in paraview with paraFoam command I face this errors:
Quote:
boundary field walls not found in object U at time=0
and it is repeated for alpha-1 and p objects too.
whenever I type interfoam command it gives me this error
Quote:
cannot find patchField entry for walls
I have checked the mesh and also I solve it in fluent with no problem.
I don't know what is wrong!!! if some one can help me I would be so appreciate.
thanks
Attached Images
File Type: jpg 1 (2).jpg (53.8 KB, 21 views)
1988 is offline   Reply With Quote

Old   April 7, 2014, 20:56
Default
  #2
Member
 
Christian Butcher
Join Date: Jul 2013
Location: Japan
Posts: 85
Rep Power: 12
chrisb2244 is on a distinguished road
In your set of case files, you have probably a "0" directory, right?
Within that directory will be a number of files with names like "U", "p", "alpha.phasename" or "alpha1" depending on your OF version.

If you open these files with a text editor, you will read the boundary conditions at the bottom.

The error you're getting is complaining that in your mesh files, you have a patch face (or faces) called "walls", and that "walls" are not specified in your "U", "p" or "alpha1" files.

Probably you're going to want something like

Code:
boundaryField
{
    walls
    {
         type zeroGradient;
    }

    other_patch_faces_here

}
On a side note, if you don't have a "0" directory, these files are going to be held in whatever your first timestep is

Edit : Correction - having looked over the question again, it occurs to me that you don't want an actual wall where your "walls" patch is/are, if "walls" is the name you gave to the ones you added and highlighted in the picture (Apologies, I have no experience with Fluent or Gambit).
As such, you probably want some sort of inletOutlet BC - maybe someone else would know more specifically here.

Last edited by chrisb2244; April 7, 2014 at 20:59. Reason: More information about BC types
chrisb2244 is offline   Reply With Quote

Old   April 8, 2014, 04:10
Default
  #3
Senior Member
 
Ali reza
Join Date: Mar 2014
Posts: 110
Rep Power: 12
1988 is on a distinguished road
thanks for protect.
unfortunately I have defined these wall boundary conditions in U ,alpha 1 and so on but the error is still staying.
let me show you one of boundary conditions which I defined.
this is alpha1 boundary condition:
Code:
dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet-1
    {
        type            fixedValue;
    value        uniform 0.0;
    }
    inlet-2
    {
        type            fixedValue;
    value        uniform 1.0;
    }
     inlet-3
    {
        type            fixedValue;
    value        uniform 1.0;
    }
    outlet
    {
        type            zeroGradient;
    }

    Walls
    {
        type            zeroGradient;
    }


}

// ************************************************************************* //
1988 is offline   Reply With Quote

Old   April 8, 2014, 04:12
Default
  #4
Member
 
Christian Butcher
Join Date: Jul 2013
Location: Japan
Posts: 85
Rep Power: 12
chrisb2244 is on a distinguished road
Quote:
Originally Posted by 1988 View Post
thanks for protect.

Code:
    Walls
    {
        type            zeroGradient;
    }
This is good, but your Walls != walls OpenFOAM is case sensitive - if you change this to walls instead of Walls you should be fine
chrisb2244 is offline   Reply With Quote

Old   April 8, 2014, 14:09
Default
  #5
Senior Member
 
Ali reza
Join Date: Mar 2014
Posts: 110
Rep Power: 12
1988 is on a distinguished road
thanks a lot.
I have done what you say but I have just one error that it is about alpha1.org .
Quote:
ERROR: In /build/buildd/paraview-3.98.0/VTK/IO/Geometry/vtkOpenFOAMReader.cxx, line 6875
vtkOpenFOAMReaderPrivate (0x9855bb0): boundaryField walls not found in object alpha1.org at time = 0
and this is my alpha1 boundary condition:
Code:
dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    water-1
    {
        type            fixedValue;
    value        uniform 0.0;
    }
    water-2
    {
        type            fixedValue;
    value        uniform 1.0;
    oil
    {
        type            fixedValue;
    value        uniform 1.0;
    }

    outlet
    {
        type            zeroGradient;
    }

   walls
    {
        type            zeroGradient;
    }

}

// ************************************************************************* //
thank you
1988 is offline   Reply With Quote

Old   April 8, 2014, 21:06
Default
  #6
Member
 
Christian Butcher
Join Date: Jul 2013
Location: Japan
Posts: 85
Rep Power: 12
chrisb2244 is on a distinguished road
If you're getting the same error for alpha1.org, it's because your alpha1.org file doesn't have a "walls" boundary field. However, it was my understanding that the alpha1.org files provided in tutorial cases are usually used by some sort of 'Allrun' script to copy out an alpha1 file before starting solutions, so I'm surprised that your solver cares about an alpha1.org file.

In particular, the default interFoam does not read an alpha1.org file. So my guess is that this error is only occurring if you try to run foamToVTK or some similar conversion utility?

I think the solution is simply to remove the alpha1.org file from your 0 directory, if you do not need it.

I have a test case I'm using to try and troubleshoot bugs in a solver/library pair I'm writing, and I find it helpful to have a script saying something like

Code:
rm -rf 0*
cp -rf backup0 0
myInterFoam 2>&1 | tee log
Here you can see that my "backup0" directory holds a set of all the field files for time 0, so if you need alpha1.org for a backup, then you might have some luck moving it to a different directory (than 0)?

Best,
Christian
chrisb2244 is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
PEM fuel cell simulation pchoopanya Mesh Generation & Pre-Processing 1 February 29, 2016 23:56
two phase flow through a micro channel ajithjec FLUENT 0 September 10, 2013 13:39
Help! How to hold turbulence in a 2d channel flow RANS simulation 1/153 Main CFD Forum 7 July 27, 2012 15:39
Modeling the mixing of air and kerosene in a flow channel StefanG CFX 3 June 11, 2012 20:21
[ICEM] Flow channel meshing problems StefanG ANSYS Meshing & Geometry 19 May 15, 2012 06:44


All times are GMT -4. The time now is 05:48.