Pulse Propagation in Optical Fibers

Size: px
Start display at page:

Download "Pulse Propagation in Optical Fibers"

Transcription

1 Pulse Propagation in Optical Fibers Arleth Manuela Gonçalves Departaento de Engenharia Electrotécnica e de Coputadores Instituto Superior Técnico Av. Rovisco Pais, 49- Lisboa, Portugal arlethgoncalves@hotail.co Abstrat This paper addresses the pulse propagation through a fiber optic syste, operating in the linear and nonlinear regies. After a brief introduction to optical fibers, we use the odal theory approach to understand the operating principle for the pulses propagating in the fiber. Then, we try to understand the effect of dispersion, as the unavoidable phenoenon to which pulses are subject, while propagating in the fiber, in the linear regie. Since dispersion strongly liits the bandwidth of the transitted signal and ay cause interference between sybols, we address several ethods for reducing the effect of dispersion, which can be very efficient: the use of dispersion copensating fibers and the propagation of solitons. In non-linear syste, the non-linear Kerr effect explains the appearance of self-phase odulation (SPM) which has a contrary action to the phenoenon of dispersion in the anoalous dispersion region, allowing the appearance of solitons, which era pulses that conserve their shape along the propagation. The propagation of pulses in nonlinear regie is governed by the nonlinear Shrondinger (NLS) equation for which, their analytical solutions are calculated. An efficient ethod is used in the nuerical siulations under the nonlinear regie: the split-step fourier Method (SSFM). Finally, a brief study of the fundaental soliton is presented. Keywords Fibre optics, dispersion, pulse propagation, solitons. T I. INTRODUCTION he appearance of optical fibers drastically changed the paradig of data transission, because of two ain advantages that ake it particularly in relation to wires: The optical fibers are copletely iune to electroagnetic interference as they are ade of dielectric aterial capable of transitting pulses of light, which eans that data is not corrupted during transission. Optical fibers do not conduct electric current, so there will be no probles with electricity, as voltage instability probles or issues with spokes. The physical phenoenon known as total internal reflection is the responsible for light pulses to travel along the fiber by successive reflections []. Although since 87 it was known that light could describe a curved path within a aterial, only in 95 the Indian physicists Narinder Kapany could coplete its experients leading to the invention the optical fibers. However, only with the appearance of devices capable of converting electronic pulses into pulses of light, it was possible to transit inforation through the fiber [9]. There are two types of fiber: Multiode fibers with core diaeters greater than the single ode, allowing light to traverse the fiber through several paths, which are strongly affected by the effect of dispersion and, because of this, are ost often used for data traffic within walking distance. Single-ode fibers with core diaeters of the order of less than one strand of hair and let light travel in the interior of fiber by a single path, and less sensitive to the dispersion effect and are therefore ost often used for long distance counications [6], [8]. With this work, we hope to have contributed to better understanding of the effects to which light pulses propagate along the fiber are subject, both in the linear and nonlinear regies. II. MODAL THEORY Because of the difference between the refractive indices of the core and the cladding, it is possible to transit of light inside the fiber [3]. The core index of refraction of the core ust be always higher than the refractive index of the cladding. The odal theory can explain the existence of an angle of incidence that enables the phenoenon of total internal reflection. A. Fibers operated in the Linear Regie Optical fibers operated in the linear regie (onoodal regie), are the fibers that allow the use of only one light signal through the fiber. They have a saller radius and greater bandwidth due to lower dispersion. It is possible to define a axiu value for the angle of incidence with respect to the fiber axis, called acceptance angle axiu θ ax, iposing that only the rays entering the fiber at an angle of less than θ ax are propagated along it [3]. Let n be the refractive index in the core and n be the refractive index in the cladding. The paraeter which describes the ability to collect light in a fiber called a nuerical aperture (NA) which is related as θ ax follows

2 NA n sin n n n ax (.) where n is the refractive index of the outer ediu in which the fiber is inserted, usually the air, where n =. Δ is the dielectric contrast given by n n n (.) For better understanding the wave propagation inside a step index fiber with core radius a, one ay consider the schee depicted in Fig... Fig. Rays propagation in fiber []. The longitudinal coponent of the electric field in cylindrical coordinates has the following for Ez( r,, z, t) EF( r)exp( i )exp i( z t) (.3) where exp[i(βz-ωt)] is the tie and space dependence of onochroatic light traveling along the fiber axis, E is the aplitude of the field, is the aziuthal variation index, F (r) represents the radial odal eigenfunction, β is the longitudinal propagation constant, which is related to the free-space propagation constant as where k ( r) n ( r) k (.4) n for r a nr () (.5) n for r a and the propagation constant in vacuu is k c. Let k q in the cladding and k h in core. Since there is a surface wave guided by the cladding sheath, suffering total internal reflection at the core-cladding interface, under these conditions, one ust have q i. Introducing noralized (diensionless) wavenubers u ha and w a, and a noralized frequency v, their relation is given by v u w (.6) Moreover, introducing the noralized odal refractive index as u w ( k ) n b v v n n (.7) there is a iniu frequency of operation (noralized cutoff frequency v ) for which a certain ode starts to propagate in c the fiber. A special value is given by v. 448, which is the cutoff value for the second propagating ode, indicating whenever the fiber is operating below that value, the syste onoodal, or ultiodal if above, i.e., for singleode fiber, for B. Propagation Modes in Fiber v vc => ultiode fiber. c v v => Light, while propagating in the fiber, can be seen sa an electroagnetic phenoenon, and the whole propagating echanis, which can be described by the electroagnetic optical fields associated to it, is governed by Maxwell equations []. Generally speaking, odes can be classified into TE (Transverse Electric) odes, where there is no coponent of electric field in the direction of propagation (Ez = ), TM (Transverse Magnetic) odes, where there is no agnetic field coponent in the direction of propagation (Hz = ), and hybrid odes, for which there have both longitudinal field coponents, i.e., they have electric field (Ez ) and agnetic field (Hz ) coponents in the direction of propagation. These latter odes can be classified into HE n and EH n, where the paraeter refers to the aziuthal field variation and the paraeter n refers to the radial field variation []. In general, the surface waves guided by an optical fiber are hybrid odes, and their field coponents are governed by the Bessel functions of first (J) and second (K) species, satisfying the following equation odal where 4 u v R( u) S( u) v uv (.8) J ( u) K( w) R ( u) (.9a) uj ( u) wk ( w) J ( u) K( w) S( u) ( ) uj ( u) wk ( w) (.9b) The resolution of the Eq. (.8), while being very coplex, turns it to be necessary to use nuerical ethods. Only when the aziuthal index is zero ( = ) it is possible to propagate transversal TE n and TM n odes. According to the Gloge approxiation (i.e., for Δ << ), in which odes propagating in optical fibers are weakly guided, one ay take R(u)= S(u), which will cause the odal Eq.(.8) to be reduced to ( ) v R u uw (.) where the plus (+) sign corresponds to the odes EHn and the inus (-) sign corresponds to the odes HEn. Through the relationship between the Bessel functions and their derivatives under the Gloge approxiation, in which the c

3 3 aziuthal index zero, the Eq (.) for both odes (EHn HEn) reduces to J( u) K( w) uj ( u) wk ( w) (.) inferring that the odes EHn are actually TMn and odes HEn are odes TEn. Moreover for weakly guided fibers, these odes are alost linearly polarized, and therefore will be tered as LPpn. One should stress that the only ode capable of propagating in a single ode fiber syste (the only ode whose cutoff corresponds to v c =), called the fundaental ode, is the HE ode. Mosr propagating odes are degenerated odes, i.e., they have the sae cut-off frequency but with different field structure. In particular, HE n odes give rise to odes LP n Modes TE n, and TM n HE n are degenerate and give rise toodes LP n HE +,n odes and EH -,n with are also giving rise to degenerate odes LP n. Therefore, one has: EH n LP pn (with p=+); HE n LP pn (with p=-). Accordingly, one has the following odal for LP odes: J p( u) K p( w) u w J ( u) K ( w) (.) p For the fundaental ode (LP ), it reduces to uj ( u) K ( w) wj ( u) K ( w) p (.3) Figure. shows the variation of the noralized odal refractive index b as a function of noralized frequency v for the first six odes propagating in the fiber. Fig. - First six LP ode fiber: diagras b(v). One can observe a general increase in the noralized refractive index b odal with increasing noralized frequency v, and the greater the noralized frequency v, ore odes propagating in the fiber. It is also interesting to analyze the influence of the dielectric contrast on the dispersion curves b(v). as in Figure.3 Fig. 3 - Influence of dielectric contrast on the dispersion curve b (v). Noting that an increase in the dielectric contrast Δ causes an increase in dispersion curves b(v), one ay conclude that is possible to reduce effects of dispersion by using fibers with a saller contrast. III. PULSE PROPAGATION IN THE LINEAR REGIME Any pulse, while propagating in an optical fiber, especially in the linear regie, suffers the effect of tie dispersion, which causes its broadening and ay create interference between sybols, which can greatly liit the bandwidth of the signal to be transitted. Let one consider an optical pulse with a finite spectral width to be launched into the fiber. Each spectral coponent of the pulse, while traveling along the fiber, has a different group velocity which depends on its wavelength according to v c g with (3.) The dispersion due to the difference between the propagation velocities of different spectral coponents is called the Group Velocity Dispersion (GVD). There are two ain types of dispersion: the Interodal Dispersion (between the various odes of propagation) and intraodal dispersion (under the sae ode of propagation). Interodal dispersion is due to the various propagation odes travelling in the fiber. Each ode has a different group velocity for the sae wavelength and the pulse width at the output of the fiber will depend on the transission ties. The tie delay between the fastest (i.e., the fundaental ode) and the slowest ode (ode of higher order, depending on how any odes can propagate in this ultiode fiber) is responsible for the broadening of the pulse at the output of the fiber. In the case of the interodal dispersion or, as coonly referred to, the chroatic dispersion, the broadening of the pulse occurs within the sae ode, ainly because of the dependence of the refractive index of the aterial with the frequency. However, the chroatic dispersion is a consequence of cobined effect of two factors: the aterial dispersion and the waveguide dispersion [6], []. In the case of the aterial dispersion, the refractive index of the constituent aterials of the fiber has a non-linear variation with the wavelength dn ng n d (3.)

4 4 Waveguide dispersion follows fro the fact that, for a given ode, the energy distribution between the core and the cladding is a function of wavelength. Generally, 8% of optical power propagating in the fiber reains confined to the core while the reaining % propagate along the cladding at a speed greater than the core causing the broadening of the pulse at the fiber output [6]. In the neighborhood of the zero dispersion wavelength (ideal frequency band) the spectral coponents at different wavelengths have alost the sae propagation velocity, in other words, propagation delay is quite constant at different wavelengths. Fig 3. - Zero dispersion [5]. Discarding higher order dispersion effects, the paraeter which describes the fiber dispersion is given by c D (3.3) with the DGV coefficient β is given by vg v ( ) (3.4) g Let one take into consideration the electroagnetic spectru as depicted in Figure 3.. Fig 3. - Electroagnetic spectru. In the anoalous dispersion region, one has D λ > => β <. The group velocity increases with the frequency, i.e., the blue spectral coponents travel faster than those red (see Figure 3.3.). In noral dispersion region, one has D λ < => β >. The group velocity decreases with the frequency, i.e., in the red spectral coponents travel faster than the blue ones. Fig Noral dispersion [5]. For a standard silicon optical fiber, the zero dispersion wavelength occurs around 3 n. Around this value, there is no pulse broadening. For this reason, the ost current optical counication systes have been developed to take advantage is this characteristic []. A. Pulse Propagation Equation in linear Regie In this Section, the differential equation that governs the propagation of pulses in the linear regie along a single ode optical fiber with sall contrast (Δ << ), is derived ignoring the effect of higher order dispersion. Let A(,t) be the pulse at the fiber input at z =. Assuing that the pulse odulates a carrier, with an angular frequency ω, and that the electric field is linearly polarized along the x, one ay write E( x, y,, t) xe ˆ F( x, y) A(, t)exp( i t) (3.5) where F (x, y) is the odal function representing the transversal variation of the fields of the LP ode. With the help of the Fourier transfor (nuerically coputed through the FFT-Fast Fourier Transfor and Inverse Fast Fourier Transfor-IFFT), the pulse envelope at a generic point of the fiber is given by A ( z, t ) (, )exp( ) A z i t d (3.6) Thus the electric field in a generic point of the fiber is given by E( x, y, z, t) xe ˆ F( x, y) A( z, t)exp i( z t) (3.7) where ( ) (3.7a) In order to obtain A(z,t)) fro A(,t), we introduce the following GVD coefficients Fig Anoalous dispersion [5]. (3.8a)

5 5 where v g ( ) vg ( ) vg (3.8b) (3.8c) The frequency deviation over the carrier is given by (3.9) As Ω <<ω, it is reasonable to disregard the dispersion coefficients above β. Moreover, discarding the losses in the fiber, the linear propagation equation, which allows calculating A(z,t) fro A(,t), can be written by A( z, t) A A( z, t) i z t t (3.) This equation can be solved using a siple algorith which allows the coputation of pulse propagation along the fiber. This algorith is herein designated RIMF algorith and is applied in three steps: First step: To copute with A(, ) FFT A(, t) A(, t)exp( it) dt Second step: Then copute A( z, ) A(, )exp i ( ) z Third step: Finally A( z, t) IFFT A( z, ) A( z, )exp( it) d ( ) (3.) Because the pulses are usually narrowband Ω <<ω, A(z,t) is a slowly varying function in tie and oscillates with exp(-iωt). In the absence of the dispersive effects, the pulse propagate without distortion with a group delay For the nuerical solution of Eq. (3.) is useful to introduce the noralized (diensionless) space and tie variables z (3.4a) L D t z So equation (3.) can be rewritten as and the Fourier pair as (3.4b) A(, ) i A(, ) sgn( ) (3.5) A(, ) A(, )exp( i ) d (3.6a) A (, ) (, )exp( ) A i d (3.6b) where ψ is a noralized frequency given by ( ) (3.7) To apply the RIMF algorith to Eq. (3.5), one has First step: To calculate A(, ) FFT A(, ) Second step: Then i A(, ) (, )exp sgn( ) A Third step: Find A(, ) IFFT A(, ) Figure 3.5 illustrates the effect of the GVD on a secanthyperbolic shaped pulse, i.e., one whose initial shape is A ( ) A (, ) sec h ( ) (3.8) z (3.) g Since it is reasonable to neglect the influence of the higher order dispersion (β 3 =), it is possible to isolate the effect of the GDV in the pulse propagation. Defining τ as the characteristic pulse tie width, the dispersion length is introduced as L D (3.3) The dispersive effects in a optical link of length L are negligible only if L D >L. Fig Coparison of the pulse input and output fiber.

6 6 A(, ) FFT A(, t) A exp ic ( ic) (3.4) The pulse spectral intensity will be B. Chirp Fig Evolution of the pulse. Since the optical signals eitted by a laser source suffer fro chirp, it is useful to introduce a chirp paraeter C which quantifies the variation in the carrier frequency C c (3.9) where β c is Henry factor, responsible for the enlargeent of the spectral line. When considering the linear propagation of a Gaussian pulse with the Chirp [3], the initial pulse can be written as A(, ) A exp C ( C ) The spectral width Δω at /e axiu will be (3.5) C (3.6) One ay conclude that the larger the value of C, the larger the pulse spectral width. ic t A(, t) A exp where A represents the pulse aplitude. (3.) Fig Evolution of Gaussian pulse with C=. Fig Initial Gaussian pulse for C=. The frequency shift δω(z,t) caused by the existence of the chirp is z t z ( z, t) C ( C ) (3.) ( z) with ( z) ( z) (3.) and η(z) is the pulse broadening factor Fig Evolution of Gaussian pulse with C=. The pulse broadening factor η(z), introduced in Eq. (3.3), also shows that a pulse ay also suffer a tie copression as it propagates, as long as β C< [3]. Figure 3. shows the variation of the broadening factor with the travelled distance for a Gaussian pulse in the anoalous dispersion region (β <) z z ( z) sgn( ) C L D L D (3.3) When applying the first step of the RIMF algorith to the initial ipulse Gaussian, one has Fig 3. - Spatial evolution of spectral width pulse for different values of C.

7 7 The pulse broadening due to the dispersion is sensitive to the sharpness of the pulse. When considering pulses with sharper steep edges, the broadening is generally greater, as is the case Super-Gaussian pulse whose initial oentu is given by A(, t) A exp ic t (3.7) Fig Evolution of Super-Gaussian pulse with C=. This expression is quite identical to the Gaussian pulse, differing only in the value of the paraeter. For a Gaussian pulse is equal to. The edges of the pulses becoe increasingly steep as increases, as shown in Figure 3., when copared with figure 3.7. Fig 3. - Initial Super-Gaussian pulse for C=, to =3. The paraeter is related to the duration t r for which the intensity of the pulse increases fro % to 9% of its peak value t r (3.8) Eq. (3.8) shows that a pulse with a saller rise tie increases faster []. For this reason the super-gaussian pulses as well as expand faster than the Gaussian, also strongly distort its original shape, as shown below, when considering =3. Fig 3. - Super-Gaussian coparison pulse input and output fiber. C. Dispersion Copensation A technique used to copensate for or iniize the proble of dispersion is the use of dispersion copensating fibers (DCF Dispersion Copensation Fiber). This technique is to cobine optical fibers with different characteristics, such that the average GDV of the link becoes quite sall, whereas the GDV of each section ay be large [7]. In each two consecutive sections, there are two types of fiber. A section of greater length L, operating in the anoalous dispersion region (with β <), and a shorter section of length L, operating in the noral region (with β >). Nevertheless, the two GDV coefficients are quite different ( β β ). This technique takes advantage of the linear nature of the syste when considering the propagating optical pulse into two sections of the fiber, whose ipulse propagation equation is given by i A( L, ) A(, )exp ( L L ) i d (3.9) where L L L. The length of dispersion copensating fiber L is chosen so that L L L L (3.3) Pick-up β << β so that L is uch lower than L. When using this technique it is coon for L to be in the order of - k, while L is in the order of k. Provided that they fulfill the condition β L +β L =, A(L,τ)=A(,τ) thus recovering its original shape once every two consecutive sections, although pulse width ay change significantly in each section. As an exaple of application, we consider L = k, β = -, L = k e β = and the pulse input of the first fiber is given by A ( ) A (, ) sec h ( ) (3.3) Fig Evolution of Super-Gaussian pulse with C=. Fig 3. 5 Coparison input and output pulse in the st fiber.

8 8 The output pulse of the first fiber is the ipulse input of the second fiber (DCF). So, taking into account the shape of the pulse to propagate Fig 4. - Power of a Gaussian pulse [7]. Fig Coparison input and output pulse in the nd fiber. IV. PULSE PROPAGATION IN NONLINEAR REGIME Pulse propagation in nonlinear regie is affected by the optical Kerr effect. The propagation of ipulses in Nonlinear Dispersive Regie (NLDR) is governed by the SPM and the GDV siultaneously [4]. To disrupt the relative dielectric constant ( xy, ) (4.) where n( x, y) n (4.) Whatever the process that led to that disturbance, the new longitudinal propagation constant is given by (4.3) Introducing the paraeter nk n n (4.4) ref wherein is the effective area and r ef is the effective radius. The optical Kerr effect provides that for certain values of P( z, t) (4.5) n P( z, t ) is the transported power, that is related to the power of the internal fiber Pin() t and attenuation coefficient of the fiber as follows P( z, t) P ( t)exp( z) (4.6) in The nonlinear phase generated by Kerr effect will be L L L ( ) ( ) NL (, ) (4.7) t dz dz P z t dz And because the effective length is exp( L ) (4.7) Coes to (t NL ) ( t) P ( t) (4.8) NL in Noting then that the nonlinear phase depends only on P in (t), hence the nae of Self-Phase Modulation. The instantaneous frequency deviation of the local (over carrier) caused by the SPM pulse is dnl dpin () () t t (4.9) dt dt It has in front of the ipulse dp dt ( t), yielding a redshift. Siilarly the tail of pulse dp dt ( t), causes a blueshift. Since this is an in effect contrary to what happens in GDV to the anoalous dispersion, thus allows the propagation of pulses that retain their shape along the propagation [4]. A. Non-Linear Equation Shrodinger In non-linear regie, the equation of propagation of pulses is governed by NLS equation, that neglecting effects of higher order and losses is given by u u i sgn( ) u u (4.) If incident peak power P, observe the relationship P (4.) nl The propagation of pulses in a fiber of length L in siplistic ters [4], is governed by the nonlinear regie for L>L NL and the linear regie for L<L NL, for a grossly way, is possible to distinguish four regies of propagation in the fiber: Non-dispersive linear regie (NDLR), when L<L NL and L<L D ( disregard GDV and APM effects) Dispersive linear regie (DLR), when L<L NL and L>L D (disregard APM, only acts GDV) Non-dispersive nonlinear regie (NDNLR), when L>L NL and L<L D (disregard GDV, only acts APM) Non-dispersive nonlinear regie (NDNLR), when L>L NL and L>L D (acts GDV and APM siultaneously, thus allowing the propagation of solitons). B. Analytical Solutions of NLS Equation By liiting the analysis to just in case the solitons as they occur (anoalous dispersion zone, where sgn(β ) = -), coes to NLS equation u u i u u (4.) Which has the analytical solution for solitary waves u(, ) sec h( q) exp i (4.3) Where is the paraeter which sets both the aplitude and pulse width, q is pulse s center in relation and is initial phase (in ). in

9 9 Treating a wave with localized surrounding ( u(, ) is independent of ), apart fro that when tends to, u(, ) approaches. Give nae of fundaental soliton the canonical for of Eq. (4.3) by aking and q to phase zero. u(, ) sec h( )exp i (4.4) When considering the loss ( L D ) and 3rd dispersion order ( 3 6 ), variables N which represents non-linearity and D which represents the dispersion, are given as a function of N i u (4.8a) 3 i D sgn( ) (4.8b) 3 What is not entirely true, but it is considered the variation of u with negligible. When considering the incident pulse u ( ) u(, ) type, Eq. s (4.7) solution will be u(, ) exp ( N D ) u ( ) (4.9) Fig 4. - Evolution of fundaental soliton. In case of periodic waves, ust be taken into account that using the inverse of the dispersion or IST (inverse scattering transfor), shows that any incident pulse shape u ( ) u (, ) N sec h ( ) (4.5) Where N represents the soliton order. For soliton order N, unlike the fundaental soliton, pulses do not retain their shape along the propagation, shows instead an evolution periodic with period in real units represents (4.6) z L D Fig Evolution solution of the 3rd order. C. Nuerical Siulation NLS Equation: Split-Step Fourier Method The disclosed ethod for solving nonlinear equations propagating pulse is SSFM, which the optical apply the nuerical [4], is to separate the equation (4.), non-linear of dispersion. Equation (4.) can be written in a ore copact for u(, ) ( D N) u(, ) (4.7) Making u( h, ) exp h( N D ) u(, ) (4.) Eq. (4.) sets up an iterative schee of longitudinal step h, which enables the start ( ) at the end of the fiber ( L LLD L ). Trying to divide the total space propagation in sall sections eleentary length h. SSFM consists of two consecutive procedures, dividing the Eq. (4.) in v(, ) exp( hn ) u(, ) (4.a) u( h, ) exp( hd ) v(, ) (4.b) According to Eq. (4.3a), v (, ) is given by h v(, ) exp exp ih u(, ) u(, ) (4.) Using the Fourier transfors (st step RIMF algorith), then defines v(, ) FFT v(, ) v(, )exp( i ) d (4.3) And these conditions D, according to Eq. (4.8b) has been D becoes na algebraic operator i 3 D sgn( ) i (4.4) By taking into account Eq. (4.b) - nd step of the RIMF algorith u( h, ) exp( hd ) v(, ) h 3 exp i sgn( ) exp( ih ) v(, ) (4.5)

10 Each iteration to finally, apply 3rd step RIMF s algorith u( h, ) IFFT u( h, ) u( h, )exp( i ) d (4.6) Thus SSMF suarizes the application of RIMF algorith in case of non-linear equations, using an iterative schee longitudinal h, showing such a ethod since it is very effective to use longitudinal steps quite sall, because the saller step size, greater iterations nuber and increasing the nuber of iterations, the ore efficient it becoes ethod. Enhancing here that for siulation of all pulses in nonlinear regies, used SSFM (as was case in Figures and ). D. Characteristics of Fundaental Soliton Area of soliton does not depend on any paraeter characteristic fiber [4] A (4.7) Energy of the soliton is given by E s (4.8) Shown to be inversely proportional to the teporal width Power spectral density is directly related to energy, as follows Es S( ) d (4.9) [6] Andrade, M. A. (9). Modelização da propagação e sisteas de counicação óptica baseados na tecnica de ultiplexage por divisão no copriento de onda. Vila Real: Universidade de Trás-os-Montes e Alto Douro. [7] Dos Santos, N. M.-D. (). Métodos variacionais aplicados ao estudo das fibras ópticas e técnicas de copensação da dispersão. Lisboa: Instituto Superior Técnico. [8] fibra.no.sapo. (5 de Outubro de 4). Obtido e Janeiro de, de [9] Hardware, C. d. (s.d.). Clube do Hardware.co.br. Obtido e Novebro de, de [] paginas.fe.up.pt. (s.d.). Obtido e Setebro de, de fibrasopticas.pdf [] PUC - Rio, C. D. (s.d.). dbd.puc-rio.br. Obtido e Novebro de, de [] RNP, R. N. ( de Abril de ). rnp.br. Obtido e Janeiro de, de [3] UFRJ, G. d.-g. (s.d.). gta.ufrj.br. Obtido e Novebro de, de CoceitosCoposio.htl V. CONCLUSIONS In this paper we have shown that an optical pulse propagating along an optical fiber, suffers the effect of dispersion which broads the pulse at the output of the fiber. This is a liiting factor for the bandwidth of the transitted signal. In linear regie, the technique of using the dispersion copensating fibers every two consecutive sections, shows to be very efficient and to solve the effect of the dispersion. In non-linear regie, pulses are subject to AMF (caused by the optical Kerr effect). In the anoalous dispersion region, this effect is contrary to the DGV, thereby allowing the propagation of pulses very interesting because it does not changes its shape over propagation in the fiber. These pulses are called soliton. For the siulation pulse schee is nonlinear regie, we have used the SSFM, which has proved to be as efficient as higher the nuber of iterations used in the ethod. REFERENCES [] Agraval, G. P. (7). Maxwell s Equations. In Nonlinear Fiber Optics, 4ª edição. Burlington, USA: Acadeic Press. [] Agrawal, G. P. (7). Super-Gaussian Pulses. In Non Fiber Optics, 4ª Edição. Burlington, USA: Acadeic Press. [3] Paiva, C. (). Fotónica. Fibras Ópticas. [4] Paiva, C. R. (). Fotónica. Solitões E Fibras Ópticas. [5] Cartaxo, A. (). Sisteas de Telecounicações. Counicações Ópticas.

Anton Bourdine. 1. Introduction. and approximate propagation constants by following simple ratio (Equation (32.22) in [1]):

Anton Bourdine. 1. Introduction. and approximate propagation constants by following simple ratio (Equation (32.22) in [1]): Matheatical Probles in Engineering olue 5, Article ID 843, pages http://dx.doi.org/.55/5/843 Research Article Fast and Siple Method for Evaluation of Polarization Correction to Propagation Constant of

More information

In this chapter we will start the discussion on wave phenomena. We will study the following topics:

In this chapter we will start the discussion on wave phenomena. We will study the following topics: Chapter 16 Waves I In this chapter we will start the discussion on wave phenoena. We will study the following topics: Types of waves Aplitude, phase, frequency, period, propagation speed of a wave Mechanical

More information

The source of THz radiation based on dielectric waveguide excited by sequence of electron bunches

The source of THz radiation based on dielectric waveguide excited by sequence of electron bunches Journal of Physics: Conference Series PAPER OPEN ACCESS The source of THz radiation based on dielectric waveguide excited by sequence of electron bunches To cite this article: A M Altark and A D Kanareykin

More information

Electromagnetic Waves

Electromagnetic Waves Electroagnetic Waves Physics 4 Maxwell s Equations Maxwell s equations suarize the relationships between electric and agnetic fields. A ajor consequence of these equations is that an accelerating charge

More information

Lecture 3 Fiber Optical Communication Lecture 3, Slide 1

Lecture 3 Fiber Optical Communication Lecture 3, Slide 1 Lecture 3 Optical fibers as waveguides Maxwell s equations The wave equation Fiber modes Phase velocity, group velocity Dispersion Fiber Optical Communication Lecture 3, Slide 1 Maxwell s equations in

More information

Proc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES

Proc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Proc. of the IEEE/OES Seventh Working Conference on Current Measureent Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Belinda Lipa Codar Ocean Sensors 15 La Sandra Way, Portola Valley, CA 98 blipa@pogo.co

More information

PH 221-2A Fall Waves - I. Lectures Chapter 16 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition)

PH 221-2A Fall Waves - I. Lectures Chapter 16 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition) PH 1-A Fall 014 Waves - I Lectures 4-5 Chapter 16 (Halliday/Resnick/Walker, Fundaentals of Physics 9 th edition) 1 Chapter 16 Waves I In this chapter we will start the discussion on wave phenoena. We will

More information

ma x = -bv x + F rod.

ma x = -bv x + F rod. Notes on Dynaical Systes Dynaics is the study of change. The priary ingredients of a dynaical syste are its state and its rule of change (also soeties called the dynaic). Dynaical systes can be continuous

More information

General Properties of Radiation Detectors Supplements

General Properties of Radiation Detectors Supplements Phys. 649: Nuclear Techniques Physics Departent Yarouk University Chapter 4: General Properties of Radiation Detectors Suppleents Dr. Nidal M. Ershaidat Overview Phys. 649: Nuclear Techniques Physics Departent

More information

On Constant Power Water-filling

On Constant Power Water-filling On Constant Power Water-filling Wei Yu and John M. Cioffi Electrical Engineering Departent Stanford University, Stanford, CA94305, U.S.A. eails: {weiyu,cioffi}@stanford.edu Abstract This paper derives

More information

Data-Driven Imaging in Anisotropic Media

Data-Driven Imaging in Anisotropic Media 18 th World Conference on Non destructive Testing, 16- April 1, Durban, South Africa Data-Driven Iaging in Anisotropic Media Arno VOLKER 1 and Alan HUNTER 1 TNO Stieltjesweg 1, 6 AD, Delft, The Netherlands

More information

OPTI510R: Photonics. Khanh Kieu College of Optical Sciences, University of Arizona Meinel building R.626

OPTI510R: Photonics. Khanh Kieu College of Optical Sciences, University of Arizona Meinel building R.626 OPTI510R: Photonics Khanh Kieu College of Optical Sciences, University of Arizona kkieu@optics.arizona.edu Meinel building R.626 Announceents HW#3 is due next Wednesday, Feb. 21 st No class Monday Feb.

More information

Multi-Scale/Multi-Resolution: Wavelet Transform

Multi-Scale/Multi-Resolution: Wavelet Transform Multi-Scale/Multi-Resolution: Wavelet Transfor Proble with Fourier Fourier analysis -- breaks down a signal into constituent sinusoids of different frequencies. A serious drawback in transforing to the

More information

Electromagnetic scattering. Graduate Course Electrical Engineering (Communications) 1 st Semester, Sharif University of Technology

Electromagnetic scattering. Graduate Course Electrical Engineering (Communications) 1 st Semester, Sharif University of Technology Electroagnetic scattering Graduate Course Electrical Engineering (Counications) 1 st Seester, 1388-1389 Sharif University of Technology Contents of lecture 5 Contents of lecture 5: Scattering fro a conductive

More information

72. (30.2) Interaction between two parallel current carrying wires.

72. (30.2) Interaction between two parallel current carrying wires. 7. (3.) Interaction between two parallel current carrying wires. Two parallel wires carrying currents exert forces on each other. Each current produces a agnetic field in which the other current is placed.

More information

PHYS 102 Previous Exam Problems

PHYS 102 Previous Exam Problems PHYS 102 Previous Exa Probles CHAPTER 16 Waves Transverse waves on a string Power Interference of waves Standing waves Resonance on a string 1. The displaceent of a string carrying a traveling sinusoidal

More information

Figure 1: Equivalent electric (RC) circuit of a neurons membrane

Figure 1: Equivalent electric (RC) circuit of a neurons membrane Exercise: Leaky integrate and fire odel of neural spike generation This exercise investigates a siplified odel of how neurons spike in response to current inputs, one of the ost fundaental properties of

More information

Lecture 4 Fiber Optical Communication Lecture 4, Slide 1

Lecture 4 Fiber Optical Communication Lecture 4, Slide 1 ecture 4 Dispersion in single-mode fibers Material dispersion Waveguide dispersion imitations from dispersion Propagation equations Gaussian pulse broadening Bit-rate limitations Fiber losses Fiber Optical

More information

The Energy Flux Method for Reverberation: Modeling and Inversion

The Energy Flux Method for Reverberation: Modeling and Inversion DISTRIBUTION STATEMENT A Approved for public release; distribution is unliited The Energy Flux Method for Reverberation: Modeling and Inversion Ji-Xun Zhou School of Mechanical Engineering Georgia Institute

More information

Distributed-Feedback Lasers

Distributed-Feedback Lasers Distributed-Feedback Lasers Class: Integrated Photonic Devices Tie: Fri. 8:00a ~ 11:00a. Classroo: 資電 06 Lecturer: Prof. 李明昌 (Ming-Chang Lee) Wavelength Dependence of Bragg Reflections Free-Space Bragg

More information

Quiz 5 PRACTICE--Ch12.1, 13.1, 14.1

Quiz 5 PRACTICE--Ch12.1, 13.1, 14.1 Nae: Class: Date: ID: A Quiz 5 PRACTICE--Ch2., 3., 4. Multiple Choice Identify the choice that best copletes the stateent or answers the question.. A bea of light in air is incident at an angle of 35 to

More information

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either

More information

Ch 12: Variations on Backpropagation

Ch 12: Variations on Backpropagation Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith

More information

CHAPTER 1: INTRODUCTION

CHAPTER 1: INTRODUCTION CHAPTER 1: INTRODUCTION 1.1 SCOPE AND CONTENT Counications and sensing systes are ubiquitous. They are found in ilitary, industrial, edical, consuer, and scientific applications eploying radio frequency,

More information

Chapter 3: Dielectric Waveguides and Optical Fibers

Chapter 3: Dielectric Waveguides and Optical Fibers Chapter 3: Dielectric Waveguides and Optical Fibers Syetric planar dielectric slab waveguides Modal and waveguide dispersion in planar waveguides Step index fibers Nuerical aperture Dispersion in single

More information

Chapter 6 Aberrations

Chapter 6 Aberrations EE90F Chapter 6 Aberrations As we have seen, spherical lenses only obey Gaussian lens law in the paraxial approxiation. Deviations fro this ideal are called aberrations. F Rays toward the edge of the pupil

More information

ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER

ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER IEPC 003-0034 ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER A. Bober, M. Guelan Asher Space Research Institute, Technion-Israel Institute of Technology, 3000 Haifa, Israel

More information

Ph 20.3 Numerical Solution of Ordinary Differential Equations

Ph 20.3 Numerical Solution of Ordinary Differential Equations Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing

More information

FIBER OPTICS. Prof. R.K. Shevgaonkar. Department of Electrical Engineering. Indian Institute of Technology, Bombay. Lecture: 07

FIBER OPTICS. Prof. R.K. Shevgaonkar. Department of Electrical Engineering. Indian Institute of Technology, Bombay. Lecture: 07 FIBER OPTICS Prof. R.K. Shevgaonkar Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture: 07 Analysis of Wave-Model of Light Fiber Optics, Prof. R.K. Shevgaonkar, Dept. of

More information

Water a) 48 o b) 53 o c) 41.5 o d) 44 o. Glass. PHYSICS 223 Exam-2 NAME II III IV

Water a) 48 o b) 53 o c) 41.5 o d) 44 o. Glass. PHYSICS 223 Exam-2 NAME II III IV PHYSICS 3 Exa- NAME. In the figure shown, light travels fro aterial I, through three layers of other aterials with surfaces parallel to one another, and then back into another layer of aterial I. The refractions

More information

Reading from Young & Freedman: For this topic, read the introduction to chapter 25 and sections 25.1 to 25.3 & 25.6.

Reading from Young & Freedman: For this topic, read the introduction to chapter 25 and sections 25.1 to 25.3 & 25.6. PHY10 Electricity Topic 6 (Lectures 9 & 10) Electric Current and Resistance n this topic, we will cover: 1) Current in a conductor ) Resistivity 3) Resistance 4) Oh s Law 5) The Drude Model of conduction

More information

PH 222-2C Fall Electromagnetic Oscillations and Alternating Current. Lectures 18-19

PH 222-2C Fall Electromagnetic Oscillations and Alternating Current. Lectures 18-19 H - Fall 0 Electroagnetic Oscillations and Alternating urrent ectures 8-9 hapter 3 (Halliday/esnick/Walker, Fundaentals of hysics 8 th edition) hapter 3 Electroagnetic Oscillations and Alternating urrent

More information

Chapter 28: Alternating Current

Chapter 28: Alternating Current hapter 8: Alternating urrent Phasors and Alternating urrents Alternating current (A current) urrent which varies sinusoidally in tie is called alternating current (A) as opposed to direct current (D).

More information

Measuring orbital angular momentum superpositions of light by mode transformation

Measuring orbital angular momentum superpositions of light by mode transformation CHAPTER 7 Measuring orbital angular oentu superpositions of light by ode transforation In chapter 6 we reported on a ethod for easuring orbital angular oentu (OAM) states of light based on the transforation

More information

In this chapter we will study sound waves and concentrate on the following topics:

In this chapter we will study sound waves and concentrate on the following topics: Chapter 17 Waves II In this chapter we will study sound waves and concentrate on the following topics: Speed of sound waves Relation between displaceent and pressure aplitude Interference of sound waves

More information

A DESIGN GUIDE OF DOUBLE-LAYER CELLULAR CLADDINGS FOR BLAST ALLEVIATION

A DESIGN GUIDE OF DOUBLE-LAYER CELLULAR CLADDINGS FOR BLAST ALLEVIATION International Journal of Aerospace and Lightweight Structures Vol. 3, No. 1 (2013) 109 133 c Research Publishing Services DOI: 10.3850/S201042862013000550 A DESIGN GUIDE OF DOUBLE-LAYER CELLULAR CLADDINGS

More information

Now multiply the left-hand-side by ω and the right-hand side by dδ/dt (recall ω= dδ/dt) to get:

Now multiply the left-hand-side by ω and the right-hand side by dδ/dt (recall ω= dδ/dt) to get: Equal Area Criterion.0 Developent of equal area criterion As in previous notes, all powers are in per-unit. I want to show you the equal area criterion a little differently than the book does it. Let s

More information

OSCILLATIONS AND WAVES

OSCILLATIONS AND WAVES OSCILLATIONS AND WAVES OSCILLATION IS AN EXAMPLE OF PERIODIC MOTION No stories this tie, we are going to get straight to the topic. We say that an event is Periodic in nature when it repeats itself in

More information

Causality and the Kramers Kronig relations

Causality and the Kramers Kronig relations Causality and the Kraers Kronig relations Causality describes the teporal relationship between cause and effect. A bell rings after you strike it, not before you strike it. This eans that the function

More information

Particle interactions with spatially localized wavepackets

Particle interactions with spatially localized wavepackets PSFC/JA-10-36 Particle interactions with spatially localized wavepackets Y. Koinis, a K. Hizanidis, a and A.K. Ra October 010 Plasa Science and Fusion Center, Massachusetts Institute of Technology Cabridge,

More information

Actuators & Mechanisms Actuator sizing

Actuators & Mechanisms Actuator sizing Course Code: MDP 454, Course Nae:, Second Seester 2014 Actuators & Mechaniss Actuator sizing Contents - Modelling of Mechanical Syste - Mechaniss and Drives The study of Mechatronics systes can be divided

More information

ANALYSIS ON RESPONSE OF DYNAMIC SYSTEMS TO PULSE SEQUENCES EXCITATION

ANALYSIS ON RESPONSE OF DYNAMIC SYSTEMS TO PULSE SEQUENCES EXCITATION The 4 th World Conference on Earthquake Engineering October -7, 8, Beijing, China ANALYSIS ON RESPONSE OF DYNAMIC SYSTEMS TO PULSE SEQUENCES EXCITATION S. Li C.H. Zhai L.L. Xie Ph. D. Student, School of

More information

Sensorless Control of Induction Motor Drive Using SVPWM - MRAS Speed Observer

Sensorless Control of Induction Motor Drive Using SVPWM - MRAS Speed Observer Journal of Eerging Trends in Engineering and Applied Sciences (JETEAS) 2 (3): 509-513 Journal Scholarlink of Eerging Research Trends Institute in Engineering Journals, 2011 and Applied (ISSN: 2141-7016)

More information

NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT

NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT PACS REFERENCE: 43.5.LJ Krister Larsson Departent of Applied Acoustics Chalers University of Technology SE-412 96 Sweden Tel: +46 ()31 772 22 Fax: +46 ()31

More information

= T. Oscillations and Waves. Example of an Oscillating System IB 12 IB 12

= T. Oscillations and Waves. Example of an Oscillating System IB 12 IB 12 Oscillation: the vibration of an object Oscillations and Waves Eaple of an Oscillating Syste A ass oscillates on a horizontal spring without friction as shown below. At each position, analyze its displaceent,

More information

III. Quantization of electromagnetic field

III. Quantization of electromagnetic field III. Quantization of electroagnetic field Using the fraework presented in the previous chapter, this chapter describes lightwave in ters of quantu echanics. First, how to write a physical quantity operator

More information

Physically Based Modeling CS Notes Spring 1997 Particle Collision and Contact

Physically Based Modeling CS Notes Spring 1997 Particle Collision and Contact Physically Based Modeling CS 15-863 Notes Spring 1997 Particle Collision and Contact 1 Collisions with Springs Suppose we wanted to ipleent a particle siulator with a floor : a solid horizontal plane which

More information

IN modern society that various systems have become more

IN modern society that various systems have become more Developent of Reliability Function in -Coponent Standby Redundant Syste with Priority Based on Maxiu Entropy Principle Ryosuke Hirata, Ikuo Arizono, Ryosuke Toohiro, Satoshi Oigawa, and Yasuhiko Takeoto

More information

Study Committee B5 Colloquium 2005 September Calgary, CANADA

Study Committee B5 Colloquium 2005 September Calgary, CANADA 36 Study oittee B olloquiu Septeber 4-6 algary, ND ero Sequence urrent opensation for Distance Protection applied to Series opensated Parallel Lines TKHRO KSE* PHL G BEUMONT Toshiba nternational (Europe

More information

Q ESTIMATION WITHIN A FORMATION PROGRAM q_estimation

Q ESTIMATION WITHIN A FORMATION PROGRAM q_estimation Foration Attributes Progra q_estiation Q ESTIMATION WITHIN A FOMATION POGAM q_estiation Estiating Q between stratal slices Progra q_estiation estiate seisic attenuation (1/Q) on coplex stratal slices using

More information

Electric Power System Transient Stability Analysis Methods

Electric Power System Transient Stability Analysis Methods 1 Electric Power Syste Transient Stability Analysis Methods João Pedro de Carvalho Mateus, IST Abstract In this paper are presented the state of the art of Electric Power Syste transient stability analysis

More information

Question 1. [14 Marks]

Question 1. [14 Marks] 6 Question 1. [14 Marks] R r T! A string is attached to the dru (radius r) of a spool (radius R) as shown in side and end views here. (A spool is device for storing string, thread etc.) A tension T is

More information

Chapter 2 General Properties of Radiation Detectors

Chapter 2 General Properties of Radiation Detectors Med Phys 4RA3, 4RB3/6R3 Radioisotopes and Radiation Methodology -1 Chapter General Properties of Radiation Detectors Ionizing radiation is ost coonly detected by the charge created when radiation interacts

More information

Spine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area

Spine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area Proceedings of the 006 WSEAS/IASME International Conference on Heat and Mass Transfer, Miai, Florida, USA, January 18-0, 006 (pp13-18) Spine Fin Efficiency A Three Sided Pyraidal Fin of Equilateral Triangular

More information

26 Impulse and Momentum

26 Impulse and Momentum 6 Ipulse and Moentu First, a Few More Words on Work and Energy, for Coparison Purposes Iagine a gigantic air hockey table with a whole bunch of pucks of various asses, none of which experiences any friction

More information

PH 221-3A Fall Waves - II. Lectures Chapter 17 (Halliday/Resnick/Walker, Fundamentals of Physics 8 th edition)

PH 221-3A Fall Waves - II. Lectures Chapter 17 (Halliday/Resnick/Walker, Fundamentals of Physics 8 th edition) PH 221-3A Fall 2010 Waves - II Lectures 27-28 Chapter 17 (Halliday/Resnick/Walker, Fundaentals of Physics 8 th edition) 1 Chapter 17 Waves II In this chapter we will study sound waves and concentrate on

More information

COS 424: Interacting with Data. Written Exercises

COS 424: Interacting with Data. Written Exercises COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well

More information

UNCERTAINTIES IN THE APPLICATION OF ATMOSPHERIC AND ALTITUDE CORRECTIONS AS RECOMMENDED IN IEC STANDARDS

UNCERTAINTIES IN THE APPLICATION OF ATMOSPHERIC AND ALTITUDE CORRECTIONS AS RECOMMENDED IN IEC STANDARDS Paper Published on the16th International Syposiu on High Voltage Engineering, Cape Town, South Africa, 2009 UNCERTAINTIES IN THE APPLICATION OF ATMOSPHERIC AND ALTITUDE CORRECTIONS AS RECOMMENDED IN IEC

More information

DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS *

DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS * IJST, Transactions of Mechanical Engineering, Vol. 39, No. M1, pp 89-100 Printed in The Islaic Republic of Iran, 2015 Shira University DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS

More information

Projectile Motion with Air Resistance (Numerical Modeling, Euler s Method)

Projectile Motion with Air Resistance (Numerical Modeling, Euler s Method) Projectile Motion with Air Resistance (Nuerical Modeling, Euler s Method) Theory Euler s ethod is a siple way to approxiate the solution of ordinary differential equations (ode s) nuerically. Specifically,

More information

Hyperbolic Horn Helical Mass Spectrometer (3HMS) James G. Hagerman Hagerman Technology LLC & Pacific Environmental Technologies April 2005

Hyperbolic Horn Helical Mass Spectrometer (3HMS) James G. Hagerman Hagerman Technology LLC & Pacific Environmental Technologies April 2005 Hyperbolic Horn Helical Mass Spectroeter (3HMS) Jaes G Hageran Hageran Technology LLC & Pacific Environental Technologies April 5 ABSTRACT This paper describes a new type of ass filter based on the REFIMS

More information

An Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period

An Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period An Approxiate Model for the Theoretical Prediction of the Velocity... 77 Central European Journal of Energetic Materials, 205, 2(), 77-88 ISSN 2353-843 An Approxiate Model for the Theoretical Prediction

More information

Fourier Series Summary (From Salivahanan et al, 2002)

Fourier Series Summary (From Salivahanan et al, 2002) Fourier Series Suary (Fro Salivahanan et al, ) A periodic continuous signal f(t), - < t

More information

Electromagnetic fields modeling of power line communication (PLC)

Electromagnetic fields modeling of power line communication (PLC) Electroagnetic fields odeling of power line counication (PLC) Wei Weiqi UROP 3 School of Electrical and Electronic Engineering Nanyang echnological University E-ail: 4794486@ntu.edu.sg Keyword: power line

More information

This model assumes that the probability of a gap has size i is proportional to 1/i. i.e., i log m e. j=1. E[gap size] = i P r(i) = N f t.

This model assumes that the probability of a gap has size i is proportional to 1/i. i.e., i log m e. j=1. E[gap size] = i P r(i) = N f t. CS 493: Algoriths for Massive Data Sets Feb 2, 2002 Local Models, Bloo Filter Scribe: Qin Lv Local Models In global odels, every inverted file entry is copressed with the sae odel. This work wells when

More information

Physics 41 Homework #2 Chapter 16. fa. Here v is the speed of the wave. 16. The speed of a wave on a massless string would be infinite!

Physics 41 Homework #2 Chapter 16. fa. Here v is the speed of the wave. 16. The speed of a wave on a massless string would be infinite! Physics 41 Hoewor # Chapter 1 Serway 7 th Conceptual: Q: 3,, 8, 11, 1, Probles P: 1, 3, 5, 9, 1, 5, 31, 35, 38, 4, 5, 57 Conceptual 3. (i) d=e, f, c, b, a (ii) Since, the saller the (the coefficient of

More information

Modulation of Harmonic Emission Spectra from Intense Laser-Plasma Interactions

Modulation of Harmonic Emission Spectra from Intense Laser-Plasma Interactions Modulation of Haronic Eission Spectra fro Intense Laser-Plasa Interactions T.J.M. Boyd and R. Ondarza-Rovira 2 Centre for Physics, University of Essex, Wivenhoe Park, Colchester CO4 3SQ, U.K. 2 ININ, A.P.

More information

3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme

3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme P-8 3D acoustic wave odeling with a tie-space doain dispersion-relation-based Finite-difference schee Yang Liu * and rinal K. Sen State Key Laboratory of Petroleu Resource and Prospecting (China University

More information

Optical Properties of Plasmas of High-Z Elements

Optical Properties of Plasmas of High-Z Elements Forschungszentru Karlsruhe Techni und Uwelt Wissenschaftlishe Berichte FZK Optical Properties of Plasas of High-Z Eleents V.Tolach 1, G.Miloshevsy 1, H.Würz Project Kernfusion 1 Heat and Mass Transfer

More information

University of Bath DOI: /OE Publication date: Link to publication

University of Bath DOI: /OE Publication date: Link to publication Citation for published version: Chen, L & Bird, DM 2011, 'Guidance in Kagoe-like photonic crystal fibres II: perturbation theory for a realistic fibre structure' Optics Express, vol. 19, no. 7, pp. 6957-6968.

More information

Department of Physics Preliminary Exam January 3 6, 2006

Department of Physics Preliminary Exam January 3 6, 2006 Departent of Physics Preliinary Exa January 3 6, 2006 Day 1: Classical Mechanics Tuesday, January 3, 2006 9:00 a.. 12:00 p.. Instructions: 1. Write the answer to each question on a separate sheet of paper.

More information

Step index planar waveguide

Step index planar waveguide N. Dubreuil S. Lebrun Exam without document Pocket calculator permitted Duration of the exam: 2 hours The exam takes the form of a multiple choice test. Annexes are given at the end of the text. **********************************************************************************

More information

Explicit Approximate Solution for Finding the. Natural Frequency of the Motion of Pendulum. by Using the HAM

Explicit Approximate Solution for Finding the. Natural Frequency of the Motion of Pendulum. by Using the HAM Applied Matheatical Sciences Vol. 3 9 no. 1 13-13 Explicit Approxiate Solution for Finding the Natural Frequency of the Motion of Pendulu by Using the HAM Ahad Doosthoseini * Mechanical Engineering Departent

More information

MODIFIED SERIES RESISTANCE MODEL - DETERMINATION OF MEAN CONCENTRATION BY INTEGRAL TRANSFORMATION

MODIFIED SERIES RESISTANCE MODEL - DETERMINATION OF MEAN CONCENTRATION BY INTEGRAL TRANSFORMATION MODIFIED SERIES RESISTANCE MODEL - DETERMINATION OF MEAN CONCENTRATION BY INTEGRAL TRANSFORMATION A. L. Venezuela a, R. F. Cantão a, R. S. Ongaratto b, and R. N. Haneda c a Universidade Federal de São

More information

P235 Midterm Examination Prof. Cline

P235 Midterm Examination Prof. Cline P235 Mier Exaination Prof. Cline THIS IS A CLOSED BOOK EXAMINATION. Do all parts of all four questions. Show all steps to get full credit. 7:00-10.00p, 30 October 2009 1:(20pts) Consider a rocket fired

More information

Plasma Modulation of Harmonic Emission Spectra from Laser-Dense Plasma Interactions

Plasma Modulation of Harmonic Emission Spectra from Laser-Dense Plasma Interactions Plasa Modulation of Haronic Eission Spectra fro Laser-Dense Plasa Interactions R. Ondarza-Rovira and T.J.M. Boyd 2 ININ, A.P. 8-027, México 80, Distrito Federal, Mexico; ondarza@nuclear.inin.x 2 Centre

More information

Non-Parametric Non-Line-of-Sight Identification 1

Non-Parametric Non-Line-of-Sight Identification 1 Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,

More information

Lecture #8-3 Oscillations, Simple Harmonic Motion

Lecture #8-3 Oscillations, Simple Harmonic Motion Lecture #8-3 Oscillations Siple Haronic Motion So far we have considered two basic types of otion: translation and rotation. But these are not the only two types of otion we can observe in every day life.

More information

Kinematics and dynamics, a computational approach

Kinematics and dynamics, a computational approach Kineatics and dynaics, a coputational approach We begin the discussion of nuerical approaches to echanics with the definition for the velocity r r ( t t) r ( t) v( t) li li or r( t t) r( t) v( t) t for

More information

The accelerated expansion of the universe is explained by quantum field theory.

The accelerated expansion of the universe is explained by quantum field theory. The accelerated expansion of the universe is explained by quantu field theory. Abstract. Forulas describing interactions, in fact, use the liiting speed of inforation transfer, and not the speed of light.

More information

Numerical Solution of the MRLW Equation Using Finite Difference Method. 1 Introduction

Numerical Solution of the MRLW Equation Using Finite Difference Method. 1 Introduction ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.1401 No.3,pp.355-361 Nuerical Solution of the MRLW Equation Using Finite Difference Method Pınar Keskin, Dursun Irk

More information

OBJECTIVES INTRODUCTION

OBJECTIVES INTRODUCTION M7 Chapter 3 Section 1 OBJECTIVES Suarize data using easures of central tendency, such as the ean, edian, ode, and idrange. Describe data using the easures of variation, such as the range, variance, and

More information

Mutual capacitor and its applications

Mutual capacitor and its applications Mutual capacitor and its applications Chun Li, Jason Li, Jieing Li CALSON Technologies, Toronto, Canada E-ail: calandli@yahoo.ca Published in The Journal of Engineering; Received on 27th October 2013;

More information

Chapter 1: Basics of Vibrations for Simple Mechanical Systems

Chapter 1: Basics of Vibrations for Simple Mechanical Systems Chapter 1: Basics of Vibrations for Siple Mechanical Systes Introduction: The fundaentals of Sound and Vibrations are part of the broader field of echanics, with strong connections to classical echanics,

More information

Optical Fiber Signal Degradation

Optical Fiber Signal Degradation Optical Fiber Signal Degradation Effects Pulse Spreading Dispersion (Distortion) Causes the optical pulses to broaden as they travel along a fiber Overlap between neighboring pulses creates errors Resulting

More information

Analytical solution for the electric potential in arbitrary anisotropic layered media applying the set of Hankel transforms of integer order

Analytical solution for the electric potential in arbitrary anisotropic layered media applying the set of Hankel transforms of integer order Geophysical Prospecting, 2006, 54, 651 661 Analytical solution for the electric potential in arbitrary anisotropic layered edia applying the set of Hankel transfors of integer order E. Pervago 1,A.Mousatov

More information

2.141 Modeling and Simulation of Dynamic Systems Assignment #2

2.141 Modeling and Simulation of Dynamic Systems Assignment #2 2.141 Modeling and Siulation of Dynaic Systes Assignent #2 Out: Wednesday Septeber 20, 2006 Due: Wednesday October 4, 2006 Proble 1 The sketch shows a highly siplified diagra of a dry-dock used in ship

More information

( ). One set of terms has a ω in

( ). One set of terms has a ω in Laptag Class Notes W. Gekelan Cold Plasa Dispersion relation Suer Let us go back to a single particle and see how it behaves in a high frequency electric field. We will use the force equation and Maxwell

More information

SOLUTIONS. PROBLEM 1. The Hamiltonian of the particle in the gravitational field can be written as, x 0, + U(x), U(x) =

SOLUTIONS. PROBLEM 1. The Hamiltonian of the particle in the gravitational field can be written as, x 0, + U(x), U(x) = SOLUTIONS PROBLEM 1. The Hailtonian of the particle in the gravitational field can be written as { Ĥ = ˆp2, x 0, + U(x), U(x) = (1) 2 gx, x > 0. The siplest estiate coes fro the uncertainty relation. If

More information

Physics 221A: HW3 solutions

Physics 221A: HW3 solutions Physics 22A: HW3 solutions October 22, 202. a) It will help to start things off by doing soe gaussian integrals. Let x be a real vector of length, and let s copute dxe 2 xt Ax, where A is soe real atrix.

More information

Scattering by a Blade on a Metallic Plane

Scattering by a Blade on a Metallic Plane UEMG #121474 Scattering by a Blade on a Metallic Plane DANILO ERRICOLO PIERGIORGIO L. E. USLENGHI BADRIA ELNOUR FRANCESCA MIOC QUERY SHEET QY: Au: Please check all artwork throughout article for correctness.

More information

EE5900 Spring Lecture 4 IC interconnect modeling methods Zhuo Feng

EE5900 Spring Lecture 4 IC interconnect modeling methods Zhuo Feng EE59 Spring Parallel LSI AD Algoriths Lecture I interconnect odeling ethods Zhuo Feng. Z. Feng MTU EE59 So far we ve considered only tie doain analyses We ll soon see that it is soeties preferable to odel

More information

In this chapter, we consider several graph-theoretic and probabilistic models

In this chapter, we consider several graph-theoretic and probabilistic models THREE ONE GRAPH-THEORETIC AND STATISTICAL MODELS 3.1 INTRODUCTION In this chapter, we consider several graph-theoretic and probabilistic odels for a social network, which we do under different assuptions

More information

m A 1 m mgd k m v ( C) AP Physics Multiple Choice Practice Oscillations

m A 1 m mgd k m v ( C) AP Physics Multiple Choice Practice Oscillations P Physics Multiple Choice Practice Oscillations. ass, attached to a horizontal assless spring with spring constant, is set into siple haronic otion. Its axiu displaceent fro its equilibriu position is.

More information

OPTI510R: Photonics. Khanh Kieu College of Optical Sciences, University of Arizona Meinel building R.626

OPTI510R: Photonics. Khanh Kieu College of Optical Sciences, University of Arizona Meinel building R.626 OPTI510R: Photonics Khanh Kieu College of Optical Sciences, University of Arizona kkieu@optics.arizona.edu Meinel building R.66 1 Announceents Hoework # is due today, HW#3 is assigned due Feb. 1 st No

More information

Supporting Information for Supression of Auger Processes in Confined Structures

Supporting Information for Supression of Auger Processes in Confined Structures Supporting Inforation for Supression of Auger Processes in Confined Structures George E. Cragg and Alexander. Efros Naval Research aboratory, Washington, DC 20375, USA 1 Solution of the Coupled, Two-band

More information

Physics 2107 Oscillations using Springs Experiment 2

Physics 2107 Oscillations using Springs Experiment 2 PY07 Oscillations using Springs Experient Physics 07 Oscillations using Springs Experient Prelab Read the following bacground/setup and ensure you are failiar with the concepts and theory required for

More information

Intelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines

Intelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes

More information

Time-of-flight Identification of Ions in CESR and ERL

Time-of-flight Identification of Ions in CESR and ERL Tie-of-flight Identification of Ions in CESR and ERL Eric Edwards Departent of Physics, University of Alabaa, Tuscaloosa, AL, 35486 (Dated: August 8, 2008) The accuulation of ion densities in the bea pipe

More information

Generalized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials

Generalized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials Copyright 7 Tech Science Press CMC vol.53 no.4 pp.37-34 7 Generalized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials S.D. Akbarov and M. Negin 3 Abstract: Dispersion of

More information