CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Negative diffusion !!

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

LinkBack Thread Tools Search this Thread Display Modes
Old   July 8, 2009, 12:16
Default Negative diffusion !!
Jagmohan Meena
Join Date: May 2009
Posts: 30
Rep Power: 15
jmmeena is on a distinguished road

I am running a case which involves diffusion. I am getting following error:
Time = 4.68

smoothSolver: Solving for Ux, Initial residual = 7.20439e-08, Final residual = 7.20439e-08, No Iterations 0
smoothSolver: Solving for Uy, Initial residual = 6.83487e-08, Final residual = 6.83487e-08, No Iterations 0
smoothSolver: Solving for Uz, Initial residual = 1.9183e-07, Final residual = 1.9183e-07, No Iterations 0
DILUPBiCG: Solving for h, Initial residual = 9.49839e-06, Final residual = 9.49839e-06, No Iterations 0
GAMG: Solving for pd, Initial residual = 5.64058e-08, Final residual = 3.39985e-09, No Iterations 1
time step continuity errors : sum local = 4.40307e-12, global = -2.87747e-13, cumulative = -5.94087e-05
rho max/min : 0.721684 .676527
DILUPBiCG: Solving for T, Initial residual = 6.34211e-07, Final residual = 6.34211e-07, No Iterations 0
#0 Foam::error:rintStack(Foam::Ostream&) in "/usr/local/OpenFOAM/OpenFOAM-1.5.x/lib/linux64GccDPOpt/"
#1 Foam::sigFpe::sigFpeHandler(int) in "/usr/local/OpenFOAM/OpenFOAM-1.5.x/lib/linux64GccDPOpt/"
#2 ?? in "/lib/"
#3 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/usr/local/OpenFOAM/OpenFOAM-1.5.x/lib/linux64GccDPOpt/"
#4 Foam::fvMatrix<double>::solve(Foam::Istream&) in "/usr/local/OpenFOAM/OpenFOAM-1.5.x/lib/linux64GccDPOpt/"
#5 Foam::fvMatrix<double>::solve() at /usr/local/OpenFOAM/OpenFOAM-1.5.x/src/finiteVolume/lnInclude/fvMatrixSolve.C:187
#6 Foam::lduMatrix::solverPerformance Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at /usr/local/OpenFOAM/OpenFOAM-1.5.x/src/finiteVolume/lnInclude/fvMatrix.C:1279
#7 main at ~/OpenFOAM/jagmohan-1.5.x/applications/solvers/Scalar_Transport_combinations/buoyantTransportSimpleFoamNew/buoyantTransportSimpleFoamNew.C:85
#8 __libc_start_main in "/lib/"
#9 _start at /build/buildd/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116

I checked my solver and at 85 line, it is scalar transport equation for diffusion which has problem. I did post processing and found that there are negative values for diffusion which is not possible. I think negative values arise due to numerical error while solving. But once the negative values occur, simulation diverges. At outlet, I am using zeroGradient condition, where it diverged first. Please tell me how to define diffusion term strictly non negative or appropriate boundary condition? I am attaching 3 pictures at different time steps which are explaining about this!!


Last edited by jmmeena; July 8, 2009 at 14:29.
jmmeena is offline   Reply With Quote

Old   July 8, 2009, 14:26
Jagmohan Meena
Join Date: May 2009
Posts: 30
Rep Power: 15
jmmeena is on a distinguished road
here are the attachments !
Attached Images
File Type: jpg time_step7.jpg (25.3 KB, 45 views)
File Type: jpg time_step8.jpg (22.7 KB, 39 views)
File Type: jpg final_time_step.jpg (22.0 KB, 34 views)
jmmeena is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
negative element volume (CFX-10.0) CFDworker CFX 8 September 27, 2011 19:16
[blockMesh] BlockMesh FOAM warning gaottino OpenFOAM Meshing & Mesh Conversion 7 July 19, 2010 15:11
[blockMesh] Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Meshing & Mesh Conversion 10 April 2, 2007 15:00
negative values for diffusion scalar cmv Siemens 0 April 18, 2005 07:21

All times are GMT -4. The time now is 21:56.