Abhijit Biswas, Laila S. Busaleh, David I. Ketcheson, Carlos Muñoz-Moncayo, Manvendra Rajvanshi
arXiv preprint 2025
We study a first-order hyperbolic approximation of the nonlinear Schr\"odinger (NLS) equation. We show that the system is strictly hyperbolic and possesses a modified Hamiltonian structure, along with at least three conserved quantities that approximate those of NLS. We provide families of explicit standing-wave solutions to the hyperbolic system, which are shown to converge uniformly to ground-state solutions of NLS in the relaxation limit. The system is formally equivalent to NLS in the relaxation limit, and we develop asymptotic preserving discretizations that tend to a consistent discretization of NLS in that limit, while also conserving mass. Examples for both the focusing and defocusing regimes demonstrate that the numerical discretization provides an accurate approximation of the NLS solution.
Abhijit Biswas, Laila S. Busaleh, David I. Ketcheson, Carlos Muñoz-Moncayo, Manvendra Rajvanshi
arXiv preprint 2025
We study a first-order hyperbolic approximation of the nonlinear Schr\"odinger (NLS) equation. We show that the system is strictly hyperbolic and possesses a modified Hamiltonian structure, along with at least three conserved quantities that approximate those of NLS. We provide families of explicit standing-wave solutions to the hyperbolic system, which are shown to converge uniformly to ground-state solutions of NLS in the relaxation limit. The system is formally equivalent to NLS in the relaxation limit, and we develop asymptotic preserving discretizations that tend to a consistent discretization of NLS in that limit, while also conserving mass. Examples for both the focusing and defocusing regimes demonstrate that the numerical discretization provides an accurate approximation of the NLS solution.
Carlos Muñoz-Moncayo
arXiv preprint 2024
In this work, we address the imposition of outflow boundary conditions for one-dimensional conservation laws. While a highly accurate numerical solution can be obtained in the interior of the domain, boundary discretization can lead to unphysical reflections. We investigate and implement some classes of relaxation methods and far-field operators to deal with this problem without significantly increasing the size of the computational domain. Relations are established between these techniques, and extensions of them are explored. In particular, we introduce a simple and robust relaxation method with a matrix-valued weight function that selectively absorbs outgoing waves. As a challenging model problem, we consider the Lagrangian formulation of the Euler equations for an isotropic gas with inflow boundary conditions determined by an oscillating piston.
Carlos Muñoz-Moncayo
arXiv preprint 2024
In this work, we address the imposition of outflow boundary conditions for one-dimensional conservation laws. While a highly accurate numerical solution can be obtained in the interior of the domain, boundary discretization can lead to unphysical reflections. We investigate and implement some classes of relaxation methods and far-field operators to deal with this problem without significantly increasing the size of the computational domain. Relations are established between these techniques, and extensions of them are explored. In particular, we introduce a simple and robust relaxation method with a matrix-valued weight function that selectively absorbs outgoing waves. As a challenging model problem, we consider the Lagrangian formulation of the Euler equations for an isotropic gas with inflow boundary conditions determined by an oscillating piston.
Carlos Muñoz-Moncayo#, David I. Ketcheson, Manuel Quezada de Luna (# corresponding author)
Communications in Applied Mathematics and Computational Science 2023
We investigate the achievable efficiency of exact solvers for the Riemann problem for two systems of first-order hyperbolic PDEs: the shallow water equations and the Euler equations of compressible gas dynamics. Many approximate solvers have been developed for these systems; exact solution algorithms have received less attention because the computation of the exact solution typically requires an iterative solution of algebraic equations, which can be expensive or unreliable. We investigate a range of iterative algorithms and initial guesses. In addition to existing algorithms, we propose simple new algorithms that are guaranteed to converge and to remain in the range of physically admissible values at all iterations. We apply the existing and new iterative schemes to an ensemble of test Riemann problems. For the shallow water equations, we find that Newton’s method with a simple modification converges quickly and reliably. For the Euler equations we obtain similar results; however, when the required precision is high, a combination of Ostrowski and Newton iterations converges faster. These solvers are slower than standard approximate solvers like Roe and HLLE, but come within a factor of two in speed. We also provide a preliminary comparison of the accuracy of a finite volume discretization using an exact solver versus standard approximate solvers.
Carlos Muñoz-Moncayo#, David I. Ketcheson, Manuel Quezada de Luna (# corresponding author)
Communications in Applied Mathematics and Computational Science 2023
We investigate the achievable efficiency of exact solvers for the Riemann problem for two systems of first-order hyperbolic PDEs: the shallow water equations and the Euler equations of compressible gas dynamics. Many approximate solvers have been developed for these systems; exact solution algorithms have received less attention because the computation of the exact solution typically requires an iterative solution of algebraic equations, which can be expensive or unreliable. We investigate a range of iterative algorithms and initial guesses. In addition to existing algorithms, we propose simple new algorithms that are guaranteed to converge and to remain in the range of physically admissible values at all iterations. We apply the existing and new iterative schemes to an ensemble of test Riemann problems. For the shallow water equations, we find that Newton’s method with a simple modification converges quickly and reliably. For the Euler equations we obtain similar results; however, when the required precision is high, a combination of Ostrowski and Newton iterations converges faster. These solvers are slower than standard approximate solvers like Roe and HLLE, but come within a factor of two in speed. We also provide a preliminary comparison of the accuracy of a finite volume discretization using an exact solver versus standard approximate solvers.
Juan Mayorga-Zambrano#, Leonardo Medina-Espinosa, Carlos Muñoz-Moncayo (# corresponding author)
Differential Equations and Dynamical Systems 2023
We consider an N-dimensional nonlinear Schrödinger equation $(P_{\varepsilon})$ with a positive Planck constant $\varepsilon>0$ and power nonlinearity $p>1$. We consider the finite case and critical frequency as described by Byeon and Wang, i.e., the continuous non-negative potential's zero level set is a singleton, around which it decays like a homogeneous positive function. In the limit $\varepsilon \to 0$, we denote the semiclassical limit problem by $(P_{\text{lim}})$. By a Ljusternik-Schnirelman scheme we get an infinite number of solutions for $(P_{\varepsilon})$ and $(P_{\text{lim}})$, $v_{k,\varepsilon}$ and $w_k$, respectively. We prove, up to a scaling, that (a) $v_{k,\varepsilon}$ subconverges to $w_k$ pointwise and in Sobolev-like norm, (b) the energy of $v_{k,\varepsilon}$ converges to that of $w_k$, and (c) $v_{k,\varepsilon}$ exponentially decays out of the zero level set of the potential.
Juan Mayorga-Zambrano#, Leonardo Medina-Espinosa, Carlos Muñoz-Moncayo (# corresponding author)
Differential Equations and Dynamical Systems 2023
We consider an N-dimensional nonlinear Schrödinger equation $(P_{\varepsilon})$ with a positive Planck constant $\varepsilon>0$ and power nonlinearity $p>1$. We consider the finite case and critical frequency as described by Byeon and Wang, i.e., the continuous non-negative potential's zero level set is a singleton, around which it decays like a homogeneous positive function. In the limit $\varepsilon \to 0$, we denote the semiclassical limit problem by $(P_{\text{lim}})$. By a Ljusternik-Schnirelman scheme we get an infinite number of solutions for $(P_{\varepsilon})$ and $(P_{\text{lim}})$, $v_{k,\varepsilon}$ and $w_k$, respectively. We prove, up to a scaling, that (a) $v_{k,\varepsilon}$ subconverges to $w_k$ pointwise and in Sobolev-like norm, (b) the energy of $v_{k,\varepsilon}$ converges to that of $w_k$, and (c) $v_{k,\varepsilon}$ exponentially decays out of the zero level set of the potential.
Juan Riofrío#, Carlos Muñoz-Moncayo, Isidro R Amaro, Israel Pineda (# corresponding author)
Artificial Intelligence, Computer and Software Engineering Advances 2021
This paper presents a new clustering algorithm to identify groups of countries. First, a layer of several clustering methods is applied to the original dataset. Then, after performing dimensionality reduction techniques like t-SNE or SOM on the resulting data, a second clustering layer (K-Means) is applied to identify the final clusters. This method is applied to a dataset from 163 countries, considering the following variables population, area, Gross Domestic Product (GDP), Gross Domestic Product adjusted for Purchase Power Parity (GDP-PPP), and COVID-19 related data (Confirmed, Recovered, and Deaths). The implementation with SOM dimensionality reduction outperformed the one with t-SNE for the considered dataset. We expect that using this information, countries can have an insight on which measures against COVID-19 replicate or avoid, based on the results in countries from the same cluster.
Juan Riofrío#, Carlos Muñoz-Moncayo, Isidro R Amaro, Israel Pineda (# corresponding author)
Artificial Intelligence, Computer and Software Engineering Advances 2021
This paper presents a new clustering algorithm to identify groups of countries. First, a layer of several clustering methods is applied to the original dataset. Then, after performing dimensionality reduction techniques like t-SNE or SOM on the resulting data, a second clustering layer (K-Means) is applied to identify the final clusters. This method is applied to a dataset from 163 countries, considering the following variables population, area, Gross Domestic Product (GDP), Gross Domestic Product adjusted for Purchase Power Parity (GDP-PPP), and COVID-19 related data (Confirmed, Recovered, and Deaths). The implementation with SOM dimensionality reduction outperformed the one with t-SNE for the considered dataset. We expect that using this information, countries can have an insight on which measures against COVID-19 replicate or avoid, based on the results in countries from the same cluster.