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

VAWT simulation with pimpleDyMFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 14, 2013, 04:07
Default VAWT simulation with pimpleDyMFoam
  #1
Member
 
Join Date: Nov 2012
Posts: 62
Rep Power: 13
Naruto is on a distinguished road
Dear all,
Hello. I have been simulating VAWT in 3-D computational domain for the last couple of months. I am using pimpleDyMFoam as my solver. My current time step is very small like 2e-05 in 3-D computational domain, when the maximum courant number is 10. I am wondering, is it normal? Because I have seen the propeller tutorial. There the time step is around 1e-04. That's why I am a little bit worried.

I have also simulated and validated my results with experimental paper in 2-D. In case of 2-D I was able to manage a time step around 1e-04. In 3-D it would be tough, as there is separation around blade tip. But still I am confused. Again I was not able to validate my results with the experimental results.

Another problem is the pressure. My pressure is fluctuating over the entire computational domain. I don't know why.
Thank you all.
__________________
Happy Foaming
Naruto is offline   Reply With Quote

Old   November 14, 2013, 05:31
Default
  #2
Senior Member
 
Jens Höpken
Join Date: Apr 2009
Location: Duisburg, Germany
Posts: 159
Rep Power: 16
jhoepken is on a distinguished road
Send a message via Skype™ to jhoepken
Have you compared your cell sizes and velocities to the ones of the propeller tutorial? These usually determine the Courantnumber, as you probably know from you CFD textbooks.

With regards to the pressure fluctuations: What are your settings in fvSolution and how large is your domain? Especially the space between your turbine and the farfield boundaries. Do you use AMI/GGI or rotate the entire mesh uniformly?

Jens
jhoepken is offline   Reply With Quote

Old   November 14, 2013, 09:08
Default
  #3
Member
 
benoit paillard
Join Date: Mar 2010
Posts: 96
Rep Power: 16
bennn is on a distinguished road
Your courant number will increase with rotor radius, rotational speed and inlet velocity, the two latters being related by TSR. Real size rotor with real inlet velocity will require very small timesteps. You might want to decrease one of these parameters, while keeping an appropriate reynolds number (you can go from 10⁷ to 10⁶ but not from 10⁶ to 10⁵ for instance). it'll give you a feasible timestep.

Fluctuating pressure is a bad sign. Check your boundary conditions, your very first timesteps, your GGI interface etc... I've had issues of fluctuating pressure with inlet mesh problems.
bennn is offline   Reply With Quote

Old   November 15, 2013, 04:33
Default
  #4
Member
 
Join Date: Nov 2012
Posts: 62
Rep Power: 13
Naruto is on a distinguished road
Sorry for my late reply. Yes I have compared my case with the propeller tutorial. It is quite similar. My turbine is small in size and the inlet velocity is high.

@Benn I am uploading my fvSchemes and fvSolutions. Would you please take a look at it?

Thanks in advance.


[ATTACH]fvSolution.txt[/ATTACH]
Attached Files
File Type: txt fvSchemes.txt (1.6 KB, 44 views)
__________________
Happy Foaming
Naruto is offline   Reply With Quote

Old   November 15, 2013, 06:29
Default
  #5
Member
 
benoit paillard
Join Date: Mar 2010
Posts: 96
Rep Power: 16
bennn is on a distinguished road
Your fvsolutions/schemes won't tell me much for that issue.

I'd need screen captures of your mesh, and your turbine characteristics.

If your inlet velocity is high then yes you'll need a very small timestep. However, you can reduce this velocity, and "extrapolate" your results to higher velocities by using coefficient of performance, and normal/tangential force coeffs. etc....
bennn is offline   Reply With Quote

Old   November 15, 2013, 06:35
Default
  #6
Member
 
Join Date: Nov 2012
Posts: 62
Rep Power: 13
Naruto is on a distinguished road
I am also attaching some snapshots of fluctuating pressure.

Capture1.JPG
Figure 1: Fluctuating pressure

Capture2.jpg
Figure 2: Meshing of a slice

Capture.JPG
Figure 3: Overall meshing

As you may see I have used structured mesh in farfield (inlet and outlet). Do you think it might be a reason?

My boundary conditions are same as propeller tutorial. I have just changed the fvSolutions and fvSchemes a little bit.
__________________
Happy Foaming
Naruto is offline   Reply With Quote

Old   November 15, 2013, 07:17
Default
  #7
Member
 
benoit paillard
Join Date: Mar 2010
Posts: 96
Rep Power: 16
bennn is on a distinguished road
What did you generate the mesh with ? is it hex or prism ? (can't really tell from your mesh) If it's prisms and cells are shared between inlet/outlet and side patches it will ruin your computation.

The pressure fluctuation looks indeed pretty nasty.
bennn is offline   Reply With Quote

Old   November 15, 2013, 08:08
Default
  #8
Member
 
Join Date: Nov 2012
Posts: 62
Rep Power: 13
Naruto is on a distinguished road
I have used pointwise to generate the mesh. It is a combination of hex and prism mesh. Sorry I did not understand the "cells are shared between inlet/outlet and side patches" part. As you may assume the prism mesh is around the blades. I am attaching some snapshots.
Capture.jpg

As you may see the prism mesh is around the blade. The other meshing elements are tetrahedra and pyramids. I have used unstructured meshing around the tip of the blade surface and structured mesh in other parts in order to save computational time. I am also attaching a snapshot.
Capture1.jpg

I have extruded the meshing on the blade surface into prism mesh. After that I used unstructured meshing only.

Do you have any advice regarding meshing? If you want I could send you my meshing file. I am afraid if I am asking for too much.
Attached Images
File Type: jpg Capture2.jpg (97.3 KB, 31 views)
__________________
Happy Foaming
Naruto is offline   Reply With Quote

Old   November 15, 2013, 08:58
Default
  #9
Member
 
benoit paillard
Join Date: Mar 2010
Posts: 96
Rep Power: 16
bennn is on a distinguished road
Post your mesh in parafoam readable format, I'll check it out.
bennn is offline   Reply With Quote

Old   November 15, 2013, 09:59
Default
  #10
Member
 
Join Date: Nov 2012
Posts: 62
Rep Power: 13
Naruto is on a distinguished road
Dear Bennn,
I have uploaded the VTK file in the following link.
https://drive.google.com/file/d/0Byl...it?usp=sharing
Please take a look if you have time.
__________________
Happy Foaming
Naruto is offline   Reply With Quote

Old   November 15, 2013, 13:31
Default
  #11
Member
 
benoit paillard
Join Date: Mar 2010
Posts: 96
Rep Power: 16
bennn is on a distinguished road
dude 1.3 gig ?
bennn is offline   Reply With Quote

Old   November 16, 2013, 05:55
Default
  #12
Member
 
Join Date: Nov 2012
Posts: 62
Rep Power: 13
Naruto is on a distinguished road
Sorry for the inconvenience and my late reply. I have uploaded another VTK format file. https://drive.google.com/file/d/0Byl...it?usp=sharing

please have a look. It is only 450MB.
__________________
Happy Foaming
Naruto is offline   Reply With Quote

Old   February 22, 2014, 11:04
Question
  #13
Senior Member
 
Himanshu Sharma
Join Date: Jul 2012
Posts: 101
Rep Power: 13
himanshu28 is on a distinguished road
Quote:
Originally Posted by Naruto View Post
Sorry for the inconvenience and my late reply. I have uploaded another VTK format file. https://drive.google.com/file/d/0Byl...it?usp=sharing

please have a look. It is only 450MB.
Hi,

I am also working on a VAWT problem being new to OpenFoam i am trying to go through the propeller tutorial in OF. Since i also intend to perform a 3D analysis can you suggest or share some reading material about defining moving body in OF.How the rotating patches are defined? also How at start a 2D case can be setup by importing geometry and meshing it.I generally use ICEM CFD for meshing.

Thank you.
himanshu28 is offline   Reply With Quote

Old   February 23, 2014, 21:37
Default
  #14
Member
 
Fengjiao Bian
Join Date: Nov 2013
Location: beijing
Posts: 30
Rep Power: 12
jiaojiao is on a distinguished road
Quote:
Originally Posted by Naruto View Post
I have used pointwise to generate the mesh. It is a combination of hex and prism mesh. Sorry I did not understand the "cells are shared between inlet/outlet and side patches" part. As you may assume the prism mesh is around the blades. I am attaching some snapshots.
Attachment 26831

As you may see the prism mesh is around the blade. The other meshing elements are tetrahedra and pyramids. I have used unstructured meshing around the tip of the blade surface and structured mesh in other parts in order to save computational time. I am also attaching a snapshot.
Attachment 26833

I have extruded the meshing on the blade surface into prism mesh. After that I used unstructured meshing only.

Do you have any advice regarding meshing? If you want I could send you my meshing file. I am afraid if I am asking for too much.
Hi! I use pointwise to generate HAWT mesh too, but not with the glf file.I am now confused about how to set the rotating zones,could you give me some advices? How do you set the rotation?with GGI or MRF or some other methods?I am not sure how to set the interface when use MRF,because when I generate mesh,the interface are disappered,I coulf only set the outer mesh boundary.So I use the topoSet/fvoption file to difine the ratating zones,but not sure is it right,plz help me.Thank you very much!
jiaojiao 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
Problem in initializing transient simulation with a finer mesh sidd CFX 8 April 29, 2016 03:25
Requiring some help regarding tornado simulation amarjogot Main CFD Forum 0 June 20, 2013 13:26
3-D VAWT simulation using pimpleDyMFOAM not working Naruto OpenFOAM Running, Solving & CFD 2 May 30, 2013 06:49
Simulation of a complex wing in solidworks flow simulation niels1900 FloEFD, FloWorks & FloTHERM 6 April 20, 2011 11:44
Continuous vs interrupted simulation sega OpenFOAM Running, Solving & CFD 4 November 3, 2008 15:29


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