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

Huangxin Chen

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

ResearcherUnclaimed source record
math.NANumerical AnalysisComputational Engineering, Finance, and Sciencephysics.comp-ph

Catalog footprint

What is connected

13works
4topics
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

13 published item(s)

preprint2024arXiv

A prediction-correction based iterative convolution-thresholding method for topology optimization of heat transfer problems

In this paper, we propose an iterative convolution-thresholding method (ICTM) based on prediction-correction for solving the topology optimization problem in steady-state heat transfer equations. The problem is formulated as a constrained minimization problem of the complementary energy, incorporating a perimeter/surface-area regularization term, while satisfying a steady-state heat transfer equation. The decision variables of the optimization problem represent the domains of different materials and are represented by indicator functions. The perimeter/surface-area term of the domain is approximated using Gaussian kernel convolution with indicator functions. In each iteration, the indicator function is updated using a prediction-correction approach. The prediction step is based on the variation of the objective functional by imposing the constraints, while the correction step ensures the monotonically decreasing behavior of the objective functional. Numerical results demonstrate the efficiency and robustness of our proposed method, particularly when compared to classical approaches based on the ICTM.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A $C^{0}$ interior penalty method for $m$th-Laplace equation

In this paper, we propose a $C^{0}$ interior penalty method for $m$th-Laplace equation on bounded Lipschitz polyhedral domain in $\mathbb{R}^{d}$, where $m$ and $d$ can be any positive integers. The standard $H^{1}$-conforming piecewise $r$-th order polynomial space is used to approximate the exact solution $u$, where $r$ can be any integer greater than or equal to $m$. Unlike the interior penalty method in [T.~Gudi and M.~Neilan, {\em An interior penalty method for a sixth-order elliptic equation}, IMA J. Numer. Anal., \textbf{31(4)} (2011), pp. 1734--1753], we avoid computing $D^{m}$ of numerical solution on each element and high order normal derivatives of numerical solution along mesh interfaces. Therefore our method can be easily implemented. After proving discrete $H^{m}$-norm bounded by the natural energy semi-norm associated with our method, we manage to obtain stability and optimal convergence with respect to discrete $H^{m}$-norm. Numerical experiments validate our theoretical estimate.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A fully discrete energy stable scheme for a phase-field moving contact line model with variable densities and viscosities

In this work, we propose a fully discrete energy stable scheme for the phase-field moving contact line model with variable densities and viscosities. The mathematical model consists of a Cahn-Hilliard equation, a Navier-Stokes equation and the generalized Navier boundary condition for the moving contact line. A scalar auxiliary variable is adopted to transform the governing system into an equivalent form, allowing the double well potential to be treated semi-explicitly. A stabilization term is added to balance the explicit nonlinear term originating from the surface energy at fluid-solid interface. A pressure stabilization method is used to decouple the computation of velocity and pressure. Some subtle implicit-explicit treatments are adopted to deal with convention and stress terms. We establish a rigorous proof of energy stability for the proposed time-marching scheme. Then a finite difference method on staggered grids is used to spatially discretize the constructed time-marching scheme. We further prove that the fully discrete scheme also satisfies the discrete energy dissipation law. Numerical results demonstrate accuracy and energy stability of the proposed scheme. Using our numerical scheme, we analyze the contact line dynamics through a shear flow driven droplet sliding case. Three-dimensional droplet spreading is also investigated on a chemically patterned surface. Our numerical simulation accurately predicts the expected energy evolutions and it successfully reproduces expected phenomena that an oil droplet contracts inwards on a hydrophobic zone and spreads outwards quickly on a hydrophilic zone.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A mixed finite element scheme for biharmonic equation with variable coefficient and von Kármán equations

In this paper, a new mixed finite element scheme using element-wise stabilization is introduced for the biharmonic equation with variable coefficient on Lipschitz polyhedral domains. The proposed scheme doesn't involve any integration along mesh interfaces. The gradient of the solution is approximated by $H({\rm div})$-conforming $BDM_{k+1}$ element or vector valued Lagrange element with order $k+1$, while the solution is approximated by Lagrange element with order $k+2$ for any $k\geq 0$.This scheme can be easily implemented and produces positive definite linear system. We provide a new discrete $H^{2}$-norm stability, which is useful not only in analysis of this scheme but also in $C^{0}$ interior penalty methods and DG methods. Optimal convergences in both discrete $H^{2}$-norm and $L^{2}$-norm are derived. This scheme with its analysis is further generalized to the von Kármán equations. Finally, numerical results verifying the theoretical estimates of the proposed algorithms are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A new physics-preserving IMPES scheme for incompressible and immiscible two-phase flow in heterogeneous porous media

In this work we consider a new efficient IMplicit Pressure Explicit Saturation (IMPES) scheme for the simulation of incompressible and immiscible two-phase flow in heterogeneous porous media with capillary pressure. Compared with the conventional IMPES schemes, the new IMPES scheme is inherently physics-preserving, namely, the new algorithm is locally mass conservative for both phases and it also enjoys another appealing feature that the total velocity is continuous in the normal direction. Moreover, the new scheme is unbiased with regard to the two phases and the saturations of both phases are bounds-preserving if the time step size is smaller than a certain value. The key ideas in the new scheme include that the Darcy flows for both phases are rewritten in the formulation based on the total velocity and an auxiliary velocity referring to as the capillary potential gradient, and the total discretized conservation equation is obtained by the summation of the discretized conservation equation for each phase. The upwind strategy is applied to update the saturations explicitly, and the upwind mixed finite element methods are used to solve the pressure-velocity systems which can be decoupled. We also present some interesting examples to demonstrate the efficiency and robustness of the new algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A superconvergent HDG method for the Maxwell equations

We present and analyze a new hybridizable discontinuous Galerkin (HDG) method for the steady state Maxwell equations. In order to make the problem well-posed, a condition of divergence is imposed on the electric field. Then a Lagrange multiplier $p$ is introduced, and the problem becomes the solution of a mixed curl-curl formulation of the Maxwell's problem. We use polynomials of degree $k+1$, $k$, $k$ to approximate $\bfu,\nabla \times \bfu$ and $p$ respectively. In contrast, we only use a non-trivial subspace of polynomials of degree $k+1$ to approximate the numerical tangential trace of the electric field and polynomials of degree $k+1$ to approximate the numerical trace of the Lagrange multiplier on the faces. On the simplicial meshes, a special choice of the stabilization parameters is applied, and the HDG system is shown to be well-posed. Moreover, we show that the convergence rates for $\boldsymbol{u}$ and $\nabla \times \boldsymbol{u}$ are independent of the Lagrange multiplier $p$. If we assume the dual operator of the Maxwell equation on the domain has adequate regularity, we show that the convergence rate for $\boldsymbol{u}$ is $O(h^{k+2})$. From the point of view of degrees of freedom of the globally coupled unknown: numerical trace, this HDG method achieves superconvergence for the electric field without postprocessing. Finally, we show that on general polyhedral elements, by a particular choice of the stabilization parameters again, the HDG system is also well-posed and the superconvergence of the HDG method is derived.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An absolutely stable $hp$-HDG method for the time-harmonic Maxwell equations with high wave number

We present and analyze a hybridizable discontinuous Galerkin (HDG) method for the time-harmonic Maxwell equations. The divergence-free condition is enforced on the electric field, then a Lagrange multiplier is introduced, and the problem becomes the solution of a mixed curl-curl formulation of the Maxwell's problem. The method is shown to be an absolutely stable HDG method for the indefinite time-harmonic Maxwell equations with high wave number. By exploiting the duality argument, the dependence of convergence of the HDG method on the wave number k, the mesh size h and the polynomial order p is obtained. Numerical results are given to verify the theoretical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A first order system least squares method for the Helmholtz equation

We present a first order system least squares (FOSLS) method for the Helmholtz equation at high wave number k, which always deduces Hermitian positive definite algebraic system. By utilizing a non-trivial solution decomposition to the dual FOSLS problem which is quite different from that of standard finite element method, we give error analysis to the hp-version of the FOSLS method where the dependence on the mesh size h, the approximation order p, and the wave number k is given explicitly. In particular, under some assumption of the boundary of the domain, the L2 norm error estimate of the scalar solution from the FOSLS method is shown to be quasi optimal under the condition that kh/p is sufficiently small and the polynomial degree p is at least O(\log k). Numerical experiments are given to verify the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

First order least squares method with weakly imposed boundary condition for convection dominated diffusion problems

We present and analyze a first order least squares method for convection dominated diffusion problems, which provides robust L2 a priori error estimate for the scalar variable even if the given data f in L2 space. The novel theoretical approach is to rewrite the method in the framework of discontinuous Petrov - Galerkin (DPG) method, and then show numerical stability by using a key equation discovered by J. Gopalakrishnan and W. Qiu [Math. Comp. 83(2014), pp. 537-552]. This new approach gives an alternative way to do numerical analysis for least squares methods for a large class of differential equations. We also show that the condition number of the global matrix is independent of the diffusion coefficient. A key feature of the method is that there is no stabilization parameter chosen empirically. In addition, Dirichlet boundary condition is weakly imposed. Numerical experiments verify our theoretical results and, in particular, show our way of weakly imposing Dirichlet boundary condition is essential to the design of least squares methods - numerical solutions on subdomains away from interior layers or boundary layers have remarkable accuracy even on coarse meshes, which are unstructured quasi-uniform.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Robust a Posteriori Error Estimates for HDG method for Convection-Diffusion Equations

We propose a robust a posteriori error estimator for the hybridizable discontinuous Galerkin (HDG) method for convection-diffusion equations with dominant convection. The reliability and efficiency of the estimator are established for the error measured in an energy norm. The energy norm is uniformly bounded even when the diffusion coefficient tends to zero. The estimators are robust in the sense that the upper and lower bounds of error are uniformly bounded with respect to the diffusion coefficient. A weighted test function technique and the Oswald interpolation are key ingredients in the analysis. Numerical results verify the robustness of the proposed a posteriori error estimator. In numerical experiments, optimal convergence is observed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Robust Multilevel Method for Hybridizable Discontinuous Galerkin Method for the Helmholtz Equation

A robust multilevel preconditioner based on the hybridizable discontinuous Galerkin method for the Helmholtz equation with high wave number is presented in this paper. There are two keys in our algorithm, one is how to choose a suitable intergrid transfer operator, and the other is using GMRES smoothing on coarse grids. The multilevel method is performed as a preconditioner in the outer GMRES iteration. To give a quantitative insight of our algorithm, we use local Fourier analysis to analyze the convergence property of the proposed multilevel method. Numerical results show that for fixed wave number, the convergence of the algorithm is mesh independent. Moreover, the performance of the algorithm depends relatively mildly on wave number.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Multilevel Preconditioner with Stable Coarse Grid Corrections for the Helmholtz Equation

In this paper we consider a class of robust multilevel precontioners for the Helmholtz equation with high wave number. The key idea in this work is to use the continuous interior penalty finite element methods (CIP-FEM) studied in \cite{Wu12,Wu12-hp} to construct the stable coarse grid correction problems. The multilevel methods, based on GMRES smoothing on coarse grids, are then served as a preconditioner in the outer GMRES iteration. In the one dimensional case, convergence property of the modified multilevel methods is analyzed by the local Fourier analysis. From our numerical results, we find that the proposed methods are efficient for a reasonable range of frequencies. The performance of the algorithms depends relatively mildly on wave number. In particular, only one GMRES smoothing step may guarantee the optimal convergence of our multilevel algorithm, which remedies the shortcoming of the multilevel algorithm in \cite{EEO01}.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Hybridizable Discontinuous Galerkin Method for the Helmholtz Equation with High Wave Number

This paper analyzes the error estimates of the hybridizable discontinuous Galerkin (HDG) method for the Helmholtz equation with high wave number in two and three dimensions. The approximation piecewise polynomial spaces we deal with are of order $p\geq 1$. Through choosing a specific parameter and using the duality argument, it is proved that the HDG method is stable without any mesh constraint for any wave number $κ$. By exploiting the stability estimates, the dependence of convergence of the HDG method on $κ,h$ and $p$ is obtained. Numerical experiments are given to verify the theoretical results.

0 citations0 reviewsNo code linkedOpen access