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

Lift and drag coefficients of a flapping wing aircraft

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

Like Tree2Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   August 5, 2012, 12:15
Arrow Lift and drag coefficients of a flapping wing aircraft
  #1
Member
 
Join Date: Aug 2012
Posts: 30
Rep Power: 4
Julian121 is on a distinguished road
Dear all,

I am simulating the flapping motion of a wing using Fluent to calculate lift and drag coefficients. I have used a UDF code and dynamic mesh to simulate the flapping motion.

I am simulating a rigid wing which has pitch and roll rotations and have defined four boundary conditions:
1- Pressure far field
2- Wall
3- Wall-solid
4- Interior-solid

I am defining the x,y,z components of the airflow using pressure far field boundary and assume that the flow is static and the wing moves.

What force direction should I use to calculate lift force?
Do I have to run the solution at steady first and then transient?

I would appreciate if anyone help me with this.
Julian121 is offline   Reply With Quote

Old   August 5, 2012, 14:08
Default
  #2
Senior Member
 
Join Date: Mar 2011
Location: Germany
Posts: 414
Rep Power: 11
cfd seeker is on a distinguished road
Quote:
Originally Posted by Julian121 View Post
Dear all,

I am simulating the flapping motion of a wing using Fluent to calculate lift and drag coefficients. I have used a UDF code and dynamic mesh to simulate the flapping motion.

I am simulating a rigid wing which has pitch and roll rotations and have defined four boundary conditions:
1- Pressure far field
2- Wall
3- Wall-solid
4- Interior-solid

I am defining the x,y,z components of the airflow using pressure far field boundary and assume that the flow is static and the wing moves.

What force direction should I use to calculate lift force?
Do I have to run the solution at steady first and then transient?

I would appreciate if anyone help me with this.
lift coefficient depends on angle of attack....in monitor panel give values of the cos and sine components at which you performing your analysis and run steady solution solution for static case without UDF so that flow develops and then switch to transient case and hook UDF to simulate flapping motion
cfd seeker is offline   Reply With Quote

Old   August 6, 2012, 03:09
Default
  #3
Member
 
Touré
Join Date: Aug 2010
Location: Montréal
Posts: 85
Rep Power: 6
Touré is on a distinguished road
I totally agree with cfd seeker. If you want to automate it, you can use parameters to change the angle but I wouldn't recommend that. The process would be:
1. Create a new project
2. Choose an angle alpha
3. Check the force Fy on Y-axis
Redo the process for several alpha until you get a curve or Lift versus alpha. It's better to do one alpha per project so that if it crash you wont loose all your other simulation.
For example, if you want to simulate alpha = 20, you can load a previous solution data of alpha = 10 so that the computing will be less difficult for the software to converge to the solution.
To find Fy (figure 1), go to Report -> forces ; Fy = Lift (Fx= Drag)
If you really want to use your UDF, you must save the data for the several times with Autosave (figure 2). Then, you will check the alpha and the corresponding Fy for each saved time data .
Attached Images
File Type: jpg Capture0.JPG (44.7 KB, 82 views)
File Type: jpg Capture00.JPG (48.0 KB, 68 views)
Touré is offline   Reply With Quote

Old   August 6, 2012, 16:26
Default
  #4
Member
 
Join Date: Aug 2012
Posts: 30
Rep Power: 4
Julian121 is on a distinguished road
Since I am simulating flapping motion of a wing and not a static airfoil, the lift force varies at different wing positions. How can I plot the lift force over time?
If I use Forces option in Fluent, it only gives lift at a particular wing position and not total lift over time.
Julian121 is offline   Reply With Quote

Old   August 6, 2012, 21:01
Default
  #5
Member
 
Touré
Join Date: Aug 2010
Location: Montréal
Posts: 85
Rep Power: 6
Touré is on a distinguished road
Fluent loads only the flow-field data of one time step and not all of them but maybe you can do this in FLUENT

Report -> forces -> Save Output Parameter

You will see the Parameter set at the bottom of the project (see figure).
Double click on it to check the values
Attached Images
File Type: jpg Capture0.JPG (21.7 KB, 33 views)
Rakeshvb likes this.
Touré is offline   Reply With Quote

Old   August 7, 2012, 05:51
Default
  #6
Member
 
Join Date: Aug 2012
Posts: 30
Rep Power: 4
Julian121 is on a distinguished road
Thank you for your reply. As the direction of lift changes during the wing motion, how can I define the direction of lift force while the direction changes at each wing position?
Julian121 is offline   Reply With Quote

Old   August 7, 2012, 08:27
Default
  #7
Member
 
Touré
Join Date: Aug 2010
Location: Montréal
Posts: 85
Rep Power: 6
Touré is on a distinguished road
Maybe you're talking about the projection of the Lift. The Drag (Fx) and Lift (Fy) which are orthogonal never change direction. The lift is always on the same direction as the gravity (vertical) and the drag is horizontal even if you wing moves. I don't think that FLUENT gives the projection of the Lift on the new position of you axis when they rotate and I don't think that you really need to compute it. I suggest simply Excel or Matlab to compute the projected Lift forces if you really want them by multiplying the Lift with cos(alpha) or sin(alpha) .
Touré is offline   Reply With Quote

Old   August 7, 2012, 09:54
Default
  #8
Senior Member
 
Join Date: Mar 2011
Location: Germany
Posts: 414
Rep Power: 11
cfd seeker is on a distinguished road
Quote:
Originally Posted by Julian121 View Post
Since I am simulating flapping motion of a wing and not a static airfoil, the lift force varies at different wing positions. How can I plot the lift force over time?
If I use Forces option in Fluent, it only gives lift at a particular wing position and not total lift over time.
I guess you are talking about the variation of lift as the wing flaps with time....its simple just go to Monitors and put the sine and cos components of angle of attack and also check plot and print tabs.....by this way you will get the variation of lift with time as it goes from positive to negative until a cyclic behavior which shows convergence of your problem....
cfd seeker is offline   Reply With Quote

Old   August 8, 2012, 12:06
Default
  #9
Member
 
Touré
Join Date: Aug 2010
Location: Montréal
Posts: 85
Rep Power: 6
Touré is on a distinguished road
Monitors -> Create (Lift)
Touré is offline   Reply With Quote

Old   August 9, 2012, 13:07
Default
  #10
Member
 
Join Date: Aug 2012
Posts: 30
Rep Power: 4
Julian121 is on a distinguished road
Quote:
Originally Posted by Touré View Post
Monitors -> Create (Lift)
I started the simulation and solved the problem in steady-state first. However, when I change the steady state to transient problem, the mesh quality decreases and I get an error:


Static mesh:
Minimum Orthogonal Quality = 0.6

After clicking on the preview motion:
Minimum Orthogonal Quality = 0.00000e+00

Warning: non-positive volumes exist.

Primitive Error at Node 0: Update-Dynamic-Mesh failed. Negative cell volume detected.

Primitive Error at Node 1: Update-Dynamic-Mesh failed. Negative cell volume detected.

Primitive Error at Node 2: Update-Dynamic-Mesh failed. Negative cell volume detected.

Primitive Error at Node 3: Update-Dynamic-Mesh failed. Negative cell volume detected.

Error: Update-Dynamic-Mesh failed. Negative cell volume detected.
Error Object: #f

Any idea what I should do?
Julian121 is offline   Reply With Quote

Old   August 10, 2012, 03:00
Default
  #11
Senior Member
 
Join Date: Mar 2011
Location: Germany
Posts: 414
Rep Power: 11
cfd seeker is on a distinguished road
Quote:
Originally Posted by Julian121 View Post
I started the simulation and solved the problem in steady-state first. However, when I change the steady state to transient problem, the mesh quality decreases and I get an error:


Static mesh:
Minimum Orthogonal Quality = 0.6

After clicking on the preview motion:
Minimum Orthogonal Quality = 0.00000e+00

Warning: non-positive volumes exist.

Primitive Error at Node 0: Update-Dynamic-Mesh failed. Negative cell volume detected.

Primitive Error at Node 1: Update-Dynamic-Mesh failed. Negative cell volume detected.

Primitive Error at Node 2: Update-Dynamic-Mesh failed. Negative cell volume detected.

Primitive Error at Node 3: Update-Dynamic-Mesh failed. Negative cell volume detected.

Error: Update-Dynamic-Mesh failed. Negative cell volume detected.
Error Object: #f

Any idea what I should do?
negative volumes are appearing during re-meshing process....decrease your time step size.........
cfd seeker is offline   Reply With Quote

Old   August 10, 2012, 03:39
Default
  #12
Member
 
Join Date: Aug 2012
Posts: 30
Rep Power: 4
Julian121 is on a distinguished road
Quote:
Originally Posted by cfd seeker View Post
negative volumes are appearing during re-meshing process....decrease your time step size.........

I have reduced the time step to 0.0001 but still does not work.
Julian121 is offline   Reply With Quote

Old   August 10, 2012, 09:42
Default
  #13
Senior Member
 
Join Date: Mar 2011
Location: Germany
Posts: 414
Rep Power: 11
cfd seeker is on a distinguished road
Quote:
Originally Posted by Julian121 View Post
I have reduced the time step to 0.0001 but still does not work.
try step size of 0.000001....what are your flow conditions?
cfd seeker is offline   Reply With Quote

Old   August 10, 2012, 11:54
Default
  #14
Member
 
Touré
Join Date: Aug 2010
Location: Montréal
Posts: 85
Rep Power: 6
Touré is on a distinguished road
It's maybe because of a coarse mesh. Maybe you need very small elements in the zone where you have your movement (trailing edge) at the beginning of the simulation. The dynamic mesh destroys the quality of your meshing, so do your best at the beginning.
Touré is offline   Reply With Quote

Old   August 10, 2012, 13:39
Default
  #15
Member
 
Join Date: Aug 2012
Posts: 30
Rep Power: 4
Julian121 is on a distinguished road
Quote:
Originally Posted by Touré View Post
It's maybe because of a coarse mesh. Maybe you need very small elements in the zone where you have your movement (trailing edge) at the beginning of the simulation. The dynamic mesh destroys the quality of your meshing, so do your best at the beginning.
I refined the mesh at the trailing edge but it did not work. Still, I get that error. I do not understand why Minimum Orthogonal Quality changes during the simulation. When the time step is 0.0001 the simulation goes up to 20 while the time step is 0.00001 it goes to 228 and then stops.

Is the error because of the mesh? I have tried tri/quad coarse/fine mesh but still does not work.
Julian121 is offline   Reply With Quote

Old   August 10, 2012, 14:33
Default
  #16
Member
 
Touré
Join Date: Aug 2010
Location: Montréal
Posts: 85
Rep Power: 6
Touré is on a distinguished road
Did you select "Remesh" when you have selected "Dynamic Mesh" ?
Use "Previewing the Dynamic Mesh" button, so that you see what is going on before running your simulation.
The "Minimum Orthogonal Quality" should change because the cells on the boundary of your airfoil moves making their angles to change also.
Touré is offline   Reply With Quote

Old   August 11, 2012, 19:07
Default
  #17
Member
 
Join Date: Aug 2012
Posts: 30
Rep Power: 4
Julian121 is on a distinguished road
Quote:
Originally Posted by Touré View Post
Did you select "Remesh" when you have selected "Dynamic Mesh" ?
Use "Previewing the Dynamic Mesh" button, so that you see what is going on before running your simulation.
The "Minimum Orthogonal Quality" should change because the cells on the boundary of your airfoil moves making their angles to change also.
Thank you very much Touré for your help.

I have got five zones in dynamic mesh zones: interior-solid, pressure_far_field, solid (fluid), wall and wall_solid

Which zone should be stationary or deforming? Should all of them be deforming except the wall-solid which does rigid body?

When I set every zone to stationary except the solid(fluid) to deforming and wall_solid to rigid body, I don't get any negative cell error.
Julian121 is offline   Reply With Quote

Old   August 11, 2012, 23:20
Default
  #18
Member
 
Touré
Join Date: Aug 2010
Location: Montréal
Posts: 85
Rep Power: 6
Touré is on a distinguished road
You 've got it.
Stationary is for the boundaries that are not moving like your far-field and the two symmetrical boundaries
Deforming is for a deforming boundary or zone such as your fluid
Rigid body is your airfoil and you have to link it with your motion after compiling the code.c file
Did you use "Peview Mesh Motion" before doing the computation?
Keep going
Touré is offline   Reply With Quote

Old   August 12, 2012, 03:28
Default
  #19
Member
 
Join Date: Aug 2012
Posts: 30
Rep Power: 4
Julian121 is on a distinguished road
Quote:
Originally Posted by Touré View Post
You 've got it.
Stationary is for the boundaries that are not moving like your far-field and the two symmetrical boundaries
Deforming is for a deforming boundary or zone such as your fluid
Rigid body is your airfoil and you have to link it with your motion after compiling the code.c file
Did you use "Peview Mesh Motion" before doing the computation?
Keep going
The problem still exists! When interior-solid is stationary, the minimum orthogonal quality remains constant and the error does not appear. However, when interior-solid is stationary, no motion is seen by using "preview motion".

What is the difference between solid (fluid) and interior-solid zones? Should it be deforming? When I change interior-solid to deforming, the error appear again.

The wing tips remain in the domain during the movement of the wing.

Last edited by Julian121; August 12, 2012 at 07:35.
Julian121 is offline   Reply With Quote

Old   August 12, 2012, 04:31
Default
  #20
Member
 
Join Date: Aug 2012
Posts: 30
Rep Power: 4
Julian121 is on a distinguished road
Quote:
Originally Posted by Touré View Post
You 've got it.
Stationary is for the boundaries that are not moving like your far-field and the two symmetrical boundaries
Deforming is for a deforming boundary or zone such as your fluid
Rigid body is your airfoil and you have to link it with your motion after compiling the code.c file
Did you use "Peview Mesh Motion" before doing the computation?
Keep going
The lift coefficient curve shows that it varies between 0.75 to 0.76 which seems to be wrong. Should not it go to minus value during downstroke? I have attached a screenshot of the lift curve. Do I have to do the calculation through "Run calculation", as when I preview the motion the lift curve starts drawing?
Attached Images
File Type: jpg lift.jpg (41.3 KB, 44 views)
Julian121 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
lift and drag coefficients around a ground vehicle Pedro CFX 3 September 5, 2012 18:31
Design of a flapping wing aircraft yahya_azizi FLUENT 0 July 11, 2012 11:26
Great lift coefficients but negative drag coefficient?! nez FLUENT 0 November 28, 2011 14:57
Lift and Drag on airplane wing Arj CFX 4 October 6, 2006 04:05
Lift and Drag Coefficients Reliability Luis FLUENT 2 December 27, 2005 15:45


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