Source author record

Wenqiang Feng

Wenqiang Feng appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

4works
1topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

4 published item(s)

preprint2016arXiv

A Second Order Energy Stable Scheme for the Cahn-Hilliard-Hele-Shaw Equations

We present a second-order-in-time finite difference scheme for the Cahn-Hilliard-Hele-Shaw equations. This numerical method is uniquely solvable and unconditionally energy stable. At each time step, this scheme leads to a system of nonlinear equations that can be efficiently solved by a nonlinear multigrid solver. Owing to the energy stability, we derive an $\ell^2 (0,T; H_h^3)$ stability of the numerical scheme. To overcome the difficulty associated with the convection term $\nabla \cdot (ϕ\boldsymbol{u})$, we perform an $\ell^\infty (0,T; H_h^1)$ error estimate instead of the classical $\ell^\infty (0,T; \ell^2)$ one to obtain the optimal rate convergence analysis. In addition, various numerical simulations are carried out, which demonstrate the accuracy and efficiency of the proposed numerical scheme.

preprint2016arXiv

An Energy Stable Finite-Difference Scheme for Functionalized Cahn-Hilliard Equation and its Convergence Analysis

We present and analyze an unconditionally energy stable and convergent finite difference scheme for the Functionalized Cahn-Hilliard equation. One key difficulty associated with the energy stability is based on the fact that one nonlinear energy functional term in the expansion appears as non-convex, non-concave. To overcome this subtle difficulty, we add two auxiliary terms to make the combined term convex, which in turns yields a convex-concave decomposition of the physical energy. As a result, an application of the convex splitting methodology assures both the unique solvability and the unconditional energy stability of the proposed numerical scheme. To deal with a 4-Laplacian solver in an $H^{-1}$ gradient flow at each time step, we apply an efficient preconditioned steepest descent algorithm to solve the corresponding nonlinear systems. In addition, a global in time $H_{\rm per}^2$ stability of the numerical scheme is established at a theoretical level, which in turn ensures the full order convergence analysis of the scheme. A few numerical results are presented, which confirm the stability and accuracy of the proposed numerical scheme.

preprint2016arXiv

Convergence Analysis and Numerical Implementation of a Second Order Numerical Scheme for the Three-Dimensional Phase Field Crystal Equation

In this paper we analyze and implement a second-order-in-time numerical scheme for the three-dimensional phase field crystal (PFC) equation. The numerical scheme was proposed in [46], with the unique solvability and unconditional energy stability established. However, its convergence analysis remains open. We present a detailed convergence analysis in this article, in which the maximum norm estimate of the numerical solution over grid points plays an essential role. Moreover, we outline the detailed multigrid method to solve the highly nonlinear numerical scheme over a cubic domain, and various three-dimensional numerical results are presented, including the numerical convergence test, complexity test of the multigrid solver and the polycrystal growth simulation.

preprint2014arXiv

A Fourier Pseudospectral Method for the "Good" Boussinesq Equation with Second Order Temporal Accuracy

In this paper, we discuss the nonlinear stability and convergence of a fully discrete Fourier pseudospectral method coupled with a specially designed second order time-stepping for the numerical solution of the "good" Boussinesq equation. Our analysis improves the existing results presented in earlier literature in two ways. First, an $l_\infty(0, T^*; H2)$ convergence for the solution and $l_\infty(0, T^*; l_2)$ convergence for the time-derivative of the solution are obtained in this paper, instead of the $l_\infty(0, T^*; l_2)$ convergence for the solution and the $l_\infty(0, T^*; H^{-2})$ convergence for the time-derivative, given in [17]. In addition, the stability and convergence of this method is shown to be unconditional for the time step in terms of the spatial grid size, compared with a severe restriction time step restriction $Δt \leq Ch^2$ reported in [17].