KMITL

Permanent URI for this communityhttps://dspace.kmitl.ac.th/handle/123456789/1

Browse

Settings

Author: search.filters.author.Sakkaplangkul, PutthaSubject: search.filters.subject.Finite differences

Search Results

Now showing 1 - 2 of 2
  • Loading...
    Thumbnail Image
    Some of the metrics are blocked by your 
    consent settings
    Item type:Publication,
    Convergence analysis of operator splitting methods for Maxwell’s equations in dispersive media of Debye type
    (2023-12-01)
    Sakkaplangkul, Puttha
    In this paper, two new effective operator splitting methods (SS-MD and SM-MD) for the Maxwell’s equations for dispersive media in two dimensions transverse electric polarization (the 2D Maxwell–Debye TE model) are presented and analyzed. The splitting schemes consist of two sub-stages in each time step, each of which requires solving a number of 1D discrete sub-problems. The Crank–Nicolson approach is used to solve each sub-problem’s time discretization. Both splitting methods satisfy the energy decay and are unconditionally stable. The convergence result of the SS-MD scheme is shown to be of first order in time and of second order in space based on the energy technique, whereas the SM-MD scheme is of second order in both time and space. We also analyze numerical dispersion analysis to obtain two identities of the discrete numerical dispersion relations of both splitting schemes. Examples and numerical experiments are provided to demonstrate and support our theoretical results.
  • Loading...
    Thumbnail Image
    Some of the metrics are blocked by your 
    consent settings
    Item type:Publication,
    Dispersion analysis of finite difference and discontinuous Galerkin schemes for Maxwell's equations in linear Lorentz media
    (2019-10-01)
    Jiang, Yan
    ;
    Sakkaplangkul, Puttha
    ;
    Bokil, Vrushali A.
    ;
    Cheng, Yingda
    ;
    Li, Fengyan
    In this paper, we consider Maxwell's equations in linear dispersive media described by a single-pole Lorentz model for electronic polarization. We study two classes of commonly used spatial discretizations: finite difference methods (FD) with arbitrary even order accuracy in space and high spatial order discontinuous Galerkin (DG) finite element methods. Both types of spatial discretizations are coupled with second order semi-implicit leap-frog and implicit trapezoidal temporal schemes. By performing detailed dispersion analysis for the semi-discrete and fully discrete schemes, we obtain rigorous quantification of the dispersion error for Lorentz dispersive dielectrics. In particular, comparisons of dispersion error can be made taking into account the model parameters, and mesh sizes in the design of the two types of schemes. This work is a continuation of our previous research on energy-stable numerical schemes for nonlinear dispersive optical media [6,7]. The results for the numerical dispersion analysis of the reduced linear model, considered in the present paper, can guide us in the optimal choice of discretization parameters for the more complicated and nonlinear models. The numerical dispersion analysis of the fully discrete FD and DG schemes, for the dispersive Maxwell model considered in this paper, clearly indicate the dependence of the numerical dispersion errors on spatial and temporal discretizations, their order of accuracy, mesh discretization parameters and model parameters. The results obtained here cannot be arrived at by considering discretizations of Maxwell's equations in free space. In particular, our results contrast the advantages and disadvantages of using high order FD or DG schemes and leap-frog or trapezoidal time integrators over different frequency ranges using a variety of measures of numerical dispersion errors. Finally, we highlight the limitations of the second order accurate temporal discretizations considered.