Numerical Investigations on the Resonance Errors of Multiscale Discontinuous Galerkin Methods for One-Dimensional Stationary Schrödinger Equation

In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al.in J Sci Comput 66:321–345,2016;Dong and Wang in J Comput Appl Math 380:1–11,2020)for a one-dimensional stationary Schrödinger equation.Previous work showed that penalty parameters were required to be positive in error analysis,but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes.In this work,by performing extensive numerical experiments,we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods,and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.

A Posteriori Error Estimate of Two Grid Mixed Finite Element Methods for Semilinear Elliptic Equations

In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the mixed method equations. Then, the averaging technique is used to construct the a posteriori error estimates of the two-grid mixed finite element method and theoretical analysis are given for the error estimators. Finally, we give some numerical examples to verify the reliability and efficiency of the a posteriori error estimator.

Error estimates of H^1-Galerkin mixed finite element method for Schrdinger equation

An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.

Error Estimates for Mixed Finite Element Methods for Sobolev Equation

The purpose of this paper is to investigate the convergence of the mixed finite element method for the initial-boundary value problem for the Sobolev equation Ut-div{aut + b1 u} = f based on the Raviart-Thomas space Vh × Wh H(div; × L2(). Optimal order estimates are obtained for the approximation of u, ut, the associated velocity p and divp respectively in L(0,T;L2()), L(0,T;L2()), L(0,T;L2()2), and L2(0, T; L2()). Quasi-optimal order estimates are obtained for the approximations of u, ut in L(0, T; L()) and p in L(0,T; L()2).

ANALYSIS AND DISCRETIZATION FOR AN OPTIMAL CONTROL PROBLEM OF A VARIABLE-COEFFICIENT RIESZ-FRACTIONAL DIFFUSION EQUATION WITH POINTWISE CONTROL CONSTRAINTS

We present a mathematical and numerical study for a pointwise optimal control problem governed by a variable-coefficient Riesz-fractional diffusion equation.Due to the impact of the variable diffusivity coefficient,existing regularity results for their constantcoefficient counterparts do not apply,while the bilinear forms of the state(adjoint)equation may lose the coercivity that is critical in error estimates of the finite element method.We reformulate the state equation as an equivalent constant-coefficient fractional diffusion equation with the addition of a variable-coefficient low-order fractional advection term.First order optimality conditions are accordingly derived and the smoothing properties of the solutions are analyzed by,e.g.,interpolation estimates.The weak coercivity of the resulting bilinear forms are proven via the Garding inequality,based on which we prove the optimal-order convergence estimates of the finite element method for the(adjoint)state variable and the control variable.Numerical experiments substantiate the theoretical predictions.

Influence of dissipation on solitary wave solution to generalized Boussinesq equation

This paper uses the theory of planar dynamic systems and the knowledge of reaction-diffusion equations,and then studies the bounded traveling wave solution of the generalized Boussinesq equation affected by dissipation and the influence of dissipation on solitary waves.The dynamic system corresponding to the traveling wave solution of the equation is qualitatively analyzed in detail.The influence of the dissipation coefficient on the solution behavior of the bounded traveling wave is studied,and the critical values that can describe the magnitude of the dissipation effect are,respectively,found for the two cases of b_3<0 and b_3>0 in the equation.The results show that,when the dissipation effect is significant(i.e.,r is greater than the critical value in a certain situation),the traveling wave solution to the generalized Boussinesq equation appears as a kink-shaped solitary wave solution;when the dissipation effect is small(i.e.,r is smaller than the critical value in a certain situation),the traveling wave solution to the equation appears as the oscillation attenuation solution.By using the hypothesis undetermined method,all possible solitary wave solutions to the equation when there is no dissipation effect(i.e.,r=0)and the partial kink-shaped solitary wave solution when the dissipation effect is significant are obtained;in particular,when the dissipation effect is small,an approximate solution of the oscillation attenuation solution can be achieved.This paper is further based on the idea of the homogenization principles.By establishing an integral equation reflecting the relationship between the approximate solution of the oscillation attenuation solution and the exact solution obtained in the paper,and by investigating the asymptotic behavior of the solution at infinity,the error estimate between the approximate solution of the oscillation attenuation solution and the exact solution is obtained,which is an infinitesimal amount that decays exponentially.The influence of the dissipation coefficient on the amplitude,frequency,period,and energy of the bounded traveling wave solution of the equation is also discussed.

Virtual Element Discretization of Optimal Control Problem Governed by Brinkman Equations

In this paper, we discuss virtual element method (VEM) approximation of optimal control problem governed by Brinkman equations with control constraints. Based on the polynomial projections and variational discretization of the control variable, we build up the virtual element discrete scheme of the optimal control problem and derive the discrete first order optimality system. A priori error estimates for the state, adjoint state and control variables in L2 and H1 norm are derived. The theoretical findings are illustrated by the numerical experiments.

Linearization of the Liouville Equation, Multiple Splits of the Chandler Frequency, Markowitz Wobbles, and Error Analysis

The rotation of the physical Earth is far more complex than the rotation of a biaxial or slightly triaxial rigid body can represent. The linearization of the Liouville equation via the Munk and MacDonal perturbation scheme has oversimplified polar excitation physics. A more conventional linearization of the Liouville equation as the generalized equation of motion for free rotation of the physical Earth reveals: 1) The reference frame is most essential, which needs to be unique and physically located in the Earth;2) Physical angular momentum perturbation arises from motion and mass redistribution to appear as relative angular momentum in a rotating Earth, which excites polar motion and length of day variations;3) At polar excitation, the direction of the rotation axis in space does not change besides nutation and precession around the invariant angular momentum axis, while the principal axes shift responding only to mass redistribution;4) Two inertia changes appear simultaneously at polar excitation;one is due to mass redistribution, and the other arises from the axial near-symmetry of the perturbed Earth;5) The Earth at polar excitation becomes slightly triaxial and axially near-symmetrical even it was originally biaxial;6) At polar excitation, the rotation of a non-rigid Earth becomes unstable;7) The instantaneous figure axis or mean excitation axis around which the rotation axis physically wobbles is not a principal axis;8) In addition to amplitude excitation, the Chandler wobble possesses also multiple frequency-splits and is slow damping;9) Secular polar drift is after the products of inertia and always associated with the Chandler wobble;both belong to polar motion;10) The Earth will reach its stable rotation only after its rotation axis, major principal axis, and instantaneous figure axis or mean excitation axis are all completely aligned with each other to arrive at the minimum energy configuration of the system;11) The observation of the multiple splits of the Chandler frequency is further examined by means of exact-bandwidth filtering and spectral analysis, which confirms the theoretical prediction of the linearized Liouville equation. After the removal of the Gibbs phenomenon from the polar motion spectra, Markowitz wobbles are also observed;12) Error analysis of the ILS data demonstrates that the incoherent noises from the Wars in 1920-1945 are separable from polar motion and removable, so the ILS data are still reliable and useful for the study of the continuation of polar motion.

A Study Regarding Solution of a Knowledge Model Based on the Containing-Type Error Matrix Equation

An error matrix equation based on error matrix theory was presented in previous research of the error-eliminating theory. The purpose of solving the error matrix equation is to create a decision support on how to switch from bad to good status. A matrix based on error logic is used to express current status u, expectant status u1 and transformation matrix T. It is u, u1, and T that are used to build error matrix equation T (u)= u1. This allows us to find a method whereby bad status “u” changes to good status “u1” by solving the equation. The conversion method that transform from current to expectant status can be obtained from the transformation matrix T. On this basis, this paper proposes a new kind of error matrix equation named “containing-type error matrix equation”. This equation is more suitable for analyzing the realistic question. The method of solving, existence and form of solution for this type of equation have been presented in this paper. This research provides a potential useful new technique for decision analysis.

Linear Error Equation on Field F_2

在这份报纸，我们将给一个方法在 F2 上解决线性错误 equationon，由使用地 F2 和分区理论上的线性代数学。

A Priori and A Posteriori Error Estimates of Streamline Diffusion Finite Element Method for Optimal Control Problem Governed by Convection Dominated Diffusion Equation

In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existence and uniqueness of the discretized scheme.Then a priori and a posteriori error estimates are derived for the state,the co-state and the control.Three numerical examples are presented to illustrate our theoretical results.

ERROR ESTIMATES OF GALERKIN METHOD FOR HIGH DIMENSIONAL STEADY STATE KURAMOTO-SIVASHINSKY EQUATION

Error estimates of Galerkin method for Kuramoto-Sivashingsky (K-S) equation in space dimension ≥3 are derived in the paper. These results furnish strong evidence for the computation of the solutions.

PARALLEL SCHWARZ ALGORITHMS FOR PARABOLIC EQUATIONS AND ERROR ESTIMATES

In this paper we introduce two kinds of parallel Schwarz domain decomposition me thods for general, selfadjoint, second order parabolic equations and study the dependence of their convergence rates on parameters of time-step and space-mesh. We prove that the, approximate solution has convergence independent of iteration times at each time-level. And the L^2 error estimates are given.

A NONLINEAR GALERKIN MIXED ELEMENT METHOD AND A POSTERIORI ERROR ESTIMATOR FOR THE STATIONARY NAVIER-STOKES EQUATIONS

A nonlinear Galerkin mixed element (NGME) method and a posteriori error exstimator based on the method are established for the stationary Navier-Stokes equations. The existence and error estimates of the NGME solution are first discussed, and then a posteriori error estimator based on the NGME method is derived.

OPTIMAL ERROR BOUND IN A SOBOLEV SPACE OF REGULARIZED APPROXIMATION SOLUTIONS FOR A SIDEWAYS PARABOLIC EQUATION

The inverse heat conduction problem （IHCP） is a severely ill-posed problem in the sense that the solution （ if it exists） does not depend continuously on the data. But now the results on inverse heat conduction problem are mainly devoted to the standard inverse heat conduction problem. Some optimal error bounds in a Sobolev space of regularized approximation solutions for a sideways parabolic equation, i. e. , a non-standard inverse heat conduction problem with convection term which appears in some applied subject are given.

THE LARGE TIME ERROR ESTIMATES OF FOURIER SPECTRAL METHOD FOR GENERALIZED BENJAMIN-BONA-MAHONY EQUATIONS

In this paper, a spectral method to analyze the generalized Benjamin Bona Mahony equations is used. The existence and uniqueness of global smooth solution of these equations are proved. The large time error estimation between the spectral approximate solution and the exact solution is obtained.

A COUPLING METHOD OF DIFFERENCE WITH HIGH ORDER ACCURACY AND BOUNDARY INTEGRAL EQUATION FOR EVOLUTIONARY EQUATION AND ITS ERROR ESTIMATES

In the present paper, a new numerical method for solving initial-boundary value problems of evolutionary equations is proposed and studied, combining difference method with high accuracy with boundary integral equation method. The numerical approximate schemes for both problems on a bounded or unbounded domain in R3 are proposed and their prior error estimates are obtained.

A POSTERIORI ERROR ESTIMATE OF THE DSD METHOD FOR FIRST-ORDER HYPERBOLIC EQUATIONS

A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illustrate the accuracy and feasibility of this method.

RESIDUAL A POSTERIORI ERROR ESTIMATE TWO-GRID METHODS FOR THE STEADY (NAVIER-STOKES) EQUATION WITH STREAM FUNCTION FORM

Residual based on a posteriori error estimates for conforming finite element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level method were derived. The posteriori error estimates contained additional terms in comparison to the error estimates for the solution obtained by the standard finite element method. The importance of these additional terms in the error estimates was investigated by studying their asymptotic behavior. For optimal scaled meshes, these bounds are not of higher order than of convergence of discrete solution.

Optimal error estimates for Fourier spectral approximation of the generalized KdV equation

A Fourier spectral method for the generalized Korteweg-de Vries equation with periodic boundary conditions is analyzed, and a corresponding optimal error estimate in L^2-norm is obtained. It improves the result presented by Maday and Quarteroni. A modified Fourier pseudospectral method is also presented, with the same convergence properties as the Fourier spectral method.