|
[Sponsors] |
[solids4Foam] FSI simulation of a Wind Turbine |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
Dinmukhamed
Join Date: Nov 2020
Posts: 3
Rep Power: 2 ![]() |
Dear all,
I am trying to perform an FSI simulation of a Wind Turbine, and, during my search for appropriate tool, I found solids4foam by Dr. Cardiff. I've build my case based on tutorials, which were included with the solids4foam, and I faced some difficulties. Here are some information about my setup: 1) The solids4foam was built on foam-extend 4.1 2) I am using full rotor model of the NREL 5MW wind turbine Further, I will explain the difficulties, which I faced: 1) Aligning fluid and solid models. In all solids4foam tutorials, the models were created using blockMesh and they are perfectly align with each other. However, I am studying blades, which have complex structure and they are much harder to align perfectly. I've used ANSYS fluent and snappyHexMesh to create a mesh, but using both methods gave me the following error: Code:
--> FOAM FATAL ERROR: ZoneA and zoneB interfaces are not conformal (zoneA patch = interface, zoneB patch = interface) directMap method requires conformal interfaces! From function directMap::checkZoneSizes() const in file fluidSolidInterfaces/fluidSolidInterface/interfaceToInterfaceMapping/directMapInterfaceToInterfaceMapping/directMapInterfaceToInterfaceMapping.C at line 65 RBF interpolation: Code:
terminate called after throwing an instance of 'std::bad_alloc' what(): std::bad_alloc Aborted Code:
--> FOAM FATAL ERROR: Master point addressing is not correct From function GGIInterpolation::masterToSlavePointInterpolate(const Field<Type> pf) in file /home/dimash/foam/foam-extend-4.1/src/foam/lnInclude/GGIInterpolation.C at line 492. FOAM aborting Aborted I attached the figure with fluid model, solid model, and an overlap of the two models: model.PNG 2) Rotating and deforming at the same time. As I understand, solids4foam uses dynamicMotionSolverFvMesh to deform the mesh, but to simulate rotation effect I should use different dynamicFvMesh parameter. Is it possible to simulate a rotation effect with two way FSI without changing the source code, or I need to modify it? Thank you in advance Best regards, Dinmukhamed |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 884
Rep Power: 28 ![]() ![]() |
Hi Dinmukhamed,
This looks like an interesting problem. I have a few comments:
Philip |
|
![]() |
![]() |
![]() |
![]() |
#3 |
New Member
Dinmukhamed
Join Date: Nov 2020
Posts: 3
Rep Power: 2 ![]() |
Thank you for your reply,
I've tried to do what you suggested: 1) I've tried to build the solids4foam on foam-extend 4.0 and both com and org versions of OpenFoam, but I did not succeed in it. Is foam-extend 4.0 better than 4.1 or I can continue to work on the current version? 2) As you suggested, now I am testing the case on a tip of the blade. I am also turning coupling after some simulation time. 3) The GGI interpolation works for me, but now I have a floating point error. I will add the log file to this post. When I set the wind speed to 0, this error disappears. As I understand, this problem may be caused by inaccurate calculation of forces, acting on a solid model. Can you suggest how to fix it? I will try to enhance both fluid and solid meshes and see will it help or not. The log file: log.solids4Foam.txt Best regards, Dinmukhamed |
|
![]() |
![]() |
![]() |
![]() |
#4 | |||
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 884
Rep Power: 28 ![]() ![]() |
Quote:
Quote:
Quote:
Philip |
||||
![]() |
![]() |
![]() |
![]() |
#5 |
New Member
Dinmukhamed
Join Date: Nov 2020
Posts: 3
Rep Power: 2 ![]() |
After performing multiple simulations and testing different models, I successful run an FSI simulation of a wind turbine, but I am still struggling with rotational part.
The progress in FSI: 1) The main problem was poor quality of the solid mesh. I used snappyHexMesh to create better solid mesh and now everything works fine with both linear and non-linear deformation.I am still struggling with implementing rotation to the FSI simulation. Can you suggest what I should do? At this moment I found following: The deformation is controlled by dynamicFvMesh dynamicMotionSolverFvMesh from the dynamicMeshDict file. The rotation is controlled by dynamicFvMesh turboFvMesh from the same file. In order to have both deformation from FSI and rotation effect, I should combine both of the dynamic solvers.Should I change the source code and to implement my own dynamic solver or I can somehow rotate the cellZone with FSI deformation with existing tools? Thank you in advance Update (20.12.2020): I've tried to a 2D simulation with fixedRotation boundary condition at solid part, however my solid model is not rotating. I will continue to search for possible solutions. Sincerely, Dinmukhamed Last edited by Dimash; December 20, 2020 at 09:30. |
|
![]() |
![]() |
![]() |
![]() |
#6 |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 884
Rep Power: 28 ![]() ![]() |
Hi Dinmukhamed,
Glad to hear that you have made progress. Thanks for your sharing your experience and giving tips. As regards rotation, can you give more details? Do you mean that you would like to apply a rotation boundary condition to the root of the blade? If so, then you need to apply a rotation boundary condition ("fixedRotation" is an option) to the solid at the blade root (the solid does not have a mesh motion solver as it uses a Lagrangian approach), and then you need to make sure the fluid mesh motion solver is capable of letting the fluid mesh be "pushed" around or even to specify the rotation of part of it. I suggest creating a simple representative test case, probably 2-D to check it works. Philip |
|
![]() |
![]() |
![]() |
Thread Tools | Search this Thread |
Display Modes | |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
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 problem | sanchovg2 | STAR-CCM+ | 12 | March 5, 2013 05:18 |
[ANSYS Meshing] 3d wind turbine mesh for multiphase simulation | mingersai | ANSYS Meshing & Geometry | 0 | January 17, 2012 19:20 |
3D simulation of wind turbine in Yaw wind(in a lateral wind) | mohammad | Main CFD Forum | 0 | December 28, 2010 04:26 |