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

Weizhang Huang

Weizhang Huang 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 Sciencemath.AP

Catalog footprint

What is connected

39works
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

39 published item(s)

preprint2022arXiv

A study on CFL conditions for the DG solution of conservation laws on adaptive moving meshes

The selection of time step plays a crucial role in improving stability and efficiency in the Discontinuous Galerkin (DG) solution of hyperbolic conservation laws on adaptive moving meshes that typically employs explicit stepping. A commonly used selection of time step is a direct extension based on Courant-Friedrichs-Levy (CFL) conditions established for fixed and uniform meshes. In this work, we provide a mathematical justification for those time step selection strategies used in practical adaptive DG computations. A stability analysis is presented for a moving mesh DG method for linear scalar conservation laws. Based on the analysis, a new selection strategy of the time step is proposed, which takes into consideration the coupling of the $α$-function (that is related to the eigenvalues of the Jacobian matrix of the flux and the mesh movement velocity) and the heights of the mesh elements. The analysis also suggests several stable combinations of the choices of the $α$-function in the numerical scheme and in the time step selection. Numerical results obtained with a moving mesh DG method for Burgers' and Euler equations are presented. For comparison purpose, numerical results obtained with an error-based time step-size selection strategy are also given.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Metric Tensor Approach to Data Assimilation with Adaptive Moving Meshes

Adaptive moving spatial meshes are useful for solving physical models given by time-dependent partial differentialequations. However, special consideration must be given when combining adaptive meshing procedures with ensemble-based data assimilation (DA) techniques. In particular, we focus on the case where each ensemble member evolvesindependently upon its own mesh and is interpolated to a common mesh for the DA update. This paper outlines aframework to develop time-dependent reference meshes using locations of observations and the metric tensors (MTs)or monitor functions that define the spatial meshes of the ensemble members. We develop a time-dependent spatiallocalization scheme based on the metric tensor (MT localization). We also explore how adaptive moving mesh tech-niques can control and inform the placement of mesh points to concentrate near the location of observations, reducingthe error of observation interpolation. This is especially beneficial when we have observations in locations that wouldotherwise have a sparse spatial discretization. We illustrate the utility of our results using discontinuous Galerkin(DG) approximations of 1D and 2D inviscid Burgers equations. The numerical results show that the MT localizationscheme compares favorably with standard Gaspari-Cohn localization techniques. In problems where the observationsare sparse, the choice of common mesh has a direct impact on DA performance. The numerical results also demonstratethe advantage of DG-based interpolation over linear interpolation for the 2D inviscid Burgers equation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A quasi-conservative DG-ALE method for multi-component flows using the non-oscillatory kinetic flux

A high-order quasi-conservative discontinuous Galerkin (DG) method is proposed for the numerical simulation of compressible multi-component flows. A distinct feature of the method is a predictor-corrector strategy to define the grid velocity. A Lagrangian mesh is first computed based on the flow velocity and then used as an initial mesh in a moving mesh method (the moving mesh partial differential equation or MMPDE method ) to improve its quality. The fluid dynamic equations are discretized in the direct arbitrary Lagrangian-Eulerian framework using DG elements and the non-oscillatory kinetic flux while the species equation is discretized using a quasi-conservative DG scheme to avoid numerical oscillations near material interfaces. A selection of one- and two-dimensional examples are presented to verify the convergence order and the constant-pressure-velocity preservation property of the method. They also demonstrate that the incorporation of the Lagrangian meshing with the MMPDE moving mesh method works well to concentrate mesh points in regions of shocks and material interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A well-balanced positivity-preserving quasi-Lagrange moving mesh DG method for the shallow water equations

A high-order, well-balanced, positivity-preserving quasi-Lagrange moving mesh DG method is presented for the shallow water equations with non-flat bottom topography. The well-balance property is crucial to the ability of a scheme to simulate perturbation waves over the lake-at-rest steady state such as waves on a lake or tsunami waves in the deep ocean. The method combines a quasi-Lagrange moving mesh DG method, a hydrostatic reconstruction technique, and a change of unknown variables. The strategies in the use of slope limiting, positivity-preservation limiting, and change of variables to ensure the well-balance and positivity-preserving properties are discussed. Compared to rezoning-type methods, the current method treats mesh movement continuously in time and has the advantages that it does not need to interpolate flow variables from the old mesh to the new one and places no constraint for the choice of an update scheme for the bottom topography on the new mesh. A selection of one- and two-dimensional examples are presented to demonstrate the well-balance property, positivity preservation, and high-order accuracy of the method and its ability to adapt the mesh according to features in the flow and bottom topography.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Anisotropic mesh quality measures and adaptation for polygonal meshes

Anisotropic mesh quality measures and anisotropic mesh adaptation are studied for polygonal meshes. Three sets of alignment and equidistribution measures are developed, one based on least squares fitting, one based on generalized barycentric mapping, and the other based on singular value decomposition of edge matrices. Numerical tests show that all three sets of mesh quality measures provide good measurements for the quality of polygonal meshes under given metrics. Based on one of these sets of quality measures and using a moving mesh partial differential equation, an anisotropic adaptive polygonal mesh method is constructed for the numerical solution of second order elliptic equations. Numerical examples are presented to demonstrate the effectiveness of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Domain Decomposition Parabolic Monge-Ampère Approach for Fast Generation of Adaptive Moving Meshes

A fast method is presented for adaptive moving mesh generation in multi-dimensions using a domain decomposition parabolic Monge-Ampère approach. The domain decomposition procedure employed here is non-iterative and involves splitting the computational domain into overlapping subdomains. An adaptive mesh on each subdomain is then computed as the image of the solution of the $L^2$ optimal mass transfer problem using a parabolic Monge-Ampère method. The domain decomposition approach allows straightforward implementation for the parallel computation of adaptive meshes which helps to reduce computational time significantly. Results are presented to show the numerical convergence of the domain decomposition solution to the single domain solution. Several numerical experiments are given to demonstrate the performance and efficiency of the proposed method. The numerical results indicate that the domain decomposition parabolic Monge-Ampère method is more efficient than the standard implementation of the parabolic Monge-Ampère method on the whole domain, in particular when computing adaptive meshes in three spatial dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Conditioning of implicit Runge-Kutta integration for finite element approximation of linear diffusion equations on anisotropic meshes

The conditioning of implicit Runge-Kutta (RK) integration for linear finite element approximation of diffusion equations on general anisotropic meshes is investigated. Bounds are established for the condition number of the resulting linear system with and without diagonal preconditioning for the implicit Euler and general implicit RK methods. Two solution strategies are considered for the linear system resulting from general implicit RK integration: the simultaneous solution (the system is solved as a whole) and a successive solution which follows the commonly used implementation of implicit RK methods to first transform the system into smaller systems using the Jordan normal form of the RK matrix and then solve them successively. For the simultaneous solution in case of a positive semidefinite symmetric part of the RK coefficient matrix and for the successive solution it is shown that . If the smallest eigenvalue of the symmetric part of the RK coefficient matrix is negative and the simultaneous solution strategy is used, an upper bound on the time step is given so that the system matrix is positive definite. The obtained bounds for the condition number have explicit geometric interpretations and take the interplay between the diffusion matrix and the mesh geometry into full consideration. They show that there are three mesh-dependent factors that can affect the conditioning: the number of elements, the mesh nonuniformity measured in the Euclidean metric, and the mesh nonuniformity with respect to the inverse of the diffusion matrix. They also reveal that the preconditioning using the diagonal of the system matrix, the mass matrix, or the lumped mass matrix can effectively eliminate the effects of the mesh nonuniformity measured in the Euclidean metric. Numerical examples are given.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

A quasi-Lagrangian moving mesh discontinuous Galerkin method for hyperbolic conservation laws

A moving mesh discontinuous Galerkin method is presented for the numerical solution of hyperbolic conservation laws. The method is a combination of the discontinuous Galerkin method and the mesh movement strategy which is based on the moving mesh partial differential equation approach and moves the mesh continuously in time and orderly in space. It discretizes hyperbolic conservation laws on moving meshes in the quasi-Lagrangian fashion with which the mesh movement is treated continuously and no interpolation is needed for physical variables from the old mesh to the new one. Two convection terms are induced by the mesh movement and their discretization is incorporated naturally in the DG formulation. Numerical results for a selection of one- and two-dimensional scalar and system conservation laws are presented. It is shown that the moving mesh DG method achieves the theoretically predicted order of convergence for problems with smooth solutions and is able to capture shocks and concentrate mesh points in non-smooth regions. Its advantage over uniform meshes and its insensitiveness to mesh smoothness are also demonstrated.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Adaptive Finite Element Solution of the Porous Medium Equation in Pressure Formulation

A lack of regularity in the solution of the porous medium equation poses a serious challenge in its theoretical and numerical studies. A common strategy in theoretical studies is to utilize the pressure formulation of the equation where a new variable called the mathematical pressure is introduced. It is known that the new variable has much better regularity than the original one and Darcy's law for the movement of the free boundary can be expressed naturally in this new variable. The pressure formulation has not been used in numerical studies. The goal of this work is to study its use in the adaptive finite element solution of the porous medium equation. The MMPDE moving mesh strategy is employed for adaptive mesh movement while linear finite elements are used for spatial discretization. The free boundary is traced explicitly by integrating Darcy's law with the Euler scheme. Numerical results are presented for three two-dimensional examples. The method is shown to be second-order in space and first-order in time in the pressure variable. Moreover, the convergence order of the error in the location of the free boundary is almost second-order in the maximum norm. However, numerical results also show that the convergence order of the method in the original variable stays between first-order and second-order in the $L^1$ norm or between 0.5th-order and first-order in the $L^2$ norm. Nevertheless, the current method can offer some advantages over numerical methods based on the original formulation for situations with large exponents or when a more accurate location of the free boundary is desired.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

An adaptive moving mesh discontinuous Galerkin method for the radiative transfer equation

The radiative transfer equation models the interaction of radiation with scattering and absorbing media and has important applications in various fields in science and engineering. It is an integro-differential equation involving time, space and angular variables and contains an integral term in angular directions while being hyperbolic in space. The challenges for its numerical solution include the needs to handle with its high dimensionality, the presence of the integral term, and the development of discontinuities and sharp layers in its solution along spatial directions. Its numerical solution is studied in this paper using an adaptive moving mesh discontinuous Galerkin method for spatial discretization together with the discrete ordinate method for angular discretization. The former employs a dynamic mesh adaptation strategy based on moving mesh partial differential equations to improve computational accuracy and efficiency. Its mesh adaptation ability, accuracy, and efficiency are demonstrated in a selection of one- and two-dimensional numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Conditioning of Finite Volume Element Method for Diffusion Problems with General Simplicial Meshes

The conditioning of the linear finite volume element discretization for general diffusion equations is studied on arbitrary simplicial meshes. The condition number is defined as the ratio of the maximal singular value of the stiffness matrix to the minimal eigenvalue of its symmetric part. This definition is motivated by the fact that the convergence rate of the generalized minimal residual method for the corresponding linear systems is determined by the ratio. An upper bound for the ratio is established by developing an upper bound for the maximal singular value and a lower bound for the minimal eigenvalue of the symmetric part. It is shown that the bound depends on three factors, the number of the elements in the mesh, the mesh nonuniformity measured in the Euclidean metric, and the mesh nonuniformity measured in the metric specified by the inverse diffusion matrix. It is also shown that the diagonal scaling can effectively eliminates the effects from the mesh nonuniformity measured in the Euclidean metric. Numerical results for a selection of examples in one, two, and three dimensions are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Moving mesh finite difference solution of non-equilibrium radiation diffusion equations

A moving mesh finite difference method based on the moving mesh partial differential equation is proposed for the numerical solution of the 2T model for multi-material, non-equilibrium radiation diffusion equations. The model involves nonlinear diffusion coefficients and its solutions stay positive for all time when they are positive initially. Nonlinear diffusion and preservation of solution positivity pose challenges in the numerical solution of the model. A coefficient-freezing predictor-corrector method is used for nonlinear diffusion while a cutoff strategy with a positive threshold is used to keep the solutions positive. Furthermore, a two-level moving mesh strategy and a sparse matrix solver are used to improve the efficiency of the computation. Numerical results for a selection of examples of multi-material non-equilibrium radiation diffusion show that the method is capable of capturing the profiles and local structures of Marshak waves with adequate mesh concentration. The obtained numerical solutions are in good agreement with those in the existing literature. Comparison studies are also made between uniform and adaptive moving meshes and between one-level and two-level moving meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Moving Mesh simulation of contact sets in two dimensional models of elastic-electrostatic deflection problems

Numerical and analytical methods are developed for the investigation of contact sets in electrostatic-elastic deflections modeling micro-electro mechanical systems. The model for the membrane deflection is a fourth-order semi-linear partial differential equation and the contact events occur in this system as finite time singularities. Primary research interest is in the dependence of the contact set on model parameters and the geometry of the domain. An adaptive numerical strategy is developed based on a moving mesh partial differential equation to dynamically relocate a fixed number of mesh points to increase density where the solution has fine scale detail, particularly in the vicinity of forming singularities. To complement this computational tool, a singular perturbation analysis is used to develop a geometric theory for predicting the possible contact sets. The validity of these two approaches are demonstrated with a variety of test cases.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

A New Functional for Variational Mesh Generation and Adaptation Based on Equidistribution and Alignment Conditions

A new functional is presented for variational mesh generation and adaptation. It is formulated based on combining the equidistribution and alignment conditions into a single condition with only one dimensionless parameter. The functional is shown to be coercive but not convex. A solution procedure using a discrete moving mesh partial differential equation is employed. It is shown that the element volumes and altitudes of a mesh trajectory of the mesh equation associated with the new functional are bounded away from zero and the mesh trajectory stays nonsingular if it is so initially. Numerical examples demonstrate that the new functional performs comparably as an existing one that is also based on the equidistribution and alignment conditions and known to work well but contains an additional parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

A study on nonnegativity preservation in finite element approximation of Nagumo-type nonlinear differential equations

Preservation of nonnengativity and boundedness in the finite element solution of Nagumo-type equations with general anisotropic diffusion is studied. Linear finite elements and the backward Euler scheme are used for the spatial and temporal discretization, respectively. An explicit, an implicit, and two hybrid explicit-implicit treatments for the nonlinear reaction term are considered. Conditions for the mesh and the time step size are developed for the numerical solution to preserve nonnegativity and boundedness. The effects of lumping of the mass matrix and the reaction term are also discussed. The analysis shows that the nonlinear reaction term has significant effects on the conditions for both the mesh and the time step size. Numerical examples are given to demonstrate the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Anisotropic mesh adaptation for 3D anisotropic diffusion problems with application to fractured reservoir simulation

Anisotropic mesh adaptation is studied for linear finite element solution of 3D anisotropic diffusion problems. The M-uniform mesh approach is used, where an anisotropic adaptive mesh is generated as a uniform one in the metric specified by a tensor. In addition to mesh adaptation, preservation of the maximum principle is also studied. Some new sufficient conditions for maximum principle preservation are developed, and a mesh quality measure is defined to server as a good indicator. Four different metric tensors are investigated: one is the identity matrix, one focuses on minimizing an error bound, another one on preservation of the maximum principle, while the fourth combines both. Numerical examples show that these metric tensors serve their purposes. Particularly, the fourth leads to meshes that improve the satisfaction of the maximum principle by the finite element solution while concentrating elements in regions where the error is large. Application of the anisotropic mesh adaptation to fractured reservoir simulation in petroleum engineering is also investigated, where unphysical solutions can occur and mesh adaptation can help improving the satisfaction of the maximum principle.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Selection of the Regularization Parameter in the Ambrosio-Tortorelli Approximation of the Mumford-Shah Functional for Image Segmentation

The Ambrosio-Tortorelli functional is a phase-field approximation of the Mumford-Shah functional that has been widely used for image segmentation. The approximation has the advantages of being easy to implement, maintaining the segmentation ability, and $Γ$-converging to the Mumford-Shah functional. However, it has been observed in actual computation that the segmentation ability of the Ambrosio-Tortorelli functional varies significantly with different values of the parameter and it even fails to $Γ$-converge to the original functional for some cases. In this paper we present an asymptotic analysis on the gradient flow equation of the Ambrosio-Tortorelli functional and show that the functional can have different segmentation behavior for small but finite values of the regularization parameter and eventually loses its segmentation ability as the parameter goes to zero when the input image is treated as a continuous function. This is consistent with the existing observation as well as the numerical examples presented in this work. A selection strategy for the regularization parameter and a scaling procedure for the solution are devised based on the analysis. Numerical results show that they lead to good segmentation of the Ambrosio-Tortorelli functional for real images.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Study on Moving Mesh Finite Element Solution of the Porous Medium Equation

An adaptive moving mesh finite element method is studied for the numerical solution of the porous medium equation with and without variable exponents and absorption. The method is based on the so-called moving mesh partial differential equation approach and employs its newly developed implementation. Three types of mesh are considered, uniform and arclength-based and Hessian-based adaptive meshes. The method shows a first order convergence for uniform and arclength-based adaptive meshes and a second-order convergence for Hessian-based adaptive meshes. It is also shown that the method can be used for situations with complex free boundaries, emerging and splitting of free boundaries, and the porous medium equation with variable exponents and absorption. Two dimensional numerical results are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stability of explicit one-step methods for P1-finite element approximation of linear diffusion equations on anisotropic meshes

We study the stability of explicit one-step integration schemes for the linear finite element approximation of linear parabolic equations. The derived bound on the largest permissible time step is tight for any mesh and any diffusion matrix within a factor of $2(d+1)$, where $d$ is the spatial dimension. Both full mass matrix and mass lumping are considered. The bound reveals that the stability condition is affected by two factors. The first one depends on the number of mesh elements and corresponds to the classic bound for the Laplace operator on a uniform mesh. The other factor reflects the effects of the interplay of the mesh geometry and the diffusion matrix. It is shown that it is not the mesh geometry itself but the mesh geometry in relation to the diffusion matrix that is crucial to the stability of explicit methods. When the mesh is uniform in the metric specified by the inverse of the diffusion matrix, the stability condition is comparable to the situation with the Laplace operator on a uniform mesh. Numerical results are presented to verify the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A comparative numerical study of meshing functionals for variational mesh adaptation

We present a comparative numerical study for three functionals used for variational mesh adaptation. One of them is a generalisation of Winslow's variable diffusion functional while the others are based on equidistribution and alignment. These functionals are known to have nice theoretical properties and work well for most mesh adaptation problems either as a stand-alone variational method or combined within the moving mesh framework. Their performance is investigated numerically in terms of equidistribution and alignment mesh quality measures. Numerical results in 2D and 3D are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A geometric discretization and a simple implementation for variational mesh generation and adaptation

We present a simple direct discretization for functionals used in the variational mesh generation and adaptation. Meshing functionals are discretized on simplicial meshes and the Jacobian matrix of the continuous coordinate transformation is approximated by the Jacobian matrices of affine mappings between elements. The advantage of this direct geometric discretization is that it preserves the basic geometric structure of the continuous functional, which is useful in preventing strong decoupling or loss of integral constraints satisfied by the functional. Moreover, the discretized functional is a function of the coordinates of mesh vertices and its derivatives have a simple analytical form, which allows a simple implementation of variational mesh generation and adaptation on computer. Since the variational mesh adaptation is the base for a number of adaptive moving mesh and mesh smoothing methods, the result in this work can be used to develop simple implementations of those methods. Numerical examples are given.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A hybrid LDG-HWENO scheme for KdV-type equations

A hybrid LDG-HWENO scheme is proposed for the numerical solution of KdV-type partial differential equations. It evolves the cell averages of the physical solution and its moments (a feature of Hermite WENO) while discretizes high order spatial derivatives using the local DG method. The new scheme has the advantages of both LDG and HWENO methods, including the ability to deal with high order spatial derivatives and the use of a small number of global unknown variables. The latter is independent of the order of the scheme and the spatial order of the underlying differential equations. One and two dimensional numerical examples are presented to show that the scheme can attain the same formal high order accuracy as the LDG method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Discrete maximum principle for the weak Galerkin method for anisotropic diffusion problems

A weak Galerkin discretization of the boundary value problem of a general anisotropic diffusion problem is studied for preservation of the maximum principle. It is shown that the direct application of the $M$-matrix theory to the stiffness matrix of the weak Galerkin discretization leads to a strong mesh condition requiring all of the mesh dihedral angles to be strictly acute (a constant-order away from 90 degrees). To avoid this difficulty, a reduced system is considered and shown to satisfy the discrete maximum principle under weaker mesh conditions. The discrete maximum principle is then established for the full weak Galerkin approximation using the relations between the degrees of freedom located on elements and edges. Sufficient mesh conditions for both piecewise constant and general anisotropic diffusion matrices are obtained. These conditions provide a guideline for practical mesh generation for preservation of the maximum principle. Numerical examples are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Monotone finite difference schemes for anisotropic diffusion problems via nonnegative directional splittings

Nonnegative directional splittings of anisotropic diffusion operators in the divergence form are investigated. Conditions are established for nonnegative directional splittings to hold in a neighborhood of an arbitrary interior point. The result is used to construct monotone finite difference schemes for the boundary value problem of anisotropic diffusion operators. It is shown that such a monotone scheme can be constructed if the underlying diffusion matrix is continuous on the closure of the physical domain and symmetric and uniformly positive definite on the domain, the mesh spacing is sufficiently small, and the size of finite difference stencil is sufficiently large. An upper bound for the stencil size is obtained, which is determined completely by the diffusion matrix. Loosely speaking, the more anisotropic the diffusion matrix is, the larger stencil is required. An exception is the situation with a strictly diagonally dominant diffusion matrix where a three-by-three stencil is sufficient for the construction of a monotone finite difference scheme. Numerical examples are presented to illustrate the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A study on anisotropic mesh adaptation for finite element approximation of eigenvalue problems with anisotropic diffusion operators

Anisotropic mesh adaptation is studied for the linear finite element solution of eigenvalue problems with anisotropic diffusion operators. The M-uniform mesh approach is employed with which any nonuniform mesh is characterized mathematically as a uniform one in the metric specified by a metric tensor. Bounds on the error in the computed eigenvalues are established for quasi-M-uniform meshes. Numerical examples arising from the Laplace-Beltrami operator on parameterized surfaces, nonlinear diffusion, thermal diffusion in a magnetic field in plasma physics, and the Laplacian operator on an L-shape domain are presented. Numerical results show that anisotropic adaptive meshes can lead to more accurate computed eigenvalues than uniform or isotropic adaptive meshes. They also confirm the second order convergence of the error that is predicted by the theoretical analysis. The effects of approximation of curved boundaries on the computation of eigenvalue problems is also studied in two dimensions. It is shown that the initial mesh used to define the geometry of the physical domain should contain at least \sqrt{N} boundary points to keep the effects of boundary approximation at the level of the error of the finite element approximation, where N is the number of the elements in the final adaptive mesh. Only about N^{1/3} boundary points in the initial mesh are needed for boundary value problems. This implies that the computation of eigenvalue problems is more sensitive to the boundary approximation than that of boundary value problems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

How a nonconvergent recovered Hessian works in mesh adaptation

Hessian recovery has been commonly used in mesh adaptation for obtaining the required magnitude and direction information of the solution error. Unfortunately, a recovered Hessian from a linear finite element approximation is nonconvergent in general as the mesh is refined. It has been observed numerically that adaptive meshes based on such a nonconvergent recovered Hessian can nevertheless lead to an optimal error in the finite element approximation. This also explains why Hessian recovery is still widely used despite its nonconvergence. In this paper we develop an error bound for the linear finite element solution of a general boundary value problem under a mild assumption on the closeness of the recovered Hessian to the exact one. Numerical results show that this closeness assumption is satisfied by the recovered Hessian obtained with commonly used Hessian recovery methods. Moreover, it is shown that the finite element error changes gradually with the closeness of the recovered Hessian. This provides an explanation on how a nonconvergent recovered Hessian works in mesh adaptation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Image Segmentation with Eigenfunctions of an Anisotropic Diffusion Operator

We propose the eigenvalue problem of an anisotropic diffusion operator for image segmentation. The diffusion matrix is defined based on the input image. The eigenfunctions and the projection of the input image in some eigenspace capture key features of the input image. An important property of the model is that for many input images, the first few eigenfunctions are close to being piecewise constant, which makes them useful as the basis for a variety of applications such as image segmentation and edge detection. The eigenvalue problem is shown to be related to the algebraic eigenvalue problems resulting from several commonly used discrete spectral clustering models. The relation provides a better understanding and helps developing more efficient numerical implementation and rigorous numerical analysis for discrete spectral segmentation methods. The new continuous model is also different from energy-minimization methods such as geodesic active contour in that no initial guess is required for in the current model. The multi-scale feature is a natural consequence of the anisotropic diffusion operator so there is no need to solve the eigenvalue problem at multiple levels. A numerical implementation based on a finite element method with an anisotropic mesh adaptation strategy is presented. It is shown that the numerical scheme gives much more accurate results on eigenfunctions than uniform meshes. Several interesting features of the model are examined in numerical examples and possible applications are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Maximum principle in linear finite element approximations of anisotropic diffusion-convection-reaction problems

A mesh condition is developed for linear finite element approximations of anisotropic diffusion-convection-reaction problems to satisfy a discrete maximum principle. Loosely speaking, the condition requires that the mesh be simplicial and $\mathcal{O}(\|\V{b}\|_\infty h + \|c\|_\infty h^2)$-nonobtuse when the dihedral angles are measured in the metric specified by the inverse of the diffusion matrix, where $h$ denotes the mesh size and $\V{b}$ and $c$ are the coefficients of the convection and reaction terms. In two dimensions, the condition can be replaced by a weaker mesh condition (an $\mathcal{O}(\|\V{b}\|_\infty h + \|c\|_\infty h^2)$ perturbation of a generalized Delaunay condition). These results include many existing mesh conditions as special cases. Numerical results are presented to verify the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Maximum principle for the finite element solution of time dependent anisotropic diffusion problems

Preservation of the maximum principle is studied for the combination of the linear finite element method in space and the $θ$-method in time for solving time dependent anisotropic diffusion problems. It is shown that the numerical solution satisfies a discrete maximum principle when all element angles of the mesh measured in the metric specified by the inverse of the diffusion matrix are non-obtuse and the time step size is bounded below and above by bounds proportional essentially to the square of the maximal element diameter. The lower bound requirement can be removed when a lumped mass matrix is used. In two dimensions, the mesh and time step conditions can be replaced by weaker Delaunay-type conditions. Numerical results are presented to verify the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Sign-preserving of principal eigenfunctions in P1 finite element approximation of eigenvalue problems of second-order elliptic operators

This paper is concerned with the P1 finite element approximation of the eigenvalue problem of second-order elliptic operators subject to the Dirichlet boundary condition. The focus is on the preservation of basic properties of the principal eigenvalue and eigenfunctions of continuous problems. It is shown that when the stiffness matrix is an irreducible $M$-matrix, the algebraic eigenvalue problem maintains those properties such as the smallest eigenvalue being real and simple and the corresponding eigenfunctions being either positive or negative inside the physical domain. Mesh conditions leading to such a stiffness matrix are also studied. A sufficient condition is that the mesh is simplicial, acute when measured in the metric specified by the inverse of the diffusion matrix, and interiorly connected. The acute requirement can be replaced by the Delaunay condition in two dimensions. Numerical results are presented to verify the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The cutoff method for the numerical computation of nonnegative solutions of parabolic PDEs with application to anisotropic diffusion and lubrication-type equations

The cutoff method, which cuts off the values of a function less than a given number, is studied for the numerical computation of nonnegative solutions of parabolic partial differential equations. A convergence analysis is given for a broad class of finite difference methods combined with cutoff for linear parabolic equations. Two applications are investigated, linear anisotropic diffusion problems satisfying the setting of the convergence analysis and nonlinear lubrication-type equations for which it is unclear if the convergence analysis applies. The numerical results are shown to be consistent with the theory and in good agreement with existing results in the literature. The convergence analysis and applications demonstrate that the cutoff method is an effective tool for use in the computation of nonnegative solutions. Cutoff can also be used with other discretization methods such as collocation, finite volume, finite element, and spectral methods and for the computation of positive solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Unconditionally stable high-order time integration for moving mesh finite difference solution of linear convection-diffusion equations

This paper is concerned with moving mesh finite difference solution of partial differential equations. It is known that mesh movement introduces an extra convection term and its numerical treatment has a significant impact on the stability of numerical schemes. Moreover, many implicit second and higher order schemes, such as the Crank-Nicolson scheme, will loss their unconditional stability. A strategy is presented for developing temporally high order, unconditionally stable finite difference schemes for solving linear convection-diffusion equations using moving meshes. Numerical results are given to demonstrate the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Adaptive finite elements with anisotropic meshes

The paper presents a numerical study for the finite element method with anisotropic meshes. We compare the accuracy of the numerical solutions on quasi-uniform, isotropic, and anisotropic meshes for a test problem which combines several difficulties of a corner singularity, a peak, a boundary layer, and a wavefront. Numerical experiment clearly shows the advantage of anisotropic mesh adaptation. The conditioning of the resulting linear equation system is addressed as well. In particular, it is shown that the conditioning with adaptive anisotropic meshes is not as bad as generally assumed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Conditioning of Finite Element Equations with Arbitrary Anisotropic Meshes

Bounds are developed for the condition number of the linear finite element equations of an anisotropic diffusion problem with arbitrary meshes. They depend on three factors. The first, factor proportional to a power of the number of mesh elements, represents the condition number of the linear finite element equations for the Laplacian operator on a uniform mesh. The other two factors arise from the mesh nonuniformity viewed in the Euclidean metric and in the metric defined by the diffusion matrix. The new bounds reveal that the conditioning of the finite element equations with adaptive anisotropic meshes is much better than what is commonly feared. Diagonal scaling for the linear system and its effects on the conditioning are also studied. It is shown that the Jacobi preconditioning, which is an optimal diagonal scaling for a symmetric positive definite sparse matrix, can eliminate the effects of mesh nonuniformity viewed in the Euclidean metric and reduce those effects of the mesh viewed in the metric defined by the diffusion matrix. Tight bounds on the extreme eigenvalues of the stiffness and mass matrices are obtained. Numerical examples are given.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An anisotropic mesh adaptation method for the finite element solution of heterogeneous anisotropic diffusion problems

Heterogeneous anisotropic diffusion problems arise in the various areas of science and engineering including plasma physics, petroleum engineering, and image processing. Standard numerical methods can produce spurious oscillations when they are used to solve those problems. A common approach to avoid this difficulty is to design a proper numerical scheme and/or a proper mesh so that the numerical solution validates the discrete counterpart (DMP) of the maximum principle satisfied by the continuous solution. A well known mesh condition for the DMP satisfaction by the linear finite element solution of isotropic diffusion problems is the non-obtuse angle condition that requires the dihedral angles of mesh elements to be non-obtuse. In this paper, a generalization of the condition, the so-called anisotropic non-obtuse angle condition, is developed for the finite element solution of heterogeneous anisotropic diffusion problems. The new condition is essentially the same as the existing one except that the dihedral angles are now measured in a metric depending on the diffusion matrix of the underlying problem. Several variants of the new condition are obtained. Based on one of them, two metric tensors for use in anisotropic mesh generation are developed to account for DMP satisfaction and the combination of DMP satisfaction and mesh adaptivity. Numerical examples are given to demonstrate the features of the linear finite element method for anisotropic meshes generated with the metric tensors.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Anisotropic Mesh Adaptation for Variational Problems Using Error Estimation Based on Hierarchical Bases

Anisotropic mesh adaptation has been successfully applied to the numerical solution of partial differential equations but little considered for variational problems. In this paper, we investigate the use of a global hierarchical basis error estimator for the development of an anisotropic metric tensor needed for the adaptive finite element solution of variational problems. The new metric tensor is completely a~posteriori and based on residual, edge jumps and the hierarchical basis error estimator. Numerical results show that it performs comparable with existing metric tensors based on Hessian recovery. A few sweeps of the symmetric Gauß-Seidel iteration for solving the global error problem prove sufficient to provide directional information necessary for successful mesh adaptation. .

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Discrete maximum principle and a Delaunay-type mesh condition for linear finite element approximations of two-dimensional anisotropic diffusion problems

The finite element solution of two-dimensional anisotropic diffusion problems is considered. A Delaunay-type mesh condition is developed for linear finite element approximations to satisfy a discrete maximum principle. The condition is shown to be weaker than the existing anisotropic non-obtuse angle condition. It reduces to the well known Delaunay condition for the special case with the identity diffusion matrix. Numerical results are presented to verify the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A posteriori error analysis for finite element solution of elliptic differential equations using equidistributing meshes

The paper is concerned with the adaptive finite element solution of linear elliptic differential equations using equidistributing meshes. A strategy is developed for defining this type of mesh based on residual-based a posteriori error estimates and rigorously analyzing the convergence of a linear finite element approximation using them. The existence and computation of equidistributing meshes and the continuous dependence of the finite element approximation on mesh are also studied. Numerical results are given to verify the theoretical findings.

0 citations0 reviewsNo code linkedOpen access