convergence in tetramesh

khaiching · October 27, 2010, 22:30

Hi colleagues:

I would like to gather comments from you all regarding the convergence performance using tetrahedral mesh.

I build up my 3D tetrahedral mesh using Delaunay method and find that the convergence is kinda irregular. For a simple incompressible flow condition in a flow distributor at a certain tetrahedral mesh setup, I could not achieve convergence at all. The pressure correction residual error is getting unstable ( residual error goes down to a certain level and it spikes up again).

Subsequently, I extract the cells (problematic cell) that contain the largest residual error and try to remesh the region of those problematic cells. In this case, the convergence can be attained.

Dear colleagues, I would like to know whether you have come across this phenomenon before? Something come across my mind: for the matrix solver of pressure correction equation, I utilise the Bicgstab method with preconditioner. Could this problem due to the characteristic of the Bicgstab method itself, whereby the error cannot be pushed out from the domain as efficiently as the AMG method, particularly in the unstructured tetrahedral environment? I have not tried the AMG yet.

Or, could it be due to any other reasons, which is unknown to me?

Thanks.

eplanting · October 28, 2010, 00:01

Variable dt or fixed dt? Also, are you using upwinding or midpoint, midpoint is a bit more accurate.

khaiching · October 28, 2010, 02:57

Hi:

I am using steady-state algorithm to solve the equations. SIMPLE algorithm in fact.

I am using 1st order upwind scheme for all the convective terms. Therefore, I think the instability occured is not due to my convection scheme.

eplanting · October 28, 2010, 06:10

ick, while upwinding is powerful, it tends to give more pretty pictures than accurate ones. Anyways, I still think it's probably something to do with a value of dt occuring that's too high for the solution, that's a pretty common problem around the graduate level at least, one gradient is pretty steep, and it looks like the rest of the solution is converging, but one area is developing an error, anyway, I'd recheck ur code for the calculation of dt, or you could post it.

October 27, 2010, 22:30	convergence in tetramesh	#1
khaiching New Member KaySee Join Date: Apr 2010 Posts: 5 Rep Power: 16	Hi colleagues: I would like to gather comments from you all regarding the convergence performance using tetrahedral mesh. I build up my 3D tetrahedral mesh using Delaunay method and find that the convergence is kinda irregular. For a simple incompressible flow condition in a flow distributor at a certain tetrahedral mesh setup, I could not achieve convergence at all. The pressure correction residual error is getting unstable ( residual error goes down to a certain level and it spikes up again). Subsequently, I extract the cells (problematic cell) that contain the largest residual error and try to remesh the region of those problematic cells. In this case, the convergence can be attained. Dear colleagues, I would like to know whether you have come across this phenomenon before? Something come across my mind: for the matrix solver of pressure correction equation, I utilise the Bicgstab method with preconditioner. Could this problem due to the characteristic of the Bicgstab method itself, whereby the error cannot be pushed out from the domain as efficiently as the AMG method, particularly in the unstructured tetrahedral environment? I have not tried the AMG yet. Or, could it be due to any other reasons, which is unknown to me? Thanks.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Convergence of CFX field in FSI analysis	nasdak	CFX	2	June 29, 2009 02:17
Problems with convergence with an easy system	franzdrs	Main CFD Forum	0	June 15, 2009 19:17
increasing mesh quality is leading to poor convergence	tippo	CFX	2	May 5, 2009 11:55
Defect correction and convergence	ganesh	Main CFD Forum	4	June 30, 2006 15:20
Convergence problems	Chetan	FLUENT	3	April 15, 2004 20:13

October 28, 2010, 00:01		#2
eplanting New Member eric planting Join Date: Oct 2010 Posts: 3 Rep Power: 16	Variable dt or fixed dt? Also, are you using upwinding or midpoint, midpoint is a bit more accurate.

October 28, 2010, 02:57		#3
khaiching New Member KaySee Join Date: Apr 2010 Posts: 5 Rep Power: 16	Hi: I am using steady-state algorithm to solve the equations. SIMPLE algorithm in fact. I am using 1st order upwind scheme for all the convective terms. Therefore, I think the instability occured is not due to my convection scheme.

October 28, 2010, 06:10		#4
eplanting New Member eric planting Join Date: Oct 2010 Posts: 3 Rep Power: 16	ick, while upwinding is powerful, it tends to give more pretty pictures than accurate ones. Anyways, I still think it's probably something to do with a value of dt occuring that's too high for the solution, that's a pretty common problem around the graduate level at least, one gradient is pretty steep, and it looks like the rest of the solution is converging, but one area is developing an error, anyway, I'd recheck ur code for the calculation of dt, or you could post it.