[Tempest]

Hello OFers,

I am trying to add a source term in the energy equation (see image) for simulating a lid driven cavity. This involves second partial derivative of Temperature.

For this, I am doing something like this:

volVectorField gradT = fvc::grad(T);
volScalarField Ty = gradT.component(1); // T_y
volVectorField gradT2 = fvc::grad(Ty); // (T_xy, T_yy, T_zy)
volScalarField Tyy = gradT2.component(1); // T_yy

Is this correct? Or what can be the correct way to do it?

[Santiago]

Originally Posted by Tempest

Hello OFers,

I am trying to add a source term in the energy equation (see image) for simulating a lid driven cavity. This involves second partial derivative of Temperature.

For this, I am doing something like this:

volVectorField gradT = fvc::grad(T);
volScalarField Ty = gradT.component(1); // T_y
volVectorField gradT2 = fvc::grad(Ty); // (T_xy, T_yy, T_zy)
volScalarField Tyy = gradT2.component(1); // T_yy

Is this correct? Or what can be the correct way to do it?

Isnt that a laplacian? OF implements a laplacian object in the fvm/fvc namespace