No abstract uploaded!

Inspired by works of Castéras (Pac J Math 276:321–345, 2015), Li and Zhu (Calc Var Partial Differ Equ 58:1–18, 2019), Sun and Zhu (Calc Var Partial Differ Equ 60:1–26, 2021), we propose a heat flow for the mean field equation on a connected finite graph G=(V,E). Namely ∂_tϕ(u)=Δu−Q+ρ\frac{e^u}{∫_V e^u dμ} u(⋅,0)=u_0, where Δ is the standard graph Laplacian, ρ is a real number, Q:V→R is a function satisfying ∫_V Qdμ=ρ, and ϕ:R→R is one of certain smooth functions including ϕ(s)=es. We prove that for any initial data u_0 and any ρ∈R, there exists a unique solution u:V×[0,+∞)→R of the above heat flow; moreover, u(x, t) converges to some function u_∞:V→R uniformly in x∈V as t→+∞, and u_∞ is a solution of the mean field equation Δu_∞−Q+ρ\frac{e^{u_∞}}{∫_V e^{u_∞}dμ}=0. Though G is a finite graph, this result is still unexpected, even in the special case Q≡0. Our approach reads as follows: the short time existence of the heat flow follows from the ODE theory; various integral estimates give its long time existence; moreover we establish a Lojasiewicz–Simon type inequality and use it to conclude the convergence of the heat flow.

We consider the Schrödinger equation for the harmonic oscillator$i$$∂$_{$t$}$u$=$Hu$, where$H$=−Δ+|$x$|^{2}, with initial data in the Hermite–Sobolev space$H$^{−$s$/2}$L$^{2}(ℝ^{$n$}). We obtain smoothing and maximal estimates and apply these to perturbations of the equation and almost everywhere convergence problems.

This paper concerns the physical behaviors of any solutions to the one dimensional compressible Navier-Stokes equations for viscous and heat conductive gases with constant viscosities and heat conductivity for fast decaying density at far fields only. First, it is shown that the specific entropy becomes not uniformly bounded immediately after the initial time, as long as the initial density decays to vacuum at the far field at the rate not slower than $O\left(\frac1{|x|^{\ell_\rho}}\right)$ with $\ell_\rho>2$. Furthermore, for faster decaying initial density, i.e., $\ell_\rho\geq4$, a sharper result is discovered that the absolute temperature becomes uniformly positive at each positive time, no matter whether it is uniformly positive or not initially, and consequently the corresponding entropy behaves as $O(-\log(\varrho_0(x)))$ at each positive time, independent of the boundedness of the initial entropy. Such phenomena are in sharp contrast to the case with slowly decaying initial density of the rate no faster than $O(\frac1{x^2})$, for which our previous works \cite{LIXINADV,LIXINCPAM,LIXIN3DK} show that the uniform boundedness of the entropy can be propagated for all positive time and thus the temperature decays to zero at the far field. These give a complete answer to the problem concerning the propagation of uniform boundedness of the entropy for the heat conductive ideal gases and, in particular, show that the algebraic decay rate $2$ of the initial density at the far field is sharp for the uniform boundedness of the entropy. The tools to prove our main results are based on some scaling transforms, including the Kelvin transform, and a Hopf type lemma for a class of degenerate equations with possible unbounded coefficients.

We consider a system coupling the incompressible Navier–Stokes equations to the Vlasov–Fokker–Planck equation. Such a problem arises in the description of particulate flows. We design a numerical scheme to simulate the behavior of the system. This scheme is asymptotic-preserving, thus efficient in both the kinetic and hydrodynamic regimes. It has a numerical stability condition controlled by the non-stiff convection operator, with an implicit treatment of the stiff drag term and the Fokker–Planck operator. Yet, consistent to a standard asymptotic-preserving Fokker–Planck solver or an incompressible Navier–Stokes solver, only the conjugate–gradient method and fast Poisson and Helmholtz solvers are needed. Numerical experiments are presented to demonstrate the accuracy and asymptotic behavior of the scheme, with several interesting applications.

In this paper, we study the transonic shock problem for the full compressible Euler system in a general two-dimensional de Laval nozzle as proposed in Courant and Friedrichs (Supersonic flow and shock waves, Interscience, New York, 1948): given the appropriately large exit pressure pe(x), if the upstream flow is still supersonic behind the throat of the nozzle, then at a certain place in the diverging part of the nozzle, a shock front intervenes and the gas is compressed and slowed down to subsonic speed so that the position and the strength of the shock front are automatically adjusted such that the end pressure at the exit becomes pe(x). We solve this problem completely for a general class of de Laval nozzles whose divergent parts are small and arbitrary perturbations of divergent angular domains for the full steady compressible Euler system. The problem can be reduced to solve a nonlinear free boundary value problem for a mixed hyperbolic–elliptic system. One of the key ingredients in the analysis is to solve a nonlinear free boundary value problem in a weighted Hölder space with low regularities for a second order quasilinear elliptic equation with a free parameter (the position of the shock curve at one wall of the nozzle) and non-local terms involving the trace on the shock of the first order derivatives of the unknown function.

Let G=(V,E) be a locally finite graph. Firstly, using calculus of variations, including a direct method of variation and the mountain-pass theory, we get sequences of solutions to several local equations on G (the Schrödinger equation, the mean field equation, and the Yamabe equation). Secondly, we derive uniform estimates for those local solution sequences. Finally, we obtain global solutions by extracting convergent sequence of solutions. Our method can be described as a variational method from local to global.

We establish the global existence and uniqueness of classical solutions to the Cauchy problem for the isentropic compressible Navier-Stokes equations in three spatial dimensions with smooth initial data that are of small energy but possibly large oscillations with constant state as far field, which could be either vacuum or nonvacuum. The initial density is allowed to vanish, and the spatial measure of the set of vacuum can be arbitrarily large; in particular, the initial density can even have compact support. These results generalize previous results on classical solutions for initial densities being strictly away from vacuum and are the first for global classical solutions that may have large oscillations and can contain vacuum states. © 2012 Wiley Periodicals, Inc

For the three-dimensional full compressible Navier–Stokes system describing the motion of a viscous, compressible, heat-conductive, and Newtonian polytropic fluid, we establish the global existence and uniqueness of classical solutions with smooth initial data which are of small energy but possibly large oscillations where the initial density is allowed to vanish. Moreover, for the initial data, which may be discontinuous and contain vacuum states, we also obtain the global existence of weak solutions. These results generalize previous ones on classical and weak solutions for initial density being strictly away from a vacuum, and are the first for global classical and weak solutions which may have large oscillations and can contain vacuum states.

This paper analyzes the random fluctuations obtained by a heterogeneous multi-scale first-order finite element method applied to solve elliptic equations with a random potential. Several multiscale numerical algorithms have been shown to correctly capture the homogenized limit of solutions of elliptic equations with coefficients modeled as stationary and ergodic random fields. Because theoretical results are available in the continuum setting for such equations, we consider here the case of a secondorder elliptic equations with random potential in two dimensions of space. We show that the random fluctuations of such solutions are correctly estimated by the heterogeneous multi-scale algorithm when appropriate fine-scale problems are solved on subsets that cover the whole computational domain. However, when the fine-scale problems are solved over patches that do not cover the entire domain, the random fluctuations may or may not be estimated accurately. In the case of random potentials with short-range interactions, the variance of the random fluctuations is amplified as the inverse of the fraction of the medium covered by the patches. In the case of random potentials with long-range interactions, however, such an amplification does not occur and random fluctuations are correctly captured independent of the (macroscopic) size of the patches. These results are consistent with those obtained in [9] for more general equations in the one-dimensional setting and provide indications on the loss in accuracy that results from using coarser, and hence computationally less intensive, algorithms.

In this paper, we consider the cutoff Boltzmann equation near Maxwellian, we proved the global existence and uniqueness for the cutoff Boltzmann equation in polynomial weighted space for all $\gamma \in (-3, 1]$. We also proved initially polynomial decay for the large velocity in $L^2$ space will induce polynomial decay rate, while initially exponential decay will induce exponential rate for the convergence. Our proof is based on newly established inequalities for the cutoff Boltzmann equation and semigroup techniques. Moreover, by generalizing the $L_x^\infty L^1_v \cap L^\infty_{x, v}$ approach, we prove the global existence and uniqueness of a mild solution to the Boltzmann equation with bounded polynomial weighted $L^\infty_{x, v}$ norm under some small condition on the initial $L^1_x L^\infty_v$ norm and entropy so that this initial data allows large amplitude oscillations.

In this paper, we study entire solutions of the Allen-Chan equation in one-dimensional Euclidean space. This equation is a scalar reaction-diffusion equation with a bistable nonlinearity. It is well-known that this equation admits three different types of traveling fronts connecting two of its three constant states. Under certain conditions on the wave speeds, the existence of entire solutions with merging these three traveling fronts is shown by constructing a suitable pair of super-sub-solutions.

In this paper, we derive apriori estimates for constant scalar curvature K\"ahler metrics on a compact K\"ahler manifold. We show that higher order derivatives can be estimated in terms of a $C^0$ bound for the K\"ahler potential. We also discuss some local versions of these estimates which can be of independent interest.

We present a simple proof on the formation of flocking in the F. Cucker and S. Smale system [Jpn. J. Math. (3) 2, No. 1, 197–227 (2007; Zbl 1166.92323)] based on the explicit construction of a Lyapunov functional. Our results also provide a unified condition on the initial states in which the exponential convergence to the flocking state will occur. For large particle systems, we give a rigorous justification for the mean-field limit from the many particle Cucker-Smale system to the Vlasov equation with flocking dissipation as the number of particles goes to infinity.

No abstract uploaded!

The existence of topological solutions for the Chern-Simons equation with two Higgs particles has been proved by Lin, Ponce and Yang [16]. However, both the uniqueness problem and the existence of non-topological solutions have been left open. In this paper, we consider the case of one vortex at origin. Among others, we prove the uniqueness of topological solutions and give a complete study of the radial solutions, in particular, the existence of some non-topological solutions

We construct a conformally invariant random family of closed curves in the plane by welding of random homeomorphisms of the unit circle. The homeomorphism is constructed using the exponential of$βX$, where X is the restriction of the 2-dimensional free field on the circle and the parameter$β$is in the “high temperature” regime $$ \beta < \sqrt {2} $$ . The welding problem is solved by studying a non-uniformly elliptic Beltrami equation with a random complex dilatation. For the existence a method of Lehto is used. This requires sharp probabilistic estimates to control conformal moduli of annuli and they are proven by decomposing the free field as a sum of independent fixed scale fields and controlling the correlations of the complex dilatation restricted to dyadic cells of various scales. For the uniqueness we invoke a result by Jones and Smirnov on conformal removability of Hölder curves. Our curves are closely related to SLE($ϰ$) for$ϰ$<4.

We prove the existence of a resonance-free region in scattering by a strictly convex obstacle $\mathcal{O}$ with the Robin boundary condition $\partial_{\nu}u+\gamma u|_{\partial\mathcal{O}}=0$ . More precisely, we show that the scattering resonances lie below a cubic curve ℑ$ζ$=−$S$|$ζ$|^{1/3}+$C$. The constant$S$is the same as in the case of the Neumann boundary condition$γ$=0. This generalizes earlier results on cubic pole-free regions obtained for the Dirichlet boundary condition.

This paper investigates infinite-time spreading and finite-time blow-up for the Keller-Segel system. For 0 < m ≤ 2 − 2 / d, the L p space for both dynamic and steady solutions are detected with ${p:=\frac{d(2-m)}{2} }$ . Firstly, the global existence of the weak solution is proved for small initial data in L p . Moreover, when m > 1 − 2 / d, the weak solution preserves mass and satisfies the hyper-contractive estimates in L q for any p < q < ∞. Furthermore, for slow diffusion 1 < m ≤ 2 − 2/d, this weak solution is also a weak entropy solution which blows up at finite time provided by the initial negative free energy. For m > 2 − 2/d, the hyper-contractive estimates are also obtained. Finally, we focus on the L p norm of the steady solutions, it is shown that the energy critical exponent m = 2d/(d + 2) is the critical exponent separating finite L p norm and infinite L p norm for the steady state solutions.

The immersed boundary method is one of the most useful computational methods in studying fluid structure interaction. On the other hand, the Immersed Boundary method is also known to require small time steps to maintain stability when solved with an explicit method. Many implicit or approximately implicit methods have been proposed in the literature to remove this severe time step stability constraint, but none of them give satisfactory performance. In this paper, we propose an efficient semi-implicit scheme to remove this stiffness from the immersed boundary method for the Navier–Stokes equations. The construction of our semi-implicit scheme consists of two steps. First, we obtain a semiimplicit discretization which is proved to be unconditionally stable. This unconditionally stable semi-implicit scheme is still quite expensive to implement in practice. Next, we apply the small scale decomposition to the unconditionally stable semi-implicit scheme to construct our efficient semi-implicit scheme. Unlike other implicit or semi-implicit schemes proposed in the literature, our semi-implicit scheme can be solved explicitly in the spectral space. Thus the computational cost of our semi-implicit schemes is comparable to that of an explicit scheme. Our extensive numerical experiments show that our semiimplicit scheme has much better stability property than an explicit scheme. This offers a substantial computational saving in using the immersed boundary method.

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.

For periodic initial data with density allowed to vanish initially, we establish the global existence of strong and weak solutions to the two-dimensional barotropic compressible Navier–Stokes equations with no restrictions on the size of initial data provided the shear viscosity is a positive constant and the bulk one is λ =ρ^β with β>4/3. These results generalize and improve the previous ones due to Vaigant–Kazhikhov [Sib. Math. J. 36 (1995) 1283–1316] who required β>3. Moreover, we also prove that the densities for both the strong and weak solutions remain bounded from above independently of time. As a consequence, it is shown that both the strong and weak solutions converge to the equilibrium state as time tends to infinity.

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].

No abstract uploaded!

The Boltzmann equation which describes the time evolution of a large number of particles through the binary collision in statistics physics has close relation to the systems of fluid dynamics, that is, Euler equations and NavierStokes equations. As for a basic wave pattern to Euler equations, we consider the nonlinear stability of contact discontinuities to the Boltzmann equation. Even though the stability of the other two nonlinear waves, i.e., shocks and rarefaction waves has been extensively studied, there are few stability results on the contact discontinuity because unlike shock waves and rarefaction waves, its derivative has no definite sign, and decays slower than a rarefaction wave. Moreover, it behaves like a linear wave in a nonlinear setting so that its coupling with other nonlinear waves reveals a complicated interaction mechanism. Based on the new definition of contact waves to the Boltzmann equation