On Tue, May 23, 2017 at 8:05 AM, Sergio Manzetti
<sergio.manze...@fjordforsk.no> wrote:
>
> Dear Daniel. In this example I was trying to find out a solution for the
> complex differential equation:
>
> i dphi/dt + d²phi/dx^2 + |phi^2|phi = 0
>
> However, I just learned from the fipy ring that complex conjugates are not
> possible to include in Fipy.
Right, FiPy knows nothing about complex numbers.
> However, if I should calculate the given eqn without the complex conjugate:
>
> dphi/dt + d²phi/dx^2 + |phi^2|phi = 0
That is a strange equation because it induces counter diffusion, which
probably won't solve. I'm not sure whether the source term does
something to counter that, but you can try.
> which example should I follow in the manuscript?
There isn't an example that matches this equation exactly, but try
http://www.ctcms.nist.gov/fipy/examples/diffusion/generated/examples.diffusion.mesh1D.html
for the TransientTerm and DiffustionTerm and
http://www.ctcms.nist.gov/fipy/examples/phase/generated/examples.phase.simple.html
for the ImplicitSourceTerm.
> I am not sure what the script does, when one sets a phi value before the
> given PDE...when I thought that the phi value was found exactly by FipY?
Setting the value before solving the PDEs provide the initial value if
it's an initial value problem. Even a non-initial value problem needs
to be start somewhere for the non-linear numerical iteration.
--
Daniel Wheeler
_______________________________________________
fipy mailing list
fipy@nist.gov
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
