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

Unexpected behavior for interFoam with fixedGradient pressure BC

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 31, 2013, 05:47
Default Unexpected behavior for interFoam with fixedGradient pressure BC
  #1
Member
 
Join Date: May 2012
Posts: 55
Rep Power: 5
styleworker is on a distinguished road
Hello,

for validation of my actual case I'm simulating the interphase flow of the capillary rise tutorial with modified pressure boundary.

In my case, the inlet pressure BC should have a neg pressure gradient. The pressure gradient is higher than the pressure gradient produced by the curvature of the free surface, so that the fluid moves downwards and out of the inlet patch. Basically it does move downwards, but the flow flux is way too low (see case 4 in the appended picture).

I also tested a case with a positive pressure gradient (case 3). As expected the fluid moves upwards, but this case is way slower than the standard capillary rise tutorial (case 2). Actually it should move faster, because there is a additional pressure gradient beside the pressure gradient produced by the interphase.

The applied BC's for pressure are:

Code:
    inlet
    {
        type            fixedGradient;
        gradient        uniform -100;
    }

    atmosphere
    {
        type            fixedValue;
        value           uniform 0;
    }

    walls
    {
        type            fixedFluxPressure;
        adjoint         no;
    }
velocity BC's:
Code:
inlet
    {
        type            pressureInletOutletVelocity;
        value           uniform (0 0 0);
    }

    atmosphere
    {
        type            pressureInletOutletVelocity;
        value           uniform (0 0 0);
    }

    walls
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
alpha BC's:
Code:
    inlet
    {
        type            inletOutlet;
        value           uniform 1;
        inletValue      uniform 1;
    }

    atmosphere
    {
        type            zeroGradient;
    }

    walls
    {
        type           constantAlphaContactAngle;
        theta0         45;
        limit          gradient;
        value          uniform 0;
    }
Attached Images
File Type: jpg capillary_rise.jpg (32.3 KB, 37 views)
Attached Files
File Type: zip cap.zip (30.4 KB, 4 views)
styleworker is offline   Reply With Quote

Old   January 31, 2013, 11:34
Default
  #2
Member
 
Join Date: May 2012
Posts: 55
Rep Power: 5
styleworker is on a distinguished road
The solution is actually clear. I assumed before, that the reference pressure is the pressure of the ambient air, which does not make sense :-)

Last edited by styleworker; January 31, 2013 at 12:36.
styleworker is offline   Reply With Quote

Old   May 20, 2013, 22:56
Default
  #3
New Member
 
Nikhil
Join Date: Sep 2011
Posts: 9
Rep Power: 5
Nikhilcfd is on a distinguished road
Hello,

I faced a similar issue with capillary rise in interFOAM. In case 3, (where you applied positive pressure gradient) the fluid should move faster as you explained? Could you please explain how I can specify reference pressure and on what basis?

Whenever I applied positive or negative (or counter) pressure gradient (across the capillary tube/plate), I did not get required rise. When I set interface compression term (Calpha) to zero, the height rise was good in comparison with the analytical solution. For cases where I don't have pressure gradient across the capillary tube (i.e p=0 at both inlet and outlet), I got good results for both Calpha = 1 and 0 (interface was diffuse for c=0 but height rise was good). It appears that the Calpha term is imposing additional pressure on interface for cases with positive/negative pressure gradients across the capillary tubes. Please let me know your suggestions.

Thanks,

NP
Nikhilcfd 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
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
Unsteady pressure differential between inlet and outlet of the pipe for single phase joshi20h FLUENT 0 September 26, 2012 12:41
How to get the total pressure in the UDF? zgzhai Fluent UDF and Scheme Programming 2 September 12, 2012 06:15
Difference between pressure, absolute pressure and Total Pressure shaswat CFX 1 September 6, 2012 06:12
Setup/monitor points of pressure and force coefficients siw CFX 3 October 22, 2010 06:07


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