0% found this document useful (0 votes)
42 views4 pages

Solucion 01 PDF

The document summarizes the derivation of the Euler equations for incompressible and compressible fluid flow in one, two, and three dimensions. For incompressible 3D flow, the Euler equations decouple into equations for conservation of mass, momentum, and energy. The document shows the derivation of these decoupled equations. For compressible 1D flow, the document expresses the Euler equations in terms of density, momentum, and total energy. It derives the Jacobian and right eigenvectors of the equations, expressing them in terms of the sound speed.

Uploaded by

librospiratas
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
42 views4 pages

Solucion 01 PDF

The document summarizes the derivation of the Euler equations for incompressible and compressible fluid flow in one, two, and three dimensions. For incompressible 3D flow, the Euler equations decouple into equations for conservation of mass, momentum, and energy. The document shows the derivation of these decoupled equations. For compressible 1D flow, the document expresses the Euler equations in terms of density, momentum, and total energy. It derives the Jacobian and right eigenvectors of the equations, expressing them in terms of the sound speed.

Uploaded by

librospiratas
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 4

CME306 / CS205B Theory Homework 8

Euler equations
For incompressible flow the inviscid 1D Euler equations decouple to:
ρt + uρx = 0
px
ut + =0
ρ
et + uex = 0
The 3D Euler equations are given by
       
ρ ρu ρv ρw
 ρu   ρu2 + p   ρuv   ρuw 
       
 ρv  +  ρuv  +  ρv 2 + p  +  ρvw  = 0 (1)
       
ρw  ρuw   ρvw   ρw2 + p 
E t (E + p)u x (E + p)v y (E + p)w z
where ρ is the density, u = (u, v, w) are the velocities, E is the total energy per unit volume and p is the
pressure. The total energy is the sum of the internal energy and the kinetic energy.
 
1 2
E = ρ e + kuk
2
= ρe + ρ(u2 + v 2 + w2 )/2
where e is the internal energy per unit mass. The assumption of incompressiblity gives
∇ · u = ux + vy + wz = 0, (2)
Show that in 3D the inviscid Euler equations with the assumption of incompressible flow decouple to:
ρt + u · ∇ρ = 0
px
ut + u · ∇u + =0
ρ
py
vt + u · ∇v + =0
ρ
pz
wt + u · ∇w + =0
ρ
et + u · ∇e = 0
The mass conservation equation takes the form:
0 = ρt + ∇ · (ρu)
= ρt + ρ∇ · u + u · ∇ρ
= ρt + u · ∇ρ = 0 .
The momentum equation along the x-axis can be condensed into
0 = (ρu)t + (ρu2 )x + (ρuv)y + (ρuw)z + px
= ρut + uρt + ρuux + u(ρu)x + ρvuy + u(ρv)y + ρwuz + u(ρw)z + px
= ρut + ρuux + ρvuy + ρwuz + px + (ρt + (ρv)y + (ρu)x + (ρw)z )
= ρut + ρu · ∇u + px + (ρt + ∇ · (ρu))
px
⇒ ut + u · ∇u + =0.
ρ

1
A similar argument reveals that the y- and z-axis momentum equations reduce to their appropriate equations,
giving (in vector form):
∇p
⇒ ut + (u · ∇)u + =0. (3)
ρ
Finally, The energy equation can be manipulated in the following way:

0 = Et + ∇ · [(E + p)u]
= Et + ∇ · (Eu) + ∇ · (pu)
= Et + E∇ · u + u · ∇E + p∇ · u + u · ∇p
     
1 1 1
= ρ e + u · u + ρt e + u · u + u · ∇ ρe + ρ u · u + u · ∇p
2 t 2 2
   
1 1
= ρet + ρu · ut + ρt e + u · u + u · ∇(ρe) + u · ∇ ρ u · u + u · ∇p
2 2
   
1 1
= ρet + ρu · ut + ρt e + u · u + ρu · ∇e + eu · ∇ρ + u · u u · ∇ρ + ρu · (u · ∇u) + u · ∇p
2 2
   
1 ∇p
= ρet + ρu · ∇e + e + u · u (ρt + u · ∇ρ) + ρu · ut + u · ∇u +
2 ρ
⇒ et + u · ∇e = 0 .

2
Compressible Flow
Find the Jacobian and the right eigenvectors for Euler’s equations in 1-D, (hint: itq
is useful, in the calculation
of the eigenvectors, to consider the enthalpy H = ρ , and the sound speed c = γp
E+p
ρ ).
   
ρ ρu
ρu +  ρu2 + p  = 0. (4)
E t Eu + pu x

You should assume the ideal gas law as your equation of state,

p(ρ, e) = (γ − 1)ρe. (5)

We begin by converting the flux term into our independent variables, x1 = ρ, x2 = ρu and x3 = E. Then
we can write the Flux term as:
  
x2 x2

 x22 2
 2
x
 x1 + (γ − 1) x3 + 12 x21  (γ − 1)x3 + 21 (3 − γ) xx2 

 = 1
(6)
2 x32
 
x3 x2 1 x2 x2 γx x2
+ 1
(1 − γ)
x 1
+ (γ − 1) x 3 + 2x x 1
3 x1 12 x 2
1

which gives our Jacobian the form:


0 1 0
 
x2
J =
 1
− 2 (3 − γ) x22 (3 − γ) xx21 (γ − 1)
1

x3 γx3 x2
−γ xx3 x2 2 + (γ − 1) x32 x1 + 32 (1 − γ) x22 γ xx21
1 1 1
 
0 1 0
=  − 12 (3 − γ)u2 (3 − γ)u (γ − 1) (7)
−γ Eu
ρ + (γ − 1)u
3
γ Eρ − 32 (γ − 1)u2 γu
q
γp
We are given the eigenvalues in lecture as λ = {u, u ± c}, where c = ρ . The first eigenvector then
simply becomes:

J~v = λ~v
⇒v2 = uv1
1
− (3 − γ)u2 v1 + (3 − γ)u2 v1 + (γ − 1)v3 = u2 v1
2
1
⇒ u 2 v1 = v3
2
1
⇒ v1 = 1 v2 = u v3 = u2
2
In order to solve the other eigenvectors, it is useful to introduce the enthalpy term ρH = E + p. Then
E+p
H=
ρ
1
= e + u2 + (γ − 1)e
2
1
= u2 + γe
2
1 E
= (1 − γ)u2 + γ .
2 ρ

3
Then we can manipulate the following to get our eigenvectors:

J~v = λ~v
⇒v2 = λv1
1
(γ − 3)u2 v1 + (3 − γ)uv2 + (γ − 1)v3 = λv2
2  
2 1 2
⇒γv3 = v3 + λ − (3 − γ)uλ − (γ − 3)u v1
2
E E 3
− γ uv1 + (γ − 1)u3 v1 + γ λv1 − (γ − 1)u2 λv1 + γuv3 = λv3
ρ ρ 2
   
1 1 3
− H + (γ − 1)u2 uv1 + (γ − 1)u3 v1 + H + (γ − 1)u2 λv1 − (γ − 1)u2 λv1 + γuv3 = λv3
2 2 2
1
− Huv1 + (γ − 1)u3 v1 + Hλv1 − (γ − 1)u2 λv1 + γuv3 = λv3
2
(λ − u)Hv1 + λ2 u − 2u2 λ + u3 v1 = (λ − u)v3


Hv1 + (λ − u)uv1 = v3
⇒v3 = H ± ucv1
⇒ v1 = 1 v2 = λ v3 = H ± uc

You might also like