Hi,
I am trying to solve an eq of this type:
A(d phi/d t) = div (D grad phi) + q
in cylindrical coordinates. phi is the temperature. (A= 3500000, D= 05, q a
source(joule heating (which is as an hyperbola))
(BCs: Left -> fixedflux =0, Right -> FixedValue = 20.)
With 2D mesh (with 1 cell only (ny=1) the results are perfect but with 1D mesh
there are problems.
I tried simplifying the equation in various ways to find where the problem is:
1) A(d phi/d t) = q
2) A(d phi/d t) = div (D grad phi) (changing BCs to both Dirichlet type)
3) div (D grad phi) + q = 0 (Poisson)
While for eq 1 and 2 in both 1D and 2D the result are exactly the same, for eq
3 as for the main equation the result are totally different between 1D and 2D
(with 2D correct).
It's as if the diffusion is much more important than what it should be....
The behavior is strange but it seems that there are problems in the computation
of DiffusionTerm + source with 1D Cylindrical mesh otherwise I may have
misunderstood something.
Thanks for your help.
PS: ( I posted a new topic as the previous object was completly wrong)
_________________________________________________________________
Hotmail: Trusted email with powerful SPAM protection.
https://signup.live.com/signup.aspx?id=60969
