Abstract
We present an efficient second-order accurate scheme to treat stiff source terms within the framework of higher order Godunov’s methods. We employ Duhamel’s formula to devise a modified predictor step which accounts for the effects of stiff source terms on the conservative fluxes and recovers the correct isothermal behavior in the limit of an infinite cooling/reaction rate. Source term effects on the conservative quantities are fully accounted for by means of a one-step, second-order accurate semi-implicit corrector scheme based on the deferred correction method of Dutt et al. We demonstrate the accurate, stable and convergent results of the proposed method through a set of benchmark problems for a variety of stiffness conditions and source types.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
