[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