I am an Assistant Professor in the Department of Mathematics and Statistics at Oakland University, Michigan. Previously, I was a Postdoctoral Researcher in the Department of Mathematics at the University of Houston where I worked with Dr. Loic Cappanera. In 2020, I obtained my Ph.D. in Applied Mathematics at the University of Waterloo in Canada under the supervision of Dr. Sander Rhebergen.
My research interests lie in the design of finite element discretizations for the solution of partial differential equations with applications to realistic and challenging problems in science and engineering. For my projects, I implement novel discretization methods using open-source parallel codes to produce accurate and realistic simulations. Furthermore, I analyze the theoretical convergence of the methods under consideration, which involves rigorous mathematical proofs.
Highlights. The modeling of multiphase flows in porous media is a computationally and theoretically challenging problem due to the high nonlinear nature of the problem. I am interested in developing and analyzing discontinuous Galerkin discretizations for the three-phase (aqueous, liquid and vapor) case. Our paper https://doi.org/10.1093/imanum/drac053 was the first one to show existence and uniqueness of a DG method for three-phase flows in porous media, and we provided error estimates. Then, in https://doi.org/10.1016/j.amc.2024.128909 we developed a second-order time stepping method for the two-phase flow problem that preserves the energy exactly. We are currently extending this method to the black oil problem, in which phase changes occur.
Collaborators. Loic Cappanera (University of Ottawa), Catalin Trenchea (University of Pittsburgh), Maicon Correa (University of Sao Paulo), Mark Simmons (graduated from University of Houston).
Highlights. The phase-field crystal equation is a sixth-order nonlinear PDE that models microstructure evolution of two- phase systems (crystal and liquid) on atomic length and diffusive time scales. The high order in the equation leads to a large linear system. Moreover, higher order in time is highly desirable to be able to capture the evolution of the phases properly, especially for long time simulations. In https://doi.org/10.1016/j.apnum.2026.03.003 we designed a novel hybridizable discontinuous Galerkin method based on reformulation of the PDE as a first order system. We proved existence and uniqueness of a solution and our numerical results demonstrate that with a significantly smaller number of degrees of freedom, our method can achieve comparable results than the state of the art literature. Currently, our Doctoral student David Hennicken is working on a second-order time stepping method that preserves energy exactly.
Collaborators. Tamas Horvath (Oakland University), Natasha Sharma (University of Texas at El Paso), David Hennicken (current grad student at Oakland University).
Highlights. The modeling of welding of plates using lasers involves a heat equation with radiation boundary conditions. I recently published a paper https://authors.elsevier.com/sd/article/S0377-0427(26)00379-1 in which I showed that a previously noted issue in which the convergence of the method slows down in some cases, is not due to loss of regularity of the solution, but rather it is an artifact of the linearization scheme. I am currently working on obtaining optimal error estimates, and also on extending the method to include advection effects and a time-dependent problem.
G. Sosa Jones, M. Shillor, “A model for contact of a rod with an obstacle using the Damped Normal Compliance condition”. Accepted for publication in the IMA Journal of Applied Mathematics. 2026. Preprint
G. Saylor, “On the linearization scheme of an embedded discontinuous Galerkin method for the second-order elliptic problem with nonlinear Newton boundary conditions”. Accepted for publication in the Journal of Computational and Applied Mathematics. 2026.
G. Saylor, T. Horváth, N. Sharma, “An unconditionally stable hybridizable-embedded discontinuous Galerkin method for the phase field crystal equation”. Applied Numerical Mathematics Journal, 225, pp 241-258. 2026. Preprint
M. Shillor, G. Saylor, “Simulations of a dynamic rod-obstacle contact problem with the Damped Normal Compliance contact condition”. SeMA Journal, 82, pp 205-218. 2025.
G. Sosa Jones, C. Trenchea, “Discrete energy balance equation via a symplectic second-order method for two-phase flow in porous media”. Applied Mathematics and Computation, 480, 128909. 2024. Preprint
G. Sosa Jones, B. Riviere, L. Cappanera, “Existence and convergence of a discontinuous Galerkin method for the incompressible three-phase flow problem in porous media”. IMA Journal of Numerical Analysis, 43:5, pp 2714–2747, 2023. Preprint
G. Sosa Jones, S. Rhebergen, “An interface-tracking space-time hybridizable/embedded discontinuous Galerkin method for nonlinear free-surface flows” . Computers & Fluids, 246, 105625. 2022. Preprint
G. Sosa Jones, J.J. Lee, S. Rhebergen, “A space-time hybridizable discontinuous Galerkin method for linear free-surface waves” . Journal of Scientific Computing, 85:61. 2020. Preprint
G. Sosa Jones, J. Arteaga, O. Jiménez, “A study of mimetic and finite difference methods for the static diffusion equation” . Computers & Mathematics with Applications, 76/3, pp 633-648. 2018.
S. Buitrago, G. Sosa Jones, O. Jiménez, “An Upwind Finite Volume Method on Non-Orthogonal Quadrilateral Meshes for the Convection Diffusion Equation in Porous Media” . Applicable Analysis: An International Journal (Taylor & Francis Production), 95/10, pp 2203-2223. 2015.
G. Sosa Jones, “Space-time hybridizable discontinuous Galerkin methods for free-surface wave problems” . Ph.D. Thesis, Applied Mathematics, University of Waterloo. 2020.
Oakland University
University of Houston
University of Waterloo
Simón Bolívar University