Anisotropic Radial Basis Function Methods for Continental Size Ice Sheet Simulations

In this paper we develop and implement anisotropic radial basis function methods for simulating the dynamics of ice sheets and glaciers. We test the methods on two problems: the well-known benchmark ISMIP-HOM B that corresponds to a glacier size ice and a synthetic ice sheet whose geometry is inspired by the EISMINT benchmark that corresponds to a continental size ice sheet. We illustrate the advantages of the radial basis function methods over a standard finite element method. We also show how the use of anisotropic radial basis functions allows for accurate simulation of the velocities on a large ice sheet, which was not possible with standard isotropic radial basis function methods due to a large aspect ratio between the ice length and the ice thickness. Additionally, we implement a partition of unity method in order to improve the computational efficiency of the radial basis function methods.Comment: The authors contributed equally to this wor

Glueball relevant study on isoscalars from Nf=2N_f=2 lattice QCD

We perform a glueball-relevant study on isoscalars based on anisotropic $N_f=2$ lattice QCD gauge configurations. In the scalar channel, we identify the ground state obtained through gluonic operators to be a single-particle state through its dispersion relation. When $q\bar{q}$ operator is included, we find the mass of this state does not change, and the $q\bar{q}$ operator couples very weakly to this state. So this state is most likely a glueball state. For pseudoscalars, along with the exiting lattice results, our study implies that both the conventional $q\bar{q}$ state $\eta_2$ (or $\eta'$ in flavor $SU(3)$) and a heavier glueball-like state with a mass of roughly 2.6 GeV exist in the spectrum of lattice QCD with dynamical quarks.Comment: 8 pages, 3 figures, 3 tables, talk presented at the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spai