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

Jia Zhao

Jia Zhao contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.RAmath.RTmath.SPNeural and Evolutionary Computing

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
6topics
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

11 published item(s)

preprint2022arXiv

Cohomology and deformations of Relative Rota-Baxter operators on Lie-Yamaguti algebras

In this paper, we establish the cohomology of relative Rota-Baxter operators on Lie-Yamaguti algebras via the Yamaguti cohomology. Then we use this type of cohomology to characterize deformations of relative Rota-Baxter operators on Lie-Yamaguti algebras. We show that if two linear or formal deformations of a relative Rota-Baxter operator are equivalent, then their infinitesimals are in the same cohomology class in the first cohomology group. Moreover, an order $n$ deformation of a relative Rota-Baxter operator can be extended to an order $n+1$ deformation if and only if the obstruction class in the second cohomology group is trivial.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical Approximations of the Allen-Cahn-Ohta-Kawasaki (ACOK) Equation with Modified Physics Informed Neural Networks (PINNs)

The physics informed neural networks (PINNs) has been widely utilized to numerically approximate PDE problems. While PINNs has achieved good results in producing solutions for many partial differential equations, studies have shown that it does not perform well on phase field models. In this paper, we partially address this issue by introducing a modified physics informed neural networks. In particular, they are used to numerically approximate Allen-Cahn-Ohta-Kawasaki (ACOK) equation with a volume constraint.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A General Framework to Derive Linear, Decoupled and Energy-stable Schemes for Reversible-Irreversible Thermodynamically Consistent Models: Part I Incompressible Hydrodynamic Models

In this paper, we present a general numerical platform for designing accurate, efficient, and stable numerical algorithms for incompressible hydrodynamic models that obeys the thermodynamical laws. The obtained numerical schemes are automatically linear in time. It decouples the hydrodynamic variable and other state variables such that only small-size linear problems need to be solved at each time marching step. Furthermore, if the classical velocity projection method is utilized, the velocity field and pressure field can be decoupled. In the end, only a few elliptic-type equations shall be solved in each time step. This strategy is made possible through a sequence of model reformulations by fully exploring the models' thermodynamic structures. The generalized Onsager principle directly guides these reformulation procedures. In the reformulated but equivalent models, the reversible and irreversible components can be identified, guiding the numerical platform to decouple the reversible and irreversible dynamics. This eventually leads to decoupled numerical algorithms, given that the coupling terms only involve irreversible dynamics. To further demonstrate the numerical platform's power, we apply it to several specific incompressible hydrodynamic models. The energy stability of the proposed numerical schemes is shown in detail. The second-order accuracy in time is verified numerically through time step refinement tests. Several benchmark numerical examples are presented to further illustrate the proposed numerical framework's accuracy, stability, and efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Eigenvalues of Laplace operators on non-bipartite graphs

This paper considers the comparison between the eigenvalues of Laplace operators with the standard conditions and the anti-standard conditions on non-bipartite graphs which are equilateral or inequilateral. First of all, we show the calculation of the eigenvalues of Laplace operators on equilateral metric graphs with arbitrary edge length. Based on this method, we use the properties of the cosine function and the arccosine function to find the comparison between the eigenvalues of Laplace operators with the standard conditions and the anti-standard conditions on equilateral non-bipartite graphs. In addition, we give the inequalities between standard and anti-standard eigenvalues on a special inequilateral non-bipartite graph.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Improving the Accuracy and Consistency of the Scalar Auxiliary Variable (SAV) Method with Relaxation

The scalar auxiliary variable (SAV) method was introduced by Shen et al. and has been broadly used to solve thermodynamically consistent PDE problems. By utilizing scalar auxiliary variables, the original PDE problems are reformulated into equivalent PDE problems. The advantages of the SAV approach, such as linearity, unconditionally energy stability, and easy-to-implement, are prevalent. However, there is still an open issue unresolved, i.e., the numerical schemes resulted from the SAV method preserve a "modified" energy law according to the auxiliary variables instead of the original variables. Truncation errors are introduced during numerical calculations so that the numerical solutions of the auxiliary variables are no longer equivalent to their original continuous definitions. In other words, even though the SAV scheme satisfies a modified energy law, it does not necessarily satisfy the energy law of the original PDE models. This paper presents one essential relaxation technique to overcome this issue, which we named the relaxed-SAV (RSAV) method. Our RSAV method penalizes the numerical errors of the auxiliary variables by a relaxation technique. In general, the RSAV method keeps all the advantages of the baseline SAV method and improves its accuracy and consistency noticeably. Several examples have been presented to demonstrate the effectiveness of the RSAV approach.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Non-uniform Time-stepping Convex Splitting Scheme for the Time-fractional Cahn-Hilliard Equation

In this paper, a non-uniform time-stepping convex-splitting numerical algorithm for solving the widely used time-fractional Cahn-Hilliard equation is introduced. The proposed numerical scheme employs the $L1^+$ formula for discretizing the time-fractional derivative and a second-order convex-splitting technique to deal with the non-linear term semi-implicitly. Then the pseudospectral method is utilized for spatial discretization. As a result, the fully discrete scheme has several advantages: second-order accurate in time, spectrally accurate in space, uniquely solvable, mass preserving, and unconditionally energy stable. Rigorous proofs are given, along with several numerical results to verify the theoretical results, and to show the accuracy and effectiveness of the proposed scheme. Also, some interesting phase separation dynamics of the time-fractional Cahn-Hilliard equation has been investigated.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Discovering Phase Field Models from Image Data with the Pseudo-spectral Physics Informed Neural Networks

In this paper, we introduce a new deep learning framework for discovering the phase field models from existing image data. The new framework embraces the approximation power of physics informed neural networks (PINN), and the computational efficiency of the pseudo-spectral methods, which we named pseudo-spectral PINN or SPINN. Unlike the baseline PINN, the pseudo-spectral PINN has several advantages. First of all, it requires less training data. A minimum of two snapshots with uniform spatial resolution would be adequate. Secondly, it is computationally efficient, as the pseudo-spectral method is used for spatial discretization. Thirdly, it requires less trainable parameters compared with the baseline PINN. Thus, it significantly simplifies the training process and assures less local minima or saddle points. We illustrate the effectiveness of pseudo-spectral PINN through several numerical examples. The newly proposed pseudo-spectral PINN is rather general, and it can be readily applied to discover other PDE-based models from image data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Solving Allen-Cahn and Cahn-Hilliard Equations using the Adaptive Physics Informed Neural Networks

Phase field models, in particular, the Allen-Cahn type and Cahn-Hilliard type equations, have been widely used to investigate interfacial dynamic problems. Designing accurate, efficient, and stable numerical algorithms for solving the phase field models has been an active field for decades. In this paper, we focus on using the deep neural network to design an automatic numerical solver for the Allen-Cahn and Cahn-Hilliard equations by proposing an improved physics informed neural network (PINN). Though the PINN has been embraced to investigate many differential equation problems, we find a direct application of the PINN in solving phase-field equations won't provide accurate solutions in many cases. Thus, we propose various techniques that add to the approximation power of the PINN. As a major contribution of this paper, we propose to embrace the adaptive idea in both space and time and introduce various sampling strategies, such that we are able to improve the efficiency and accuracy of the PINN on solving phase field equations. In addition, the improved PINN has no restriction on the explicit form of the PDEs, making it applicable to a wider class of PDE problems, and shedding light on numerical approximations of other PDEs in general.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Arbitrarily High-order Linear Schemes for Gradient Flow Models

We present a paradigm for developing arbitrarily high order, linear, unconditionally energy stable numerical algorithms for gradient flow models. We apply the energy quadratization (EQ) technique to reformulate the general gradient flow model into an equivalent gradient flow model with a quadratic free energy and a modified mobility. Given solutions up to $t_n=n Δt$ with $Δt$ the time step size, we linearize the EQ-reformulated gradient flow model in $(t_n, t_{n+1}]$ by extrapolation. Then we employ an algebraically stable Runge-Kutta method to discretize the linearized model in $(t_n, t_{n+1}]$. Then we use the Fourier pseudo-spectral method for the spatial discretization to match the order of accuracy in time. The resulting fully discrete scheme is linear, unconditionally energy stable, uniquely solvable, and may reach arbitrarily high order. Furthermore, we present a family of linear schemes based on prediction-correction methods to complement the new linear schemes. Some benchmark numerical examples are given to demonstrate the accuracy and efficiency of the schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Arbitrarily High-order Unconditionally Energy Stable Schemes for Thermodynamically Consistent Gradient Flow Models

We present a systematical approach to developing arbitrarily high order, unconditionally energy stable numerical schemes for thermodynamically consistent gradient flow models that satisfy energy dissipation laws. Utilizing the energy quadratization (EQ) method, We formulate the gradient flow model into an equivalent form with a corresponding quadratic free energy functional. Based on the equivalent form with a quadratic energy, we propose two classes of energy stable numerical approximations. In the first approach, we use a prediction-correction strategy to improve the accuracy of linear numerical schemes. In the second approach, we adopt the Gaussian collocation method to discretize the equivalent form with a quadratic energy, arriving at an arbitrarily high-order scheme for gradient flow models. Schemes derived using both approaches are proved rigorously to be unconditionally energy stable. The proposed schemes are then implemented in four gradient flow models numerically to demonstrate their accuracy and effectiveness. Detailed numerical comparisons among these schemes are carried out as well. These numerical strategies are rather general so that they can be readily generalized to solve any thermodynamically consistent PDE models.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

An Accurate and Efficient Algorithm for The Time-fractional Molecular Beam Epitaxy Model with Slope Selection

In this paper, we propose a time-fractional molecular beam epitaxy (MBE) model with slope selection and its efficient, accurate, full discrete, linear numerical approximation. The numerical scheme utilizes the fast algorithm for the Caputo fractional derivative operator in time discretization and Fourier spectral method in spatial discretization. Refinement tests are conducted to verify the $2-α$ order of time convergence, with $α\in (0, 1]$ the fractional order of derivative. Several numerical simulations are presented to demonstrate the accuracy and efficiency of our newly proposed scheme. By exploring the fast algorithm calculating the Caputo fractional derivative, our numerical scheme makes it practice for long time simulation of MBE coarsening, which is essential for MBE model in practice. With the proposed fractional MBE model, we observe that the scaling law for the energy decays as $ O(t^{-\fracα{3}})$ and the roughness increases as $O(t^{\fracα{3}})$, during the coarsening dynamics with random initial condition. That is to say, the coarsening rate of MBE model could be manipulated by the fractional order $α$, and it is linearly proportional to $α$. This is the first time in literature to report/discover such scaling correlation. It provides a potential application field for fractional differential equations. Besides, the numerical approximation strategy proposed in this paper can be readily applied to study many classes of time-fractional and high dimensional phase field models.

0 citations0 reviewsNo code linkedOpen access