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

Negative diffusion !!

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 8, 2009, 11:16
Default Negative diffusion !!
  #1
Member
 
Jagmohan Meena
Join Date: May 2009
Posts: 30
Rep Power: 8
jmmeena is on a distinguished road
Hello,

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/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/usr/local/OpenFOAM/OpenFOAM-1.5.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib/libc.so.6"
#3 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/usr/local/OpenFOAM/OpenFOAM-1.5.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 Foam::fvMatrix<double>::solve(Foam::Istream&) in "/usr/local/OpenFOAM/OpenFOAM-1.5.x/lib/linux64GccDPOpt/libfiniteVolume.so"
#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/libc.so.6"
#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!!

Regard,
JM

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

Old   July 8, 2009, 13:26
Default
  #2
Member
 
Jagmohan Meena
Join Date: May 2009
Posts: 30
Rep Power: 8
jmmeena is on a distinguished road
here are the attachments !
Attached Images
File Type: jpg time_step7.jpg (25.3 KB, 32 views)
File Type: jpg time_step8.jpg (22.7 KB, 32 views)
File Type: jpg final_time_step.jpg (22.0 KB, 29 views)
jmmeena 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
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20
negative element volume (CFX-10.0) CFDworker CFX 8 September 27, 2011 18:16
BlockMesh FOAM warning gaottino OpenFOAM Native Meshers: blockMesh 7 July 19, 2010 14:11
Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Native Meshers: blockMesh 10 April 2, 2007 14:00
negative values for diffusion scalar cmv CD-adapco 0 April 18, 2005 06:21


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