EDITED
Hello,
I am a beginning user of COMSOL.
Currently, I am trying to create a preliminary model of Marangoni convection, like this one available on the COMSOL site: www.comsol.com/model/marangoni-convection-110
The main difference is that mine is axisymmetric, rather than two-dimenional, and instead of a temperature difference, the fluid is heated with a point heat source on the upper, free surface.
I was planning on doing both a stationary and time dependent study, but both are running into errors.
EDIT: The problem with the strange temperature distribution has been solved. The heat conductivity was improperly defined. However, now it reports errors in user defined functions.
The error message is as follows:
Error in user-defined function.
- Function: rho
Failed to evaluate variable.
- Variable: comp1.ht.rho
- Defined as: ((unit_kg_cf*rho((comp1.T/unit_K_cf)))/(unit_m_cf^3))
Failed to evaluate expression.
- Expression: d((-2*comp1.ht.rho*comp1.ht.Cp*(comp1.ht.ur*comp1.Tr+comp1.ht.uz*comp1.Tz)*test(comp1.T)*pi*r)*(dvol),{test@15})
- Feature: Stationary Solver 1 (sol2/s1)
Could someone help me get this right?
I've attached my model to this post.
With regards,
Aaron Ayal
Hello,
I am a beginning user of COMSOL.
Currently, I am trying to create a preliminary model of Marangoni convection, like this one available on the COMSOL site: www.comsol.com/model/marangoni-convection-110
The main difference is that mine is axisymmetric, rather than two-dimenional, and instead of a temperature difference, the fluid is heated with a point heat source on the upper, free surface.
I was planning on doing both a stationary and time dependent study, but both are running into errors.
EDIT: The problem with the strange temperature distribution has been solved. The heat conductivity was improperly defined. However, now it reports errors in user defined functions.
The error message is as follows:
Error in user-defined function.
- Function: rho
Failed to evaluate variable.
- Variable: comp1.ht.rho
- Defined as: ((unit_kg_cf*rho((comp1.T/unit_K_cf)))/(unit_m_cf^3))
Failed to evaluate expression.
- Expression: d((-2*comp1.ht.rho*comp1.ht.Cp*(comp1.ht.ur*comp1.Tr+comp1.ht.uz*comp1.Tz)*test(comp1.T)*pi*r)*(dvol),{test@15})
- Feature: Stationary Solver 1 (sol2/s1)
Could someone help me get this right?
I've attached my model to this post.
With regards,
Aaron Ayal