This paper examines a linear extrapolation time-discretization of a phase relaxation model with temperature dependent convection and reaction. The model consists of a diffusion-advection PDE for temperature and an ODE with double obstacle for phase variable. Under a stability constraint, this scheme is shown to converge with optimal orders for temperature and enthalpy, and for heat flux as time-step .
A family of nonlinear conservative finite difference schemes for the multidimensional Boussinesq Paradigm Equation is considered. A second order of convergence and a preservation of the discrete energy for this approach are proved. Existence and boundedness of the discrete solution on an appropriate time interval are established. The schemes have been numerically tested on the models of the propagation of a soliton and the interaction of two solitons. The numerical experiments demonstrate that the...
An energy conservative scheme is proposed for the regularized long wave (RLW) equation. The integral method with variational limit is used to discretize the spatial derivative and the finite difference method is used to discretize the time derivative. The energy conservation of the scheme and existence of the numerical solution are proved. The convergence of the order and unconditional stability are also derived. Numerical examples are carried out to verify the correctness of the theoretical analysis....
Finite element methods with piecewise polynomial spaces in space for solving the nonstationary heat equation, as a model for parabolic equations are considered. The discretization in time is performed using the Crank-Nicolson method. A new a priori estimate is proved. Thanks to this new a priori estimate, a new error estimate in the discrete norm of is proved. An -error estimate is also shown. These error estimates are useful since they allow us to get second order time accurate approximations...
First–order accurate monotone conservative schemes have good convergence and stability properties, and thus play a very important role in designing modern high resolution shock-capturing schemes. Do the monotone difference approximations always give a good numerical solution in sense of monotonicity preservation or suppression of oscillations? This note will investigate this problem from a numerical point of view and show that a -point monotone scheme may give an oscillatory solution even though...
First–order accurate monotone conservative schemes have good convergence and stability properties, and thus play a very important role in designing modern high resolution shock-capturing schemes. Do the monotone difference approximations always give a good numerical solution in sense of monotonicity preservation or suppression of oscillations? This note will investigate this problem from a numerical point of view and show that a (2K+1)-point monotone scheme may give an oscillatory solution even...
We prove a posteriori error estimates of optimal order for linear Schrödinger-type equations in the L∞(L2)- and the L∞(H1)-norm. We discretize only in time by the Crank-Nicolson method. The direct use of the reconstruction technique, as it has been proposed by Akrivis et al. in [Math. Comput. 75 (2006) 511–531], leads to a posteriori upper bounds that are of optimal order in the L∞(L2)-norm, but of suboptimal order in the L∞(H1)-norm. The optimality in the case of L∞(H1)-norm is recovered by using...
We prove a posteriori error estimates of optimal order for linear Schrödinger-type equations in the L∞(L2)- and the L∞(H1)-norm. We discretize only in time by the Crank-Nicolson method. The direct use of the reconstruction technique, as it has been proposed by Akrivis et al. in [Math. Comput.75 (2006) 511–531], leads to a posteriori upper bounds that are of optimal order in the L∞(L2)-norm, but of suboptimal order in the L∞(H1)-norm. The optimality in the case of L∞(H1)-norm is recovered by using...
The paper presents an a posteriori error estimator for a (piecewise linear) nonconforming finite element approximation of the heat equation in , or 3, using backward Euler’s scheme. For this discretization, we derive a residual indicator, which use a spatial residual indicator based on the jumps of normal and tangential derivatives of the nonconforming approximation and a time residual indicator based on the jump of broken gradients at each time step. Lower and upper bounds form the main results...