In this paper, we develop an energy stable fully discrete discontinuous Galekin (DG) finite element method for the thin film epitaxy problem. Based on the method of lines, we construct and prove the energy stability of the spatial semi-discrete DG scheme firstly. To avoid the strict time step restriction of the explicit time integration method, the first order convex splitting method is used to get an unconditionally stable fully discrete DG method, which is a linearly implicit scheme for this nonlinear problem. The energy stability of the fully discrete convex splitting DG scheme is also proved. To improve the temporal accuracy, spectral deferred correction (SDC) method is adapted to achieve the high order accuracy in both time and space. Combining with the convex splitting method, the SDC method can be linearly solvable, high order accurate and stable in our numerical tests. These advantages ensure that the resulting fully discrete DG scheme is efficient to perform the long time simulation of the thin film epitaxy model. Numerical experiments of the accuracy and long time simulation show the capability and efficiency of the method.

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.

No abstract uploaded!

No abstract uploaded!

A two-dimensional chiral conformal field theory can be viewed mathematically as the representation theory of its chiral algebra, a vertex operator algebra. Vertex operator algebras are especially well suited for studying logarithmic conformal field theory (in which correlation functions have logarithmic singularities arising from non-semisimple modules for the chiral algebra) because of the logarithmic tensor category theory of Huang, Lepowsky, and Zhang. In this paper, we study not-necessarily-semisimple or rigid braided tensor categories C of modules for the fixed-point vertex operator subalgebra V^G of a vertex operator (super)algebra V with finite automorphism group G. The main results are that every V^G-module in C with a unital and associative V-action is a direct sum of g-twisted V-modules for possibly several g∈G, that the category of all such twisted V-modules is a braided G-crossed (super)category, and that the G-equivariantization of this braided G-crossed (super)category is braided tensor equivalent to the original category C of V^G-modules. This generalizes results of Kirillov and Müger proved using rigidity and semisimplicity. We also apply the main results to the orbifold rationality problem: whether V^G is strongly rational if V is strongly rational. We show that V^G is indeed strongly rational if V is strongly rational, G is any finite automorphism group, and V^G is C_2-cofinite.