BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Weizhang Huang

Weizhang Huang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical AnalysisComputational Engineering, Finance, and Sciencemath.AP

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
21works
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

21 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
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

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