A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first ...A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.展开更多
A new numerical technique based on the wavelet derivative operator is presented as an alternative to BPM to study the integrated optical waveguide. The wavelet derivative operator is used instead of FFT/IFFT or finite...A new numerical technique based on the wavelet derivative operator is presented as an alternative to BPM to study the integrated optical waveguide. The wavelet derivative operator is used instead of FFT/IFFT or finite difference to calculate the derivatives of the transverse variable in the Helmholtz equation. Results of numerically simulating the injected field at z =0 are exhibited with Gaussian distribution in transverse direction propagating through the two dimensional waveguides (with linear and/or nonlinear refractive index) , which are similar to those in the related publications. Consequently it is efficient and needs not absorbing boundary by introducing the interpolation operator during calculating the wavelet derivative operator. The iterative process needs fewer steps to be stable. Also, when the light wave meets the changes of mediums, the wavelet derivative operator has the adaptive property to adjust those changes at the boundaries.展开更多
A modified wide-angle beam propagation based on the Douglas operators is presented.The truncation error in the modified wide-angle beam propagation is reduced to o (Δ x ) 4 in the transverse direction nearly wi...A modified wide-angle beam propagation based on the Douglas operators is presented.The truncation error in the modified wide-angle beam propagation is reduced to o (Δ x ) 4 in the transverse direction nearly without any increase of the computation time,whereas the error in the ordinary wide-angle beam propagation method is typically o (Δ x ) 2.With trivial programming changes,the accuracy is higher,especially in wide-angle propagation.展开更多
A beam propagation method based on the Galerkin method with Hermite-Gauss basis functions for studying optical field propagation in weakly guiding dielectric structures is described. The selected basis functions natur...A beam propagation method based on the Galerkin method with Hermite-Gauss basis functions for studying optical field propagation in weakly guiding dielectric structures is described. The selected basis functions naturally satisfy the required boundary conditions at infinity so that the boundary truncation is avoided. The paraxial propagation equation is converted into a set of first-order ordinary differential equations, which are solved by means of standard numerical library routines. Besides, the calculation is efficient due to its small resulted matrix. The evolution of the injected field and its normalized power along the propagation distance in an asymmetric slab waveguide and directional coupler are presented, and the solutions are good agreement with those obtained by finite difference BPM, which tests the validity of the present approach.展开更多
Based on the dynamic governing equation of propagating buckle on a beam on a nonlinear elastic foundation, this paper deals with an important problem of buckle arrest by combining the FEM with a time integration techn...Based on the dynamic governing equation of propagating buckle on a beam on a nonlinear elastic foundation, this paper deals with an important problem of buckle arrest by combining the FEM with a time integration technique. A new conclusion completely different from that by the quasi-static analysis about the buckle arrestor design is drawn. This shows that the inertia of the beam cannot be ignored in the analysis under consideration, especially when the buckle propagation is suddenly stopped by the arrestors.展开更多
For exact estimation of efficiency of a buckle arrestor, it is necessary to take the effect of structural inertia into account in the analysis of buckle propagation on elastic structures after meeting arrestors. Under...For exact estimation of efficiency of a buckle arrestor, it is necessary to take the effect of structural inertia into account in the analysis of buckle propagation on elastic structures after meeting arrestors. Under this consideration, this paper deals with the dynamics of buckle arrest and its numerical simulation on the basis of the beam system model used by Chater and Hutchinson (1983). The FEM combined with an improving are-length control method is adopted to solve the dynamic equations describing the arresting of buckle propagation. A new group of parameters for arrestor design which differs greatly from that by the quasi-static analysis is obtained. The present results support the conclusion that the inertia of the beam cannot be neglected in such analysis.展开更多
A modified alternating direction implicit approach is proposed to discretize the three-dimensional full-vectorial beam propagation method （3D-FV-BPM） formulation along the longitudinal direction. The cross-coupling ...A modified alternating direction implicit approach is proposed to discretize the three-dimensional full-vectorial beam propagation method （3D-FV-BPM） formulation along the longitudinal direction. The cross-coupling terms （CCTs） are neglected at the first substep, and then double used at the second substep. The order of two substeps is reversed for each transverse electric field component so that the CCTs are always expressed in an implicit form, thus the calculation is efficient and stable. Based on the multinomial interpolation, a universal finite difference scheme with a high accuracy is developed to approximate the 3D-FV-BPM formulation along the transverse directions, in which the discontinuities of the normal components of the electric field across the abrupt dielectric interfaces are taken into account and can be applied to both uniform and non-uniform grids. The corresponding imaginary-distance procedure is first applied to a buried rectangular and a GaAs-based deeply-etched rib waveguide. The field patterns and the normalized propagation constants of the fundamental and the first order modes are presented and the hybrid nature of the full-vectorial guided-modes is demonstrated, which shows the validity and utility of the present approach. Then the modal characteristics of the deeply- and shallow-etched rib waveguides based on the InGaAsp/InGaAsP strained multiple quantum wells in InP substrate are investigated in detail. The results are necessary for modeling and the design of the planar lightwave circuits or photonic integrated circuits based on these waveguides.展开更多
The propagation of an optical vortex in a hexagonally arranged single mode multicore fiber structure is investigated for possible generation of additional vortices and their spread dynamics. Fields are separated into ...The propagation of an optical vortex in a hexagonally arranged single mode multicore fiber structure is investigated for possible generation of additional vortices and their spread dynamics. Fields are separated into a slowly varying paraxial envelope and a rapidly changing exponential component. Solutions are derived from the paraxial inhomogeneous Schrodinger equation in two dimensions along with the index of refraction of the proposed structure. Numerical analyses are based on the beam propagation method and transparent boundary conditions in matrix form with different parameters to represent the intensity and phase of all derived fields. Vortices are numerically identified by their points of zero intensity and their phase change or polarity. The optical interferogram with a plane wave reference is also employed to distinguish the dislocation points in the transverse directions of the propagating fields.展开更多
In this paper the analytical solutions of the impact of a particle on Timoshenko beams with four kinds of different boundary conditions are obtained according to Navier's idea, which is further developed. The init...In this paper the analytical solutions of the impact of a particle on Timoshenko beams with four kinds of different boundary conditions are obtained according to Navier's idea, which is further developed. The initial values of the impact forces are exactly determined by the momentum conservation law. The propagation of the longitudinal and transverse waves along the beam, especially, the effects of boundary conditions on the characteristics of the reflected waves, are investigated in detail. Some results are compared with those by MSC/NASTRAN.展开更多
A novel three-dimensional beam propagation method (BPM) based on the variable transformed Galerkin's method is introduced for simulating optical field propagation in three-dimensional dielectric structures. The in...A novel three-dimensional beam propagation method (BPM) based on the variable transformed Galerkin's method is introduced for simulating optical field propagation in three-dimensional dielectric structures. The infinite Cartesian x-y plane is mapped into a unit square by a tangent-type function transformation. Consequently, the infinite region problem is converted into the finite region problem. Thus, the boundary truncation is eliminated and the calculation accuracy is promoted. The three-dimensional BPM basic equation is reduced to a set of first-order ordinary differential equations through sinusoidal basis function, which fits arbitrary cladding optical waveguide, then direct solution of the resulting equations by means of the Runge-Kutta method. In addition, the calculation is efficient due to the small matrix derived from the present technique. Both z-invariant and z-variant examples are considered to test both the accuracy and utility of this approach.展开更多
Based on the vectorial structure of electromagnetic beam and the method of stationary phase, the analytical structure of Hermite Gaussian beam in far field is presented. The structural energy flux distributions are al...Based on the vectorial structure of electromagnetic beam and the method of stationary phase, the analytical structure of Hermite Gaussian beam in far field is presented. The structural energy flux distributions are also investigated in the far field. The structural pictures of some Hermite Gaussian beams are depicted in the far field. As the structure of Hermite Gaussian beam is dominated by the transverse mode numbers and the initial transverse Gaussian half width, it is more complex than that of Gaussian beam. The ratios of the structural energy fluxes to the whole energy flux are independent of the transverse mode numbers and the initial transverse Gaussian half width. The present research reveals the internal vectorial structure of Hermite Gaussian beam from other viewpoint.展开更多
This paper presents a systematic method to derive Beam Propagation Models for optical waveguides.The technique is based on the use of the symbolic calculus rules for pseudodifferential operators.The cases of straight ...This paper presents a systematic method to derive Beam Propagation Models for optical waveguides.The technique is based on the use of the symbolic calculus rules for pseudodifferential operators.The cases of straight and bent optical waveguides are successively considered.展开更多
In this article a new approach is considered for implementing operator splitting methods for transport problems, influenced by electric fields. Our motivation came to model PE-CVD (plasma-enhanced chemical vapor depos...In this article a new approach is considered for implementing operator splitting methods for transport problems, influenced by electric fields. Our motivation came to model PE-CVD (plasma-enhanced chemical vapor deposition) processes, means the flow of species to a gas-phase, which are influenced by an electric field. Such a field we can model by wave equations. The main contributions are to improve the standard discretization schemes of each part of the coupling equation. So we discuss an improvement with implicit Runge- Kutta methods instead of the Yee’s algorithm. Further we balance the solver method between the Maxwell and Transport equation.展开更多
In this paper, a split step ABCD matrix method is suggested to investigate Gaussian beam propagation in a Kerr type metamaterial medium. This method is based on dividing the medium interval into subsequent steps.Meanw...In this paper, a split step ABCD matrix method is suggested to investigate Gaussian beam propagation in a Kerr type metamaterial medium. This method is based on dividing the medium interval into subsequent steps.Meanwhile, Gaussian beam profile in every step is obtained by finding the ABCD matrix of that particular step, and is used to find the ABCD matrix of the next step. Results of the suggested matrix method have been compared with the results of numerical split-step Fourier method for a Kerr medium, which indicates a good agreement. Then, we use the ABCD matrix to investigate Gaussian beams propagation in a Kerr type metamaterial, which is also in agreement with pervious results by other methods.展开更多
文摘A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.
文摘A new numerical technique based on the wavelet derivative operator is presented as an alternative to BPM to study the integrated optical waveguide. The wavelet derivative operator is used instead of FFT/IFFT or finite difference to calculate the derivatives of the transverse variable in the Helmholtz equation. Results of numerically simulating the injected field at z =0 are exhibited with Gaussian distribution in transverse direction propagating through the two dimensional waveguides (with linear and/or nonlinear refractive index) , which are similar to those in the related publications. Consequently it is efficient and needs not absorbing boundary by introducing the interpolation operator during calculating the wavelet derivative operator. The iterative process needs fewer steps to be stable. Also, when the light wave meets the changes of mediums, the wavelet derivative operator has the adaptive property to adjust those changes at the boundaries.
文摘A modified wide-angle beam propagation based on the Douglas operators is presented.The truncation error in the modified wide-angle beam propagation is reduced to o (Δ x ) 4 in the transverse direction nearly without any increase of the computation time,whereas the error in the ordinary wide-angle beam propagation method is typically o (Δ x ) 2.With trivial programming changes,the accuracy is higher,especially in wide-angle propagation.
文摘A beam propagation method based on the Galerkin method with Hermite-Gauss basis functions for studying optical field propagation in weakly guiding dielectric structures is described. The selected basis functions naturally satisfy the required boundary conditions at infinity so that the boundary truncation is avoided. The paraxial propagation equation is converted into a set of first-order ordinary differential equations, which are solved by means of standard numerical library routines. Besides, the calculation is efficient due to its small resulted matrix. The evolution of the injected field and its normalized power along the propagation distance in an asymmetric slab waveguide and directional coupler are presented, and the solutions are good agreement with those obtained by finite difference BPM, which tests the validity of the present approach.
基金This project is supported by the National Natural Science Foundation of China(NNSF 18572029).
文摘Based on the dynamic governing equation of propagating buckle on a beam on a nonlinear elastic foundation, this paper deals with an important problem of buckle arrest by combining the FEM with a time integration technique. A new conclusion completely different from that by the quasi-static analysis about the buckle arrestor design is drawn. This shows that the inertia of the beam cannot be ignored in the analysis under consideration, especially when the buckle propagation is suddenly stopped by the arrestors.
基金This work was financially supported by the National Natural Science Foundation of China(No.19572029)
文摘For exact estimation of efficiency of a buckle arrestor, it is necessary to take the effect of structural inertia into account in the analysis of buckle propagation on elastic structures after meeting arrestors. Under this consideration, this paper deals with the dynamics of buckle arrest and its numerical simulation on the basis of the beam system model used by Chater and Hutchinson (1983). The FEM combined with an improving are-length control method is adopted to solve the dynamic equations describing the arresting of buckle propagation. A new group of parameters for arrestor design which differs greatly from that by the quasi-static analysis is obtained. The present results support the conclusion that the inertia of the beam cannot be neglected in such analysis.
文摘A modified alternating direction implicit approach is proposed to discretize the three-dimensional full-vectorial beam propagation method （3D-FV-BPM） formulation along the longitudinal direction. The cross-coupling terms （CCTs） are neglected at the first substep, and then double used at the second substep. The order of two substeps is reversed for each transverse electric field component so that the CCTs are always expressed in an implicit form, thus the calculation is efficient and stable. Based on the multinomial interpolation, a universal finite difference scheme with a high accuracy is developed to approximate the 3D-FV-BPM formulation along the transverse directions, in which the discontinuities of the normal components of the electric field across the abrupt dielectric interfaces are taken into account and can be applied to both uniform and non-uniform grids. The corresponding imaginary-distance procedure is first applied to a buried rectangular and a GaAs-based deeply-etched rib waveguide. The field patterns and the normalized propagation constants of the fundamental and the first order modes are presented and the hybrid nature of the full-vectorial guided-modes is demonstrated, which shows the validity and utility of the present approach. Then the modal characteristics of the deeply- and shallow-etched rib waveguides based on the InGaAsp/InGaAsP strained multiple quantum wells in InP substrate are investigated in detail. The results are necessary for modeling and the design of the planar lightwave circuits or photonic integrated circuits based on these waveguides.
文摘The propagation of an optical vortex in a hexagonally arranged single mode multicore fiber structure is investigated for possible generation of additional vortices and their spread dynamics. Fields are separated into a slowly varying paraxial envelope and a rapidly changing exponential component. Solutions are derived from the paraxial inhomogeneous Schrodinger equation in two dimensions along with the index of refraction of the proposed structure. Numerical analyses are based on the beam propagation method and transparent boundary conditions in matrix form with different parameters to represent the intensity and phase of all derived fields. Vortices are numerically identified by their points of zero intensity and their phase change or polarity. The optical interferogram with a plane wave reference is also employed to distinguish the dislocation points in the transverse directions of the propagating fields.
文摘In this paper the analytical solutions of the impact of a particle on Timoshenko beams with four kinds of different boundary conditions are obtained according to Navier's idea, which is further developed. The initial values of the impact forces are exactly determined by the momentum conservation law. The propagation of the longitudinal and transverse waves along the beam, especially, the effects of boundary conditions on the characteristics of the reflected waves, are investigated in detail. Some results are compared with those by MSC/NASTRAN.
文摘A novel three-dimensional beam propagation method (BPM) based on the variable transformed Galerkin's method is introduced for simulating optical field propagation in three-dimensional dielectric structures. The infinite Cartesian x-y plane is mapped into a unit square by a tangent-type function transformation. Consequently, the infinite region problem is converted into the finite region problem. Thus, the boundary truncation is eliminated and the calculation accuracy is promoted. The three-dimensional BPM basic equation is reduced to a set of first-order ordinary differential equations through sinusoidal basis function, which fits arbitrary cladding optical waveguide, then direct solution of the resulting equations by means of the Runge-Kutta method. In addition, the calculation is efficient due to the small matrix derived from the present technique. Both z-invariant and z-variant examples are considered to test both the accuracy and utility of this approach.
文摘Based on the vectorial structure of electromagnetic beam and the method of stationary phase, the analytical structure of Hermite Gaussian beam in far field is presented. The structural energy flux distributions are also investigated in the far field. The structural pictures of some Hermite Gaussian beams are depicted in the far field. As the structure of Hermite Gaussian beam is dominated by the transverse mode numbers and the initial transverse Gaussian half width, it is more complex than that of Gaussian beam. The ratios of the structural energy fluxes to the whole energy flux are independent of the transverse mode numbers and the initial transverse Gaussian half width. The present research reveals the internal vectorial structure of Hermite Gaussian beam from other viewpoint.
文摘This paper presents a systematic method to derive Beam Propagation Models for optical waveguides.The technique is based on the use of the symbolic calculus rules for pseudodifferential operators.The cases of straight and bent optical waveguides are successively considered.
文摘In this article a new approach is considered for implementing operator splitting methods for transport problems, influenced by electric fields. Our motivation came to model PE-CVD (plasma-enhanced chemical vapor deposition) processes, means the flow of species to a gas-phase, which are influenced by an electric field. Such a field we can model by wave equations. The main contributions are to improve the standard discretization schemes of each part of the coupling equation. So we discuss an improvement with implicit Runge- Kutta methods instead of the Yee’s algorithm. Further we balance the solver method between the Maxwell and Transport equation.
文摘In this paper, a split step ABCD matrix method is suggested to investigate Gaussian beam propagation in a Kerr type metamaterial medium. This method is based on dividing the medium interval into subsequent steps.Meanwhile, Gaussian beam profile in every step is obtained by finding the ABCD matrix of that particular step, and is used to find the ABCD matrix of the next step. Results of the suggested matrix method have been compared with the results of numerical split-step Fourier method for a Kerr medium, which indicates a good agreement. Then, we use the ABCD matrix to investigate Gaussian beams propagation in a Kerr type metamaterial, which is also in agreement with pervious results by other methods.
