Conservative finite-difference scheme for the problem of propagation of a femtosecond pulse in a nonlinear photonic crystal with nonreflecting boundary conditions

Conservative finite-difference schemes are constructed for the problems of self-action of a femtosecond laser pulse and of second-harmonic generation in a one-dimensional nonlinear photonic crystal with nonreflecting boundary conditions. The invariants of the governing equations are found taking into account these conditions. Nonreflecting conditions substantially improve the efficiency of conservative finite-difference schemes used in the modeling of complex nonlinear effects in photonic crystals, which require much smaller steps in space and time than those used in the case of linear propagation. The numerical experiments performed show that the boundary reflects no more than 0.01% of the transmitted energy, which corresponds to the truncation error in the boundary conditions. The amplitude of the reflected pulse is less than that of the pulse transmitted through the boundary by two (and more) orders of magnitude. The simulation is based on the approach proposed by the authors for the given class of problems.     

Conservative finite-difference scheme for the problem of a femtosecond laser pulse with an axially symmetric profile propagating in a medium with cubic nonlinearity

Conservative finite-difference schemes are constructed for the problem of a femtosecond laser pulse propagating in a cubically nonlinear medium in the axially symmetric case with allowance for temporal dispersion of the nonlinear response of the medium. The process is governed by the nonlinear Schrödinger equation involving the time derivative of the nonlinear term. The invariants of the differential problem are presented. It is shown that the difference analogues of these invariants hold for the solution to the finite-difference schemes proposed for the problem. As an example, the numerical results obtained for the self-focusing of a femtosecond light beam are presented.     

Two-dimensional high-resolution schemes and their application in the modeling of ionizing waves in gas discharges

The method for constructing upwind high-resolution schemes is proposed in application to the modeling of ionizing waves in gas discharges. The flux-limiting criterion for continuity equations is derived using the proposed partial monotony property of a finite difference scheme. For two-dimensional extension, the cone transport upwind approach for constructing genuinely two-dimensional difference schemes is used. It is shown that when calculating rotations of symmetric profiles by using this scheme, a circular form of isolines is not distorted in a distinct from the coordinate splitting method. The conservative second order finite-difference scheme is proposed for solving the equations system of the drift-diffusion model of electric discharge; this scheme implies finite-difference conservation laws of electric charge and full electric current (fully conservative scheme). Computations demonstrate absence of numeric oscillations and good resolution of two-dimensional ionizing fronts in simulations of streamer and barrier discharges     

The splitting finite-difference time-domain methods for Maxwell's equations in two dimensions

In this paper, we consider splitting methods for Maxwell's equations in two dimensions. A new kind of splitting finite-difference time-domain methods on a staggered grid is developed. The corresponding schemes consist of only two stages for each time step, which are very simple in computation. The rigorous analysis of the schemes is given. By the energy method, it is proved that the scheme is unconditionally stable and convergent for the problems with perfectly conducting boundary conditions. Numerical dispersion analysis and numerical experiments are presented to show the efficient performance of the proposed methods. Furthermore, the methods are also applied to solve a scattering problem successfully.     

Two difference schemes for the numerical solution of Maxwell’s equations as applied to extremely and super low frequency signal propagation in the Earth-ionosphere waveguide

Two explicit two-time-level difference schemes for the numerical solution of Maxwell’s equations are proposed to simulate propagation of small-amplitude extremely and super low frequency electromagnetic signals (200 Hz and lower) in the Earth-ionosphere waveguide with allowance for the tensor conductivity of the ionosphere. Both schemes rely on a new approach to time approximation, specifically, on Maxwell’s equations represented in integral form with respect to time. The spatial derivatives in both schemes are approximated to fourth-order accuracy. The first scheme uses field equations and is second-order accurate in time. The second scheme uses potential equations and is fourth-order accurate in time. Comparative test computations show that the schemes have a number of important advantages over those based on finite-difference approximations of time derivatives. Additionally, the potential scheme is shown to possess better properties than the field scheme.     

Consistent two-sided estimates for the solutions of quasilinear parabolic equations and their approximations

For a linearized finite-difference scheme approximating the Dirichlet problem for a multidimensional quasilinear parabolic equation with unbounded nonlinearity, we establish pointwise two-sided solution estimates consistent with similar estimates for the differential problem. These estimates are used to prove the convergence of finite-difference schemes in the grid L ₂ norm.     

A linearized finite-difference scheme for the numerical solution of the nonlinear cubic Schrödinger equation

A linearized finite-difference scheme is used to transform the initial/boundary-value problem associated with the nonlinear Schrödinger equation into a linear algebraic system. This method is developed by re placing the time and the space partial derivatives by parametric finite-difference re placements and the nonlinear term by an appropriate parametric linearized scheme based on Taylor’s expansion. The resulting finite-difference method is analysed for stability and convergence. The results of a number of numerical experiments for the single-soliton wave are given.     

Simulation of soliton solutions for the propagation of a femtosecond pulse in a medium with a cubic nonlinearity

Soliton solutions are constructed numerically for the problem of propagation of a femtosecond pulse in a medium with a cubic nonlinearity. The problem is posed as an eigenvalue problem with an operator nonlinear in the eigenfunctions. For given values of the propagation parameter we find the real eigenvalue λ and the corresponding eigenvector. This eigenvector is a soliton, i.e., a solution that does not vary in the coordinate of propagation of the light pulse. An algorithm is proposed to find the minimum eigenvalue and the corresponding eigenfunctions that satisfy given conditions. Translated from Prikladnaya Matematika i Informatika, No. 2, pp. 63–68, 1999.     

Propagation of TM waves in a layer with arbitrary nonlinearity

A boundary value problem for Maxwell’s equations describing propagation of TM waves in a nonlinear dielectric layer with arbitrary nonlinearity is considered. The layer is located between two linear semi-infinite media. The problem is reduced to a nonlinear boundary eigenvalue problem for a system of second-order nonlinear ordinary differential equations. A dispersion equation for the eigenvalues of the problem (propagation constants) is derived. For a given nonlinearity function, the dispersion equation can be studied both analytically and numerically. A sufficient condition for the existence of at least one eigenvalue is formulated.     

A new averaging scheme for the riemann problem in pure water

The numerical investigation of shock phenomena in gas or liquid media where enthalpy is the preferred thermodynamic variable poses special problems. When an expression for internal energy is available, the usual procedure is to employ a splitting scheme to remove source terms from the Euler equations, then upwind-biased shock capturing algorithms are built around the Riemann problem for the conservative system which remains. However, when the governing equations arc formulated in terms of total enthalpy, treatment of a pressure time derivative as a source term leads to a Riemann problem for a system where one equation is not a conservation law. The present research establishes that successful upwind-biased shock capturing schemes can be based upon the pseudo-conservative system. A new averaging scheme for solving the associated Riemann problem is developed. The method is applied to numerical simulations of shock wave propagation in pure water.     

Construction of a novel finite-difference scheme for a nonlinear diffusion equation

A novel finite-difference scheme is constructed for a particular nonlinear diffusion equation. This scheme does not correspond to any of those that might be expected to follow from application of the standard procedures. A major advantage is that in the appropriate limit the difference equation model reduces to schemes that reproduce exactly the behaviors of the solution to the corresponding differential equation.     

二维非线性对流扩散方程的非振荡特征差分方法

１．引言 近十几年来,双曲守恒律问题的高分辨率格式已取得很大发展,具有局部自适应选取节点的非振荡插值算法（如 ＵＮＯ［１］, ＥＮＯ［２］等）在这些格式的构造中起着重要的作用．特征差分法是求解对流扩散问题的一种较为有效方法,但在求解具有陡峭前线问题时,也会产生非物理振荡阻（见４）．本文将把特征差分法与非振荡插值算法相结合构造对流扩散问题的高分辨率差分格式． ［１］中的 ＵＮＯ及［２］中的 ＥＮＯ插值都是一维的,有关讨论二维 ＵＮＯ及ＥＮＯ插值的文章还不多见,本文将构造二维基于六节点的二次非振荡插值以及…     

Partially Coherent Waves in Nonlinear Periodic Lattices

We study the propagation of partially coherent (random-phase) waves in nonlinear periodic lattices. The dynamics in these systems is governed by the threefold interplay between the nonlinearity, the lattice properties, and the statistical (coherence) properties of the waves. Such dynamic interplay is reflected in the characteristic properties of nonlinear wave phenomena (e.g., solitons) in these systems. While the propagation of partially coherent waves in nonlinear periodic systems is a universal problem, we analyze it in the context of nonlinear photonic lattices, where recent experiments have proven their existence.     

CABARET scheme for the numerical solution of aeroacoustics problems: Generalization to linearized one-dimensional Euler equations

A generalization of the CABARET finite difference scheme is proposed for linearized one-dimensional Euler equations based on the characteristic decomposition into local Riemann invariants. The new method is compared with several central finite difference schemes that are widely used in computational aeroacoustics. Numerical results for the propagation of an acoustic wave in a homogeneous field and the refraction of this wave through a contact discontinuity obtained on a strongly nonuniform grid are presented.     

三维对流扩散方程的三种高精度分裂格式

在算子分裂法思想的基础上,将两种高精度的离散格式推广应用于三维对流扩散方程,同时对经典ADI格式的对流项做了改进,改进后的格式的对流项对空间具有4阶精度,而经典ADI格式对空间只有2阶精度,由此可见,提高了该格式的实用性．最后对两种典型的浓度场进行了数值模拟,将3种格式的计算结果与解析解以及其它传统差分格式的计算结果进行了对比,得出当Peclet数不大于5时,3种格式均获得了令人满意的数值结果,说明推广的这三种方法具有很高的准确性和可靠性．     

On exact dimensional splitting for a multidimensional scalar quasilinear hyperbolic conservation law

A dimensional splitting scheme is applied to a multidimensional scalar homogeneous quasilinear hyperbolic equation (conservation law). It is proved that the splitting error is zero. The proof is presented for the above partial differential equation in an arbitrary number of dimensions. A numerical example is given that illustrates the proved accuracy of the splitting scheme. In the example, the grid convergence of split (locally one-dimensional) compact and bicompact difference schemes and unsplit bicompact schemes combined with high-order accurate time-stepping schemes (namely, Runge–Kutta methods of order 3, 4, and 5) is analyzed. The errors of the numerical solutions produced by these schemes are compared. It is shown that the orders of convergence of the split schemes remain high, which agrees with the conclusion that the splitting error is zero.     

Some numerical experiments on the splitting of Burgers' equation

The combined approach of linearization and splitting up is used for devising new algorithms to solve a one-dimensional Burgers' equation. Two schemes are discussed and the computed solutions are compared with the exact solution. For this problem it is found that the schemes proposed yield excellent numerical results for Reynolds number R ranges from 50 up to 1500. The schemes were also tested for another problem whose R = 10000. In this case a filtering technique is used to overcome the nonlinear instability.     

The Convergence and Stability of Splitting Finite-Difference Schemes for Nonlinear Evolutionary Equations

We consider a splitting finite-difference scheme for an initial-boundary value problem for a two-dimensional nonlinear evolutionary equation. The problem is split into nonlinear and linear parts. The linear part is also split into locally one-dimensional equations. We prove the convergence and stability of the scheme in L ₂ and C norms. Printed in Lietuvos Matematikos Rinkinys, Vol. 45, No. 3, pp. 413–434, July–September, 2005.     

Efficient energy-stable schemes for the hydrodynamics coupled phase-field model

In this article, several efficient and energy-stable semi–implicit schemes are presented for the Cahn–Hilliard phase-field model of two-phase incompressible flows. A scalar auxiliary variable (SAV) approach is implemented to solve the Cahn–Hilliard equation, while a splitting method based on pressure stabilization is used to solve the Navier–Stokes equation. At each time step, the schemes involve solving only a sequence of linear elliptic equations, and computations of the phase-field variable, velocity, and pressure are totally decoupled. A finite-difference method on staggered grids is adopted to spatially discretize the proposed time-marching schemes. We rigorously prove the unconditional energy stability for the semi-implicit schemes and the fully discrete scheme. Numerical results in both two and three dimensions are obtained, which demonstrate the accuracy and effectiveness of the proposed schemes. Using our numerical schemes, we compare the SAV, invariant energy quadratization (IEQ), and stabilization approaches. Bubble rising dynamics and coarsening dynamics are also investigated in detail. The results demonstrate that the SAV approach is more accurate than the IEQ approach and that the stabilization approach is the least accurate among the three approaches. The energy stability of the SAV approach appears to be better than that of the other approaches at large time steps.