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.

No abstract uploaded!

It is well-known that sparse grid algorithm has been widely accepted as an efficient tool to overcome the “curse of dimensionality" in some degree. In this note, we give the error estimate of hyperbolic cross (HC) approximations with all sorts of Askey polynomials. These polynomials are useful in generalized polynomial chaos (gPC) in the field of uncertainty quantification. The exponential convergences in both regular and optimized HC approximations have been shown under the condition that the random variable depends on the random inputs smoothly in some degree. Moreover, we apply gPC to numerically solve the ordinary differential equations with slightly higher dimensional random inputs. Both regular and optimized HC have been investigated with Laguerre-chaos, Charlier-chaos and Hermite-chaos in the numerical experiment. The discussion of the connection between the standard ANOVA approximation and Galerkin approximation is in the appendix.

In this paper we study the generalization of the Brocard congurations. The denition of q-period broken line is given and then the new necessary and sucient condition about the existence of two conjugate Brocard points of the q􀀀period loop broken line is investigated. Based on the past research, a geometric method for the construction of all possible loop broken lines with conjugate Brocard points is developed.

No abstract uploaded!