Roe-Solver Instability at Mach=1

Dwayne · July 6, 2013, 09:34

Hey Guys,

can somebody tell me the reasons why Roes solver become instable at Mach=1? The literature says that the solver becomes a central-differences-method at Mach=1 because the right side of the Roe-flux vanishes, is this true? Isn't it only one eigenvalue of the jacobian that becomes zero?

Thanks for help!

immortality · July 7, 2013, 05:00

hi
What software you use?

Dwayne · July 7, 2013, 07:36

Quote:

Originally Posted by immortality

hi
What software you use?

I use Matlab to calculate a one-dimensional shock-tube-problem for the Euler-Equations.

immortality · July 7, 2013, 08:50

how?is there such program in Matlab or you have written it?

Dwayne · July 7, 2013, 09:06

Quote:

Originally Posted by immortality

how?is there such program in Matlab or you have written it?

I have written it.

Dwayne · July 8, 2013, 14:00

Nobody any idea?

praveen · July 9, 2013, 22:13

At sonic point one eigenvalue becomes zero. This leads to loss of dissipation in one of the eigenfields. It manifests as an entropy violating jump inside a sonic rarefaction. Look for "entropy fix" to overcome this problem.

Dwayne · July 10, 2013, 06:57

Quote:

Originally Posted by praveen

At sonic point one eigenvalue becomes zero. This leads to loss of dissipation in one of the eigenfields. It manifests as an entropy violating jump inside a sonic rarefaction. Look for "entropy fix" to overcome this problem.

Ok i understand, thank you!

July 6, 2013, 09:34	Roe-Solver Instability at Mach=1	#1
Dwayne New Member Join Date: Jul 2013 Posts: 11 Rep Power: 13	Hey Guys, can somebody tell me the reasons why Roes solver become instable at Mach=1? The literature says that the solver becomes a central-differences-method at Mach=1 because the right side of the Roe-flux vanishes, is this true? Isn't it only one eigenvalue of the jacobian that becomes zero? Thanks for help! immortality likes this.

July 7, 2013, 05:00		#2
immortality Senior Member Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,208 Rep Power: 27	hi What software you use? __________________ Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked.

July 7, 2013, 08:50		#4
immortality Senior Member Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,208 Rep Power: 27	how?is there such program in Matlab or you have written it? __________________ Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked.

July 9, 2013, 22:13		#7
praveen Super Moderator Praveen. C Join Date: Mar 2009 Location: Bangalore Posts: 343 Blog Entries: 6 Rep Power: 18	At sonic point one eigenvalue becomes zero. This leads to loss of dissipation in one of the eigenfields. It manifests as an entropy violating jump inside a sonic rarefaction. Look for "entropy fix" to overcome this problem. Dwayne likes this. __________________ http://cpraveen.github.io http://twitter.com/cfdlab http://github.com/cpraveen

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
OpenCL linear solver for OpenFoam 1.7 (alpha) will come out very soon	qinmaple	OpenFOAM Announcements from Other Sources	4	August 10, 2012 12:00
roe solver	morteza08	Main CFD Forum	3	January 4, 2011 12:03
roe solver	morteza08	Main CFD Forum	0	May 25, 2010 16:41
ROE Riemann Solver and MUSCL	ares	Main CFD Forum	10	May 7, 2010 05:47
Who has solve RichtmyerMeshkov RM instability by using VOF solver	luckyluke	OpenFOAM Running, Solving & CFD	0	May 31, 2005 04:52

July 8, 2013, 14:00		#6
Dwayne New Member Join Date: Jul 2013 Posts: 11 Rep Power: 13	Nobody any idea?