CFD Online Logo CFD Online URL
Home > Forums > General Forums > Main CFD Forum

Maths problem when deriving conservation of energy equations.

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   June 15, 2009, 11:08
Default Maths problem when deriving conservation of energy equations.
New Member
Join Date: Jun 2009
Posts: 4
Rep Power: 16
Demonwolf is on a distinguished road
Ive been reading John Andersons book 'Computational Fluid Dynamics the basics with applications'. Im upto the part where the energy equation is being derived.
The energy equation is derived in nonconservation form and is in terms of the total energy e+V^{2}/2. This is on page 70.
It then goes on to describe other ways that the equation can be written in different energy terms and begins to do some manipulations to rewrite the energy equations in terms of internal energy only.

It starts by multiplying the momentum equations in nonconversation form by u, v, and w respectively. The original momentum equations are given below:

\rho\frac{Du}{Dt} = -\frac{\partial p}{\partial x}+\frac{\partial \tau_{xx}}{\partial x}+\frac{\partial \tau_{yx}}{\partial y}+\frac{\partial \tau_{zx}}{\partial z}+\rhof_x -- (2.50a)

\rho\frac{Dv}{Dt} = -\frac{\partial p}{\partial y}+\frac{\partial \tau_{xy}}{\partial x}+\frac{\partial \tau_{yy}}{\partial y}+\frac{\partial \tau_{zy}}{\partial z}+\rhof_y -- (2.50b)

\rho\frac{Dw}{Dt} = -\frac{\partial p}{\partial z}+\frac{\partial \tau_{xz}}{\partial x}+\frac{\partial \tau_{yz}}{\partial y}+\frac{\partial \tau_{zz}}{\partial z}+\rhof_z -- (2.50c)

The first part of the manipulation involves multiplying the equations 2.50a, 2.50b, and 2.50c by u, v and w respectively. The book gives the results as shown below:

\rho\frac{D(u^{2}/2)}{Dt} = -u\frac{\partial p}{\partial x}+u\frac{\partial \tau_{xx}}{\partial x}+u\frac{\partial \tau_{yx}}{\partial y}+u\frac{\partial \tau_{zx}}{\partial z}+\rhouf_x -- (2.67)

\rho\frac{D(v^{2}/2)}{Dt} = -v\frac{\partial p}{\partial y}+v\frac{\partial \tau_{xy}}{\partial x}+v\frac{\partial \tau_{yy}}{\partial y}+v\frac{\partial \tau_{zy}}{\partial z}+\rhovf_y -- (2.68)

\rho\frac{D(w^{2}/2)}{Dt} = -w\frac{\partial p}{\partial z}+w\frac{\partial \tau_{xz}}{\partial x}+w\frac{\partial \tau_{yz}}{\partial y}+w\frac{\partial \tau_{zz}}{\partial z}+\rhowf_z -- (2.69)

I can understand the right hand sides of the above three equations. But I dont understand how to get the left hand sides. I am confused how the following is equal:

\rho\frac{Du}{Dt}\times u = \rho\frac{D(u^{2}/2)}{Dt}

Can someone show me how this is done, and also tell me what area of differential mathematics I need to read to plug the gap in my knowledge. Thank you.
Demonwolf is offline   Reply With Quote

Old   June 15, 2009, 14:21
New Member
Join Date: Jun 2009
Posts: 4
Rep Power: 16
Demonwolf is on a distinguished road
Ive just worked out the solution to this myself, which Ill give below.. so no need to post. It was pretty simple in the end and Im surprised I didnt get it sooner.

Use substitution and then the chain rule to show how \frac{Du}{Dt}\times u = \frac{D(u^2/2)}{Dt}

let y = u^2/2 Therefore
\frac{Dy}{Du} = u

Then the chain rule states :

\frac{Dy}{Dt} =

\Rightarrow \frac{Dy}{Dt} = u\times\frac{Du}{Dt}

Demonwolf is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
turbulent energy spectrum problem cfd.newbie Main CFD Forum 6 September 24, 2015 16:13
SIMPLE and energy equation convergence Fabio Main CFD Forum 0 June 1, 2007 06:06
finite volume daniel Main CFD Forum 17 April 11, 2005 20:36
Energy equation convergence problem Reza Main CFD Forum 0 August 27, 2003 13:09
extremely simple problem... can you solve it properly? Mikhail Main CFD Forum 40 September 9, 1999 09:11

All times are GMT -4. The time now is 13:15.