A Nonoscillatory Second-Order Time-Stepping Procedure for Reaction-Diffusion Equations
After a theory of morphogenesis in chemical cells was introduced in the 1950s, much attention had been devoted to the numerical solution of reaction-diffusion (RD) partial differential equations (PDEs). The Crank–Nicolson (CN) method has been a common second-order time-stepping procedure. However, t...
Saved in:
| Main Authors: | , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Wiley
2020-01-01
|
| Series: | Complexity |
| Online Access: | http://dx.doi.org/10.1155/2020/5163704 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1832567445730099200 |
|---|---|
| author | Philku Lee George V. Popescu Seongjai Kim |
| author_facet | Philku Lee George V. Popescu Seongjai Kim |
| author_sort | Philku Lee |
| collection | DOAJ |
| description | After a theory of morphogenesis in chemical cells was introduced in the 1950s, much attention had been devoted to the numerical solution of reaction-diffusion (RD) partial differential equations (PDEs). The Crank–Nicolson (CN) method has been a common second-order time-stepping procedure. However, the CN method may introduce spurious oscillations for nonsmooth data unless the time step size is sufficiently small. This article studies a nonoscillatory second-order time-stepping procedure for RD equations, called a variable-θmethod, as a perturbation of the CN method. In each time level, the new method detects points of potential oscillations to implicitly resolve the solution there. The proposed time-stepping procedure is nonoscillatory and of a second-order temporal accuracy. Various examples are given to show effectiveness of the method. The article also performs a sensitivity analysis for the numerical solution of biological pattern forming models to conclude that the numerical solution is much more sensitive to the spatial mesh resolution than the temporal one. As the spatial resolution becomes higher for an improved accuracy, the CN method may produce spurious oscillations, while the proposed method results in stable solutions. |
| format | Article |
| id | doaj-art-4ad2f4795fed4775a09cf3f64f1e1a39 |
| institution | Kabale University |
| issn | 1076-2787 1099-0526 |
| language | English |
| publishDate | 2020-01-01 |
| publisher | Wiley |
| record_format | Article |
| series | Complexity |
| spelling | doaj-art-4ad2f4795fed4775a09cf3f64f1e1a392025-02-03T01:01:30ZengWileyComplexity1076-27871099-05262020-01-01202010.1155/2020/51637045163704A Nonoscillatory Second-Order Time-Stepping Procedure for Reaction-Diffusion EquationsPhilku Lee0George V. Popescu1Seongjai Kim2Department of Mathematics and Statistics, Mississippi State University, Starkville, Mississippi State, MS 39762, USAInstitute for Genomics, Biocomputing and Biotechnology, Mississippi State University, Starkville, Mississippi State, MS 39762, USADepartment of Mathematics and Statistics, Mississippi State University, Starkville, Mississippi State, MS 39762, USAAfter a theory of morphogenesis in chemical cells was introduced in the 1950s, much attention had been devoted to the numerical solution of reaction-diffusion (RD) partial differential equations (PDEs). The Crank–Nicolson (CN) method has been a common second-order time-stepping procedure. However, the CN method may introduce spurious oscillations for nonsmooth data unless the time step size is sufficiently small. This article studies a nonoscillatory second-order time-stepping procedure for RD equations, called a variable-θmethod, as a perturbation of the CN method. In each time level, the new method detects points of potential oscillations to implicitly resolve the solution there. The proposed time-stepping procedure is nonoscillatory and of a second-order temporal accuracy. Various examples are given to show effectiveness of the method. The article also performs a sensitivity analysis for the numerical solution of biological pattern forming models to conclude that the numerical solution is much more sensitive to the spatial mesh resolution than the temporal one. As the spatial resolution becomes higher for an improved accuracy, the CN method may produce spurious oscillations, while the proposed method results in stable solutions.http://dx.doi.org/10.1155/2020/5163704 |
| spellingShingle | Philku Lee George V. Popescu Seongjai Kim A Nonoscillatory Second-Order Time-Stepping Procedure for Reaction-Diffusion Equations Complexity |
| title | A Nonoscillatory Second-Order Time-Stepping Procedure for Reaction-Diffusion Equations |
| title_full | A Nonoscillatory Second-Order Time-Stepping Procedure for Reaction-Diffusion Equations |
| title_fullStr | A Nonoscillatory Second-Order Time-Stepping Procedure for Reaction-Diffusion Equations |
| title_full_unstemmed | A Nonoscillatory Second-Order Time-Stepping Procedure for Reaction-Diffusion Equations |
| title_short | A Nonoscillatory Second-Order Time-Stepping Procedure for Reaction-Diffusion Equations |
| title_sort | nonoscillatory second order time stepping procedure for reaction diffusion equations |
| url | http://dx.doi.org/10.1155/2020/5163704 |
| work_keys_str_mv | AT philkulee anonoscillatorysecondordertimesteppingprocedureforreactiondiffusionequations AT georgevpopescu anonoscillatorysecondordertimesteppingprocedureforreactiondiffusionequations AT seongjaikim anonoscillatorysecondordertimesteppingprocedureforreactiondiffusionequations AT philkulee nonoscillatorysecondordertimesteppingprocedureforreactiondiffusionequations AT georgevpopescu nonoscillatorysecondordertimesteppingprocedureforreactiondiffusionequations AT seongjaikim nonoscillatorysecondordertimesteppingprocedureforreactiondiffusionequations |