inviscid 2d computation for M=0.3 flow past cylinder on fine mesh can not converge

boubalos · March 17, 2011, 09:16

I try to use very fine mesh to test my Euler solver, however I get the following surprising results. I wish someone can give me a explain.

Both my solver and Fluent can get a correct result on a coarse mesh. The Mach contour shows symmetrical. wake is very weak.

Both my solver and Fluent can not get a correct result on a fine mesh. The Mach contour shows unsymmetrical, and wake is strong and long. Fluent even can not get a converged result!

The diameter of the cylinder is 2.0. The mesh is structured. The first layer mesh is 0.001 for fine mesh; and 0.02 for coarse mesh. Far field boundary is 40 times of the diameter away. And I am quite sure that my solver and Fluent is correctly setted.

boubalos · March 20, 2011, 01:27

Can't any one help me, on such a simple problem ...?

praveen · March 20, 2011, 02:41

Describe the numerical scheme, post pictures of grid, solution. Then somebody might be able to help.

Martin Hegedus · March 20, 2011, 04:26

The issue may be that a subsonic Euler solution has very little dissipation. When the solution starts up a high pressure compression zone will occur at the front and a low pressure expansion zone occurs at the back side. The solution will then oscillate between the front and the back. Eventually these oscillations may die down. Maybe. The solution methodology may keep it going. If you are just testing out your code by comparing it to Fluent you can run it at a Mach number high enough to get a shock. Maybe Mach 0.7. The shock will isolate the back from the front. But it still may take time to converge.

March 17, 2011, 09:16	inviscid 2d computation for M=0.3 flow past cylinder on fine mesh can not converge	#1
boubalos New Member Zhu Xiaofeng Join Date: Sep 2009 Posts: 4 Rep Power: 16	I try to use very fine mesh to test my Euler solver, however I get the following surprising results. I wish someone can give me a explain. Both my solver and Fluent can get a correct result on a coarse mesh. The Mach contour shows symmetrical. wake is very weak. Both my solver and Fluent can not get a correct result on a fine mesh. The Mach contour shows unsymmetrical, and wake is strong and long. Fluent even can not get a converged result! The diameter of the cylinder is 2.0. The mesh is structured. The first layer mesh is 0.001 for fine mesh; and 0.02 for coarse mesh. Far field boundary is 40 times of the diameter away. And I am quite sure that my solver and Fluent is correctly setted.

Thread Tools	Search this Thread
Show Printable Version Email this Page	Search this Thread: Advanced Search
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
FLUENT, flow past a circular cylinder	Paul Reichl	Main CFD Forum	16	December 8, 2012 15:21
Flow past a cylinder	ckl	Main CFD Forum	2	October 27, 1999 12:36
computation about flow around a yawed cone	Tylor Xie	Main CFD Forum	0	June 9, 1999 07:33
Free-surface flow past a cylinder, help !	Tony Chen	Main CFD Forum	0	February 5, 1999 02:06
Simulation of the Flow past a circular cylinder using STAR-CD	M. S. GUEROUACHE	Main CFD Forum	0	October 1, 1998 10:51

March 20, 2011, 01:27		#2
boubalos New Member Zhu Xiaofeng Join Date: Sep 2009 Posts: 4 Rep Power: 16	Can't any one help me, on such a simple problem ...?

March 20, 2011, 02:41		#3
praveen Super Moderator Praveen. C Join Date: Mar 2009 Location: Bangalore Posts: 342 Blog Entries: 6 Rep Power: 18	Describe the numerical scheme, post pictures of grid, solution. Then somebody might be able to help.

March 20, 2011, 04:26		#4
Martin Hegedus Senior Member Martin Hegedus Join Date: Feb 2011 Posts: 500 Rep Power: 19	The issue may be that a subsonic Euler solution has very little dissipation. When the solution starts up a high pressure compression zone will occur at the front and a low pressure expansion zone occurs at the back side. The solution will then oscillate between the front and the back. Eventually these oscillations may die down. Maybe. The solution methodology may keep it going. If you are just testing out your code by comparing it to Fluent you can run it at a Mach number high enough to get a shock. Maybe Mach 0.7. The shock will isolate the back from the front. But it still may take time to converge.