Simulation of third order nonlinearities such as four wave mixing in silicon waveguides

Hello;

I myself think that simulation of nonlinearities is possible by COMSOL. However, I have not tried it until now.
The question in my mind is possibility of launching several fields (several wavelength) to a waveguide and simulating this case for propagation and interaction of them. What's your view? How can I do this?

Thanks for your help.
Sincerely,
Uldouz