Source author record

Svetlana Roudenko

Svetlana Roudenko appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2025arXiv

Soliton profiles: Classical Numerical Schemes vs. Neural Network - Based Solvers

We present a comparative study of classical numerical solvers, such as Petviashvili's method or finite difference with Newton iterations, and neural network-based methods for computing ground states or profiles of solitary-wave solutions to the one-dimensional dispersive PDEs that include the nonlinear Schrödinger, the nonlinear Klein-Gordon and the generalized KdV equations. We confirm that classical approaches retain high-order accuracy and strong computational efficiency for single-instance problems in the one-dimensional setting. Physics-informed neural networks (PINNs) are also able to reproduce qualitative solutions but are generally less accurate and less efficient in low dimensions than classical solvers due to expensive training and slow convergence. We also investigate the operator-learning methods, which, although computationally intensive during training, can be reused across many parameter instances, providing rapid inference after pretraining, making them attractive for applications involving repeated simulations or real-time predictions. For single-instance computations, however, the accuracy of operator-learning methods remains lower than that of classical methods or PINNs, in general.

preprint2022arXiv

Interaction with an obstacle in the 2d focusing nonlinear Schrödinger equation

We present a numerical study of solutions to the $2d$ cubic and quintic focusing nonlinear Schrödinger equation in the exterior of a smooth, compact and strictly convex obstacle (a disk) with Dirichlet boundary condition. We first investigate the effect of the obstacle on the behavior of solutions traveling toward the obstacle at different angles and with different velocities directions. We introduce a new concept of weak and strong interactions of the solutions with the obstacle. Next, we study the existence of blow-up solutions depending on the type of the interaction and show how the presence of the obstacle changes the overall behavior of solutions (e.g., from blow-up to global existence), especially in the strong interaction case, as well as how it affects the shape of solutions compared to their initial data, (e.g., splitting into transmitted and reflected parts). We also investigate the influence of the size of the obstacle on the eventual existence of blow-up solutions in the strong interaction case in terms of the transmitted and the reflected parts of the mass. Moreover, we show that the sharp threshold for global existence vs. finite time blow-up solutions in the mass critical case in the presence of the obstacle is the same as the one given by Weinstein for {\rm{NLS}} in the whole Euclidean space $\R^d$. Finally, we construct new Wall-type initial data that blows up in finite time after a strong interaction with an obstacle and having a very distinct dynamics compared with all other blow-up scenarios and dynamics for the {\rm{NLS}} in the whole Euclidean space $\R^d$.

preprint2022arXiv

Well-posedness and dynamics of solutions to the generalized KdV with low power nonlinearity

We consider two types of the generalized Korteweg - de Vries equation, where the nonlinearity is given with or without absolute values, and, in particular, including the low powers of nonlinearity, an example of which is the Schamel equation. We first prove the local well-posedness of both equations in a weighted subspace of $H^1$ that includes functions with polynomial decay, extending the result of Linares et al [39] to fractional weights. We then investigate solutions numerically, confirming the well-posedness and extending it to a wider class of functions that includes exponential decay. We include a comparison of solutions to both types of equations, in particular, we investigate soliton resolution for the positive and negative data with different decay rates. Finally, we study the interaction of various solitary waves in both models, showing the formation of solitons, dispersive radiation and even breathers, all of which are easier to track in nonlinearities with lower power.

preprint2020arXiv

Asymptotic stability of solitary waves of the 3D quadratic Zakharov-Kuznetsov equation

We consider the quadratic Zakharov-Kuznetsov equation $$ \partial_t u + \partial_x Δu + \partial_x u^2 =0 $$ on $\mathbb{R}^3$. A solitary wave solution is given by $Q(x-t,y,z)$, where $Q$ is the ground state solution to $-Q + ΔQ + Q^2 =0$. We prove the asymptotic stability of these solitary wave solutions. Specifically, we show that initial data close to $Q$ in the energy space, evolves to a solution that, as $t\to\infty$, converges to a rescaling and shift of $Q(x-t,y,z)$ in $L^2$ in a rightward shifting region $x> δt -\tan θ\sqrt{y^2+z^2} $ for $0 \leq θ\leq \fracπ{3}-δ$.

preprint2020arXiv

Global behavior of solutions to the focusing generalized Hartree equation

We study the global behavior of solutions to the nonlinear generalized Hartree equation, where the nonlinearity is of the non-local type and is expressed as a convolution, $$ i u_t + Δu + (|x|^{-(N-γ)} \ast |u|^p)|u|^{p-2}u=0, \quad x \in \mathbb{R}^N, t\in \mathbb{R}. $$ Our main goal is to understand behavior of $H^1$ (finite energy) solutions of this equation in various settings. In this work we make an initial attempt towards this goal. We first investigate the $H^1$ local wellposedness and small data theory. We then, in the intercritical regime ($0<s<1$), classify the behavior of $H^1$ solutions under the mass-energy assumption $\mathcal{ME}[u_0]<1$, identifying the sharp threshold for global versus finite time solutions via the sharp constant of the corresponding convolution type Gagliardo-Nirenberg interpolation inequality (note that the uniqueness of a ground state is not known in the general case). In particular, depending on the size of the initial mass and gradient, solutions will either exist for all time and scatter in $H^1$, or blow up in finite time or diverge along an infinity time sequence. To either obtain $H^1$ scattering or divergence to infinity, in this paper we employ the well-known concentration compactness and rigidity method of Kenig-Merle [36] with the novelty of studying the nonlocal nonlinear potential given via convolution with negative powers of $|x|$ and different, including fractional, powers of nonlinearities.

preprint2020arXiv

Stable blow-up dynamics in the $L^2$-critical and $L^2$-supercritical generalized Hartree equation

We study stable blow-up dynamics in the generalized Hartree equation with radial symmetry, a Schrödinger-type equation with a nonlocal, convolution-type nonlinearity: $iu_t+Δu +\left(|x|^{-(d-2)} \ast |u|^{p} \right) |u|^{p-2}u = 0, x \in \mathbb{R}^d$. First, we consider the $L^2$-critical case in dimensions d=3, 4, 5, 6, 7 and obtain that a generic blow-up has a self-similar structure and exhibits not only the square root blowup rate but also the log-log correction (via asymptotic analysis and functional fitting). In this setting we also study blow-up profiles and show that generic blow-up solutions converge to the rescaled $Q$, a ground state solution of the elliptic equation $-ΔQ+Q- \left(|x|^{-(d-2)} \ast |Q|^p \right) |Q|^{p-2} Q =0$. We also consider the $L^2$-supercritical case in dimensions d=3,4. We derive the profile equation for the self-similar blow-up and establish the existence and local uniqueness of its solutions. As in the NLS $L^2$-supercritical regime, the profile equation exhibits branches of non-oscillating, polynomially decaying (multi-bump) solutions. A numerical scheme of putting constraints into solving the corresponding ODE is applied during the process of finding the multi-bump solutions. Direct numerical simulation of solutions to the generalized Hartree equation by the dynamic rescaling method indicates that the $Q_{1,0}$ is the profile for the stable blow-up. In this supercritical case, we obtain the blow-up rate without any correction. This blow-up happens at the focusing level $10^{-5}$, and thus, numerically observable (unlike the $L^2$-critical case). In summary, we find that the results are similar to the behavior of stable blowup dynamics in the corresponding NLS settings. Consequently, one may expect that the form of the nonlinearity in the Schrödinger-type equations is not essential in stable blow-up.

preprint2014arXiv

Going beyond the threshold: scattering and blow-up in the focusing NLS equation

We study the focusing nonlinear Schrödinger equation in the $L^2$-supercritical regime with finite energy and finite variance initial data. We investigate solutions above the energy (or mass-energy) threshold. In our first result, we extend the known scattering versus blow-up dichotomy above that threshold for finite variance solutions in the energy-subcritical and energy-critical regimes, obtaining scattering and blow-up criteria for solutions with arbitrarily large mass and energy. As a consequence, we characterize the behavior of the ground state initial data modulated by a quadratic phase. Our second result gives two blow up criteria, which are also applicable in the energy-supercritical NLS setting. We finish with various examples illustrating our results.

preprint2014arXiv

Modeling the Influence of Environment and Intervention on Cholera in Haiti

We propose a simple model with two infective classes in order to model the cholera epidemic in Haiti. We include the impact of environmental events (rainfall, temperature and tidal range) on the epidemic in the Artibonite and Ouest regions by introducing terms in the transmission rate that vary with environmental conditions. We fit the model on weekly data from the beginning of the epidemic until December 2013, including the vaccination programs that were recently undertaken in the Ouest and Artibonite regions. We then modified these projections excluding vaccination to assess the programs' effectiveness. Using real-time daily rainfall, we found lag times between precipitation events and new cases that range from 3.4 to 8.4 weeks in Artibonite and 5.1 to 7.4 in Ouest. In addition, it appears that, in the Ouest region, tidal influences play a significant role in the dynamics of the disease. Intervention efforts of all types have reduced case numbers in both regions; however, persistent outbreaks continue. In Ouest, where the population at risk seems particularly besieged and the overall population is larger, vaccination efforts seem to be taking hold more slowly than in Artibonite, where a smaller core population was vaccinated. The models including the vaccination programs predicted that a year and six months later, the mean number of cases in Artibonite would be reduced by about two thousand cases, and in Ouest by twenty four hundred cases below that predicted by the models without vaccination. We also found that vaccination is best when done in the early spring, and as early as possible in the epidemic. Comparing vaccination between the first spring and the second, there is a drop of about 40% in the case reduction due to the vaccine and about 10% per year after that.

preprint2012arXiv

A solution to the focusing 3d NLS that blows up on a contracting sphere

We rigorously construct radial $H^1$ solutions to the 3d cubic focusing NLS equation $i\partial_t ψ+ Δψ+ 2 |ψ|^2ψ=0$ that blow-up along a contracting sphere. With blow-up time set to $t=0$, the solutions concentrate on a sphere at radius $\sim t^{1/3}$ but focus towards this sphere at the faster rate $\sim t^{2/3}$. Such dynamics were originally proposed heuristically by Degtyarev-Zakharov-Rudakov (1975) and independently later in Holmer-Roudenko (2007), where it was demonstrated to be consistent with all conservation laws of this equation. In the latter paper, it was proposed as a solution that would yield divergence of the $L_x^3$ norm within the "wide" radius $\sim |\nabla u(t)|_{L_x^2}^{-1/2}$ but not within the "tight" radius $\sim |\nabla u(t)|_{L_x^2}^{-2}$, the second being the rate of contraction of self-similar blow-up solutions observed numerically and described in detail in Chapter 7 of Sulem-Sulem (1998).

preprint2010arXiv

A class of solutions to the 3d cubic nonlinear Schroedinger equation that blow-up on a circle

We consider the 3d cubic focusing nonlinear Schroedinger equation (NLS) i\partial_t u + Δu + |u|^2 u=0, which appears as a model in condensed matter theory and plasma physics. We construct a family of axially symmetric solutions, corresponding to an open set in H^1_{axial}(R^3) of initial data, that blow-up in finite time with singular set a circle in xy plane. Our construction is modeled on Raphaël's construction \cite{R} of a family of solutions to the 2d quintic focusing NLS, i\partial_t u + Δu + |u|^4 u=0, that blow-up on a circle.

preprint2010arXiv

Blow-up solutions on a sphere for the 3d quintic NLS in the energy space

We prove that if $u(t)$ is a log-log blow-up solution, of the type studied by Merle-Raphaël (2001-2005), to the $L^2$ critical focusing NLS equation $i\partial_t u +Δu + |u|^{4/d} u=0$ with initial data $u_0\in H^1(\mathbb{R}^d)$ in the cases $d=1, 2$, then $u(t)$ remains bounded in $H^1$ away from the blow-up point. This is obtained without assuming that the initial data $u_0$ has any regularity beyond $H^1(\mathbb{R}^d)$. As an application of the $d=1$ result, we construct an open subset of initial data in the radial energy space $H^1_{rad}(\mathbb{R}^3)$ with corresponding solutions that blow-up on a sphere at positive radius for the 3d quintic ($\dot H^1$-critical) focusing NLS equation $i\partial_tu + Δu + |u|^4u=0$. This improves Raphaël-Szeftel (2009), where an open subset in $H^3_{rad}(\mathbb{R}^3)$ is obtained. The method of proof can be summarized as follows: on the whole space, high frequencies above the blow-up scale are controlled by the bilinear Strichartz estimates. On the other hand, outside the blow-up core, low frequencies are controlled by finite speed of propagation.

preprint2010arXiv

Divergence of infinite-variance nonradial solutions to the 3d NLS equation

We consider solutions $u(t)$ to the 3d NLS equation $i\partial_t u + Δu + |u|^2u=0$ such that $\|xu(t)\|_{L^2} = \infty$ and $u(t)$ is nonradial. Denoting by $M[u]$ and $E[u]$, the mass and energy, respectively, of a solution $u$, and by $Q(x)$ the ground state solution to $-Q+ΔQ+|Q|^2Q=0$, we prove the following: if $M[u]E[u]<M[Q]E[Q]$ and $\|u_0\|_{L^2}\|\nabla u_0\|_{L^2}>\|Q\|_{L^2}\|\nabla Q\|_{L^2}$, then either $u(t)$ blows-up in finite positive time or $u(t)$ exists globally for all positive time and there exists a sequence of times $t_n\to +\infty$ such that $\|\nabla u(t_n)\|_{L^2} \to \infty$. Similar statements hold for negative time.

preprint2010arXiv

Maximizers for the Strichartz norm for small solutions of mass-critical NLS

Consider the mass-critical nonlinear Schrödinger equations in both focusing and defocusing cases for initial data in $L^2$ in space dimension N. By Strichartz inequality, solutions to the corresponding linear problem belong to a global $L^p$ space in the time and space variables, where $p=2+4/N$. In 1D and 2D, the best constant for the Strichartz inequality was computed by D.~Foschi who has also shown that the maximizers are the solutions with Gaussian initial data. Solutions to the nonlinear problem with small initial data in $L^2$ are globally defined and belong to the same global $L^p$ space. In this work we show that the maximum of the $L^p$ norm is attained for a given small mass. In addition, in 1D and 2D, we show that the maximizer is unique and obtain a precise estimate of the maximum. In order to prove this we show that the maximum for the linear problem in 1D and 2D is nondegenerated.