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

Does someone has the same problem

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 20, 2008, 21:30
Default hi, i am simulating wind pres
  #1
Senior Member
 
weihong yao
Join Date: Mar 2009
Posts: 117
Rep Power: 8
ivanyao is on a distinguished road
hi,
i am simulating wind pressure on a cube.everything is fine,but when i calculated,it appeared sth wrong. like that;
Create mesh for time = 0
> Reading field p
> Reading field U
> Reading/calculating face flux field phi
> Selecting incompressible transport model Newtonian
> Selecting turbulence model kEpsilon
> Starting time loop
> Time = 0.005
> Courant Number mean: 0 max: 0.0774444
> DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual =
2.04207e-09, No Iterations 1
> DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual =
4.95291e-09, No Iterations 1
> DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual =
4.73817e-09, No Iterations 1
> --> FOAM FATAL ERROR : Continuity error cannot be removed by
adjusting the outflow.
> Please check the velocity boundary conditions and/or run
potentialFoam to initialise the outflow.
> Total flux : 971.7
> Specified mass inflow : 50184
> Specified mass outflow : 0
> Adjustable mass outflow : 0
> From function adjustPhi(surfaceScalarField& phi, const
volVectorField& U,const volScalarField& p
> in file cfdTools/general/adjustPhi/adjustPhi.C at line 111.
> FOAM exiting
does someone has the same problem?how to fix it?
ivanyao is offline   Reply With Quote

Old   October 20, 2008, 23:58
Default hi, there must be something wr
  #2
Member
 
roy fokker
Join Date: Mar 2009
Posts: 44
Rep Power: 8
dbxmcf is on a distinguished road
hi, there must be something wrong with your boundary condition, maybe no fixedValue pressure boundary for both inlet and outlet, search the error "Continuity error cannot be removed by adjusting the outflow", I remember this has been discussed several times.
dbxmcf is offline   Reply With Quote

Old   October 21, 2008, 10:16
Default hi, thank you!i will check it
  #3
Senior Member
 
weihong yao
Join Date: Mar 2009
Posts: 117
Rep Power: 8
ivanyao is on a distinguished road
hi,
thank you!i will check it.
ivanyao is offline   Reply With Quote

Old   October 22, 2008, 21:04
Default hi, when i set fixedValue pre
  #4
Senior Member
 
weihong yao
Join Date: Mar 2009
Posts: 117
Rep Power: 8
ivanyao is on a distinguished road
hi,
when i set fixedValue pressure boundary for both inlet and outlet.and then calculate.it also appeared some errors:
Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model kEpsilon

Starting time loop

Time = 0.005

Courant Number mean: 0 max: 0.0774444
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 2.04207e-09, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 4.95291e-09, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 4.73817e-09, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 9.37246e-07, No Iterations 300
time step continuity errors : sum local = 1.70147e-10, global = 1.95661e-12, cumulative = 1.95661e-12
DICPCG: Solving for p, Initial residual = 7.43961e-05, Final residual = 9.77997e-07, No Iterations 120
time step continuity errors : sum local = 1.61161e-10, global = 1.48469e-11, cumulative = 1.68035e-11
#0 Foam::error::printStack(Foam:stream&) in "/home/ivan/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/ivan/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xb7fae420]
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/ivan/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<foam::fvpatchfield,>(Foam::GeometricF ield<double,>&, Foam::GeometricField<double,> const&, Foam::GeometricField<double,> const&) in "/home/ivan/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libincompressibleTurbulenc eModels.so"
#5 Foam::tmp<foam::geometricfield<double,> > Foam::operator/<foam::fvpatchfield,>(Foam::tmp<foam::geometricfie ld<double,> > const&, Foam::GeometricField<double,> const&) in "/home/ivan/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libincompressibleTurbulenc eModels.so"
#6 Foam::turbulenceModels::kEpsilon::correct() in "/home/ivan/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libincompressibleTurbulenc eModels.so"
#7 main in "/home/ivan/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/turbFoam"
#8 __libc_start_main in "/lib/libc.so.6"
#9 Foam::regIOobject::readIfModified() in "/home/ivan/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/turbFoam"

where did me make mistake?
ivanyao is offline   Reply With Quote

Old   October 23, 2008, 16:00
Default hi weihong This may be caus
  #5
Senior Member
 
su_junwei's Avatar
 
su junwei
Join Date: Mar 2009
Location: Xi'an China
Posts: 151
Rep Power: 10
su_junwei is on a distinguished road
Send a message via MSN to su_junwei
hi weihong

This may be caused by the exception of divide operator. Divide by zero for instance.

Junwei Su
su_junwei is offline   Reply With Quote

Old   October 26, 2008, 09:20
Default hi, I am sorry, but I can n
  #6
Senior Member
 
weihong yao
Join Date: Mar 2009
Posts: 117
Rep Power: 8
ivanyao is on a distinguished road
hi,
I am sorry, but I can not understand what's the problem exactly。would you tell me in detail.it is important to me.thank you.
ivanyao is offline   Reply With Quote

Old   October 27, 2008, 05:16
Default Dear weihong Run the a
  #7
Senior Member
 
su_junwei's Avatar
 
su junwei
Join Date: Mar 2009
Location: Xi'an China
Posts: 151
Rep Power: 10
su_junwei is on a distinguished road
Send a message via MSN to su_junwei
Dear weihong

Run the application in gdb, and using the command "bt" (backtrace) to judge where such a exception occurred.

Also, try setting the fixed value boundary for inlet and zero gradient for outlet for the pressure p.

Junwei Su
su_junwei is offline   Reply With Quote

Old   October 27, 2008, 06:14
Default hi,junwei thank you,i will tr
  #8
Senior Member
 
weihong yao
Join Date: Mar 2009
Posts: 117
Rep Power: 8
ivanyao is on a distinguished road
hi,junwei
thank you,i will try it.
ivanyao is offline   Reply With Quote

Old   October 27, 2008, 08:49
Default Hi Weihong, I believe your
  #9
Member
 
John Wang
Join Date: Mar 2009
Location: Singapore
Posts: 71
Rep Power: 8
cwang5 is on a distinguished road
Hi Weihong,

I believe your problem is caused by your boundary condition setup. Please refer to the PitzDaily tutorial in simpleFoam for the proper boundary condition for inlet and outlet.

btw, I think this post should be in the general OpenFoam forum instead of the bug forum.

Hope that helps.

John
cwang5 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
VOF diffuse problem on two fluids problem Fang Jin FLUENT 5 February 17, 2011 06:38
Problem in Modelling Heat Transfer Problem Deepak R FLUENT 1 December 6, 2007 10:37
Problem in cavity flow problem saad Main CFD Forum 4 November 1, 2007 08:45
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 06:29
problem in solving "wave generation" problem san FLUENT 2 April 3, 2006 23:37


All times are GMT -4. The time now is 02:49.