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

CL history for pitching airfoil

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By star

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 14, 2014, 05:16
Default CL history for pitching airfoil
  #1
Member
 
Harry
Join Date: May 2013
Posts: 68
Rep Power: 13
star is on a distinguished road
Hi,
In pitching airfoil case the angle of attack changes with time. In Solve/Monitors, we can plot CL history but that is for some fixed angle of attack (as there is no udf option there). Can anybody tell me how to get CL history in this case?
Thanks in advance
star is offline   Reply With Quote

Old   February 14, 2014, 05:35
Smile
  #2
Senior Member
 
Anna Tian's Avatar
 
Meimei Wang
Join Date: Jul 2012
Posts: 494
Rep Power: 16
Anna Tian is on a distinguished road
Quote:
Originally Posted by star View Post
Hi,
In pitching airfoil case the angle of attack changes with time. In Solve/Monitors, we can plot CL history but that is for some fixed angle of attack (as there is no udf option there). Can anybody tell me how to get CL history in this case?
Thanks in advance
Monitor -> edit 'Cl' -> click on 'write' and fill in the file name

Then after running the simulation, you can find the cl-history file in the case folder. That is the same for the changing AoA.

Good luck!
__________________
Best regards,
Meimei
Anna Tian is offline   Reply With Quote

Old   February 14, 2014, 11:45
Default
  #3
Member
 
Harry
Join Date: May 2013
Posts: 68
Rep Power: 13
star is on a distinguished road
Thanks Meimei for your reply. I am not sure whether i should write the X and Y components in CL because it depends on angle of attack (in my case it's changing) but i will try it.
star is offline   Reply With Quote

Old   February 22, 2014, 05:44
Default
  #4
New Member
 
Musa
Join Date: Mar 2012
Posts: 14
Rep Power: 14
fkanyako is on a distinguished road
Hi Star,
I am also running simulation pitching airfoil
using the equation a(t) = a0 +a1 * sin(wt).
where a0 = 15, a1 = 10, (ao is the mean angle of attack). I took the time derivative of this equation to get the angular velocity and I wrote the udf like this omega(2) = w *a1* cos(wt). As for the a0 i set it in fluent that the velocity is coming at 15 degrees which is the mean angle of attack.
x = -5sin(15)
y = 5cos(15) let say free stream velocity is 5.
I also set the force vectors for lift and drag under solve/monitor. at this same angle 15 deg.

But the lift coefficient history when plotted against the angle of attack its not comparable to experiment. Please let me know if am writing the udf correctly or do you mind just to show me your udf... thanks
fkanyako is offline   Reply With Quote

Old   February 22, 2014, 07:10
Default
  #5
Member
 
Harry
Join Date: May 2013
Posts: 68
Rep Power: 13
star is on a distinguished road
Quote:
Originally Posted by fkanyako View Post
Hi Star,
I am also running simulation pitching airfoil
using the equation a(t) = a0 +a1 * sin(wt).
where a0 = 15, a1 = 10, (ao is the mean angle of attack). I took the time derivative of this equation to get the angular velocity and I wrote the udf like this omega(2) = w *a1* cos(wt). As for the a0 i set it in fluent that the velocity is coming at 15 degrees which is the mean angle of attack.
x = -5sin(15)
y = 5cos(15) let say free stream velocity is 5.
I also set the force vectors for lift and drag under solve/monitor. at this same angle 15 deg.

But the lift coefficient history when plotted against the angle of attack its not comparable to experiment. Please let me know if am writing the udf correctly or do you mind just to show me your udf... thanks
Before running, you must check pitch motion through Display Zone Motion. If it works as you want then udf is ok. I am not expert in this. you can find more udf's in this forum. Here i have attached mine. It works fine.

#include <stdio.h>
#include <math.h>
#include "udf.h"

#define pi 3.141592653589793 //

DEFINE_CG_MOTION(airfoil_pitch, dt, vel, omega, time, dtime)
{
real pitchampl, pitchdispl, omg, count;

/* reset velocities */
NV_S (vel, =, 0.0);
NV_S (omega, =, 0.0);

/* motion */
pitchampl = 4 * pi / 180; /* angular pitch amplitude */
omg = 2 * pi * 40; /* angular frequency */
pitchdispl = pitchampl * sin(omg * time); /* resultant displacement */

omega[2] = pitchampl * omg * cos(omg * time);
blgypeng likes this.
star is offline   Reply With Quote

Old   February 22, 2014, 07:50
Default
  #6
New Member
 
Musa
Join Date: Mar 2012
Posts: 14
Rep Power: 14
fkanyako is on a distinguished road
Thanks a lot Star.
Looking your UDF, we basically have the same. When I preview mesh motion its working fine. But the result is not so good. The problem am having maybe its with the setting especially the boundary condition. are you ptiching your airfoil at any mean angle of attack ?.

Also I get the lift coefficient between 0.4 and -0.1 vs time, sinusoidal. But the experimental is between 3.0 and -0.5 vs time. If you manage to validate your result please explain the settings and also how u computed the lift coefficient vs angle of attack. Thanks
fkanyako is offline   Reply With Quote

Old   February 22, 2014, 09:00
Default
  #7
Member
 
Harry
Join Date: May 2013
Posts: 68
Rep Power: 13
star is on a distinguished road
It depends on no. of factors. This question is not specific. Quality of your mesh, Reynolds number, time step size and other boundary conditions. You should search previous posts related to your problem.
[Be sure you are using correct pitching amplitude and frequency]
star is offline   Reply With Quote

Old   August 4, 2017, 08:20
Default
  #8
New Member
 
flanag
Join Date: Jun 2017
Posts: 12
Rep Power: 9
lfla is on a distinguished road
I am currently doing something similar and i have set inlet velocity in x and y components. x=30cos8 and y:30sin8.
I applied the udf to the fluid using mesh motion and selecting the udf

#include <stdio.h>
#include <math.h>
#include "udf.h"

DEFINE_CG_MOTION(airfoil_pitch, dt, vel, omega, time, dtime)
{
real pitchampl, pitchdispl, omg, pi, count;

pi = 3.14159265;


/* reset velocities */
NV_S(vel, =, 0.0);
NV_S(omega, =, 0.0);

/* motion */
pitchampl = 5.5 * pi / 180; /* angular pitch amplitude */
omg = 2 * pi * 1.2; /* angular frequency */
pitchdispl = pitchampl * sin(omg * time); /* resultant displacement */

omega[2] = pitchampl * omg * cos(omg * time);
}

Do you have any advice for, thank you i would reallu appreciate it
Attached Images
File Type: jpg mesh1.jpg (25.3 KB, 9 views)
File Type: jpg inlet.jpg (117.9 KB, 12 views)
File Type: jpg cell zone.jpg (88.7 KB, 13 views)
lfla is offline   Reply With Quote

Old   August 4, 2017, 08:24
Default
  #9
New Member
 
flanag
Join Date: Jun 2017
Posts: 12
Rep Power: 9
lfla is on a distinguished road
Quote:
Originally Posted by fkanyako View Post
Hi Star,
As for the a0 i set it in fluent that the velocity is coming at 15 degrees which is the mean angle of attack.
x = -5sin(15)
y = 5cos(15) let say free stream velocity is 5.
I also set the force vectors for lift and drag under solve/monitor. at this same angle 15 deg.
Hi fkanyako,
Could i ask you why did you not use x=5cos15 and y 5sin15 Thank you
lfla is offline   Reply With Quote

Reply


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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Low Speed Airfoil Mancusi FLUENT 7 April 3, 2014 07:11
Problem with restart solution in shape_optimization.py robyTKD SU2 Shape Design 21 May 29, 2013 10:26
[GAMBIT] Meshing airfoil using .dat file problem creggie ANSYS Meshing & Geometry 10 June 27, 2010 20:24
Modeling Backflow for a 3D Airfoil (Wing of Finite Span) Josh CFX 9 August 18, 2009 12:31
Airfoil boundary condition Frank Main CFD Forum 1 April 21, 2008 19:36


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