An autocatalytic biochemical system in the presence of recycling enzyme is solved numerically using two numerical methods based on finite difference schemes. The first method is the well known Euler method which is an...An autocatalytic biochemical system in the presence of recycling enzyme is solved numerically using two numerical methods based on finite difference schemes. The first method is the well known Euler method which is an explicit method, whereas the second method is implicit. Although the implicit method, method 2, is first-order accurate in time it converges to the fixed point(s) for large time step, L Numerical results show the existence of hard excitation and birhythmicity.展开更多
In this work, the MMC-TDGL equation, a stochastic Cahn-Hilliard equation, is solved numerically by using the finite difference method in combination with a convex splitting technique of the energy functional.For the n...In this work, the MMC-TDGL equation, a stochastic Cahn-Hilliard equation, is solved numerically by using the finite difference method in combination with a convex splitting technique of the energy functional.For the non-stochastic case, we develop an unconditionally energy stable difference scheme which is proved to be uniquely solvable. For the stochastic case, by adopting the same splitting of the energy functional, we construct a similar and uniquely solvable difference scheme with the discretized stochastic term. The resulted schemes are nonlinear and solved by Newton iteration. For the long time simulation, an adaptive time stepping strategy is developed based on both first- and second-order derivatives of the energy. Numerical experiments are carried out to verify the energy stability, the efficiency of the adaptive time stepping and the effect of the stochastic term.展开更多
文摘An autocatalytic biochemical system in the presence of recycling enzyme is solved numerically using two numerical methods based on finite difference schemes. The first method is the well known Euler method which is an explicit method, whereas the second method is implicit. Although the implicit method, method 2, is first-order accurate in time it converges to the fixed point(s) for large time step, L Numerical results show the existence of hard excitation and birhythmicity.
基金supported by the Hong Kong General Research Fund (Grant Nos. 202112, 15302214 and 509213)National Natural Science Foundation of China/Research Grants Council Joint Research Scheme (Grant Nos. N HKBU204/12 and 11261160486)+1 种基金National Natural Science Foundation of China (Grant No. 11471046)the Ministry of Education Program for New Century Excellent Talents Project (Grant No. NCET-12-0053)
文摘In this work, the MMC-TDGL equation, a stochastic Cahn-Hilliard equation, is solved numerically by using the finite difference method in combination with a convex splitting technique of the energy functional.For the non-stochastic case, we develop an unconditionally energy stable difference scheme which is proved to be uniquely solvable. For the stochastic case, by adopting the same splitting of the energy functional, we construct a similar and uniquely solvable difference scheme with the discretized stochastic term. The resulted schemes are nonlinear and solved by Newton iteration. For the long time simulation, an adaptive time stepping strategy is developed based on both first- and second-order derivatives of the energy. Numerical experiments are carried out to verify the energy stability, the efficiency of the adaptive time stepping and the effect of the stochastic term.