[babakflame]

Dear Fellows

In my solver, by solving transport equations, I am achieving which is bulk charge density and which is velocity. I have added the following snippet in my createFields.H file:

I want to find the distribution of from my code in which is the component of velocity in x direction.

I have added the following snippet in my createFields.H file:

 volScalarField Ux
        (
            IOobject
            (
                "Ux",
                runTime.timeName(),
                mesh,
                IOobject::NO_READ,
                IOobject::NO_WRITE
            ),
            U.component(vector::X)
        );

volScalarField rhoE      
    (
        IOobject
        (
            "rhoE",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT,
            IOobject::AUTO_WRITE
        ),
        mesh
//    dimensionedScalar("rhoE", dimensionSet(0, -3, 1, 0, 0, 1, 0), scalar(0.0))
    );

  volScalarField rhoEUx
    ( 
     IOobject
        (
            "rhoEUx",
             runTime.timeName(),
             mesh,
             IOobject::NO_READ,
             IOobject::AUTO_WRITE
         ),
         rhoE*Ux
     );

And I have added this line to my Ueqn.H

     rhoEUx = rhoE * Ux;
    
 //     volVectorField temptemp = rhoE*E;


//    volVectorField temptemp = 3e-06*(fvc::grad(UeExt));
    fvVectorMatrix UEqn
    (
        fvm::ddt(rho, U)
      + fvm::div(rhoPhi, U)
      - fvm::laplacian(muEff, U)
      - (fvc::grad(U) & fvc::grad(muEff))
      - tmp2
      - tmp1 
    );

    UEqn.relax();

Everything with the solver is fine, However, although the values for and are calculated correctly, the is not correct from the code,


However if I use calculator filter in paraview, I would be able to create the distribution of . Does any body know, where did I make mistake?

Regards