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.展开更多
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.展开更多
文摘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.
文摘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.
