Sorry to be asking so many questions, but I've hit yet another hurdle. This is sort of a continuation of my laser heating problem. Thanks to the mailing list I was able to obtain the correct temperature distribution, but the next part of the problem is to use that temperature distribution to calculate a thermo-elastic expansion given by



where u is a displacement vector of the sample surface. The first term looks like a diffusion term and the last term looks like a source, but the middle term has the divergence and gradient reversed from the convection case so I don't know how to represent it. Divergence of a gradient should give a scalar, but gradient of a divergence should give a vector. I think I should define u as
CellVariable(mesh=m, hasOld=True, value=[[0.], [0.]], elementshape=(2,))
And the right side as
2*a*(1+v)*T.grad 

or
 
2*a*(1+v)*T.faceGrad
But I can't find a way to write out the full equation that doesn't produce errors. How would you set up the left side?

Thanks again,
Adam


_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
  [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]

Reply via email to