CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Wind Turbine from SRFSimpleFoam to pimpleDymFoam (http://www.cfd-online.com/Forums/openfoam/120404-wind-turbine-srfsimplefoam-pimpledymfoam.html)

be_inspired July 6, 2013 06:44

Wind Turbine from SRFSimpleFoam to pimpleDymFoam?How?
 
1 Attachment(s)
Hi all,

I am going to reformulate the question with more information so in this way, it could be more interesting for others in future.

I would like to simulate a wind turbine using the sliding mesh aproach into OF2.1.1. My starting point is a mesh previously used for SRFSimpleFoam. It is a 120 grades sector where one blade is included and it was meshed using IcemCFD.

I suppose that the first task would be to triple the mesh, rotate each zone and merge the zones. The next step that I do not know if it is necesary to do is apply a "stitchMesh" process between each pair of patches.

I have tried to apply the stitchMesh ( without -perfect flag, because when I use it, it tells me that the distance between points is higher than a specific tolerance defined by the minor edge length) and it seems that I can do the process two times but when I want to do the last one, it tells me:


Code:

--> FOAM FATAL ERROR:
Face 23340818 reduced to less than 3 points.  Topological/cutting error B.
Old face: 2(72 73) new face: 2(72 73)

Could it be possible to impose a cylicBC to each pair of patches?

I have the impresion that because I do not have a perfect 120 grades sector and because the minimum edge length is very low due to the boundary layer refinement, I can not use the -perfect flag.

What do you think is the best way to obtain the mesh required for pimpleDymFoam solver?

Thank you very much!!

be_inspired July 24, 2013 12:04

3 Attachment(s)
I have tried to progress and now I am facing the next issue.
All the work about the mesh was done into Fluent and now a have a msh with two regions and the patches that will be the AMI patches.

I have run a simpleFoam simulation to check that all goes ok and it seems that there is no problem. When I have run pimpleDyMFoam solver, the courant number goes 5.86138e+07....

There are some guidelines http://143.248.98.6/ExCFD/dynamicMesh.pdf where they do not say anything about baffles or createAMIDict files.

Can I check where max Courant number happens in the field?
Could someone explain me why or why not it is neccesary to follow the propeller tutorial even if the mesh has been made using ICEMCFD/Fluent?

I already have defined the regions and the patches into ICEM/Fluent

Thank you very much

bennn October 8, 2013 04:25

This was a long time ago but here are a few very general guidelines.

- With unsteady AMI the courant number needs to be really low. your value will make the computation crash for sure. You need to reduce timestep size.

- Across an AMI you want the mesh size to be as constant as possible. That will make the computation crash as well, or at least create weird oscillation at the boundary.

- Your external mesh looks rather small. Expand or maybe get rid of it (you'll get rid of the AMI then which might be a good option).

- Best thing you can do is follow the propeller tutorial. The origin of the mesh has nothing to do with that.

Last but not least : With just 3 blades there's no point in running unsteady simulations (as long as the blades are not stalled). It would make sense with a pile or something. And furthermore the periodic BC will reduce your computation time a lot. I would use that in the unsteady case if possible.


All times are GMT -4. The time now is 04:35.