[{"data":1,"prerenderedAt":-1},["ShallowReactive",2],{"doc-detail-83254-en":3,"doc-seo-83254-105":30,"detail-sidebar-cat-0-en-105":91},{"code":4,"msg":5,"data":6},0,"success",{"doc_id":7,"user_id":8,"nickname":9,"user_avatar":10,"doc_module":4,"category_id":11,"category_name":12,"doc_title":13,"doc_description":14,"doc_content":15,"file_id":16,"file_url":17,"file_type":18,"file_size":19,"view_count":20,"is_deleted":4,"is_public":21,"is_downloadable":21,"audit_status":21,"page_count":22,"language":23,"language_code":24,"site_id":25,"html_lang":24,"table_of_contents":26,"faqs":27,"seo_title":13,"seo_description":14,"update_tm":28,"read_time":29},83254,962075114101,"Seraphina","https://ap-avatar.wpscdn.com/avatar/e000253a75eb197efd?x-image-process=image/resize,m_fixed,w_180,h_180&k=1780044092746381165",8,"Research & Report","Explicit stabilized implementation of singly diagonally implicit Runge-Kutta methods","Implicit Runge–Kutta methods provide an effective way to integrate stiff differential equations by avoiding restrictive time-step limits of standard explicit schemes. Explicit stabilized integrators offer an alternative, especially for high-dimensional advection–diffusion–reaction problems with diffusive terms. A new explicit stabilized implementation for a class of diagonally implicit Runge–Kutta methods is presented, enabling high-order singly-diagonally-implicit schemes for PDEs with a computational cost analogous to classical explicit stabilized methods. The method rewrites the implicit update as a steady state of a modified auxiliary system, solved via a partitioned Runge–Kutta–Chebyshev procedure inspired by optimization.","arXiv :2607 .07497v1 [math .NA] 8 Jul 2026  \nExplicit stabilized implementation of singly diagonally implicit Runge-Kutta methods  \nIbrahim Almuslimani∗, Gilles Vilmart†, and Konstantinos Zygalakis‡  \nJuly 9, 2026  \nAbstract  \nImplicit methods are a natural approach for the integration of stiff differential equations, to avoid time-step restrictions faced by standard explicit integrators. Explicit stabilised integrators are an alternative to implicit methods, which can be particularly efficient in highdimensional applications with diffusive terms. Towards the best of both worlds, we introduce a new explicit stabilised implementation of a class of diagonally implicit Runge-Kutta methods. This allows us to implement high-order singly-diagonally-implicit Runge-Kutta methods for advection-diffusion-reaction PDEs with a provable computational cost analogous to that of standard explicit stabilised methods. The main ingredient is to recast the implicit Runge– Kutta update as the steady state of a modified auxiliary system, which is then computed using a partitioned Runge–Kutta–Chebyshev method inspired by optimisation techniques.  \nKeywords: diagonally implicit Runge-Kutta methods, explicit stabilised methods, partitioned Runge-Kutta methods, stiff problems, advection-diffusion-reaction problems.  \nAMS subject classification (2010): 65L20, 65M12, 65M20  \n1 Introduction  \nIn this paper, we consider systems of ordinary differential equations (ODEs) representing space discretisations of partial differential equations (PDEs) of the form  \ny˙ = F (y), (1a)  \nF (y) := FD (y) + FA (y) + FR (y) . (1b)  \nwhere FD (y), FA (y), FR (y) 2 Rd represent diffusion, advection and reaction terms respectively. Typically (1) arises from the spatial discretisation of an advection-diffusion-reaction problem in N spatial dimensions, and assuming a spatial grid of size ∆x, the dimension of the ODE grows like d = O((∆x)−N ) .  \nThe numerical time integration of (1) is very challenging since this is a stiff ODE. Following [13], we call a differential equation stiff when a standard explicit numerical integrator such asa standard explicit Runge-Kutta method faces a severe time-step restriction. In particular, in the case of (1) the eigenvalues of the Jacobian of FD are typically distributed along the negative  \n∗ École Polytechnique Fédérale de Lausanne (EPFL), Swiss Plasma Center (SPC), CH-1015 Lausanne, Switzer[land. Ibrahim.Almuslimani@epfl.ch](land. Ibrahim.Almuslimani@epfl.ch)  \n†Section de Mathématiques, Université de Genève, CP 64, 1211 Genève 4, [Switzerland. Gilles.Vilmart@unige.ch](Switzerland. Gilles.Vilmart@unige.ch)  \n‡School of Mathematics and the Maxwell Institute for Mathematical Sciences, University of Edinburgh, Edinburgh, EH9 3FD, [UK. K.Zygalakis@ed.ac.uk](UK. K.Zygalakis@ed.ac.uk)  \nreal axis in an interval that grows as [􀀀O(∆x−2) , 0] (for a symmetric diffusion operator), which implies that the time step ∆t that can be used by the explicit Euler method suffers from the CFL ∆t 􀀔 C∆x2 . Furthermore, the eigenvalues of the Jacobian of FA are typically distributed along the imaginary axis in an interval that grows as [􀀀iO(∆x−1) , iO(∆x−1)], while the eigenvalues of the Jacobian of the reaction term FR even though are usually not related to ∆x can have a ratio maxj jRλj j/ minj jRλj j that can be very large or vary over several orders of magnitude when modeling multiscale reaction terms.  \nThere are two different approaches to dealing with the issue of stiffness within the class of Runge-Kutta methods. The first relates to the design of implicit methods, with a prime example being the implicit Euler method. These methods have the advantage of achieving stability unconditionally with respect to the stepsize, which comes with the cost of solving a system of nonlinear equations. When the dimension of the corresponding differential equation is low, the nonlinear system can be solved efficiently using Newton’s method [17] . As the dimension of the dif","cbCaiiY6gOQ8q67Y","https://ap.wps.com/l/cbCaiiY6gOQ8q67Y","pdf",771206,2,1,25,"English","en",105,"# Introduction\n## Stiff ODEs from PDE discretization\n## Implicit Runge-Kutta and diagonally implicit structure\n## Explicit stabilized methods and Runge-Kutta-Chebyshev families","[{\"question\":\"What problem motivates explicit stabilized methods in this paper?\",\"answer\":\"The paper targets stiff ODE systems arising from spatial discretizations of advection–diffusion–reaction PDEs, where standard explicit integrators face severe time-step restrictions.\"},{\"question\":\"How does the proposed approach relate implicit updates to an auxiliary steady-state system?\",\"answer\":\"It recasts the implicit Runge–Kutta update as the steady state of a modified auxiliary system, which is then computed numerically.\"},{\"question\":\"What numerical technique is used to compute the auxiliary system’s steady state?\",\"answer\":\"A partitioned Runge–Kutta–Chebyshev method is employed, inspired by optimization techniques, to obtain an explicit stabilized implementation with provable cost.\"}]",1784186280,63,{"code":4,"msg":31,"data":32},"ok",{"site_id":25,"language":24,"slug":33,"title":13,"keywords":34,"description":14,"schema_data":35,"social_meta":86,"head_meta":88,"extra_data":90,"updated_unix":28},"explicit-stabilized-implementation-of-singly-diagonally-implicit-runge-kutta-methods","",{"@graph":36,"@context":85},[37,53,68],{"@type":38,"itemListElement":39},"BreadcrumbList",[40,44,47,50],{"item":41,"name":42,"@type":43,"position":21},"https://docshare.wps.com","Home","ListItem",{"item":45,"name":46,"@type":43,"position":20},"https://docshare.wps.com/document/","Document",{"item":48,"name":12,"@type":43,"position":49},"https://docshare.wps.com/document/research-report/",3,{"item":51,"name":13,"@type":43,"position":52},"https://docshare.wps.com/document/explicit-stabilized-implementation-of-singly-diagonally-implicit-runge-kutta-methods/83254/",4,{"url":51,"name":13,"@type":54,"author":55,"headline":13,"publisher":57,"fileFormat":60,"inLanguage":24,"description":14,"dateModified":61,"datePublished":62,"encodingFormat":60,"isAccessibleForFree":63,"interactionStatistic":64},"DigitalDocument",{"name":9,"@type":56},"Person",{"url":41,"name":58,"@type":59},"DocShare","Organization","application/pdf","2026-07-23","2026-07-16",true,{"@type":65,"interactionType":66,"userInteractionCount":20},"InteractionCounter",{"@type":67},"ViewAction",{"@type":69,"mainEntity":70},"FAQPage",[71,77,81],{"name":72,"@type":73,"acceptedAnswer":74},"What problem motivates explicit stabilized methods in this paper?","Question",{"text":75,"@type":76},"The paper targets stiff ODE systems arising from spatial discretizations of advection–diffusion–reaction PDEs, where standard explicit integrators face severe time-step restrictions.","Answer",{"name":78,"@type":73,"acceptedAnswer":79},"How does the proposed approach relate implicit updates to an auxiliary steady-state system?",{"text":80,"@type":76},"It recasts the implicit Runge–Kutta update as the steady state of a modified auxiliary system, which is then computed numerically.",{"name":82,"@type":73,"acceptedAnswer":83},"What numerical technique is used to compute the auxiliary system’s steady state?",{"text":84,"@type":76},"A partitioned Runge–Kutta–Chebyshev method is employed, inspired by optimization techniques, to obtain an explicit stabilized implementation with provable cost.","https://schema.org",{"og:url":51,"og:type":87,"og:title":13,"og:site_name":58,"og:description":14},"article",{"robots":89,"canonical":51},"index,follow",{"doc_id":7,"site_id":25},{"code":4,"msg":5,"data":92},[93,97,101,105,110,115,120,123,128,131,135],{"id":21,"doc_module":4,"doc_module_name":46,"category_name":94,"show_sort_weight":95,"slug":96},"Story & Novel",90,"story-novel",{"id":20,"doc_module":4,"doc_module_name":46,"category_name":98,"show_sort_weight":99,"slug":100},"Literature",80,"literature",{"id":52,"doc_module":4,"doc_module_name":46,"category_name":102,"show_sort_weight":103,"slug":104},"Exam",70,"exam",{"id":106,"doc_module":4,"doc_module_name":46,"category_name":107,"show_sort_weight":108,"slug":109},5,"Comic",60,"comic",{"id":111,"doc_module":4,"doc_module_name":46,"category_name":112,"show_sort_weight":113,"slug":114},6,"Technology",50,"technology",{"id":116,"doc_module":4,"doc_module_name":46,"category_name":117,"show_sort_weight":118,"slug":119},7,"Healthcare",40,"healthcare",{"id":11,"doc_module":4,"doc_module_name":46,"category_name":12,"show_sort_weight":121,"slug":122},30,"research-report",{"id":124,"doc_module":4,"doc_module_name":46,"category_name":125,"show_sort_weight":126,"slug":127},9,"Religion & Spirituality",20,"religion-spirituality",{"id":126,"doc_module":4,"doc_module_name":46,"category_name":129,"show_sort_weight":126,"slug":130},"World Cup","world-cup",{"id":132,"doc_module":4,"doc_module_name":46,"category_name":133,"show_sort_weight":132,"slug":134},10,"Lifestyle","lifestyle",{"id":136,"doc_module":4,"doc_module_name":46,"category_name":137,"show_sort_weight":106,"slug":138},19,"General","general"]