CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

Wind turbine simulation issues

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 19, 2017, 18:52
Default Wind turbine simulation issues
  #1
New Member
 
Carlos Muñoz
Join Date: Apr 2017
Posts: 5
Rep Power: 6
Carlos_lboro is on a distinguished road
Hello everybody, this is my first post in this forum so I'll try to be concise:

I'm simulating the near wake of a wind turbine by designing an actuator disc, in the simulation I have made the actuator disc as a subdomain porous medium with a momentum loss.
https://goo.gl/photos/1MSUACnsA8D1vZCt9
https://goo.gl/photos/e9bsYd3Si4SH4aBp7
That model worked perfectly fine, no issues there, but then I tried to make it a bit more complex by creating several concentric discs (or rings) each of those being as well a porous subdomain and with different loss coefficients because the energy taken form the wind is different in different parts of the wing.
I found out that the program was not able to calculate de solution, throwing at me an error 2.
Pictures of the second model:
https://goo.gl/photos/TLbbWVTYMZVSeXxTA
https://goo.gl/photos/szG4kNFt66aW75B16
https://goo.gl/photos/3dpiw3YEVVB5JQPe8

Thank you very much for any help or hints that you can give me.
Carlos_lboro is offline   Reply With Quote

Old   April 19, 2017, 19:08
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,877
Rep Power: 132
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Your images are not attached.

But I do not recommend you use a porous model to model a turbine. I have never seen a wind turbine which looked like a sponge. You should use a general momentum source term for this, then you will have much more control over the source term.

And finally: If you are trying to model different amounts of momentum being applied to different parts of the wind turbine disk then I think you are pushing the source term approach too far. If you want to see the effects of different loading across the turbine disk then I recommend you model the wind turbine properly and get the full, complete 3D flow field.
ghorrocks is offline   Reply With Quote

Old   April 19, 2017, 20:03
Default
  #3
New Member
 
Carlos Muñoz
Join Date: Apr 2017
Posts: 5
Rep Power: 6
Carlos_lboro is on a distinguished road
momenumloss.jpg

Is this what you mean by general momentum source?

Thank you!
Carlos_lboro is offline   Reply With Quote

Old   April 19, 2017, 20:40
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,877
Rep Power: 132
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Yes, but do not use a loss model and a general momentum source at the same time. Remove the loss model.

Also note you need to correctly set the source term coefficient to get convergence - read the documentation on how to configure that.
ghorrocks is offline   Reply With Quote

Old   April 19, 2017, 21:58
Default
  #5
New Member
 
Carlos Muñoz
Join Date: Apr 2017
Posts: 5
Rep Power: 6
Carlos_lboro is on a distinguished road
I have read this:

https://www.sharcnet.ca/Software/Ans...fxBasiMomeIsot

but it doesn't really clarify anything about convergence or correct values, it randomly throws a 10^5 coefficient that I cannot understand.

In this other one they talk about the convergence:

https://www.sharcnet.ca/Software/Ans.../i1300063.html

But I still can't understand the partial derivation of Sm vs U...

I'm sorry, I'm sure this must sound super easy for you but for me is a lot of fluids theory in a really short amout of time and with a very brief explanation by anyone.
Carlos_lboro is offline   Reply With Quote

Old   April 20, 2017, 21:55
Default
  #6
New Member
 
Carlos Muñoz
Join Date: Apr 2017
Posts: 5
Rep Power: 6
Carlos_lboro is on a distinguished road
Doing what you said of avoiding the loss coefficient and just using the general momentum loss does not "take" any energy from the flow, therefore the velocity at both sides of my "actuator disc" are the same. Any ideas of what could have gone wrong?
Carlos_lboro is offline   Reply With Quote

Old   April 20, 2017, 22:53
Default
  #7
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,877
Rep Power: 132
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
You will not get a velocity increase from one side of the actuator to the other as this does not conserve mass. Rather you get the flow accelerates into the actuator, passes through the actuator at a high velocity and then dissipates into the surrounding lower velocity fluid.

The source term coefficient is just there for convergence. If you are converging OK then you might not need it. Note the page you quote explains source term coefficients, but in section 1.3.2.1.1.

Have a look at the units of the source term and the units of the momentum your turbine removes from the flow. That might help you understand how they are linked together.
ghorrocks is offline   Reply With Quote

Old   April 23, 2017, 10:04
Default
  #8
New Member
 
Carlos Muñoz
Join Date: Apr 2017
Posts: 5
Rep Power: 6
Carlos_lboro is on a distinguished road
Hi Glenn, thank you for your answer, it wasn't showing to me.

I want an actuator disc as a sink of momentum, I need it to drain momentum from the flow thats going through it and therefore change the pressure (and velocity).

I understood from the documentation and other forum posts that I have been looking at that the source term is only for convergence and that usually is a large number, in several posts I've seen used 10000-100000.

The units I get on ANSYS of the momentum source are [kg m^-2 s^-2] and the ones from the turbine are [kg m s^-1] or that is what I've found... But still dont get how they relate to each other.

I need to eventually get the thrust so I would need to calculate:
F = (m dot * V)e - (m dot * V)0 + (pe - p0) * Ae
being e the exit and 0 the undisturbed flow before the actuator.

Again thank you very much for your help, my apologies for being so slow learning this kind of things.
Carlos_lboro is offline   Reply With Quote

Old   April 23, 2017, 20:17
Default
  #9
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,877
Rep Power: 132
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Quote:
Hi Glenn, thank you for your answer, it wasn't showing to me.
Can you PM me with what happened? I would like to know why my post was not showing.


Quote:
I understood from the documentation and other forum posts that I have been looking at that the source term is only for convergence and that usually is a large number, in several posts I've seen used 10000-100000.
You are getting confused with the source term coefficient, but you mean the source term. I recommend for now you do not worry about the source term coefficient. Get the source term right first and see if you need a source term coefficient. Not all cases need a source term coefficient.

Quote:
The units I get on ANSYS of the momentum source are [kg m^-2 s^-2] and the ones from the turbine are [kg m s^-1] or that is what I've found... But still dont get how they relate to each other.
If you divide the momentum source by the turbine figure you quote you get [m^-3 s]. That is the inverse of the volume flow rate. Hopefully that is a clue for you.

Quote:
F = (m dot * V)e - (m dot * V)0 + (pe - p0) * Ae
You can simplify this. m(dot) and V is the same at the inlet and exit plane as mass is conserved. This means it just comes down the pressure difference.
ghorrocks is offline   Reply With Quote

Old   September 14, 2021, 01:24
Cool Actuator disc solver
  #10
New Member
 
Edgar Alejandro Martínez Ojeda
Join Date: Jul 2019
Posts: 20
Rep Power: 4
Edgar Alejandro Martínez is on a distinguished road
Look at my YT video in which I generalized Erik Svenning's actuator disc. You can calculate the pressure drop in ParaView and the solver gives you the thrust.
https://www.youtube.com/watch?v=MbqbgGLJRso

Link to the solver: https://github.com/EdgarAMO/multipleDiskSimpleFoam

Link to the case file: https://github.com/EdgarAMO/actuator-disk-farm
Edgar Alejandro Martínez is offline   Reply With Quote

Reply

Tags
cfx, moment loss, porous, turbine

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
wind turbine blade flow simulation rsskarthikeyan FLUENT 0 May 27, 2015 06:40
Torque of wind turbine simulation caohan FLUENT 8 August 12, 2014 00:01
Simulation of fan Vs wind turbine mohammad Main CFD Forum 0 November 5, 2013 09:43
wind turbine simulation, Ansys Post results question Laions CFX 7 September 20, 2011 06:13
3D simulation of wind turbine in Yaw wind(in a lateral wind) mohammad Main CFD Forum 0 December 28, 2010 04:26


All times are GMT -4. The time now is 03:07.