[verboomj]

Hello,

For the purpose of my solver, I want to refine the mesh at a certain edge, that is found by my topology optimization. This can be seen in the attached figure.

After some research, dynamicRefineFvMesh seems to be the 'go to' class. I've implemented the mesh.update() command in my simplFoam loop, similar to the pimpleDyMFoam implementation:

        fvc::makeAbsolute(phi, U);

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

        magAlphaGrad = mag(fvc::grad(alpha));
        
        counter += 1;
        if(counter>activateVolume)
        {
        mesh.update();
        }

It compiles and appears to work, but when I run a case it crashes at the moment it wants to refine the mesh.

Time = 2

Selected 468 cells for refinement out of 120000.
Refined from 120000 to 123276 cells.
Selected 0 split points out of a possible 468.
We're sorry, but the application crashed and safe stack tracing isn't available in this current implementation of blueCFD-Core patches for OpenFOAM.

The dynamicMeshDict file is rather standard as compared to exisiting tutorials for this dynamicRefineFvMesh.

Any ideas where I should look, regarding this crash?

[wyldckat]

Greetings verboomj,

I'm really late to give an answer to you on this, but I do have two suggestions for you:

1. If possible, upgrade to blueCFD-Core 2017, which provides stack tracing, so that you can then see why it crashed, instead of giving a cryptic message.
2. Old-school debugging is always an option, namely to use calls to "Info" to help isolate the location of the problem in the source code. For example:
   

   Info << "Got here 01" << endl;
   
   //...
   
   Info << "Got here 02" << endl;

   And so on.

Best regards,
Bruno

[Cambridge]

Originally Posted by verboomj

Hello,

For the purpose of my solver, I want to refine the mesh at a certain edge, that is found by my topology optimization. This can be seen in the attached figure.

After some research, dynamicRefineFvMesh seems to be the 'go to' class. I've implemented the mesh.update() command in my simplFoam loop, similar to the pimpleDyMFoam implementation:

        fvc::makeAbsolute(phi, U);

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

        magAlphaGrad = mag(fvc::grad(alpha));
        
        counter += 1;
        if(counter>activateVolume)
        {
        mesh.update();
        }

It compiles and appears to work, but when I run a case it crashes at the moment it wants to refine the mesh.

Time = 2

Selected 468 cells for refinement out of 120000.
Refined from 120000 to 123276 cells.
Selected 0 split points out of a possible 468.
We're sorry, but the application crashed and safe stack tracing isn't available in this current implementation of blueCFD-Core patches for OpenFOAM.

The dynamicMeshDict file is rather standard as compared to exisiting tutorials for this dynamicRefineFvMesh.

Any ideas where I should look, regarding this crash?

Hi verboomj

I am new in topology optimization of fluid and thermal by using OpenFOAM, could you please share any basic codes or tutorials of this? It would better for me to start to replicate some cases. Thank you so much.

regards