CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

CLSVOF in InterFoam

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

Like Tree4Likes
  • 4 Post By feddy

LinkBack Thread Tools Display Modes
Old   February 11, 2014, 13:54
Default CLSVOF in InterFoam
New Member
Deepak Kalaikadal
Join Date: Feb 2014
Posts: 6
Rep Power: 3
kalaikdr is on a distinguished road

I am trying to implement the CLSVOF scheme in interFoam. as given to Albadawi et al, International Journal of Multiphase Flow, 2013.

In his paper, it is given that the redistancing of the level set function is solved by using the equation

d(phi)/dt' = sign(phi)*(1-|grad(phi)|)

where t' is an artificial time step chosen as delta t' = 0.1*delta x

when I tried it on interFoam as

+(sign(lsPhi) * mag(fvc::grad(lsPhi)))

Foam compiles, but crashes on running due to dimension mismatch.

Any idea how to set it right? Does t' have the units of distance, as it is based on the mesh size?

Also, how do I solve this in a loop of artificial time until it converges to the steady state solution?
kalaikdr is offline   Reply With Quote

Old   March 16, 2014, 11:11
New Member
Join Date: Jul 2012
Posts: 1
Rep Power: 0
feddy is on a distinguished road
Try this

1- you define the field d and d0 (d is your lsPhi)

volScalarField d

volScalarField d0
dimensionedScalar("d0",dimensionSet(0,0,0,0,0,0,0) ,scalar(0)),

2- d have no dimension, for do not have one error of dimension you define:

dimensionedScalar one = dimensionedScalar("one", dimensionSet(0,1,0,0,0,0,0), 1.0);

3- you define the quantity present in the article:

const surfaceScalarField deltaCoeff = d.mesh().deltaCoeffs();
const scalar DX = (1.0/min(deltaCoeff).value());

scalar Dtau = 0.1*DX;
scalar Gamma = 0.75*DX;
scalar EPS = 0.1*DX;
scalar NITER = EPS/Dtau;

4- the pseudo iteration

d == Gamma*(2*alpha1-1.0);
d0 == d;

for (int loop = 0; loop < NITER; loop++)
// calculate distance function field

fvScalarMatrix dEqn
fvm::Sp(scalar(1),d) + Dtau*sign(d0)*mag(one*fvc::grad(d)) == d + Dtau*sign(d0)
feddy is offline   Reply With Quote

Old   November 27, 2014, 12:50
Default Clsvof
New Member
praveen kumar sharma
Join Date: Jun 2014
Posts: 14
Rep Power: 3
P Sharma is on a distinguished road
Hii everyone,
I am new in openFoam. I want to know how to use CLSVOF model in openFoam (or in foam extend) . Or I have to implement clsvof model separately. Is their any chance to install that model direct and i can use.
P Sharma is offline   Reply With Quote


clsvof, openfoam

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
interFoam in parallel gooya_kabir OpenFOAM Running, Solving & CFD 0 December 9, 2013 06:09
Problem of InterFoam with LES SpalartAllmarasIDDES keepfit OpenFOAM 3 August 29, 2013 11:21
InterFoam stops after deltaT goes to 1e14 francesco_b OpenFOAM Running, Solving & CFD 8 July 31, 2013 02:29
interFoam vs. simpleFoam channel flow comparison DanM OpenFOAM Running, Solving & CFD 11 January 5, 2013 07:21
Open Channel Flow using InterFoam type solver sxhdhi OpenFOAM Running, Solving & CFD 3 May 5, 2009 21:58

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