BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Zhonghua Qiao

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

ResearcherUnclaimed source record
math.NANumerical Analysismath.PRComputer Visionmath-phmath.APmath.MP

Catalog footprint

What is connected

16works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

16 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Acceleration for Microflow Simulations of High-Order Moment Models by Using Lower-Order Model Correction

We study the acceleration of steady-state computation for microflow, which is modeled by the high-order moment models derived recently from the steady-state Boltzmann equation with BGK-type collision term. By using the lower-order model correction, a novel nonlinear multi-level moment solver is developed. Numerical examples verify that the resulting solver improves the convergence significantly thus is able to accelerate the steady-state computation greatly. The behavior of the solver is also numerically investigated. It is shown that the convergence rate increases, indicating the solver would be more efficient, as the total levels increases. Three order reduction strategies of the solver are considered. Numerical results show that the most efficient order reduction strategy would be $m_{l-1} = \lceil m_{l} / 2 \rceil$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Characterizing the stabilization size for semi-implicit Fourier-spectral method to phase field equations

Recent results in the literature provide computational evidence that stabilized semi-implicit time-stepping method can efficiently simulate phase field problems involving fourth-order nonlinear dif- fusion, with typical examples like the Cahn-Hilliard equation and the thin film type equation. The up-to-date theoretical explanation of the numerical stability relies on the assumption that the deriva- tive of the nonlinear potential function satisfies a Lipschitz type condition, which in a rigorous sense, implies the boundedness of the numerical solution. In this work we remove the Lipschitz assumption on the nonlinearity and prove unconditional energy stability for the stabilized semi-implicit time-stepping methods. It is shown that the size of stabilization term depends on the initial energy and the perturba- tion parameter but is independent of the time step. The corresponding error analysis is also established under minimal nonlinearity and regularity assumptions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gradient bounds for a thin film epitaxy equation

We consider a gradient flow modeling the epitaxial growth of thin films with slope selection. The surface height profile satisfies a nonlinear diffusion equation with biharmonic dissipation. We establish optimal local and global wellposedness for initial data with critical regularity. To understand the mechanism of slope selection and the dependence on the dissipation coefficient, we exhibit several lower and upper bounds for the gradient of the solution in physical dimensions $d\le 3$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An unconditionally energy stable finite difference scheme for a stochastic Cahn-Hilliard equation

In this work, the MMC-TDGL equation, a stochastic Cahn-Hilliard equation is solved numerically by using the finite difference method in combination with a convex splitting technique of the energy functional. For the non-stochastic case, we develop an unconditionally energy stable difference scheme which is proved to be uniquely solvable. For the stochastic case, by adopting the same splitting of the energy functional, we construct a similar and uniquely solvable difference scheme with the discretized stochastic term. The resulted schemes are nonlinear and solved by Newton iteration. For the long time simulation, an adaptive time stepping strategy is developed based on both first- and second-order derivatives of the energy. Numerical experiments are carried out to verify the energy stability, the efficiency of the adaptive time stepping and the effect of the stochastic term.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Convergence of a Fast Explicit Operator Splitting Method for the Molecular Beam Epitaxy Model

A fast explicit operator splitting (FEOS) method for the molecular beam epitaxy model has been presented in [Cheng, et al., Fast and stable explicit operator splitting methods for phase-field models, J. Comput. Phys., submitted]. The original problem is split into linear and nonlinear subproblems. For the linear part, the pseudo-spectral method is adopted; for the nonlinear part, a 33-point difference scheme is constructed. Here, we give a compact center-difference scheme involving fewer points for the nonlinear subproblem. Besides, we analyze the convergence rate of the algorithm. The global error order $\mathcal{O}(τ^2+h^4)$ in discrete $L^2$-norm is proved theoretically and verified numerically. Some numerical experiments show the robustness of the algorithm for small coefficients of the fourth-order term for the one-dimensional case. Besides, coarsening dynamics are simulated in large domains and the $1/3$ power laws are observed for the two-dimensional case.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

NRxx Simulation of Microflows with Shakhov Model

In this paper, we propose a method to simulate the microflows with Shakhov model using the NRxx method developed in [4, 5, 6]. The equation under consideration is the Boltzmann equation with force terms and the Shakhov model is adopted to achieve the correct Prandtl number. As the focus of this paper, we derive a uniform framework for different order moment systems on the wall boundary conditions, which is a major difficulty in the moment methods. Numerical examples for both steady and unsteady problems are presented to show the convergence in the number of moments.

0 citations0 reviewsNo code linkedOpen access