2025-07-04 | | Total: 22
We introduce two multiscale numerical schemes for the time integration of weakly nonlinear Schrödinger equations, built upon the discretization of Picard iterates of the solution. These high-order schemes are designed to achieve high precision with respect to the small nonlinearity parameter under particular CFL condition. By exploiting the scattering properties of these schemes thanks to a low-frequency projected linear flow, we also establish its uniform accuracy over long time horizons. Numerical simulations are provided to illustrate the theoretical results, and these schemes are further applied to investigate dynamics in the framework of wave turbulence.
We address the initial source identification problem for the heat equation, a notably ill-posed inverse problem characterized by exponential instability. Departing from classical Tikhonov regularization, we propose a novel approach based on moment analysis of the heat flow, transforming the problem into a more stable inverse moment formulation. By evolving the measured terminal time moments backward through their governing ODE system, we recover the moments of the initial distribution. We then reconstruct the source by solving a convex optimization problem that minimizes the total variation of a measure subject to these moment constraints. This formulation naturally promotes sparsity, yielding atomic solutions that are sums of Dirac measures. Compared to existing methods, our moment-based approach reduces exponential error growth to polynomial growth with respect to the terminal time. We provide explicit error estimates on the recovered initial distributions in terms of moment order, terminal time, and measurement errors. In addition, we develop efficient numerical discretization schemes and demonstrate significant stability improvements of our approach through comprehensive numerical experiments.
We study the continuous model of the localized wave propagation corresponding to the one-dimensional diatomic crystal lattice. From the mathematical point of view the problem can be described in terms of the Cauchy problem with localized initial data for a system of two pseudo-differential equations. We assume two small parameters in this formulation -- the lattice step and the size if the initial perturbation. We construct the asymptotic solution of the continuous Cauchy problem with respect to the size of perturbation. The ratio of the small parameters drastically affects the form of the solution. We consider two situations -- when the size of the perturbation is sufficiently large and when it is comparable with the lattice step. In each situations we provide analytical formulae for the asymptotic solution via Airy function.
We consider the setting of two disks in a domain in R2 which are almost touching and have finite and positive conductivities, giving rise to a divergence form elliptic equation with discontinuous coefficients. We use the maximum principle to give a new proof of a gradient bound of Li-Vogelius.
The linearization principle states that the stability (or instability) of solutions to a suitable linearization of a nonlinear problem implies the stability (or instability) of solutions to the original nonlinear problem. In this work, we prove this principle for solutions of abstract fractional reaction-diffusion equations with a fractional derivative in time of order α∈(0,1). Then, we apply these results to particular fractional reaction-diffusion equations, obtaining, for example, the counterpart of the classical Turing instability in the case of fractional equations.
We consider the Cauchy problem for the modified Camassa-Holm equation ut+((u2−u2x)m)x=0,m\coloneqqu−uxx,t>0, −∞<x<+∞ subject to the step-like initial data: u(x,0)→A1 as x→−∞ and u(x,0)→A2 as x→+∞, where 0<A1<A2. The goal is to establish the global existence of the solution of this problem.
The Brenner-Navier-Stokes-Fourier (BNSF) system, introduced by Howard Brenner, was developed to address some deficiencies in the classical Navier-Stokes-Fourier system, based on the concept of volume velocity. We consider the one-dimensional BNSF system in Lagrangian mass coordinates, incorporating temperature-dependent transport coefficients, which yields a more physically realistic framework. We establish the existence and uniqueness of monotone traveling wave solutions (or viscous shocks) to the BNSF system with any positive C2 dissipation coefficients, provided that the shock amplitude is sufficiently small. We utilize geometric singular perturbation theory as in the constant coefficient case [13]; however, due to the arbitrary nonlinearities of the coefficients, we employ the implicit function theorem, which grants robustness to our approach. This work is motivated by [12], which proves a contraction property of any large solutions to the BNSF system around the traveling wave solutions. Thus, we also derive some quantitative estimates on the traveling wave solutions that play a fundamental role in [12].
We consider a Navier-Stokes fluid-plate interaction (FSI) system which describes the evolutions of the fluid contained within a 3D cavity, as it interacts with a deformable elastic membrane on the ``free" upper boundary of the cavity. These models arise in various aeroelastic and biomedical applications as well as in the control of ocular pressure, and sloshing phenomena. We analyze the well-posedness of weak solutions to the stationary (λ-parametrized) coupled PDE system by way of invoking the nonlinear generalization of the abstract variational formulations which was introduced in \cite{girault2012finite}, wherein an inf-sup approach is followed to show existence-uniqueness of solutions under a small data assumption. In addition, we provide a numerical approximation scheme of the infinite dimensional coupled system via a finite element method approximation (FEM). The numerical results use a standard conforming scheme and handle the introduced nonlinearities via Picard iterations. Numerical results are obtained for an appropriate test problem satisfying the necessary boundary conditions and coupling. Moreover, error bounds between the FEM and theoretical solution in terms of the characteristic mesh size are supplied in appropriate Sobolev norms which agree with the established literature. These FEM approximations of the coupled system with their associated error bounds validate the theoretical findings.
In this paper, we consider the Cauchy problem to the incompressible Euler and Navie-Stokes equations on the d-dimensional torus.Our aim of this paper is two fold. Firstly, we construct a new initial data and present a simple proof of the ill-posedness of the Euler equations in different senses: (1) the solution map of the Euler equations starting from u0 is discontinuous at t=0 in Bsp,∞ with s>0 and 1≤p≤∞, which covers the result obtained by Cheskidov and Shvydkoy in ;(2) the solution map of the Euler equations is not continuous as a map from Bsp,∞ to L∞T(Bsp,∞);(3) the solution map of the Euler equations cannot be Holder continuous in time variable in Besov spaces Bsp,r.
We establish the existence and provide explicit expressions for the stationary states of the one-dimensional Schrödinger equation with a repulsive delta-prime potential and a focusing nonlinearity of power type. Furthermore, we prove that, if the nonlinearity is subcritical, then ground states exist for any strength of the delta-prime interaction and for every positive value of the mass. This result supplies an example of ground states arising from a repulsive potential, a counterintuitive phenomenon explained by the emergence of an additional dimension in the energy space, induced by the delta-prime interaction. This new dimension contains states of lower energy and is thus responsible for the existence of nonlinear ground states that do not originate from linear eigenfunctions. The explicit form of the ground states is derived by addressing the ancillary problem of minimizing the action functional on the Nehari manifold. We solve such problem in the subcritical, critical, and supercritical regimes.
In this paper we show the non-uniqueness of mild solutions to the two-dimensional forced Navier-Stokes equations in the half space under the noslip boundary condition, following the program established by Albritton, Brué, and Colombo in 2022. Our construction of non-unique solutions is based on the instability of self-similar vorticity at high Reynolds numbers which concentrates around the boundary at the initial time. In our construction, therefore, a kind of boundary layer has to be taken into account in the analysis, contrasting to the known results where the unstable self-similar vorticity is located away from the boundary with O(1) distance around the initial time.
The aim of this paper is to define a field-road system in 2D where the road is a merely 1D-rectifiable set. For this purpose we introduce a general setting in order to define a parabolic problem onto a rectifiable set, which is coupled with another more classical parabolic problem outside this set, with transmission conditions.
We investigate the barotropic compressible Navier-Stokes equations with the Navier-slip boundary conditions in a general two-dimensional bounded simply connected domain. For initial density that is allowed to vanish, we establish the global existence and exponential decay of weak, strong, and classical solutions when the bulk viscosity coefficient is suitably large, without any restrictions on the size of the initial data. Furthermore, we prove that when the bulk viscosity coefficient tends to infinity, the solutions of the compressible Navier-Stokes equations converge to those of the inhomogeneous incompressible Navier-Stokes equations. The key idea is to utilize the logarithmic interpolation inequality on general bounded domains and apply the compensated compactness lemma.
Let Ω⊂R2 be a bounded domain on which Hardy's inequality holds. We prove that [exp(u2)−1]/δ2∈L1(Ω) if u∈H10(Ω), where δ denotes the distance to ∂Ω. The corresponding higher-dimensional result is also given. These results contain both Hardy's and Trudinger's inequalities, and yield a new variational characterization of the maximal solution of the Liouville equation on smooth domains, in terms of a renormalized functional. A global H1 bound on the difference between the maximal solution and the first term of its asymptotic expansion follows.
This paper investigates the Cauchy problem for the barotropic compressible Navier-Stokes equations in R2 with the constant state as far field, which could be vacuum or nonvacuum. Under the assumption of a sufficiently large bulk viscosity coefficient, we establish the long-time existence of weak, strong, and classical solutions, without imposing any extra restrictions on the initial velocity divergence. Moreover, we demonstrate that the solutions of the compressible Navier-Stokes equations converge to solutions of the inhomogeneous incompressible Navier-Stokes equations, as the bulk viscosity coefficient tends to infinity. The incompressible limit of the weak solutions holds even without requiring the initial velocity to be divergence-free. The key obstacle in the Cauchy problem is the failure of the Poincaré's inequality. This could be resolved by introducing a time-dependent Poincaré's type inequality, but it needs imposing weighted-integrability conditions on the initial density.
In this paper, we prove generalizations to the L^p setting of the Hardy-Rellich inequalities on domains of R^N with singularity given by the distance function to the boundary. The inequalities we obtain are either sharp in bounded domains, where we provide concrete minimizing sequences, or give a new bound for the sharp constant, while also depending on the geometric properties of the domain and its boundary. We also give applications to the existence and non-existence of solutions for a singular problem using variational methods and a Pohozaev identity.
The paper deals with homogenisation problems for high-contrast symmetric convolution-type operators with integrable kernels in media with a periodic microstructure. We adapt the two-scale convergence method to nonlocal convolution-type operators and obtain the homogenisation result both for problems stated in the whole space and in bounded domains with the homogeneous Dirichlet boundary condition. Our main focus is on spectral analysis. We describe the spectrum of the limit two-scale operator and characterize the limit behaviour of the spectrum of the original problem as the microstructure period tends to zero. It is shown that the spectrum of the limit operator is a subset the limit of the spectrum of the original operator, and that they need not coincide.
This paper is concerned with a 2D channel flow that is periodic horizontally but bounded above and below by hard walls. We assume the presence of horizontal viscosity only. We study the well-posedness, large-time behavior, and stability of solutions. For global well-posedness, we aim to assume less differentiability on initial velocity (u0,v0): in particular, we assume u0,v0∈L2(Ω) and ∂yu0∈L2(Ω).
This paper studies the inviscid limit problem for the two-dimensional Navier-Stokes equations with anisotropic viscosity. The fluid is assumed to be bounded above and below by impenetrable walls, with a no-slip boundary condition imposed on the bottom wall. For H2 initial velocity, we establish strong convergence in the Lp norm to the limiting problem as the vertical viscosity approaches zero, for any 2≤p<∞. The main challenge lies in the mismatch of boundary conditions - specifically, the no-slip condition in the original problem versus the slip condition in the limiting problem.
We provide the classical Boundary Harnack principle in Lipschitz domains for solutions to two different linear uniformly elliptic equations with the same principal part.
We develop an existence and regularity theory for solutions to a geometric free boundary problem motivated by models of tumor growth. In this setting, the tumor invades an accessible region D, its motion is directed along a constant vector V, and it cannot penetrate another region K acting as an obstacle to the spread of the tumor. Due to the non variational structure of the problem, we show existence of viscosity solutions via Perron's method. Subsequently, we prove interior regularity for the free boundary near regular points by means of an improvement of flatness argument. We further analyze the boundary regularity and we prove that the free boundary meets the obstacle as a C1,α graph. A key step in the analysis of the boundary regularity involves the study of a thin obstacle problem with oblique boundary conditions, for which we establish C1,α estimates.