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

how to write a scalar value to a file ????

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 7, 2013, 12:00
Default how to write a scalar value to a file ????
  #1
Member
 
Samer
Join Date: Jan 2013
Posts: 31
Rep Power: 4
SamerAli is on a distinguished road
Dear All
I am trying to calculate the lift and drag forces on a (cylinder + plate) and i would love to extract the total lift and the total drag and write them in a file at each writing interval, how would i do that ???
I am using Openfoam 1.6-ext, and this is the header file for lift and drag calculations :

scalar lift = 0;
scalar drag = 0;

const vectorField& Sfp = mesh.boundary()[plateID].Sf();
forAll(p.boundaryField()[plateID], facei)
{
vector faceForce = p.boundaryField()[plateID][facei] * Sfp[facei];
lift += vector(0,1,0) & faceForce;
drag += vector(1,0,0) & faceForce;
}

const vectorField& Sfc = mesh.boundary()[cylinderID].Sf();
forAll(p.boundaryField()[cylinderID], facei)
{
vector faceForce = p.boundaryField()[cylinderID][facei] * Sfc[facei];
lift += vector(0,1,0) & faceForce;
drag += vector(1,0,0) & faceForce;
}

how do i extract and write the scalar lift and the scalar drag to a file ??

Best regards
SamerAli is offline   Reply With Quote

Old   May 7, 2013, 13:19
Default
  #2
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,608
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi Samer,

It sounds like you are trying to accomplish the functionality in the forces and forceCoeffs functionObjects, which can be added at runTime.

Search the forum for its use and you can also find an example on how to extract the forces in this tutorial:

Code:
<path_to_tutorials>/incompressible/pimpleDyMFoam/wingMotion/wingMotion2D_simpleFoam/system/controlDict
There are also other tutorials using the forces, and you can find these by typing the following in the command line (with a sourced OF-version)

Code:
tut; find ./ -name "controlDict" -type f | xargs grep forces
Kind regards,

Niels
ngj 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
Helyx-OS (GUI for SnappyHexMesh elvis OpenFOAM Native Meshers: snappyHexMesh and Others 176 July 10, 2015 16:50
which tutorial can find that include Riemann boundary condition? immortality OpenFOAM Running, Solving & CFD 12 May 3, 2013 18:21
build problem swak4Foam OF 2.2.0 mcathela OpenFOAM Installation 14 April 23, 2013 13:59
mesh file for flow over a circular cylinder Ardalan Main CFD Forum 6 April 17, 2010 23:40
OpenFOAM on MinGW crosscompiler hosted on Linux allenzhao OpenFOAM Installation 127 January 30, 2009 20:08


All times are GMT -4. The time now is 17:44.