ترغب بنشر مسار تعليمي؟ اضغط هنا

High-order numerical method for the nonlinear Helmholtz equation with material discontinuities in one space dimension

192   0   0.0 ( 0 )
 نشر من قبل Guy Baruch
 تاريخ النشر 2007
  مجال البحث فيزياء
والبحث باللغة English
 تأليف Guy Baruch




اسأل ChatGPT حول البحث

The nonlinear Helmholtz equation (NLH) models the propagation of electromagnetic waves in Kerr media, and describes a range of important phenomena in nonlinear optics and in other areas. In our previous work, we developed a fourth order method for its numerical solution that involved an iterative solver based on freezing the nonlinearity. The method enabled a direct simulation of nonlinear self-focusing in the nonparaxial regime, and a quantitative prediction of backscattering. However, our simulations showed that there is a threshold value for the magnitude of the nonlinearity, above which the iterations diverge. In this study, we numerically solve the one-dimensional NLH using a Newton-type nonlinear solver. Because the Kerr nonlinearity contains absolute values of the field, the NLH has to be recast as a system of two real equations in order to apply Newtons method. Our numerical simulations show that Newtons method converges rapidly and, in contradistinction with the iterations based on freezing the nonlinearity, enables computations for very high levels of nonlinearity. In addition, we introduce a novel compact finite-volume fourth order discretization for the NLH with material discontinuities.The one-dimensional results of the current paper create a foundation for the analysis of multi-dimensional problems in the future.



قيم البحث

اقرأ أيضاً

We present a novel computational methodology for solving the scalar nonlinear Helmholtz equation (NLH) that governs the propagation of laser light in Kerr dielectrics. The methodology addresses two well-known challenges in nonlinear optics: Singular behavior of solutions when the scattering in the medium is assumed predominantly forward (paraxial regime), and the presence of discontinuities in the % linear and nonlinear optical properties of the medium. Specifically, we consider a slab of nonlinear material which may be grated in the direction of propagation and which is immersed in a linear medium as a whole. The key components of the methodology are a semi-compact high-order finite-difference scheme that maintains accuracy across the discontinuities and enables sub-wavelength resolution on large domains at a tolerable cost, a nonlocal two-way artificial boundary condition (ABC) that simultaneously facilitates the reflectionless propagation of the outgoing waves and forward propagation of the given incoming waves, and a nonlinear solver based on Newtons method. The proposed methodology combines and substantially extends the capabilities of our previous techniques built for 1Dand for multi-D. It facilitates a direct numerical study of nonparaxial propagation and goes well beyond the approaches in the literature based on the augmented paraxial models. In particular, it provides the first ever evidence that the singularity of the solution indeed disappears in the scalar NLH model that includes the nonparaxial effects. It also enables simulation of the wavelength-width spatial solitons, as well as of the counter-propagating solitons.
151 - Chusei Kiumi , Kei Saito 2021
Localization is a characteristic phenomenon of space-inhomogeneous quantum walks in one dimension, where particles remain localized at their initial position. Eigenvectors of time evolution operators are deeply related to the amount of trapping. In t his paper, we introduce the analytical method for the eigenvalue problem using a transfer matrix to quantitatively evaluate localization by deriving the time-averaged limit distribution and reveal the condition of strong trapping.
156 - Kun Xu , Quanhua Sun , 2010
Due to the limited cell resolution in the representation of flow variables, a piecewise continuous initial reconstruction with discontinuous jump at a cell interface is usually used in modern computational fluid dynamics methods. Starting from the di scontinuity, a Riemann problem in the Godunov method is solved for the flux evaluation across the cell interface in a finite volume scheme. With the increasing of Mach number in the CFD simulations, the adaptation of the Riemann solver seems introduce intrinsically a mechanism to develop instabilities in strong shock regions. Theoretically, the Riemann solution of the Euler equations are based on the equilibrium assumption, which may not be valid in the non-equilibrium shock layer. In order to clarify the flow physics from a discontinuity, the unsteady flow behavior of one-dimensional contact and shock wave is studied on a time scale of (0~10000) times of the particle collision time. In the study of the non-equilibrium flow behavior from a discontinuity, the collision-less Boltzmann equation is first used for the time scale within one particle collision time, then the direct simulation Monte Carlo (DSMC) method will be adapted to get the further evolution solution. The transition from the free particle transport to the dissipative Navier-Stokes (NS) solutions are obtained as an increasing of time. The exact Riemann solution becomes a limiting solution with infinite number of particle collisions. For the high Mach number flow simulations, the points in the shock transition region, even though the region is enlarged numerically to the mesh size, should be considered as the points inside a highly non-equilibrium shock layer.
158 - Alvaro Salas 2008
In this paper we show some exact solutions for the general fifth order KdV equation. These solutions are obtained by the extended tanh method.
245 - P.Cruz , E.L. Lakshtanov 2008
We have constructed a sequence of solutions of the Helmholtz equation forming an orthogonal sequence on a given surface. Coefficients of these functions depend on an explicit algebraic formulae from the coefficient of the surface. Moreover, for exter ior Helmholtz equation we have constructed an explicit normal derivative of the Dirichlet Green function. In the same way the Dirichlet-to-Neumann operator is constructed. We proved that normalized coefficients are uniformly bounded from zero.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا