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

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

 May 7, 2013, 12:00 how to write a scalar value to a file ???? #1 Member   Samer Join Date: Jan 2013 Posts: 31 Rep Power: 4 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

 May 7, 2013, 13:19 #2 Senior Member   Niels Gjoel Jacobsen Join Date: Mar 2009 Location: Rotterdam, The Netherlands Posts: 1,595 Rep Power: 24 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: `/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

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post elvis OpenFOAM Native Meshers: snappyHexMesh and Others 175 March 20, 2015 12:32 immortality OpenFOAM Running, Solving & CFD 12 May 3, 2013 18:21 mcathela OpenFOAM Installation 14 April 23, 2013 13:59 Ardalan Main CFD Forum 6 April 17, 2010 23:40 allenzhao OpenFOAM Installation 127 January 30, 2009 20:08

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