Bird Wing simulation suggestions
Hi!!!
As stated in the title, for my third year undergraduate project I would like to simulate the flight of an hummingbird. The complexity involved is greater than the one required for an undergrad thesis, therefore, my supervisor recommended me to concentrate on the wing alone. He suggested to reduce it to a small, thin ellipse and assign to it the special rotation (infinity sign alike?!). He said OpenFoam is a powerlful tool so I am getting used to it over summer in order to be ready for october. I went through the user guide provided and installed the 1.6ext where there is a rotating motion solver which may fit my case. How would you proceed? Is there a specific solver/tool which may help me in simulating the wing and extract useful data such as the variation of Lift and Drag with the movement etc? What solver or utility do you reckon is required for my case? How will I be able to produce the mesh and give to it the rotation? The part is currently a Solidworks CAD file. Any suggestion is more than welcome!! Thanks! :D |
Hi, welcome to the forum,
I'm not into the topic myself, but there was someone asking similar questions just a few weeks ago concerning the flight of mosquitos. Search the forum for 'mosquito', this should be a good starting point. cutter |
Thanks, for the reply, I had already seen that and I am currently going through Frank Bos thesis. However my case is slightly different because I do not have any coordinate for the different positions assumed by the wing and I have no idea how to proceed in this sense :confused:
|
Greetings to all!
@Andrea: I'll quote your question from the other thread: Quote:
Beyond this, in case you don't want to use a CAD software, basically what you're looking for is a deformed circle. Check the tutorials that use the feature "arc" that blockMesh has got. For example, the 2nd tutorial on the User Guide is one of those cases: http://www.openfoam.org/docs/user/plateHole.php Best regards, Bruno |
1 Attachment(s)
Quote:
I have been going through the Blender tutorial, however a cannot get any good blockMeshdict or snappyHex one. In the image my simple geometry is shown (.001m thick). I am probably doing something wrong while setting the patches because in the BMdict I only get a vast number of vertices and no actual patches coordinates. I have tried both autodetect features and selecting manually the edges. Note: the geometry has been imported as an stl file from Solidworks |
Hi Andrea,
There isn't much I can do to help, without at least an example case of what you're doing. By the way, did you read this tutorial: http://openfoamwiki.net/index.php/Tut_andrew_blender :confused: Best regards, Bruno |
I had a look at that but for the moment I decided to focus on a 2D case by applying the RBF function to the basic airfoil in the tutorial folder (pimpleDyMFoam).
BTW can RBF motion solver be applied to 3d models? Assuming that I would like to create a full 3D model, and make 2 patches out of the wings, can I apply that movement? Thank you so much!!! Basically, how do I do something similar to this http://www.youtube.com/watch?v=i6ECzy0jYxE ? :( |
Hi Andrea,
According to what I briefly saw in Frank Bos' thesis - http://repository.tudelft.nl/assets/...daa1b2/bos.pdf - it does seem to work in 3D as well. Best regards, Bruno |
1 Attachment(s)
Quote:
I read his work as well, I did get decent result in the 2D case. Im now into the 3D case but cant get the simple mesh. I keep getting the error Code:
--> FOAM FATAL ERROR: Code:
castellatedMesh true; :( Edit: I changed the blockMeshDict to a different one and now I dont get that message anymore but snappyHexMesh run. It finisher very very soon and effectively doesnt do anything as the only thing I can see in paraview is the box. I can open the geometry itself and it is fine and fits conpletely in the box. I attached snappy's Log. Andrea |
Hi Andrea,
Have a look into the tutorials listed here: http://openfoamwiki.net/index.php/Sn...als_and_Guides Best regards, Bruno |
4 Attachment(s)
I went trough those but we no success. Same result, no visible mesh.
I also tried using HelyxOS: 1)created new case 2) Left base mesh settings unchanged 3) in the volume option i did import .stl and imported my stl Tried both outside and inside but it says it is not supported (cause the geometry might not be closed??) 4) Proceeded to write the mesh. 5) I can see only the box, there is clearly something wrong with my execution of snappy but, due to my limited knowledge, I have no idea. 6) I also tried with a sphere and got nothing as well, so it is me. Did I miss any step? I left the levels of refinement unchanged cause it did not let me do neither inside or outside (with the sphere I can see that the box gets refined but cant see any sphere inside). I attached some picture which also show the hummingbird.stl inside the searchable box. Apologies :(:( Andrea |
Hi Andrea,
The quickest way to check if the problem is with the STL file is to run this command: Code:
surfaceCheck hummingbird.stl
Bruno |
Sorted thanks, I realised I was just being an idiot and messed up with the inside/outside refinement boxes etc. I got the geometry in now and I will start soon. :D:D
By the way, why if I interrupt a solver while is running from HelyxOs then is compromised and wont run a second time? :mad: Regards Andrea |
Hi Andrea,
Quote:
And please detail for them the specific steps you take, so that they can reproduce the error! ;) Best regards, Bruno |
Dear Bruno,
I am getting along fairly well!! However, now the meshes are getting coarser and my machine cannot handle them any more. Consequently, I got access to the University HPC. The cluster has OpenFoam installed but sadly no extension! My supervisor said that we can get the solver and the function installed somehow. Now, is there any way in which I can simply compile the function before every simulation or do I have to install it before? Is it the same for the solver? I need icoDyMFoam and RBFMotion Funtion. What would you recommend? Thanks!!! :D:D Regards Andrea |
Hi Andrea,
Quote:
Quote:
"Coarser" means that there are a lesser number of cells and therefore less precision in the discretization ;) Hence "coarser". Quote:
Quote:
In theory, it's just a matter of similarly following the steps from chapter 2 on this tutorial: http://openfoamwiki.net/index.php/Ho...ure_to_icoFoam Which specific version(s) of OpenFOAM do they have on their cluster? I can try to have a look into this next weekend. Best regards, Bruno |
1 Attachment(s)
Quote:
However, I did encounter a major issue. As showed in the picture, the plate is moving as it should and it did work fine on the modified geometry of the movingblockRBF case. However, I had to enlarge the domain and reduce the dimension of the block. I am using the same kinematic model I tested on the tutorial but after approx .2 seconds (62000 cells, approx 3 minutes after starting) of simulation, it crashes since the velocity goes to an outrageous order of magnitude (10^34 etc). My hypothesis is that since the mesh gets to squashed I get numerical errors. I noticed that by modifying some parameters in the dynamic mesh dict I can delay the problem but i can stop it. Do you know where I can find some tips on how to avoid this? Thanks!! :confused::confused: |
Hi Andrea,
I think the problem is because the cells are getting extremely compressed on the right side of the mesh. You need to increase the region through which the mesh can expand, so that it can distribute the cell compression along the cells to its right. Problem is that I'm not familiar with this specific mesh motion, so I don't know which parameters need tuning :( Can you share the contents of the file "constant/dynamicMeshDict"? And indicate which parameters you've been changing? Best regards, Bruno |
Dear Bruno, Thanks for the quick reply!!
Here is my dict Code:
dynamicFvMesh dynamicMotionSolverFvMesh; |
Hi Andrea,
The two most suspicious parameters are: Code:
outerRadius 1.5; Code:
IMQBCoeffs But as I said, I'm not familiar with this dynamic mesh motion, so this is just an educated guess :). Best regards, Bruno |
All times are GMT -4. The time now is 09:17. |