In this paper we consider the well-posedness of the compressible Prandtl boundary layer for the Navier–Stokes–Poisson equations with nonslip boundary condition and obtain the linear stability of the approximate solution constructed by boundary layer expansion.

The multi-scale zero relaxation singular limit for gas dynamics in thermal non-equilibrium with multiple non-equilibrium modes in multi-dimensions with physical boundaries from non-equilibrium to thermal equilibrium of compressible Euler flow is proved in this paper for analytical data by establishing the uniform local-in-time estimates. A cancellation mechanism is utilized to deal with the nonlinear singular terms that cause the increase in both time and space derivatives in energy estimates. The rates of the relaxations corresponding to different non-equilibrium modes tending to zero discussed in this paper can be arbitrarily different.

This work is concerned with our recently developed formalism of non-equilibrium thermodynamics. This formalism extends the classical irreversible thermodynamics which leads to classical thermodynamics and can not describe physical phenomena with long relaxations. With the extended theory, we obtain a generalized hydrodynamic system which can be used to describe the non-Newtonian fluids, heat transfer beyond Fourier's law and fluid flow in very short temporal and spatial scales. We study the mathematical structure of the generalized hydrodynamic system and show that it possesses a nice conservation-dissipation structure and therefore is symmetrizable hyperbolic. Moreover, we rigorously justify that this generalized hydrodynamic system tends to the classical hydrodynamics when the relaxation times approach to zero. This shows that the classical hydrodynamics can be derived as an approximation of our generalized one.

Edge states attract more and more research interests owing to the novel topologically protected properties. In this work, we studied edge modes and traveling edge states via the linear Dirac equation with so-called edge-admissible masses. The unidirectional edge state provides a heuristic approach to more general traveling edge states through the localized behavior along slowly varying edges. We show the dominated asymptotic solutions of two typical edge states that follow circular and curved edges with small curvature by the analytic and quantitative arguments.

In this article, we study wave dynamics in the fractional nonlinear Schrödinger equation with a modulated honeycomb potential. This problem arises from recent research interests in the interplay between topological materials and nonlocal governing equations. Both are current focuses in scientific research fields. We first develop the Floquet-Bloch spectral theory of the linear fractional Schrödinger operator with a honeycomb potential. Especially, we prove the existence of conical degenerate points, i.e., Dirac points, at which two dispersion band functions intersect. We then investigate the dynamics of wave packets spectrally localized at a Dirac point and derive the leading effective envelope equation. It turns out the envelope can be described by a nonlinear Dirac equation with a varying mass. With rigorous error estimates, we demonstrate that the asymptotic solution based on the effective envelope equation approximates the true solution well in the weighted-H^s space.

Mathematical analysis on electromagnetic waves in photonic graphene, a photonic topological material which has a honeycomb structure, is one of the most important current research topics. By modulating the honeycomb structure, numerous topological phenomena have been observed recently. The electromagnetic waves in such a media are generally described by the 2-dimensional wave equation. It has been shown that the corresponding elliptic operator with a honeycomb material weight has Dirac points in its dispersion surfaces. In this paper, we study the time evolution of the wave packets spectrally concentrated at such Dirac points in a modulated honeycomb material weight. We prove that such wave packet dynamics is governed by the Dirac equation with a varying mass in a large but finite time. Our analysis provides mathematical insights to those topological phenomena in photonic graphene.

In this paper, we establish the global C^{2,α} and W^{2,p} regularity for the Monge-Amp`ere equation det D^2u = f subject to boundary condition Du(Ω) = Ω^∗, where Ω and Ω^∗ are bounded convex domains in the Euclidean space R^n with C^{1,1} boundaries, and f is a Ho ̈lder continuous function. This boundary value problem arises naturally in optimal transportation and many other applications.

In this paper we study the distribution of scattering resonances for a multi-dimensional semi-classical Schrödinger operator, associated to a potential well in an island at energies close to the maximal one that limits the separation of the well and the surrounding sea.

In this article, we re-examine some of the classical pointwise multiplication theorems in Sobolev–Slobodeckij spaces, in part motivated by a simple counter-example that illustrates how certain multiplication theorems fail in Sobolev–Slobodeckij spaces when a bounded domain is replaced by R^n. We identify the source of the failure, and examine why the same failure is not encountered in Bessel potential spaces. To analyze the situation, we begin with a survey of the classical multiplication results stated and proved in the 1977 article of Zolesio, and carefully distinguish between the case of spaces defined on the all of R^n and spaces defined on a bounded domain (with e.g. a Lipschitz boundary). However, the survey we give has a few new wrinkles; the proofs we include are based almost exclusively on interpolation theory rather than Littlewood–Paley theory and Besov spaces, and some of the results we give and their proofs, including the results for negative exponents, do not appear in the literature in this form. We also include a particularly important variation of one of the multiplication theorems that is relevant to the study of nonlinear PDE systems arising in general relativity and other areas. The conditions for multiplication to be continuous in the case of Sobolev–Slobodeckij spaces are somewhat subtle and intertwined, and as a result, the multiplication theorems of Zolesio in 1977 have been cited (more than once) in the standard literature in slightly more generality than what is actually proved by Zolesio, and in cases that allow for construction of counter-examples such as the one included here.

For the class of graph-directed self-similar measures on R, which could have overlaps but are essentially of finite type, we set up a framework for deriving a closed formula for the spectral dimension of the Laplacians defined by these measures. For the class of finitely ramified graph-directed self-similar sets, the spectral dimension of the associated Laplace operators has been obtained by Hambly and Nyberg [6]. The main novelty of our results is that the graph-directed self-similar measures we consider do not need to satisfy the graph open set condition.

In [9] Anderson’s conjecture was proven by comparing values of Bloch functions with the variation of the function. We extend that result on Bloch functions from two to arbitrary dimension and prove that ∫_{[0,x]}|∇b(ζ)|e^{b(ζ)}d|ζ|<∞.. In the second part of the paper, we show that the area or volume integral ∫_{B^d}|∇u(w)|p(w,θ)dA(w) for positive harmonic functions u is bounded by the value cu(0) for at least one θ. The integral is also transferred to simply connected domains and interpreted from the point of view of stochastics. Several emerging open problems are presented.

We look for bounded periodic solutions for a parametric fractional problem involving a continuous nonlinearity with subcritical growth. By using a variant of Caffarelli and Silvestre extension method adapted to the periodic case and variational tools we prove the existence of at least three bounded periodic solutions when the parameter varies in an appropriate range.

In this paper we prove the following long-standing conjecture: stable solutions to semi-linear elliptic equations are bounded (and thus smooth) in dimension n⩽9. This result, that was only known to be true for n⩽4, is optimal: log(1/|x|^2) is a W^{1,2} singular stable solution for n⩾10. The proof of this conjecture is a consequence of a new universal estimate: we prove that, in dimension n⩽9, stable solutions are bounded in terms only of their L^1 norm, independently of the non-linearity. In addition, in every dimension we establish a higher integrability result for the gradient and optimal integrability results for the solution in Morrey spaces. As one can see by a series of classical examples, all our results are sharp. Furthermore, as a corollary, we obtain that extremal solutions of Gelfand problems are W^{1,2} in every dimension and they are smooth in dimension n⩽9. This answers to two famous open problems posed by Brezis and Brezis–Vázquez.

A key difficulty in the analysis and numerical approximation of the shallow water equations is the nonconservative product of measures due to the gravitational force acting on a sloped bottom. Solutions may be nonunique, and numerical schemes are not only consistent discretizations of the shallow water equations, but they also determine how to model the physics. Our derivation is based on a subcell reconstruction using infinitesimal singular layers at the cell boundaries, as inspired by S. Noelle, Y. Xing, and C.-W. Shu [J. Comput. Phys., 226 (2007), pp. 29-58]. One key step is to separate the singular measures. Another aspect is the reconstruction of the solution variables in the singular layers. We study three reconstructions. The first leads to the well-known scheme of Audusse et al., [SIAM J. Sci. Comput., 25 (2004), pp. 2050-2065], which introduces the hydrostatic reconstruction. The second is a modification proposed in [T. Morales de Luna, M. J. Castro D´ıaz, and C. Par´es, Appl. Math. Comput., 219 (2013), pp. 9012-9032], which analyzes whether a wave has enough energy to overcome a step. The third is our new scheme, which borrows its structure from the wet-dry front. For a number of cases discussed in recent years, where water runs down a hill, Audusse’s scheme converges slowly or fails. Morales’ scheme gives a visible improvement. Both schemes are clearly outperformed by our new scheme.

A novel discontinuous Galerkin (DG) method is developed to solve timedependent bi-harmonic type equations involving fourth derivatives in one and multiple space dimensions. We present the spatial DG discretization based on a mixed formulation and central interface numerical fluxes so that the resulting semi-discrete schemes are L^2 stable even without interior penalty. For time discretization, we use Crank–Nicolson so that the resulting scheme is unconditionally stable and second order in time. We present the optimal L^2 error estimate of O(h^{k+1}) for polynomials of degree k for semi-discrete DG schemes, and the L2 error of O(h^{k+1} + (\Delta t)^2) for fully discrete DG schemes. Extensions to more general fourth order partial differential equations and cases with non-homogeneous boundary conditions are provided. Numerical results are presented to verify the stability and accuracy of the schemes. Finally, an application to the one-dimensional Swift–Hohenberg equation endowed with a decay free energy is presented.

In this paper, we establish the local existence and uniqueness of multi-dimensional contact discontinuities for the ideal compressible magnetohydrodynamics (MHD) in Sobolev spaces, which are most typical interfacial waves for astrophysical plasmas and prototypical fundamental waves for hyperbolic systems of conservation laws. Such waves are characteristic discontinuities for which there is no flow across the discontinuity surface while the magnetic field crosses transversely, which lead to a two-phase free boundary problem where the pressure, velocity and magnetic field are continuous across the interface whereas the entropy and density may have jumps. To overcome the difficulties of possible nonlinear Rayleigh--Taylor instability and loss of derivatives, here we use crucially the Lagrangian formulation and the Cauchy formula (1882) for the magnetic field. This enables us to capture the boundary regularizing effect of the transversal magnetic field on the flow map, and on the other hand, it allows us to define a special good unknown of the magnetic field to get around the troublesome boundary integrals due to the transversality of the magnetic field. In particular, our result removes the additional assumption of the Rayleigh--Taylor sign condition required by Morando, Tarkhinin and Trebeschi (\emph{J. Differential Equations} \textbf{258} (2015), no. 7, 2531--2571; \emph{Arch. Ration. Mech. Anal.} \textbf{228} (2018), no. 2, 697--742) and holds for both 2D and 3D and hence gives a complete answer to the two open questions raised therein. Moreover, there is {\it no loss of derivatives} in our well-posedness theory. The solution is constructed as {\it the inviscid limit} of solutions to well-chosen nonlinear approximate problems for the two-phase compressible viscous non-resistive MHD.

Due to the highly degeneracy and singularities of the entropy equation, the physical entropy for viscous and heat conductive polytropic gases behave singularly in the presence of vacuum and it is thus a challenge to study its dynamics. It is shown in this paper that the uniform boundedness of the entropy and the inhomogeneous Sobolev regularities of the velocity and temperature can be propagated for viscous and heat conductive gases in $\mathbb R^3$, provided that the initial vacuum occurs only at far fields with suitably slow decay of the initial density. Precisely, it is proved that for any strong solution to the Cauchy problem of the heat conductive compressible Navier--Stokes equations, the corresponding entropy keeps uniformly bounded and the $L^2$ regularities of the velocity and temperature can be propagated, up to the existing time of the solution, as long as the initial density vanishes only at far fields with a rate no more than $O(\frac{1}{|x|^2})$. The main tools are some singularly weighted energy estimates and an elaborate De Giorgi type iteration technique. The De Giorgi type iterations are carried out to different equations in establishing the lower and upper bounds of the entropy.

In this paper, we consider the Cauchy problem to the planar non-resistive magnetohydrodynamic equations without heat conductivity, and establish the global well-posedness of strong solutions with large initial data. The key ingredient of the proof is to establish the a priori estimates on the effective viscous flux and a newly introduced ``transverse effective viscous flux" vector field inducted by the transverse magnetic field. The initial density is assumed only to be uniformly bounded and of finite mass and, in particular, the vacuum and discontinuities of the density are allowed.

In this paper, we provide rigorous justification of the hydrostatic approximation and the derivation of primitive equations as the small aspect ratio limit of the incompressible three-dimensional Navier-Stokes equations in the anisotropic horizontal viscosity regime. Setting $\varepsilon >0$ to be the small aspect ratio of the vertical to the horizontal scales of the domain, we investigate the case when the horizontal and vertical viscosities in the incompressible three-dimensional Navier-Stokes equations are of orders $O(1)$ and $O(\varepsilon^\alpha)$, respectively, with $\alpha>2$, for which the limiting system is the primitive equations with only horizontal viscosity as $\varepsilon$ tends to zero. In particular we show that for ``well prepared" initial data the solutions of the scaled incompressible three-dimensional Navier-Stokes equations converge strongly, in any finite interval of time, to the corresponding solutions of the anisotropic primitive equations with only horizontal viscosities, as $\varepsilon$ tends to zero, and that the convergence rate is of order $O\left(\varepsilon^\frac\beta2\right)$, where $\beta=\min\{\alpha-2,2\}$. Note that this result is different from the case $\alpha=2$ studied in [Li, J.; Titi, E.S.: \emph{The primitive equations as the small aspect ratio limit of the Navier-Stokes equations: Rigorous justification of the hydrostatic approximation}, J. Math. Pures Appl., \textbf{124} \rm(2019), 30--58], where the limiting system is the primitive equations with full viscosities and the convergence is globally in time and its rate of order $O\left(\varepsilon\right)$.

In this paper, we consider the initial boundary value problem in a cylindrical domain to the three dimensional primitive equations with full eddy viscosity in the momentum equations but with only horizontal eddy diffusivity in the temperature equation. Global well-posedness of $z$-weak solution is established for any such initial datum that itself and its vertical derivative belong to $L^2$. This not only extends the results in \cite{Cao5} from the spatially periodic case to general cylindrical domains but also weakens the regularity assumptions on the initial data which are required to be $H^2$ there.

In this paper, we investigate the initial-boundary value problem to the heat conductive compressible Navier-Stokes equations. Local existence and uniqueness of strong solutions is established with any such initial data that the initial density $\rho_0$, velocity $u_0$, and temperature $\theta_0$ satisfy $\rho_0\in W^{1,q}$, with $q\in(3,6)$, $u_0\in H^1$, and $\sqrt{\rho_0}\theta_0\in L^2$. The initial density is assumed to be only nonnegative and thus the initial vacuum is allowed. In addition to the necessary regularity assumptions, we do not require any initial compatibility conditions such as those proposed in (Y. Cho and H. Kim, \emph{Existence results for viscous polytropic fluids with vacuum}, J. Differential Equations {\bf 228} (2006), no.~2, 377--411.), which although are widely used in many previous works but put some inconvenient constraints on the initial data. Due to the weaker regularities of the initial data and the absence of the initial compatibility conditions, leading to weaker regularities of the solutions compared with those in the previous works, the uniqueness of solutions obtained in the current paper does not follow from the arguments used in the existing literatures. Our proof of the uniqueness of solutions is based on the following new idea of two-stages argument: (i) showing that the difference of two solutions (or part of their components) with the same initial data is controlled by some power function of the time variable; (ii) carrying out some singular-in-time weighted energy differential inequalities fulfilling the structure of the Gr\"onwall inequality. The existence is established in the Euler coordinates, while the uniqueness is proved in the Lagrangian coordinates first and then transformed back to the Euler coordinates.

We study linear and nonlinear Schrodinger equations defined by fractal measures. Under the assumption that the Laplacian has compact resolvent, we prove that there exists a uniqueness weak solution for a linear Schrodinger equation, and then use it to obtain the existence and uniqueness of weak solution of a nonlinear Schrodinger equation. We prove that for a class of self-similar measures on R with overlaps, the Schrodinger equations can be discretized so that the finite element method can be applied to obtain approximate solutions. We also prove that the numerical solutions converge to the actual solution and obtain the rate of convergence.

The weak well-posedness of strong damping wave equations defined by fractal Laplacians is proved by using Galerkin method. These fractal Laplacians are defined by self-similar measures with overlaps, such as the well-known infinite Bernoulli convolution associated with the golden ratio, the three-fold convolution of the Cantor measure, and a class of self-similar measures that we call essentially of finite type. In general, the structure of self-similar measures with overlap are complicated and intractable. However, some important information about the structure of the three measures above can be obtained. We make use of these information to set up a framework for one-dimensional measures to discretize the equations, and use the finite element and central difference methods to obtain numerical approximations of the weak solutions. We also show that the numerical solutions converge to the actual solution and obtain the rate of convergence.

No abstract uploaded!

We solve the Fu-Yau equation for arbitrary dimension and arbitrary slope $\alpha'$. Actually we obtain at the same time a solution of the open case $\alpha'>0$, an improved solution of the known case $\alpha'<0$, and solutions for a family of Hessian equations which includes the Fu-Yau equation as a special case. The method is based on the introduction of a more stringent ellipticity condition than the usual $\Gamma_k$ admissible cone condition, and which can be shown to be preserved by precise estimates with scale.