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

User-defined k or epsilon?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 12, 2009, 13:40
Default User-defined k or epsilon?
  #1
Member
 
Join Date: Apr 2009
Posts: 78
Rep Power: 8
tstorm is on a distinguished road
Has anyone written successful UDF source equations for turbulence kinetic energy (k) or turbulence dissipation rate (epsilon)?

I've tried writing mine based on the k-e equations here (http://www.cfd-online.com/Wiki/Standard_k-epsilon_model) but velocity profiles in the boundary layer don't match the Blaussius curve.

Thanks!
tstorm is offline   Reply With Quote

Old   October 12, 2009, 20:52
Default
  #2
New Member
 
Javier Larrondo
Join Date: Jun 2009
Posts: 17
Blog Entries: 1
Rep Power: 8
jalarron is on a distinguished road
Currently I'm implementing a splitting technique for k-e model, so I have to add two k-e extra equations. I've done it successfully.
Are you trying to modify existing Fluent k-e equations? or you are trying to add k-e from scratch?
Cheers,
Javier
jalarron is offline   Reply With Quote

Old   October 12, 2009, 22:13
Default
  #3
Member
 
Join Date: Apr 2009
Posts: 78
Rep Power: 8
tstorm is on a distinguished road
I'm trying to start from scratch. I'm making several v2-f models, which use the standard k and e equations.
Would you mind posting your k and e source code? I'd appreciate it.
Thanks!
tstorm is offline   Reply With Quote

Old   October 12, 2009, 23:47
Default
  #4
New Member
 
Javier Larrondo
Join Date: Jun 2009
Posts: 17
Blog Entries: 1
Rep Power: 8
jalarron is on a distinguished road
Which viscous model are you using for your implementation from scratch of k-e? laminar model + 2 UDS equation for k-e?
About the source code, I have to check first with my supervisor, just for IP reasons. I think it shouldn't be a problem, because is a method already published and I'm just implementing it in FLUENT.
Currently I haven't tried to start from scratch but it is one option that we were thinking because patching the in-built k-e it's a bit tricky specially in the near-wall region.
Cheers,
Javier
jalarron is offline   Reply With Quote

Old   October 13, 2009, 01:11
Default
  #5
Member
 
Join Date: Apr 2009
Posts: 78
Rep Power: 8
tstorm is on a distinguished road
I'm actually using the k-e model with turbulence equations turned off. I hadn't considered running it with the laminar model. I might try that tomorrow, it may solve some problems.

Just to give you some background, the v2-f model is based on the k-e model, and uses the same k and e equations, with some minor differences. The constants are calibrated differently - in particular, Ce2 is a function of turbulent quantities, and isn't constant. Also, the model is integratable to solid walls and doesn't use damping functions or wall functions. So, FLUENT's k and e equations won't work.

I hope you can share your code - I can't imagine that the k and e sources would be proprietary, but I understand if you can't share them.

Thanks!
tstorm is offline   Reply With Quote

Old   October 13, 2009, 11:35
Default
  #6
Member
 
Join Date: Apr 2009
Posts: 78
Rep Power: 8
tstorm is on a distinguished road
I just remembered why I haven't used the laminar model with 4 UDSs (k, e, v2, f). I need to use a user-defined turbulent viscosity, which is made easy when you use the k-e model with turbulence equations turned off.

I'm assuming if you plan to start from scratch with k and e UDSs you'll need to store turbulent viscosity as well. Had you figured out how you were going to hook that to the solver?
tstorm is offline   Reply With Quote

Old   October 13, 2009, 21:21
Default
  #7
New Member
 
Javier Larrondo
Join Date: Jun 2009
Posts: 17
Blog Entries: 1
Rep Power: 8
jalarron is on a distinguished road
I haven't tried anything with the laminar option, but I think that you could store the turbulent viscosity in a UDM variable, and updated in every iteration.
I'm sending you a document of what I did through a private message. OK?
Cheers,
Javier
jalarron is offline   Reply With Quote

Old   October 14, 2009, 11:55
Default
  #8
Member
 
Join Date: Apr 2009
Posts: 78
Rep Power: 8
tstorm is on a distinguished road
You could store mu_t as a UDM, but I don't see how you'd hook it to the solver in the laminar model. The laminar model should have mu_t = 0. If you figure that out I'd love to hear how it turns out.
tstorm is offline   Reply With Quote

Reply

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
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 13 November 4, 2013 15:13
user defined turbulence model manuutin STAR-CD 5 October 14, 2009 05:29
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 1 November 25, 2008 21:21
Gradient of a User defined Variable Ramadas CFX 2 August 21, 2007 09:19
User Defined Scalars - Returning Values Carlos V. FLUENT 0 April 19, 2006 18:18


All times are GMT -4. The time now is 22:57.