RP Fiber Power – Simulation and Design Software
for Optical Fibers, Amplifiers and Fiber Lasers
| Start | Purpose | Model | Fiber Data | Interface | Demos | Licenses |
Spatial Resolution
RP Fiber Power can take into account the variation of optical intensities and the excitation of laser-active ions along the fiber and also in the radial and azimuthal direction. The user can define the concentration profile of the laser-active dopant.
It is also possible to model cases where different types of ions (e.g., Er3+ and Yb3+) have different doping profiles.
Finally, optical propagation modes (see below) with different spatial profiles can be included, which then interact with the dopant.
Mode Solver
RP Fiber Power contains a powerful mode solver, which can very quickly calculate all guided modes (LP modes) for a fiber from a given radially symmetric refractive index profile. This works even for fibers with several hundred modes. The mode solver also calculates propagation constants and the fraction of power guided in the core. The resulting intensity profiles can be very easily used in the definition of optical channels for the calculation of the power evolution in the fiber.
The refractive index profile can be specified with a formula or from tabulated values as result from measurements. A wavelength dependence can also be included, which also allows one to calculate the chromatic dispersion (material and waveguide dispersion).
Level Systems and Transitions
RP Fiber Power can be used to model cases where two or even three different types of ions participate in an amplification or lasing process. Each ion may have a substantial number of metastable excited levels.
The user can freely define various kinds of transitions between electronic levels:
- spontaneous transitions (radiative or non-radiative)
- optically induced transitions (absorption or stimulated emission)
- quenching processes
- energy transfers (between ions of the same type, or of different types, e.g. between Er3+ and Yb3+ ions)
These inputs determine the rate equation system, which is then very efficiently solved with an advanced algorithm.
The figure to the right shows an example case in the context of a thulium-doped upconversion fiber laser. Here, one is dealing with three metastable levels in addition to the ground state. There are four absorbing transitions (including the quasi-three-level laser transition) and two transitions with stimulated emission (pump and signal), apart from various spontaneous transitions.
Optical Channels
RP Fiber Power allows one to take into account a large number of so-called optical channels, which may represent single propagation modes or some (possibly large) collection of propagation modes. For modeling amplified spontaneous emission (ASE), channels can have a finite optical bandwidth. For each channel, the transverse intensity distribution can be freely defined. For fiber amplifiers, the noise figure of any signal channel can be calculated.
The figure below illustrates the optical channels in an erbium-doped fiber amplifier: a pump channel, two signal channels, and 16 ASE channels.
Dynamic Simulations
The software also allows you to simulate the temporal evolution of your system. A highly flexible approach has been implemented, allowing both the precise simulation of the dynamics of lasers, fully taking into account the propagation time in the fiber, and very fast simulations of single-pass or double-pass pulse amplification.
As an example, the figure on the right side shows how a pulse shape gets distorted by gain saturation during amplification in a high-gain fiber amplifier, and how much of the stored energy is extracted by such a pulse. Another example corresponds to a Q-switched fiber laser model, which is described on a separate page.
Ultrashort Pulse Propagation
An ultrashort pulse can be numerically represented with an array of complex amplitudes in the time or frequency domain. Fourier transforms are used to relate time and frequency to each other. Wavelength-dependent amplification and chromatic dispersion are modeled in the frequency domain, whereas nonlinear effects are treated in the time domain. A refined numerical split-step algorithm with automatic step size control ensures a high accuracy combined with high speed.
Arbitrary chromatic dispersion can be included by specifying a wavelength-dependent function for the GVD (group velocity dispersion) or the effective refractive index. If the integrated mode solver is used, this can also calculate the mode dispersion.
Not only a simple Kerr nonlinearity can be modeled, but also a nonlinearity with a delayed nonlinear response, leading to stimulated Raman scattering. The user can specify an arbitrary response function, not only parameters of a given functional form, so that Raman gain with arbitrary gain spectrum can be treated.
Over 80 functions of the script language are available for defining model inputs and retrieving all sorts of pulse properties. The results can also be inspected using the convenient interactive pulse display window.
Although various competing products can also simulate ultrashort pulse propagation, RP Fiber Power is probably unique in terms of its high flexibility, as results from the full control with a powerful script language. Therefore, you may use it for a wide range of simulations – even going far beyond what the developer has anticipated.
Speed and Reliability
The numerical solution of the equations for the steady state in such a model is a challenging task. It requires an iterative procedure. If this procedure is not very well worked out, its convergence properties can strongly depend on the situation. (See the Spotlight article of 2008-11-08 with a discussion of the validation of such simulations.)
RP Fiber Power has been very carefully designed, using an iterative procedure with automatically adapted parameters. As a result, the software achieves a surprisingly high speed and reliable convergence in a very wide range of cases. This includes cases with laser operation very close to threshold or far above threshold, amplification with backward pumping and strong ASE, etc., where other software may have extreme problems in finding the correct solution, or doing this within a reasonable time.
The diagram above shows with a color scale how the number of iterations required in a fiber laser model varies when the pump power and the output coupler transmission are systematically varied, and the configuration is each time calculated from scratch. There is a slight increase of computation time near threshold, but the computation for one configuration (with 3 optical channels in that case) typically requires only a tiny fraction second on an ordinary PC, despite the good spatial resolution (with 100 points in z direction).
Ultrashort pulse propagation is also simulated with a very robust and efficient algorithm, automatically taking finer numerical steps only where needed, based on a monitoring of numerical errors.
