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

Strange Hydrostatic Boundary Backflow with interFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 15, 2016, 21:32
Default Strange Hydrostatic Boundary Backflow with interFoam
  #1
KTG
Senior Member
 
Abe
Join Date: May 2016
Posts: 115
Rep Power: 9
KTG is on a distinguished road
Hello Foamers,

I am having a strange issue when using interFoam with a hydrostatic boundary condition. The general setting is that I have two plates with a narrow gap between them. A thick, laminar liquid is flowing into the gap under a constant head, hydrostatic boundary condition. So I am using the interFoam solver. The liquid that is displacing the air is orders of magnitude denser, so I just set the boundaries to zero pressure where the air is exiting the gap.

flow.png shows the simulation after 40 seconds.

backflow.png is the same, but with forward streamlines. You can clearly see that the section that is assigned the uniformDensityHydrostaticPressure condition has split itself between inflow and outflow. I can't figure this out, because the pressure there is greater than zero.

Here are the boundary conditions for pgh (the "freeinlet" patch is located directly above the uniformDensityHydrostaticPressure patch):

Code:
boundaryField
{
    inlet
  
     {
        type            uniformDensityHydrostaticPressure;
        rho             1263.9;
        pRefValue       0;
        pRefPoint       (0 0 .16);
        value           uniform 0; // optional initial value
    }

    freeinlet
    {
        type            fixedValue;
        value            uniform 0;

    }
    

    outlet
    {
        type            fixedValue;
        value            uniform 0;

    }

    wedgewalls
    {
        type            fixedFluxPressure;    
    }    

    lid
    {
        type            fixedValue;
        value            uniform 0;

    }

}
for U:

Code:
boundaryField
{

    inlet
    {
        type        zeroGradient;
    }    

    freeinlet
    {
        type        zeroGradient;
    }    

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

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

    lid
    {
        type            slip;
    }


}
and alpha.water:

Code:
 boundaryField
{
    inlet
    {
        type            inletOutlet;
        inletValue      uniform 1;
        value           uniform 0;
    }
    freeinlet
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
    wedgewalls
    {
        type            zeroGradient;
    }
    lid
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
    defaultFaces
    {
        type            empty;
    }
}
Any ideas would be greatly appreciated! I have been messing with the geometry of this to try to make it a little more stable, but I think the real issue is in my boundary conditions. I can't figure out why the zero pressure boundaries seem to have zero flux, but the hydrostatic one has some flow that seems to defy the pressure gradient.

Thanks!
Attached Images
File Type: png flow.png (16.9 KB, 29 views)
File Type: png backflow.png (25.6 KB, 24 views)
KTG is offline   Reply With Quote

Reply

Tags
boundary conditions, hydrostatic pressure, interfoam

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Centrifugal fan j0hnny CFX 13 October 1, 2019 13:55
Domain Imbalance HMR CFX 5 October 10, 2016 05:57
Wrong flow in ratating domain problem Sanyo CFX 17 August 15, 2015 06:20
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 07:00
[Commercial meshers] Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Meshing & Mesh Conversion 2 July 15, 2005 04:15


All times are GMT -4. The time now is 13:25.