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.展开更多
Several numerical methods of differential equations and their applications in ballistic calculation are discussed for the purpose of simplification of the dynamic differential equations of projectile trajectory.Progra...Several numerical methods of differential equations and their applications in ballistic calculation are discussed for the purpose of simplification of the dynamic differential equations of projectile trajectory.Program simulations of Euler method,Heun method,lassic fourth-order Runge Kutta（RK4）method,ABM method and Hamming method are achieved based on Matlab.In addtion,the approximate solutions,local truncation errors and calculation time of the dynamic differential equations are obtained.By analyzing the simultaion results,the advantages and disadvantages of these methods are compared,which provides a basis for choice of ballistic calculation methods.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Malaria,a devastating disease caused by the Plasmodium parasite and transmitted through the bites of female Anopheles mosquitoes,remains a significant public health concern,claiming over 600,000 lives annually,predomi...Malaria,a devastating disease caused by the Plasmodium parasite and transmitted through the bites of female Anopheles mosquitoes,remains a significant public health concern,claiming over 600,000 lives annually,predominantly among children.Novel tools,including the application of Wolbachia,are being developed to combat malaria-transmitting mosquitoes.This study presents a modified susceptible-exposed-infectious-recovered-susceptible(SEIRS)compartmental mathematical model to evaluate the impact of awareness-based control measures on malaria transmission dynamics,incorporating mosquito interactions and seasonality.Employing the next-generation matrix approach,we calculated a basic reproduction number(R0)of 2.4537,indicating that without robust control measures,the disease will persist in the human population.The model equations were solved numerically using fourth and fifth-order Runge-Kutta methods.The model was fitted to malaria incidence data from Kenya spanning 2000 to 2021 using least squares curve fitting.The fitting algorithm yielded a mean absolute error(MAE)of 2.6463 when comparing the actual data points to the simulated values of infectious human population(Ih).This finding indicates that the proposed mathematical model closely aligns with the recorded malaria incidence data.The optimal values of the model parameters were estimated from the fitting algorithm,and future malaria dynamics were projected for the next decade.The research findings suggest that social media-based awareness campaigns,coupled with specific optimization control measures and effective management methods,offer the most cost-effective approach to managing malaria.展开更多
Many theoretical studies have been developed to study the spectral response of a fiber Bragg grating (FBG) under non-uniform strain distribution along the length of FBG in recent years. However, almost no experiments ...Many theoretical studies have been developed to study the spectral response of a fiber Bragg grating (FBG) under non-uniform strain distribution along the length of FBG in recent years. However, almost no experiments were designed to obtain the evolution of the spectrum when a FBG is subjected to non-uniform strain. In this paper, the spectral responses of a FBG under non-uniform strain distributions are given and a numerical simulation based on the Runge-Kutta method is introduced to investigate the responses of the FBG under some typical non-uniform transverse strain fields, including both linear strain gradient and quadratic strain field. Experiment is carried out by using loads applied at different locations near the FBG. Good agreements between experimental results and numerical simulations are obtained.展开更多
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.展开更多
A wavelet method is proposed to solve the Burgers’equation.Following this method,this nonlinear partial differential equation is first transformed into a system of ordinary differential equations using the modified w...A wavelet method is proposed to solve the Burgers’equation.Following this method,this nonlinear partial differential equation is first transformed into a system of ordinary differential equations using the modified wavelet Galerkin method recently developed by the authors.Then,the classical fourth-order explicit Runge–Kutta method is employed to solve the resulting system of ordinary differential equations.Such a wavelet-based solution procedure has been justified by solving two test examples:results demonstrate that the proposed method has a much better accuracy and efficiency than many other existing numerical methods,and whose order of convergence can go up to 5.Most importantly,our results also indicate that the present wavelet method can readily deal with those fluid dynamics problems with high Reynolds numbers.展开更多
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.展开更多
基金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.
文摘Several numerical methods of differential equations and their applications in ballistic calculation are discussed for the purpose of simplification of the dynamic differential equations of projectile trajectory.Program simulations of Euler method,Heun method,lassic fourth-order Runge Kutta（RK4）method,ABM method and Hamming method are achieved based on Matlab.In addtion,the approximate solutions,local truncation errors and calculation time of the dynamic differential equations are obtained.By analyzing the simultaion results,the advantages and disadvantages of these methods are compared,which provides a basis for choice of ballistic calculation methods.
文摘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.
文摘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.
文摘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.
基金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.
文摘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.
文摘Malaria,a devastating disease caused by the Plasmodium parasite and transmitted through the bites of female Anopheles mosquitoes,remains a significant public health concern,claiming over 600,000 lives annually,predominantly among children.Novel tools,including the application of Wolbachia,are being developed to combat malaria-transmitting mosquitoes.This study presents a modified susceptible-exposed-infectious-recovered-susceptible(SEIRS)compartmental mathematical model to evaluate the impact of awareness-based control measures on malaria transmission dynamics,incorporating mosquito interactions and seasonality.Employing the next-generation matrix approach,we calculated a basic reproduction number(R0)of 2.4537,indicating that without robust control measures,the disease will persist in the human population.The model equations were solved numerically using fourth and fifth-order Runge-Kutta methods.The model was fitted to malaria incidence data from Kenya spanning 2000 to 2021 using least squares curve fitting.The fitting algorithm yielded a mean absolute error(MAE)of 2.6463 when comparing the actual data points to the simulated values of infectious human population(Ih).This finding indicates that the proposed mathematical model closely aligns with the recorded malaria incidence data.The optimal values of the model parameters were estimated from the fitting algorithm,and future malaria dynamics were projected for the next decade.The research findings suggest that social media-based awareness campaigns,coupled with specific optimization control measures and effective management methods,offer the most cost-effective approach to managing malaria.
基金supported by the National High Technology Research and Development Program of China (No.2007AA03Z117)the Key Program of National Natural Science Foundation of China (No.50830201)
文摘Many theoretical studies have been developed to study the spectral response of a fiber Bragg grating (FBG) under non-uniform strain distribution along the length of FBG in recent years. However, almost no experiments were designed to obtain the evolution of the spectrum when a FBG is subjected to non-uniform strain. In this paper, the spectral responses of a FBG under non-uniform strain distributions are given and a numerical simulation based on the Runge-Kutta method is introduced to investigate the responses of the FBG under some typical non-uniform transverse strain fields, including both linear strain gradient and quadratic strain field. Experiment is carried out by using loads applied at different locations near the FBG. Good agreements between experimental results and numerical simulations are obtained.
基金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.
基金supported by the National Natural Science Foundation of China(Grant Nos.11032006,11072094,and 11121202)the Ph.D.Program Foundation of Ministry of Education of China(Grant No.20100211110022)+2 种基金the National Key Project of Magneto-Constrained Fusion Energy Development Program(Grant No.2013GB110002)the Fundamental Research Funds for the Central Universities(Grant No.lzujbky-2013-1)the Scholarship Award for Excellent Doctoral Student granted by the Lanzhou University
文摘A wavelet method is proposed to solve the Burgers’equation.Following this method,this nonlinear partial differential equation is first transformed into a system of ordinary differential equations using the modified wavelet Galerkin method recently developed by the authors.Then,the classical fourth-order explicit Runge–Kutta method is employed to solve the resulting system of ordinary differential equations.Such a wavelet-based solution procedure has been justified by solving two test examples:results demonstrate that the proposed method has a much better accuracy and efficiency than many other existing numerical methods,and whose order of convergence can go up to 5.Most importantly,our results also indicate that the present wavelet method can readily deal with those fluid dynamics problems with high Reynolds numbers.
文摘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.
