CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > ANSYS > CFX

Boundary Conditions for liquid Drainage in a foam

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

LinkBack Thread Tools Search this Thread Display Modes
Old   August 27, 2019, 05:56
Default Boundary Conditions for liquid Drainage in a foam
New Member
Join Date: Mar 2018
Posts: 15
Rep Power: 4
dickes is on a distinguished road
Hi everyone,

I am having problems implementing the boundary conditions for my openings properly.

First of all: I am simulating liquid drainage in a foam with a transient Multiphase approach. As I am not planning to resolve the interfaces between liquid and gas bubbles I am using an Euler-Euler approach in CFX. Instead I am using momentum sources for each phase to model the flow. Right now my domain is quasi 1-dimensional (1 element in x and z which are Symmetry planes, as well as 100 elements in y) and gravity is acting in the y-direction. I initiliase the domain with liquid=0.05 (continuous Fluid) and gas=0.95 (Dispersed Fluid)

My momentum Sources are added by a subdomain for the liquid Phase as follows:
1. Isotropic Loss Model with a permeability expression, which is dependent on liquid Volume Fraction
Foam permeability = (edge length)^2*6.3*10^-3
edge length = (1-liquid.Volume Fraction)^(-1/3)*2^(-5/3)*3^(1/6)*pi^(1/3)*d gas

2. Momentum source for the capillarity:
capillary y = -surface tension/(2*0.33*diameter_gas/2)*liquid.Volume Fraction^(-3/2)*(liquid.Volume Fraction.Gradient Y)*liquid.Volume Fraction^(2)

On top of the domain I have a pressure Opening with an Opening Pressure = 0Pa; liquid.Volume fraction = 0; gas.Volume Fraction = 1

Below I also have an opening where the liquid is able to drain through, and the boundary conditions are given by liquid.Volume Fraction= 0.05; gas Volume Fraction = 0.95 and I enter the velocity in the form of

Foam permeability *(liquid.density*gravity y+capillary y)/liquid.viscosity*liquid.Volume Fraction

Which must be the velocity according to my momentum sources.

I expect the liquid to drain from the top but stay at the liquid.Volume Fraction of 0.05 at the lower boundary, which I defined as a Boundary Condition. Through the increase of the gradient at the boundary my capillary term should be growing and reduce the velocity for the outlet, untill an equilibrium is reached. However, my liquid Volume Fraction distribution will not reach an expected equilibrium. The Volume Fraction decrease below the given values set by my boundary condition and if I let the simulation run long enough, liquid will completely drain out of the Domain.

Thanks in advance for all helpful tips!
dickes is offline   Reply With Quote


boundary conditions

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-reverse flow in outlet lesds to a mass in flow field xiexing CFX 3 March 29, 2017 11:00
Problem in setting Boundary Condition Madhatter92 CFX 12 January 12, 2016 05:39
Overflow Error in Multiphase Modelling with Two Continuous Fluids ashtonJ CFX 6 August 11, 2014 15:32
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
[blockMesh] BlockMesh FOAM warning gaottino OpenFOAM Meshing & Mesh Conversion 7 July 19, 2010 15:11

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