[Eller_OF]

Hello,

I have a question on how to change the contact angle using the fixedFluxPressure BC on the Walls. The simulation has been sent to me and I am supposed to only change the angle from now (I’ve been told) 45 degrees to 70 degrees.

Usually I would use

Walls
type constantAlphaContactAngle;
theta0 45;
limit gradient;
value uniform 0;

In 0/U folder and then just change the theta0.

But the simulation i work with has a

Walls
type fixedFluxPressure;

So now i am unfortunately unsure about what to do.

To sum up: The simulation currently uses contact angle = 45 degrees and I cannot find where this has been set.

So if anyone has an idea i would be super thankful.

(Using interfoam)

[sadsid]

Originally Posted by Eller_OF

Hello,

I have a question on how to change the contact angle using the fixedFluxPressure BC on the Walls. The simulation has been sent to me and I am supposed to only change the angle from now (I’ve been told) 45 degrees to 70 degrees.

Usually I would use

Walls
type constantAlphaContactAngle;
theta0 45;
limit gradient;
value uniform 0;

In 0/U folder and then just change the theta0.

But the simulation i work with has a

Walls
type fixedFluxPressure;

So now i am unfortunately unsure about what to do.

To sum up: The simulation currently uses contact angle = 45 degrees and I cannot find where this has been set.

So if anyone has an idea i would be super thankful.

(Using interfoam)

Solved your issue?
I am using this condition and want to know one thing. I have to set GDL contact angle. How to set this? Any idea

[Ash Kot]

Originally Posted by Eller_OF

Hello,

I have a question on how to change the contact angle using the fixedFluxPressure BC on the Walls. The simulation has been sent to me and I am supposed to only change the angle from now (I’ve been told) 45 degrees to 70 degrees.

Usually I would use

Walls
type constantAlphaContactAngle;
theta0 45;
limit gradient;
value uniform 0;

In 0/U folder and then just change the theta0.

But the simulation i work with has a

Walls
type fixedFluxPressure;

So now i am unfortunately unsure about what to do.

To sum up: The simulation currently uses contact angle = 45 degrees and I cannot find where this has been set.

So if anyone has an idea i would be super thankful.

(Using interfoam)

The following code should provide you a brief idea of setting up contact angle along with fixedFluxPressure

The code in alpha.water file:

wall
    {
        type            constantAlphaContactAngle;
        theta0         45;
        limit            gradient;
        value           uniform 0;
    }

The code in p_rgh file for that same wall patch:

wall
    {
        type            fixedFluxPressure;
        value           uniform 0;
    }

The code in U file, the same wall patch:
Either:

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

or

wall
    {
        type            noSlip;
    }

I hope this solves your question:
As you changes the contact angle the liquid behavior also changes, you can verify based upon the value computed in paraview for wetted area, and interfacial area using integrateVariables filter.

[sadsid]

Originally Posted by Ash Kot

The following code should provide you a brief idea of setting up contact angle along with fixedFluxPressure

The code in alpha.water file:

wall
    {
        type            constantAlphaContactAngle;
        theta0         45;
        limit            gradient;
        value           uniform 0;
    }

The code in p_rgh file for that same wall patch:

wall
    {
        type            fixedFluxPressure;
        value           uniform 0;
    }

The code in U file, the same wall patch:
Either:

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

or

wall
    {
        type            noSlip;
    }

I hope this solves your question:
As you changes the contact angle the liquid behavior also changes, you can verify based upon the value computed in paraview for wetted area, and interfacial area using integrateVariables filter.

Actually, I am applying the same condition. I am solving channel flow and the article for which I am trying to get the results used: 1) Hydrophobic GDL (contact angle 140) and 2) hydrophilic channel walls with a contact angle 60. The bottom wall is considered as GDL wall. Now I want to know how to set these two angles in alpha.water

[Ash Kot]

Trying to understand what you have mentioned:
so bottom wall is Hydrophobic and rest all walls are Hydrophilic is that correct?

In that case, the bottom wall will have the above mentioned code related to constantContactAngle with value defined 140.
and rest channel walls with have same code of constantContactAngle with value as 60
and you define that in alpha.water file.

Only one constantContactAngle boundary condition can be defined for one wall with one specific value of whatever the angle is.
multiple constant contact angles are not possible for one single patch wall, in one single boundary condition.
if you want the value for that constant contact angle needed to be changed then you'll probably need to use swak4Foam to setup time dependent contact angle changing BC.

[sadsid]

Originally Posted by Ash Kot

Trying to understand what you have mentioned:
so bottom wall is Hydrophobic and rest all walls are Hydrophilic is that correct?

In that case, the bottom wall will have the above mentioned code related to constantContactAngle with value defined 140.
and rest channel walls with have same code of constantContactAngle with value as 60
and you define that in alpha.water file.

Only one constantContactAngle boundary condition can be defined for one wall with one specific value of whatever the angle is.
multiple constant contact angles are not possible for one single patch wall, in one single boundary condition.
if you want the value for that constant contact angle needed to be changed then you'll probably need to use swak4Foam to setup time dependent contact angle changing BC.

Yes, the bottom wall of the channel corresponds to the GDL, which was defined as a hydrophobic surface with a contact angle of 140. So for this patch I can use:
type constantAlphaContactAngle;
theta0 140;
limit gradient;
value uniform 0;

and for top I can use:
type constantAlphaContactAngle;
theta0 60;
limit gradient;
value uniform 0;

This is what you are saying?

[sadsid]

Originally Posted by Ash Kot

Trying to understand what you have mentioned:
so bottom wall is Hydrophobic and rest all walls are Hydrophilic is that correct?

In that case, the bottom wall will have the above mentioned code related to constantContactAngle with value defined 140.
and rest channel walls with have same code of constantContactAngle with value as 60
and you define that in alpha.water file.

Only one constantContactAngle boundary condition can be defined for one wall with one specific value of whatever the angle is.
multiple constant contact angles are not possible for one single patch wall, in one single boundary condition.
if you want the value for that constant contact angle needed to be changed then you'll probably need to use swak4Foam to setup time dependent contact angle changing BC.

The goemetry is like the one I attached here.

[Ash Kot]

After looking your geometry, and the code that you have defined post #6, yes that is what I'm saying.

its pretty simple, although I don't know how much variation you might be able to visualize, based upon the contact angles.

The hydrophilic and hydrophobic energies are too small, they define their significance on 100 molecular to 1,000,000 molecular level (if we assume a drop of water is 1,000,000 molecules).
However, 10,000,000 to 100,000,000 molecular level (that is ice cube size blob of water), gravitational energy along with density and viscosity show more predominant effect and over constrain the hydrophilic and hydrophobic energies. Thus, in your CFD analysis, there is significant possibility that even at 140 or more than that angle, you won't be able to see much difference unless your mesh is finest!

Apologies, discussion in molecules and other chemical concepts, being chemical engineer it becomes evident to explain in molecules and moles :-D

[sadsid]

Originally Posted by Ash Kot

After looking your geometry, and the code that you have defined post #6, yes that is what I'm saying.

its pretty simple, although I don't know how much variation you might be able to visualize, based upon the contact angles.

The hydrophilic and hydrophobic energies are too small, they define their significance on 100 molecular to 1,000,000 molecular level (if we assume a drop of water is 1,000,000 molecules).
However, 10,000,000 to 100,000,000 molecular level (that is ice cube size blob of water), gravitational energy along with density and viscosity show more predominant effect and over constrain the hydrophilic and hydrophobic energies. Thus, in your CFD analysis, there is significant possibility that even at 140 or more than that angle, you won't be able to see much difference unless your mesh is finest!

Apologies, discussion in molecules and other chemical concepts, being chemical engineer it becomes evident to explain in molecules and moles :-D

Thank you for your detailed reply. At the moment I am only trying to obtaint he benchmark results of the paper "Residence time of water film and slug flow features in fuel cell gas
channels and their effect on instantaneous area coverage ratio".

For my own case, I will follow what you said. If I need more guidance I will post here Thank you for your response.

[Ash Kot]

No problem, your welcome.

Just for the side-note:
In their article, table 2 the summary of settings for VoF model.
The interface compression/representation scheme they used: Piecewise Linear Interface Calculation (PLIC)
Thus, if you planning to use that scheme then you'll have to use OpenFOAM V8 version (Release 2020) from OpenFOAM.org. (I have less idea about ESI OpenFOAM having that scheme).

Only OpenFOAM V8 has PLIC scheme otherwise you'll probably use MULES for alpha.water (the default scheme)

[sadsid]

Originally Posted by Ash Kot

No problem, your welcome.

Just for the side-note:
In their article, table 2 the summary of settings for VoF model.
The interface compression/representation scheme they used: Piecewise Linear Interface Calculation (PLIC)
Thus, if you planning to use that scheme then you'll have to use OpenFOAM V8 version (Release 2020) from OpenFOAM.org. (I have less idea about ESI OpenFOAM having that scheme).

Only OpenFOAM V8 has PLIC scheme otherwise you'll probably use MULES for alpha.water (the default scheme)

You are right. I am using OF-v8.