Hi,
Thanks for the previous help. Regarding your previous answer on those
expressions, I checked the three different methods: my original version,
and both of your suggestions, and I see that they match. Nevertheless, I am
taking your suggestion.
I have a couple of follow-ups regarding a
Sorry for the slow reply.
Expressions like [(1.0/density),] and [(temperature / density**2),]
unfortunately don't work reliably in FiPy. They might work, but I wouldn't
count on it.
The issue is that FiPy sees a Python list that holds the instantaneous value of
a rank-0 CelVariable, rather
Hi,
I've adjusted the system a little bit, and I wanted some verification and
more criticism.
I took your suggestions about (Z.grad)**2 and changed it to
numerix.dot(Z.grad, Z.grad) to give me a scalar.
In addition, the terms that could not be easily written as convection terms:
Kevin -
- Should I model the last term in the temperature (T) equations as a convection
term or explicit source term?
The best coupling would be obtained by treating this term as a convection term
on n with a velocity proportional to \partial T/\partial x. Unfortunately,
(D/n) sits outside
