This study presents numerical algorithms for solving a class of equations that partly consists of derivatives of the unknown state at previous certain times, as well as an integro-differential term containing a weakly...This study presents numerical algorithms for solving a class of equations that partly consists of derivatives of the unknown state at previous certain times, as well as an integro-differential term containing a weakly singular kernel. These equations are types of integro-differential equation of the second kind and were originally obtained from an aeroelasticity problem. One of the main contributions of this study is to propose numerical algorithms that do not involve transforming the original equation into the corresponding Volterra equation, but still enable the numerical solution of the original equation to be determined. The feasibility of the proposed numerical algorithm is demonstrated by applying examples in measuring the maximum errors with exact solutions at every computed nodes and calculating the corresponding numerical rates of convergence thereafter.展开更多
In this paper, a modified algorithm is proposed for solving linear integro-differential equations of the second kind. The main idea is based on applying Romberg extrapolation algorithm (REA), on Trapezoidal rule. In a...In this paper, a modified algorithm is proposed for solving linear integro-differential equations of the second kind. The main idea is based on applying Romberg extrapolation algorithm (REA), on Trapezoidal rule. In accordance with the computational perspective, the comparison has shown that Adomian decomposition approach is more effective to be utilized. The numerical results show that the modified algorithm has been successfully applied to the linear integro-differential equations and the comparisons with some existing methods appeared in the literature reveal that the modified algorithm is more accurate and convenient.展开更多
In this paper, the authors show that the general linear second order ordinary Differential Equation can be formulated as an optimization problem and that evolutionary algorithms for solving optimization problems can a...In this paper, the authors show that the general linear second order ordinary Differential Equation can be formulated as an optimization problem and that evolutionary algorithms for solving optimization problems can also be adapted for solving the formulated problem. The authors propose a polynomial based scheme for achieving the above objectives. The coefficients of the proposed scheme are approximated by an evolutionary algorithm known as Differential Evolution (DE). Numerical examples with good results show the accuracy of the proposed method compared with some existing methods.展开更多
The idea of AC = BD was applied to solve the nonlinear differential equations. Suppose that Au = 0 is a given equation to he solved and Dv = 0 is an equation to be easily solved. If the transformation u = Cv is obtain...The idea of AC = BD was applied to solve the nonlinear differential equations. Suppose that Au = 0 is a given equation to he solved and Dv = 0 is an equation to be easily solved. If the transformation u = Cv is obtained so that v satisfies Dv = 0, then the solutions for Au = 0 can be found. In order to illustrate this approach, several examples about the transformation C are given.展开更多
An accurate and rapid method for solving radiative transfer equation is presented in this paper. According to the fact that the multiple scattering component of radiance is less sensitive to the error of phase functio...An accurate and rapid method for solving radiative transfer equation is presented in this paper. According to the fact that the multiple scattering component of radiance is less sensitive to the error of phase function than the single scattering component is,we calculate the multiple scattering component by using delta-Eddington approximation and the single scattering component by solving radiative transfer equation. On the ground, when multiple sattering component is small, for example, when the total optical depth T is small, the accurate radiance can be obtained with this method. For the need of the space remote sensing, the upward radiance at the top of the atmosphere is mainly studied, and an approximate expression is presented to correct the multiple scattering component. Compared with the more precise Gauss-Seidel method.the results from this method show an accuracy of better than 10% when zenith angle 0 < 50 掳 and T < 1. The computational speed of this method is, however, much faster than that of Gauss-Seidel method.展开更多
When soldering electronic components onto circuit boards,the temperature curves of the reflow ovens across different zones and the conveyor belt speed significantly influence the product quality.This study focuses on ...When soldering electronic components onto circuit boards,the temperature curves of the reflow ovens across different zones and the conveyor belt speed significantly influence the product quality.This study focuses on optimizing the furnace temperature curve under varying settings of reflow oven zone temperatures and conveyor belt speeds.To address this,the research sequentially develops a heat transfer model for reflow soldering,an optimization model for reflow furnace conditions using the differential evolution algorithm,and an evaluation and decision model combining the differential evolution algorithm with the Technique for Order Preference by Similarity to Ideal Solution(TOPSIS)method.This approach aims to determine the optimal furnace temperature curve,zone temperatures of the reflow oven,and the conveyor belt speed.展开更多
A variation of the direct Taylor expansion algorithm is suggested and applied to several linear and nonlinear differential equations of interest in physics and engineering, and the results are compared with those obta...A variation of the direct Taylor expansion algorithm is suggested and applied to several linear and nonlinear differential equations of interest in physics and engineering, and the results are compared with those obtained from other algorithms. It is shown that the suggested algorithm competes strongly with other existing algorithms, both in accuracy and ease of application, while demanding a shorter computation time.展开更多
This paper proposed the scheme of transmission lines distance protection based on differential equation algorithms (DEA) and Hilbert-Huang transform (HHT). The measured impedance based on EDA is affected by various fa...This paper proposed the scheme of transmission lines distance protection based on differential equation algorithms (DEA) and Hilbert-Huang transform (HHT). The measured impedance based on EDA is affected by various factors, such as the distributed capacitance, the transient response characteristics of current transformer and voltage transformer, etc. In order to overcome this problem, the proposed scheme applies HHT to improve the apparent impedance estimated by DEA. Empirical mode decomposition (EMD) is used to decompose the data set from DEA into the intrinsic mode functions (IMF) and the residue. This residue has monotonic trend and is used to evaluate the impedance of faulty line. Simulation results show that the proposed scheme improves significantly the accuracy of the estimated impedance.展开更多
The multisplitting algorithm for solving large systems of ordinary differential equations on parallel computers was introduced by Jeltsch and Pohl in [1]. On fixed time intervals conver gence results could be derived ...The multisplitting algorithm for solving large systems of ordinary differential equations on parallel computers was introduced by Jeltsch and Pohl in [1]. On fixed time intervals conver gence results could be derived if the subsystems are solving exactly.Firstly,in theis paper,we deal with an extension of the waveform relaxation algorithm by us ing multisplittin AOR method based on an overlapping block decomposition. We restricted our selves to equidistant timepoints and dealed with the case that an implicit integration method was used to solve the subsystems numerically in parallel. Then we have proved convergence of multi splitting AOR waveform relaxation algorithm on a fixed window containing a finite number of timepoints.展开更多
Lie algorithm combined with differential form Wu's method is used to complete the symmetry classification of partial differential equations(PDEs)containing arbitrary parameter.This process can be reduced to solve ...Lie algorithm combined with differential form Wu's method is used to complete the symmetry classification of partial differential equations(PDEs)containing arbitrary parameter.This process can be reduced to solve a large system of determining equations,which seems rather difficult to solve,then the differential form Wu's method is used to decompose the determining equations into a series of equations,which are easy to solve.To illustrate the usefulness of this method,we apply it to some test problems,and the results show the performance of the present work.展开更多
In order to solve the problem of vibration bounce caused by the contact between moving and stationary contacts in the process of switching on,two-degree-of-freedom motion differential equation of the contact system is...In order to solve the problem of vibration bounce caused by the contact between moving and stationary contacts in the process of switching on,two-degree-of-freedom motion differential equation of the contact system is established.Genetic algorithm is used to optimize the pull in process of AC contactor.The whole process of contact bounce was observed and analyzed by high-speed photography experiment.The theory and experimental results were very similar.The iron core has collided before the contact is separated,which further aggravates the contact bounce.When the iron core bounces collided again,the bounce of the contact was not affected.During the operation of the contactor,the movement of the moving iron core will cause slight vibration of the system.The contact bounce time and the maximum amplitude are reduced.The research results provide a theoretical basis for further control and reduction of contact bounce.展开更多
Many engineering problems can be reduced to the solution of a variable coefficient differential equation. In this paper, the exact analytic method is suggested to solve variable coefficient differential equations unde...Many engineering problems can be reduced to the solution of a variable coefficient differential equation. In this paper, the exact analytic method is suggested to solve variable coefficient differential equations under arbitrary boundary condition. By this method, the general computation formal is obtained. Its convergence in proved. We can get analytic expressions which converge to exact solution and its higher order derivatives uniformy Four numerical examples are given, which indicate that satisfactory results can he obtanedby this method.展开更多
New form of necessary conditions for optimality (NCO) is considered. They can be useful for design the direct infinite- dimensional optimization algorithms for systems described by partial differential equations (PDE)...New form of necessary conditions for optimality (NCO) is considered. They can be useful for design the direct infinite- dimensional optimization algorithms for systems described by partial differential equations (PDE). Appropriate algo-rithms for unconstrained minimizing a functional are considered and tested. To construct the algorithms, new form of NCO is used. Such approach demonstrates fast uniform convergence at optimal solution in infinite-dimensional space.展开更多
We extend LeVeque's wave propagation algorithm,a widely used finite volume method for hyperbolic partial differential equations,to a third-order accurate method.The resulting scheme shares main properties with the...We extend LeVeque's wave propagation algorithm,a widely used finite volume method for hyperbolic partial differential equations,to a third-order accurate method.The resulting scheme shares main properties with the original method,i.e.,it is based on a wave decomposition at grid cell interfaces,it can be used to approximate hyperbolic problems in divergence form as well as in quasilinear form and limiting is introduced in the form of a wave limiter.展开更多
Acceptable glycemic control when examining the effects of meals was </span></span><span><span><span style="font-family:"">achieved when combining basal insulin therapy and ...Acceptable glycemic control when examining the effects of meals was </span></span><span><span><span style="font-family:"">achieved when combining basal insulin therapy and high concentration insulin injection before a meal, when using a PID controller (Proportionality, Integrity and Derivative actions) alone, when using a PID controller with basal insulin therapy and when combining the three methods of insulin delivery. Naturally, a type 1 diabetic must inject himself with insulin in well-measured doses. Thus, the management and control of diabetes become a complex task when one must be considered the disturbance due to nutrition and sports activity. This concern has been at the center of much research through different approaches through mathematical methods and Artificial Intelligence methods. This article simulates a physiological model of glycemic control in type 1 diabetics by a PID regulatory mechanism, in the context of disturbances caused by the patient’s meals and athletic activity.展开更多
A heuristic technique is developed for a nonlinear magnetohydrodynamics （MHD） Jeffery-Hamel problem with the help of the feed-forward artificial neural net- work （ANN） optimized with the genetic algorithm （GA） a...A heuristic technique is developed for a nonlinear magnetohydrodynamics （MHD） Jeffery-Hamel problem with the help of the feed-forward artificial neural net- work （ANN） optimized with the genetic algorithm （GA） and the sequential quadratic programming （SQP） method. The twodimensional （2D） MHD Jeffery-Hamel problem is transformed into a higher order boundary value problem （BVP） of ordinary differential equations （ODEs）. The mathematical model of the transformed BVP is formulated with the ANN in an unsupervised manner. The training of the weights of the ANN is carried out with the evolutionary calculation based on the GA hybridized with the SQP method for the rapid local convergence. The proposed scheme is evaluated on the variants of the Jeffery-Hamel flow by varying the Reynold number, the Hartmann number, and the an- gles of the walls. A large number of simulations are performed with an extensive analysis to validate the accuracy, convergence, and effectiveness of the scheme. The comparison of the standard numerical solution and the analytic solution establishes the correctness of the proposed designed methodologies.展开更多
Most problems in uncertainty quantification,despite their ubiquitousness in scientific computing,applied mathematics and data science,remain formidable on a classical computer.For uncertainties that arise in partial d...Most problems in uncertainty quantification,despite their ubiquitousness in scientific computing,applied mathematics and data science,remain formidable on a classical computer.For uncertainties that arise in partial differential equations(PDEs),large numbers M>>1 of samples are required to obtain accurate ensemble averages.This usually involves solving the PDE M times.In addition,to characterise the stochasticity in a PDE,the dimension L of the random input variables is high in most cases,and classical algorithms suffer from the curse-of-dimensionality.We propose new quantum algorithms for PDEs with uncertain coefficients that are more efficient in M and L in various important regimes,compared to their classical counterparts.We introduce transformations that convert the original d-dimensional equation(with uncertain coefficients)into d+L(for dissipative equations)or d+2L(for wave type equations)dimensional equations(with certain coefficients)in which the uncertainties appear only in the initial data.These transformations also allow one to superimpose the M different initial data,so the computational cost for the quantum algorithm to obtain the ensemble average from M different samples is independent of M,while also showing potential advantage in d,L and precisionεin computing ensemble averaged solutions or physical observables.展开更多
Based on a new linear, continuous and bounded operator (PGOPO), a more effective approach and optimal control algorithm than by the block-pulse functions and Walsh functions to design the linear servomechanism of time...Based on a new linear, continuous and bounded operator (PGOPO), a more effective approach and optimal control algorithm than by the block-pulse functions and Walsh functions to design the linear servomechanism of time-varying systems with time-delay is proposed in the paper. By means of the operator, the differential equation is transferred to a more explicit algebraic form which is much easier than the numerical integration of nonlinear TPBVP derived from Pantryagin's maximum principle method. Furthermore, the method is established strictly based on the theory of convergence in the mean square and it is convenient and simple in computation. So the method can be applied to industry control and aeronautics and astronautics field which is frequently mixed with time varying and time delay. Some illustrative numerical examples are interpreted to support the technique.展开更多
We introduce a new parallel evolutionary algorithm in modeling dynamic systems by nonlinear higher-order ordinary differential equations (NHODEs). The NHODEs models are much more universal than the traditional linear ...We introduce a new parallel evolutionary algorithm in modeling dynamic systems by nonlinear higher-order ordinary differential equations (NHODEs). The NHODEs models are much more universal than the traditional linear models. In order to accelerate the modeling process, we propose and realize a parallel evolutionary algorithm using distributed CORBA object on the heterogeneous networking. Some numerical experiments show that the new algorithm is feasible and efficient.展开更多
文摘This study presents numerical algorithms for solving a class of equations that partly consists of derivatives of the unknown state at previous certain times, as well as an integro-differential term containing a weakly singular kernel. These equations are types of integro-differential equation of the second kind and were originally obtained from an aeroelasticity problem. One of the main contributions of this study is to propose numerical algorithms that do not involve transforming the original equation into the corresponding Volterra equation, but still enable the numerical solution of the original equation to be determined. The feasibility of the proposed numerical algorithm is demonstrated by applying examples in measuring the maximum errors with exact solutions at every computed nodes and calculating the corresponding numerical rates of convergence thereafter.
文摘In this paper, a modified algorithm is proposed for solving linear integro-differential equations of the second kind. The main idea is based on applying Romberg extrapolation algorithm (REA), on Trapezoidal rule. In accordance with the computational perspective, the comparison has shown that Adomian decomposition approach is more effective to be utilized. The numerical results show that the modified algorithm has been successfully applied to the linear integro-differential equations and the comparisons with some existing methods appeared in the literature reveal that the modified algorithm is more accurate and convenient.
文摘In this paper, the authors show that the general linear second order ordinary Differential Equation can be formulated as an optimization problem and that evolutionary algorithms for solving optimization problems can also be adapted for solving the formulated problem. The authors propose a polynomial based scheme for achieving the above objectives. The coefficients of the proposed scheme are approximated by an evolutionary algorithm known as Differential Evolution (DE). Numerical examples with good results show the accuracy of the proposed method compared with some existing methods.
文摘The idea of AC = BD was applied to solve the nonlinear differential equations. Suppose that Au = 0 is a given equation to he solved and Dv = 0 is an equation to be easily solved. If the transformation u = Cv is obtained so that v satisfies Dv = 0, then the solutions for Au = 0 can be found. In order to illustrate this approach, several examples about the transformation C are given.
文摘An accurate and rapid method for solving radiative transfer equation is presented in this paper. According to the fact that the multiple scattering component of radiance is less sensitive to the error of phase function than the single scattering component is,we calculate the multiple scattering component by using delta-Eddington approximation and the single scattering component by solving radiative transfer equation. On the ground, when multiple sattering component is small, for example, when the total optical depth T is small, the accurate radiance can be obtained with this method. For the need of the space remote sensing, the upward radiance at the top of the atmosphere is mainly studied, and an approximate expression is presented to correct the multiple scattering component. Compared with the more precise Gauss-Seidel method.the results from this method show an accuracy of better than 10% when zenith angle 0 < 50 掳 and T < 1. The computational speed of this method is, however, much faster than that of Gauss-Seidel method.
文摘When soldering electronic components onto circuit boards,the temperature curves of the reflow ovens across different zones and the conveyor belt speed significantly influence the product quality.This study focuses on optimizing the furnace temperature curve under varying settings of reflow oven zone temperatures and conveyor belt speeds.To address this,the research sequentially develops a heat transfer model for reflow soldering,an optimization model for reflow furnace conditions using the differential evolution algorithm,and an evaluation and decision model combining the differential evolution algorithm with the Technique for Order Preference by Similarity to Ideal Solution(TOPSIS)method.This approach aims to determine the optimal furnace temperature curve,zone temperatures of the reflow oven,and the conveyor belt speed.
文摘A variation of the direct Taylor expansion algorithm is suggested and applied to several linear and nonlinear differential equations of interest in physics and engineering, and the results are compared with those obtained from other algorithms. It is shown that the suggested algorithm competes strongly with other existing algorithms, both in accuracy and ease of application, while demanding a shorter computation time.
文摘This paper proposed the scheme of transmission lines distance protection based on differential equation algorithms (DEA) and Hilbert-Huang transform (HHT). The measured impedance based on EDA is affected by various factors, such as the distributed capacitance, the transient response characteristics of current transformer and voltage transformer, etc. In order to overcome this problem, the proposed scheme applies HHT to improve the apparent impedance estimated by DEA. Empirical mode decomposition (EMD) is used to decompose the data set from DEA into the intrinsic mode functions (IMF) and the residue. This residue has monotonic trend and is used to evaluate the impedance of faulty line. Simulation results show that the proposed scheme improves significantly the accuracy of the estimated impedance.
文摘The multisplitting algorithm for solving large systems of ordinary differential equations on parallel computers was introduced by Jeltsch and Pohl in [1]. On fixed time intervals conver gence results could be derived if the subsystems are solving exactly.Firstly,in theis paper,we deal with an extension of the waveform relaxation algorithm by us ing multisplittin AOR method based on an overlapping block decomposition. We restricted our selves to equidistant timepoints and dealed with the case that an implicit integration method was used to solve the subsystems numerically in parallel. Then we have proved convergence of multi splitting AOR waveform relaxation algorithm on a fixed window containing a finite number of timepoints.
基金National Natural Science Foundation of China(No.61862048)。
文摘Lie algorithm combined with differential form Wu's method is used to complete the symmetry classification of partial differential equations(PDEs)containing arbitrary parameter.This process can be reduced to solve a large system of determining equations,which seems rather difficult to solve,then the differential form Wu's method is used to decompose the determining equations into a series of equations,which are easy to solve.To illustrate the usefulness of this method,we apply it to some test problems,and the results show the performance of the present work.
基金Natural Science Foundation of Shaanxi Province(No.2011J2009)。
文摘In order to solve the problem of vibration bounce caused by the contact between moving and stationary contacts in the process of switching on,two-degree-of-freedom motion differential equation of the contact system is established.Genetic algorithm is used to optimize the pull in process of AC contactor.The whole process of contact bounce was observed and analyzed by high-speed photography experiment.The theory and experimental results were very similar.The iron core has collided before the contact is separated,which further aggravates the contact bounce.When the iron core bounces collided again,the bounce of the contact was not affected.During the operation of the contactor,the movement of the moving iron core will cause slight vibration of the system.The contact bounce time and the maximum amplitude are reduced.The research results provide a theoretical basis for further control and reduction of contact bounce.
文摘Many engineering problems can be reduced to the solution of a variable coefficient differential equation. In this paper, the exact analytic method is suggested to solve variable coefficient differential equations under arbitrary boundary condition. By this method, the general computation formal is obtained. Its convergence in proved. We can get analytic expressions which converge to exact solution and its higher order derivatives uniformy Four numerical examples are given, which indicate that satisfactory results can he obtanedby this method.
文摘New form of necessary conditions for optimality (NCO) is considered. They can be useful for design the direct infinite- dimensional optimization algorithms for systems described by partial differential equations (PDE). Appropriate algo-rithms for unconstrained minimizing a functional are considered and tested. To construct the algorithms, new form of NCO is used. Such approach demonstrates fast uniform convergence at optimal solution in infinite-dimensional space.
基金This work was supported by the DFG through HE 4858/4-1
文摘We extend LeVeque's wave propagation algorithm,a widely used finite volume method for hyperbolic partial differential equations,to a third-order accurate method.The resulting scheme shares main properties with the original method,i.e.,it is based on a wave decomposition at grid cell interfaces,it can be used to approximate hyperbolic problems in divergence form as well as in quasilinear form and limiting is introduced in the form of a wave limiter.
文摘Acceptable glycemic control when examining the effects of meals was </span></span><span><span><span style="font-family:"">achieved when combining basal insulin therapy and high concentration insulin injection before a meal, when using a PID controller (Proportionality, Integrity and Derivative actions) alone, when using a PID controller with basal insulin therapy and when combining the three methods of insulin delivery. Naturally, a type 1 diabetic must inject himself with insulin in well-measured doses. Thus, the management and control of diabetes become a complex task when one must be considered the disturbance due to nutrition and sports activity. This concern has been at the center of much research through different approaches through mathematical methods and Artificial Intelligence methods. This article simulates a physiological model of glycemic control in type 1 diabetics by a PID regulatory mechanism, in the context of disturbances caused by the patient’s meals and athletic activity.
文摘A heuristic technique is developed for a nonlinear magnetohydrodynamics （MHD） Jeffery-Hamel problem with the help of the feed-forward artificial neural net- work （ANN） optimized with the genetic algorithm （GA） and the sequential quadratic programming （SQP） method. The twodimensional （2D） MHD Jeffery-Hamel problem is transformed into a higher order boundary value problem （BVP） of ordinary differential equations （ODEs）. The mathematical model of the transformed BVP is formulated with the ANN in an unsupervised manner. The training of the weights of the ANN is carried out with the evolutionary calculation based on the GA hybridized with the SQP method for the rapid local convergence. The proposed scheme is evaluated on the variants of the Jeffery-Hamel flow by varying the Reynold number, the Hartmann number, and the an- gles of the walls. A large number of simulations are performed with an extensive analysis to validate the accuracy, convergence, and effectiveness of the scheme. The comparison of the standard numerical solution and the analytic solution establishes the correctness of the proposed designed methodologies.
基金supported by the National Natural Science Foundation of China(Grant Nos.12031013,12341104,and 12050410230)the National Natural Science Foundation of China International Young Scientists Project(Grant No.12050410230)+6 种基金the Shanghai Municipal Science and Technology Major Project(Grant No.2021SHZDZX0102)the Innovation Program of Shanghai Municipal Education Commission(Grant No.2021-01-07-00-02-E00087)the Science and Technology Program of ShanghaiChina(Grant No.21JC1402900)the Shanghai Pujiang Talent Grant(Grant No.20PJ1408400)the Shanghai Jiao Tong University 2030 Initiativethe Fundamental Research Funds for the Central Universities。
文摘Most problems in uncertainty quantification,despite their ubiquitousness in scientific computing,applied mathematics and data science,remain formidable on a classical computer.For uncertainties that arise in partial differential equations(PDEs),large numbers M>>1 of samples are required to obtain accurate ensemble averages.This usually involves solving the PDE M times.In addition,to characterise the stochasticity in a PDE,the dimension L of the random input variables is high in most cases,and classical algorithms suffer from the curse-of-dimensionality.We propose new quantum algorithms for PDEs with uncertain coefficients that are more efficient in M and L in various important regimes,compared to their classical counterparts.We introduce transformations that convert the original d-dimensional equation(with uncertain coefficients)into d+L(for dissipative equations)or d+2L(for wave type equations)dimensional equations(with certain coefficients)in which the uncertainties appear only in the initial data.These transformations also allow one to superimpose the M different initial data,so the computational cost for the quantum algorithm to obtain the ensemble average from M different samples is independent of M,while also showing potential advantage in d,L and precisionεin computing ensemble averaged solutions or physical observables.
基金National Natural Science Foundation of China(69934010)
文摘Based on a new linear, continuous and bounded operator (PGOPO), a more effective approach and optimal control algorithm than by the block-pulse functions and Walsh functions to design the linear servomechanism of time-varying systems with time-delay is proposed in the paper. By means of the operator, the differential equation is transferred to a more explicit algebraic form which is much easier than the numerical integration of nonlinear TPBVP derived from Pantryagin's maximum principle method. Furthermore, the method is established strictly based on the theory of convergence in the mean square and it is convenient and simple in computation. So the method can be applied to industry control and aeronautics and astronautics field which is frequently mixed with time varying and time delay. Some illustrative numerical examples are interpreted to support the technique.
基金the National Natural Science Foundation of China(No.70 0 710 42 and No.60 0 73 0 43 )
文摘We introduce a new parallel evolutionary algorithm in modeling dynamic systems by nonlinear higher-order ordinary differential equations (NHODEs). The NHODEs models are much more universal than the traditional linear models. In order to accelerate the modeling process, we propose and realize a parallel evolutionary algorithm using distributed CORBA object on the heterogeneous networking. Some numerical experiments show that the new algorithm is feasible and efficient.
