amir
(amir)
#1
Dear friends,

I have a problem in coding to solve a system of thermoelastic equations.

I have written for first equation in my code as follows.

problem thel1([u,v,T],[uh,vh,Tbar])= int2d(TD)(lambda*div(u,v)*div(uh,vh)+2.*mu*(epsilon(u,v)'*epsilon(uh,vh))*

-kappa(T-Told)*div(uh,vh))- int2d(TD)(f*vh) + on(4,u=0,v=0);

However, after running, it shows me the following error!

36 : -kappa*(T-Told) error operator - <10LinearCombI7MGauche4C_F0E>, <St4pairIP6FEbaseId5v_fesEiE>

Could you please help me??

1 Like

zhaog6
#2
Note: the `-kappa * Told * div(uh, vh)`

is a right-hand side. the weak form can be written as follows

```
problem thel1([u,v,T],[uh,vh,Tbar])= int2d(TD)(lambda*div(u,v)*div(uh,vh)
+ 2.*mu*(epsilon(u,v)'*epsilon(uh,vh)) - kappa*T*div(uh,vh))
- int2d(TD)(f*vh - kappa*Told*div(uh, vh))
+ on(4,u=0,v=0);
```

1 Like

amir
(amir)
#3
Thank you very much.

Would you please help me about the second equation??

I’ve written it as follows.

problem thel2([u,v,T],[uh,vh,Tbar])= int2d(TD)(rho*c*idt*T*Tbar)

+int2d(TD)(kappa*T0*idt*div(u,v))*

+int2d(TD)(kgrad(T)'*grad(Tbar))*

-int2d(TD)(rhoc*idt*Told*Tbar+kappa*T0*idt*div(uold,vold)*Tbar)

+on(1,T=100);

I’m learning Free FEM, then I’m not ok how to code.