Computational & Applied Mathematics (CAM) Seminar
Student No.:40
Time:Tue 15:35, 2017.12.26
Instructor:Shi Zuoqiang  [Tsinghua University]
Place:Conference Room 1, Jinchun Yuan West Bldg.
Starting Date:2017-9-18
Ending Date:2018-1-5






Speaker: Chi-Wang Shu [Division of Applied Mathematics, Brown University]

Time:  10:00-11:00, 10th Dec. 2017 (Sunday) 周日上午 10:00
Location: Conference Room 1, Jin Chun Yuan West Bldg. 近春园西楼一楼第一会议室

Title: IMEX time marching for discontinuous Galerkin methods

 For discontinuous Galerkin methods approximating convection diffusion equations, explicit time marching is expensive since the time step is restricted by the square of the spatial mesh size. Implicit methods, however, would require the solution of non-symmetric, non-positive definite and nonlinear systems, which could be difficult.
The high order accurate implicit-explicit (IMEX) Runge-Kutta or multi-step time marching, which treats the diffusion term implicitly (often linear, resulting in a linear positive-definite solver) and the convection term (often nonlinear) explicitly, can greatly improve computational efficiency. We prove that certain IMEX time discretizations, up to third order accuracy, coupled with local discontinuous Galerkin method for the diffusion term treated implicitly, and regular discontinuous Galerkin method for the convection term treated explicitly, are unconditionally stable (the time step is upper bounded only by a constant depending on the diffusion coefficient but not on the spatial mesh size) and optimally convergent. The results have been generalized to multi-dimensional unstructured meshes, to other types of DG methods such as the embedded DG methods, to fourth order PDEs, and to incompressible fluid flow. The method has been applied to the drift-diffusion model in semiconductor device simulations, where a convection diffusion equation is coupled with an electrical potential equation. Numerical experiments confirm the good performance of such schemes. This is a joint work with Haijin Wang, Qiang Zhang, Yunxian Liu, Shiping Wang and Guosheng Fu.




Speaker: 曹外香 Cao Waixiang [北京师范大学 Beijing Normal University]

Time:  3:15pm, Nov. 28 (周二下午3:15)
Location: Conference Room 3, Jin Chun Yuan West Bldg. (近春园西楼二楼第三会议室)

Title: Superconvergence of Discontinuous Galerkin methods for linear hyperbolic equations

Abstract: In this talk, we will introduce superconvergence properties of discontinuous Galerkin (DG) methods linear hyperbolic conservation laws when upwind fluxes are used. We prove, under some suitable initial and boundary discretizations, the ($2k+1$)-th order superconvergence rate of the DG approximation at the downwind points and for the cell averages, Moreover, we prove that the gradient of the DG solution is superconvergent with a rate of ($k+1$)-th order at all interior left Radau points; and the function value approximation is superconvergent at all right Radau points with a rate of ($k+2$)-th order. Numerical experiments indicate that the aforementioned superconvergence rates are sharp.







Speaker: Li Guo, 中山大学

Time:  3:15pm, Nov. 14 (周二下午3:15)
Location: Conference Room 3, Jin Chun Yuan West Bldg. (近春园西楼二楼第三会议室)

Title: Energy stable local discontinuous Galerkin methods

Abstract: Energy stable local discontinuous Galerkin (LDG) method and its two applications have been investigated in this talk, including the Keller-Segel (KS) chemotaxis model and the nonlinear Schrodinger equation with wave operator (NLSW). The KS chemotaxis model may exhibit
blow-up patterns with certain initial conditions, and is not easy to approximate numerically. A free energy which is decreasing during time evolution has been constructed. Even though the energy can be negative, it is always positive if the blow-up will occur. Meanwhile, the NLSW
problem also has a conservative energy which is important to simulate long time behavior and eliminate oscillations. In this talk, we will construct a special energy stable LDG method to approximate the KS chemotaxis model with blow-up solution and numerically preserve the energy as well as a fully discrete energy conserving scheme utilizing the LDG method in space and the Crank-Nicholson algorithm in time to simulate the NLSW problem. Some numerical experiments for these two problems will be given to demonstrate the validity and performance of the energy stable LDG method.







Speaker: Su Xifeng [Beijing Normal University]


Time: 15:20-16:55, Nov. 2nd (周四下午)
Location: Conference Room 3, Jin Chun Yuan West Bldg. (近春园西楼二层第三会议室)



Title: Convergence of discrete Aubry-Mather model in the continuous limit and related topics


Abstract: We will consider the Frenkel-Kontorova models and their higher dimensional generalizations and talk about the corresponding discrete weak KAM theory. The existence of the discrete weak KAM solutions is related to the additive eigenvalue problem in ergodic optimization. In particular, I will show that the discrete weak KAM solutions converge to the weak KAM solutions of the autonomous Tonelli Hamilton-Jacobi equations as the time step goes to zero. This is a joint work with Prof. Philippe Thieullen.




Time:  3:15pm, Oct. 10 (周二下午3:15)
Location: Conference Room 3, Jin Chun Yuan West Bldg. (近春园西楼二楼第三会议室)


Title: Staggered Discontinuous Galerkin Methods for Stokes problem and elastodynamics


Speaker: Du Jie, YMSC


Abstract: Staggered discontinuous Galerkin (SDG) method is a new class of discontinuous Galerkin (DG) methods, which uses staggered mesh. It combines some good properties of finite element methods and standard DG methods through the use of staggered grid. In this talk, we develop an adaptive SDG method for the Stokes system. A computable error indicator is constructed and its reliability and efficiency are proved. Moreover, by combing the features of SDG method and traditional hybridization method, we present a staggered hybridization technique for DG methods to discretize linear elastodynamic equations. Our new approach offers several advantages, namely energy conservation, high-order optimal convergence, preservation of symmetry for the stress tensor, block diagonal mass matrices as well as low dispersion error.





Title: Self-Organized Hydrodynamic models for nematic alignment and the application to myxobacteria

Speaker: Yu Hui, YMSC

Abstract: A continuum model for a population of self-propelled particles interacting through nematic alignment is derived from an individual-based model. The methodology consists of introducing a hydrodynamic scaling of the corresponding mean field kinetic equation. The resulting perturbation problem is solved thanks to the concept of generalized collision invariants. It yields a hyperbolic but non-conservative system of equations for the nematic mean direction of the flow and the densities of particles flowing parallel or antiparallel to this mean direction. An application to myxobacteria is presented.



For the detail of CAM seminar, please refer to










Speaker: 卢朓 Lu Tiao [北京大学 Peiking University]

Time:  3:15pm, Nov. 21 (周二下午3:15)
Location: Conference Room 3, Jin Chun Yuan West Bldg. (近春园西楼二楼第三会议室)

Title: Well-posedness of the initial value problem for a class of time-dependent partial integro-differential equations

Abstract: The abstract initial value problem for the system of evolution equations for a real-valued function and a function-valued function was considered. The existence and uniqueness of classical solution to the evolution system were proved in a Banach space under assumptions on the boundedness and smoothness of data. An isomorphism between the solution
space of the evolution system and a special subspace of the Schwartz space is established. It is verified that the solution of the evolution system can be mapped to a function which is the solution of the initial value problem of an integral-differential equation. The theoretical finding has a potential application in studying the well-posedness of the stationary Wigner equation with inflow boundary conditions.


Speaker: Sijian Wang [Rutgers University]

Time:  3:15pm, Dec. 26
Location: Conference Room 3, Jin Chun Yuan West Bldg. 近春园西楼第三会议室

Title: Statistical Views on Precision Medicine in Diabetes

For details, please see

Abstract:  Diabetes is one of the most prevalent and costly chronic medical conditions worldwide, incurring significant burdens on
individuals, society, and the health care system. Currently, there is a lack of clear clinical guideline on treatment and care for each individual diabetes patient. This is mainly because diabetes patients are heterogeneous and response treatments differently. The usage of specific characteristics of a patient can be helpful for directing diagnostic or treatment strategies that are most effective for that individual. The spectrum of information that can guide personalized
decisions on diabetes treatment and care includes individual behavioral and clinical features, standard clinical laboratory findings and molecular markers. In this talk, we present some statistical methods and consideration of individualized treatments for diabetes patients in the following four problems: individualized therapy recommendation,
personalized treatment initiation, cost-constraint decision rule construction and dynamic risk (cardiovascular event) prediction. Some analysis results using data from Eli Lilly and Company diabetes research group will be presented.