Researcher profile

Zhonghua Qiao

Zhonghua Qiao contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
4topics
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

10 published item(s)

preprint2022arXiv

An Active Contour Model with Local Variance Force Term and Its Efficient Minimization Solver for Multi-phase Image Segmentation

In this paper, we propose an active contour model with a local variance force (LVF) term that can be applied to multi-phase image segmentation problems. With the LVF, the proposed model is very effective in the segmentation of images with noise. To solve this model efficiently, we represent the regularization term by characteristic functions and then design a minimization algorithm based on a modification of the iterative convolution-thresholding method (ICTM), namely ICTM-LVF. This minimization algorithm enjoys the energy-decaying property under some conditions and has highly efficient performance in the segmentation. To overcome the initialization issue of active contour models, we generalize the inhomogeneous graph Laplacian initialization method (IGLIM) to the multi-phase case and then apply it to give the initial contour of the ICTM-LVF solver. Numerical experiments are conducted on synthetic images and real images to demonstrate the capability of our initialization method, and the effectiveness of the local variance force for noise robustness in the multi-phase image segmentation.

preprint2022arXiv

An implicit--explicit second order BDF numerical scheme with variable steps for gradient flows

In this paper, we propose and analyze an efficient implicit--explicit (IMEX) second order in time backward differentiation formulation (BDF2) scheme with variable time steps for gradient flow problems using the scalar auxiliary variable (SAV) approach. We prove the unconditional energy stability of the scheme for a modified discrete energy with the adjacent time step ratio $γ_{n+1}:=\Dt_{n+1}/\Dt_{n}\leq 4.8645$. The uniform $H^{2}$ bound for the numerical solution is derived under a mild regularity restriction on the initial condition, that is $ϕ(\x,0)\in H^{2}$. Based on this uniform bound, a rigorous error estimate of the numerical solution is carried out on the temporal nonuniform mesh. Finally, serval numerical tests are provided to validate the theoretical claims. With the application of an adaptive time-stepping strategy, the efficiency of our proposed scheme can be clearly observed in the coarsening dynamics simulation.

preprint2022arXiv

Double stabilizations and convergence analysis of a second-order linear numerical scheme for the nonlocal Cahn-Hilliard equation

In this paper, we study a second-order accurate and linear numerical scheme for the nonlocal Cahn-Hilliard equation. The scheme is established by combining a modified Crank-Nicolson approximation and the Adams-Bashforth extrapolation for the temporal discretization, and by applying the Fourier spectral collocation to the spatial discretization. In addition, two stabilization terms in different forms are added for the sake of the numerical stability. We conduct a complete convergence analysis by using the higher-order consistency estimate for the numerical scheme, combined with the rough error estimate and the refined estimate. By regarding the numerical solution as a small perturbation of the exact solution, we are able to justify the discrete $\ell^\infty$ bound of the numerical solution, as a result of the rough error estimate. Subsequently, the refined error estimate is derived to obtain the optimal rate of convergence, following the established $\ell^\infty$ bound of the numerical solution. Moreover, the energy stability is also rigorously proved with respect to a modified energy. The proposed scheme can be viewed as the generalization of the second-order scheme presented in an earlier work, and the energy stability estimate has greatly improved the corresponding result therein.

preprint2022arXiv

Generalized SAV-exponential integrator schemes for Allen-Cahn type gradient flows

The energy dissipation law and the maximum bound principle (MBP) are two important physical features of the well-known Allen-Cahn equation. While some commonly-used first-order time stepping schemes have turned out to preserve unconditionally both energy dissipation law and MBP for the equation, restrictions on the time step size are still needed for existing second-order or even higher-order schemes in order to have such simultaneous preservation. In this paper, we develop and analyze novel first- and second-order linear numerical schemes for a class of Allen-Cahn type gradient flows. Our schemes combine the generalized scalar auxiliary variable (SAV) approach and the exponential time integrator with a stabilization term, while the standard central difference stencil is used for discretization of the spatial differential operator. We not only prove their unconditional preservation of the energy dissipation law and the MBP in the discrete setting, but also derive their optimal temporal error estimates under fixed spatial mesh. Numerical experiments are also carried out to demonstrate the properties and performance of the proposed schemes.

preprint2022arXiv

Multi-phase image segmentation by the Allen--Cahn Chan--Vese model

This paper proposes an Allen-Cahn Chan-Vese model to settle the multi-phase image segmentation. We first integrate the Allen--Cahn term and the Chan--Vese fitting energy term to establish an energy functional, whose minimum locates the segmentation contour. The subsequent minimization process can be attributed to variational calculation on fitting intensities and the solution approximation of several Allen-Cahn equations, wherein $n$ Allen-Cahn equations are enough to partition $m = 2^n$ segments. The derived Allen-Cahn equations are solved by efficient numerical solvers with exponential time integrations and finite difference space discretization. The discrete maximum bound principle and energy stability of the proposed numerical schemes are proved. Finally, the capability of our segmentation method is verified in various experiments for different types of images.

preprint2022arXiv

Stabilized exponential-SAV schemes preserving energy dissipation law and maximum bound principle for the Allen-Cahn type equations

It is well-known that the Allen-Cahn equation not only satisfies the energy dissipation law but also possesses the maximum bound principle (MBP) in the sense that the absolute value of its solution is pointwise bounded for all time by some specific constant under appropriate initial/boundary conditions. In recent years, the scalar auxiliary variable (SAV) method and many of its variants have attracted much attention in numerical solution for gradient flow problems due to their inherent advantage of preserving certain discrete analogues of the energy dissipation law. However, existing SAV schemes usually fail to preserve the MBP when applied to the Allen-Cahn equation. In this paper, we develop and analyze new first- and second-order stabilized exponential-SAV schemes for a class of Allen-Cahn type equations, which are shown to simultaneously preserve the energy dissipation law and MBP in discrete settings. In addition, optimal error estimates for the numerical solutions are rigorously obtained for both schemes. Extensive numerical tests and comparisons are also conducted to demonstrate the performance of the proposed schemes.

preprint2022arXiv

Strong Approximation of Monotone Stochastic Partial Differential Equations Driven by Multiplicative Noise

We establish a general theory of optimal strong error estimation for numerical approximations of a second-order parabolic stochastic partial differential equation with monotone drift driven by a multiplicative infinite-dimensional Wiener process. The equation is spatially discretized by Galerkin methods and temporally discretized by drift-implicit Euler and Milstein schemes. By the monotone and Lyapunov assumptions, we use both the variational and semigroup approaches to derive a spatial Sobolev regularity under the $L_ω^p L_t^\infty \dot H^{1+γ}$-norm and a temporal Hölder regularity under the $L_ω^p L_x^2$-norm for the solution of the proposed equation with an $\dot H^{1+γ}$-valued initial datum for $γ\in [0,1]$. Then we make full use of the monotonicity of the equation and tools from stochastic calculus to derive the sharp strong convergence rates $O(h^{1+γ}+τ^{1/2})$ and $O(h^{1+γ}+τ^{(1+γ)/2})$ for the Galerkin-based Euler and Milstein schemes, respectively.

preprint2020arXiv

Convergence analysis for a stabilized linear semi-implicit numerical scheme for the nonlocal Cahn-Hilliard equation

In this paper, we provide a detailed convergence analysis for a first order stabilized linear semi-implicit numerical scheme for the nonlocal Cahn-Hilliard equation, which follows from consistency and stability estimates for the numerical error function. Due to the complicated form of the nonlinear term, we adopt the discrete $H^{-1}$ norm for the error function to establish the convergence result. In addition, the energy stability obtained in [Du et al., J. Comput. Phys., 363:39--54, 2018] requires an assumption on the uniform $\ell^\infty$ bound of the numerical solution and such a bound is figured out in this paper by conducting the higher order consistency analysis. Taking the view that the numerical solution is indeed the exact solution with a perturbation, the error function is $\ell^\infty$ bounded uniformly under a loose constraint of the time step size, which then leads to the uniform maximum-norm bound of the numerical solution.

preprint2020arXiv

Maximum bound principles for a class of semilinear parabolic equations and exponential time differencing schemes

The ubiquity of semilinear parabolic equations has been illustrated in their numerous applications ranging from physics, biology, to materials and social sciences. In this paper, we consider a practically desirable property for a class of semilinear parabolic equations of the abstract form $u_t=\mathcal{L}u+f[u]$ with $\mathcal{L}$ being a linear dissipative operator and $f$ being a nonlinear operator in space, namely a time-invariant maximum bound principle, in the sense that the time-dependent solution $u$ preserves for all time a uniform pointwise bound in absolute value imposed by its initial and boundary conditions. We first study an analytical framework for some sufficient conditions on $\mathcal{L}$ and $f$ that lead to such a maximum bound principle for the time-continuous dynamic system of infinite or finite dimensions. Then, we utilize a suitable exponential time differencing approach with a properly chosen generator of contraction semigroup to develop first- and second-order accurate temporal discretization schemes, that satisfy the maximum bound principle unconditionally in the time-discrete setting. Error estimates of the proposed schemes are derived along with their energy stability. Extensions to vector- and matrix-valued systems are also discussed. We demonstrate that the abstract framework and analysis techniques developed here offer an effective and unified approach to study the maximum bound principle of the abstract evolution equation that cover a wide variety of well-known models and their numerical discretization schemes. Some numerical experiments are also carried out to verify the theoretical results.

preprint2020arXiv

Strong Approximation of Stochastic Allen-Cahn Equation with White Noise

We establish an optimal strong convergence rate of a fully discrete numerical scheme for second order parabolic stochastic partial differential equations with monotone drifts, including the stochastic Allen-Cahn equation, driven by an additive space-time white noise. Our first step is to transform the original stochastic equation into an equivalent random equation whose solution possesses more regularity than the original one. Then we use the backward Euler in time and spectral Galerkin in space to fully discretize this random equation. By the monotone assumption, in combination with the factorization method and stochastic calculus in martingale-type 2 Banach spaces, we derive a uniform maximum norm estimation and a Hölder-type regularity for both stochastic and random equations. Finally, the strong convergence rate of the proposed fully discrete scheme under the $l_t^\infty L^2_ωL^2_x \cap l_t^q L^q_ωL^q_x$-norm is obtained. Several numerical experiments are carried out to verify the theoretical result.