Modeling and Simulation of Piecewise Regular Multimode Fiber Links Operating in a Few-Mode Regime
This work presents an alternative model of multimode fiber links with conventional silica weakly-guiding graded-index irregular multimode fibers under a few-mode optical signal propagation generated by laser source. The proposed model is based on the piecewise regular representation. It takes into account launch conditions, differential mode delay, both lower- and higher-order mode chromatic dispersion, differential mode attenuation, and mode mixing and power diffusion occurring due to real fiber irregularity and micro- and macrobends. We present some results of introduced model approbation with following pulse propagation simulations. A close matching with measured pulse responses at the output of test fibers is noticed.
Silica multimode graded-index fibers are used in a wide variety of applications ranging from on-board to in-premises networks links with length not more than 1–2 km. Since IEEE 802.3z standard was ratified and commercial available SFP transceivers with Vertical Cavity Surface Emitting Lasers (VCSELs) appeared on telecommunication market, multimode fibers became very popular both for in-building structural cabling systems (SCS) and typical distributed networks with backbone/vertical fiber cabling systems, private networks in premises, and campus environments. Nowadays the most top applications of multimode fibers are associated with data center SCS, high bit-rate storage area networks, and radio-over-fiber (RoF) techniques over already installed multimode fiber infrastructure inside building [1–3].
Modern commercial multi-Gigabit transmission transceivers are realized on Vertical Cavity Surface Emitting Lasers (VCSELs) or single-mode Fabry-Perot laser diodes (LD) [1–4]. Because emission from the conventional VCSEL usually consists of about 5 or 6 transversal modes with maximal azimuthal order not more than 3, and single-mode LD injects just fundamental and lower-order modes and , only several guided modes are excited in multimode fiber link [5–7]. Therefore optical signal propagates over multimode fiber link in a so-called few-mode regime, and passage to the simulation of a few-mode pulse transmission over multimode fiber requires taking into account both “individual” dispersion parameters of mode component with particular order (amplitude, attenuation, delay, chromatic dispersion, etc.) and mode coupling.
A variety of methods have been developed for modeling and simulation of laser-based multi-Gigabit data transmission over multimode fibers. Monograph  can be considered as fundamental complete basis work. The following groups of publications—[8–15]—should also be noted: their authors are affiliated to IEEE 802.3z, 802.3ae, 802.3aq, and 802.3ba Task Forces. Other works are also concerned with the design and simulation of multi-Gigabit multimode optical communication systems. Thus model [16–18] was realized in commercial RSoft tool ModeSYS, [19, 20] are based on vector methods, and  applies ray-tracing analysis. All mentioned models [4–20] are primarily assigned for the estimation of differential mode delay (DMD), which is the main issue of pulse distortion for laser-based multi-Gigabit data transmission over multimode fiber links and is concerned with theoretical research of DMD dependence on source parameters, launch conditions, and real fiber refractive index profile defects and deviations from the optimum graded-index profile form. They also provide accurate simulation technique of signal generation, conversion and processing in transceivers of multi-Gigabit communication systems.
A generalized chromatic dispersion parameter is applied in some models mentioned above. It is based on well-known empiric formula, which is cited in the most of fiber specifications and data sheets . However this approach does not differ of lower- and higher-order modes and does not solve the problem of “individual” chromatic dispersion estimation for guided mode with particular order. Also mode coupling is neglected in most of the listed above. It is justified by the short length of multimode fiber link. In some works correction factors or modal noise penalty are applied, but only mode power redistribution at fiber interconnections in fiber optic distribution systems, consolidation points, patch-panels, and so forth, is recalculated without taking into account both real fiber irregularities and fiber bends.
Here we introduce an alternative model of multimode fiber links with conventional silica weakly-guiding graded-index irregular multimode fibers under a few-mode optical signal propagation generated by laser source. Proposed model is based on the piecewise regular representation. It takes into account launch conditions, differential mode delay, both lower- and higher-order mode chromatic dispersion, differential mode attenuation, and mode mixing and power diffusion occurring due to real fiber irregularity and micro and macro-bends. We also present some results of introduced model approbation with following pulse propagation simulations and matching with measured pulse responses at the output of test fibers.
2. Brief Overview of Methods for Simulation of Pulse Propagation over Irregular Optical Fibers
A several methods have been developed for computing electromagnetic wave fields during propagation over irregular waveguides. Monograph  is one of the fundamental works of this direction. Effects of wave interaction in irregular waveguides are exact described by the coupled mode theory. Complete overview, description, classification, and list of applications of coupled mode theory methods can be found in monograph .
Generally coupled mode theory is based on representation of mode fields along longitudinal cross-section of irregular optical waveguide in the kind of mode field superposition of complete system of guided and leaky modes in regular optical waveguide. This expansion varies along the fiber longitudinal axis and it is described by system of coupled mode equations. The system of equations is formulated by substitution of mentioned irregular optical waveguide mode field representation to Maxwell equations  or Helmholtz equation [26, 27]. Result of substitution is infinite system of differential equations of the first or second order (it depends on initial equations—Maxwell or Helmholtz) for each mode amplitude depending on longitudinal coordinate .
The system of equations exactly describes mode fields in longitudinal cross-section of irregular optical waveguide, but it has not general analytical solution. Problem can be much simplified for the limited number of modes—just one or two [28–30]. For example coupled mode theory is effectively applied for analysis of the fundamental mode propagation over irregular optical fiber under taking into account the birefringence effect [25–27, 31]. However few-mode regime of signal transmission over multimode fibers requires increase of mode components that much complicates solutions of the system of equations. This fact limits coupled mode theory application in the case of random fluctuations of fiber core geometry parameters along the longitudinal axis .
In [32, 33], mode components of few-mode signal are combined in mode groups with roughly equal propagation constants. In works [34–36], a passage to infinite number of mode components (mode continuum) is proposed followed by replacing the sum of discrete mode component fields by integration. By neglecting core geometry variation, random fiber macrobends is considered as the main factor of mode coupling in [37–39]. Also in [40–44] a special curvilinear coordinates are applied to describe mode coupling via fiber bends. In [45–47], mode mixing and power diffusion is simulated by combination of coupled mode theory and ray-tracing analysis.
The system of incoherently coupled nonlinear Schrodinger equations is applied for mathematical description of pulse propagation over multimode fiber in [48–54]. Although the generalized system of equations was written for infinite number of modes by taking into account fiber nonlinearity, the main application of this method is used just for simulation of coupling between two orthogonally polarized modes in single mode fibers occurring due to birefringence.
Piecewise regular representation [25, 26, 55] of irregular optical waveguide is an alternative method to coupled mode theory. It is basic for direct methods like finite-difference time-domain (FDTD) method—one of the most accurate, universal, and “popular” numerical methods, which is widely used for analysis of electromagnetic wave propagation over waveguides with arbitrarily complicated configuration (e.g., [56–62], etc.). Also well-known direct rigorous methods are method of moments [58, 62, 63], finite-difference frequency-domain (FDFD) method [55, 58, 62], finite element method [55, 58, 62], wavelets [64–66] and their modifications, and combinations with elements of coupled mode theory [41–44, 67–73].
Rigorous methods provide accurately computing of the transversal mode field distribution at the particular distance of emission source and time interval under the lowest error. The main issue of their application is complicated formulation of boundary and initial conditions according to individual problem and waveguide properties and configuration. They are very effective for the modeling of waveguide structures with quite short dimensions—not more than —like integral optics components. However simulation of pulse propagation over conventional fiber link (with length of hundreds or even tens meters) by the mentioned direct rigorous methods requires extra computing resources and is almost impossible.
In terms of piecewise regular representation of both LAN, MAN (metropolitan area network) and long-haul fiber links, matrix methods [62, 74–80], recurrent methods [81–85] and quasi-analytic time-domain models based on quasianalytic methods [86–89] and spatial transform methods [55, 86, 90–95] are widely used. The most of them were developed for simulation of pulse propagation over single mode fibers.
From this point of view, method described in [96–98] can be considered as an exception. It combined coupled mode theory with piecewise regular representation and was developed for modeling of signal transmission processes over long irregular multimode fiber links. Here a random fiber bends are preset and supposed as the main reason of mode coupling occurring. Also it takes into account birefringence. However the method neglects core geometry variation along longitudinal fiber axis and considers only multimode fibers with ideal infinite -profiles.
Therefore the most considered methods above for modeling and analysis of pulse propagation over optical fibers should be adapted for irregular multimode fiber link operating in a few-mode regime. Model has to take into account the whole total selected excited guided mode composition of transferred few-mode optical signal and their “individual” parameters during propagation over multimode fibers with real graded-index profiles differing form ideal -profiles by local defects and refractive index fluctuation. Moreover, in this case division of mode composition on lower- and higher-order mode groups with averaged propagation constant is unacceptable. Also model has to take into account launch conditions and mode coupling due to both variation of core geometry parameter along fiber length and micro- and macrobends of fiber occurring in real cable links.
3. Model of Piecewise Regular Multimode Fiber Link
In this work we also apply piecewise regular representation combined with general approach of split-step method  to simulate a processes of mode mixing and power diffusion due to mode coupling (Figure 1). Here single silica weakly guiding circular multimode fibers with arbitrary axial symmetric refractive index profile with single continue outer claddings are considered. According to piecewise regular representation the fiber is divided into regular spans with length . Inside the span fiber geometry parameters are considered as constant, and modes propagate independently without interaction and mixing. It is supposed that each excited guided mode with propagation constant, varying from one regular span to another span, satisfies to cut-off condition for whole regular spans composing fiber. Also it is assumed that at the link transmitter end each excited guided mode starts transferring single optical pulse in particular form identical to input signal (e.g., Gaussian). During pulse propagation over regular span its amplitude decreases due to mode attenuation. The signal is mainly distorted due to difference between group velocities and amplitudes of modes, that is, DMD effect. Also transferred by each mode pulse, spreading due to chromatic dispersion for particular mode is taken into account.
Boundaries of regular spans can be represented generally as ideal axially alignment splice of two almost similar optical fibers with mismatching parameters. However it is correct only for “straight” fibers. Therefore it is propose to simulate fiber bends by representation of boundaries in the form of splice of two mismatched fibers with some low angular misalignment. Mode power redistribution between the amplitudes of signal components as a result of mode interaction is estimated by mode coupling coefficient matrix computing at the joints of regular spans. Here only guided modes are considered as the main issue under pulse dynamics research during propagation over multimode fiber link in a few-mode regime. However, in fact, power loss due lo component transformation from guided to leaky mode, and reflections are also indirectly taken into account.
At the receiver end the resulting pulse envelope is considered as superposition of all mode components of signal. Here it is propose to apply a well-known expression for the frequency response of the signal transferred by mode components over regular multimode fiber with length [4, 8–15, 18, 99–101]: where is direct Fourier transform; is initial pulse at the transmitter end; and are starting amplitude and mode attenuation of th guided mode ; and are first- and second-order dispersion parameters. These dispersion parameters are elements of well-known Taylor series expansion approximation of propagation constant frequency dependence [2, 55, 99]: where Here is mode delay and is group velocity dispersion (GVD) associated with chromatic dispersion of th guided mode .
Therefore, according to introduced piecewise regular representation of irregular multimode fiber link, the frequency response of a few-mode optical signal transferred by guided modes over irregular multimode fiber with length under given particular length of regular span by taking into account expression (1) can be written in the following form: where ; is the integer part of real number .
The resulting pulse response at the receiver end of irregular multimode link is computed by the following simple expression: where is inverse Fourier transform; is the complex conjugate of .
It is obvious that the proposed model requires fast and simple method for the analysis of the optical fibers with complicated refractive index profile form close to profiles of real fibers, which provides evaluation of parameters not only for the fundamental but also higher-order guided modes. For this purpose we propose to apply an extension of modified Gaussian approximation (EMGA) generalized for estimation of any order guided mode parameters, propagating over silica weakly guiding optical fibers with an arbitrary axial-symmetric refractive index profile . EMGA is based on combination of Gaussian approximation and stratification method. It permits a much reduce computing time, especially for the calculations of higher-order mode parameters under the low error  by making a passage to analytical expressions both for the variational expression and characteristic equation. Below a detailed description of EMGA and its application for the proposed model is presented.
4. Extension of Modified Gaussian Approximation
EMGA is based on classical Gaussian approximation  of radial mode field distribution in the weakly guiding optical waveguide with an arbitrary refractive index profile by the well-known Laguerre-Gauss function expression [4, 25, 100], describing a mode field distribution in weakly guiding optical waveguide with ideal infinite parabolic index profile: where is Laguerre polynomial.
EMGA leads to equivalent normalized mode field radius estimation by solving a characteristic equation, which is derived from propagation constant variational expression under following passage to square core mode parameter variational expression, written for analyzed weakly guiding optical waveguide with given refractive index profile. Parameter is basic for this method and completely defines mode transmission parameters.
Unlike known methods and their modifications [25, 104–110], based on conventional Gaussian approximation, in EMGA optical fiber with an arbitrary graded axial-symmetric index profile is considered as multicladding optical fiber. Therefore, refractive index profile inside core region can be represented in the form of the set of layers in which the refractive stays a constant [102, 103]: and any profile function related to refractive index profile as where is profile height parameter, can be written in terms of profile parameter : where where is refractive index of layer (); is the maximal core refractive index; is cladding refractive index.
This refractive index profile representation permits writing the variational expression for core mode parameter and characteristic equation for normalized equivalent mode field radius in the form of finite nested sums as follows [15, 16]: whereand is coefficient of polynomial representation in the form of power series [111, 112]: is normalized frequency; is wavenumber; is wavelength.
Therefore analysis of weakly guiding single-cladding optical fiber with an arbitrary profile leads to the following. Refractive index profile is represented by the profile function (9a) and (9b) in the form of layers. Fiber parameters and mode orders and are substituted to the characteristic equation (11). By means of numerical solution (11), the normalized equivalent mode field radius will be obtained. Then is substituted to the expression (10), and mode core parameter is estimated; that permits evaluating the propagation constant for guided mode by the well-known expression [4, 25–27, 55, 100]:
Optical confinement factor can be considered as the second criterion for the identification of the ghost solutions: If we take into account Gaussian approximation, parameter is defined by analytical expression derived in  from the generalized integral form for weakly guiding optical fibers presented in : where
5. Mode Group Velocity
According to what is mentioned above, realization of introduced model of piecewise regular multimode fiber link requires a passage to analytical expressions for mode group velocity. Guided mode group velocity is related to propagation constant by the well-known formula [4, 25–27, 55, 100]:
By squaring and following differentiation of the expression (14) with respect to wavelength, the first partial derivative of the square of propagation constant will be obtained:
6. Mode Chromatic Dispersion Parameter
The same technique is applied to derive the analytical expressions for mode chromatic dispersion parameter associated with GVD. Chromatic dispersion parameter is related to propagation constant by the well-known formula [4, 25–27, 55, 100]:
By differentiating (20), a second partial derivative of the square of propagation constant will be obtained:
Therefore expression for leads to following form:
By differentiating expressions (22) with respect to wavelength, required second derivatives of square core mode parameter and normalized equivalent mode field radius will be obtained:
7. Material Dispersion and Profile Parameters
We shall apply a well-known Sellmeier equation to take into account the refractive index dependence on wavelength [4, 25–27, 55, 100]: where and are Sellmeier’s coefficients ( is also denoted as the resonance wavelength) which has been empirically measured for GeO2-SiO2 glasses [113, 114] under the several particular dopant concentration. Here we shall apply the method described in  to estimate Sellmeier coefficients at the graded-index profile points.
A passage from the differentiation operator to its square is used by analogously to propagation constant derivatives. This passage will much simplify expressions for the first and second derivatives of the refractive index:
Therefore the first and second derivatives of the profile height parametercan be expressed as follows:
Then derivatives of profile parameter defined by formula (9b) are determined by the following expressions:
8. Differential Mode Attenuation
We shall apply a simple empirical relation proposed in work  using experimental data from  to estimate differential mode attenuation depending on mode order: where is principal mode number; is the attenuation of lower-order modes (it is supposed equal to the attenuation at the correspondence wavelength mentioned in fiber specification); is total number of modes satisfying the cutoff condition for analyzed fiber: where is gradient factor of smoothed -profile.
9. Mode Coupling
Nowadays calculation of mode coupling coefficients is still one of the most well-known and simple methods for optical waveguide junction analysis. It is widely applied for estimation of insertion loss and reflection at the splices of the same type optical fibers with nonidentical technologic parameters [117–121], modeling, and research of the influence of launch conditions on optical waveguide mode composition excitation [25, 122–124] including the modeling of a few-mode optical signal propagation over multimode fibers [4, 8, 10].
Therefore substitution of radial mode field distribution expressions to integrals (34) is required for estimation of mode coupling coefficients. The simplest mode field structure corresponds to the fundamental mode . As a result, it permits obtaining analytical formulas for the overlap integral (34) by applying simple approximation expressions or by a particular function series expansion for the fundamental mode field representation. These algorithms are widely used for analysis of the single mode fiber splices [25, 117–123].
Contrariwise, higher-order mode field structure is more complicated. For this reason in certain papers concerning multimode fiber analysis [4, 8, 10] it is propose to integrate expression (34) numerically or to substitute simple approximation expressions, which are solutions of wave equation for the fibers with ideal step or infinite parabolic refractive index profile. However this substitution is not correct entirely: as it was mentioned above, real refractive index profiles differ much from the model ideal smoothed -profiles  even for the last generation of commercial silica optical fibers both single mode and multimode.
Here we propose to apply the overlap integral method combined with introduced EMGA. Under the EMPG application, the use of the algorithm mentioned above for the following mode coupling estimation becomes acceptable: here local features of real silica optical fiber refractive index profile are taken into account that provides decreasing the error of calculations. Thereby the fiber splice analysis will led to substitution of simple approximation expression (6) to the overlap integral (34) under preliminary estimation of equivalent mode field radiuses of injected and excited modes with given azimuthal and radial orders by EMPG. In terms of introduced model of piecewise regular multimode fiber link, this approach provides both simulation of mode mixing and power diffusion by recalculation of the matrix of mode coupling coefficients at the junctions of regular span boundaries and modeling of initial mode composition excitation in multimode fiber by laser source at the transmitter end under taking into account launch condition (single mode pigtail, central launching, offset, air gap, special matching fiber, etc.).
Another advantage of proposed method combining the overlap integral method and EMGA is the ability of passage to analytical expressions both for mode coupling coefficients at the ideal centralized splice and splice with offset or tilt for modes with arbitrary order. In recent work  the derivation of them is described in detail. Therefore here only final heavy analytical expressions are presented.
First of all let us consider the idealized model of diverse optical fiber splice: jointed mismatched type fibers have ideal smooth ends cleaved under 90° to the core axes, and there are no any longitudinal, axial, or angular misalignments between core centers. On the one hand it is very simplified fiber splice description; on the other hand this ideal splice representation corresponds to central launching during the modeling of fiber excitation by laser source.
A passage from the generalized form of the overlap integral (1) to analytical expression for the arbitrary order mode coupling coefficient estimation at the central splice of the optical fibers with mismatched parameters and without any misalignments was proposed in work . Analytical expression deriving is described in detail in . Finally the formula for arbitrary order mode coupling coefficient at the central splice is written as follows: Here mode coupling occurs only for modes with the same azimuthal order ; Γ is Gamma function; and are injected and excited mode field radiuses.
The analytical expression for arbitrary order mode coupling coefficient at the optical fiber splice under low angle misalignment has the following form : where is confluent hypergeometric function of the 1st kind [111, 112]: whereand are Laguerre polynomial expansion factors of (13a), (13b); is refractive index of launching medium (air gap, core of adjusting/exiting fiber, etc.).
The final expression for mode coupling coefficients for 0th azimuth order modes and at the splice with transverse offset has the form : where the lower and upper limits of second and forth sums are and and the “odd” factorial is :
The analytic formula mode coupling coefficients for the 1st azimuth order modes and at the splice with transverse offset is written in the following form :
Finally the expression for mode coupling coefficients for the arbitrary higher-order () modes and at the splice with transverse offset leads to the following form: where
10. Results and Discussion
Commercial GeO2-doped graded-index multimode fiber samples of different generations of categories OM2 and OM2+/OM3 with length about 300–500 m packed in coils and conventional spools (Figure 2) were selected for experimental approbations of introduced model.
The refractive index profiles of selected samples were preliminary measured  by optical fiber analyzer EXFO NR-9200 via refracted near field scan technique . Profiles reconstructed from the protocol data are presented in Figure 3. Also each fiber sample passed an OTDR test with following fiber length measurement by backscattering method. The OTDR Hewlett Packard E6000A with single mode plug-in unit 1310/1550 was used, and test was produced at the wavelength nm [128, 129].
We applied a typical optical pulse response workstation R2D2 for DMD diagnostics [128, 129]. Simplified block scheme of R2D2 is presented in Figure 4. The experimental setup and software R2D2 screenshot are presented in Figures 5 and 6.
The FP LD emits light with a wavelength nm. The input pulse is shown in Figure 7. It has the full width at a half maximum (FWHM) of about 340 ps. The first series of measurements did not apply special launching conditions. Here only conventional FC/PC adaptors were applied. Therefore an angular misalignment was taking into account according to  during the modeling of transmitter end of link.
We simulated fiber irregularity by core diameter variation. It was set from the data of fiber diameter measurements been produced during fiber drawing  (Figure 8). Micro- and macrobends were modeled by equivalent low angular misalignment at the boundaries of regular spans with length m . Measured and simulated pulse responses of fibers OM2+/OM3 and OM2 with refractive index profiles shown in Figure 3 are presented in Figures 9(a) and 9(b). Here a close matching can be noticed. Thus for OM2+/OM3 fiber the error of the FWHM is 1.16%, and output pulse dispersion error is 7.42%. For OM2 fiber the errors of second peak amplitude and time scale position is 1.63% and 0.43%, and DMD error is 1.39%.
Another series of tests were produced under the control of launching conditions. Simplified block-scheme of the experimental setup is shown in Figure 10. Here emission from the LD launched to tested multimode fiber via single mode pigtail (ITU-T Rec. G.652) which was spliced with multimode fiber with precision offset. Refractive index profile of the single mode fiber is presented in Figure 11. The offset was realized by Ericsson FSU-975 program number 8 “Attenuator” . Here parameter “attenuation” was set to 0, “offset”—to particular value, and “ECF factor” (surface tension compensation factor) also was set to 0 to avoid deformations in the heat zone of the offset fusion splice (Figure 12). Some measured and simulated pulse responses of fibers OM2+/OM3 and OM2 under particular offset launching conditions are presented in Figure 13. A close matching also can be noticed.
(a) Central launching
|(b) Offset m|
|(c) Offset m|
(d) Central launching
|(e) Offset m|
|(f) Offset m|
An alternative model of multimode fiber links with conventional silica weakly-guiding graded-index irregular multimode fibers under a few-mode optical signal propagation generated by laser source is introduced. The model is based on the piecewise regular representation. It takes into account launch conditions, differential mode delay, both lower- and higher-order mode chromatic dispersion, differential mode attenuation, and mode mixing and power diffusion occurring due to real fiber irregularity and micro- and macrobends. Results of theoretical approbation are presented. A good agreement between developed model and measurements of pulse response were obtained.
“BSRIA update: structured cabling market, data centers, IIMS, and potential new studies,” BSRIA Presentation Materials, January 2009.View at: Google Scholar
Sh. Bois, “Next generation fibers and standards,” FOLS Presentation Materials, October 2009.View at: Google Scholar
S. Bottacchi, Multi-Gigabit Transmission over Multimode Optical Fibre. Theory and Design Methods for 10GbE Systems, John Wiley & Sons, 2006.
C. Degen, I. Fischer, and W. Elsäßer, “Transverse modes in oxide confined VCSELs: influence of pump profile, spatial hole burning, and thermal effects,” Optics Express, vol. 5, no. 3, pp. 38–47, 1999.View at: Google Scholar
L. Raddatz, I. H. White, D. G. Cunningham, and M. C. Nowell, “An experimental and theoretical study of the offset launch technique for the enhancement of the bandwidth of multimode fiber links,” IEEE Journal of Lightwave Technology, vol. 16, no. 3, pp. 324–331, 1998.View at: Publisher Site | Google Scholar
J. S. Abbot, M. J. Hackert, D. E. Harshbarger, D. G. Cunningham, T. Ch. di Minico, and I. H. White, “Analysis of multimode fiber behavior with laser sources in the development of the Gigabit Ethernet fiber optic specifications,” in Proceedings of the 47th International Wire & Cable Symposium (IWCS '98), pp. 897–907, Philadelphia, Pa, USA, 1998.View at: Google Scholar
P. K. Pepeljugoski and D. M. Kuchta, “Design of optical communications data links,” IBM Journal of Research and Development, vol. 47, no. 2-3, pp. 223–237, 2003.View at: Google Scholar
A. Gholami, D. Molin, P. Matthijsse, G. Kuyt, and P. A. Sillard, “A complete physical model for Gigabit Ethernet optical communication systems,” in Proceedings of the 57th International Wire & Cable Symposium (IWCS '08), pp. 289–294, Providence, RI, USA, 2008.View at: Google Scholar
J. Morikuni, P. Mena, B. K. Whitlock, and R. Scarmozzino, “Simulation-based prediction of multimode fiber bandwidth for 10 Gb/s systems,” in Proceedings of the 15th Annual Meeting of the IEEE Laser and Electro-Optics Society (LEOS '02), vol. 2, pp. 604–605, Glasgow, Scotland, November 2002.View at: Publisher Site | Google Scholar
J. Morikuni, P. Mena, B. K. Whitlock, and R. Scarmozzino, “Link-level design, analysis, and simulation of multimode data communication systems,” in Proceedings of the 19th Annual National Fiber Optic Engineers Conference (NFOEC '03), pp. 858–867, Orlando, Fla, USA, 2003.View at: Google Scholar
G. Shaulov and B. Whitlock, “Multimode fiber communication system simulation,” IEEE 802.3aq 10GBase-LRM Task Force. Presentation Materials, July 2004.View at: Google Scholar
F. J. Achten, M. Bingle, and B. P. de Hon, “DMD measurement and simulation on short length, 850-nm laser-optimized, 50 mm core diameter, GIMM optical fibres,” in Proceedings of the 12th Symposium on Optical Fiber Measurements, pp. 101–104, Boulder, Colo, USA, 2002.View at: Google Scholar
A. Sengupta, “Simulation of 10GbE multimode optical communication systems,” in Proceedings of the 57th International Wire & Cable Symposium (IWCS '08), pp. 320–326, Providence, RI, USA, 2008.View at: Google Scholar
N. Sharma, A. S. Chitambar, H. Ch. Bianchi, and K. Sivaprasad, “Model for differential mode delay distribution in multimode fibers,” Progress in Electromagnetic Research, vol. 1, no. 4, pp. 469–472, 2005.View at: Google Scholar
A. Listvin, V. Listvin, and D. Shvyrkov, Fibers for Telecommunication Lines, LESARart, Moscow, Russia, 2003.
B. Katsenelenbaum, Theory of Irregular Waveguides with Slow Varying Parameters, Mir, Moscow, Russia, 1961.
A. Baribin, Electrodynamics of Waveguide Structures. Theory of Perturbation and Coupled Mode Theory, Fizmalit, Moscow, Russia, 2007.
A. Snyder and J. Love, Optical Waveguide Theory, Chapman & Hall, London, UK, 1983.
H. G. Unger, Planar Optical Waveguides and Fibres, Clarendon Press, Oxford, UK, 1977.
D. Marcuse, Light Transmission Optics, Van Nostrand Reinhold Company, London, UK, 1972.
S. Geckeler, “Pulse broadening in optical fibers with mode mixing,” Applied Optics, vol. 18, no. 13, pp. 2192–2198, 1979.View at: Google Scholar
I. O. Zolotovski and D. I. Sementsev, “Dispersion influence to Gaussian pulse transformation in periodically irregular fiber,” Optics and Spectroscopy, vol. 85, no. 2, pp. 304–308, 1998.View at: Google Scholar
F. Tian, Y.-Z. Wu, and P.-D. Ye, “Analysis of polarization fluctuation in single-mode optical fibers with continuous random coupling,” IEEE Journal of Lightwave Technology, vol. 5, no. 9, pp. 1165–1168, 1987.View at: Google Scholar
L. Jeunhomme and J. P. Pocholle, “Mode coupling in a multimode optical fiber with microbends,” Applied Optics, vol. 14, no. 10, pp. 2400–2405, 1975.View at: Google Scholar
M. Miyagi, S. Kawakami, M. Ohashi, and S. Nishida, “Measurement of mode conversion coefficients an mode dependent losses in a multimode fibers,” Applied Optics, vol. 17, no. 20, pp. 3238–3244, 1978.View at: Google Scholar
D. Świerk and A. Heimrath, “Coupled power equations in multimode waveguides: several, random, correlated imperfections,” Journal of Modern Optics, vol. 43, no. 11, pp. 2301–2310, 1996.View at: Google Scholar
R. Vacek, “Electromagnetic wave propagation in general multimode waveguide structures exhibiting simultaneously attenuation, dispersion and coupling phenomena,” Journal of Microwaves and Optoelectronics, vol. 1, no. 3, pp. 42–52, 1998.View at: Google Scholar
R. Olshansky, “Mode coupling effects in graded-index optical fibers,” Applied Optics, vol. 14, no. 4, pp. 935–945, 1975.View at: Google Scholar
M. Skorobogatiy, M. Ibanescu, S. G. Johnson et al., “Analysis of general geometric scaling perturbations in a transmitting waveguide: fundamental connection between polarization-mode dispersion and group-velocity dispersion,” Journal of the Optical Society of America B, vol. 19, no. 12, pp. 2867–2875, 2002.View at: Google Scholar
M. Skorobogatiy, S. A. Jacobs, S. G. Johnson, and Y. Fink, “Geometric variations in high index-contrast waveguides, coupled mode theory in curvilinear coordinates,” Optics Express, vol. 10, no. 21, pp. 1227–1243, 2002.View at: Google Scholar
D. Gloge, “Bending loss in multimode fibers with graded and ungraded core index,” Applied Optics, vol. 11, no. 11, pp. 2506–2513, 1972.View at: Google Scholar
B. Crosignani, A. Cutolo, and P. di Porto, “Coupled-mode theory of nonlinear propagation in multimode and single-mode fibers: envelope solitons and self-confinement,” Journal of Optical Society of America, vol. 72, no. 9, pp. 1136–1141, 1982.View at: Google Scholar
S. Solimeno, B. Crosignani, and P. di Porto, Guiding, Diffraction and Confinement of Optical Radiation, Academic Press, San Diego, Calif, USA, 1986.
Y. S. Kivshar and G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals, Academic Press, San Diego, Calif, USA, 2003.
K. Okamoto, Fundamentals of Optical Waveguides, Academic Press, San Diego, Calif, USA, 2000.
S. T. Chu and S. K. Chaudhuri, “Finite-difference time-domain method for optical waveguide analysis,” Progress in Electromagnetic Research, vol. 11, pp. 255–300, 1995.View at: Google Scholar
K. Sawaya, “Numerical techniques for analysis of electromagnetic problems,” IEICE Transactions on Communications, vol. 83, no. 3, pp. 444–452, 2000.View at: Google Scholar
G. Bao, L. Cowsar, and W. Masters, Mathematical Modeling in Optical Science, SIAM, Philadelphia, Pa, USA, 2001.
A. B. Shvartsburg, Time Domain Optics of Ultrashort Waveforms, Clarendon Press, Oxford, UK, 1996.
A. B. Shvartsburg, Impulse Time-Domain Electromagnetics of Continues Media, Birkhäauser, Boston, Mass, USA, 1999.
D. M. Sullivan, Electromagnetic Simulation Using the FDTD Method, IEEE Press, New York, NY, USA, 2000.
V. A. Neganov, O. V. Osipov, S. B. Rayevskiy, and G. P. Yarovoy, Electrodynamics and Wave Propagation, Radiotechnika, Moscow, Russia, 2007.
R. Mitra, Computational Methods for Electrodynamics, Mir, Moscow, Russia, 1977.
G. W. Wornell, Signal Processing With Fractals: A Wavelet-Based Approach, Prentice-Hall, Upper Saddle River, NJ, USA, 1996.
T. K. Sarkar, C. Su, R. Adve, M. Salazar-Palma, L. Garcia-Castillo, and R. R. Boix, “A tutorial on wavelets from an electrical engineering perspective, part 1: discrete wavelet techniques,” IEEE Antennas and Propagation Magazine, vol. 40, no. 5, pp. 49–68, 1998.View at: Google Scholar
T. K. Sarkar and C. Su, “A tutorial on wavelets from an electrical engineering perspective, part 2: the continuous case,” IEEE Antennas and Propagation Magazine, vol. 40, no. 6, pp. 36–48, 1998.View at: Google Scholar
A. Vukovic, P. Sewell, T. M. Benson, and P. C. Kendall, “Spectral method applied to design of spotsize converters,” Electronics Letters, vol. 33, no. 25, pp. 2121–2123, 1997.View at: Google Scholar
M. Reed, T. M. Benson, P. Sewell, P. C. Kendall, G. M. Berry, and S. V. Dewar, “Free space radiation mode analysis of rectangular dielectric waveguides,” Optical and Quantum Electronics, vol. 28, no. 9, pp. 1175–1179, 1996.View at: Google Scholar
A. Niiyama and M. Koshiba, “3-dimensional eam propagation analysis of nonlinear optical fibers,” IEICE Transactions on Communications, vol. E80-B, no. 4, pp. 522–527, 1997.View at: Google Scholar
R. Magnanini and F. Santosa, “Wave propagation in a 2-D optical waveguide,” SIAM Journal on Applied Mathematics, vol. 61, no. 4, pp. 1237–1252, 2000.View at: Google Scholar
G.-W. Chern, J.-F. Chang, and L. A. Wang, “Modeling of nonlinear pulse propagation in periodic and quasi-periodic binary long-period fiber gratings,” Journal of the Optical Society of America B, vol. 19, no. 7, pp. 1497–1508, 2002.View at: Google Scholar
G. B. Malykin, V. I. Pozdnyakova, and I. A. Shereshevskii, “Random groups in the optical waveguides theory,” Journal of Nonlinear Mathematical Physics, vol. 8, no. 4, pp. 491–517, 2001.View at: Google Scholar
N. A. Galchenko, “Matrix theory of excitation of electromagnetic waves in irregular waveguides,” Radiophysics and Quantum Electronics, vol. 40, no. 6, pp. 744–751, 1997.View at: Google Scholar
O. O. Silitchev, “Matrix method for computing of coherent laser pulse propagation,” Quantum Electronics, vol. 20, no. 10, pp. 983–990, 1993.View at: Google Scholar
T. A. Martynova, “On the computing of irregular waveguides with local layered inhomogeneity,” Radiophysics and Quantum Electronics, vol. 18, no. 5, pp. 1178–1188, 1975.View at: Google Scholar
T. A. Ramadan, “A recurrence technique for computing the effective indexes of the guided modes of coupled single-mode waveguides,” Progress in Electromagnetic Research, vol. 4, pp. 33–46, 2008.View at: Google Scholar
V. V. Merkoulov and I. S. Sineva, “On the repeated reflections in the non-homogeneous transmission line,” Journal of Radio Electronics, no. 5, paper 8, 2000.View at: Google Scholar
S. G. Akopov, V. N. Korshunov, B. S. Solov'ev, and B. N. Fomichev, “Energetic characteristics of distribution functions of irregularities of geometric structure of optic fiber,” Elektrosvyaz, no. 9, pp. 9–11, 1992.View at: Google Scholar
S. G. Akopov, V. N. Korshunov, B. S. Solovyov, and V. N. Fomichev, “Estimation of characteristics of random spatial process of irregularity geometry structure distribution in optical fiber,” Electrosvyaz, no. 10, pp. 20–22, 1992.View at: Google Scholar
Y. Y. Lu, “Some techniques for computing wave propagation in optical waveguides,” Communications in Computational Physics, vol. 1, no. 6, pp. 1056–1075, 2006.View at: Google Scholar
U. Osterberg, “Signal processing in optical fibers,” Modern Signal Processing, vol. 46, pp. 301–316, 2003.View at: Google Scholar
D. Arakaki, W. Yu, and R. Mittra, “On the solution of a class of large body problems with partial circular symmetry (multiple asymmetries) by using a hybrid-dimensional Finite-Difference Time-Domain (FDTD) method,” IEEE Transactions on Antennas and Propagation, vol. 49, no. 3, pp. 354–360, 2001.View at: Publisher Site | Google Scholar
K. A. Amaratunga, “A wavelet-based approach for compressing kernel data in large-scale simulations of 3D integral problems,” Computing in Science and Engineering, no. 7-8, pp. 34–35, 2000.View at: Google Scholar
G. P. Agrawal, Nonlinear Fiber Optics, Academic Press, San Diego, Calif, USA, 1989.
M. J. Adams, An Introduction to Optical Waveguides, John Wiley & Sons, New York, NY, USA, 1981.
J. J. R. Hansen and E. Nicolaisen, “Propagation in graded-index fibers: comparison between experiment and three theories,” Applied Optics, vol. 17, no. 17, pp. 2831–2835, 1978.View at: Google Scholar
A. V. Bourdine and O. R. Delmukhametov, “Calculation of transmission parameters of the launched higher-order modes based on the combination of a modified Gaussian approximation and a finite element method,” Telecommunications and Radio Engineering, vol. 72, no. 2, pp. 111–123, 2013.View at: Publisher Site | Google Scholar
H. Meher and S. I. Hosain, “Variational approximations for single-mode graded-index fibers: some interesting applications,” Journal of Optical Communications, vol. 24, no. 1, pp. 25–30, 2003.View at: Google Scholar
R. Tewari, S. I. Hosain, and K. Thyagarajan, “Scalar variational analysis of single mode fibers with Gaussian and smoothed-out profiles,” Optics Communications, vol. 48, no. 3, pp. 176–180, 1983.View at: Google Scholar
I. Gradstein and I. Ryjik, Tables of Integrals, GIFML, Moscow, Russia, 1963.
M. Abramovitz and I. Stegan, Handbook of Mathematical Functions with Formulas, Graphs and Mathematical Tables, Nauka, Moscow, Russia, 1979.
L. N. Binh, “Design guidelines for ultra-broadmand dispersion-flattened optical fibers with segmented-core index profile,” Tech. Rep. MECSE-14-2003, Monash University, Clayton, Australia, 2003.View at: Google Scholar
J. W. Fleming, “Dispersion in GeO2-SiO2 glasses,” Applied Optics, vol. 23, no. 24, pp. 4886–4493, 1984.View at: Google Scholar
V. A. Burdin, “Methods for computation of Sellmeier coefficients for dispersion analysis of silica optical fibers,” Infocommunication Technologies, vol. 4, no. 2, pp. 30–34, 2006.View at: Google Scholar
M. Rousseau and L. Jeunhomme, “Optimum index profile in multimode optical fiber with respect to mode coupling,” Optics Communications, vol. 23, no. 2, pp. 275–278, 1977.View at: Google Scholar
S. P. Gurji and V. B. Katok, “Loss calculation at the splices of single mode fibers with complicated refractive index profile,” Elektrosvyaz, no. 10, pp. 25–27, 1990.View at: Google Scholar
V. A. Srapionov, “Mode coupling at the splices of optical fibers with mismatched parameters,” Elektrosvyaz, no. 10, pp. 10–12, 1985.View at: Google Scholar
I. A. Avrutski, V. A. Sychugov, and A. V. Tishenko, “The study of excitation, radiation and reflection processes in corrugated waveguides,” Proceedings of IOFAN, no. 34, pp. 3–99, 1991.View at: Google Scholar
R. Chandra, K. Thyagarajan, and A. K. Ghatak, “Mode excitation by tilted and offset Gaussian beams in W-type fibers,” Applied Optics, vol. 17, no. 17, pp. 2842–2847, 1978.View at: Google Scholar
A. V. Bourdine, “Mode coupling at the splice of diverse optical fibers,” in Optical Technologies for Telecommunications, vol. 8787 of Proceedings of SPIE, pp. 878706-1–878706-12, 2012.View at: Google Scholar
A. V. Bourdine and K. A. Yablochkin, “Investigations of refractive index profile defects of silica graded-index multimode fibers of telecommunication cables,” Infocommunication Technologies, vol. 8, no. 2, pp. 22–27, 2010.View at: Google Scholar
“Optical Fiber Analyzer EXFO NR-9200/NR-9200HR,” EXFO, 2006.View at: Google Scholar
A. V. Bourdine, “Differential mode delay of different generations of silica multimode optical fibers,” Photon Express, vol. 5-6, no. 69-70, pp. 20–22, 2008.View at: Google Scholar
A. V. Bourdine, “On multimode optical fiber DMD diagnostics,” Infocommunication Technologies, vol. 6, no. 4, pp. 33–38, 2008.View at: Google Scholar
K. J. Lyytikainen, Control of complex structural geometry in optical fibre drawing [Ph.D. thesis], School of Physics and Optical Fibre Technology Centre, University of Sydney, 2004, A Thesis Submitted for the Degree of Doctor of Philosophy.
“Ericsson FSU-975. Instruction Manual,” Ericsson, 2001.View at: Google Scholar