fvSchemes in FOAM extend 3.1

Shenan · December 12, 2014, 11:34

I tried running the LRR turbulence model in OF 2.2.1 but ran into difficulty with time step continuity errors blowing up. I tightened the pressure convergence criteria but still it blew up; so I thought I'd give FOAM 3.1 a try.

I thought most things would be portable but some of the fvSchemes that I was using in OF are not recognised by FOAM. For instance, in OF I had:

Code:

div(phi,U)      bounded Gauss linearUpwindV grad(U);

In FOAM,

Code:

bounded

was not recognised; so I removed it and then

Code:

grad(U)

was not recognised. The problem is that I don't know what to replace it with.

Is there a "dictionary" to "translate" between OF and FOAM?

Tushar@cfd · December 16, 2014, 23:26

Quote:

Originally Posted by Shenan

I tried running the LRR turbulence model in OF 2.2.1 but ran into difficulty with time step continuity errors blowing up. I tightened the pressure convergence criteria but still it blew up; so I thought I'd give FOAM 3.1 a try.

I thought most things would be portable but some of the fvSchemes that I was using in OF are not recognised by FOAM. For instance, in OF I had:

Code:

div(phi,U)      bounded Gauss linearUpwindV grad(U);

In FOAM,

Code:

bounded

was not recognised; so I removed it and then

Code:

grad(U)

was not recognised. The problem is that I don't know what to replace it with.

Is there a "dictionary" to "translate" between OF and FOAM?

Hello Shenan,

If I am correct then the "bounded" term was included in OpenFOAM from OF ver. 2.2 onwards. This term will be beneficial if you are running steady-state solver. It will bound your solution and help achieve better convergence.

-
Best Regards!

hjasak · December 17, 2014, 11:39

Try

Gauss linearUpwindV Gauss linear;

Shenan · December 22, 2014, 09:45

Thanks for the tip, Hrv. Unfortunately, the solver did not converge.

The time step continuity errors (sum local) went from 0.01 in the second iteration to 2e+05 in the third.

Of the other parameters, only bounding epsilon (max) showed signs of divergence: from 470 in the second iteration to 2e+08 in the third.

I will have a look at this again after Xmas. It might not be a problem with fvSchemes.

December 12, 2014, 11:34	fvSchemes in FOAM extend 3.1	#1
Shenan New Member Shenan Join Date: Nov 2014 Posts: 5 Rep Power: 11	I tried running the LRR turbulence model in OF 2.2.1 but ran into difficulty with time step continuity errors blowing up. I tightened the pressure convergence criteria but still it blew up; so I thought I'd give FOAM 3.1 a try. I thought most things would be portable but some of the fvSchemes that I was using in OF are not recognised by FOAM. For instance, in OF I had: Code: div(phi,U) bounded Gauss linearUpwindV grad(U); In FOAM, Code: bounded was not recognised; so I removed it and then Code: grad(U) was not recognised. The problem is that I don't know what to replace it with. Is there a "dictionary" to "translate" between OF and FOAM?

December 17, 2014, 11:39		#3
hjasak Senior Member Hrvoje Jasak Join Date: Mar 2009 Location: London, England Posts: 1,905 Rep Power: 33	Try Gauss linearUpwindV Gauss linear; __________________ Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[Other] mesh airfoil NACA0012	anand_30	OpenFOAM Meshing & Mesh Conversion	13	March 7, 2022 17:22
[blockMesh] error message with modeling a cube with a hold at the center	hsingtzu	OpenFOAM Meshing & Mesh Conversion	2	March 14, 2012 09:56
[blockMesh] BlockMesh FOAM warning	gaottino	OpenFOAM Meshing & Mesh Conversion	7	July 19, 2010 14:11
[blockMesh] Axisymmetrical mesh	Rasmus Gjesing (Gjesing)	OpenFOAM Meshing & Mesh Conversion	10	April 2, 2007 14:00
[Gmsh] Import gmsh msh to Foam	adorean	OpenFOAM Meshing & Mesh Conversion	24	April 27, 2005 08:19

December 22, 2014, 09:45		#4
Shenan New Member Shenan Join Date: Nov 2014 Posts: 5 Rep Power: 11	Thanks for the tip, Hrv. Unfortunately, the solver did not converge. The time step continuity errors (sum local) went from 0.01 in the second iteration to 2e+05 in the third. Of the other parameters, only bounding epsilon (max) showed signs of divergence: from 470 in the second iteration to 2e+08 in the third. I will have a look at this again after Xmas. It might not be a problem with fvSchemes.