This paper discusses the stability of theoretical solutions for nonlinear multi-variable delay perturbation problems (MVDPP) of the form x′(t)=f(x(t),x(t-τ 1(t)),...,x(t-τ m(t)),y(t),y(t-τ 1(t)),...,y(t-τ m(t...This paper discusses the stability of theoretical solutions for nonlinear multi-variable delay perturbation problems (MVDPP) of the form x′(t)=f(x(t),x(t-τ 1(t)),...,x(t-τ m(t)),y(t),y(t-τ 1(t)),...,y(t-τ m(t))), and εy′(t)=g(x(t),x(t-τ 1(t)),...,x(t-τ m(t)),y(t),y(t-τ 1(t)),...,y(t-τ m(t))), where 0<ε1. A sufficient condition of stability for the systems is obtained. Additionally we prove the numerical solutions of the implicit Euler method are stable under this condition.展开更多
The numerical solution of large scale multi-dimensional convection diffusion equations often requires efficient parallel algorithms.In this work,we consider the extension of a recently proposed non-overlapping domain ...The numerical solution of large scale multi-dimensional convection diffusion equations often requires efficient parallel algorithms.In this work,we consider the extension of a recently proposed non-overlapping domain decomposition method for two dimensional time dependent convection diffusion equations with variable coefficients. By combining predictor-corrector technique,modified upwind differences with explicitimplicit coupling,the method under consideration provides intrinsic parallelism while maintaining good stability and accuracy.Moreover,for multi-dimensional problems, the method can be readily implemented on a multi-processor system and does not have the limitation on the choice of subdomains required by some other similar predictor-corrector or stabilized schemes.These properties of the method are demonstrated in this work through both rigorous mathematical analysis and numerical experiments.展开更多
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.展开更多
An algorithm composed of an iterative modified approximate factorization(MAF(k)) method with Navier-Stokes characteristic boundary conditions(NSCBC) is proposed for solving subsonic viscous flows.A transformation on t...An algorithm composed of an iterative modified approximate factorization(MAF(k)) method with Navier-Stokes characteristic boundary conditions(NSCBC) is proposed for solving subsonic viscous flows.A transformation on the matrix equation in MAF(k) is made in order to impose the implicit boundary conditions properly.To be in consistent with the implicit solver for the interior domain,an implicit scheme for NSCBC is formulated.The performance of the developed algorithm is investigated using spatially evolving zero pressure gradient boundary layer over a flat plate and a wall jet mixing with a cross flow over a flat plate with a square hole as the test cases.The numerical results are compared to the existing experimental datasets and a number of general correlations,together with other available numerical solutions,which demonstrate that the developed algorithm possesses promising capacity for simulating the subsonic viscous flows with large CFL number.展开更多
文摘This paper discusses the stability of theoretical solutions for nonlinear multi-variable delay perturbation problems (MVDPP) of the form x′(t)=f(x(t),x(t-τ 1(t)),...,x(t-τ m(t)),y(t),y(t-τ 1(t)),...,y(t-τ m(t))), and εy′(t)=g(x(t),x(t-τ 1(t)),...,x(t-τ m(t)),y(t),y(t-τ 1(t)),...,y(t-τ m(t))), where 0<ε1. A sufficient condition of stability for the systems is obtained. Additionally we prove the numerical solutions of the implicit Euler method are stable under this condition.
基金the National Natural Science Foundation of China(No.10571017)supported in part by the National Natural Science Foundation of China(No.60533020)supported in part by NSF DMS 0712744
文摘The numerical solution of large scale multi-dimensional convection diffusion equations often requires efficient parallel algorithms.In this work,we consider the extension of a recently proposed non-overlapping domain decomposition method for two dimensional time dependent convection diffusion equations with variable coefficients. By combining predictor-corrector technique,modified upwind differences with explicitimplicit coupling,the method under consideration provides intrinsic parallelism while maintaining good stability and accuracy.Moreover,for multi-dimensional problems, the method can be readily implemented on a multi-processor system and does not have the limitation on the choice of subdomains required by some other similar predictor-corrector or stabilized schemes.These properties of the method are demonstrated in this work through both rigorous mathematical analysis and numerical experiments.
文摘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.
文摘An algorithm composed of an iterative modified approximate factorization(MAF(k)) method with Navier-Stokes characteristic boundary conditions(NSCBC) is proposed for solving subsonic viscous flows.A transformation on the matrix equation in MAF(k) is made in order to impose the implicit boundary conditions properly.To be in consistent with the implicit solver for the interior domain,an implicit scheme for NSCBC is formulated.The performance of the developed algorithm is investigated using spatially evolving zero pressure gradient boundary layer over a flat plate and a wall jet mixing with a cross flow over a flat plate with a square hole as the test cases.The numerical results are compared to the existing experimental datasets and a number of general correlations,together with other available numerical solutions,which demonstrate that the developed algorithm possesses promising capacity for simulating the subsonic viscous flows with large CFL number.