- #1

DianeLR

- 7

- 0

Hello,

I want to model the thermal behaviour of a moving heat transfer fluid in 1D, with convective exchanges with the walls. I have obtained the following equation (1 on the figure). I have performed a second order spatial discretization with decentred schemes at the extremities (y = 0 and H). After spatial discretisation, equations (2 to 4) are obtained.

By scoring these equations in OpenModelica (a software with a DASSL time integrator), I obtain consistent results at the extremities but not at the centre. I think this is due to the discretization, especially the mcp term.

Do you have any idea how to correct this problem?

I want to model the thermal behaviour of a moving heat transfer fluid in 1D, with convective exchanges with the walls. I have obtained the following equation (1 on the figure). I have performed a second order spatial discretization with decentred schemes at the extremities (y = 0 and H). After spatial discretisation, equations (2 to 4) are obtained.

By scoring these equations in OpenModelica (a software with a DASSL time integrator), I obtain consistent results at the extremities but not at the centre. I think this is due to the discretization, especially the mcp term.

Do you have any idea how to correct this problem?