Adaptive Mesh refinement for steady state solver

yhaomin2007 · July 17, 2013, 13:35

Hey, all,

I have a problem about adaptive mesh refinement. I looked around all the solver, it seems all the solver that can do AMR is transient solver. Does anyone know if there is a existing AMR steady solver?
I also tried to implement AMR into rhoSimplecFoam solver. The procedure was simple, but I met error. I added mesh.update() in the solver.
The error info is :

Quote:

Time = 10

Selected 455 cells for refinement out of 112000.
Refined from 112000 to 115185 cells.
Selected 0 split points out of a possible 455.
Execution time for mesh.update() = 8.01 s
time step continuity errors : sum local = 379.765, global = -189.418, cumulative = -212.852

--> FOAM FATAL ERROR:
field does not correspond to level 0 sizes: field = 115185 level = 112000

From function void GAMGAgglomeration::restrictField(Field<Type>& cf, const Field<Type>& ff, const label fineLevelIndex) const
in file lnInclude/GAMGAgglomerationTemplates.C at line 47.

FOAM aborting

It seems the solver can find the cells that need to refine, but it fails to find the split point. The refinement is not success. Does anyone have experience on this?

thank you in advance~

Akshay · July 18, 2013, 01:55

Hey!

Is your cacheAgglomeration off? Try that.

yhaomin2007 · July 18, 2013, 12:07

Yes, this is the problem. I worked it out.
thank you

atoof · August 20, 2013, 04:11

Dear Haomin,

Did you just add mesh.update() in your solver? Is it sufficient to have adaptive mesh refinement for any steady state solver?

Regards,

Hossein

yhaomin2007 · August 20, 2013, 11:47

Hi, you also have to use "dynamicMesh" class to construct mesh. You can easily find examples in solvers that used dynamicMesh.
And, yes, the sentence that make difference is "mesh.update()".

verboomj · April 12, 2017, 07:35

Hi,

I'm also trying to implement adaptive meshing in my simpleFoam solver, but I seem to get the following error when I include mesh.update(); in my Simple Loop.

error: 'class Foam::fvMesh' has no member named 'update'

I've been trying to look at pimpleDyMFoam and use this as a starting point for my simpleFoam solver.

Any ideas how I can tackle this problem?

July 18, 2013, 01:55		#2
Akshay Member Akshay Kumar Join Date: Aug 2010 Location: India Posts: 84 Rep Power: 16	Hey! Is your cacheAgglomeration off? Try that. hua1015 likes this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Adding adaptive mesh refinement to the forwardStep tutorial	AlphaSheep	OpenFOAM Running, Solving & CFD	10	August 3, 2019 15:40
[snappyHexMesh] snappyHexMesh refinement regions ignored	guitarbren	OpenFOAM Meshing & Mesh Conversion	2	April 9, 2013 04:59
adaptive mesh refinement in STARCCM+	abdullahkarimi	Siemens	1	November 16, 2010 19:26
[snappyHexMesh] external flow with snappyHexMesh	chelvistero	OpenFOAM Meshing & Mesh Conversion	11	January 15, 2010 20:43
Mesh generator and CFD solver	Gennady Kireyko	Main CFD Forum	0	May 6, 2001 12:13

July 18, 2013, 12:07		#3
yhaomin2007 Member Haomin Yuan Join Date: Jan 2012 Location: Madison, Wisconsin, USA Posts: 59 Rep Power: 14	Yes, this is the problem. I worked it out. thank you

August 20, 2013, 04:11		#4
atoof Member Hossein Join Date: Apr 2010 Posts: 65 Rep Power: 16	Dear Haomin, Did you just add mesh.update() in your solver? Is it sufficient to have adaptive mesh refinement for any steady state solver? Regards, Hossein

August 20, 2013, 11:47		#5
yhaomin2007 Member Haomin Yuan Join Date: Jan 2012 Location: Madison, Wisconsin, USA Posts: 59 Rep Power: 14	Hi, you also have to use "dynamicMesh" class to construct mesh. You can easily find examples in solvers that used dynamicMesh. And, yes, the sentence that make difference is "mesh.update()".

April 12, 2017, 07:35		#6
verboomj New Member Jeroen Join Date: Oct 2016 Posts: 21 Rep Power: 10	Hi, I'm also trying to implement adaptive meshing in my simpleFoam solver, but I seem to get the following error when I include mesh.update(); in my Simple Loop. error: 'class Foam::fvMesh' has no member named 'update' I've been trying to look at pimpleDyMFoam and use this as a starting point for my simpleFoam solver. Any ideas how I can tackle this problem?