Advances in Numerical Analysis
https://www.hindawi.com
The latest articles from Hindawi
© 2017 , Hindawi Limited . All rights reserved.

Numerical Simulation of Barite Sag in Pipe and Annular Flow
Thu, 10 Aug 2017 00:00:00 +0000
http://www.hindawi.com/journals/ana/2017/2672438/
With the ever increasing global energy demand and diminishing petroleum reserves, current advances in drilling technology have resulted in numerous directional wells being drilled as operators strive to offset the everrising operating costs. In as much as deviatedwell drilling allows drillers to exploit reservoir potential by penetrating the pay zone in a horizontal, rather than vertical, fashion, it also presents conditions under which the weighting agents can settle out of suspension. The present work is categorized into two parts. In the first part, governing equations were built inside a twodimensional horizontal pipe geometry and the finite element method utilized to solve the equationsets. In the second part, governing equations were built inside a threedimensional horizontal annular geometry and the finite volume method utilized to solve the equationsets. The results of the first part of the simulation are the solid concentration, mixture viscosity, and a prediction of the barite bed characteristics. For the second part, simulation results show that the highest occurrence of barite sag is at low annular velocities, nonrotating drill pipe, and eccentric drill pipe. The CFD approach in this study can be utilized as a research study tool in understanding and managing the barite sag problem.
Patrick Kabanda and Mingbo Wang
Copyright © 2017 Patrick Kabanda and Mingbo Wang. All rights reserved.

Numerical Analysis on Flow Behavior of Molten Iron and Slag in Main Trough of Blast Furnace during Tapping Process
Thu, 09 Mar 2017 00:00:00 +0000
http://www.hindawi.com/journals/ana/2017/6713160/
The threedimensional model was developed according to number 4 of the main trough of blast furnace at China Steel Co. (CSC BF4). The equations and volume of fluid (VOF) were used for describing the turbulent flow at the impinging zone of trough, indicating fluids of liquid iron, molten slag, and air in the governing equation, respectively, in this paper. The pressure field and velocity profile were then obtained by the finite volume method (FVM) and the pressure implicit with splitting of operators (PISO), respectively, followed by calculating the wall shear stress through Newton’s law of viscosity for validation. Then, the operation conditions and the main trough geometry were numerically examined for the separation efficiency of iron from slag stream. As shown in the results, the molten iron losses associated with the slag can be reduced by increasing the height difference between the slag and iron ports, reducing the tapping rate, and increasing the height of the opening under the skimmer.
Li Wang, ChienNan Pan, and WenTung Cheng
Copyright © 2017 Li Wang et al. All rights reserved.

The Convergence of a Class of Parallel NewtonType Iterative Methods
Sun, 05 Mar 2017 07:02:28 +0000
http://www.hindawi.com/journals/ana/2017/4281684/
A general iterative process is proposed, from which a class of parallel Newtontype iterative methods can be derived. A unified convergence theorem for the general iterative process is established. The convergence of these Newtontype iterative methods is obtained from the unified convergence theorem. The results of efficiency analyses and numerical example are satisfactory.
Qinglong Huang
Copyright © 2017 Qinglong Huang. All rights reserved.

Analysis of Tensegrity Structures with Redundancies, by Implementing a Comprehensive Equilibrium Equations Method with Force Densities
Sun, 25 Dec 2016 13:55:36 +0000
http://www.hindawi.com/journals/ana/2016/4815289/
A general approach is presented to analyze tensegrity structures by examining their equilibrium. It belongs to the class of equilibrium equations methods with force densities. The redundancies are treated by employing Castigliano’s second theorem, which gives the additional required equations. The partial derivatives, which appear in the additional equations, are numerically replaced by statically acceptable internal forces which are applied on the structure. For both statically determinate and indeterminate tensegrity structures, the properties of the resulting linear system of equations give an indication about structural stability. This method requires a relatively small number of computations, it is direct (there is no iteration procedure and calculation of auxiliary parameters) and is characterized by its simplicity. It is tested on both 2D and 3D tensegrity structures. Results obtained with the method compare favorably with those obtained by the Dynamic Relaxation Method or the Adaptive Force Density Method.
Miltiades Elliotis, Petros Christou, and Antonis Michael
Copyright © 2016 Miltiades Elliotis et al. All rights reserved.

Revisited Optimal Error Bounds for Interpolatory Integration Rules
Wed, 16 Nov 2016 05:45:03 +0000
http://www.hindawi.com/journals/ana/2016/3170595/
We present a unified way to obtain optimal error bounds for general interpolatory integration rules. The method is based on the Peano form of the error term when we use Taylor’s expansion. These bounds depend on the regularity of the integrand. The method of integration by parts “backwards” to obtain bounds is also discussed. The analysis includes quadrature rules with nodes outside the interval of integration. Best error bounds for composite integration rules are also obtained. Some consequences of symmetry are discussed.
François Dubeau
Copyright © 2016 François Dubeau. All rights reserved.

Analysis of Subgrid Stabilization Method for StokesDarcy Problems
Thu, 22 Sep 2016 14:19:10 +0000
http://www.hindawi.com/journals/ana/2016/7389102/
A number of techniques, used as remedy to the instability of the Galerkin finite element formulation for Stokes like problems, are found in the literature. In this work we consider a coupled StokesDarcy problem, where in one part of the domain the fluid motion is described by Stokes equations and for the other part the fluid is in a porous medium and described by Darcy law and the conservation of mass. Such systems can be discretized by heterogeneous mixed finite elements in the two parts. A better method, from a computational point of view, consists in using a unified approach on both subdomains. Here, the coupled StokesDarcy problem is analyzed using equalorder velocity and pressure approximation combined with subgrid stabilization. We prove that the obtained finite element solution is stable and converges to the classical solution with optimal rates for both velocity and pressure.
Kamel Nafa
Copyright © 2016 Kamel Nafa. All rights reserved.

Analysis of a Decoupled TimeStepping Scheme for Evolutionary Micropolar Fluid Flows
Mon, 19 Sep 2016 13:01:27 +0000
http://www.hindawi.com/journals/ana/2016/7010645/
Micropolar fluid model consists of NavierStokes equations and microrotational velocity equations describing the dynamics of flows in which microstructure of fluid is important. In this paper, we propose and analyze a decoupled timestepping algorithm for the evolutionary micropolar flow. The proposed method requires solving only one uncoupled NavierStokes and one microrotation subphysics problem per time step. We derive optimal order error estimates in suitable norms without assuming any stability condition or time step size restriction.
S. S. Ravindran
Copyright © 2016 S. S. Ravindran. All rights reserved.

Lebesgue Constant Using Sinc Points
Mon, 19 Sep 2016 09:21:25 +0000
http://www.hindawi.com/journals/ana/2016/6758283/
Lebesgue constant for Lagrange approximation at Sinc points will be examined. We introduce a new barycentric form for Lagrange approximation at Sinc points. Using Thiele’s algorithm we show that the Lebesgue constant grows logarithmically as the number of interpolation Sinc points increases. A comparison between the obtained upper bound of Lebesgue constant using Sinc points and other upper bounds for different set of points, like equidistant and Chebyshev points, is introduced.
Maha Youssef, Hany A. ElSharkawy, and Gerd Baumann
Copyright © 2016 Maha Youssef et al. All rights reserved.

A Finite VolumeComplete Flux Scheme for a Polluted Groundwater Site
Wed, 24 Aug 2016 13:44:03 +0000
http://www.hindawi.com/journals/ana/2016/1945958/
We present a model of a polluted groundwater site. The model consists of a coupled system of advectiondiffusionreaction equations for the groundwater level and the concentration of the pollutant. We use the complete flux scheme for the space discretization in combination with the method for time integration and we prove a new stability result for the scheme. Numerical results are computed for the Guarani Aquifer in South America and they show good agreement with results in literature.
M. F. P. ten Eikelder, J. H. M. ten Thije Boonkkamp, M. P. T. Moonen, and B. V. Rathish Kumar
Copyright © 2016 M. F. P. ten Eikelder et al. All rights reserved.

Selection of an Interval for Variable Shape Parameter in Approximation by Radial Basis Functions
Thu, 18 Aug 2016 11:55:30 +0000
http://www.hindawi.com/journals/ana/2016/1397849/
In radial basis function approximation, the shape parameter can be variable. The values of the variable shape parameter strategies are selected from an interval which is usually determined by trial and error. As yet there is not any algorithm for determining an appropriate interval, although there are some recipes for optimal values. In this paper, a novel algorithm for determining an interval is proposed. Different variable shape parameter strategies are examined. The results show that the determined interval significantly improved the accuracy and is suitable enough to count on in variable shape parameter strategies.
Jafar Biazar and Mohammad Hosami
Copyright © 2016 Jafar Biazar and Mohammad Hosami. All rights reserved.

Analysis of Linear Piecewise Constant Delay Systems Using a Hybrid Numerical Scheme
Sun, 07 Aug 2016 13:29:54 +0000
http://www.hindawi.com/journals/ana/2016/8903184/
An efficient computational technique for solving linear delay differential equations with a piecewise constant delay function is presented. The new approach is based on a hybrid of blockpulse functions and Legendre polynomials. A key feature of the proposed framework is the excellent representation of smooth and especially piecewise smooth functions. The operational matrices of delay, derivative, and product corresponding to the mentioned hybrid functions are implemented to transform the original problem into a system of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the proposed numerical scheme.
H. R. Marzban and S. M. Hoseini
Copyright © 2016 H. R. Marzban and S. M. Hoseini. All rights reserved.

Remarks on Numerical Experiments of the AllenCahn Equations with Constraint via Yosida Approximation
Tue, 07 Jun 2016 12:12:45 +0000
http://www.hindawi.com/journals/ana/2016/1492812/
We consider a onedimensional AllenCahn equation with constraint from the viewpoint of numerical analysis. The constraint is provided by the subdifferential of the indicator function on the closed interval, which is the multivalued function. Therefore, it is very difficult to perform a numerical experiment of our equation. In this paper we approximate the constraint by the Yosida approximation. Then, we study the approximating system of the original model numerically. In particular, we give the criteria for the standard forward Euler method to give the stable numerical experiments of the approximating equation. Moreover, we provide the numerical experiments of the approximating equation.
Tomoyuki Suzuki, Keisuke Takasao, and Noriaki Yamazaki
Copyright © 2016 Tomoyuki Suzuki et al. All rights reserved.

EighthOrder Compact Finite Difference Scheme for 1D Heat Conduction Equation
Mon, 16 May 2016 10:05:30 +0000
http://www.hindawi.com/journals/ana/2016/8376061/
The purpose of this paper is to develop a highorder compact finite difference method for solving onedimensional (1D) heat conduction equation with Dirichlet and Neumann boundary conditions, respectively. A parameter is used for the direct implementation of Dirichlet and Neumann boundary conditions. The introduced parameter adjusts the position of the neighboring nodes very next to the boundary. In the case of Dirichlet boundary condition, we developed eighthorder compact finite difference method for the entire domain and fourthorder accurate proposal is presented for the Neumann boundary conditions. In the case of Dirichlet boundary conditions, the introduced parameter behaves like a free parameter and could take any value from its defined domain but for the Neumann boundary condition we obtained a particular value of the parameter. In both proposed compact finite difference methods, the order of accuracy is the same for all nodes. The time discretization is performed by using CrankNicholson finite difference method. The unconditional convergence of the proposed methods is presented. Finally, a set of 1D heat conduction equations is solved to show the validity and accuracy of our proposed methods.
Asma Yosaf, Shafiq Ur Rehman, Fayyaz Ahmad, Malik Zaka Ullah, and Ali Saleh Alshomrani
Copyright © 2016 Asma Yosaf et al. All rights reserved.

Preconditioning and Uniform Convergence for ConvectionDiffusion Problems Discretized on ShishkinType Meshes
Sun, 28 Feb 2016 10:40:31 +0000
http://www.hindawi.com/journals/ana/2016/2161279/
A onedimensional linear convectiondiffusion problem with a perturbation parameter multiplying the highest derivative is considered. The problem is solved numerically by using the standard upwind scheme on special layeradapted meshes. It is proved that the numerical solution is uniform accurate in the maximum norm. This is done by a new proof technique in which the discrete system is preconditioned in order to enable the use of the principle where “uniform stability plus uniform consistency implies uniform convergence.” Without preconditioning, this principle cannot be applied to convectiondiffusion problems because the consistency error is not uniform in . At the same time, the condition number of the discrete system becomes independent of due to the same preconditioner; otherwise, the condition number of the discrete system before preconditioning increases when tends to 0. We obtained such results in an earlier paper, but only for the standard Shishkin mesh. In a nontrivial generalization, we show here that the same proof techniques can be applied to the whole class of Shishkintype meshes.
Thái Anh Nhan and Relja Vulanović
Copyright © 2016 Thái Anh Nhan and Relja Vulanović. All rights reserved.

The Exponential Cubic BSpline Algorithm for Kortewegde Vries Equation
Tue, 17 Feb 2015 13:07:04 +0000
http://www.hindawi.com/journals/ana/2015/367056/
The exponential cubic Bspline algorithm is presented to find the numerical solutions of the Kortewegde Vries (KdV) equation. The problem is reduced to a system of algebraic equations, which is solved by using a variant of Thomas algorithm. Numerical experiments are carried out to demonstrate the efficiency of the suggested algorithm.
Ozlem Ersoy and Idris Dag
Copyright © 2015 Ozlem Ersoy and Idris Dag. All rights reserved.

Approximation Properties of a Gradient Recovery Operator Using a Biorthogonal System
Tue, 03 Feb 2015 06:34:06 +0000
http://www.hindawi.com/journals/ana/2015/187604/
A gradient recovery operator based on projecting the discrete gradient onto the standard finite element space is considered. We use an oblique projection, where the test and trial spaces are different, and the bases of these two spaces form a biorthogonal system. Biorthogonality allows efficient computation of the recovery operator. We analyze the approximation properties of the gradient recovery operator. Numerical results are presented in the twodimensional case.
Bishnu P. Lamichhane and Adam McNeilly
Copyright © 2015 Bishnu P. Lamichhane and Adam McNeilly. All rights reserved.

A FourthOrder Collocation Scheme for TwoPoint Interface Boundary Value Problems
Mon, 24 Nov 2014 12:55:44 +0000
http://www.hindawi.com/journals/ana/2014/875013/
A fourthorder accurate orthogonal spline collocation scheme is formulated to approximate linear twopoint boundary value problems with interface conditions. The coefficients of the differential operator may have jump discontinuities at the interface point, a nodal point of the scheme. Existence and uniqueness of the numerical solution are proved. Optimal order error estimates in the maximum norm are obtained, and a superconvergence property of the numerical solution in the maximal nodal norm is proved. Numerical results are presented confirming the theoretical estimates.
Rakhim Aitbayev and Nazgul Yergaliyeva
Copyright © 2014 Rakhim Aitbayev and Nazgul Yergaliyeva. All rights reserved.

A Method to Accelerate the Convergence of the Secant Algorithm
Wed, 19 Nov 2014 09:08:14 +0000
http://www.hindawi.com/journals/ana/2014/321592/
We present an acceleration technique for the Secant method. The Secant method is a rootsearching algorithm for a general function . We exploit the fact that the combination of two Secant steps leads to an improved, socalled firstorder approximant of the root. The original Secant algorithm can be modified to a firstorder accelerated algorithm which generates a sequence of firstorder approximants. This process can be repeated: two th order approximants can be combined in a th order approximant and the algorithm can be modified to an th order accelerated algorithm which generates a sequence of such approximants. We show that the sequence of th order approximants converges to the root with the same order as methods using polynomial fits of of degree .
M. J. P. Nijmeijer
Copyright © 2014 M. J. P. Nijmeijer. All rights reserved.

Iterative Algorithms for a Finite Family of Multivalued QuasiNonexpansive Mappings
Thu, 13 Nov 2014 07:18:11 +0000
http://www.hindawi.com/journals/ana/2014/181049/
Let be a nonempty closed and convex subset of a uniformly convex real Banach space and let be multivalued quasinonexpansive mappings. A new iterative algorithm is constructed and the corresponding sequence is proved to be an approximating fixed point sequence of each ; that is, . Then, convergence theorems are proved under appropriate additional conditions. Our results extend and improve some important recent results (e.g., Abbas et al. (2011)).
C. Diop, M. Sene, and N. Djitté
Copyright © 2014 C. Diop et al. All rights reserved.

Identifying an Unknown Coefficient in the ReactionDiffusion Equation Using He’s VIM
Sun, 19 Oct 2014 11:20:05 +0000
http://www.hindawi.com/journals/ana/2014/547973/
An inverse heat problem of finding an unknown parameter p(t) in the parabolic initialboundary value problem is solved with variational iteration method (VIM). For solving the discussed inverse problem, at first we transform it into a nonlinear direct problem and then use the proposed method. Also an error analysis is presented for the method and prior and posterior error bounds of the approximate solution are estimated. The main property of the method is in its flexibility and ability to solve nonlinear equation accurately and conveniently. Some examples are given to illustrate the effectiveness and convenience of the method.
F. Parzlivand and A. M. Shahrezaee
Copyright © 2014 F. Parzlivand and A. M. Shahrezaee. All rights reserved.

A ThirdOrder NewtonType Method for Finding Polar Decomposition
Tue, 30 Sep 2014 12:08:45 +0000
http://www.hindawi.com/journals/ana/2014/576325/
It is attempted to present an iteration method for finding polar decomposition. The approach is categorized in the scope of Newtontype methods. Error analysis and rate of convergence are studied. Some illustrations are also given to disclose the numerical behavior of the proposed method.
F. Khaksar Haghani
Copyright © 2014 F. Khaksar Haghani. All rights reserved.

The Perturbation Bound for the Spectral Radius of a Nonnegative Tensor
Sun, 28 Sep 2014 08:06:00 +0000
http://www.hindawi.com/journals/ana/2014/109525/
We study the perturbation bound for the spectral radius of an mthorder ndimensional nonnegative tensor . The main contribution of this paper is to show that when is perturbed to a nonnegative tensor by , the absolute difference between the spectral radii of and is bounded by the largest magnitude of the ratio of the ith component of and the ith component , where is an eigenvector associated with the largest eigenvalue of in magnitude and its entries are positive. We further derive the bound in terms of the entries of only when is not known in advance. Based on the perturbation analysis, we make use of the NQZ algorithm
to estimate the spectral radius of a nonnegative tensor in general. On the other hand, we study the backward error matrix and obtain its smallest error bound for its perturbed largest eigenvalue and associated eigenvector of an irreducible nonnegative tensor. Based on the backward error analysis, we can estimate the stability of computation of the largest eigenvalue of an irreducible nonnegative tensor by the NQZ algorithm. Numerical examples are presented to illustrate the theoretical results of our perturbation analysis.
Wen Li and Michael K. Ng
Copyright © 2014 Wen Li and Michael K. Ng. All rights reserved.

On the Stream FunctionVorticity Finite Element Formulation for Incompressible Flow in Porous Media
Thu, 18 Sep 2014 06:14:19 +0000
http://www.hindawi.com/journals/ana/2014/459761/
Stream functionvorticity finite element formulation for incompressible flow in porous media is presented. The model consists of the heat equation, the equation for the concentration, and the equations of motion under the Darcy law. The existence of solution for the discrete problem is established. Optimal a priori error estimates are given.
Abdellatif Agouzal, Karam Allali, and Siham Binna
Copyright © 2014 Abdellatif Agouzal et al. All rights reserved.

Error Analysis for a Noisy Lacunary Cubic Spline Interpolation and a Simple Noisy Cubic Spline Quasi Interpolation
Wed, 10 Sep 2014 09:10:39 +0000
http://www.hindawi.com/journals/ana/2014/353194/
We mainly present the error analysis for two new cubic spline based methods; one is a lacunary interpolation method and the other is a very simple
quasi interpolation method. The new methods are able to reconstruct a function and its first two derivatives from noisy function data. The explicit error bounds for the methods
are given and proved. Numerical tests and comparisons are performed. Numerical results verify the efficiency of our methods.
FengGong Lang and XiaoPing Xu
Copyright © 2014 FengGong Lang and XiaoPing Xu. All rights reserved.

A Meshless Method for the Numerical Solution of a TwoDimension IHCP
Sun, 03 Aug 2014 06:45:04 +0000
http://www.hindawi.com/journals/ana/2014/456142/
This paper uses the collocation method and radial basis functions (RBFs) to analyze the solution of a twodimension inverse heat conduction problem (IHCP). The accuracy of the method is tested in terms of Error and RMS errors. Also, the stability of the technique is investigated by perturbing the additional specification data by increasing the amounts of random noise. The results of numerical experiments are compared with the analytical solution in illustrative examples to confirm the accuracy and efficiency of the presented scheme.
F. Parzlivand and A. M. Shahrezaee
Copyright © 2014 F. Parzlivand and A. M. Shahrezaee. All rights reserved.

An Efficient Family of TraubSteffensenType Methods for Solving Systems of Nonlinear Equations
Wed, 02 Jul 2014 09:14:36 +0000
http://www.hindawi.com/journals/ana/2014/152187/
Based on TraubSteffensen method, we present a derivative free threestep family of sixthorder methods for solving systems of nonlinear equations. The local convergence order of the family is determined using firstorder divided difference operator for functions of several variables and direct computation by Taylor's expansion. Computational efficiency is discussed, and a comparison between the efficiencies of the proposed techniques with the existing ones is made. Numerical tests are performed to compare the methods of the proposed family with the existing methods and to confirm the theoretical results. It is shown that the new family is especially efficient in solving large systems.
Janak Raj Sharma and Puneet Gupta
Copyright © 2014 Janak Raj Sharma and Puneet Gupta. All rights reserved.

On Nonlinear Inverse Problems of Heat Transfer with Radiation Boundary Conditions: Application to Dehydration of Gypsum Plasterboards Exposed to Fire
Sun, 13 Apr 2014 13:39:37 +0000
http://www.hindawi.com/journals/ana/2014/634712/
The paper investigates boundary optimal controls and parameter estimates to the wellposedness nonlinear model of dehydration of thermic problems. We summarize the general formulations for the boundary control for initialboundary value problem for nonlinear partial differential equations modeling the heat transfer and derive necessary optimality conditions, including the adjoint equation, for the optimal set of parameters minimizing objective functions . Numerical simulations illustrate several numerical optimization methods, examples, and realistic cases, in which several interesting phenomena are observed. A large amount of computational effort is required to solve the coupled state equation and the adjoint equation (which is backwards in time), and the algebraic gradient equation (which implements the coupling between the adjoint and control variables). The state and adjoint equations are solved using the finite element method.
A. Belmiloudi and F. Mahé
Copyright © 2014 A. Belmiloudi and F. Mahé. All rights reserved.

General Family of Third Order Methods for Multiple Roots of Nonlinear Equations and Basin Attractors for Various Methods
Mon, 31 Mar 2014 16:39:23 +0000
http://www.hindawi.com/journals/ana/2014/963878/
A general scheme of third order convergence is described for finding multiple roots of nonlinear equations. The proposed scheme requires one evaluation of , and each per iteration and contains several known onepoint third order methods for finding multiple roots, as particular cases. Numerical examples are included to confirm the theoretical results and demonstrate convergence behavior of the proposed methods. In the end, we provide the basins of attraction for some methods to observe their dynamics in the complex plane.
Rajni Sharma and Ashu Bahl
Copyright © 2014 Rajni Sharma and Ashu Bahl. All rights reserved.

Monotonicity Preserving Rational Quadratic Fractal Interpolation Functions
Thu, 30 Jan 2014 07:25:34 +0000
http://www.hindawi.com/journals/ana/2014/504825/
Fractal interpolation is an advanced technique for analysis and synthesis of scientific and engineering data. We introduce the rational quadratic fractal interpolation functions (FIFs) through a suitable rational quadratic iterated function system (IFS). The novel notion
of shape preserving fractal interpolation without any shape parameter is introduced through the rational fractal interpolation model in the literature for the first time. For a prescribed set of monotonic data, we derive the sufficient conditions by restricting the scaling factors for shape preserving rational quadratic FIFs. A local modification pertaining to any subinterval is possible in this model if the scaling factors are chosen appropriately. We establish the convergence results of a monotonic rational quadratic FIF to the original function in . For given data with derivatives at grids, our approach generates several monotonicity preserving rational quadratic FIFs, whereas this flexibility is not available in the classical approach. Finally, numerical experiments support the importance of the developed rational quadratic IFS scheme through construction of visually pleasing monotonic rational fractal curves including
the classical one.
A. K. B. Chand and N. Vijender
Copyright © 2014 A. K. B. Chand and N. Vijender. All rights reserved.

Nonlinear DoubleLayer Bragg Waveguide: Analytical and Numerical Approaches to Investigate Waveguiding Problem
Wed, 22 Jan 2014 08:08:52 +0000
http://www.hindawi.com/journals/ana/2014/231498/
The paper is concerned with propagation of surface TE waves in a circular nonhomogeneous twolayered dielectric waveguide filled with nonlinear medium. The problem is reduced to the analysis of a
nonlinear integral equation with a kernel in the form of the Green function. The existence of propagating TE waves for chosen nonlinearity (the Kerr law) is proved using the contraction mapping method. Conditions under which k waves can propagate are obtained, and intervals of localization of the corresponding propagation constants are found. For numerical solution of the problem, a method based on solving an auxiliary Cauchy problem (the shooting method) is proposed. In numerical experiment, two types of nonlinearities are considered and compared: the Kerr nonlinearity and nonlinearity with saturation. New propagation regime is found.
Yury G. Smirnov, Eugenii Yu. Smol’kin, and Dmitry V. Valovik
Copyright © 2014 Yury G. Smirnov et al. All rights reserved.