Как мне представить диффузный источник, который включает координаты?

Question

C (h) \ frac {\ частичный h} {\ частичный t} = \ nabla \ cdot (K (h) \ nabla (h + z)) У меня есть два вопроса:

1. Как мне представить диффузный источник, который включает координаты? Основное уравнение для h;k (h) является функцией h;z была координатой z. Я знал, как написать диффузный источник без координаты z, но я не знал, как представить диффузионный источник, включающий координату.
2. Сделать "phi [0] .updateOld ()" в "... \ fipy \ examples \ diffusion \ mesh1D.py" обновить "eq = transientTerm () == DiffusionTerm (coeff = D0 * (1 -phi [0])) "когда просто использовать" res = eq.sweep (var = phi [0], dt = timeStepDuration) "

Спасибо за вашу помощь.

Answer 1

1. 

и

т.

Второй термин в RHS записывается в явном виде как (K(h).faceValue * [[0], [0], [1]]).divergence.

Если K(h) является линейной функцией h, то вы можете выразить это как ConvectionTerm. Например, если K(h) = h * K1(h), то вы могли бы написать ConvectionTerm(coeff=K1(h).faceValue * [[0], [0], [1]]).

Вы имеете в виду, обновляется ли coeff=D0 * (1 - phi[0])? Да.