CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

How to simulate on geometry with curvature

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 14, 2012, 06:51
Default How to simulate on geometry with curvature
  #1
Member
 
Join Date: Aug 2011
Posts: 33
Rep Power: 5
seboxx is on a distinguished road
Hi Foamers,

I have created a bump using blockMesh.
Visually everything works out nice, I mean I can see the bump in ParaFoam and it looks as i intended it to look. (see attached jpg)
But the problem is, the simulations i run are not correct. Actually it is just the same as if I would just run it on a disc with equal radius.

For example I need to calculate the volume of the domain at some point during the simulation. No matter whether I use domainIntegrate, volumeIntegrate or sum(mesh.V().field()) I get the volume 1256.62, which is exactly pi * r^2
At first I thought this is a problem of the domainIntegrate function or whatever.
But doing simulations of localized structures that should definitely feel the curvature it just goes through as if it would be flat.

I am using pisoFoam, i think I have just modified the pisoFoam source code to my needs... Is it possible that I need additional libraries to calculate on curved meshes? Do I need another solver as a base code?

Also checkMesh reports that the Mesh is ok and the solver doesn't give any error message...

I am happy about any suggestion


Here is the code:
Code:
#include "fvCFD.H"


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])
{
    #include "setRootCase.H"

    #include "createTime.H"
    #include "createMesh.H"
    #include "createFields.H"
  //#include "initContinuityErrs.H"
  //#include "fvcVolumeIntegrate.H"    
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

    Info<< "\nStarting time loop\n" << endl;


    while (runTime.loop())
      {
        Info<< "Time = " << runTime.timeName() << nl << endl;

	#include "readPISOControls.H"
	
        // --- PISO loop 
	
        for (int corr=0; corr<nCorr; corr++)
        {
	  // #include "uEqn.H"
	  // #include "vEqn.H"  	  
	}
	   
        Info<< "sum(mesh.V().field())= " << sum(mesh.V().field()) << endl;
	Info<< "domainIntegrate= " << fvc::domainIntegrate(pos(u-(1e-12))) <<endl;
	Info<< "volumeIntegrate= " << sum(fvc::volumeIntegrate(pos(u-(1e-12)))) <<endl;
	
	//	outputDict.set("S",scale_spaceunit*fvc::volumeIntegrate(pos(u-(1e-12))));
        runTime.write();
	
	    Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
            << "  ClockTime = " << runTime.elapsedClockTime() << " s"
            << nl << endl;
      }

    Info<< "End\n" << endl;

    return 0;
}
Best regards,

seb
Attached Images
File Type: jpg bump.jpg (70.3 KB, 514 views)
seboxx is offline   Reply With Quote

Old   January 18, 2012, 05:43
Default
  #2
Member
 
Join Date: Aug 2011
Posts: 33
Rep Power: 5
seboxx is on a distinguished road
Hmm, or is it a problem with my blockMeshDict?

I have tried the same geometry with modified Versions of icoFoam and simpleFoam. The modification was only
Code:
sum(mesh.V().field())
and the result was the same. The volume of a flat plate...

Hmm checkMesh says Mesh Ok. But I realised that the wrong total volume already appears there...

Please find my blockMeshDict attached. I would be happy about any suggestions...

Regards

seb
Attached Files
File Type: zip blockMeshDict.zip (1.3 KB, 7 views)
seboxx is offline   Reply With Quote

Reply

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
Problem Importing Geometry ProE to CFX fatb0y CFX 3 January 14, 2012 19:42
2D geometry to 3D geometry Sgonzalezg ANSYS Meshing & Geometry 11 April 8, 2011 11:02
importing geometry into starccm krishna086 STAR-CCM+ 0 November 28, 2010 21:54
car geometry Arjun FLUENT 4 September 1, 2007 15:08
urgent help! car geometry Arjun FLUENT 0 August 13, 2007 13:07


All times are GMT -4. The time now is 20:41.