|
[Sponsors] |
March 1, 2011, 13:16 |
Curvature of the interface using interFoam
|
#1 |
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 16 |
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 |
|
March 2, 2011, 04:19 |
|
#2 |
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 16 |
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_) ); thanks andrea |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Wind turbine simulation | Saturn | CFX | 60 | July 17, 2024 05:45 |
interFoam - stratified flow - problem with shear stress at interface | AnjaMiehe | OpenFOAM Running, Solving & CFD | 8 | June 14, 2010 06:49 |
RPM in Wind Turbine | Pankaj | CFX | 9 | November 23, 2009 04:05 |
Convective Heat Transfer - Heat Exchanger | Mark | CFX | 6 | November 15, 2004 15:55 |
Replace periodic by inlet-outlet pair | lego | CFX | 3 | November 5, 2002 20:09 |