Hi, Alshikh,
from what I understand, you shall introduce another callback function
like your eps_bragg(const vec &pt), which will return the value of
the Kerr coefficient. It should be called in a similar manner as e.g.
the susceptibilities are.
I always found it useful to test out the Kerr nonlinearity in a
homogeneous medium first; when I got an obvious anharmonic wave
distortion after a while, I could be sure the nonlinearity is well
defined.
Hope this helps a bit,
Filip
2015-06-15 16:02 GMT+02:00, alsh...@netcourrier.com <alsh...@netcourrier.com>:
>
> Hello,
>
>
> I have a bragg structure like the one defined in the
> example bragg_transmission.cpp
>
>
>
> double eps_bragg(const vec &pt) {
>
> double z = pt.z() - zsize * 0.5;
>
>
>
> if (fabs(z)*2 > Nperiods)
>
> return nlo*nlo;
>
> else {
>
> double zi;
>
> double zf = modf(z, &zi);
>
> if (zf < 0) zf += 1;
>
> if (zf < wlo)
>
> return (nlo*nlo);
>
> else
>
> return (nhi*nhi);
>
> }
>
> }
>
>
>
> how can I add (kerr) nonlinear property to the material which have high
> refrective index? so I can do s.chi3_set(value).......
>
>
>
> Thank you .....
>
>
>
>
>
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
