From a variational perspective, we derive a series of magnetization and quantum spin current systems coupled via an “s-d” potential term, including the Schrödinger--Landau--Lifshitz--Maxwell system, the Pauli--Landau--Lifshitz system, and the Schrödinger--Landau--Lifshitz system with successive simplifications. For the latter two systems, we propose using the time splitting spectral method for the quantum spin current and the Gauss--Seidel projection method for the magnetization. Accuracy of the time splitting spectral method applied to the Pauli equation is analyzed and verified by numerous examples. Moreover, behaviors of the Schrödinger--Landau--Lifshitz system in different “s-d” coupling regimes are explored numerically.

In this paper, we discuss error estimates associated with three different aggregation-diffusion splitting schemes for the Keller-Segel equations. We start with one algorithm based on the Trotter product formula, and we show that the convergence rate is CΔt, where Δt is the time-step size. Secondly, we prove the convergence rate CΔt2 for the Strang's splitting. Lastly, we study a splitting scheme with the linear transport approximation, and prove the convergence rate CΔt.

In this paper we study the existence assertion of the initial boundary value problem for the equation $\frac{\partial u}{\partial t} = \Delta e^{-\Delta u}$. This problem arises in the mathematical description of the evolution of crystal surfaces. Our investigations reveal that the exponent in the equation can have a singular part in the sense of the Lebesgue decomposition theorem, and the exponential nonlinearity somehow “cancels” it out. The net result is that we obtain a solution $u$ that satisfies the equation and the initial boundary conditions in the almost everywhere (a.e.) sense.

Abstract In this paper, we provide an alternative proof for the classical Sz. Nagy inequality in one dimension by a variational method and generalize it to higher dimensions J(h):=(∫Rd|h| dx)a611∫Rd|63h|2 dx(∫Rd|h|m+1 dx)a+1m+1≥β0, where m > 0 for d = 1, 2, for , and . The Euler–Lagrange equation for critical points of in the non-negative radial decreasing function space is given by a free boundary problem for a generalized Lane–Emden equation, which has a unique solution (denoted by h c ) and the solution determines the best constant for the above generalized Sz. Nagy inequality. The connection between the critical mass for the thin-film equation and the best constant of the Sz. Nagy inequality in one dimension was first noted by Witelski et al (2004 Eur. J. Appl. Math. 15 223–56). For the following critical thin film equation in multi-dimension ht+6366(h 63Δh)+6366(h 63hm)=0,x∈Rd, where m = 1 + 2/d, the critical mass is also given by . A finite time blow-up occurs for solutions with the initial mass larger than M c . On the other hand, if the initial mass is less than M c and a global non-negative entropy weak solution exists, then the second moment goes to infinity as or in for some subsequence . This shows that a part of the mass spreads to infinity.

This paper investigates the existence of a uniform in time $L^{\infty}$ bounded weak entropy solution for the quasilinear parabolic-parabolic Keller-Segel model with the supercritical diffusion exponent $0 < m < 2-\frac{2}{d}$ in the multi-dimensional space ${\mathbb{R}}^d$ under the condition that the $L^{\frac{d(2-m)}{2}}$ norm of initial data is smaller than a universal constant. Moreover, the weak entropy solution $u(x,t)$ satisfies mass conservation when $m > 1-\frac{2}{d}$. We also prove the local existence of weak entropy solutions and a blow-up criterion for general $L^1\cap L^{\infty}$ initial data.

In this paper, we introduce a random particle blob method for the Keller-Segel equation (with dimension ) and establish a rigorous convergence analysis.

We study the stability of traveling waves of the nonlinear Schrödinger equation with nonzero condition at infinity obtained via a constrained variational approach. Two important physical models for this are the Gross–Pitaevskii (GP) equation and the cubic-quintic equation. First, under a non-degeneracy condition we prove a sharp instability criterion for 3D traveling waves of (GP), which had been conjectured in the physical literature. This result is also extended for general nonlinearity and higher dimensions, including 4D (GP) and 3D cubic-quintic equations. Second, for cubic-quintic type nonlinearity, we construct slow traveling waves and prove their nonlinear instability in any dimension. For dimension two, the non-degeneracy condition is also proved for these slow traveling waves. For general traveling waves without vortices (that is nonvanishing) and with general nonlinearity in any dimension, we find a sharp condition for linear instability. Third, we prove that any 2D traveling wave of (GP) is transversally unstable, and we find the sharp interval of unstable transversal wave numbers. Near unstable traveling waves of all of the above cases, we construct unstable and stable invariant manifolds.

We investigate the long-time dynamics of an opinion formation model inspired by a work by Borghesi, Bouchaud and Jensen. First, we derive a Fokker–Planck-type equation under the assumption that interactions between individuals produce little consensus of opinion (grazing collision approximation). Second, we study conditions under which the Fokker–Planck equation has non-trivial equilibria and derive the macroscopic limit (corresponding to the long-time dynamics and spatially localized interactions) for the evolution of the mean opinion. Finally, we compare two different types of interaction rates: the original one given in the work of Borghesi, Bouchaud and Jensen (symmetric binary interactions) and one inspired from works by Motsch and Tadmor (non-symmetric binary interactions). We show that the first case leads to a conservative model for the density of the mean opinion whereas the second case leads to a non-conservative equation. We also show that the speed at which consensus is reached asymptotically for these two rates has fairly different density dependence.

In this paper, we study existence of global entropy weak solutions to a critical-case unstable thin lm equation in one-dimensional case ht + @x(h n @xxxh) + @x(h n+2@xh) = 0; where n  1. There exists a critical mass Mc = 2 p 6 3 found by Witelski etal. (2004 Euro. J. of Appl. Math. 15, 223{256) for n = 1. We obtain global existence of a non-negative entropy weak solution if initial mass is less than Mc. For n  4, entropy weak solutions are positive and unique. For n = 1,a nite time blow-up occurs for solutions with initial mass larger than Mc.For the Cauchy problem with n = 1 and initial mass less than Mc, we show that at least one of the following long-time behavior holds: the second moment goes to innity as the time goes to innity or h(
; tk) * 0 in L1 (R) for some subsequence tk ! 1.

We establish an optimal error estimate for a random particle blob method for the Keller-Segel equation in Rd (d  2). With a blob size " = h

No abstract uploaded!

No abstract uploaded!

ABSTRACT We study coagulation-fragmentation equations inspired by a simple model proposed in fisheries science to explain data for the size distribution of schools of pelagic fish. Although the equations lack detailed balance and admit no $H$-theorem, we are able to develop a rather complete description of equilibrium profiles and large-time behavior, based on recent developments in complex function theory for Bernstein and Pick functions. In the large-population continuum limit, a scaling-invariant regime is reached in which all equilibria are determined by a single scaling profile. This universal profile exhibits power-law behavior crossing over from exponent $-\frac23$ for small size to $-\frac32$ for large size, with an exponential cut-off.

This work considers the rigorous derivation of continuum models of step motion starting from a mesoscopic Burton-Cabrera-Frank (BCF) type model following the work [Xiang, SIAM J. Appl. Math. 2002]. We prove that as the lattice parameter goes to zero, for a finite time interval, a modified discrete model converges to the strong solution of the limiting PDE with first order convergence rate.

We study in this work a continuum model derived from 1D attachment-detachment-limited (ADL) type step flow on vicinal surface, $ u_t=-u^2(u^3)_{hhhh}$, where $u$, considered as a function of step height $h$, is the step slope of the surface. We formulate a notion of weak solution to this continuum model and prove the existence of a global weak solution, which is positive almost everywhere. We also study the long time behavior of weak solution and prove it converges to a constant solution as time goes to infinity. The space-time H\"older continuity of the weak solution is also discussed as a byproduct.

As V. I. Arnold observed in the 1960s, the Euler equations of incompressible fluid flow correspond formally to geodesic equations in a group of volume-preserving diffeomorphisms. Working in an Eulerian framework, we study incompressible flows of shapes as critical paths for action (kinetic energy) along transport paths constrained to be shape densities (characteristic functions). The formal geodesic equations for this problem are Euler equations for incompressible, inviscid potential flow of fluid with zero pressure and surface tension on the free boundary. The problem of minimizing this action exhibits an instability associated with microdroplet formation, with the following outcomes: Any two shapes of equal volume can be approximately connected by an Euler spray---a countable superposition of ellipsoidal geodesics. The infimum of the action is the Wasserstein distance squared, and is almost never attained except in dimension 1. Every Wasserstein geodesic between bounded densities of compact support provides a solution of the (compressible) pressureless Euler system that is a weak limit of (incompressible) Euler sprays. Each such Wasserstein geodesic is also the unique minimizer of a relaxed least-action principle for a two-fluid mixture theory corresponding to incompressible fluid mixed with vacuum.

We study convolution groups generated by completely monotone sequences and completely monotone functions. Using a convolution group, we define a fractional calculus for a certain class of distributions. When acting on causal functions, this definition agrees with the traditional Riemann-Liouville definition for t>0 but includes some singularities at t=0 so that the group property holds. Using this group, we are able to extend the definition of Caputo derivatives of order in (0,1) to a certain class of locally integrable functions without using the first derivative. The group property allows us to de-convolve the fractional differential equations to integral equations with completely monotone kernels, which then enables us to prove the general Gronwall inequality (or comparison principle) with the most general conditions. This then opens the door of a priori energy estimates of fractional PDEs. Some other fundamental results for fractional ODEs are also established within this frame under very weak conditions. Besides, we also obtain some interesting results about completely monotone sequences.

Abstract We study a class of fourth order nonlinear parabolic equations which include the thin-film equation and the quantum drift-diffusion model as special cases. We investigate these equations by first developing functional inequalities of the type $$ \int_\Omega u^{2\gamma-\alpha-\beta}\Delta\ua\Delta\ub dx \geq c\int_\Omega|\Delta \ur |^2dx, $$ which seem to be of interest on their own right.

A longstanding question in the dual Brunn–Minkowski theory is “What are the dual analogues of Federer’s curvature measures for convex bodies?” The answer to this is provided. This leads naturally to dual versions of Minkowski-type problems: What are necessary and sufficient conditions for a Borel measure to be a dual curvature measure of a convex body? Sufficient conditions, involving measure concentration, are established for the existence of solutions to these problems.

In this paper we are concerned with the existence of a weak solution to the initial boundary value problem for the equation ∂ t u=Δ(Δu) −3 . This problem arises in the mathematical modeling of the evolution of a crystal surface. Existence of a weak solution u with Δu≥0 is obtained via a suitable substitution. Our investigations reveal the close connection between this problem and the equation ∂ t ρ+ρ 2 Δ 2 ρ 3 =0 , another crystal surface model first proposed by H. Al Hajj Shehadeh, R. V. Kohn and J. Weare in [1].

Motivated by a model proposed by Peng et al. [Advances in Coll. and Interf. Sci. 206 (2014)] for break-up of tear films on human eyes, we study the dynamics of a generalized thin film model. The governing equations form a fourth-order coupled system of nonlinear parabolic PDE for the film thickness and salt concentration subject to nonconservative effects representing evaporation. We analytically prove the global existence of solutions to this model with mobility exponents in several different ranges and the results are then validated against PDE simulations. We also numerically capture other interesting dynamics of the model, including finite-time rupture-shock phenomenon due to the instabilities caused by locally elevated evaporation rates, convergence to equilibrium and infinite-time thinning.

In this paper we study classification of homogeneous solutions to the stationary Euler equation with locally finite energy. Written in the form $\bu = \n^\perp \Psi$, $\Psi(r,\th) = r^{\l} \psi(\th)$, for $\l >0$, we show that only trivial solutions exist in the range $0<\l<1/2$, i.e. parallel shear and rotational flows. In other cases many new solutions are exhibited that have hyperbolic, parabolic and elliptic structure of streamlines. In particular, for $\l>9/2$ the number of different non-trivial elliptic solutions is equal to the cardinality of the set $(2,\sqrt{2\l}) \cap \N$. The case $\l = 2/3$ is relevant to Onsager's conjecture. We underline the reasons why no anomalous dissipation of energy occurs for such solutions despite their critical Besov regularity $1/3$.

For autonomous Tonelli systems on $\R^n$, we develop an intrinsic proof of the existence of generalized characteristics using sup-convolutions. This approach, together with convexity estimates for the fundamental solution, leads to new results such as the global propagation of singularities along generalized characteristics.

We address the topology of the set of singularities of a solution of a Hamilton-Jacobi equation. For this we will apply the idea of the first two authors (Generalized characteristic and Lax-Oleinik operators: global result, preprint, arXiv:1605.07581, 2016.) to use the positive Lax-Oleinik semi-group to propagate singularities. {\it To cite this article: P. Cannarsa, W. Cheng, A. Fathi, C. R. Acad. Sci. Paris, Ser. I 340 (2005).

It is well-known that in Banach spaces with finite cotype, the $R$ -bounded and $\gamma $ -bounded families of operators coincide. If in addition $X$ is a Banach lattice, then these notions can be expressed as square function estimates. It is also clear that $R$ -boundedness implies $\gamma $ -boundedness. In this note we show that all other possible inclusions fail. Furthermore, we will prove that $R$ -boundedness is stable under taking adjoints if and only if the underlying space is $K$ -convex.