Researcher profile

Hai-Wei Sun

Hai-Wei Sun contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2021arXiv

A preconditioner based on sine transform for two-dimensional Riesz space factional diffusion equations in convex domains

In this paper, we develop a fast numerical method for solving the time-dependent Riesz space fractional diffusion equations with a nonlinear source term in the convex domain. An implicit finite difference method is employed to discretize the Riesz space fractional diffusion equations with a penalty term in a rectangular region by the volume-penalization approach. The stability and the convergence of the proposed method are studied. As the coefficient matrix is with the Toeplitz-like structure, the generalized minimum residual method with a preconditioner based on the sine transform is exploited to solve the discretized linear system, where the preconditioner is constructed in view of the combination of two approximate inverse $τ$ matrices, which can be diagonalized by the sine transform. The spectrum of the preconditioned matrix is also investigated. Numerical experiments are carried out to demonstrate the efficiency of the proposed method.

preprint2021arXiv

Exponential-sum-approximation technique for variable-order time-fractional diffusion equations

In this paper, we study the variable-order (VO) time-fractional diffusion equations. For a VO function $α(t)\in(0,1)$, we develop an exponential-sum-approximation (ESA) technique to approach the VO Caputo fractional derivative. The ESA technique keeps both the quadrature exponents and the number of exponentials in the summation unchanged at the different time levels. Approximating parameters are properly selected to achieve efficient accuracy. Compared with the general direct method, the proposed method reduces the storage requirement from $\mathcal{O}(n)$ to $\mathcal{O}(\log^2 n)$ and the computational cost from $\mathcal{O}(n^2)$ to $\mathcal{O}(n\log^2 n)$, respectively, with $n$ being the number of the time levels. When this fast algorithm is exploited to construct a fast ESA scheme for the VO time-fractional diffusion equations, the computational complexity of the proposed scheme is only of $\mathcal{O}(mn\log^2 n)$ with $\mathcal{O}(m\log^2n)$ storage requirement, where $m$ denotes the number of spatial grids. Theoretically, the unconditional stability and error analysis of the fast ESA scheme are given. The effectiveness of the proposed algorithm is verified by numerical examples.

preprint2021arXiv

Fast second-order evaluation for variable-order Caputo fractional derivative with applications to fractional sub-diffusion equations

In this paper, we propose a fast second-order approximation to the variable-order (VO) Caputo fractional derivative, which is developed based on $L2$-$1_σ$ formula and the exponential-sum-approximation technique. The fast evaluation method can achieve the second-order accuracy and further reduce the computational cost and the acting memory for the VO Caputo fractional derivative. This fast algorithm is applied to construct a relevant fast temporal second-order and spatial fourth-order scheme ($FL2$-$1_σ$ scheme) for the multi-dimensional VO time-fractional sub-diffusion equations. Theoretically, $FL2$-$1_σ$ scheme is proved to fulfill the similar properties of the coefficients as those of the well-studied $L2$-$1_σ$ scheme. Therefore, $FL2$-$1_σ$ scheme is strictly proved to be unconditionally stable and convergent. A sharp decrease in the computational cost and the acting memory is shown in the numerical examples to demonstrate the efficiency of the proposed method.

preprint2021arXiv

Robust fast method for variable-order time-fractional diffusion equations without regularity assumptions

In this paper, we develop a robust fast method for mobile-immobile variable-order (VO) time-fractional diffusion equations (tFDEs), superiorly handling the cases of small or vanishing lower bound of the VO function. The valid fast approximation of the VO Caputo fractional derivative is obtained using integration by parts and the exponential-sum-approximation method. Compared with the general direct method, the proposed algorithm ($RF$-$L1$ formula) reduces the acting memory from $\mathcal{O}(n)$ to $\mathcal{O}(\log^2 n)$ and computational cost from $\mathcal{O}(n^2)$ to $\mathcal{O}(n \log^2 n)$, respectively, where $n$ is the number of time levels. Then $RF$-$L1$ formula is applied to construct the fast finite difference scheme for the VO tFDEs, which sharp decreases the memory requirement and computational complexity. The error estimate for the proposed scheme is studied only under some assumptions of the VO function, coefficients, and the source term, but without any regularity assumption of the true solutions. Numerical experiments are presented to verify the effectiveness of the proposed method.

preprint2021arXiv

Spectral analysis for preconditioning of multi-dimensional Riesz fractional diffusion equations

In this paper, we analyze the spectra of the preconditioned matrices arising from discretized multi-dimensional Riesz spatial fractional diffusion equations. The finite difference method is employed to approximate the multi-dimensional Riesz fractional derivatives, which will generate symmetric positive definite ill-conditioned multi-level Toeplitz matrices. The preconditioned conjugate gradient method with a preconditioner based on the sine transform is employed to solve the resulting linear system. Theoretically, we prove that the spectra of the preconditioned matrices are uniformly bounded in the open interval (1/2,3/2) and thus the preconditioned conjugate gradient method converges linearly. The proposed method can be extended to multi-level Toeplitz matrices generated by functions with zeros of fractional order. Our theoretical results fill in a vacancy in the literature. Numerical examples are presented to demonstrate our new theoretical results in the literature and show the convergence performance of the proposed preconditioner that is better than other existing preconditioners.