[phani45]

May be it's a simple thing , but ı'm confused and seek your help.

I was trying to solve a 1D diffusion equation:




fvScalarMatrix CEqn

(

fvm::ddt(C)-fvm::laplacian(DT1,C)

);

CEqn.solve();




and when I look at the flux using:




volVectorField C_flux("C_flux",(fvc::reconstruct(CEqn.flux())));




if (runTime.outputTime())

{

C_flux.write();

}




I'm getting non uniform distribution.




I'm getting a Gaussian distribution where as the flux is supposed to be constant. Isn't it?




Can any one help me to understand it