We employ the Duan-Rach-Wazwaz modified Adomian decomposition method for solving initial value problems for the systems of nonlinear ordinary differential equations numerically. In order to confirm practicality, robus...We employ the Duan-Rach-Wazwaz modified Adomian decomposition method for solving initial value problems for the systems of nonlinear ordinary differential equations numerically. In order to confirm practicality, robustness and reliability of the method, we compare the results from the modified Adomian decomposition method with those from the MATHEMATICA solutions and also from the fourth-order Runge Kutta method solutions in some cases. Furthermore, we apply Padé approximants technique to improve the solutions of the modified decomposition method whenever the exact solutions exist.展开更多
In this paper, the multi-scale Runge-Kutta_Galerkin method is developed for solving the evolution equations, with the spatial variables of the equations being discretized by the multi-scale Galerkin method based on th...In this paper, the multi-scale Runge-Kutta_Galerkin method is developed for solving the evolution equations, with the spatial variables of the equations being discretized by the multi-scale Galerkin method based on the multi-scale orthogonal bases in Ho (a, b) and then the classical fourth order explicit Runge-Kutta method being applied to solve the resulting initial problem of the ordinary differential equations for the coefficients of the approximate solution. The proposed numerical scheme is validated by applications to the Burgers equation (nonlinear convection-diffusion problem), the KdV equation (single solitary and 2-solitary wave problems) and the KdV-Burgers equation, where analytical solutions are available for estimating the errors. Numerical results show that using the algorithm we can solve these equations stably without the need for extra stabilization processes and obtain accurate solutions that agree very well with the corresponding exact solutions in all cases.展开更多
With the advantages of noncontact,high accuracy,and high flexibility,optical tweezers hold huge potential for micro-manipulation and force measurement.However,the majority of previous research focused on the state of ...With the advantages of noncontact,high accuracy,and high flexibility,optical tweezers hold huge potential for micro-manipulation and force measurement.However,the majority of previous research focused on the state of the motion of particles in the optical trap,but paid little attention to the early dynamic process between the initial state of the particles and the optical trap.Note that the viscous forces can greatly affect the motion of micro-spheres.In this paper,based on the equations of Newtonian mechanics,we investigate the dynamics of laser-trapped micro-spheres in the surrounding environment with different viscosity coefficients.Through the calculations,over time the particle trajectory clearly reveals the subtle details of the optical capture process,including acceleration,deceleration,turning,and reciprocating oscillation.The time to equilibrium mainly depends on the corresponding damping coefficient of the surrounding environment and the oscillation frequency of the optical tweezers.These studies are essential for understanding various mechanisms to engineer the mechanical motion behavior of molecules or microparticles in liquid or air.展开更多
Nonlinear methods are combined with Runge Kutta methods to develop A stable explicit nonlinear Runge Kutta methods for solving stiff differential equations and a class of the third order formulae are constructed.I...Nonlinear methods are combined with Runge Kutta methods to develop A stable explicit nonlinear Runge Kutta methods for solving stiff differential equations and a class of the third order formulae are constructed.It avoids solving the nonlinear equations which implicit methods must solve. Implementation is very simple and the computation cost for each step is small.This paper uses a shift transformation to avoid the order reduction of nonlinear methods at y = 0 . Thus the method is very practicable. Numerical tests show that the method is more efficient than explicit methods or implicit methods of the same order.展开更多
A finite-volume high-resolution numerical model for coupling the shallow water flows and pollutant diffusions was presented based on using a hybrid TVD scheme in space discretization and a Runge-Kutta method in time d...A finite-volume high-resolution numerical model for coupling the shallow water flows and pollutant diffusions was presented based on using a hybrid TVD scheme in space discretization and a Runge-Kutta method in time discretization. Numerical simulations for modelling dam-break, enlarging open channel flow and pollutant dispersion were implemented and compared with experimental data or other published computations. The validation of this method shows that it can not only deal with the problem involving discontinuities and unsteady flows, but also solve the general shallow water flows and pollutant diffusions.展开更多
Error analysis methods in frequency domain are developed in this paper for determining the characteristic root and transfer function errors when the linear multipass algorithms are used to solve linear differential eq...Error analysis methods in frequency domain are developed in this paper for determining the characteristic root and transfer function errors when the linear multipass algorithms are used to solve linear differential equations. The relation between the local truncation error in time domain and the error in frequency domain is established, which is the basis for developing the error estimation methods. The error estimation methods for the digital simulation model constructed by using the Runge-Kutta algorithms and the linear multistep predictor-corrector algorithms are also given.展开更多
This work is concerned with Magnetohydrodynamic viscous flow due to a shrinking sheet in the presence of suction. The cases of two dimensional and axisymmetric shrinking are discussed. The governing boundary layer equ...This work is concerned with Magnetohydrodynamic viscous flow due to a shrinking sheet in the presence of suction. The cases of two dimensional and axisymmetric shrinking are discussed. The governing boundary layer equations are written into a dimensionless form by similarity transformations. The transformed coupled nonlinear ordinary differential equations are numerically solved by using an advanced numeric technique. Favorability comparisons with previously published work are presented. Numerical results for the dimensionless velocity, temperature and concentration profiles as well as for the skin friction, heat and mass transfer and deposition rate are obtained and displayed graphically for pertinent parameters to show interesting aspects of the solution.展开更多
An Euler solver is developed for calculating the flowfield around the propeller. The Euler equations are recast in the blade attached rotating coordinate system with absolute flow variables.The finite volume explicit...An Euler solver is developed for calculating the flowfield around the propeller. The Euler equations are recast in the blade attached rotating coordinate system with absolute flow variables.The finite volume explicit Runge Kutta time stepping scheme is employed for solving the equations.The presented method is applied for a two bladed propeller NACA10 (3) (066) 033.An algebraic O O grid is generated. Comparison of the mumerical results shows good agreement with the experimental data.展开更多
The numerical approximations of the dynamical systems governed by semilinear parabolic equations are considered. An abstract framework for long time error estimates is established. When applied to reaction diffusion...The numerical approximations of the dynamical systems governed by semilinear parabolic equations are considered. An abstract framework for long time error estimates is established. When applied to reaction diffusion equation, Navier Stokes equations and Chan Hilliard equation, approximated by Galerkin and nonlinear Galerkin methods in space and by Runge Kutta method in time, our framework yields error estimates uniform in time.展开更多
A new transfer matrix method for long-period fiber gratings with coupled multiple cladding modes is proposed and numerically characterized. The transmission spectra of uniform and non-uniform longperiod fiber gratings...A new transfer matrix method for long-period fiber gratings with coupled multiple cladding modes is proposed and numerically characterized. The transmission spectra of uniform and non-uniform longperiod fiber gratings are numerically characterized. The theoretical results excellently agree with the experimental measurements, Compared with commonly used methods, such as using the fourth-order adaptive step size control of the Runge-Kutta algorithm in solving the coupled mode equation, the new transfer matrix method exhibits a faster calculation speed.展开更多
The aim of this work is to realize a new numerical program based on the development of a mathematical model allowing determining the parameters of the supersonic flow through a conical shock under hypothesis at high t...The aim of this work is to realize a new numerical program based on the development of a mathematical model allowing determining the parameters of the supersonic flow through a conical shock under hypothesis at high temperature, in the context of correcting the perfect gas model. In this case, the specific heat at constant pressure does not remain constant and varies with the increase of temperature. The stagnation temperature becomes an important parameter in the calculation.The mathematical model is presented by the numerical resolution of a system of first-order nonlinear differential equations with three coupled unknowns for initial conditions. The numerical resolution is made by adapting the higher order Runge Kutta method. The parameters through the conical shock can be determined by considering a new model of an oblique shock at high temperature. All isentropic parameters of after the shock flow depend on the deviation of the flow from the transverse direction. The comparison of the results is done with the perfect gas model for low stagnation temperatures, upstream Mach number and cone deviation angle. A calculation of the error is made between our high temperature model and the perfect gas model. The application is made for air.展开更多
Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin(DG) discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta(IRK) scheme...Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin(DG) discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta(IRK) scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.展开更多
文摘We employ the Duan-Rach-Wazwaz modified Adomian decomposition method for solving initial value problems for the systems of nonlinear ordinary differential equations numerically. In order to confirm practicality, robustness and reliability of the method, we compare the results from the modified Adomian decomposition method with those from the MATHEMATICA solutions and also from the fourth-order Runge Kutta method solutions in some cases. Furthermore, we apply Padé approximants technique to improve the solutions of the modified decomposition method whenever the exact solutions exist.
基金supported by the National Marine Public Welfare Research Projects of China(Grant No.201005002)the National Natural Science Foundation of China(Grant No.11071264)the Fundamental Research Funds for the Cen-tral Universities
文摘In this paper, the multi-scale Runge-Kutta_Galerkin method is developed for solving the evolution equations, with the spatial variables of the equations being discretized by the multi-scale Galerkin method based on the multi-scale orthogonal bases in Ho (a, b) and then the classical fourth order explicit Runge-Kutta method being applied to solve the resulting initial problem of the ordinary differential equations for the coefficients of the approximate solution. The proposed numerical scheme is validated by applications to the Burgers equation (nonlinear convection-diffusion problem), the KdV equation (single solitary and 2-solitary wave problems) and the KdV-Burgers equation, where analytical solutions are available for estimating the errors. Numerical results show that using the algorithm we can solve these equations stably without the need for extra stabilization processes and obtain accurate solutions that agree very well with the corresponding exact solutions in all cases.
基金Project supported by the National Natural Science Foundation of China(Grant No.11804399)the Special Funds for Basic Scientific Research at the Central University of South-Central University for Nationalities(Grant No.CZQ20018)Special Funds for Basic Scientific Research at Central Universities(Grant No.YZZ17005)。
文摘With the advantages of noncontact,high accuracy,and high flexibility,optical tweezers hold huge potential for micro-manipulation and force measurement.However,the majority of previous research focused on the state of the motion of particles in the optical trap,but paid little attention to the early dynamic process between the initial state of the particles and the optical trap.Note that the viscous forces can greatly affect the motion of micro-spheres.In this paper,based on the equations of Newtonian mechanics,we investigate the dynamics of laser-trapped micro-spheres in the surrounding environment with different viscosity coefficients.Through the calculations,over time the particle trajectory clearly reveals the subtle details of the optical capture process,including acceleration,deceleration,turning,and reciprocating oscillation.The time to equilibrium mainly depends on the corresponding damping coefficient of the surrounding environment and the oscillation frequency of the optical tweezers.These studies are essential for understanding various mechanisms to engineer the mechanical motion behavior of molecules or microparticles in liquid or air.
基金Supported by the National Natural Science Foundationof China(No.195 710 4 6 )
文摘Nonlinear methods are combined with Runge Kutta methods to develop A stable explicit nonlinear Runge Kutta methods for solving stiff differential equations and a class of the third order formulae are constructed.It avoids solving the nonlinear equations which implicit methods must solve. Implementation is very simple and the computation cost for each step is small.This paper uses a shift transformation to avoid the order reduction of nonlinear methods at y = 0 . Thus the method is very practicable. Numerical tests show that the method is more efficient than explicit methods or implicit methods of the same order.
文摘A finite-volume high-resolution numerical model for coupling the shallow water flows and pollutant diffusions was presented based on using a hybrid TVD scheme in space discretization and a Runge-Kutta method in time discretization. Numerical simulations for modelling dam-break, enlarging open channel flow and pollutant dispersion were implemented and compared with experimental data or other published computations. The validation of this method shows that it can not only deal with the problem involving discontinuities and unsteady flows, but also solve the general shallow water flows and pollutant diffusions.
基金This project was supported by the National Natural Science Foundation of China (No. 19871080).
文摘Error analysis methods in frequency domain are developed in this paper for determining the characteristic root and transfer function errors when the linear multipass algorithms are used to solve linear differential equations. The relation between the local truncation error in time domain and the error in frequency domain is established, which is the basis for developing the error estimation methods. The error estimation methods for the digital simulation model constructed by using the Runge-Kutta algorithms and the linear multistep predictor-corrector algorithms are also given.
文摘This work is concerned with Magnetohydrodynamic viscous flow due to a shrinking sheet in the presence of suction. The cases of two dimensional and axisymmetric shrinking are discussed. The governing boundary layer equations are written into a dimensionless form by similarity transformations. The transformed coupled nonlinear ordinary differential equations are numerically solved by using an advanced numeric technique. Favorability comparisons with previously published work are presented. Numerical results for the dimensionless velocity, temperature and concentration profiles as well as for the skin friction, heat and mass transfer and deposition rate are obtained and displayed graphically for pertinent parameters to show interesting aspects of the solution.
文摘An Euler solver is developed for calculating the flowfield around the propeller. The Euler equations are recast in the blade attached rotating coordinate system with absolute flow variables.The finite volume explicit Runge Kutta time stepping scheme is employed for solving the equations.The presented method is applied for a two bladed propeller NACA10 (3) (066) 033.An algebraic O O grid is generated. Comparison of the mumerical results shows good agreement with the experimental data.
文摘The numerical approximations of the dynamical systems governed by semilinear parabolic equations are considered. An abstract framework for long time error estimates is established. When applied to reaction diffusion equation, Navier Stokes equations and Chan Hilliard equation, approximated by Galerkin and nonlinear Galerkin methods in space and by Runge Kutta method in time, our framework yields error estimates uniform in time.
基金supported by the National Natural Science Foundation of China (No. 61040016)the Open Foundation for Henan Provincial Open Laboratory for Control Engineering Key Disciplines China (No. KG2009-16)the Doctor Foundation for Henan Polytechnic University China (No. 648393)
文摘A new transfer matrix method for long-period fiber gratings with coupled multiple cladding modes is proposed and numerically characterized. The transmission spectra of uniform and non-uniform longperiod fiber gratings are numerically characterized. The theoretical results excellently agree with the experimental measurements, Compared with commonly used methods, such as using the fourth-order adaptive step size control of the Runge-Kutta algorithm in solving the coupled mode equation, the new transfer matrix method exhibits a faster calculation speed.
文摘The aim of this work is to realize a new numerical program based on the development of a mathematical model allowing determining the parameters of the supersonic flow through a conical shock under hypothesis at high temperature, in the context of correcting the perfect gas model. In this case, the specific heat at constant pressure does not remain constant and varies with the increase of temperature. The stagnation temperature becomes an important parameter in the calculation.The mathematical model is presented by the numerical resolution of a system of first-order nonlinear differential equations with three coupled unknowns for initial conditions. The numerical resolution is made by adapting the higher order Runge Kutta method. The parameters through the conical shock can be determined by considering a new model of an oblique shock at high temperature. All isentropic parameters of after the shock flow depend on the deviation of the flow from the transverse direction. The comparison of the results is done with the perfect gas model for low stagnation temperatures, upstream Mach number and cone deviation angle. A calculation of the error is made between our high temperature model and the perfect gas model. The application is made for air.
文摘Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin(DG) discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta(IRK) scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.