Search

In this paper, we present two multiple scalar auxiliary variable (MSAV)-based, finite element numerical schemes for the Abels-Garcke-Gr{\"u}n (AGG) model, which is a thermodynamically consistent phase field model of two-phase incompressible flows with different densities. Both schemes are decoupled, linear, second-order in time, and the numerical implementation turns out to be straightforward. The first scheme solves the Navier-Stokes equations in a saddle point formulation, while the second one employs the artificial compressibility method, leading to a fully decoupled structure with a time-independent pressure update equation. In terms of computational cost, only a sequence of independent elliptic or saddle point systems needs to be solved at each time step. At a theoretical level, the unique solvability and unconditional energy stability (with respect to a modified energy functional) of the proposed schemes are established. In addition, comprehensive numerical simulations are performed to verify the effectiveness and robustness of the proposed schemes.

This paper studies a non-singular coupling scheme for solving the acoustic and elastic wave scattering problems and its extension to the problems of Laplace and Lam\'e equations and the problem with a compactly supported inhomogeneity is also briefly discussed. Relying on the solution representation of the wave scattering problem, a Robin-type artificial boundary condition in terms of layer potentials whose kernels are non-singular, is introduced to obtain a reduced problem on a bounded domain. The wellposedness of the reduced problems and the a priori error estimates of the corresponding finite element discretization are proved. Numerical examples are presented to demonstrate the accuracy and efficiency of the proposed method.

In this work, we investigate the propagation of electromagnetic waves in the Cole-Cole dispersive medium by using the discontinuous Galerkin (DG) method to solve the coupled time-domain Maxwell's equations and polarization equation. We define a new and sharpened total energy function for the Cole-Cole model, which better describes the behaviors of the energy than what is available in the current literature. A major theme in the time-domain numerical modeling of this problem has been tackling the difficulty of handling the nonlocal term involved in the time-domain polarization equation. Based on the diffusive representation and the quadrature formula, we derive an approximate system, where the convolution kernel is replaced by a finite number of auxiliary variables that satisfy local-in-time ordinary differential equations. To ensure the resulted approximate system is stable, a nonlinear constrained optimization numerical scheme is established to determine the quadrature coefficients. By a special choice of the numerical fluxes and projections, we obtain {for the constant coefficient case } an optimal-order convergence result for the semi-discrete DG scheme. The temporal discretization is achieved by the standard two-step backward difference formula and a fast algorithm with linear complexity is constructed. Numerical examples are provided for demonstrating the efficiency of the proposed algorithm, validating the theoretical results and illustrating the behaviors of the energy., Comment: 24 pages, 6 figures, 3 tables

It is known that both linear and nonlinear optical phenomena can be produced when the plasmon in metallic nanostructures are excited by the external electromagnetic waves. In this work, a coupled system of Maxwell equations and a gas dynamic model including a quantum pressure term is employed to simulate the plasmon dynamics of free electron fluid in different metallic nanostructures using a discontinuous Galerkin method in two dimensions. Numerical benchmarks demonstrate that the proposed numerical method can simulate both the high order harmonic generation and the nonlocal effect from metallic nanostructures. Based on the switch-on-and-off investigation, we can conclude that the quantum pressure term in gas dynamics is responsible for the bulk plasmon resonance. In addition, for the dielectric-filled nano-cavity, a coupled effective polarization model is further adopted to investigate the optical behavior of bound electrons. Concerning the numerical setting in this work, a strengthened influence of bound electrons on the generation of high order harmonic waves has been observed., Comment: 17 pages, 11 figures

In this work, we investigate numerical solutions of the two-dimensional shallow water wave using a fully nonlinear Green-Naghdi model with an improved dispersive effect. For the purpose of numerics, the Green-Naghdi model is rewritten into a formulation coupling a pseudo-conservative system and a set of pseudo-elliptic equations. Since the pseudo-conservative system is no longer hyperbolic and its Riemann problem can only be approximately solved, we consider the utilization of the central discontinuous Galerkin method which possesses an important feature of needlessness of Riemann solvers. Meanwhile, the stationary elliptic part will be solved using the finite element method. Both the well-balanced and the positivity-preserving features which are highly desirable in the simulation of the shallow water wave will be embedded into the proposed numerical scheme. The accuracy and efficiency of the numerical model and method will be illustrated through numerical tests., Comment: 27 pages, 35 figures