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.

Weak Form of a transient heat problem in solid

Please login with a confirmed email address before reporting spam

Hi to everyone. I am writing here because I have created two files, one is built using heat transfer module and the second is built using the weak form module.

I do not understand why the results are so different.

these are the caractheristics:

Module: heat transfer in solids Geometry: 0.5m x 0.5m Mesh normal Initial temperature 293.15 left edge=293.15 Right edge=303.15 transient: time step size 0.001; from 0 to 0.1 k=cp=rho=1

weak file: dependent variabile: T Dirichlet boundary conditions: left 293.15; right 293.15 initial value:293.15 the other characteristics are the same. weak formula: -rhocp(Tttest(T)-(kTxx)test(T)-(kTyy)test(T))

Can someone help me with this issue?


I am struggling.

Best Regards


2 Replies Last Post 11.06.2018, 07:10 GMT-4

Please login with a confirmed email address before reporting spam

Posted: 6 years ago 11.06.2018, 07:06 GMT-4
Updated: 6 years ago 11.06.2018, 07:11 GMT-4

Hello,

I have two suggestions for you :

Hello, I have two suggestions for you :

Please login with a confirmed email address before reporting spam

Posted: 6 years ago 11.06.2018, 07:10 GMT-4

Hello,

It probally comes from your Weak Form. Firstly, you have a parenthesis problem. You souldn't multiply rho x Cp with k.

Secondly, in ordre to calculate the weak form, you should :

1) Multiplie your PDE by the test function

2) Intergrate the expression on the volume

3) Use the Integration By Parts and the Green Theorem to reduce the derivation degree (avoid the Txx and transform them into Tx) on the dependent variable. For the state-stady part of the heat equation, the method is given step by step in this COMSOL article. You will notice that an intergral on the boundaries appears. Do not enter it in the Weak Form PDE node, it is taken care of by the Zero Flux node.

In tested this equation and it worked :

Hope it helps !

Hello, It probally comes from your Weak Form. Firstly, you have a parenthesis problem. You souldn't multiply rho x Cp with k. Secondly, in ordre to calculate the weak form, you should : 1) Multiplie your PDE by the test function 2) Intergrate the expression on the volume 3) Use the Integration By Parts and the Green Theorem to reduce the derivation degree (avoid the Txx and transform them into Tx) on the dependent variable. For the state-stady part of the heat equation, the method is given step by step in this [COMSOL article](http://www.comsol.com/blogs/strength-weak-form/). You will notice that an intergral on the boundaries appears. Do not enter it in the Weak Form PDE node, it is taken care of by the Zero Flux node. In tested this equation and it worked : - Tt*test(T) -test(Tx)*Tx-test(Ty)*Ty-test(Tz)*Tz Hope it 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.