Solving PDE involving Gaussian fields.

  • Thread starter vivek.iitd
  • Start date
  • #1
54
0
I need to solve three coupled differential equations. The equations are as follows:

dE1/dz = f(E2,E3)
dE2/dz = f(E1,E3)
dE3/dz = f(E1,E2)

Where E1,E2,E3 represents field amplitudes. In case of plane waves these amplitudes will be constant in the transverse direction therefore i can directly solve these equation given their initial values. Although, i would like to, how should i solve these equations if i have a Gaussian distribution in the transverse direction? Do i need to solve it by creating a mesh and solving these equations for each pairs of points?

I thank you for your time.
 

Answers and Replies

  • #2
54
0
sorry, equations are as follows

dE1/dz = f(E2,E3)
dE2/dz = g(E1,E3)
dE3/dz = h(E1,E2).
 
  • #3
54
0
I am solving it numerically.
 
  • #4
320
49
Do they represent laser beams?
 
  • #5
54
0
yes they represent laser beams.
 
  • #6
54
0
Actually it is a nonlinear phenomenon (sum frequency generation). To give a physical aspect to this problem it is like this, two beams (with amplitudes E1 and E2) are incident on a nonlinear crystal and they mix up with each other to give third beam (E3) this whole phenomenon is represented by these coupled equations. The direction of propagation is 'z'. In case of plane waves the amplitude will be constant in the transverse direction (plane perpendicular to 'z'), therefore i can directly solve these equations using numerical method. Although, I want to get more realistic, so i am taking Gaussian distribution. Therefore in the transverse direction, amplitude distribution is Gaussian. I am confused as how to solve these equations in that case.
 
  • #7
320
49
Many times, the sum-difference problem is analyzed assuming plane waves.
If you want to consider Gaussian beams you must propagate them but the algorithm is more complex.
As an intermediate approach, I'd suggest assuming the gaussian beams have constant q paremeters. This isn't strictly true but simplifies the job. This way, at every radial coordinate r you have a gaussian amplitude and you apply the original differential equations at a given radial coordinate.
 
  • #8
54
0
Thank you for your reply. You are right that usually these problems are solved assuming plane waves because they are simplest to solve.
As you suggested I will assume q parameter to be constant as of now i.e. neglecting diffraction and other effects. After that i may incorporate that also using split-step method.
Also to solve for Gaussian profile, i guess i will have to discretize along the radial direction. And after that i can apply these ODEs by taking one step at a time along radial direction.

I think this is the correct way to solve these problem for Gaussian profile right?
 
  • #9
320
49
Even though it isn´t a "physically blameless" approach, I'd say it's a fair solution
 

Related Threads on Solving PDE involving Gaussian fields.

Replies
2
Views
1K
Replies
2
Views
3K
Replies
2
Views
776
Replies
14
Views
860
  • Last Post
Replies
8
Views
6K
  • Last Post
Replies
3
Views
3K
  • Last Post
Replies
5
Views
2K
  • Last Post
Replies
3
Views
695
  • Last Post
Replies
2
Views
620
Top