Jacobian matrix and Navier Stokes equation

Click For Summary
SUMMARY

The discussion centers on the application of the Jacobian matrix in the linearization of the Navier-Stokes equations using the Beam-Warming method. A participant struggles to derive specific terms in the A and B matrices, particularly when differentiating the energy vector E with respect to density and velocity components. The conversation also touches on the relationship between pressure, energy, and density, emphasizing the need to express E in terms of the components of U. Additionally, the discussion introduces interpolation schemes, specifically the QUICK scheme and the Van Leer/Van Albada schemes, and questions the consistency of slopes in these methods.

PREREQUISITES
  • Understanding of the Navier-Stokes equations
  • Familiarity with the Jacobian matrix in fluid dynamics
  • Knowledge of the Beam-Warming method for linearization
  • Basic concepts of interpolation schemes in numerical methods
NEXT STEPS
  • Study the derivation of the Jacobian matrix in the context of the Navier-Stokes equations
  • Learn about the Beam-Warming method and its applications in computational fluid dynamics
  • Investigate the QUICK scheme and its properties in numerical interpolation
  • Explore the Van Leer and Van Albada schemes to understand their consistency in slope calculations
USEFUL FOR

Researchers, graduate students, and professionals in fluid dynamics, computational physics, and numerical analysis who are working with the Navier-Stokes equations and interpolation methods.

mertcan
Messages
343
Reaction score
6
upload_2018-2-15_22-52-31.png
upload_2018-2-15_22-57-3.png


Hi, in first attachment/picture you can see the generalized navier stokes equation in general form. In order to linearize these equation we use Beam Warming method and for the linearization process we deploy JACOBİAN MATRİX as in the second attachment/picture. But on my own I can ONLY obtain the third row in A matrix and second row in B matrix. I can not handle the OTHER terms, by the way GAMMA is the ratio of specific heats as it is written in picture 2.
Could you help me about my situation?

ıf it helps, pictures has been cut off that link
<< Mentor Note -- link to copyrighted textbook deleted >>
 

Attachments

  • upload_2018-2-15_22-52-31.png
    upload_2018-2-15_22-52-31.png
    18 KB · Views: 7,490
  • upload_2018-2-15_22-57-3.png
    upload_2018-2-15_22-57-3.png
    37.7 KB · Views: 5,549
Last edited by a moderator:
Engineering news on Phys.org
to be more explicit : when I take the derivative of E vector with respect to density (first element of U vector) I can obtain almost the first column of A matrix in attachment 2 BUT when I take the derivative with respect to pu( density*velocity on x axis, second element of U) or pv or pE(density*energy) I can not obtain the other elements of A matrix entirely. For instance if I take the derivative of pu(density*velocity on x axis, first element of E vector) with respect to pv(density*velocity on y axis, third element of U vector) it should NOT BE 0 because there is a common variable between pu and pv and it is the density variable or if I take the derivative of pu(density*velocity on x axis, first element of E vector) with respect to p(density, first element of U vector) it should BE u NOT zero. Could you help me about that??

Also there is a pressure term which equals (gamma-1)*density*(E-0.5*(u^2+v^2)
 
pleas use latex to help format equations:

https://www.physicsforums.com/help/latexhelp/

It helps to introduce new variables for the components of U and rewrite E in terms of U:
##U= (u_1, u_2, u_3,u_4)##
Now rewrite E in terms of these U and note that the energy and pressure are functions of U.
As you said, the pressure can be rewritten using the equation of state as: ##p=(\gamma-1)\rho(E-\frac{1}{2}(u^2+v^2))## and you need to substitute it into the vector E before determining the derivatives.
Then note that the total energy should be rewritten as the sum of the internal and kinetic energy to make clear the dependence on ##u_1..u_4##.

Can you rewrite the E-vector in terms of the components ##u_1..u_4##?
 
bigfooted said:
pleas use latex to help format equations:

https://www.physicsforums.com/help/latexhelp/

It helps to introduce new variables for the components of U and rewrite E in terms of U:
##U= (u_1, u_2, u_3,u_4)##
Now rewrite E in terms of these U and note that the energy and pressure are functions of U.
As you said, the pressure can be rewritten using the equation of state as: ##p=(\gamma-1)\rho(E-\frac{1}{2}(u^2+v^2))## and you need to substitute it into the vector E before determining the derivatives.
Then note that the total energy should be rewritten as the sum of the internal and kinetic energy to make clear the dependence on ##u_1..u_4##.

Can you rewrite the E-vector in terms of the components ##u_1..u_4##?
Initially, I got the logic @bigfooted thank you for valuable return. By the way I do not want to create another thread so I would like to ask DIFFERENT question related to interpolation schemes here.
First of all I am aware of the fact that QUICK SCHEME has consistent slope (for instance at the left side of node 3 in my attachment same slopes exist) as you can see in my picture/attachment. But I must express that I can prove slopes at the left side of node(like in picture) are equal in QUICK SCHEME thus it is consistent but I know there are another schemes like VAN LEER VAN ALBADA SCHEME which are non linear and I can NOT prove how those SCHEMES may be consistent in terms of slopes at the left side of node like QUICK SCHEME. At the centre of length (length between node 2 and node 3 in my attachment) which means left side of node 3 QUICK SCHEME always ensure the consistency of slope and I can prove but HOW DO WE KNOW THAT VAN LEER VAN ALBADA SCHEMES MAY ENSURE THE CONSISTENCY OF slope at the left side of node 3?? How can we PROVE it??
 

Attachments

  • picture:attachment.png
    picture:attachment.png
    32.5 KB · Views: 710
Last edited:
I recommend that you start a new post on this. can you be a bit more precise? Which van Leer/Albada scheme do you mean? Could you please show what you tried to do?
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
3K
Replies
9
Views
3K
  • · Replies 8 ·
Replies
8
Views
3K
  • · Replies 10 ·
Replies
10
Views
3K
Replies
3
Views
2K
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 8 ·
Replies
8
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
5K