Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Simple heat exchanger model error

Please login with a confirmed email address before reporting spam

Hi everyone,

I'm working on a heat transfer problem and could use some help. It concerns a large horizontal pipe containing an elongated u-shaped pipe (the 'loop'). The pipe contains air and the loop contains water. The water inside the u-shaped pipe is flowing and is laminar. The air in the larger pipe is initially at rest. There is heat being produced inside the u-shape. So heat is being transported from the loop through the pipe by convection and out of the model (using the external natural convection option inside the heat flux node). I'm interested in the steady state temperature in both pipe and u-shaped pipe.

Initially I started working on a simpler model which only contained 2 cylinders (ie. I used a cylinder instead of the loop). This model worked perfectly, but when I switched the smaller cylinder for the loop I started getting errors.

1. The direct solver gave the error of a singular matrix
2. the iterative solver returned a error message saying it ran into an undefined value.

I was wondering if these errors can be caused by poor meshing, since my boundary conditions are exactly the same for the working model and the non working model? Any help on solving my problem would be greatly appreciated!

5 Replies Last Post 2017年5月31日 GMT-4 07:18

Please login with a confirmed email address before reporting spam

Posted: 9 years ago 2015年11月11日 GMT-5 09:33
In my experience, and without knowing anything else about your model, these errors result from incompletely or incorrectly specified boundary conditions.

--
Steven Conrad, MD PhD
LSU Health
In my experience, and without knowing anything else about your model, these errors result from incompletely or incorrectly specified boundary conditions. -- Steven Conrad, MD PhD LSU Health

Please login with a confirmed email address before reporting spam

Posted: 9 years ago 2015年11月11日 GMT-5 09:57
Hi Steve,

Thanks for your reply. That's exactly why I'm puzzled. The only difference between the working and non working model is the path of the loop inside the pipe. The BC's are identically for both models, the only thing that is (slightly) different is the geometry.

kind regards
Jurriaan
Hi Steve, Thanks for your reply. That's exactly why I'm puzzled. The only difference between the working and non working model is the path of the loop inside the pipe. The BC's are identically for both models, the only thing that is (slightly) different is the geometry. kind regards Jurriaan


Please login with a confirmed email address before reporting spam

Posted: 9 years ago 2015年11月11日 GMT-5 11:18
I changed to a direct solver and received more information on the error:

Division by zero.
- Function: /
Failed to evaluate variable.
- Variable: comp1.ht.rho
- Defined as: (Mw_a*comp1.spf.pA)/(comp1.ht.fluid1.minput_temperature*R_const)
Failed to evaluate expression.
- Expression: d((-comp1.ht.rho*comp1.ht.Cp*(comp1.ht.ux*comp1.Tx+comp1.ht.uy*comp1.Ty+comp1.ht.uz*comp1.Tz)*test(comp1.T)*comp1.ht.d)*(dvol),{test@3})
- Feature: Stationary Solver 1 (sol1/s1)

It seems that either comp1.ht.fluid1.minput_temperature or R_const is hitting zero. Don't have time to fully investigate why at this time, but hope that provides some clues to move forward.

--
Steven Conrad, MD PhD
LSU Health
I changed to a direct solver and received more information on the error: Division by zero. - Function: / Failed to evaluate variable. - Variable: comp1.ht.rho - Defined as: (Mw_a*comp1.spf.pA)/(comp1.ht.fluid1.minput_temperature*R_const) Failed to evaluate expression. - Expression: d((-comp1.ht.rho*comp1.ht.Cp*(comp1.ht.ux*comp1.Tx+comp1.ht.uy*comp1.Ty+comp1.ht.uz*comp1.Tz)*test(comp1.T)*comp1.ht.d)*(dvol),{test@3}) - Feature: Stationary Solver 1 (sol1/s1) It seems that either comp1.ht.fluid1.minput_temperature or R_const is hitting zero. Don't have time to fully investigate why at this time, but hope that provides some clues to move forward. -- Steven Conrad, MD PhD LSU Health

Please login with a confirmed email address before reporting spam

Posted: 7 years ago 2017年5月30日 GMT-4 13:36
Hi,

I am having a similar error. Did you solve it? How?

Thank you!

José
Hi, I am having a similar error. Did you solve it? How? Thank you! José

Please login with a confirmed email address before reporting spam

Posted: 7 years ago 2017年5月31日 GMT-4 07:18
Hi José,

Yep, turned out my boundary conditions were ill-defined. It was necessary to add a pressure reference point to the model (laminar flow>point>pressure point constraint). This solved my problem.

Hope this helps!
Hi José, Yep, turned out my boundary conditions were ill-defined. It was necessary to add a pressure reference point to the model (laminar flow>point>pressure point constraint). This solved my problem. Hope this helps!

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.