Discrete Vector Calculus and Helmholtz Hodge Decomposition for Classical Finite Difference Summation by Parts Operators

In this article,discrete variants of several results from vector calculus are studied for clas-sical finite difference summation by parts operators in two and three space dimensions.It is shown that existence theorems for scalar/vector potentials of irrotational/solenoidal vector fields cannot hold discretely because of grid oscillations,which are characterised explicitly.This results in a non-vanishing remainder associated with grid oscillations in the discrete Helmholtz Hodge decomposition.Nevertheless,iterative numerical methods based on an interpretation of the Helmholtz Hodge decomposition via orthogonal projections are pro-posed and applied successfully.In numerical experiments,the discrete remainder vanishes and the potentials converge with the same order of accuracy as usual in other first-order partial differential equations.Motivated by the successful application of the Helmholtz Hodge decomposition in theoretical plasma physics,applications to the discrete analysis of magnetohydrodynamic(MHD) wave modes are presented and discussed.

High-Order Method with Moving Frames to Compute the Covariant Derivatives of Vectors on General 2D Curved Surfaces

The covariant derivative is a generalization of differentiating vectors.The Euclidean derivative is a special case of the covariant derivative in Euclidean space.The covariant derivative gathers broad attention,particularly when computing vector derivatives on curved surfaces and volumes in various applications.Covariant derivatives have been computed using the metric tensor from the analytically known curved axes.However,deriving the global axis for the domain has been mathematically and computationally challenging for an arbitrary two-dimensional(2D)surface.Consequently,computing the covariant derivative has been difficult or even impossible.A novel high-order numerical scheme is proposed for computing the covariant derivative on any 2D curved surface.A set of orthonormal vectors,known as moving frames,expand vectors to compute accurately covariant derivatives on 2D curved surfaces.The proposed scheme does not require the construction of curved axes for the metric tensor or the Christoffel symbols.The connectivity given by the Christoffel symbols is equivalently provided by the attitude matrix of orthonormal moving frames.Consequently,the proposed scheme can be extended to the general 2D curved surface.As an application,the Helmholtz‐Hodge decomposition is considered for a realistic atrium and a bunny.

ON VERY WEAK SOLUTIONS OF A-HARMONICEQUATION WITH VERY WEAK BOUNDARYVALUES

In this paper,the following result is given by using Hodge decomposition: There exists r(0) = r(0)(n,p,a,b), such that if u is an element of W-loc(1,r)(Omega) is a very weak solution of (1.1),with C max{1,p - 1} < r < p and u is an element of W-0(1,r)(Omega;partial derivativeOmega\E) where E subset of partial derivativeOmega is a closed set and small in an appropriate capacity sense, then u = 0, a.e. in Omega provided that r(0) < r < p.

Constraint-induced restriction and extension operators with applications

The Stokes operator is a differential-integral operator induced by the Stokes equations. In this paper, we analyze the Stokes operator from the point of view of the Helmholtz minimum dissipation principle. We show that, through the Hodge orthogonal decomposition, a pair of bounded linear operators, a restriction operator and an extension operator, are induced from the divergence-free constraint. As a consequence, we use it to calculate the eigenvalues of the Stokes operator.

REGULARITY FOR VERY WEAK SOLUTIONS TO A-HARMONIC EQUATION

In this paper, the following result is given by using Hodge decomposition and weak reverse Holder inequality： For every r1 with P-（2^n＋1 100^n^2 p（2^3＋n/（P-1）＋1）b/a）^-1〈r1〈p,there exists the exponent r2 = r2（n, r1,p） 〉 p, such that for every very weak solution u∈W^1r1,loc（Ω） to A-harmonic equation, u also belongs to W^1r2,loc（Ω） . In particular, u is the weak solution to A-harmonic equation in the usual sense.

Partial Regularity for Two Dimensional Landau-Lifshitz Equations

It is proved that any weak solution to the initial value problem of two dimensional Landau-Lifshitz equation is unique and is smooth with the exception of at most finitely many points, provided that the weak solution has finite energy.

Regularity for weakly(K_(1),K_(2))-quasiregular mappings

In this paper, we first give the definition of weakly (K1,K2)-quasiregular mappings, and then by using the Hodge decomposition and the weakly reverse Holder inequality, we obtain their regularity property: For any q1 that satisfies 0＜K1n(n+4)/22n+1×100n2[23n/2(25n+1)](n - q1) ＜ 1, there exists p1 = p1(n,q1,K1,K2)＞n, such that any (K1,K2)-quasiregular mapping f ∈ W1,q1loc(Ω,Rn) is in fact in W1n,p1loc (Ω, Rn). That is, f is (K1, K2)-quasiregular in the usual sense.

Regularity Results for the Generalized Beltrami System

For the generalized Beltrami system with two characteristic matrices, we deal with the regularity of its very weak solutions in the Sobolev class (1 < r < n). By changing the generalized Beltrami system into a class of a divergent elliptic system with nonhomogeneous items, we obtain that each of its very weak solutions is essentially a classical weak solution of a usual Sobolev class. Furthermore, we also establish a higher regularity of its weak solution if the regularity hypotheses of two characteristic matrices are improved.