Existence and stability of TE modes in a stratified non-linear dielectric

Starting from Maxwell's equations for a stratified optical mediumwith a non-linear refractive index, we derive the equationsfor monochromatic planar TE modes. It is then shown that TEmodes in which the electromagnetic fields are travelling wavescorrespond to solutions of these reduced equations in the formof standing waves. The equations of the paraxial approximationare then formulated and the stability of the travelling wavesis investigated in that context.     

On the problem of propagation of nonlinear coupled TE-TM waves in a layer

The problem of simultaneous propagation of two types of electromagnetic waves (TE and TM) in a plane dielectric waveguide filled with a nonlinear medium is considered. These polarized waves have different frequencies and different propagation constants. The physical problem is reduced to a nonlinear two-parameter transmission eigenvalue problem for Maxwell’s equations in a layer. The coupled eigenvalues are coupled propagation constants. A theorem on the existence and localization of coupled eigenvalues corresponding to coupled polarized electromagnetic waves is proved.     

Optical Soliton Bullets in (2+1)D Nonlinear Bragg Resonant Periodic Geometries

We consider light propagation in a Kerr-nonlinear 2D waveguide with a Bragg grating in the propagation direction and homogeneous in the transverse direction. Using Newton's iteration method we construct both stationary and travelling solitary wave solutions of the corresponding mathematical model, the 2D nonlinear coupled mode equations (2D CME). We call these solutions 2D gap solitons due to their similarity with the gap solitons of 1D CME (fiber grating). Long-time stable evolution preserving the solitary fashion is demonstrated numerically despite the fact that, as we show, for the 2D CME no local constrained minima of the Hamiltonian functional exist. Building on the 1D study of [ 1 ], we demonstrate trapping of slow enough 2D gap solitons at localized defects. We explain the mechanism of trapping as resonant transfer of energy from the soliton to one or more nonlinear defect modes. For a special class of defects, we construct a family of nonlinear defect modes by numerically following a bifurcation curve starting at analytically or numerically known linear defect modes. Compared to 1D the dynamics of trapping are harder to fully analyze and the existence of many defect modes for a given defect potential causes that slow solitons store a part of their energy for virtually all of the studied attractive defects.     

Computer simulation of a twisted nanotube buckling

We develop procedures for solving the problems of dynamic nanostructure deformation and buckling numerically. The procedures are based on discretization with respect to time of the nonlinear equations of molecular mechanics whose matrices and vectors are determined using the Morse potential for the central forces of interaction between atoms and fictitious truss elements accounting for the variations of the angle between atomic bonds. To determine the critical values of deformation parameters and the shapes of buckling nanostructures we use a stability loss criterion for solutions to nonlinear ordinary differential equations on a finite time interval. We implemented our procedures in the PIONER code, using which we solve the problem of a twisted nanotube buckling in the conditions of a quasistatic deformation. To determine the postcritical equilibrium modes we solve the same problem in a dynamic formulation. We show that the modes of equilibrium configurations of the nanotube in the initial postcritical deformation correspond to a buckling mode obtained both at the bifurcation point of quasistatic solutions and at the quasibifurcation point of dynamic solutions.     

A numerical study of electromagnetic waves in periodic waveguides

Here are considered time‐harmonic electromagnetic waves in a quadratic waveguide consisting of a periodic dielectric core enclosed by conducting walls. The permittivity function may be smooth or have jumps. The electromagnetic field is given by a magnetic vector potential in Lorenz gauge, and defined on a Floquet cell. The Helmholtz operator is approximated by a Chebyshev collocation, Fourier–Galerkin method. Laurent's rule and the inverse rule are employed for the representation of Fourier coefficients of products of functions. The computations yield, for known wavenumbers, values of the first few eigenfrequencies of the field. In general, the dispersion curves exhibit band gaps. Field patterns are identified as transverse electric, TE, transverse magnetic, TM, or hybrid modes. Maxwell's equations are fulfilled. A few trivial solutions appear when the permittivity varies in the guiding direction and across it. The results of the present method are consistent with exact results and with those obtained by a low‐order finite element software. The present method is more efficient than the low‐order finite element method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 490–513, 2014     

Analysis of quadruple corner-cut ridged square waveguide using a scaled boundary finite element method

This paper describes the development of an efficient semi-analytical method, namely scaled boundary finite-element method (SBFEM) for a quadruple corner-cut ridged square waveguide. Thinking about its symmetry, only a quarter of its cross-section needs to be considered and divided into a few sub-domains. Only the boundaries of the sub-domains are discretized with line elements leading to great flexibility in mesh generation. The singularities in the re-entrant corners are represented analytically by locating the scaling center in those points. Variational principle approach is presented to formulate the basis SBFE equations for the sub-domains. Then, an equation of the ‘stiffness matrix’ on the discretized boundary is established. Finally, by using the continued-fraction solution and introducing auxiliary variables, a generalized eigenvalue equation with respect to the cutoff wave number is obtained without introducing an internal mesh. Numerical results are presented to verify the accuracy and efficiency of the present technique. Variations of the cutoff wave numbers of the dominant and higher-order modes for both TE and TM cases with the corner-cut ridge dimensions are investigated in details. Simple approximate equations are found to accurately predict the cutoff wave number of TE_20U, TE₂₂, TM₁₁ and TM_13L modes. The single mode bandwidth of the waveguide is also calculated.     

Analysis of microwave induced natural convection in a single mode cavity (Influence of sample volume, placement, and microwave power level)

The heating of water layer using microwave oven with a rectangular waveguide has been studied both numerically and experimentally. The mathematical model is validated with the experimental data. The transient Maxwell’s equations are solved by using the Finite Difference Time Domain (FDTD) method to describe the electromagnetic field inside the waveguide and sample. The temperature profile and velocity field within sample are determined by the solutions of the momentum, energy and Maxwell’s equations. In this study, the effects of physical parameters, e.g. microwave power level, placement of sample inside the waveguide, volume of sample, are studied. The distribution of electric field, temperature profile and velocity field are presented in details. The results show good agreement between simulation results and experimental data. Conclusively, the mathematical model presented here correctly explains the phenomena of microwave heating of water layer.     

Permanent Wave Structures and Resonant Triads in a Layered Fluid

A closed three layer fluid with small density differences between the layers has two closely related modes of gravity wave propagation. The nonlinear interactions between the wave modes are investigated, particularly the nearly resonant or significant interactions. Permanent wave solutions are calculated, and it is shown that a permanent wave of the slower mode can generate resonantly a wave harmonic of the faster mode. The equations governing resonant triads of the two modes are derived, and solutions having a permanent structure are calculated from them. It is found that some resonant triad solutions vanish when the triad is embedded in the set of all harmonics with wavenumbers in its neighborhood     

Nonlinear transmission eigenvalue problem describing TE wave propagation in two-layered cylindrical dielectric waveguides

The paper is concerned with propagation of surface TE waves in a circular nonhomogeneous two-layered dielectric waveguide filled with a Kerr nonlinear medium. The problem is reduced to the analysis of a nonlinear integral equation with a kernel in the form of a Green’s function. The existence of propagating TE waves is proved using the contraction mapping method. For the numerical solution of the problem, two methods are proposed: an iterative algorithm (whose convergence is proved) and a method based on solving an auxiliary Cauchy problem (the shooting method). The existence of roots of the dispersion equation (propagation constants of the waveguide) is proved. Conditions under which k waves can propagate are obtained, and regions of localization of the corresponding propagation constants are found.     

Modelling Axi-symmetric Travelling Waves in a Dielectric with Nonlinear Refractive Index

We consider an isotropic dielectric with a nonlinear refractive index. The medium may be inhomogeneous but its spatial variation has an axial symmetry. We characterize all monochromatic axi-symmetric travelling waves as solution of a system of six second order differential equations on (0, ). Boundary conditions at 0 ensure the regularity of the fields on the axis. Guided waves satisfy additional conditions at . Special solutions of this system correspond to what are normally referred to as TE and TM modes.Lecture held in the Seminario Matematico e Fisico on June 13, 2003Received: February, 2004     

On the evolution equation for magnetic geodesics

Orbits of charged particles under the effect of a magnetic field are mathematically described by magnetic geodesics. They appear as solutions to a system of (nonlinear) ordinary differential equations of second order. But we are only interested in periodic solutions. To this end, we study the corresponding system of (nonlinear) parabolic equations for closed magnetic geodesics and, as a main result, eventually prove the existence of long time solutions. As generalization one can consider a system of elliptic nonlinear partial differential equations (PDEs) whose solutions describe the orbits of closed p-branes under the effect of a “generalized physical force”. For the corresponding evolution equation, which is a system of parabolic nonlinear PDEs associated to the elliptic PDE, we can establish existence of short time solutions.     

Influence of a wall impedance insertion on the behavior of the acoustic field in a waveguide

A rigid waveguide with an impedance covering occupying a bounded portion of its surface is considered. The influence exerted by the impedance insertion on the behavior of the waveguide field downstream of the insertion is studied depending on the impedance value, type, and the relative size of the impedance insertion. Special waveguide field propagation modes are found, and the areas of local maxima in the field suppression level are examined in the case of a layer, circular, and rectangular waveguide.     

Resonances, radiation damping and instabilitym in Hamiltonian nonlinear wave equations

We consider a class of nonlinear Klein-Gordon equations which are Hamiltonian and are perturbations of linear dispersive equations. The unperturbed dynamical system has a bound state, a spatially localized and time periodic solution. We show that, for generic nonlinear Hamiltonian perturbations, all small amplitude solutions decay to zero as time tends to infinity at an anomalously slow rate. In particular, spatially localized and time-periodic solutions of the linear problem are destroyed by generic nonlinear Hamiltonian perturbations via slow radiation of energy to infinity. These solutions can therefore be thought of as metastable states. The main mechanism is a nonlinear resonant interaction of bound states (eigenfunctions) and radiation (continuous spectral modes), leading to energy transfer from the discrete to continuum modes. This is in contrast to the KAM theory in which appropriate nonresonance conditions imply the persistence of invariant tori. A hypothesis ensuring that such a resonance takes place is a nonlinear analogue of the Fermi golden rule, arising in the theory of resonances in quantum mechanics. The techniques used involve: (i) a time-dependent method developed by the authors for the treatment of the quantum resonance problem and perturbations of embedded eigenvalues, (ii) a generalization of the Hamiltonian normal form appropriate for infinite dimensional dispersive systems and (iii) ideas from scattering theory. The arguments are quite general and we expect them to apply to a large class of systems which can be viewed as the interaction of finite dimensional and infinite dimensional dispersive dynamical systems, or as a system of particles coupled to a field. Oblatum: 6-XI-1998 & 12-VI-1998 / Published online: 14 January 1999     

非线性波方程准确孤立波解的符号计算

该文将机械化数学方法应用于偏微分方程领域，建立了构造一类非线性发展方程孤立波解的一种统一算法，并在计算机数学系统上加以实现，推导出了一批非线性发展方程的精确孤立波解．算法的基本原理是利用非线性发展方程孤立波解的局部性特点，将孤立波表示为双曲正切函数的多项式．从而将非线性发展方程（组）的求解问题转化为非线性代数方程组的求解问题．利用吴文俊消元法在计算机代数系统上求解非线性代数方程组，最终获得非线性发展方程（组）的准确孤立波解.     

Exact solutions for a class of nonlinear evolution equations: A unified ansätze approach

In this paper, we propose a modified generalized transformation for constructing analytic solutions to nonlinear differential equations. This improved unified ansätze is utilized to acquire exact solutions that are general solutions of simpler equations that are either integrable or possess special solutions. The ansätze is constructed via the choice of an integrable differential operator or a basis set of functions. The technique is implemented to obtain several families of exact solutions for a class of nonlinear evolution equations with nonlinear term of any order. In particular, the Klein–Gordon, the Sine–Gordon and Landau–Ginburg–Higgs equations are chosen as examples to illustrate the method.     

Strong Interactions between Solitary Waves Belonging to Different Wave Modes

Strong interactions between weakly nonlinear long waves are studied. Strong interactions occur when the linear long wave phase speeds are nearly equal although the waves belong to different modes. Specifically we study this situation in the context of internal wave modes propagating in a density stratified fluid. The interaction is described by two coupled Korteweg-deVries equations, which possess both dispersive and nonlinear coupling terms. It is shown that the coupled equations possess an exact analytical solution involving the characteristic “sech²” profile of the Korteweg-deVries equation. It is also shown that when the coefficients satisfy some special conditions, the coupled equations possess an n-solition solution analogous to the Korteweg-deVries n-solition solution. In general though the coupled equations are found not to be amenable to solution by the inverse scattering transform technique, and thus a numerical method has been employed in order to find solutions. This method is described in detail in Appendix A. Several numerical solutions of the coupled equations are presented.     

Asymptotic Behavior of Solutions of a Strongly Nonlinear Model of a Crystal Lattice

We consider a system of hyperbolic nonlinear equations describing the dynamics of interaction between optical and acoustic modes of a complex crystal lattice (without a symmetry center) consisting of two sublattices. This system can be considered a nonlinear generalization of the well-known Born-Huang Kun model to the case of arbitrarily large sublattice displacements. For a suitable choice of parameters, the system reduces to the sine-Gordon equation or to the classical equations of elasticity theory. If we introduce physically natural dissipative forces into the system, then we can prove that a compact attractor exists and that trajectories converge to equilibrium solutions. In the one-dimensional case, we describe the structure of equilibrium solutions completely and obtain asymptotic solutions for the wave propagation. In the presence of inhomogeneous perturbations, this system is reducible to the well-known Hopfield model describing the attractor neural network and having complex behavior regimes.__________Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 143, No. 3, pp. 357–367, June, 2005.     

Nonlinear Hyperbolic Equations in Infinite Homogeneous Waveguides

ABSTRACT

In this paper we prove global and almost global existence theorems for nonlinear wave equations with quadratic nonlinearities in infinite homogeneous waveguides. We can handle both the case of Dirichlet boundary conditions and Neumann boundary conditions. In the case of Neumann boundary conditions we need to assume a natural nonlinear Neumann condition on the quasilinear terms. The results that we obtain are sharp in terms of the assumptions on the dimensions for the global existence results and in terms of the lifespan for the almost global results. For nonlinear wave equations, in the case where the infinite part of the waveguide has spatial dimension three, the hypotheses in the theorem concern whether or not the Laplacian for the compact base of the waveguide has a zero mode or not.     

Traveling wave solutions of the n-dimensional coupled Yukawa equations

We discuss traveling wave solutions to the Yukawa equations, a system of nonlinear partial differential equations which has applications to meson–nucleon interactions. The Yukawa equations are converted to a six-dimensional dynamical system, which is then studied for various values of the wave speed and mass parameter. The stability of the solutions is discussed, and the methods of competitive modes is used to describe parameter regimes for which chaotic behaviors may appear. Numerical solutions are employed to better demonstrate the dependence of traveling wave solutions on the physical parameters in the Yukawa model. We find a variety of interesting behaviors in the system, a few of which we demonstrate graphically, which depend upon the relative strength of the mass parameter to the wave speed as well as the initial data.     

A method of solving a nonlinear boundary value problem with a parameter for a loaded differential equation

A nonlinear loaded differential equation with a parameter on a finite interval is studied. The interval is partitioned by the load points, at which the values of the solution to the equation are set as additional parameters. A nonlinear boundary value problem for the considered equation is reduced to a nonlinear multipoint boundary value problem for the system of nonlinear ordinary differential equations with parameters. For fixed parameters, we obtain the Cauchy problems for ordinary differential equations on the subintervals. Substituting the values of the solutions to these problems into the boundary condition and continuity conditions at the partition points, we compose a system of nonlinear algebraic equations in parameters. A method of solving the boundary value problem with a parameter is proposed. The method is based on finding the solution to the system of nonlinear algebraic equations composed.