CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

How to compute the viscous dissipation term

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 17, 2013, 03:14
Question How to compute the viscous dissipation term
  #1
Member
 
yehanyu
Join Date: Mar 2012
Location: Beijing, China
Posts: 36
Rep Power: 5
yhy20081016 is on a distinguished road
Hi, everybody. I am performing 2D unsteady CFD simulations on a cross-flow fan. I am using ANSYS FLUENT 12.1. I adopt the standard k-omega model to simulate turbulence.
I want to compute the viscous dissipation term. I cannot find this variable in FLUENT. Hence I exported the solution to Tecplot, and then write some equations to compute it.
I use the following equation to compute the viscous dissipation term: (This equation can be found in ANSYS FLUENT 12.0 Theory Guide, page 4-24 and 4-25)
viscous dissipation=\frac{\partial}{\partial x_{j}}(u_{i}(\tau_{ij})_{eff})
where
(\tau_{ij})_{eff}=\mu_{eff}(\frac{\partial u_{i}}{\partial x_{j}}+\frac{\partial u_{j}}{\partial x_{i}})-\frac{2}{3}\mu_{eff}\frac{\partial u_{k}}{\partial x_{k}}\delta_{ij}

I use [Data]|[Alter]|[Specify Equations…] in Tecplot to implement these equations. The expressions which I type into Tecplot are listed as follows:
{duxdx}=ddx({x-velocity})
{duxdy}=ddy({x-velocity})
{duydx}=ddx({y-velocity})
{duydy}=ddy({y-velocity})
{viscous_heating2}=ddx({x-velocity}*({viscosity-eff}*({duxdx}*2)-2.0/3.0*{viscosity-eff}*({duxdx}+{duydy}))+{y-velocity}*({viscosity-eff}*({duxdy}+{duydx}))) + ddy({x-velocity}*({viscosity-eff}*({duydx}+{duxdy}))+{y-velocity}*({viscosity-eff}*({duydy}*2)-2.0/3.0*{viscosity-eff}*({duxdx}+{duydy})))
Finally, I display contour of the variable “viscous_heating2”. However, to my surprise, the minimum value of this variable is -679152, and the maximum value is 352578. This is weird – because the dissipation should be always positive. How could negative values occur?
yhy20081016 is offline   Reply With Quote

Old   December 17, 2013, 03:33
Default
  #2
Senior Member
 
cfdnewbie
Join Date: Mar 2010
Posts: 551
Rep Power: 11
cfdnewbie is on a distinguished road
One possible cause of your troubles is that you are using the ddx functions in tecplot to compute derivatives. That's among the worst methods you can use, they are first oder accurate (at most second), and depending on your grid, that can be way to inaccurate.
Can you compute your derivatives in your code? On the original mesh maybe? Then there's at least the chance that you are not hunting tecplot ghosts with this....
cfdnewbie is offline   Reply With Quote

Old   December 17, 2013, 23:26
Default
  #3
Member
 
yehanyu
Join Date: Mar 2012
Location: Beijing, China
Posts: 36
Rep Power: 5
yhy20081016 is on a distinguished road
Thank you very much. I'll try it.
yhy20081016 is offline   Reply With Quote

Reply

Tags
dissipation, tecplot, turbulence, viscous

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Viscous term expansion santiagomarquezd OpenFOAM Programming & Development 1 December 13, 2011 21:44
Viscous work term query siw CFX 1 December 6, 2010 07:41
Adding viscous dissipation term to energy equation newbee OpenFOAM Running, Solving & CFD 4 August 25, 2006 14:33
VISCOUS heat dissipation Joseph CFX 0 October 1, 2004 08:52
Is the viscous dissipation term fully included? yann Phoenics 0 June 5, 2003 10:23


All times are GMT -4. The time now is 18:20.