It is shown that the conforming Q2,1.1.2-Q'1.mixed element is stable, and provides optimal order of approximation for the Stokes equations on rectangular grids. Here, Q2,1.1.2 = Q2,1.1.; Q1.2, and Q2,1.denotes ...
详细信息
It is shown that the conforming Q2,1.1.2-Q'1.mixed element is stable, and provides optimal order of approximation for the Stokes equations on rectangular grids. Here, Q2,1.1.2 = Q2,1.× Q1.2, and Q2,1.denotes the space of continuous piecewise-polynomials of degree 2 or less in the x direction but of degree 1.in the y direction. Q'1.is the space of discontinuous bilinear polynomials, with spurious modes filtered. To be precise, Q'1.is the divergence of the discrete velocity space Q2,1.1.2. Therefore, the resulting finite element solution for the velocity is divergence-free pointwise, when solving the Stokes equations. This element is the lowest order one in a family of divergence-free element, similar to the families of the Bernardi-Raugel element and the RaviartThomas element.
In this paper, we present an arbitrarily high-order numerical scheme for the two-component Camassa–Holm system, ensuring the preservation of three invariants: energy and two Casimir functions. The spatial discretizat...
In this paper, we present an arbitrarily high-order numerical scheme for the two-component Camassa–Holm system, ensuring the preservation of three invariants: energy and two Casimir functions. The spatial discretization is achieved using Fourier–Galerkin methods, resulting in a semi-discrete system which retains a Hamiltonian structure and approximates the invariants of the original continuous system. Subsequently, an energy-preserving Runge–Kutta method, such as Hamiltonian boundary value methods, is applied to the semi-discrete system, yielding a fully discrete scheme of arbitrarily high order. The proposed scheme is validated through numerical simulations, demonstrating its accuracy and capability in capturing different types of solutions of the two-component Camassa–Holm equation.
In this paper, we provide a theoretical method(PUFEM), which belongs to the analysis of the partition of unity finite element family of meshfree methods. The usual error analysis only shows the order of error estima...
详细信息
In this paper, we provide a theoretical method(PUFEM), which belongs to the analysis of the partition of unity finite element family of meshfree methods. The usual error analysis only shows the order of error estimate to the same as the local approximations[1.]. Using standard linear finite element base functions as partition of unity and polynomials as local approximation space, in l-d case, we derive optimal order error estimates for PUFEM interpolants. Our analysis show that the error estimate is of one order higher than the local approximations. The interpolation error estimates yield optimal error estimates for PUFEM solutions of elliptic boundary value problems.
In this paper,we investigate the dependence of the solutions on the parameters(order,initial function,right-hand function)of fractional delay differential equations(FDDEs)with the Caputo fractional *** results includi...
详细信息
In this paper,we investigate the dependence of the solutions on the parameters(order,initial function,right-hand function)of fractional delay differential equations(FDDEs)with the Caputo fractional *** results including an estimate of the solutions of FDDEs are given *** results are verified by some numerical examples.
The method of recovering a low-rank matrix with an unknown fraction whose entries are arbitrarily corrupted is known as the robust principal component analysis (RPCA). This RPCA problem, under some conditions, can b...
详细信息
The method of recovering a low-rank matrix with an unknown fraction whose entries are arbitrarily corrupted is known as the robust principal component analysis (RPCA). This RPCA problem, under some conditions, can be exactly solved via convex optimization by minimizing a combination of the nuclear norm and the 1. norm. In this paper, an algorithm based on the Douglas-Rachford splitting method is proposed for solving the RPCA problem. First, the convex optimization problem is solved by canceling the constraint of the variables, and ~hen the proximity operators of the objective function are computed alternately. The new algorithm can exactly recover the low-rank and sparse components simultaneously, and it is proved to be convergent. Numerical simulations demonstrate the practical utility of the proposed algorithm.
We present and analyze a robust preconditioned conjugate gradient method for the higher order Lagrangian finite element systems of a class of elliptic problems. An auxiliary linear element stiffness matrix is chosen t...
详细信息
We present and analyze a robust preconditioned conjugate gradient method for the higher order Lagrangian finite element systems of a class of elliptic problems. An auxiliary linear element stiffness matrix is chosen to be the preconditioner for higher order finite elements. Then an algebraic multigrid method of linear finite element is applied for solving the preconditioner. The optimal condition number which is independent of the mesh size is obtained. Numerical experiments confirm the efficiency of the algorithm.
In this paper,we apply an a posteriori error control theory that we develop in a very recent paper to three families of the discontinuous Galerkin methods for the Reissner-Mindlin plate *** derive robust a posteriori ...
详细信息
In this paper,we apply an a posteriori error control theory that we develop in a very recent paper to three families of the discontinuous Galerkin methods for the Reissner-Mindlin plate *** derive robust a posteriori error estimators for them and prove their reliability and efficiency.
This paper considers a class of discontinuous Galerkin method,which is constructed by Wong-Zakai approximation with the orthonormal Fourier basis,for numerically solving nonautonomous Stratonovich stochastic delay dif...
详细信息
This paper considers a class of discontinuous Galerkin method,which is constructed by Wong-Zakai approximation with the orthonormal Fourier basis,for numerically solving nonautonomous Stratonovich stochastic delay differential *** prove that the discontinuous Galerkin scheme is strongly convergent:globally stable and analogously asymptotically stable in mean square *** addition,this method can be easily extended to solve nonautonomous Stratonovich stochastic pantograph differential *** tests indicate that the method has first-order and half-order strong mean square convergence,when the diffusion term is without delay and with delay,respectively.
In this paper,we study the a posteriori error estimator of SDG method for variable coefficients time-harmonic Maxwell's *** propose two a posteriori error estimators,one is the recovery-type estimator,and the othe...
详细信息
In this paper,we study the a posteriori error estimator of SDG method for variable coefficients time-harmonic Maxwell's *** propose two a posteriori error estimators,one is the recovery-type estimator,and the other is the residual-type *** first propose the curl-recovery method for the staggered discontinuous Galerkin method(SDGM),and based on the super-convergence result of the postprocessed solution,an asymptotically exact error estimator is *** residual-type a posteriori error estimator is also proposed,and it's reliability and effectiveness are proved for variable coefficients time-harmonic Maxwell's *** efficiency and robustness of the proposed estimators is demonstrated by the numerical experiments.
For fractional Volterra integro-differential equations(FVIDEs)with weakly singular kernels,this paper proposes a generalized Jacobi spectral Galerkin *** basis functions for the provided method are selected generalize...
详细信息
For fractional Volterra integro-differential equations(FVIDEs)with weakly singular kernels,this paper proposes a generalized Jacobi spectral Galerkin *** basis functions for the provided method are selected generalized Jacobi functions(GJFs),which can be utilized as natural basis functions of spectral methods for weakly singular FVIDEs when appropriately *** developed method's spectral rate of convergence is determined using the L^(∞)-norm and the weighted L^(2)-*** results indicate the usefulness of the proposed method.
暂无评论