Nonlinear Wave Equation (Nonlinear Helmholtz)

  • #1
1
0

Summary:

Looking to solve a nonlinear wave equation using some sort of numerical method. Struggling with convergence issues.
I am trying to solve a PDE (which I believe can be approximated as an ODE). I have tried to solve it using 4th Order Runge-Kutta in MATLAB, but have struggled with convergence, even at an extremely high number of steps (N=100,000,000). The PDE is:
[tex] \frac{\partial^2 E(z)}{\partial z^2} + \frac{\omega^2}{c^2}\epsilon(z)E(z)=0 [/tex]
Where [itex] \epsilon(z) = \epsilon[ 1 + i \gamma |E(z)|^2] [/itex]. This is the "nonlinear" part where the imaginary part of the permittivity ([itex] \epsilon [/itex]) is dependent on the intensity of the field ([itex] ~ |E(z)|^2 [/itex]). The coefficient [itex] \gamma [/itex] can be considered the Loss Coefficient. With low [itex] \gamma [/itex] the intensity dependent term is negligible and the RK4 approach will converge, but with high intensities or higher [itex] \gamma [/itex] then it blows up.

I am using the numerical outputs from another simulation to give the field and the value of it's derivative for initial conditions.

Any suggestions on how to solve this type of differential equation? Analytically or numerically? I do not have much experience in numerical methods.
 

Attachments

Answers and Replies

  • #2
253
117
Your d.e. describes non-linear optics. Please see https://en.wikipedia.org/wiki/Soliton_(optics) for the closed form derivation that results in a soliton. You will have to account for absorption due to your complex index of refraction.
 
  • Like
Likes jim mcnamara

Related Threads on Nonlinear Wave Equation (Nonlinear Helmholtz)

Replies
1
Views
6K
Replies
0
Views
2K
Replies
5
Views
5K
Replies
1
Views
4K
Replies
4
Views
3K
Replies
11
Views
2K
Replies
0
Views
2K
Replies
3
Views
3K
  • Last Post
Replies
13
Views
1K
Replies
5
Views
951
Top