The purpose of this work is to introduce a new idea how avoid the factorization large matrices during solution stiff systems ODEs. Starting from general form an explicit linear multistep method we suggest adaptively choose its coefficients on each integration step in order minimize norm residual implicit BDF formula. Thereby reduce number unknowns n O(1), where dimension ODE system. We call typ...

This work presents two high-order, semi-discrete, central-upwind schemes for computing approximate solutions of 1D systems of conservation laws. We propose a central weighted essentially non-oscillatory (CWENO) reconstruction, also we apply a fourth-order reconstruction proposed by Peer et al., and afterwards, we combine these reconstructions with a semi-discrete central-upwind numerical flux ...

Abstract In the context of coupling hyperbolic problems, maximum stable time step an explicit numerical scheme may depend on design procedure. If this is case, procedure sensitive to changes in model parameters independent Courant–Friedrichs–Levy condition. This sensitivity can cause artificial stiffness that degrades performance a scheme. To overcome problem, we present systematic and general ...

The spatial discretization of the magnetic vector potential formulation of magnetoquasistatic field problems results in an infinitely stiff differential-algebraic equation system. It is transformed into a finitely stiff ordinary differential equation system by applying a generalized Schur complement. Applying the explicit Euler time integration scheme to this system results in a small maximum s...

In this paper we propose a class of second derivative multistep methods for solving some well-known classes of LaneEmden type equations which are nonlinear ordinary differential equations on the semi-infinite domain. These methods, which have good stability and accuracy properties, are useful in deal with stiff ODEs. We show superiority of these methods by applying them on the some famous Lane-...

For reaction-diffusion-advection equations, the stiffness from the reaction and diffusion terms often requires very restricted time step size, while the nonlinear advection term may lead to a sharp gradient in localized spatial regions. It is challenging to design numerical methods that can efficiently handle both difficulties. For reaction-diffusion systems with both stiff reaction and diffusi...