RP Photonics

Software … powerful tools for your research & development!

deutsch

RP Fiber Power – Simulation and Design Software
for Fiber Optics, Amplifiers and Fiber Lasers

Example Case: Launching Light into a Single-mode Fiber

Description of the Model

We simulate the following case:

Some technical details of the simulation:

Note that while the mode solver is limited to cases with radial symmetry of the refractive index profile, them beam propagation could be calculated for arbitrary index profiles, as long as they are weakly guiding (which is essentially the case for all all-glass fibers).

Setting up the beam propagation with a few lines of script code is simple, after some parameters have been defined:

lambda := 1 um
 
; Define the refractive index profile:
n_cl := 1.45 { cladding index }
NA := 0.08
n_co := sqrt(n_cl^2 + NA^2) { core index }
r_co := 4 um { core radius }
n_f(r) := if r <= r_co then n_co else n_cl
 
; Grid parameters for beam propagation:
r_max := 30 um
N := 2^6
dr := 2 * r_max / N
z_max := 10 mm
dz := 10 um
N_z := z_max / dz
 
; Incident beam profile: Gaussian beam from a laser
w0 := 3 um { beam radius }
d := 2 um { vertical position error }
theta_in := 0 deg { angle error }
A0%(x, y) := exp(-(x^2 + (y - d)^2) / w0^2) / sqrt(0.5 * pi * w0^2)
  * expi((2pi / lambda) * x * sin(theta_in))
 
loss(x, y) := 1 * ((x^2 + y^2) / (10 um)^2)^6
 
calc 
  begin
    bp_set_grid(r_max, N, r_max, N, z_max, N_z, 2);
    bp_set_channel(lambda);
    bp_set_n('n_f(sqrt(x^2 + y^2))'); { index profile }
    bp_set_loss('loss(x, y)'); { loss profile }
    bp_set_A0('A0%(x, y)'); { initial amplitude }
  end

Results

Figure 1 shows how the amplitude distribution evolves in the fiber. We see some wiggles, which result from a vertical offset of the incident beam position. Also, one can see how some of the light escapes from the core; this is the non-guided part which is lost after some length of fiber. The red curve shows how the power in the fiber drops. The final value is in good agreement with the analytical estimate.

intensity profiles at the end of a multimode fiber
Figure 1: Field amplitude profile in the yz plane, showing how the radiation evolves in the fiber for a slightly offset input beam.

We can also systematically vary the initial beam radius and see how the launch efficiency varies. This is shown by Figure 2. Here, the beam is assumed to be perfectly aligned to the core.

launch efficiency as a function of the initial beam radius
Figure 2: Launch efficiency as a function of the initial beam radius.

The numerical beam propagation results nicely agree with the analytical calculation. Note that they could also easily be carried out in more complicated cases. For example, one might have a fiber core which has no radial symmetry, any other kind of input beam profile, or a bent fiber.

(back to the list of example cases)

arrow