CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Why sovle e p U but not T p U in sonicFoam?

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 16, 2010, 10:43
Default Why sovle e p U but not T p U in sonicFoam?
  #1
New Member
 
yu
Join Date: Nov 2009
Posts: 26
Rep Power: 16
universez is on a distinguished road
Hi, I am just curious why variable solved in sonicFoam are e p U. Then, calcualte T using function

inline scalar T (scalar f, scalar T0, scalar (specieThermo::*F)(const scalar) const, scalar (specieThermo::*dFdT)(const scalar) const ) const

of class specieThermo.

What is the benifit comparing to solve T p U.


Thanks, yu.
universez is offline   Reply With Quote

Old   July 19, 2010, 09:11
Default
  #2
Senior Member
 
Stefan Herbert
Join Date: Dec 2009
Location: Darmstadt, Germany
Posts: 129
Rep Power: 17
herbert is on a distinguished road
Hi Yu,

I think it is easier to include temperature-dependent thermodynamic properties (especially the heat capacity) when solving for e instead of T.

Regards,
Stefan
herbert is offline   Reply With Quote

Old   July 19, 2010, 09:28
Default
  #3
New Member
 
yu
Join Date: Nov 2009
Posts: 26
Rep Power: 16
universez is on a distinguished road
Stefan,

Thank you for your reply.

I agree. It will make a implict couple between heat capacity and temperature.

Is that possible that solving e will give more smoother result than solving T? So, solving e will be more numerically stable>

Thanks

Yu

Quote:
Originally Posted by herbert View Post
Hi Yu,

I think it is easier to include temperature-dependent thermodynamic properties (especially the heat capacity) when solving for e instead of T.

Regards,
Stefan
universez is offline   Reply With Quote

Old   July 19, 2010, 17:25
Default
  #4
New Member
 
yu
Join Date: Nov 2009
Posts: 26
Rep Power: 16
universez is on a distinguished road
I guess that there are several difficulties to solve T instead of e.

1, solving e is to solve Cv*T, but involves only one variable, which makes code easier.

2, in openfoam, solving Cv*T is not straight forward. I tried different ways to put Cv in the equation, this one works,

{
solve
(
Cv*fvm::ddt(rho, T)
+ Cv*fvm::div(phi, T)
- Cv*fvm::laplacian(turbulence->alphaEff(), T)
==
- p*fvc::div(phi/fvc::interpolate(rho))
);
}

However, it separate Cv and T, and make the scheme less accurate.
universez is offline   Reply With Quote

Reply

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
energy in sonicFoam joern OpenFOAM Running, Solving & CFD 1 September 24, 2019 01:15
Crash when using sonicFoam Horus OpenFOAM 1 June 16, 2010 13:57
sonicFoam - Total Temperature not constant cboss OpenFOAM Running, Solving & CFD 0 April 21, 2010 12:11
Using Compressible Solver (sonicFoam) to solve subsonic flows ezsoal OpenFOAM 0 October 27, 2009 10:13
SonicFoam forwardStepTutorial doesnbt complete the run alberto OpenFOAM Bugs 1 June 10, 2007 16:35


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