A novel evolving surface finite element method, based on a novel equivalent formulation of the continuous problem, is proposed for computing the evolution of a closed hypersurface moving under a prescribed velocity field in two- and threedimensional spaces. The method improves the mesh quality of the approximate surface by minimizing the rate of deformation using an artificial tangential motion. The transport evolution equations of the normal vector and the extrinsic Weingarten matrix are derived and coupled with the surface evolution equations to ensure stability and convergence of the numerical approximations. Optimal-order convergence of the semi-discrete evolving surface finite element method is proved for finite elements of degree k ≥ 2. Numerical examples are provided to illustrate the convergence of the proposed method and its effectiveness in improving mesh quality on the approximate evolving surface.
@article{BHL2024SJNA,title={A convergent evolving finite element method with artificial tangential motion for surface evolution under a prescribed velocity field},author={Bai, Genming and Hu, Jiashun and Li, Buyang},date={2024},volume={62},number={5},pages={2172--2195},doi={10.1137/23M156968X},langid={english},journal={SIAM Journal on Numerical Analysis},year={2024},dimensions={true},}
A class of high-order mass- and energy-conserving methods is proposed for the nonlinear Schr\o"dinger equation based on Gauss collocation in time and finite element discretization in space, by introducing a mass- and energy-correction post-process at every time level. The existence, uniqueness, and high-order convergence of the numerical solutions are proved. In particular, the error of the numerical solution is O(τk+1 + hp) in the L∞(0, T ; H1) norm after incorporating the accumulation errors arising from the post-processing correction procedure at all time levels, where k and p denote the degrees of finite elements in time and space, respectively, which can be arbitrarily large. Several numerical examples are provided to illustrate the performance of the proposed new method, including the conservation of mass and energy and the high-order convergence in simulating solitons and bi-solitons.
@article{BHL2024SJSC,title={High-order mass- and energy-conserving methods for the nonlinear {{Schr\"odinger}} equation},author={Bai, Genming and Hu, Jiashun and Li, Buyang},date={2024-04-30},volume={46},number={2},pages={A1026-A1046},doi={10.1137/22M152178X},langid={english},journal={SIAM Journal on Scientific Compututing},year={2024},dimensions={true},}
An artificial tangential velocity is introduced into the evolving finite element methods for mean curvature flow and Willmore flow proposed by Kovács et al. (Numer Math 143(4), 797-853, 2019, Numer Math 149, 595-643, 2021) in order to improve the mesh quality in the computation. The artificial tangential velocity is constructed by considering a limiting situation in the method proposed by Barrett et al. (J Comput Phys 222(1), 441-467, 2007, J Comput Phys 227(9), 4281-4307, 2008, SIAM J Sci Comput 31(1), 225-253, 2008) . The stability of the artificial tangential velocity is proved. The optimal-order convergence of the evolving finite element methods with artificial tangential velocity are proved for both mean curvature flow and Willmore flow. Extensive numerical experiments are presented to illustrate the convergence of the method and the performance of the artificial tangential velocity in improving the mesh quality.
@article{HL2022NM,title={Evolving finite element methods with an artificial tangential velocity for mean curvature flow and {{Willmore}} flow},author={Hu, Jiashun and Li, Buyang},date={2022},volume={152},number={1},pages={127--181},doi={10.1007/s00211-022-01309-9},langid={english},journal={Numerische Mathematik},year={2022},dimensions={true},}
For the Cauchy problem of the high frequency wave-type equation with Wentzel–Kramers–Brillouin (WKB) type initial data, the extended Wentzel–Kramers–Brillouin (E-WKB) ansatz is an asymptotically valid solution. This ansatz is globally defined and formulated as an integral of coherent states over the displaced Lagrangian submanifold. This paper proves the optimal first order error estimate of the proposed E-WKB ansatz in \Ł^2\ L 2 norm for the wave-type equation in the semi-classical regime. The key ingredients in the proof are the moving frame technique developed in Zheng (Commun Math Sci 11:105–140, 2013) and the deep relations between the E-WKB analysis and the classical WKB analysis. Numerical results on the linear KdV equation verify the theoretical analysis.
2018
Front. Math. China
Global geometrical optics method for vector-valued Schrödinger problems
@article{HMZ2018FMC,title={Global geometrical optics method for vector-valued {{Schr\"odinger}} problems},author={Hu, Jiashun and Ma, Xiang and Zheng, Chunxiong},date={2018-06},volume={13},number={3},pages={579--606},doi={10.1007/s11464-018-0704-1},eissn={1673-3576},journal={Frontiers of Mathematics in China},year={2018}}
2017
J. Comput. Appl. Math.
Fast and stable evaluation of the exact absorbing boundary condition for the semi-discrete linear Schrödinger equation in unbounded domains
@article{HZ2017JCAM,title={Fast and stable evaluation of the exact absorbing boundary condition for the semi-discrete linear {{Schr\"odinger}} equation in unbounded domains},author={Hu, Jiashun and Zheng, Chunxiong},date={2017-12-15},volume={326},pages={116--125},doi={10.1016/j.cam.2017.05.018},eissn={1879-1778},journal={Journal of Computational and Applied Mathematics},year={2017}}
Numerical computation of a nonlocal diffusion equation on the real axis is considered in this paper. We first apply an extensively studied quadrature scheme to obtain a discrete nonlocal diffusion system on an unbounded domain. Then we derive an alternative formulation of the discrete problem based on the spectral analysis of the \z\-transform. This new formulation can be seen as a system defined on a bounded domain with an artificial boundary condition, and it allows us to reformulate the original infinite domain problem into an equivalent bounded domain problem. To numerically implement the exact artificial boundary condition, we apply the trapezoidal quadrature rule to approximate the contour integral induced by the inverse \z\-transform. Numerical examples are presented to demonstrate the effectiveness of our approach.
@article{ZHD+2017SJSC,title={Numerical solution of the nonlocal diffusion equation on the real line},author={Zheng, Chunxiong and Hu, Jiashun and Du, Qiang and Zhang, Jiwei},date={2017-01-01},volume={39},number={5},pages={A1951-A1968},publisher={{Society for Industrial and Applied Mathematics}},doi={10.1137/16M1090107},journal={SIAM Journal on Scientific Computing},year={2017},dimensions={true},}
