[ayushrai]

Dear Foamers,

I am very new to OpenFOAM and hence, kindly pardon me if my doubt appears too trivial and dumb. I had read couple of previous OpenFOAM threads about implementation of temperature in simpleFoam, but none of them worked in my case. I am using OpenFOAM-2.4.0-MNF version.

Problem Description : I am simulating a case in which a cold fluid (293K) flows inside a pipe, whose wall is at constant temperature (353K). The pipe is linearly converging in nature with inlet diameter of 10 cms and outlet diameter of 6 cms. The fluid enters the pipe at 5 m/seconds. The flow is incompressible, steady-state and laminar in nature due to which I chose simpleFoam to model this problem.

Since, there are only pressure and velocity fields in simpleFoam, I decided to add temperature equation to it and name it as thermalSimpleFoam. In case of geometry, I created a wedge with small angle of 4 degrees and sufficiently long length (such that it exceeds entrance length). Moreover, I also modified one of the boundary conditions, totalTemperature, and renamed it to convectiveHeatTransfer to implement boundary condition for the Wall.

Issue : My simulation ran successfully without any issue but the results are not satisfying. Analytically in such cases, I expect a parabolic profile for temperature variation from center towards the wall, but I got very weird results.

I have attached thermalSimpleFoam solver codes, case files with time directories (for usage by ParaView to visualize the flow) and boundary condition, convectiveHeatTransfer files as well for your kind reference. The files can be downloaded from the link mentioned below:
https://we.tl/t-oKMPnAS3R1

I, therefore, request you to please take a look at those files and settings, and let me know possible mistakes and corresponding rectifications. I would be highly obliged for your act of kindness.

Thanks