PMatrix in Numerics Structure

atelcikti1 · November 6, 2020, 06:13

Hi people,

I am trying to understand SU2 jacobian matrix structure and components in depth. In numerics structure cpp file there is a Pmatrix and its inverse function. Why it is implemented? I checked where those functions called and see used in implicit application of jacobian matrix. I am aware that those are related to eigen-values and vectors of jacobian as in stated and derived in Hirsh book however what is the main reason to calculate? Any reference beside Hirsh will be appreciated.

Have a nice day.

-Ali

aerosjc · May 29, 2022, 04:01

You have to use the diagonalized form instead of the Roe matrix. It may looks like that the Roe matrix A = P lambda P^{-1}. However, it is not in fact. Because you have the entropy correction which modifies the engenvalues lambda.

November 6, 2020, 06:13	PMatrix in Numerics Structure	#1
atelcikti1 Member Join Date: Mar 2017 Posts: 61 Rep Power: 9	Hi people, I am trying to understand SU2 jacobian matrix structure and components in depth. In numerics structure cpp file there is a Pmatrix and its inverse function. Why it is implemented? I checked where those functions called and see used in implicit application of jacobian matrix. I am aware that those are related to eigen-values and vectors of jacobian as in stated and derived in Hirsh book however what is the main reason to calculate? Any reference beside Hirsh will be appreciated. Have a nice day. -Ali

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Lattice Structure Wind Load – STAR-CCM+	Fluidy	STAR-CCM+	1	May 17, 2020 14:33
Number of time step in Transient structure for Fluid Structure interaction	YANNAPOL	Structural Mechanics	1	July 5, 2015 09:41
New OpenFOAM Forum Structure	jola	OpenFOAM	2	October 19, 2011 07:55
New ANSYS forum structure, what do you think?	Peter	FLUENT	10	February 6, 2009 05:07
New ANSYS forum structure, what do you think?	Peter	CFX	5	February 4, 2009 12:59

May 29, 2022, 04:01		#2
aerosjc Member Jingchang.Shi Join Date: Aug 2012 Location: Hang Zhou, China Posts: 78 Rep Power: 14	You have to use the diagonalized form instead of the Roe matrix. It may looks like that the Roe matrix A = P lambda P^{-1}. However, it is not in fact. Because you have the entropy correction which modifies the engenvalues lambda.