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

Lift, Drag Vs time chart,calculations

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 15, 2005, 07:18
Default Lift, Drag Vs time chart,calculations
  #1
Jamesd69climber
Guest
 
Posts: n/a
Hi

I am interested in calculating the lift and drag coefficients over time for a bluff body.

Is it possible to have such a graph made in post processor or solver, or do I need to calculate the coefficients manually for each timestep?

Regards

James
  Reply With Quote

Old   February 15, 2005, 08:00
Default Re: Lift, Drag Vs time chart,calculations
  #2
James Date
Guest
 
Posts: n/a
You need to create expressions for Cl, Cd and then monitor them in the output menu.

This should give you an idea of what to do.

The expressions go something like this:

LIBRARY:

CEL:

EXPRESSIONS:

Angle = 0 [deg]

Span = 0.1 [m]

Chord = 0.18574 [m]

SectionArea = Chord*Span

Drag = ((force_x()@Section)*cos(Angle))+((force_y()@Secti on)*sin(Angle))

Free Stream = 4.1 [m s^-1]

CD = Drag/(0.5*1027.6*SectionArea*Free Stream^2)

Lift = (-1.0*(force_x()@Section)*sin(Angle))+((force_y()@Se ction)*cos(Angle))

CL = Lift/(0.5*1027.6*SectionArea*Free Stream^2)

CLCD = Lift/Drag

CM = torque_z()@Section/(0.5*1027.6*SectionArea*Free Stream^2)

NormalF = (Lift*cos(Angle))+(Drag*sin(Angle))

CNN = NormalF/(0.5*1027.6*SectionArea*Free Stream^2)

CPX = torque_z()@Section/NormalF

FX = force_x()@Section

FY = force_y()@Section

XTorque = torque_x()@Section

YTorque = torque_y()@Section

ZTorque = torque_z()@Section

uu = Free Stream*cos(Angle)

vv = Free Stream*sin(Angle)

ww = 0 [m s^-1]

END

END

The monitor points go something like this:

OUTPUT CONTROL:

MONITOR OBJECTS:

MONITOR BALANCES:

Option = Full

END

MONITOR FORCES:

Option = Full

END

MONITOR PARTICLES:

Option = Full

END

MONITOR POINT: CLCD Ratio

Expression Value = CLCD

Option = Expression

END

MONITOR POINT: Centre of Pressure

Expression Value = CPX

Option = Expression

END

MONITOR POINT: Drag Coefficient

Expression Value = CD

Option = Expression

END

MONITOR POINT: Drag Force

Expression Value = Drag

Option = Expression

END

MONITOR POINT: Force X

Expression Value = FX

Option = Expression

END

MONITOR POINT: Force Y

Expression Value = FY

Option = Expression

END

MONITOR POINT: Lift Coefficient

Expression Value = CL

Option = Expression

END

MONITOR POINT: Lift Force

Expression Value = Lift

Option = Expression

END

MONITOR POINT: Moment Coefficient

Expression Value = CM

Option = Expression

END

MONITOR POINT: Normal Coefficient

Expression Value = CNN

Option = Expression

END

MONITOR POINT: Normal Force

Expression Value = NormalF

Option = Expression

END

MONITOR POINT: TorqueX

Expression Value = XTorque

Option = Expression

END

MONITOR POINT: TorqueY

Expression Value = YTorque

Option = Expression

END

MONITOR POINT: TorqueZ

Expression Value = ZTorque

Option = Expression

END

MONITOR RESIDUALS:

Option = Full

END

MONITOR TOTALS:

Option = Full

END

END

RESULTS:

File Compression Level = Default

Option = Full

Output Boundary Flows = All

Output Equation Residuals = All

Output Variable Operators = All

END

END

James

  Reply With Quote

Old   February 15, 2005, 08:19
Default Re: Lift, Drag Vs time chart,calculations
  #3
Jamesd69climber
Guest
 
Posts: n/a
James

Thank you for your reply, most appreciated!

James
  Reply With Quote

Old   February 15, 2005, 09:43
Default Re: Lift, Drag Vs time chart,calculations
  #4
Ogbeni
Guest
 
Posts: n/a
James, Excellent !, any chance you can define some of the variables ? and some explantion of the commands , for instance what's 1027.6 related to or I am asking too much,
  Reply With Quote

Old   February 15, 2005, 15:25
Default Re: Lift, Drag Vs time chart,calculations
  #5
James Date
Guest
 
Posts: n/a
You're best bet is to have a look at the CCL guide in the CFX user manual to work out any of the commands which I've used which aren't familiar.

Angle = Angle of Incidence

Span = Span of Section

Chord = Chord of Section

Free Stream = Free Stream Flow Velocity

CD = Drag Coefficient

CL = Lift Coefficient

CNN = Normal Force Coefficient

CLCD = Lift/Drag Ratio

CM = Moment Coefficient

CPX = Centre of Pressure in the x-direction

uu, vv & ww = Flow Boundary Condition Velocities

1027.6 is the density of salt water at 5 deg C.

James
  Reply With Quote

Old   February 16, 2005, 15:29
Default Re: Lift, Drag Vs time chart,calculations
  #6
Akin
Guest
 
Posts: n/a
James, I am not sure CFX produce a force_x,force_y or z for that matter. Also the statement @section, refers to the wall in question, please clarify.
  Reply With Quote

Old   February 16, 2005, 18:58
Default Re: Lift, Drag Vs time chart,calculations
  #7
James Date
Guest
 
Posts: n/a
The Statement:

force_x()@Section

Gives the total force in the x-direction acting on the wall defined with the name "Section".

James
  Reply With Quote

Old   February 17, 2005, 10:43
Default Re: Lift, Drag Vs time chart,calculations
  #8
Rui
Guest
 
Posts: n/a
Hi,

Take a look at Tutorial 20 page 414 (CFX-5.7 Documentation), where the force on the ball due to fluid flow, FFlow, is defined as force_y()Ball
  Reply With Quote

Old   February 17, 2005, 18:23
Default Re: Lift, Drag Vs time chart,calculations
  #9
Akin
Guest
 
Posts: n/a
Cheers Janes and Rui
  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
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
Multiple floating objects CKH OpenFOAM 10 September 21, 2011 23:13
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03
Convergence moving mesh lr103476 OpenFOAM Running, Solving & CFD 30 November 19, 2007 15:09


All times are GMT -4. The time now is 07:33.