2,939 research outputs found

    An accurate block solver for stiff initial value problems

    Get PDF
    New implicit block formulae that compute solution of stiff initial value problems at two points simultaneously are derived and implemented in a variable step size mode. The strategy for changing the step size for optimum performance involves halving, increasing by a multiple of 1.7, or maintaining the current step size. The stability analysis of the methods indicates their suitability for solving stiff problems. Numerical results are given and compared with some existing backward differentiation formula algorithms. The results indicate an improvement in terms of accuracy

    A new variable step size block backward differentiation formula for solving stiff initial value problems

    Get PDF
    A new block backward differentiation formula of order 4 with variable step size is formulated. By varying a parameter in the formula, different sets of formulae with A-stability property can be generated. At the cost of an additional function evaluation, the accuracy of the method is seen to outperform some existing backward differentiation formula algorithms. The strategy involved in controlling the step size ratio is also described. The problems tested with the method show its efficiency in solving stiff initial value problems

    Numerical solution for stiff initial value problems using 2-point block multistep method

    Get PDF
    This paper focuses on the derivation of an improved 2-point Block Backward Differentiation Formula of order five (I2BBDF(5)) for solving stiff first order Ordinary Differential Equations (ODEs). The I2BBDF(5) method is derived by using Taylor's series expansion to obtain the coefficients of the formula. To verify the efficiency of the I2BBDF(5) method, stiff problems from the literature are tested and compared with the existing solver for stiff ODEs. From the numerical results, we conclude that the I2BBDF(5) method can be an alternative solver for solving stiff ODE

    An implicit 2-point block extended backward differentiation formula for integration of stiff initial value problems

    Get PDF
    A class of implicit 2-point block extended backward differentiation formula (BEBDF) of order 4 is presented. The stability region of the method is constructed and shown to be A – stable. Results obtained are compared with an existing block backward differentiation formula (BBDF). The comparison shows that using constant step size and the same number of integration steps, our method achieves greater accuracy than the 2-point BBDF and is suitable for solving stiff initial value problems

    Fourth-order time-stepping for stiff PDEs on the sphere

    Full text link
    We present in this paper algorithms for solving stiff PDEs on the unit sphere with spectral accuracy in space and fourth-order accuracy in time. These are based on a variant of the double Fourier sphere method in coefficient space with multiplication matrices that differ from the usual ones, and implicit-explicit time-stepping schemes. Operating in coefficient space with these new matrices allows one to use a sparse direct solver, avoids the coordinate singularity and maintains smoothness at the poles, while implicit-explicit schemes circumvent severe restrictions on the time-steps due to stiffness. A comparison is made against exponential integrators and it is found that implicit-explicit schemes perform best. Implementations in MATLAB and Chebfun make it possible to compute the solution of many PDEs to high accuracy in a very convenient fashion

    Block Milne’s Implementation For Solving Fourth Order Ordinary Differential Equations

    Get PDF
    Block predictor-corrector method for solving non-stiff ordinary differential equations (ODEs) started with Milne’s device. Milne’s device is an extension of the block predictor corrector method providing further benefits and better results. This study considers Milne’s devise for solving fourth order ODEs. A combination of Newton’s backward difference interpolation polynomial and numerical integration method are applied and integrated at some selected grid points to formulate the block predictor-corrector method. Moreover, Milne’s devise advances the computational efficiency by applying the chief local truncation error] (CLTE) of the block predictor-corrector method after establishing the order. The numerical results were exhibited to attest the functioning of Milne’s devise in solving fourth order ODEs. The complete results were obtained with the aid of Mathematica 9 kernel for Microsoft Windows. Numerical results showcase that Milne’s device is more effective than existent methods in terms of design new step size, determining the convergence criteria and maximizing errors at all examined convergence levels

    Emergent properties of the G1/S network

    Get PDF
    Tato práce se zabývá buněčným cyklem kvasinky Saccgaromyces cerevisiae. Oblastí našeho zájmu je přechod mezi G1 a S fází, kde je naším cílem identifikovat velikosti buňky v době počátku DNA replikace. Nejprve se věnujeme nedávno publikovanému matematickému modelu, který popisuje mechanismy vedoucí k S fázi. Práce poskytuje detailní popis tohoto modelu, stejně jako časový průběh některých důležitých proteinů či jejich sloučenin. Dále se zabýváme pravděpodobnostním modelem aktivace replikačních počátků DNA. Nově uvažujeme vliv šíření DNA replikace mezi sousedícími počátky a analyzujeme jeho důsledky. Poskytujeme také senzitivní analýzu kritické velikosti buňky vzhledem ke konstantám popisujícím dynamiku reakcí v modelu G1/S přechodu.In this thesis we deal with the cell cycle of the yeast Saccharomyces cerevisiae. We are interested in its G1 to S transition, and our main goal is to determine the cell size at the onset of its DNA replication. At first, we study a recent mathematical model describing the mechanisms leading to the S phase, we provide its detailed description and present the dynamics of some significant protein and protein complexes. Further, we take a closer look at the probabilistic model for firing of DNA replication origins. We newly consider the influence of DNA replication spreading among neighboring origins, and we analyze its consequences. We also provide a sensitivity analysis of the critical cell size with respect to rate constants of G1 to S transition model.

    Partitioning Techniques and Their Parallelization for Stiff System of Ordinary Differential Equations

    Get PDF
    A new code based on variable order and variable stepsize component wise partitioning is introduced to solve a system of equations dynamically. In previous partitioning technique researches, once an equation is identified as stiff, it will remain in stiff subsystem until the integration is completed. In this current technique, the system is treated as nonstiff and any equation that caused stiffness will be treated as stiff equation. However, should the characteristics showed the elements of nonstiffness, and then it will be treated again with Adam method. This process will continue switching from stiff to nonstiff vice versa whenever it is necessary until the interval of integration is completed.Next, a block method with R-points generate R new approximate solution values;is a strategy for solving a system and also for parallelizing ODEs. Partitioning this block method to solve stiff differential equations is a new strategy; it is more efficient and takes less computational time compared to the sequential methods. Two partitioning techniques are constructed, Intervalwise Block Partitioning (IBP) and Componentwise Block Partitioning (CBP). Numerical results are compared as validation of its effectiveness. Intervalwise block partitioning will initially treat the systems of equations as nonstiff and solve them using Adams method, by switching to the Backward Differentiation formula when there is a step failure and indication of stiffness. Componentwise block partitioning will place the necessary equations that cause instability and stiffness into the stiff subsystem and solve using Backward Differentiation Formula, while all other equations will still be treated as non-stiff and solved using Adams formula. Parallelizing the partitioning strategies using Message Passing Interface (MPI) is the most appropriate method to solve large system of equations. Parallelizing the right algorithm in the partitioning code will give a better perfonnance with shorter execution times. The graphs of its performance and execution time, visualize the advantages of parallelizing

    A new fifth order implicit block method for solving first order stiff ordinary differential equations

    Get PDF
    A new implicit block backward differentiation formula that computes 3–points simultaneously is derived. The method is of order 5 and solves system of stiff ordinary differential equations (ODEs). The stability analysis indicates that the method is A–stable. Numerical results show that the method outperformed some existing block and non-block methodsfor solving stiff ODEs

    Programming codes of block-Milne's device for solving fourth-order ODEs

    Get PDF
    Block-Milne’s device is an extension of block-predictor-corrector method and specifically developed to design a worthy step size, resolve the convergence criteria and maximize error. In this study, programming codes of block- Milne’s device (P-CB-MD) for solving fourth order ODEs are considered. Collocation and interpolation with power series as the basic solution are used to devise P-CB-MD. Analysing the P-CB-MD will give rise to the principal local truncation error (PLTE) after determining the order. The P-CB-MD for solving fourth order ODEs is written using Mathematica which can be utilized to evaluate and produce the mathematical results. The P-CB-MD is very useful to demonstrate speed, efficiency and accuracy compare to manual computation applied. Some selected problems were solved and compared with existing methods. This was made realizable with the support of the named computational benefit
    corecore