We propose a method to construct numerical solutions of parabolic equations on the unit sphere. The time discretization uses Laplace transforms and quadrature. The spatial approximation of the solution employs radial basis functions restricted to the sphere. The method allows us to construct high accuracy numerical solutions in parallel. We establish L2 error estimates for smooth and nonsmooth ...

Phase fields models offer a systematic physical approach for investigating complex multiphase systems behaviors such as near-critical interfacial phenomena, phase separation under shear, and microstructure evolution during solidification. However, because interfaces are replaced by thin transition regions (diffuse interfaces), phase field simulations require resolution of very thin layers to ca...

This paper presents a new numerical method for the compressible Navier-Stokes equations governing the flow of an ideal isentropic gas. To approximate the continuity equation, the method utilizes a discontinuous Galerkin discretization on piecewise constants and a basic upwind flux. For the momentum equation, the method is a new combined discontinuous Galerkin and finite element method approxima...

A new realization of a finite element level set method for simulation of immiscible fluid flows is introduced and validated on numerical benchmarks. The new method involves a mixed discretization of the dependent variables, discretizing the flow variables with non-conforming RannacherTurek finite elements while keeping a simple first order conforming discretization of the level set field. A thr...

In this paper, an alternating direction Galerkin finite element method is presented for solving 2D time fractional reaction sub-diffusion equation with nonlinear source term. Firstly, one order implicit-explicit method is used for time discretization, then Galerkin finite element method is adopted for spatial discretization and obtain a fully discrete linear system. Secondly, Galerkin alternati...

We deal with the numerical solution of time-dependent convection-diffusion-reaction equations. We combine the local projection stabilization method for the space discretization with two different time discretization schemes: the continuous Galerkin-Petrov (cGP) method and the discontinuous Galerkin (dG) method of polynomial of degree k. We establish the optimal error estimates and present numer...

We propose a simple fast spectral method for the Boltzmann collision operator with general collision kernels. Compared with the direct spectral method [8] which requires O(N6) memory to store precomputed weights and has O(N6) numerical complexity, the new method has complexity O(MN4 logN), where N is the number of discretization points in each velocity dimension, M is the total number discretiz...

To model periodic structures with oblique incident waves/scan angles in FDTD, the field transformation method is successfully used to analyze their characteristics. In the field transformation method, Maxwell’s equations are Floquet-transformed so that only a single period of infinite periodic structure can be modeled in FDTD by using periodic boundary conditions (PBCs). A new discretization me...

Efficient and accurate numerical methods are presented for computing ground states and dynamics of the three-dimensional (3D) nonlinear relativistic Hartree equation both without and with an external potential. This equation was derived recently for describing the mean field dynamics of boson stars. In its numerics, due to the appearance of pseudodifferential operator which is defined in phase ...

the railway transportation planning under the fuzzy environment is investigated in this paper. as a main result, a new modeling method, called minimum risk chance-constrained model, is presented based on the credibility measure. for the convenience ofs olving the mathematical model, the crisp equivalents ofc hance functions are analyzed under the condition that the involved fuzzy parameter...