Re: inertial terms

2016-05-17 Thread Daniel Wheeler
Hi Francisco, Again, it's the same issue in FiPy, this type of term needs to be split into a convection term and a source. Cheers, Daniel On Mon, May 16, 2016 at 1:45 PM, Francisco Vega Reyes wrote: > A related question is there is away to express terms of the form: > > u(x) \nabla f(x) > > wh

Re: transient term with coefficient

2016-05-17 Thread Daniel Wheeler
Hi Francisco, In FiPy, it is necessary to add both a TransientTerm and a source term for terms such as, density (times) (\partial f[x]/\partial t) So, the code will look something like this, TransientTerm(var=var, coeff=density) - field * (density - density.old) / dt There may be ways