krystian, how did you solve
 

krystian,

how did you solve this error:

#0 Foam::error::printStack(Foam:stream&)
#1 Foam::sigFpe::sigFpeHandler(int)
#2 __restore_rt
#3 Foam::momentCoefficient(Foam::GeometricField<foam: :vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::dimensioned<double> const&, Foam::word const&, Foam::Vector<double> const&, double const&, double const&)
#4 main
#5 __libc_start_main
#6 __gxx_personality_v0
Floating point exception

i am getting the same message using v1.4 on fedora for the cavity example ....
am looking now in the forum for any clues but would really appreciate it if you could explain how you solved it

cheers

Nalini, It would be good if
 

Nalini,

It would be good if you could also include what appears before these lines which you posted.

Generally, the answer for that is, that liftDrag is not suited for every solver. Current stage only allows its use for one-phase flows or at least it doesn't work for two-phase flows which I tested.

Try to run it for simpleFoam, that should work. However for icoFoam it doesn't.

As you see it is kind of limited. If anyone managed to run it for rasInterFoam, any info will be appreciated.

Thanks,
K

thanks i was using icoFoam w
 

thanks
i was using icoFoam
will have a play and see what i can come up with

let us know what you've got.
 

let us know what you've got.

HI All, I am getting follow
 

HI All,

I am getting following error at the last step http://www.cfd-online.com/OpenFOAM_D...lipart/sad.gif
of
Any clues on what is going on here.
thanks in advance

g++ -m32 -Dlinux -DDP -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -O3 -DNoRepository -ftemplate-depth-40 -I/home/raj/OpenFOAM/OpenFOAM-1.4.1/src/postProcessing/incompressible/lnInclude -I/home/raj/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude -I/home/raj/OpenFOAM/OpenFOAM-1.4.1/src/transportModels -I/home/raj/OpenFOAM/OpenFOAM-1.4.1/src/turbulenceModels -I/home/raj/OpenFOAM/OpenFOAM-1.4.1/src/LESmodels -I/home/raj/OpenFOAM/OpenFOAM-1.4.1/src/LESmodels/LESdeltas/lnInclude -IlnInclude -I. -I/home/raj/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude -fPIC -pthread Make/linuxGccDPOpt/liftDrag.o -L/home/raj/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt \
-lincompressiblePostProcessing -lfoamUtil -lOpenFOAM -ldl -lm -o /home/raj/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/liftDrag
/usr/bin/ld: cannot find -lfoamUtil
collect2: ld returned 1 exit status
make: *** [/home/raj/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/liftDrag] Error 1

Hi I have a dreadfully slow c
 

Hi
I have a dreadfully slow connection and downloading OF1.2 will take 4 hrs or more. Can anyone just give me the liftdrag sources in a tarball??

Thanks.

You can find all the files acc
 

You can find all the files accessible individually in the

dev SVN repository.

Try:

src/postProcessing/incompressible/liftDrag/

Enjoy,

Hrv

Hello Forum, I just instal
 

Hello Forum,

I just installed the utility liftDrag as suggested in:

http://www.cfd-online.com/OpenFOAM_D...es/1/2299.html

I am running a very simple case in simpleFoam. I have a naca0012 at 0 angle of attack, and k-epsilon for modeling turbulence. I get the analysis to converge nicely, pressure looks good, velocity, k, epsolin, etc. look very good. When I run liftDrag . naca0012 -time 700 for example, I get this output:

[gtg627eOpenFOAM@ruzzene03 simpleFoam]$ liftDrag . naca0012 -time 700
/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4.1 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : liftDrag . naca0012 -time 700
Date : Oct 29 2007
Time : 03:42:42
Host : ruzzene03
PID : 9647
Root : /home/gtg627eOpenFOAM/OpenFOAM/gtg627eOpenFOAM-1.4.1/run/tutorials/simpleFoam
Case : naca0012
Nprocs : 1
Create time

Create mesh for time = 700

Time = 700
Reading U

Reading p

Inlet velocity: (50 0 0)
Wall patch 1 named airfoiln :
Reference area: 0.297421 Reference length: 1 Drag coefficient: 0.00116109 Lift coefficient: (0 8.21932e-05 -1.32776e-21) Moment coefficient: (-2.05483e-06 2.90273e-05 -1.95907e-05)

This result seems to not include viscous forces. The drag for a naca0012 at 0 deg and Re ~ 3 million should be approximately 0.0058, from "theory of wing sections."

Thank you in advance for any comments,

Alessandro

Hello Forum, In directory "
 

Hello Forum,

In directory "applications/utilities/postProcessing/miscellaneous/liftDrag/liftDrag.C

I see the follwing call:

-------------------------------------------------------------

Info<< "Wall patch " << patchI
<< " named " << patches[patchI].name() << " : " << nl
<< "Reference area: " << Aref
<< " Reference length: " << Lref
<< " Drag coefficient: "
<< dragCoefficient
(
U,
p,
nu,
patches[patchI].name(),
Uav,
Aref

)

------------------------------------------------------------------

but in file

../src/postProcessing/incompressible/liftDrag/liftDrag.C

I see that I can also compute turbDragCoefficient.

Following the above example for the drag coefficient I can call:


scalar turbDragCoefficient
(
const autoPtr<foam::turbulencemodel>& turbulence,
const volVectorField& U,
const volScalarField& p,
const dimensionedScalar& mu,
const word& patchName,
const vector& Uinf,
const scalar& Aref
);

My questions are:

1. How do I specify the input const autoPtr<foam::turbulencemodel>& turbulence?

2. is mu the same as nu? It seems so in the first file I mentioned above.

Thank you again,

Alessandro

Hello Forum, Sorry for the
 

Hello Forum,

Sorry for the impulsive message writing. The solution to the problem I pointed out in the previous two messages can be found in:

http://www.cfd-online.com/OpenFOAM_D...es/1/5067.html

Thank you,

Alessandro

Hi Forum, i am a new OpenFO
 

Hi Forum,
i am a new OpenFOAM user and have a problem with compieling the lift and drag force tool.
At step 6 II the following error showed up, and an error occurs if i try to run step 7.
Would be perfect if somebody would help me.
The following lines is only the end of the error.


/tmp/cchPMiak.s:11202: Error: unknown pseudo-op: `.weakref'
/tmp/cchPMiak.s:11203: Error: unknown pseudo-op: `.weakref'
/tmp/cchPMiak.s:11204: Error: unknown pseudo-op: `.weakref'
/tmp/cchPMiak.s:11205: Error: unknown pseudo-op: `.weakref'
/tmp/cchPMiak.s:11206: Error: unknown pseudo-op: `.weakref'
/tmp/cchPMiak.s:11207: Error: unknown pseudo-op: `.weakref'
/tmp/cchPMiak.s:11208: Error: unknown pseudo-op: `.weakref'
/tmp/cchPMiak.s:11209: Error: unknown pseudo-op: `.weakref'
/tmp/cchPMiak.s:11210: Error: unknown pseudo-op: `.weakref'
/tmp/cchPMiak.s:11211: Error: unknown pseudo-op: `.weakref'
/tmp/cchPMiak.s:11212: Error: unknown pseudo-op: `.weakref'
/tmp/cchPMiak.s:11213: Error: unknown pseudo-op: `.weakref'
/tmp/cchPMiak.s:11214: Error: unknown pseudo-op: `.weakref'
/tmp/cchPMiak.s:11215: Error: unknown pseudo-op: `.weakref'
/tmp/cchPMiak.s:11917: Warning: setting incorrect section attributes for .data.DW.ref.__gxx_personality_v0
/tmp/cchPMiak.s:11924: Warning: setting incorrect section attributes for .text.__i686.get_pc_thunk.bx
make: *** [Make/linuxGccDPOpt/liftDrag.o] Error 1

Thanks a lot.

This is a compiler problem - y
 

This is a compiler problem - you need to rebuild gcc on your machine and it will go away.

Hrv

It may not be that bad. This w
 

It may not be that bad. This weakref thing is usually a problem with the assembler. Usually this means that the binutils (where the assembler belongs) are waaaay older than the compiler. So recompiling the binutils might help as well (it is as difficult as recompiling gcc, but takes less time)

Bernhard

Hi again, thank you for th
 

Hi again,

thank you for the help, now step 6 compiled well, but i still get an error during step 7.
Sorry for asking proberbly stupid questions buti am quiet new with linux and OpenFOAM

Here is the error massage.
ccDPOpt/liftDrag.o -L/home/fabian/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt \
-lincompressiblePostProcessing -lfoamUtil -lOpenFOAM -ldl -lm -o /home/fabian/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/liftDrag
/usr/bin/ld: cannot find -lfoamUtil
collect2: ld returned 1 exit status
make: *** [/home/fabian/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/liftDrag] Error 1

Is it still a compiler problem? I compiled gcc-4.2.1 on my computer:

Thanks for helping me again.

Everything what i tried before
 

Everything what i tried before was executed with OpenFOAM 1.4.1,
When i try it with OpenFOAM 1.4, then i get the same massage exept that it cannot find -lincompressiblePostProcessing

Could anyone help me with my l
 

Could anyone help me with my liftDrag to 1.4.1 conversion?

I followed the directions posted in this thread until I ran in to the same problem as Fabian and Rajneesh ie. "foamUtil not found".

I then copied the foamUtil directory from 1.2 and ran "wmake libso foamUtil".

After that liftDrag compiled without errors but when I test it on the icoFoam/cavity:

/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4.1 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : liftDrag . cavity -latestTime
Date : Feb 26 2008
Time : 16:37:31
Host : LinuxLoota
PID : 27737
Root : /home/juho/OpenFOAM/peltol20-1.4.1/run/tutorials/icoFoam
Case : cavity
Nprocs : 1
Create time

Create mesh for time = 0.5

Time = 0.5
Reading U

Reading p

Inlet velocity: (1 0 0)
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/home/juho/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/juho/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xffffe420]
#3 Foam::momentCoefficient(Foam::GeometricField<foam: :vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::dimensioned<double> const&, Foam::word const&, Foam::Vector<double> const&, double const&, double const&) in "/home/juho/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libincompressiblePostProce ssing.so"
#4 main in "/home/juho/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/liftDrag"
#5 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#6 Foam::regIOobject::readIfModified() in "/home/juho/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/liftDrag"
Floating point exception (core dumped)

==================================================


It gives a similar error message when I try it on my automotive test case but in this case it also lists the inlet velocity as 60 instead of 30 it really is.

I'm using a 32-bit Ubuntu Gutsy on a Intel machine.

Any guidance is greatly appreciated!

Juho

There is a division by zero so
 

There is a division by zero somewhere. Please run the liftDrag utility on a different case to confirm this. Alternately, you can find several threads in this forum that help you to modify your solver to print lift and drag info directly in the log file for every time step.

Thank you. I ran it on my s
 

Thank you.

I ran it on my simpleFoam automotive test case with a similar result.

I found your (?) turbFoam_1 solver with built in lift and drag calculation on the forum. I'll look into it next.

I included the lift and drag c
 

I included the lift and drag calculations to the simpleFoam but in a very crude way. I compiled in the reference velocity and kinematic viscosity.

Seems to be working...

Hello!! I'm running some ca
 

Hello!!

I'm running some cases of the tutorial... I've found some differences between my results and the results that appear in the User Guide (v.OF-1.4). Specially in the cavityHighRe case. (in spite of we're using the same set-up) The Guide says that U should converge at Time=1.63, but when I run it, it converged at Time=1.805. Then we run it in other machine and other version of OpenFoam (1.4.1) and the results where the same of v.1.4 (Time=1.805), but different to the Manual: Our results at Time=1.63 are:


>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

Time = 1.63

Courant Number mean: 0.11671 max: 0.839843
DILUPBiCG: Solving for Ux, Initial residual = 2.18463e-05, Final residual = 2.98021e-07, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 2.52903e-05, Final residual = 2.92166e-07, No Iterations 1
DICPCG: Solving for p, Initial residual = 2.94105e-05, Final residual = 8.12416e-07, No Iterations 8
time step continuity errors : sum local = 6.12008e-09, global = 4.01513e-19, cumulative = 1.31549e-18
DICPCG: Solving for p, Initial residual = 1.08552e-05, Final residual = 8.55134e-07, No Iterations 5
time step continuity errors : sum local = 6.10618e-09, global = -6.73987e-19, cumulative = 6.41498e-19
ExecutionTime = 0.84 s ClockTime = 1 s

Time = 1.635

Courant Number mean: 0.116714 max: 0.839846
DILUPBiCG: Solving for Ux, Initial residual = 2.12558e-05, Final residual = 2.96012e-07, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 2.46569e-05, Final residual = 2.84806e-07, No Iterations 1
DICPCG: Solving for p, Initial residual = 2.85749e-05, Final residual = 9.04189e-07, No Iterations 9
time step continuity errors : sum local = 9.14192e-09, global = -3.49815e-19, cumulative = 2.91684e-19
DICPCG: Solving for p, Initial residual = 1.07851e-05, Final residual = 9.91138e-07, No Iterations 6
time step continuity errors : sum local = 8.95587e-09, global = -4.25171e-20, cumulative = 2.49167e-19
ExecutionTime = 0.84 s ClockTime = 1 s

<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< <


THE QUESTIONS ARE: WHY ARE THESE DIFFERENCES??
has this differences happened to someone else???
Is this because of the OpenFoam Version????
Which is the OpenFOAM version used in the icoFoam Tutorial??????


I will appreciate any comment or suggestion !!! ;-) .... Thanks in advance

Gaby