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

Smearing interface for capillary flows using interfoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 12, 2016, 12:47
Default Smearing interface for capillary flows using interfoam
  #1
New Member
 
Join Date: Sep 2015
Posts: 2
Rep Power: 0
m.b. is on a distinguished road
Dear Foamers,

Iím completely new at OpenFoam and also new in this Forum, so first I want to say hello to everyone!

Iím working on capillary flows in nano structures and for this I am using the InterFoam solver.
Unfortunately I get strange results for my liquid-Air interface and for this reason I would like to ask for some explanations or hints for this.

My simulation case is describing a 2D rectangular horizontal capillary tube with a T-junction in it and also with different contact angles at the top and bottom wall. Unfortunately my liquid interface is smearing with the time so I am not able to detect the position of the ďrealĒ interface.

I have tried different discretization schemes to getting better results but as before the result is not satisfying.

It would be very nice if someone take a look at my case and give me some helpful hints to fix this.

I will post my simulation files below this message and also in this link:

https://www.dropbox.com/sh/swneqv0un...RRPWqyP4a?dl=0

Regards,

M.



My initial conditions:
-----------------------------------------------------------------------------------

For p:

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
Inlet
{
type fixedValue;
value uniform 0;
}

Outlet
{
type fixedValue;
value uniform 0;
}

bottom
{
type fixedFluxPressure;
}

top
{
type fixedFluxPressure;
}

Front
{
type empty;
}

Back
{
type empty;
}
}
----------------------------------------------------------------------------

For u:

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
Inlet
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}

Outlet
{
type zeroGradient;
}

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

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

Front
{
type empty;
}

Back
{
type empty;
}
}
--------------------------------------------------------------------

For alpha:

dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{
Inlet
{
type inletOutlet;
value uniform 1;
inletValue uniform 1;
}

bottom
{
type constantAlphaContactAngle;
theta0 50;
limit gradient;
value uniform 0;
}

top
{
type constantAlphaContactAngle;
theta0 35;
limit gradient;
value uniform 0;
}

Outlet
{
type zeroGradient;
}

Front
{
type empty;

Back
{
type empty;
}
}

----------------------------------------------------------

My fvSchemes file:

ddtSchemes
{
default Euler;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
div(rhoPhi,U) Gauss vanLeer;
div(phi,alpha) Gauss vanLeer01;
div(phirb,alpha) Gauss interfaceCompression;
div((muEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p_rgh;
pcorr;
alpha.water;
}


----------------------------------------------------------

My fvSolution file:

solvers
{
alpha.water
{
nAlphaCorr 3;
nAlphaSubCycles 3;
cAlpha 1;
}

pcorr
{
solver PCG;
preconditioner DIC;
tolerance 1e-12;
relTol 0;
}

p_rgh
{
solver PCG;
preconditioner DIC;
tolerance 1e-12;
relTol 0.05;
}

p_rghFinal
{
$p_rgh;
tolerance 1e-12;
relTol 0;
}

U
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-12;
relTol 0;
}
}

PIMPLE
{
momentumPredictor no;
nCorrectors 7;
nNonOrthogonalCorrectors 0;
}
Attached Images
File Type: png alpha1.png (19.8 KB, 51 views)
File Type: png alpha2.png (23.6 KB, 37 views)
File Type: png alpha3.png (24.4 KB, 43 views)
m.b. is offline   Reply With Quote

Old   March 14, 2016, 15:05
Default
  #2
New Member
 
Join Date: Sep 2015
Posts: 2
Rep Power: 0
m.b. is on a distinguished road
Does somebody have any idea?
m.b. is offline   Reply With Quote

Old   April 19, 2017, 13:29
Default
  #3
New Member
 
rasoul
Join Date: Mar 2015
Posts: 24
Rep Power: 10
rasool_soofi is on a distinguished road
Hi m.b
I have this problem, too and I cant solve it yet.
As far as I know, this problem refers to two factors:
1- the version of OpenFOAM or Foam-ext software.
2- the size of your channel and meshes.
rasool_soofi is offline   Reply With Quote

Reply

Tags
capillary flow, 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
Wind turbine simulation Saturn CFX 58 July 3, 2020 01:13
Implementing a stress on a portion of interface when using interFoam Kleiber OpenFOAM Programming & Development 0 May 18, 2014 00:55
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 04:05
Smoothing interface of two-phase flows aiqch Main CFD Forum 0 June 18, 2009 13:31
Smoothing interface of two-phase flows aiqch Main CFD Forum 0 June 18, 2009 13:26


All times are GMT -4. The time now is 12:35.