# CLSVOF in InterFoam

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

 February 11, 2014, 13:54 CLSVOF in InterFoam #1 New Member   Deepak Kalaikadal Join Date: Feb 2014 Posts: 6 Rep Power: 3 Hi, 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 solve ( fvm::ddt(lsPhi) +(sign(lsPhi) * mag(fvc::grad(lsPhi))) -sign(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?

 March 16, 2014, 11:11 #2 New Member   adong Join Date: Jul 2012 Posts: 1 Rep Power: 0 Hi Try this 1- you define the field d and d0 (d is your lsPhi) volScalarField d ( IOobject ( "d", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); volScalarField d0 ( IOobject ( "d0", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE ), mesh, dimensionedScalar("d0",dimensionSet(0,0,0,0,0,0,0) ,scalar(0)), d.boundaryField().types() ); 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) ); dEqn.solve(); } nimasam, Martin_K_lalelu, kalaikdr and 1 others like this.

 November 27, 2014, 12:50 Clsvof #3 New Member   praveen kumar sharma Join Date: Jun 2014 Posts: 14 Rep Power: 3 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.

 Tags clsvof, openfoam

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post gooya_kabir OpenFOAM Running, Solving & CFD 0 December 9, 2013 06:09 keepfit OpenFOAM 3 August 29, 2013 11:21 francesco_b OpenFOAM Running, Solving & CFD 8 July 31, 2013 02:29 DanM OpenFOAM Running, Solving & CFD 11 January 5, 2013 07:21 sxhdhi OpenFOAM Running, Solving & CFD 3 May 5, 2009 21:58

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