Ivar KJELBERG
COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)
Please login with a confirmed email address before reporting spam
Posted:
9 years ago
31 août 2015, 05:34 UTC−4
Hi
there are a few possibilities why you get unexpected results in these types of conjugated heat transfer.
1) what often happens, one miss a BC between the SPF and the HT section, you need to define your SPF inflow and outflow, as well as the HT outflow, and Temperature constant at the inflow area, as well as at the side boundary of your tube
2) most common, the meshing for the very anisotropic shape of your long and narrow tube. You need to be sure your mesh resolves the gradient of both T and the velocity along the wall. If you check the oil heat diffusivity between 20 and 150°C you will see its harly more than some 0.06-0.09 mm^2/s and the flow velocity is betwen 0.3-1.5 m/s. Oil is viscous, hence the no-slip layer is quite thick and the heat will only propagate here by conduction. Resulting all in the requirement of a high density mesh along the warm wall, as your heat and flow velocity changes only in the first 3-5 mm along the wall
I would propose you start with a 2D-axi model, with structured mapped mesh, with a distribution along the axis of some 200 elements, and across (radial) some 15 elements and a ratio of 14 with a geometric sequence
Then you need to check that this mesh is sufficient, by gently densifying/coarsening it and checking the effect on the results, I havent don ethat in the model below.
To help the solver you could use a Poiseuil type pressure drop as initial conditions, I see some 140[Pa/m] along the axial length at 0.3[m/s] or 0.5[kg/s] flow.
Finally you flow needs a meter to stabilise in velocity, so you could perhaps add a meter and not heat this first meter of the pipe wall.
I would also suggest that you use the stationary solver with an auxialliary Sweep and NOT the Parametry Sweep node, as for me in my v5.1.0.180 I notice the Parametric Sweep node did not detect correctly the continuation condition, and restarted the initial conditions from scratch for each parametric step, what takes twice as long to solve, this is not normal ;(
--
Good luck
Ivar
Hi
there are a few possibilities why you get unexpected results in these types of conjugated heat transfer.
1) what often happens, one miss a BC between the SPF and the HT section, you need to define your SPF inflow and outflow, as well as the HT outflow, and Temperature constant at the inflow area, as well as at the side boundary of your tube
2) most common, the meshing for the very anisotropic shape of your long and narrow tube. You need to be sure your mesh resolves the gradient of both T and the velocity along the wall. If you check the oil heat diffusivity between 20 and 150°C you will see its harly more than some 0.06-0.09 mm^2/s and the flow velocity is betwen 0.3-1.5 m/s. Oil is viscous, hence the no-slip layer is quite thick and the heat will only propagate here by conduction. Resulting all in the requirement of a high density mesh along the warm wall, as your heat and flow velocity changes only in the first 3-5 mm along the wall
I would propose you start with a 2D-axi model, with structured mapped mesh, with a distribution along the axis of some 200 elements, and across (radial) some 15 elements and a ratio of 14 with a geometric sequence
Then you need to check that this mesh is sufficient, by gently densifying/coarsening it and checking the effect on the results, I havent don ethat in the model below.
To help the solver you could use a Poiseuil type pressure drop as initial conditions, I see some 140[Pa/m] along the axial length at 0.3[m/s] or 0.5[kg/s] flow.
Finally you flow needs a meter to stabilise in velocity, so you could perhaps add a meter and not heat this first meter of the pipe wall.
I would also suggest that you use the stationary solver with an auxialliary Sweep and NOT the Parametry Sweep node, as for me in my v5.1.0.180 I notice the Parametric Sweep node did not detect correctly the continuation condition, and restarted the initial conditions from scratch for each parametric step, what takes twice as long to solve, this is not normal ;(
--
Good luck
Ivar