CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Mosquito Analysis in OpenFoam

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

Like Tree15Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   June 20, 2013, 03:11
Default Mosquito Analysis in OpenFoam
  #1
psk
Member
 
david
Join Date: Jun 2013
Location: Montreal, Canada
Posts: 62
Rep Power: 3
psk is an unknown quantity at this point
Greetings OF friends,

I am trying to do cfd analysis in mosquito model, by using cameras we capture frame by frames and extract the coordinates (X Y Z) of model.

One Wing Beat = 10 frames. Each frame we have the coordinates of model.I am trying to understand the physics and to find lift & drag.

I have only little experience in open foam. I did tutorial problems from that I understand how to use and I hope I get enough confident.

Generally, it covers unsteady, Dynamic & moving mesh,Fluid solid Interaction.

The problem is I don't have STL file format, I have only coordinates.What i have to do, How i can start with these coordinates, whether i have to write any script file to import? I install OF 2.1 kindly secure and guide me.


Many Regards,
david
psk is offline   Reply With Quote

Old   June 20, 2013, 07:56
Default
  #2
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
I am unsure what your problem is. Is the main issue to convert coordinates to STL? Google supplied me with this script: http://www.mathworks.com/matlabcentr...l-surface-mesh
Bernhard is offline   Reply With Quote

Old   June 20, 2013, 08:56
Default
  #3
Senior Member
 
Elvis
Join Date: Mar 2009
Location: Sindelfingen, Germany
Posts: 577
Blog Entries: 5
Rep Power: 13
elvis is on a distinguished road
Hi,

I remember that Frank Bos developed Radial Basis Function (RBF) with Insect flight in mind

http://ftp.jaist.ac.jp/pub/sourcefor...h/frankBos.pdf
Moving and deforming meshes for flapping flight at low Reynolds numbers

You find Radial Basis Function (RBF) included with Open FOAM 1.6-extend
Warning you will not find Radial Basis Function (RBF) included on a non EXTEND-Version of OpenFOAM!. So if you want to use GGI of RBF or Block Coupled Matrix Solver or Finite Area Method (FAM) you have to use OF-1.6-EXTEND-Version

The PHD-Thesis from Frank Bos is also available online
http://repository.tudelft.nl/assets/...daa1b2/bos.pdf (size 11MB)
psk likes this.
elvis is offline   Reply With Quote

Old   June 20, 2013, 09:08
Default
  #4
psk
Member
 
david
Join Date: Jun 2013
Location: Montreal, Canada
Posts: 62
Rep Power: 3
psk is an unknown quantity at this point
Hai Bernhard, thanks for coming to help me. My interest is flapping & flexible wing MAV.I will tell you step by step what I did.


·My goal is to find lift for flexible flapping model of my wing (see the Image).



·For that what i did is, I fabricate the wing and I initially consider 5 points on the wing (1, 2, 3, 4, 5) by using fluorescent ink to mark the points on the wing.


·By using the camera I capture the wing movement by frame by frame. Camera captures 100 frames per second. One complete wing beat it takes 10 frames. See the link [ http://www.unc.edu/~thedrick/software1.html]. I am

following the same procedure of Hedrick lab what they did. Please see the link.


·From the digitilization, I extract the XYZ coordinates of my wing by frame by frame. Afterwards I go for mat lab to plot the XYZ coordinates to see the replication my wing.



From the mat lab how I can incorporate open foam to find the lift of wing and how I start with Pre-Processing (Meshing)?




Many thanks,
David
Attached Images
File Type: jpg Image.jpg (8.4 KB, 66 views)
psk is offline   Reply With Quote

Old   June 20, 2013, 09:30
Default
  #5
psk
Member
 
david
Join Date: Jun 2013
Location: Montreal, Canada
Posts: 62
Rep Power: 3
psk is an unknown quantity at this point
HAI ELVIS,

Thanks for coming to help me,Please see my above post what i did up to right now. currently I have mat lab file, I don't know how to incorporate Open Foam with mat lab for Meshing,kindly give your feedback.


Many thanks,
david
psk is offline   Reply With Quote

Old   June 20, 2013, 10:28
Default
  #6
Senior Member
 
Elvis
Join Date: Mar 2009
Location: Sindelfingen, Germany
Posts: 577
Blog Entries: 5
Rep Power: 13
elvis is on a distinguished road
Hi David,

i am not an expert on that topic. I just knew that OF was used for insect wing CFD.
And that RBF was implemented to handle mesh deformation as comfortable as possible for such a tricky application.

If RBF is the right mesh motion technigue you want to use for your simulation, I guess that the marker dot you filmed and extracted with mat lab are the stearing points for RBF to do mesh deformation. I do not know what the best tool to do meshing
elvis is offline   Reply With Quote

Old   June 20, 2013, 12:29
Default Meshing Mosquito in Openfoam
  #7
psk
Member
 
david
Join Date: Jun 2013
Location: Montreal, Canada
Posts: 62
Rep Power: 3
psk is an unknown quantity at this point
Thanks elvis,Correct I extract the xyz coordinates by using digitizing [ DLTdv5 MATLAB code].

I am looking for only meshing part not solving, first thing I have to mesh properly.I don't know how to start in open foam for meshing.Anybody have idea/suggestions kindly help me.


Many Thanks
David

psk is offline   Reply With Quote

Old   June 20, 2013, 14:40
Default
  #8
Senior Member
 
Elvis
Join Date: Mar 2009
Location: Sindelfingen, Germany
Posts: 577
Blog Entries: 5
Rep Power: 13
elvis is on a distinguished road
Dear David,

I would start to get OF1.6-ext installed! If RBF is the tool you want to use for your simulation. You are lucky OF1.6-extend is back online I would study the RBF tutorials that are added, and when you came to that point, I would watch out for which tool is the right Meshingsoftware
elvis is offline   Reply With Quote

Old   June 21, 2013, 02:21
Default Meshing Mosquito in Openfoam
  #9
psk
Member
 
david
Join Date: Jun 2013
Location: Montreal, Canada
Posts: 62
Rep Power: 3
psk is an unknown quantity at this point
Hai elvis, thanks for coming again,

I go through Frank Bos post in OF forum, I think he used OF 1.2 & OF 1.3 [
icoFoamAutoMotion and interFoam] and I thought,finally he goes OF 1.6 ext.I hope OF 1.2, OF 1.3 and OF 1.6 ext may be help this problem.But I concentrate only on meshing part,I am totally confused about Pre- processing of OF which one is useful for this kind of problem. I don't exactly but still I am searching for good meshing software.

For my problem I have to use any one of the following Meshing.

1.Moving
mesh
2.Dynamic mesh
3.ALE mesh
4.Arbitrary mesh interface


Kindly give your feedback.

Many Thanks,
david
psk is offline   Reply With Quote

Old   June 21, 2013, 05:44
Default
  #10
Senior Member
 
Elvis
Join Date: Mar 2009
Location: Sindelfingen, Germany
Posts: 577
Blog Entries: 5
Rep Power: 13
elvis is on a distinguished road
Hi David,

if I remember correctly RBF was introduced with OF-1.5-dev and it was transfered to latest "community Version" called OF1.6-extend. So you do have the choice to use OF1.5-dev or OF1.6-extend => I would try the supported Version 1.6-extend.

To tell you a little bit about dynamic mesh handling ;-)
Arbitrary mesh interface (AMI) is a mesh technology that was introduced with OpenFOAM® v2.1.0:
Arbitrary Mesh Interface (AMI) it is "OpenCFD Ltd´s" "answer" to the Generalized Grid Interface (GGI) mesh technology that was introduced by Hrv Jasak and Martin Beaudoin in OF1.5-dev & OF1.6-extend

You find s.th similar http://www.openfoam.org/version2.0.0/surface-film.php and (Finite Area Method (FAM)
looking for a thin film model

http://powerlab.fsb.hr/ped/kturbo/Op...Manuscript.pdf
I would suggest explore http://powerlab.fsb.hr/ped/kturbo/OpenFOAM/ folder papers & docs
immortality and psk like this.

Last edited by elvis; June 21, 2013 at 08:05.
elvis is offline   Reply With Quote

Old   June 21, 2013, 10:02
Default
  #11
psk
Member
 
david
Join Date: Jun 2013
Location: Montreal, Canada
Posts: 62
Rep Power: 3
psk is an unknown quantity at this point
Thanks a lot elvis for your kind help & information,

It takes time for me to read and get some idea.Keep in touch


Many thanks,
David
psk is offline   Reply With Quote

Old   June 21, 2013, 12:35
Default
  #12
Senior Member
 
Elvis
Join Date: Mar 2009
Location: Sindelfingen, Germany
Posts: 577
Blog Entries: 5
Rep Power: 13
elvis is on a distinguished road
Hi David,

I think tutorial flappingConsoleSmall is what you should try
/OpenFOAM/OpenFOAM-1.6-ext/tutorials/stressAnalysis/icoFsiFoam/flappingConsoleSmall
elvis is offline   Reply With Quote

Old   June 21, 2013, 12:57
Default
  #13
psk
Member
 
david
Join Date: Jun 2013
Location: Montreal, Canada
Posts: 62
Rep Power: 3
psk is an unknown quantity at this point
Hai elvis thanks for kind information.

Currently I am go through OF 1.6 ext what you mention earlier post, For meshing I thought there is no need to go in-depth of Open Foam. Instead of that, may be use some structural stuff meshing software for meshing. In case meshing is tough in open Foam.

For solving part to use Open Foam,Please give your feedback about my idea, whether its correct or wrong approach?

Many Thanks,
David
psk is offline   Reply With Quote

Old   June 24, 2013, 08:55
Default icoDyMFoam error
  #14
psk
Member
 
david
Join Date: Jun 2013
Location: Montreal, Canada
Posts: 62
Rep Power: 3
psk is an unknown quantity at this point
Hai friends,
I install open foam 1.6 ext in Ubuntu 12.04 LTS I try to run icoDyMFoam tutorials it shows error like this,

Case : /home/david/OpenFOAM/OpenFOAM-1.6-ext/tutorials/incompressible/icoDyMFoam/movingBlockRBF
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Creating block mesh from
"/home/david/OpenFOAM/OpenFOAM-1.6-ext/tutorials/incompressible/icoDyMFoam/movingBlockRBF/constant/polyMesh/blockMeshDict"
Creating curved edges
Creating topology blocks
Creating topology patches

Reading patches section

Creating block mesh topology

Reading physicalType from existing boundary file

Default patch type set to empty
--> FOAM Warning :
From function blockMesh::createTopology(IOdictionary&)
in file blockMesh/blockMeshTopology.C at line 517
Reading "/home/david/OpenFOAM/OpenFOAM-1.6-ext/tutorials/incompressible/icoDyMFoam/movingBlockRBF/constant/polyMesh/blockMeshDict" from line 17 to line 150
For patch top overriding type 'patch' with 'wall' (read from boundary file)

Check topology

Basic statistics
Number of internal faces : 8
Number of boundary faces : 32
Number of defined boundary faces : 32
Number of undefined boundary faces : 0
Checking patch -> block consistency

Creating block offsets
Creating merge list .

Creating polyMesh from blockMesh
Creating patches
Creating cells
Creating points with scale 1

Writing polyMesh
----------------
Mesh Information
----------------
boundingBox: (-12.5 -12.5 -0.1) (12.5 12.5 0.1)
nPoints: 1352
nCells: 620
nFaces: 2536
nInternalFaces: 1184
----------------
Patches
----------------
patch 0 (start: 1184 size: 12) name: block
patch 1 (start: 1196 size: 25) name: top
patch 2 (start: 1221 size: 25) name: bottom
patch 3 (start: 1246 size: 25) name: left
patch 4 (start: 1271 size: 25) name: right
patch 5 (start: 1296 size: 1240) name: frontAndBack

End

david@ubuntu:~/OpenFOAM/OpenFOAM-1.6-ext/tutorials/incompressible/icoDyMFoam/movingBlockRBF$ icoDyMFoam
icoDyMFoam: command not found
david@ubuntu:~/OpenFOAM/OpenFOAM-1.6-ext/tutorials/incompressible/icoDyMFoam/movingBlockRBF$


Kindly guide me


Many Thanks,
david
psk is offline   Reply With Quote

Old   June 24, 2013, 14:55
Default
  #15
Senior Member
 
Elvis
Join Date: Mar 2009
Location: Sindelfingen, Germany
Posts: 577
Blog Entries: 5
Rep Power: 13
elvis is on a distinguished road
Dear David,

did you follow http://openfoamwiki.net/index.php/In...1.6-ext/Ubuntu ? in thread OpenFOAM-1.6-ext on Ubuntu 12.04 install errors there used to be a link to Martin Beaudoin OF-Ext Install Guideline that is not existing anymore.

I would suggest to try a live DVD-Image http://rcc.its.psu.edu/ofw6/ofw6_v8.iso from 6th OF Workshop to have a 100% working Extend-version
immortality and psk like this.
elvis is offline   Reply With Quote

Old   June 24, 2013, 16:50
Default
  #16
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,253
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings to all!

@elvis:
Quote:
Originally Posted by elvis View Post
did you follow http://openfoamwiki.net/index.php/In...1.6-ext/Ubuntu ? in thread OpenFOAM-1.6-ext on Ubuntu 12.04 install errors there used to be a link to Martin Beaudoin OF-Ext Install Guideline that is not existing anymore.
In the section "2.3 Ubuntu 12.04" at http://openfoamwiki.net/index.php/In...u#Ubuntu_12.04 , it gives the links to 3 copies of Martin Beaudoin's guides on how to install OpenFOAM 1.6-ext on Ubuntu 12.04!

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   June 25, 2013, 10:36
Default
  #17
psk
Member
 
david
Join Date: Jun 2013
Location: Montreal, Canada
Posts: 62
Rep Power: 3
psk is an unknown quantity at this point
Thanks elvis and wyldckat


regards,
david
psk is offline   Reply With Quote

Old   June 25, 2013, 14:38
Default
  #18
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
Hi psk
is the usage of mosquite in optimization of geometry or only for moving mesh?
and whats its difference from Dakota if it can be used for geometry optimization?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   June 26, 2013, 01:40
Default
  #19
psk
Member
 
david
Join Date: Jun 2013
Location: Montreal, Canada
Posts: 62
Rep Power: 3
psk is an unknown quantity at this point
Thanks for your kind response immortality, its only Moving mesh / Dynamic mesh.

I didn't understand [whats its difference from Dakota if it can be used for geometry optimization] what it means , If you don't mind kindly let me know.

kindly go through my post if you have any idea, let me know. As per elvis and other OF friends instructions i am working in Open Foam 1.6 ext.

many thanks,
david
psk is offline   Reply With Quote

Old   July 11, 2013, 14:26
Default
  #20
psk
Member
 
david
Join Date: Jun 2013
Location: Montreal, Canada
Posts: 62
Rep Power: 3
psk is an unknown quantity at this point
Good Morning folks and elvis, this is david , last two weeks i struck with my personal commitments, I install OF 1.6 and i run IcoDyMFoam it works fine.

Before start my work, I would like to start with rigid flat wing with flapping amplitude of -10 to + 10

Mosquito wing- FSI problem, I thought its little complicated, Initially i would like to start with rigid then go for flexible.

If you have any tutorials related rigid flat wing,kindly let me know

How I can Proceed?

Many Thanks,
david
psk is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
OpenFOAM stress analysis results are showing 8-10% error with analytical solution. Sargam05 OpenFOAM 2 September 10, 2012 15:06
2D Mesh Generation Tutorial for GMSH aeroslacker Open Source Meshers: Gmsh, Netgen, CGNS, ... 12 January 19, 2012 04:52
use of openfoam in attic analysis kirankarki OpenFOAM 1 September 12, 2010 21:02
OpenFOAM 1.5.x package - CentOS 5.3 x86_64 linnemann OpenFOAM Installation 7 July 30, 2009 03:14
OpenFOAM Install problem masb OpenFOAM 3 May 25, 2009 11:32


All times are GMT -4. The time now is 10:40.