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

Huazhong Tang

Huazhong Tang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisastro-ph.IMhep-thmath-phmath.MPnlin.SIphysics.comp-phphysics.flu-dynphysics.plasm-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2022arXiv

A general class of linear unconditionally energy stable schemes for the gradient flows

This paper studies a class of linear unconditionally energy stable schemes for the gradient flows. Such schemes are built on the SAV technique and the general linear time discretization (GLTD) as well as the linearization based on the extrapolation for the nonlinear term, and may be arbitrarily high-order accurate and very general, containing many existing SAV schemes and new SAV schemes. It is shown that the semi-discrete-in-time schemes are unconditionally energy stable when the GLTD is algebraically stable, and are convergent with the order of $\min\{\hat{q},ν\}$ under the diagonal stability and some suitable regularity and accurate starting values, where $\hat{q}$ is the generalized stage order of the GLTD and $ν$ denotes the number of the extrapolation points in time. The energy stability results can be easily extended to the fully discrete schemes, for example, if the Fourier spectral method is employed in space when the periodic boundary conditions are specified. Some numerical experiments on the Allen-Cahn, Cahn-Hilliard, and phase field crystal models are conducted to validate those theories as well as the effectiveness, the energy stability and the accuracy of our schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

High-order accurate entropy stable adaptive moving mesh finite difference schemes for (multi-component) compressible Euler equations with the stiffened equation of state

This paper extends the high-order entropy stable (ES) adaptive moving mesh finite difference schemes developed in [14] to the two- and three-dimensional (multi-component) compressible Euler equations with the stiffened equation of state. The two-point entropy conservative (EC) flux is first constructed in the curvilinear coordinates. The high-order semi-discrete EC schemes are given with the aid of the two-point EC flux and the high-order discretization of the geometric conservation laws, and then the high-order semi-discrete ES schemes satisfying the entropy inequality are derived by adding the high-order dissipation term based on the multi-resolution weighted essentially non-oscillatory (WENO) reconstruction for the scaled entropy variables to the EC schemes. The explicit strong-stability-preserving Runge-Kutta methods are used for the time discretization and the mesh points are adaptively redistributed by iteratively solving the mesh redistribution equations with an appropriately chosen monitor function. Several 2D and 3D numerical tests are conducted on the parallel computer system with the MPI programming to validate the accuracy and the ability to capture effectively the localized structures of the proposed schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Positivity-Preserving Well-Balanced Central Discontinuous Galerkin Schemes for the Euler Equations under Gravitational Fields

This paper designs and analyzes positivity-preserving well-balanced (WB) central discontinuous Galerkin (CDG) schemes for the Euler equations with gravity. A distinctive feature of these schemes is that they not only are WB for a general known stationary hydrostatic solution, but also can preserve the positivity of the fluid density and pressure. The standard CDG method does not possess this feature, while directly applying some existing WB techniques to the CDG framework may not accommodate the positivity and keep other important properties at the same time. In order to obtain the WB and positivity-preserving properties simultaneously while also maintaining the conservativeness and stability of the schemes, a novel spatial discretization is devised in the CDG framework based on suitable modifications to the numerical dissipation term and the source term approximation. The modifications are based on a crucial projection operator for the stationary hydrostatic solution, which is proposed for the first time in this work. This novel projection has the same order of accuracy as the standard $L^2$-projection, can be explicitly calculated, and is easy to implement without solving any optimization problems. More importantly, it ensures that the projected stationary solution has the same cell averages on both the primal and dual meshes, which is a key to achieve the desired properties of our schemes. Based on some convex decomposition techniques, rigorous positivity-preserving analyses for the resulting WB CDG schemes are carried out. Several one- and two-dimensional numerical examples are performed to illustrate the desired properties of these schemes, including the high-order accuracy, the WB property, the robustness for simulations involving the low pressure or density, high resolution for the discontinuous solutions and the small perturbations around the equilibrium state.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An analytical solution of the isentropic vortex problem in the special relativistic magnetohydrodynamics

The isentropic vortex problem is frequently solved to test the accuracy of numerical methods and verify corresponding code. Unfortunately, its existing solution was derived in the relativistic magnetohydrodynamics by numerically solving an ordinary differential equation. This note provides an analytical solution of the 2D isentropic vortex problem with explicit algebraic expressions in the special relativistic hydrodynamics and magnetohydrodynamics and extends it to the 3D case.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

High-order accurate entropy stable adaptive moving mesh finite difference schemes for special relativistic (magneto)hydrodynamics

This paper develops high-order accurate entropy stable (ES) adaptive moving mesh finite difference schemes for the two- and three-dimensional special relativistic hydrodynamic (RHD) and magnetohydrodynamic (RMHD) equations, which is the high-order accurate extension of [J.M. Duan and H.Z. Tang, Entropy stable adaptive moving mesh schemes for 2D and 3D special relativistic hydrodynamics, J. Comput. Phys., 426(2021), 109949]. The key point is the derivation of the higher-order accurate entropy conservative (EC) and ES finite difference schemes in the curvilinear coordinates by carefully dealing with the discretization of the temporal and spatial metrics and the Jacobian of the coordinate transformation and constructing the high-order EC and ES fluxes with the discrete metrics. The spatial derivatives in the source terms of the symmetrizable RMHD equations and the geometric conservation laws are discretized by using the linear combinations of the corresponding second-order case to obtain high-order accuracy. Based on the proposed high-order accurate EC schemes and the high-order accurate dissipation terms built on the WENO reconstruction, the high-order accurate ES schemes are obtained for the RHD and RMHD equations in the curvilinear coordinates. The mesh iteration redistribution or adaptive moving mesh strategy is built on the minimization of the mesh adaption functional. Several numerical tests are conducted to validate the shock-capturing ability and high efficiency of our high-order accurate ES adaptive moving mesh methods on the parallel computer system with the MPI communication. The numerical results show that the high-order accurate ES adaptive moving mesh schemes outperform both their counterparts on the uniform mesh and the second-order ES adaptive moving mesh schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Entropy stable adaptive moving mesh schemes for 2D and 3D special relativistic hydrodynamics

This paper develops entropy stable (ES) adaptive moving mesh schemes for the 2D and 3D special relativistic hydrodynamic (RHD) equations. They are built on the ES finite volume approximation of the RHD equations in curvilinear coordinates, the discrete geometric conservation laws, and the mesh adaptation implemented by iteratively solving the Euler-Lagrange equations of the mesh adaption functional in the computational domain with suitably chosen monitor functions. First, a sufficient condition is proved for the two-point entropy conservative (EC) flux, by mimicking the derivation of the continuous entropy identity in curvilinear coordinates and using the discrete geometric conservation laws given by the conservative metrics method. Based on such sufficient condition, the EC fluxes for the RHD equations in curvilinear coordinates are derived and the second-order accurate semi-discrete EC schemes are developed to satisfy the entropy identity for the given convex entropy pair. Next, the semi-discrete ES schemes satisfying the entropy inequality are proposed by adding a suitable dissipation term to the EC scheme and utilizing linear reconstruction with the minmod limiter in the scaled entropy variables in order to suppress the numerical oscillations of the above EC scheme. Then, the semi-discrete ES schemes are integrated in time by using the second-order strong stability preserving explicit Runge-Kutta schemes. Finally, several numerical results show that our 2D and 3D ES adaptive moving mesh schemes effectively capture the localized structures, such as sharp transitions or discontinuities, and are more efficient than their counterparts on uniform mesh.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the explicit two-stage fourth-order accurate time discretizations

This paper continues to study the explicit two-stage fourth-order accurate time discretiza- tions [5, 7]. By introducing variable weights, we propose a class of more general explicit one-step two-stage time discretizations, which are different from the existing methods, such as the Euler methods, Runge-Kutta methods, and multistage multiderivative methods etc. We study the absolute stability, the stability interval, and the intersection between the imaginary axis and the absolute stability region. Our results show that our two-stage time discretizations can be fourth-order accurate conditionally, the absolute stability region of the proposed methods with some special choices of the variable weights can be larger than that of the classical explicit fourth- or fifth-order Runge-Kutta method, and the interval of absolute stability can be almost twice as much as the latter. Several numerical experiments are carried out to demonstrate the performance and accuracy as well as the stability of our proposed methods

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

High-order accurate entropy stable finite difference schemes for one- and two-dimensional special relativistic hydrodynamics

This paper develops the high-order accurate entropy stable finite difference schemes for one- and two-dimensional special relativistic hydrodynamic equations. The schemes are built on the entropy conservative flux and the weighted essentially non-oscillatory (WENO) technique as well as explicit Runge-Kutta time discretization. The key is to technically construct the affordable entropy conservative flux of the semi-discrete second-order accurate entropy conservative schemes satisfying the semi-discrete entropy equality for the found convex entropy pair. As soon as the entropy conservative flux is derived, the dissipation term can be added to give the semi-discrete entropy stable schemes satisfying the semi-discrete entropy inequality with the given convex entropy function. The WENO reconstruction for the scaled entropy variables and the high-order explicit Runge-Kutta time discretization are implemented to obtain the fully-discrete high-order schemes. Several numerical tests are conducted to validate the accuracy and the ability to capture discontinuities of our entropy stable schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Model reduction of a kinetic swarming model by operator projection

This paper derives the arbitrary order globally hyperbolic moment system for a non-linear kinetic description of the Vicsek swarming model by using the operator projection. It is built on our careful study of a family of the complicate Grad type orthogonal functions depending on a parameter (angle of macroscopic velocity). We calculate their derivatives with respect to the independent vari- able, and projection of those derivatives, the product of velocity and basis, and collision term. The moment system is also proved to be hyperbolic, rotational invariant, and mass-conservative. The relationship between Grad type expansions in different parameter is also established. A semi-implicit numerical scheme is presented to solve a Cauchy problem of our hyperbolic moment system in order to verify the convergence behavior of the moment method. It is also compared to the spectral method for the kinetic equation. The results show that the solutions of our hyperbolic moment system converge to the solutions of the kinetic equation for the Vicsek model as the order of the moment system increases, and the moment method can capture key features such as rarefaction and shock waves, contact discontinuity, and vortex formation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Multi-hump solitary waves of nonlinear Dirac equation

This paper concentrates on a (1+1)-dimensional nonlinear Dirac (NLD) equation with a general self-interaction, being a linear combination of the scalar, pseudoscalar, vector and axial vector self-interactions to the power of the integer $k+1$. The solitary wave solutions to the NLD equation are analytically derived, and the upper bounds of the hump number in the charge, energy and momentum densities for the solitary waves are proved in theory. The results show that: (1) for a given integer $k$, the hump number in the charge density is not bigger than $4$, while that in the energy density is not bigger than $3$; (2) those upper bounds can only be achieved in the situation of higher nonlinearity, namely, $k\in\{5,6,7,\cdots \}$ for the charge density and $k\in\{3,5,7,\cdots\}$ for the energy density; (3) the momentum density has the same multi-hump structure as the energy density; (4) more than two humps (resp. one hump) in the charge (resp. energy) density can only happen under the linear combination of the pseudoscalar self-interaction and at least one of the scalar and vector (or axial vector) self-interactions. Our results on the multi-hump structure will be interesting in the interaction dynamics for the NLD solitary waves.

0 citations0 reviewsNo code linkedOpen access