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

Interfoam boundary setup problem tank with dip tube in tank

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 14, 2022, 14:10
Default Interfoam boundary setup problem tank with dip tube in tank
  #1
New Member
 
Daniel Norton
Join Date: Feb 2015
Posts: 17
Rep Power: 11
Dano62 is on a distinguished road
I am trying to model flow through a vented dip tube into a tank, and I'm currently having difficulty setting stable boundary conditions. Then tank has liquid feeding through a dip tube which is vented back to the tank. The tank is maintained at a constant level through a level controller, and I've set prgh to correspond with this height. I'm using interfoam with local timestepping to first solve for an average free surface.

The liquid level stays appropriately stable, and the liquid flow appears to develop properly but the air velocity field explodes. I suspect that there's some sort of feedback loop developing between the outlet and air inlet causing issues.

I've attached images summarizing the setup and main issue for reference as included my U and p_rgh setup below.

The figure with contours is a simulation stopped early. The blue crop is velocities exceeding 10m/s of air, while the green crop shows the current liquid levels. The air velocity through the air inlet should be near negligible instead of consistently accelerating. It ranges from 10m/s to 80m/s at the time of this image.

The issue I'm having is in properly specifying the inletAir conditions, or possibly the outlet conditions. It seems no matter what I specify the system becomes unstable and the velocity of air through the inlet continuously increases. I've moved the air inlet away from the vent line, increased mesh refinement of the vent, and tried a variety of boundary conditions for the inlet and outlet but always seem to have stability issues.

The U and p_rgh conditions currently match those in the waterChannel tutorial.

Any help at all would be much appreciated. Thanks!

U:
Code:
boundaryField
{
    #includeEtc "caseDicts/setConstraintTypes"

    inletAcid
    {
        type                flowRateInletVelocity;
        volumetricFlowRate  constant 0.278931;
    }

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

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

    "walls.*"
    {
        type            noSlip;
    }
}
p_rgh:
Code:
boundaryField
{
    #includeEtc "caseDicts/setConstraintTypes"

    inletAcid
    {
        type            fixedFluxPressure;
        value           uniform -64948;
    }

    inletAir
    {
        type            totalPressure;
        p0              uniform -64948;
    }

    outlet
    {
        type            fixedFluxPressure;
        value           uniform -64948;
    }

    "walls.*"
    {
        type            fixedFluxPressure;
        value           uniform -64948;
    }
}
Attached Images
File Type: png setup.png (178.4 KB, 14 views)
File Type: jpg air_problem.jpg (35.9 KB, 13 views)
Dano62 is offline   Reply With Quote

Old   August 23, 2022, 12:50
Default
  #2
New Member
 
Daniel Norton
Join Date: Feb 2015
Posts: 17
Rep Power: 11
Dano62 is on a distinguished road
I found a working solution to this, so I'll post here in case anyone is looking for a reference later.

In the model, air was consistently accelerating from the air inlet through the combined air/liquid outlet. It's essentially a pressure-pressure boundary condition pairing with negligible line losses. The solution was to split the outlet region into a dedicated liquid outlet below the free surface and walls. The air inlet then acted as an atmospheric vent and the only air opening in the system.

I redefined the p_rgh pressure at the air inlet to a totalPressure of 0, and set the new liquid outlet P_rgh to a fixedValue of -64948 which corresponds to the static head of liquid surface beneath the origin.

The end result was a consistent free surface at the elevation I expected, and the air inlet acting as an atmospheric vent. The downside is that I've introduced an artificial wall into the model, but it's far enough away from my area of interest to not be a significant problem.
Dano62 is offline   Reply With Quote

Old   September 13, 2022, 06:07
Default
  #3
New Member
 
Eric Segalerba
Join Date: Dec 2021
Location: Italy
Posts: 24
Rep Power: 4
EricS is on a distinguished road
Hi Daniel, I'm facing a similar issue and I'd like to solve it in the same way. Could I ask you how were you able to split the outlet patch in two? Because it's a flat patch so I don't really know how to divide it in two.


Thank you in advance
EricS is offline   Reply With Quote

Old   September 13, 2022, 15:50
Default
  #4
New Member
 
Daniel Norton
Join Date: Feb 2015
Posts: 17
Rep Power: 11
Dano62 is on a distinguished road
I modified the geometry to have a dedicated outlet that was removed from my area of interest. I built the geometry in Salome, so I made the change there. You can also split a face in Salome, which may help you in modifying your geometry too if you just want to split an existing face. Be sure to allow free standing edges in snappyHexMesh if you do this and are using SHM.
Dano62 is offline   Reply With Quote

Old   September 14, 2022, 06:51
Default
  #5
New Member
 
Eric Segalerba
Join Date: Dec 2021
Location: Italy
Posts: 24
Rep Power: 4
EricS is on a distinguished road
Ok thank you, I was struggling to make it in snappy. In your way it looks much easier!
EricS is offline   Reply With Quote

Reply

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
Table bounds warnings at: END OF TIME STEP CFXer CFX 4 July 17, 2020 00:44
Basic Nozzle-Expander Design karmavatar CFX 20 March 20, 2016 09:44
Problem in setting Boundary Condition Madhatter92 CFX 12 January 12, 2016 05:39
Problem on boundary conditions wiiketi FLUENT 1 March 25, 2014 00:06
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 08:00


All times are GMT -4. The time now is 19:33.