A New Linearized Crank-Nicolson Mixed Element Scheme for the Extended Fisher-Kolmogorov Equation

We present a new mixed finite element method for solving the extended Fisher-Kolmogorov (EFK) equation. We first decompose the EFK equation as the two second-order equations, then deal with a second-order equation employing finite element method, and handle the other second-order equation using a new mixed finite element method. In the new mixed finite element method, the gradient ∇u belongs to the weaker (L2(Ω))2 space taking the place of the classical H(div;Ω) space. We prove some a priori bounds for the solution for semidiscrete scheme and derive a fully discrete mixed scheme based on a linearized Crank-Nicolson method. At the same time, we get the optimal a priori error estimates in L2 and H1-norm for both the scalar unknown u and the diffusion term w=−Δu and a priori error estimates in (L2)2-norm for its gradient χ=∇u for both semi-discrete and fully discrete schemes

A Time Two-Mesh Compact Difference Method for the One-Dimensional Nonlinear Schrödinger Equation

The nonlinear Schrödinger equation is an important model equation in the study of quantum states of physical systems. To improve the computing efficiency, a fast algorithm based on the time two-mesh high-order compact difference scheme for solving the nonlinear Schrödinger equation is studied. The fourth-order compact difference scheme is used to approximate the spatial derivatives and the time two-mesh method is designed for efficiently solving the resulting nonlinear system. Comparing to the existing time two-mesh algorithm, the novelty of the new algorithm is that the fine mesh solution, which becomes available, is also used as the initial guess of the linear system, which can improve the calculation accuracy of fine mesh solutions. Compared to the two-grid finite element methods (or finite difference methods) for nonlinear Schrödinger equations, the numerical calculation of this method is relatively simple, and its two-mesh algorithm is implemented in the temporal direction. Taking advantage of the discrete energy, the result with O(τC4+τF2+h4) in the discrete L2-norm is obtained. Here, τC and τF are the temporal parameters on the coarse and fine mesh, respectively, and h is the space step size. Finally, some numerical experiments are conducted to demonstrate its efficiency and accuracy. The numerical results show that the new algorithm gives highly accurate results and preserves conservation laws of charge and energy. Furthermore, by comparing with the standard nonlinear implicit compact difference scheme, it can reduce the CPU time without loss of accuracy

Nonlinear Vibrations of a Rotor-Active Magnetic Bearing System with 16-Pole Legs and Two Degrees of Freedom

The asymptotic perturbation method is used to analyze the nonlinear vibrations and chaotic dynamics of a rotor-active magnetic bearing (AMB) system with 16-pole legs and the time-varying stiffness. Based on the expressions of the electromagnetic force resultants, the influences of some parameters, such as the cross-sectional area Aα of one electromagnet and the number N of windings in each electromagnet coil, on the electromagnetic force resultants are considered for the rotor-AMB system with 16-pole legs. Based on the Newton law, the governing equation of motion for the rotor-AMB system with 16-pole legs is obtained and expressed as a two-degree-of-freedom system with the parametric excitation and the quadratic and cubic nonlinearities. According to the asymptotic perturbation method, the four-dimensional averaged equation of the rotor-AMB system is derived under the case of 1 : 1 internal resonance and 1 : 2 subharmonic resonances. Then, the frequency-response curves are employed to study the steady-state solutions of the modal amplitudes. From the analysis of the frequency responses, both the hardening-type nonlinearity and the softening-type nonlinearity are observed in the rotor-AMB system. Based on the numerical solutions of the averaged equation, the changed procedure of the nonlinear dynamic behaviors of the rotor-AMB system with the control parameter is described by the bifurcation diagram. From the numerical simulations, the periodic, quasiperiodic, and chaotic motions are observed in the rotor-active magnetic bearing (AMB) system with 16-pole legs, the time-varying stiffness, and the quadratic and cubic nonlinearities

An expanded mixed covolume element method for integro-differential equation of Sobolev type on triangular grids

Abstract The expanded mixed covolume Element (EMCVE) method is studied for the two-dimensional integro-differential equation of Sobolev type. We use a piecewise constant function space and the lowest order Raviart-Thomas ( RT 0 $\mathit{RT}_{0}$ ) space as the trial function spaces of the scalar unknown u and its gradient σ and flux λ, respectively. The semi-discrete and backward Euler fully-discrete EMCVE schemes are constructed, and the optimal a priori error estimates are derived. Moreover, numerical results are given to verify the theoretical analysis

A Time Two-Mesh Finite Difference Numerical Scheme for the Symmetric Regularized Long Wave Equation

The symmetric regularized long wave (SRLW) equation is a mathematical model used in many areas of physics; the solution of the SRLW equation can accurately describe the behavior of long waves in shallow water. To approximate the solutions of the equation, a time two-mesh (TT-M) decoupled finite difference numerical scheme is proposed in this paper to improve the efficiency of solving the SRLW equation. Based on the time two-mesh technique and two time-level finite difference method, the proposed scheme can calculate the velocity u(x,t) and density ρ(x,t) in the SRLW equation simultaneously. The linearization process involves a modification similar to the Gauss-Seidel method used for linear systems to improve the accuracy of the calculation to obtain solutions. By using the discrete energy and mathematical induction methods, the convergence results with O(τC2+τF+h2) in the discrete L∞-norm for u(x,t) and in the discrete L2-norm for ρ(x,t) are proved, respectively. The stability of the scheme was also analyzed. Finally, some numerical examples, including error estimate, computational time and preservation of conservation laws, are given to verify the efficiency of the scheme. The numerical results show that the new method preserves conservation laws of four quantities successfully. Furthermore, by comparing with the original two-level nonlinear finite difference scheme, the proposed scheme can save the CPU time

A New Mixed Element Method for a Class of Time-Fractional Partial Differential Equations

A kind of new mixed element method for time-fractional partial differential equations is studied. The Caputo-fractional derivative of time direction is approximated by two-step difference method and the spatial direction is discretized by a new mixed element method, whose gradient belongs to the simple L2Ω2 space replacing the complex H(div;Ω) space. Some a priori error estimates in L2-norm for the scalar unknown u and in L22-norm for its gradient σ. Moreover, we also discuss a priori error estimates in H1-norm for the scalar unknown u