# Curvature of the interface using interFoam

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

March 1, 2011, 14:16
Curvature of the interface using interFoam
#1
Senior Member

Andrea Ferrari
Join Date: Dec 2010
Posts: 313
Rep Power: 9
Hi all,

I am using interFoam to simulate multiphase flow in different geometries and I tried to write a little piece of code (writeCurvature.C that i call as a post processing) to calculate the curvature of the interface following the same method that openFoam normally use in the VOF method. (the calculations that i have done are identical to those in interfaceProperties.C). The problem is that I get strange result. I expected to find high values near the interface and values close to zero in the bulk but it is not!

i attached 2 jpg. one is the colour plot of alpha1 and the other one is the colour plot of the curvature that i got.
Any suggestion?

thanks
andrea
Attached Images
 alpha.png (21.7 KB, 122 views) curvature.jpg (41.9 KB, 130 views)

 March 2, 2011, 05:19 #2 Senior Member   Andrea Ferrari Join Date: Dec 2010 Posts: 313 Rep Power: 9 someone who has any idea why I find this curvature? These are the calculations I've done: Code: ```reading alpha1 first then // now calc the gradients volVectorField gradAlpha1 = fvc::grad(alpha1); Info<< "Writing gradAlpha1 to " << gradAlpha1.name() << " in " << runTime.timeName() << endl; gradAlpha1.write(); // Interpolated face-gradient of alpha surfaceVectorField gradAlphaf = fvc::interpolate(gradAlpha1); dimensionedScalar deltaN_ = 1e-8/pow(average(alpha1.mesh().V()), 1.0/3.0); // Face unit interface normal surfaceVectorField nHatfv = gradAlphaf/(mag(gradAlphaf) + deltaN_); // Face unit interface normal flux const surfaceVectorField& Sf = mesh.Sf(); surfaceScalarField nHatf_ = nHatfv & Sf; volScalarField curvature ( IOobject ( "Curvature", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), -fvc::div(nHatf_) );``` The thing that I can not understand is if I'm doing something wrong or if it is normal to find such values thanks andrea

 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 Saturn CFX 45 February 8, 2016 05:42 AnjaMiehe OpenFOAM Running, Solving & CFD 8 June 14, 2010 06:49 Pankaj CFX 9 November 23, 2009 05:05 Mark CFX 6 November 15, 2004 16:55 lego CFX 3 November 5, 2002 21:09

All times are GMT -4. The time now is 18:33.