hi,
i am simulating wind pres
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? |
hi, there must be something wr
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.
|
hi,
thank you!i will check it
hi,
thank you!i will check it. |
hi,
when i set fixedValue pre
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:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) 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? |
hi weihong
This may be caus
hi weihong
This may be caused by the exception of divide operator. Divide by zero for instance. Junwei Su |
hi,
I am sorry, but I can n
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. |
Dear weihong
Run the a
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 |
hi,junwei
thank you,i will tr
hi,junwei
thank you,i will try it. |
Hi Weihong,
I believe your
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 |
All times are GMT -4. The time now is 18:52. |